Skip to content
Snippets Groups Projects
Commit 8a0cf88a authored by Dmitri Naumov's avatar Dmitri Naumov
Browse files

Merge branch 'deviated_dbhe' into 'master'

[BHE/coaxial] Calculate advection vectors based on the element direction for deviated coaxial BHE.

See merge request ogs/ogs!3360
parents 4ce034c2 6c54debe
No related branches found
No related tags found
No related merge requests found
...@@ -88,18 +88,18 @@ BHECommonCoaxial::pipeHeatConductions() const ...@@ -88,18 +88,18 @@ BHECommonCoaxial::pipeHeatConductions() const
std::array<Eigen::Vector3d, BHECommonCoaxial::number_of_unknowns> std::array<Eigen::Vector3d, BHECommonCoaxial::number_of_unknowns>
BHECommonCoaxial::pipeAdvectionVectors( BHECommonCoaxial::pipeAdvectionVectors(
Eigen::Vector3d const& /*elem_direction*/) const Eigen::Vector3d const& elem_direction) const
{ {
double const rho_r = refrigerant.density; double const rho_r = refrigerant.density;
double const Cp_r = refrigerant.specific_heat_capacity; double const Cp_r = refrigerant.specific_heat_capacity;
auto v = velocities(); auto const v = velocities();
return {{// pipe i, Eq. 26 and Eq. 23 return {// pipe i, Eq. 26 and Eq. 23
{0, 0, -rho_r * Cp_r * std::abs(v[0])}, rho_r * Cp_r * std::abs(v[0]) * elem_direction,
// pipe o, Eq. 27 and Eq. 24 // pipe o, Eq. 27 and Eq. 24
{0, 0, rho_r * Cp_r * std::abs(v[1])}, -rho_r * Cp_r * std::abs(v[1]) * elem_direction,
// grout g, Eq. 28 and Eq. 25 // grout g, Eq. 28 and Eq. 25
{0, 0, 0}}}; {0, 0, 0}};
} }
std::array<double, BHECommonCoaxial::number_of_unknowns> std::array<double, BHECommonCoaxial::number_of_unknowns>
......
...@@ -74,6 +74,8 @@ Currently there are 4 types of BHE available. Following the convention in Diersc ...@@ -74,6 +74,8 @@ Currently there are 4 types of BHE available. Following the convention in Diersc
* CXA: coaxial pipe with annular space as the inlet downwards flow and the centre part as outlet upwards flow; * CXA: coaxial pipe with annular space as the inlet downwards flow and the centre part as outlet upwards flow;
* CXC: coaxial pipe with a reversed flow direction to CXA type. * CXC: coaxial pipe with a reversed flow direction to CXA type.
Especially in CXA and CXC type, the direction of the borehole itself could be deviated by any angle, which is defined by mesh. The inflow direction will be in accordance with the direction of the line element (represents the BHE borehole) in the mesh. And the outlet direction is the opposite of the inflow direction.
The cross-sections of these 4 types of BHEs are illustrated in the following figures. The cross-sections of these 4 types of BHEs are illustrated in the following figures.
{{< img src="../u_type.png" width="50">}} {{< img src="../u_type.png" width="50">}}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment