Detailed Description

Definition at line 20 of file BHECommonCoaxial.h.

#include <BHECommonCoaxial.h>

Inheritance diagram for ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial:

Collaboration diagram for ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial:

Public Member Functions
	BHECommonCoaxial (BoreholeGeometry const &borehole, RefrigerantProperties const &refrigerant, GroutParameters const &grout, FlowAndTemperatureControl const &flowAndTemperatureControl, PipeConfigurationCoaxial const &pipes, bool const use_python_bcs)
double	updateFlowRateAndTemperature (double T_out, double current_time)
std::vector< double >	calcThermalResistances (double const Nu_inner_pipe, double const Nu_annulus_pipe, int const section_index=0) const
std::array< double, number_of_unknowns >	pipeHeatCapacities () const
std::array< double, number_of_unknowns >	pipeHeatConductions (int const section_index=0) const
std::array< Eigen::Vector3d, number_of_unknowns >	pipeAdvectionVectors (Eigen::Vector3d const &, int const section_index=0) const
void	updateHeatTransferCoefficients (double const flow_rate)
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
int	getNumberOfSections () const
	Get number of sections in the borehole geometry.
std::vector< double > const &	getSectionBoundaries () const
	Get section boundaries (cumulative distance from wellhead).
double	thermalResistanceAtSection (int const unknown_index, int const section_index) const
	Get thermal resistance for a specific section and unknown.

Static Public Member Functions
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 = 3
static constexpr int	number_of_grout_zones = 1
static constexpr std::pair< int, int >	inflow_outflow_bc_component_ids []

Protected Member Functions
std::array< double, 2 >	velocities () const
virtual std::vector< double >	getThermalResistances (double const &R_gs, double const &R_ff, double const &R_fg) const =0
Protected Member Functions inherited from ProcessLib::HeatTransportBHE::BHE::BHECommon
template<typename Fn>
void	recomputeSectionalResistances (Fn &&calcForSection)
	Recompute _sectional_thermal_resistances for all sections.
double	getClampedFlowVelocity (int const section_index) const
	Get velocity for a section, clamping to last section if index exceeds the number of velocity entries. Returns 0 when empty.
std::vector< double > const &	getThermalResistancesAtSection (int const section_index) const
	Get the thermal resistance vector for a given section.

Protected Attributes
PipeConfigurationCoaxial const	_pipes
double	cross_section_area_inner_pipe
double	cross_section_area_annulus
Protected Attributes inherited from ProcessLib::HeatTransportBHE::BHE::BHECommon
std::vector< double >	_flow_velocities
std::vector< std::vector< double > >	_sectional_thermal_resistances

Private Member Functions
virtual void	assignVelocities (double inner_vel, double annulus_vel)=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

Constructor & Destructor Documentation

◆ BHECommonCoaxial()

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

Definition at line 15 of file BHECommonCoaxial.cpp.

    : BHECommon{borehole, refrigerant, grout, flowAndTemperatureControl,
                use_python_bcs},
      _pipes(pipes)
{
    cross_section_area_inner_pipe = _pipes.inner_pipe.area();
    cross_section_area_annulus =
        _pipes.outer_pipe.area() - _pipes.inner_pipe.outsideArea();
}

References ProcessLib::HeatTransportBHE::BHE::BHECommon::BHECommon(), _pipes, cross_section_area_annulus, cross_section_area_inner_pipe, ProcessLib::HeatTransportBHE::BHE::BHECommon::flowAndTemperatureControl, ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant, and ProcessLib::HeatTransportBHE::BHE::BHECommon::use_python_bcs.

Referenced by ProcessLib::HeatTransportBHE::BHE::BHE_CXA::BHE_CXA(), and ProcessLib::HeatTransportBHE::BHE::BHE_CXC::BHE_CXC().

Member Function Documentation

◆ assignVelocities()

virtual void ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::assignVelocities	(	double	inner_vel,
		double	annulus_vel )

privatepure virtual

Assigns _flow_velocities from the physical channel velocities. Subclasses encode which channel maps to which unknown.

Implemented in ProcessLib::HeatTransportBHE::BHE::BHE_CXA, and ProcessLib::HeatTransportBHE::BHE::BHE_CXC.

Referenced by updateHeatTransferCoefficients().

◆ calcThermalResistances()

std::vector< double > ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::calcThermalResistances	(	double const	Nu_inner_pipe,
		double const	Nu_annulus_pipe,
		int const	section_index = 0 ) const

Definition at line 95 of file BHECommonCoaxial.cpp.

{
    // thermal resistances due to advective flow of refrigerant in the pipes
    auto const R_advective = calculateAdvectiveThermalResistance(
        _pipes.inner_pipe, _pipes.outer_pipe, refrigerant, Nu_inner_pipe,
        Nu_annulus_pipe);
 
    // thermal resistance due to thermal conductivity of the pipe wall material
    auto const R_conductive = calculatePipeWallThermalResistance(
        _pipes.inner_pipe, _pipes.outer_pipe);
 
    // thermal resistance due to the grout transition and grout-soil exchange.
    auto const R = calculateGroutAndGroutSoilExchangeThermalResistance(
        _pipes.outer_pipe, grout,
        borehole_geometry.sections.diameterAtSection(section_index));
 
    // thermal resistance due to grout-soil exchange
    double const R_gs = R.grout_soil;
 
    // Eq. 56 and 57
    double const R_ff = R_advective.inner_pipe_coaxial + R_advective.a_annulus +
                        R_conductive.inner_pipe_coaxial;
    double const R_fg =
        R_advective.b_annulus + R_conductive.annulus + R.conductive_b;
 
    return getThermalResistances(R_gs, R_ff, R_fg);
}

References _pipes, ProcessLib::HeatTransportBHE::BHE::BHECommon::borehole_geometry, ProcessLib::HeatTransportBHE::BHE::calculateAdvectiveThermalResistance(), ProcessLib::HeatTransportBHE::BHE::calculateGroutAndGroutSoilExchangeThermalResistance(), ProcessLib::HeatTransportBHE::BHE::calculatePipeWallThermalResistance(), getThermalResistances(), ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, and ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant.

Referenced by updateHeatTransferCoefficients().

◆ getBHEBottomDirichletBCNodesAndComponents()

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

static

Definition at line 137 of file BHECommonCoaxial.cpp.

{
    return {{std::make_pair(bottom_node_id, in_component_id),
             std::make_pair(bottom_node_id, out_component_id)}};
}

◆ getBHEInflowDirichletBCNodesAndComponents()

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

static

Definition at line 126 of file BHECommonCoaxial.cpp.

{
    return {std::make_pair(top_node_id, in_component_id),
            std::make_pair(top_node_id, in_component_id + 1)};
}

◆ getThermalResistances()

virtual std::vector< double > ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::getThermalResistances	(	double const &	R_gs,
		double const &	R_ff,
		double const &	R_fg ) const

protectedpure virtual

Implemented in ProcessLib::HeatTransportBHE::BHE::BHE_CXA, and ProcessLib::HeatTransportBHE::BHE::BHE_CXC.

Referenced by calcThermalResistances().

◆ pipeAdvectionVectors()

std::array< Eigen::Vector3d, BHECommonCoaxial::number_of_unknowns > ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::pipeAdvectionVectors	(	Eigen::Vector3d const &	elem_direction,
		int const	section_index = 0 ) const

Definition at line 79 of file BHECommonCoaxial.cpp.

{
    // Pipe dimensions are constant; velocity does not vary by section.
    double const rho_r = refrigerant.density;
    double const Cp_r = refrigerant.specific_heat_capacity;
    auto const v = velocities();
 
    return {// pipe i, Eq. 26 and Eq. 23
            rho_r * Cp_r * std::abs(v[0]) * elem_direction,
            // pipe o, Eq. 27 and Eq. 24
            -rho_r * Cp_r * std::abs(v[1]) * elem_direction,
            // grout g, Eq. 28 and Eq. 25
            {0, 0, 0}};
}

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

◆ pipeHeatCapacities()

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

Definition at line 32 of file BHECommonCoaxial.cpp.

{
    double const rho_r = refrigerant.density;
    double const specific_heat_capacity = refrigerant.specific_heat_capacity;
    double const rho_g = grout.rho_g;
    double const porosity_g = grout.porosity_g;
    double const heat_cap_g = grout.heat_cap_g;
 
    return {{/*i*/ rho_r * specific_heat_capacity,
             /*o*/ rho_r * specific_heat_capacity,
             /*g*/ (1.0 - porosity_g) * rho_g * heat_cap_g}};
}

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

◆ pipeHeatConductions()

std::array< double, BHECommonCoaxial::number_of_unknowns > ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::pipeHeatConductions ( int const section_index = 0 ) const

Definition at line 55 of file BHECommonCoaxial.cpp.

{
    // Pipe dimensions are constant; velocity does not vary by section.
    double const lambda_r = refrigerant.thermal_conductivity;
    double const rho_r = refrigerant.density;
    double const Cp_r = refrigerant.specific_heat_capacity;
    double const alpha_L = _pipes.longitudinal_dispersion_length;
    double const porosity_g = grout.porosity_g;
    double const lambda_g = grout.lambda_g;
 
    auto v = velocities();
    // Here we calculate the laplace coefficients in the governing
    // equations of BHE. These governing equations can be found in
    // 1) Diersch (2013) FEFLOW book on page 952, M.120-122, or
    // 2) Diersch (2011) Comp & Geosci 37:1122-1135, Eq. 26-28, 23-25
    return {{// pipe i, Eq. 26 and Eq. 23
             (lambda_r + rho_r * Cp_r * alpha_L * std::abs(v[0])),
             // pipe o, Eq. 27 and Eq. 24
             (lambda_r + rho_r * Cp_r * alpha_L * std::abs(v[1])),
             // pipe g, Eq. 28 and Eq. 25
             (1.0 - porosity_g) * lambda_g}};
}

References _pipes, ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant, and velocities().

◆ updateFlowRateAndTemperature()

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

Definition at line 45 of file BHECommonCoaxial.cpp.

{
    auto values =
        visit([&](auto const& control) { return control(T_out, current_time); },
              flowAndTemperatureControl);
    updateHeatTransferCoefficients(values.flow_rate);
    return values.temperature;
}

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

◆ updateHeatTransferCoefficients()

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

Definition at line 145 of file BHECommonCoaxial.cpp.

{
    auto const tm_flow_inner = calculateThermoMechanicalFlowPropertiesPipe(
        _pipes.inner_pipe, borehole_geometry.length, refrigerant, flow_rate);
 
    auto const tm_flow_annulus = calculateThermoMechanicalFlowPropertiesAnnulus(
        _pipes.inner_pipe, _pipes.outer_pipe, borehole_geometry.length,
        refrigerant, flow_rate);
 
    assignVelocities(tm_flow_inner.velocity, tm_flow_annulus.velocity);
 
    recomputeSectionalResistances(
        [&](int i)
        {
            return calcThermalResistances(tm_flow_inner.nusselt_number,
                                          tm_flow_annulus.nusselt_number, i);
        });
}

References _pipes, assignVelocities(), ProcessLib::HeatTransportBHE::BHE::BHECommon::borehole_geometry, calcThermalResistances(), ProcessLib::HeatTransportBHE::BHE::calculateThermoMechanicalFlowPropertiesAnnulus(), ProcessLib::HeatTransportBHE::BHE::calculateThermoMechanicalFlowPropertiesPipe(), ProcessLib::HeatTransportBHE::BHE::BHECommon::recomputeSectionalResistances(), and ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant.

Referenced by ProcessLib::HeatTransportBHE::BHE::BHE_CXA::BHE_CXA(), ProcessLib::HeatTransportBHE::BHE::BHE_CXC::BHE_CXC(), and updateFlowRateAndTemperature().

◆ velocities()

std::array< double, 2 > ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::velocities ( ) const

inlineprotected

Returns fluid velocities indexed by unknown: {v_for_unknown_0, v_for_unknown_1}. The mapping from physical channels (inner pipe / annulus) to unknowns differs between CXA and CXC and is set once per flow-rate update via assignVelocities().

Definition at line 75 of file BHECommonCoaxial.h.

    {
        return {_flow_velocities[0], _flow_velocities[1]};
    }

References ProcessLib::HeatTransportBHE::BHE::BHECommon::_flow_velocities.

Referenced by pipeAdvectionVectors(), and pipeHeatConductions().

Member Data Documentation

◆ _pipes

PipeConfigurationCoaxial const ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::_pipes

protected

Definition at line 67 of file BHECommonCoaxial.h.

Referenced by BHECommonCoaxial(), calcThermalResistances(), ProcessLib::HeatTransportBHE::BHE::BHE_CXA::crossSectionAreas(), ProcessLib::HeatTransportBHE::BHE::BHE_CXC::crossSectionAreas(), pipeHeatConductions(), updateHeatTransferCoefficients(), ProcessLib::HeatTransportBHE::BHE::BHE_CXA::withGeometry(), and ProcessLib::HeatTransportBHE::BHE::BHE_CXC::withGeometry().

◆ cross_section_area_annulus

double ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::cross_section_area_annulus

protected

Definition at line 69 of file BHECommonCoaxial.h.

Referenced by BHECommonCoaxial(), ProcessLib::HeatTransportBHE::BHE::BHE_CXA::crossSectionAreas(), and ProcessLib::HeatTransportBHE::BHE::BHE_CXC::crossSectionAreas().

◆ cross_section_area_inner_pipe

double ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::cross_section_area_inner_pipe

protected

Definition at line 69 of file BHECommonCoaxial.h.

Referenced by BHECommonCoaxial(), ProcessLib::HeatTransportBHE::BHE::BHE_CXA::crossSectionAreas(), and ProcessLib::HeatTransportBHE::BHE::BHE_CXC::crossSectionAreas().

◆ inflow_outflow_bc_component_ids

std::pair<int, int> ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::inflow_outflow_bc_component_ids[]

staticconstexpr

Initial value:

= {

{0, 1}}

Definition at line 42 of file BHECommonCoaxial.h.

42 {

43 {0, 1}};

◆ number_of_grout_zones

int ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::number_of_grout_zones = 1

staticconstexpr

Definition at line 31 of file BHECommonCoaxial.h.

◆ number_of_unknowns

int ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::number_of_unknowns = 3

staticconstexpr

Definition at line 30 of file BHECommonCoaxial.h.

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

Detailed Description

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members

Constructor & Destructor Documentation

◆ BHECommonCoaxial()

Member Function Documentation

◆ assignVelocities()

◆ calcThermalResistances()

◆ getBHEBottomDirichletBCNodesAndComponents()

◆ getBHEInflowDirichletBCNodesAndComponents()

◆ getThermalResistances()

◆ pipeAdvectionVectors()

◆ pipeHeatCapacities()

◆ pipeHeatConductions()

◆ updateFlowRateAndTemperature()

◆ updateHeatTransferCoefficients()

◆ velocities()

Member Data Documentation

◆ _pipes

◆ cross_section_area_annulus

◆ cross_section_area_inner_pipe

◆ inflow_outflow_bc_component_ids

◆ number_of_grout_zones

◆ number_of_unknowns