23916574395

Committed 02 Apr 2026 06:50PM UTC coverage: 81.336% (+0.01%) from 81.324%

Build # 23916574395

Build Type

Pull #3734

github

Committed by

web-flow

Commit Message

Merge 3420117d8 into d9b30bbbd

Pull Request Pull Request #3734: Specify temperature from a field (structured mesh only)

Coverage Stats

17720 of 25601 branches covered (69.22%)

Branch coverage included in aggregate %.

183 of 204 new or added lines in 15 files covered. (89.71%)

69 existing lines in 3 files now uncovered.

58284 of 67843 relevant lines covered (85.91%)

44753955.49 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

83.33

/include/openmc/field.h

#ifndef OPENMC_FIELD_H
#define OPENMC_FIELD_H

#include "openmc/mesh.h"
#include "openmc/vector.h"

namespace openmc {

class ScalarField {
public:
  //----------------------------------------------------------------------------
  // Constructors
  ScalarField() = default;
  ScalarField(
    Mesh* mesh_ptr, vector<double> values, const std::string& field_name);
  ScalarField(Mesh* mesh_ptr, vector<double> values)
    : ScalarField(mesh_ptr, values, "ScalarField") {};

  //----------------------------------------------------------------------------
  // Methods

  //! Returns the distance to the next mesh boundary gien a particle position
  //! and direction. If the particle is initially outside, the distance will
  //! correspond to the nearest distance to the outer boundaries of the mesh.
  //
  //! \param[in] r Position of the particle
  //! \param[in] u Direction of the particle
  //! \return The distance in cm to the next mesh boundary
  double distance_to_next_boundary(const Position& r, const Direction& d);

  //----------------------------------------------------------------------------
  // Accessors

  // Field type
  const std::string& field_type() const { return this->field_type_; }

  // Mesh pointer
  Mesh* mesh_ptr() const
  {
    if (this->mesh_ptr_ == nullptr) {
      fatal_error(fmt::format("No mesh found for {}!", this->field_type_));
    } else {
      return this->mesh_ptr_;
    }
  }

  // Values
  double& value(int i) { return values_[i]; }
  const double& value(int i) const { return values_[i]; }
  const vector<double>& values() const { return values_; }

private:
  //----------------------------------------------------------------------------
  // Data members
  std::string field_type_; //! Name of field type
  Mesh* mesh_ptr_;         //!< Pointer to the geometric mesh
  vector<double> values_;  //!< Values associated with each mesh cell
};

class TemperatureField : public ScalarField {
public:
  //----------------------------------------------------------------------------
  // Constructors
  TemperatureField() = default;
  TemperatureField(Mesh* mesh_ptr, vector<double> values)
    : ScalarField(mesh_ptr, values, "TemperatureField") {};

  //----------------------------------------------------------------------------
  // Methods

  //! Returns the temperature in Kelvin corresponding to the given position.
  //
  //! \param[in] r Position of the particle
  //! \return Temperature in Kelvin
  double get_temperature(const Position& r);

  //! Returns the square root of the temperature multiplied by the Boltzmann
  //! constant in eV for the given position.
  //
  //! \param[in] r Position of the particle
  //! \return Sqrt(k_Boltzmann * temperature) in eV
  double get_sqrtkT(const Position& r);

  //! Update the temperature of a particle based on its position and direction.
  //! If the particle is inside the temperature field, its temperature is
  //! updated. If outside, the particle takes the temperature value
  //! associated with the current cell instance.
  //
  //! \param[inout] p Particle
  void update_particle_temperature(Particle& p);
};

} //  namespace openmc
#endif // OPENMC_FIELD_H

1	#ifndef OPENMC_FIELD_H
2	#define OPENMC_FIELD_H
3
4	#include "openmc/mesh.h"
5	#include "openmc/vector.h"
6
7	namespace openmc {
8
9	class ScalarField {	8,350✔
10	public:
11	//----------------------------------------------------------------------------
12	// Constructors
13	ScalarField() = default;	8,335✔
14	ScalarField(
15	Mesh* mesh_ptr, vector<double> values, const std::string& field_name);
16	ScalarField(Mesh* mesh_ptr, vector<double> values)
17	: ScalarField(mesh_ptr, values, "ScalarField") {};
18
19	//----------------------------------------------------------------------------
20	// Methods
21
22	//! Returns the distance to the next mesh boundary gien a particle position
23	//! and direction. If the particle is initially outside, the distance will
24	//! correspond to the nearest distance to the outer boundaries of the mesh.
25	//
26	//! \param[in] r Position of the particle
27	//! \param[in] u Direction of the particle
28	//! \return The distance in cm to the next mesh boundary
29	double distance_to_next_boundary(const Position& r, const Direction& d);
30
31	//----------------------------------------------------------------------------
32	// Accessors
33
34	// Field type
35	const std::string& field_type() const { return this->field_type_; }
36
37	// Mesh pointer
38	Mesh* mesh_ptr() const	2,027,157✔
39	{
40	if (this->mesh_ptr_ == nullptr) {	2,027,157!
NEW 41	fatal_error(fmt::format("No mesh found for {}!", this->field_type_));	×
42	} else {
43	return this->mesh_ptr_;	2,027,157✔
44	}
45	}
46
47	// Values
48	double& value(int i) { return values_[i]; }	663,080✔
49	const double& value(int i) const { return values_[i]; }
50	const vector<double>& values() const { return values_; }	663,080!
51
52	private:
53	//----------------------------------------------------------------------------
54	// Data members
55	std::string field_type_; //! Name of field type
56	Mesh* mesh_ptr_; //!< Pointer to the geometric mesh
57	vector<double> values_; //!< Values associated with each mesh cell
58	};
59
60	class TemperatureField : public ScalarField {	16,700✔
61	public:
62	//----------------------------------------------------------------------------
63	// Constructors
64	TemperatureField() = default;	8,335✔
65	TemperatureField(Mesh* mesh_ptr, vector<double> values)	15✔
66	: ScalarField(mesh_ptr, values, "TemperatureField") {};	30✔
67
68	//----------------------------------------------------------------------------
69	// Methods
70
71	//! Returns the temperature in Kelvin corresponding to the given position.
72	//
73	//! \param[in] r Position of the particle
74	//! \return Temperature in Kelvin
75	double get_temperature(const Position& r);
76
77	//! Returns the square root of the temperature multiplied by the Boltzmann
78	//! constant in eV for the given position.
79	//
80	//! \param[in] r Position of the particle
81	//! \return Sqrt(k_Boltzmann * temperature) in eV
82	double get_sqrtkT(const Position& r);
83
84	//! Update the temperature of a particle based on its position and direction.
85	//! If the particle is inside the temperature field, its temperature is
86	//! updated. If outside, the particle takes the temperature value
87	//! associated with the current cell instance.
88	//
89	//! \param[inout] p Particle
90	void update_particle_temperature(Particle& p);
91	};
92
93	} // namespace openmc
94	#endif // OPENMC_FIELD_H

openmc-dev / openmc / 23916574395

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous