OGS
ProcessLib::HeatTransportBHE::BHE::BHE_1U Class Referencefinal

Detailed Description

The BHE_1U class is the realization of 1U type of Borehole Heat Exchanger. In this class, the pipe heat capacity, pipe heat conduction, pipe advection vectors are initialized according to the geometry of 1U type of BHE. For 1U type of BHE, 4 primary unknowns are assigned on the 1D BHE elements. They are the temperature in inflow pipe T_in, temperature in outflow pipe T_out temperature of the two grout zones surrounding the inflow and outflow pipe T_g1 and T_g2. These primary variables are solved according to heat convection and conduction equations on the pipes and also in the grout zones. The interaction of the 1U type of BHE and the surrounding soil is regulated through the thermal resistance values, which are calculated specifically during the initialization of the class.

Definition at line 39 of file BHE_1U.h.

#include <BHE_1U.h>

Inheritance diagram for ProcessLib::HeatTransportBHE::BHE::BHE_1U:
[legend]
Collaboration diagram for ProcessLib::HeatTransportBHE::BHE::BHE_1U:
[legend]

Public Member Functions

 BHE_1U (BoreholeGeometry const &borehole, RefrigerantProperties const &refrigerant, GroutParameters const &grout, FlowAndTemperatureControl const &flowAndTemperatureControl, PipeConfigurationUType const &pipes, bool const use_python_bcs)
std::array< double, number_of_flow_legs > flowLegs () const
std::array< double, number_of_unknowns > pipeHeatCapacities () const
std::array< double, number_of_unknowns > pipeHeatConductions () const
std::array< Eigen::Vector3d, number_of_unknowns > pipeAdvectionVectors (Eigen::Vector3d const &elem_direction) const
std::vector< double > thermalResistances (ParameterLib::SpatialPosition const &pos) const
double updateFlowRateAndTemperature (double T_out, double current_time)
 Return the inflow temperature for the boundary condition.
std::array< double, number_of_unknowns > crossSectionAreas (ParameterLib::SpatialPosition const &pos) const
void updateHeatTransferCoefficients (double const flow_rate)
Public Member Functions inherited from ProcessLib::HeatTransportBHE::BHE::BHECommonUType
 BHECommonUType (BoreholeGeometry const &borehole, RefrigerantProperties const &refrigerant, GroutParameters const &grout, FlowAndTemperatureControl const &flowAndTemperatureControl, PipeConfigurationUType const &pipes, bool const use_python_bcs)
Public Member Functions inherited from ProcessLib::HeatTransportBHE::BHE::BHECommon
 BHECommon (BoreholeGeometry const &borehole_geometry_, RefrigerantProperties const &refrigerant_, GroutParameters const &grout_, FlowAndTemperatureControl const &flowAndTemperatureControl_, bool const use_python_bcs_)
constexpr bool isPowerBC () const

Static Public Member Functions

template<int NPoints, typename SingleUnknownMatrixType, typename RMatrixType, typename RPiSMatrixType, typename RSMatrixType>
static void assembleRMatrices (int const idx_bhe_unknowns, Eigen::MatrixBase< SingleUnknownMatrixType > const &matBHE_loc_R, Eigen::MatrixBase< RMatrixType > &R_matrix, Eigen::MatrixBase< RPiSMatrixType > &R_pi_s_matrix, Eigen::MatrixBase< RSMatrixType > &R_s_matrix)
static std::array< std::pair< std::size_t, int >, 2 > getBHEInflowDirichletBCNodesAndComponents (std::size_t const top_node_id, std::size_t const, int const in_component_id)
static std::optional< std::array< std::pair< std::size_t, int >, 2 > > getBHEBottomDirichletBCNodesAndComponents (std::size_t const bottom_node_id, int const in_component_id, int const out_component_id)

Static Public Attributes

static constexpr int number_of_unknowns = 4
static constexpr int number_of_grout_zones = 2
static constexpr int number_of_flow_legs = 2
static constexpr std::pair< int, int > inflow_outflow_bc_component_ids []

Private Member Functions

std::vector< double > calcThermalResistances (double const Nu, ParameterLib::SpatialPosition const &pos) const
 Nu is the Nusselt number.

Private Attributes

double cached_nu_ = 0.0

Additional Inherited Members

Public Attributes inherited from ProcessLib::HeatTransportBHE::BHE::BHECommon
BoreholeGeometry const borehole_geometry
RefrigerantProperties const refrigerant
GroutParameters const grout
FlowAndTemperatureControl const flowAndTemperatureControl
bool const use_python_bcs
Protected Attributes inherited from ProcessLib::HeatTransportBHE::BHE::BHECommonUType
PipeConfigurationUType const _pipes
double flow_velocity = 0.0

Constructor & Destructor Documentation

◆ BHE_1U()

ProcessLib::HeatTransportBHE::BHE::BHE_1U::BHE_1U ( BoreholeGeometry const & borehole,
RefrigerantProperties const & refrigerant,
GroutParameters const & grout,
FlowAndTemperatureControl const & flowAndTemperatureControl,
PipeConfigurationUType const & pipes,
bool const use_python_bcs )

Definition at line 20 of file BHE_1U.cpp.

26 : BHECommonUType{borehole, refrigerant, grout, flowAndTemperatureControl,
27 pipes, use_python_bcs}
28{
29 checkEqualPipeOutsideDiameters(_pipes.inlet, _pipes.outlet, "BHE 1U");
30 // Initialize thermal resistances.
31 auto values = visit(
32 [&](auto const& control)
33 {
34 return control(refrigerant.reference_temperature,
35 0. /* initial time */);
36 },
37 flowAndTemperatureControl);
38 updateHeatTransferCoefficients(values.flow_rate);
39}
ProcessLib::HeatTransportBHE::BHE::BHECommonUType::BHECommonUType
BHECommonUType(BoreholeGeometry const &borehole, RefrigerantProperties const &refrigerant, GroutParameters const &grout, FlowAndTemperatureControl const &flowAndTemperatureControl, PipeConfigurationUType const &pipes, bool const use_python_bcs)
Definition BHECommonUType.h:19
ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant
RefrigerantProperties const refrigerant
Definition BHECommon.h:51
ProcessLib::HeatTransportBHE::BHE::BHECommon::flowAndTemperatureControl
FlowAndTemperatureControl const flowAndTemperatureControl
Definition BHECommon.h:53
ProcessLib::HeatTransportBHE::BHE::BHE_1U::updateHeatTransferCoefficients
void updateHeatTransferCoefficients(double const flow_rate)
Definition BHE_1U.cpp:134
ProcessLib::HeatTransportBHE::BHE::checkEqualPipeOutsideDiameters
void checkEqualPipeOutsideDiameters(Pipe const &inlet, Pipe const &outlet, std::string_view const context)
Definition BHEParameterValidation.cpp:92

References ProcessLib::HeatTransportBHE::BHE::BHECommonUType::BHECommonUType(), ProcessLib::HeatTransportBHE::BHE::BHECommonUType::_pipes, ProcessLib::HeatTransportBHE::BHE::checkEqualPipeOutsideDiameters(), ProcessLib::HeatTransportBHE::BHE::BHECommon::flowAndTemperatureControl, ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant, updateHeatTransferCoefficients(), and ProcessLib::HeatTransportBHE::BHE::BHECommon::use_python_bcs.

Member Function Documentation

◆ assembleRMatrices()

template<int NPoints, typename SingleUnknownMatrixType, typename RMatrixType, typename RPiSMatrixType, typename RSMatrixType>
void ProcessLib::HeatTransportBHE::BHE::BHE_1U::assembleRMatrices ( int const idx_bhe_unknowns,
Eigen::MatrixBase< SingleUnknownMatrixType > const & matBHE_loc_R,
Eigen::MatrixBase< RMatrixType > & R_matrix,
Eigen::MatrixBase< RPiSMatrixType > & R_pi_s_matrix,
Eigen::MatrixBase< RSMatrixType > & R_s_matrix )
inlinestatic

Definition at line 80 of file BHE_1U.h.

86 {
87 switch (idx_bhe_unknowns)
88 {
89 case 0: // PHI_fig
90 R_matrix.block(0, 2 * NPoints, NPoints, NPoints) +=
91 -1.0 * matBHE_loc_R;
92 R_matrix.block(2 * NPoints, 0, NPoints, NPoints) +=
93 -1.0 * matBHE_loc_R;
94
95 R_matrix.block(0, 0, NPoints, NPoints) +=
96 1.0 * matBHE_loc_R; // K_i1
97 R_matrix.block(2 * NPoints, 2 * NPoints, NPoints, NPoints) +=
98 1.0 * matBHE_loc_R; // K_ig
99 return;
100 case 1: // PHI_fog
101 R_matrix.block(NPoints, 3 * NPoints, NPoints, NPoints) +=
102 -1.0 * matBHE_loc_R;
103 R_matrix.block(3 * NPoints, NPoints, NPoints, NPoints) +=
104 -1.0 * matBHE_loc_R;
105
106 R_matrix.block(NPoints, NPoints, NPoints, NPoints) +=
107 1.0 * matBHE_loc_R; // K_o1
108 R_matrix.block(3 * NPoints, 3 * NPoints, NPoints, NPoints) +=
109 1.0 * matBHE_loc_R; // K_og
110 return;
111 case 2: // PHI_gg
112 R_matrix.block(2 * NPoints, 3 * NPoints, NPoints, NPoints) +=
113 -1.0 * matBHE_loc_R;
114 R_matrix.block(3 * NPoints, 2 * NPoints, NPoints, NPoints) +=
115 -1.0 * matBHE_loc_R;
116
117 R_matrix.block(2 * NPoints, 2 * NPoints, NPoints, NPoints) +=
118 1.0 * matBHE_loc_R; // K_ig // notice we only have
119 // 1 PHI_gg term here.
120 R_matrix.block(3 * NPoints, 3 * NPoints, NPoints, NPoints) +=
121 1.0 * matBHE_loc_R; // K_og // see Diersch 2013 FEFLOW
122 // book page 954 Table M.2
123 return;
124 case 3: // PHI_gs
125 R_s_matrix.template block<NPoints, NPoints>(0, 0).noalias() +=
126 1.0 * matBHE_loc_R;
127
128 R_pi_s_matrix.block(2 * NPoints, 0, NPoints, NPoints) +=
129 -1.0 * matBHE_loc_R;
130 R_pi_s_matrix.block(3 * NPoints, 0, NPoints, NPoints) +=
131 -1.0 * matBHE_loc_R;
132 R_matrix.block(2 * NPoints, 2 * NPoints, NPoints, NPoints) +=
133 1.0 * matBHE_loc_R; // K_ig
134 R_matrix.block(3 * NPoints, 3 * NPoints, NPoints, NPoints) +=
135 1.0 * matBHE_loc_R; // K_og
136 return;
137 default:
138 OGS_FATAL(
139 "BHE_1U::assembleRMatrices: unknown index {:d} "
140 "out of range.",
141 idx_bhe_unknowns);
142 }
143 }
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:10

References OGS_FATAL.

◆ calcThermalResistances()

std::vector< double > ProcessLib::HeatTransportBHE::BHE::BHE_1U::calcThermalResistances ( double const Nu,
ParameterLib::SpatialPosition const & pos ) const
private

Nu is the Nusselt number.

Definition at line 150 of file BHE_1U.cpp.

152{
153 constexpr double pi = std::numbers::pi;
154
155 double const lambda_r = refrigerant.thermal_conductivity;
156 double const lambda_g = grout.lambda_g;
157 // t=0.0: borehole properties are physically time-invariant; genuinely
158 // time-varying parameter types are rejected in createPipe.
159 double const lambda_p_inlet =
160 sampleStrictPositive(_pipes.inlet.wall_thermal_conductivity, 0.0, pos,
161 "inlet wall_thermal_conductivity");
162 double const lambda_p_outlet =
163 sampleStrictPositive(_pipes.outlet.wall_thermal_conductivity, 0.0, pos,
164 "outlet wall_thermal_conductivity");
165
166 // thermal resistances due to advective flow of refrigerant in the _pipes
167 // Eq. 36 in Diersch_2011_CG
168 double const R_adv_i1 = 1.0 / (Nu * lambda_r * pi);
169 double const R_adv_o1 = 1.0 / (Nu * lambda_r * pi);
170
171 // thermal resistance due to thermal conductivity of the pipe wall material
172 // Eq. 49
173 double const inlet_outside_diameter = _pipes.inlet.outsideDiameter();
174 double const R_con_a_inlet =
175 pipeWallThermalResistance(_pipes.inlet, lambda_p_inlet);
176 double const R_con_a_outlet =
177 pipeWallThermalResistance(_pipes.outlet, lambda_p_outlet);
178
179 // the average outer diameter of the _pipes
180 double const d0 = inlet_outside_diameter;
181 double const D = sampleStrictPositive(borehole_geometry.diameter, 0.0, pos,
182 "borehole_diameter");
183 // Context prefix shared by the acosh-argument checks below. `cause` names
184 // the physical degeneracy that drives the argument to <= 1 so the fatal
185 // message is actionable.
186 auto const geometry_context =
187 [&](std::string_view const equation, std::string_view const cause)
188 {
189 return uTypeGeometryContext(pos, D, d0, _pipes.distance, equation,
190 cause);
191 };
192 // Eq. 51: chi requires D > sqrt(2) * d0 so that log(D/(sqrt(2)*d0)) > 0.
193 checkBoreholeVsPipeDiameter(D, std::sqrt(2.0) * d0, pos,
194 "BHE 1U chi formula (Eq. 51)");
195 double const chi = std::log(std::sqrt(D * D + 2 * d0 * d0) / 2 / d0) /
196 std::log(D / std::sqrt(2) / d0);
197 // Eq. 52
198 // thermal resistances of the grout
199 double const acosh_arg_R_g =
200 (D * D + d0 * d0 - _pipes.distance * _pipes.distance) / (2 * D * d0);
201 checkAcoshArg(
202 acosh_arg_R_g,
203 geometry_context(
204 "BHE 1U R_g (Eq. 52)",
205 "The argument drops to <= 1 when a pipe reaches the borehole wall "
206 "(distance >= D - d0), giving a zero grout resistance R_g and an "
207 "infinite (1/R) assembly coefficient."));
208 double const R_g = std::acosh(acosh_arg_R_g) / (2 * pi * lambda_g) *
209 (1.601 - 0.888 * _pipes.distance / D);
210
211 // thermal resistance due to inter-grout exchange
212 double const acosh_arg_R_ar =
213 (2.0 * _pipes.distance * _pipes.distance - d0 * d0) / d0 / d0;
214 checkAcoshArg(
215 acosh_arg_R_ar,
216 geometry_context(
217 "BHE 1U R_ar",
218 "The argument drops to <= 1 when the pipes touch (distance <= d0), "
219 "giving a zero inter-grout resistance R_ar and an infinite (1/R) "
220 "assembly coefficient."));
221 double const R_ar = std::acosh(acosh_arg_R_ar) / (2.0 * pi * lambda_g);
222
223 auto const [chi_new, R_gg, R_gs] =
224 thermalResistancesGroutSoil(chi, R_ar, R_g);
225
226 // thermal resistance due to the grout transition.
227 double const R_con_b = chi_new * R_g;
228 // Eq. 29 and 30
229 double const R_fig = R_adv_i1 + R_con_a_inlet + R_con_b;
230 double const R_fog = R_adv_o1 + R_con_a_outlet + R_con_b;
231
232 return {R_fig, R_fog, R_gg, R_gs};
233}
ProcessLib::HeatTransportBHE::BHE::checkAcoshArg
void checkAcoshArg(double const arg, std::string_view const context)
Definition BHEParameterValidation.cpp:53
ProcessLib::HeatTransportBHE::BHE::sampleStrictPositive
double sampleStrictPositive(ParameterLib::Parameter< double > const &param, double const t, ParameterLib::SpatialPosition const &pos, std::string_view const param_role)
Definition BHEParameterValidation.cpp:24
ProcessLib::HeatTransportBHE::BHE::thermalResistancesGroutSoil
std::array< double, 3 > thermalResistancesGroutSoil(double const chi, double const R_ar, double const R_g)
Definition BHE_1U.cpp:102
ProcessLib::HeatTransportBHE::BHE::checkBoreholeVsPipeDiameter
void checkBoreholeVsPipeDiameter(double const D, double const min_diameter, ParameterLib::SpatialPosition const &pos, std::string_view const context)
Definition BHEParameterValidation.cpp:62
ProcessLib::HeatTransportBHE::BHE::uTypeGeometryContext
std::string uTypeGeometryContext(ParameterLib::SpatialPosition const &pos, double const D, double const d0, double const distance, std::string_view const equation, std::string_view const cause)
Definition BHEParameterValidation.cpp:40
ProcessLib::HeatTransportBHE::BHE::pipeWallThermalResistance
double pipeWallThermalResistance(Pipe const &pipe, double const wall_thermal_conductivity)
Definition ThermalResistanceHelpers.h:22

References ProcessLib::HeatTransportBHE::BHE::BHECommonUType::_pipes, ProcessLib::HeatTransportBHE::BHE::BHECommon::borehole_geometry, ProcessLib::HeatTransportBHE::BHE::checkAcoshArg(), ProcessLib::HeatTransportBHE::BHE::checkBoreholeVsPipeDiameter(), ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, ProcessLib::HeatTransportBHE::BHE::pipeWallThermalResistance(), ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant, ProcessLib::HeatTransportBHE::BHE::sampleStrictPositive(), ProcessLib::HeatTransportBHE::BHE::thermalResistancesGroutSoil(), and ProcessLib::HeatTransportBHE::BHE::uTypeGeometryContext().

Referenced by thermalResistances().

◆ crossSectionAreas()

std::array< double, BHE_1U::number_of_unknowns > ProcessLib::HeatTransportBHE::BHE::BHE_1U::crossSectionAreas ( ParameterLib::SpatialPosition const & pos) const

Definition at line 256 of file BHE_1U.cpp.

258{
259 double const D = sampleStrictPositive(borehole_geometry.diameter, 0.0, pos,
260 "borehole_diameter");
261 double const borehole_area = Pipe::circleArea(D);
262 double const half_borehole_area = borehole_area / number_of_grout_zones;
263 return {
264 {_pipes.inlet.area(), _pipes.outlet.area(),
265 checkedGroutArea(half_borehole_area, _pipes.inlet.outsideArea(), pos),
266 checkedGroutArea(half_borehole_area, _pipes.outlet.outsideArea(),
267 pos)}};
268}
ProcessLib::HeatTransportBHE::BHE::BHE_1U::number_of_grout_zones
static constexpr int number_of_grout_zones
Definition BHE_1U.h:50
ProcessLib::HeatTransportBHE::BHE::checkedGroutArea
double checkedGroutArea(double const borehole_area_fraction, double const pipe_outside_area, ParameterLib::SpatialPosition const &pos)
Definition BHEParameterValidation.cpp:77
ProcessLib::HeatTransportBHE::BHE::Pipe::circleArea
static double circleArea(double const d)
Definition Pipe.h:46

References ProcessLib::HeatTransportBHE::BHE::BHECommonUType::_pipes, ProcessLib::HeatTransportBHE::BHE::BHECommon::borehole_geometry, ProcessLib::HeatTransportBHE::BHE::checkedGroutArea(), ProcessLib::HeatTransportBHE::BHE::Pipe::circleArea(), number_of_grout_zones, and ProcessLib::HeatTransportBHE::BHE::sampleStrictPositive().

◆ flowLegs()

std::array< double, number_of_flow_legs > ProcessLib::HeatTransportBHE::BHE::BHE_1U::flowLegs ( ) const
inline

Signed fluid velocity per flow leg. Positive = flow along +elem_direction; negative = against it. Leg ordering matches the BHE unknowns: { i1 (down), o1 (up) }.

Definition at line 56 of file BHE_1U.h.

57 {
58 return {+flow_velocity, -flow_velocity};
59 }

References ProcessLib::HeatTransportBHE::BHE::BHECommonUType::flow_velocity.

Referenced by pipeAdvectionVectors().

◆ getBHEBottomDirichletBCNodesAndComponents()

std::optional< std::array< std::pair< std::size_t, int >, 2 > > ProcessLib::HeatTransportBHE::BHE::BHE_1U::getBHEBottomDirichletBCNodesAndComponents ( std::size_t const bottom_node_id,
int const in_component_id,
int const out_component_id )
static

Definition at line 247 of file BHE_1U.cpp.

251{
252 return {{std::make_pair(bottom_node_id, in_component_id),
253 std::make_pair(bottom_node_id, out_component_id)}};
254}

◆ getBHEInflowDirichletBCNodesAndComponents()

std::array< std::pair< std::size_t, int >, 2 > ProcessLib::HeatTransportBHE::BHE::BHE_1U::getBHEInflowDirichletBCNodesAndComponents ( std::size_t const top_node_id,
std::size_t const ,
int const in_component_id )
static

Definition at line 236 of file BHE_1U.cpp.

240{
241 return {std::make_pair(top_node_id, in_component_id),
242 std::make_pair(top_node_id, in_component_id + 1)};
243}

◆ pipeAdvectionVectors()

std::array< Eigen::Vector3d, BHE_1U::number_of_unknowns > ProcessLib::HeatTransportBHE::BHE::BHE_1U::pipeAdvectionVectors ( Eigen::Vector3d const & elem_direction) const

Definition at line 82 of file BHE_1U.cpp.

83{
84 double const& rho_r = refrigerant.density;
85 double const& Cp_r = refrigerant.specific_heat_capacity;
86
87 auto const legs = flowLegs();
88 auto leg_adv = [&](double const v_signed) -> Eigen::Vector3d
89 { return rho_r * Cp_r * v_signed * elem_direction; };
90
91 return {{leg_adv(legs[0]), // i1
92 leg_adv(legs[1]), // o1
93 {0, 0, 0}, // g1
94 {0, 0, 0}}}; // g2
95}
ProcessLib::HeatTransportBHE::BHE::BHE_1U::flowLegs
std::array< double, number_of_flow_legs > flowLegs() const
Definition BHE_1U.h:56

References flowLegs(), and ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant.

◆ pipeHeatCapacities()

std::array< double, BHE_1U::number_of_unknowns > ProcessLib::HeatTransportBHE::BHE::BHE_1U::pipeHeatCapacities ( ) const

Definition at line 41 of file BHE_1U.cpp.

43{
44 double const rho_r = refrigerant.density;
45 double const specific_heat_capacity = refrigerant.specific_heat_capacity;
46 double const rho_g = grout.rho_g;
47 double const porosity_g = grout.porosity_g;
48 double const heat_cap_g = grout.heat_cap_g;
49
50 return {{/*i1*/ rho_r * specific_heat_capacity,
51 /*o1*/ rho_r * specific_heat_capacity,
52 /*g1*/ (1.0 - porosity_g) * rho_g * heat_cap_g,
53 /*g2*/ (1.0 - porosity_g) * rho_g * heat_cap_g}};
54}

References ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, and ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant.

◆ pipeHeatConductions()

std::array< double, BHE_1U::number_of_unknowns > ProcessLib::HeatTransportBHE::BHE::BHE_1U::pipeHeatConductions ( ) const

Definition at line 56 of file BHE_1U.cpp.

58{
59 double const lambda_r = refrigerant.thermal_conductivity;
60 double const rho_r = refrigerant.density;
61 double const Cp_r = refrigerant.specific_heat_capacity;
62 double const alpha_L = _pipes.longitudinal_dispersion_length;
63 double const porosity_g = grout.porosity_g;
64 double const lambda_g = grout.lambda_g;
65
66 double const velocity_norm = std::abs(flow_velocity);
67
68 // Here we calculate the laplace coefficients in the governing
69 // equations of BHE. These governing equations can be found in
70 // 1) Diersch (2013) FEFLOW book on page 952, M.120-122, or
71 // 2) Diersch (2011) Comp & Geosci 37:1122-1135, Eq. 19-22.
72 auto const pipe_conduction =
73 lambda_r + rho_r * Cp_r * alpha_L * velocity_norm;
74 auto const grout_conduction = (1.0 - porosity_g) * lambda_g;
75 return {{pipe_conduction, // i1
76 pipe_conduction, // o1
77 grout_conduction, // g1
78 grout_conduction}}; // g2
79}

References ProcessLib::HeatTransportBHE::BHE::BHECommonUType::_pipes, ProcessLib::HeatTransportBHE::BHE::BHECommonUType::flow_velocity, ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, and ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant.

◆ thermalResistances()

std::vector< double > ProcessLib::HeatTransportBHE::BHE::BHE_1U::thermalResistances ( ParameterLib::SpatialPosition const & pos) const

Return the full vector of thermal resistances for the element at pos, computed once using the cached Nusselt number. Intended to be called once per element by the assembler so that the resistance computation is not repeated per unknown and per integration point.

Definition at line 143 of file BHE_1U.cpp.

145{
146 return calcThermalResistances(cached_nu_, pos);
147}
ProcessLib::HeatTransportBHE::BHE::BHE_1U::calcThermalResistances
std::vector< double > calcThermalResistances(double const Nu, ParameterLib::SpatialPosition const &pos) const
Nu is the Nusselt number.
Definition BHE_1U.cpp:150
ProcessLib::HeatTransportBHE::BHE::BHE_1U::cached_nu_
double cached_nu_
Definition BHE_1U.h:169

References cached_nu_, and calcThermalResistances().

◆ updateFlowRateAndTemperature()

double ProcessLib::HeatTransportBHE::BHE::BHE_1U::updateFlowRateAndTemperature ( double T_out,
double current_time )

Return the inflow temperature for the boundary condition.

Definition at line 270 of file BHE_1U.cpp.

272{
273 auto values =
274 visit([&](auto const& control) { return control(T_out, current_time); },
275 flowAndTemperatureControl);
276 updateHeatTransferCoefficients(values.flow_rate);
277 return values.temperature;
278}

References ProcessLib::HeatTransportBHE::BHE::BHECommon::flowAndTemperatureControl, and updateHeatTransferCoefficients().

◆ updateHeatTransferCoefficients()

void ProcessLib::HeatTransportBHE::BHE::BHE_1U::updateHeatTransferCoefficients ( double const flow_rate)

Definition at line 134 of file BHE_1U.cpp.

135{
136 auto const tm_flow = calculateThermoMechanicalFlowPropertiesPipe(
137 _pipes.inlet, borehole_geometry.length, refrigerant, flow_rate);
138
139 flow_velocity = tm_flow.velocity;
140 cached_nu_ = tm_flow.nusselt_number;
141}
ProcessLib::HeatTransportBHE::BHE::calculateThermoMechanicalFlowPropertiesPipe
ThermoMechanicalFlowProperties calculateThermoMechanicalFlowPropertiesPipe(Pipe const &pipe, double const length, RefrigerantProperties const &fluid, double const flow_rate)
Definition ThermoMechanicalFlowProperties.h:23

References ProcessLib::HeatTransportBHE::BHE::BHECommonUType::_pipes, ProcessLib::HeatTransportBHE::BHE::BHECommon::borehole_geometry, cached_nu_, ProcessLib::HeatTransportBHE::BHE::calculateThermoMechanicalFlowPropertiesPipe(), ProcessLib::HeatTransportBHE::BHE::BHECommonUType::flow_velocity, and ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant.

Referenced by BHE_1U(), and updateFlowRateAndTemperature().

Member Data Documentation

◆ cached_nu_

double ProcessLib::HeatTransportBHE::BHE::BHE_1U::cached_nu_ = 0.0
private

Definition at line 169 of file BHE_1U.h.

Referenced by thermalResistances(), and updateHeatTransferCoefficients().

◆ inflow_outflow_bc_component_ids

std::pair<int, int> ProcessLib::HeatTransportBHE::BHE::BHE_1U::inflow_outflow_bc_component_ids[]
staticconstexpr
Initial value:
= {
{0, 1}}

Definition at line 148 of file BHE_1U.h.

148 {
149 {0, 1}};

◆ number_of_flow_legs

int ProcessLib::HeatTransportBHE::BHE::BHE_1U::number_of_flow_legs = 2
staticconstexpr

Definition at line 51 of file BHE_1U.h.

◆ number_of_grout_zones

int ProcessLib::HeatTransportBHE::BHE::BHE_1U::number_of_grout_zones = 2
staticconstexpr

Definition at line 50 of file BHE_1U.h.

Referenced by crossSectionAreas().

◆ number_of_unknowns

int ProcessLib::HeatTransportBHE::BHE::BHE_1U::number_of_unknowns = 4
staticconstexpr

Definition at line 49 of file BHE_1U.h.

The documentation for this class was generated from the following files: