Detailed Description

Definition at line 25 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	thermalResistance (int const unknown_index) const

double	updateFlowRateAndTemperature (double T_out, double current_time)

std::array< double, number_of_unknowns >	calcThermalResistances (double const Nu_inner_pipe, double const Nu_annulus_pipe) 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 &) const

void	updateHeatTransferCoefficients (double const flow_rate)

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)

Public Attributes
double	cross_section_area_inner_pipe

double	cross_section_area_annulus

double	cross_section_area_grout

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

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
virtual std::array< double, 2 >	velocities () const =0

virtual std::array< double, number_of_unknowns >	getThermalResistances (double const &R_gs, double const &R_ff, double const &R_fg) const =0

Protected Attributes
PipeConfigurationCoaxial const	_pipes

std::array< double, number_of_unknowns >	_thermal_resistances

double	_flow_velocity_inner = std::numeric_limits<double>::quiet_NaN()
	Flow velocity inside the pipes and annulus. Depends on the flow_rate.

double	_flow_velocity_annulus = std::numeric_limits<double>::quiet_NaN()

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 23 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();
    cross_section_area_grout =
        borehole_geometry.area() - _pipes.outer_pipe.outsideArea();
 
    _thermal_resistances.fill(std::numeric_limits<double>::quiet_NaN());
}

References _pipes, _thermal_resistances, ProcessLib::HeatTransportBHE::BHE::BoreholeGeometry::area(), ProcessLib::HeatTransportBHE::BHE::Pipe::area(), ProcessLib::HeatTransportBHE::BHE::BHECommon::borehole_geometry, cross_section_area_annulus, cross_section_area_grout, cross_section_area_inner_pipe, ProcessLib::HeatTransportBHE::BHE::PipeConfigurationCoaxial::inner_pipe, ProcessLib::HeatTransportBHE::BHE::PipeConfigurationCoaxial::outer_pipe, and ProcessLib::HeatTransportBHE::BHE::Pipe::outsideArea().

Member Function Documentation

◆ calcThermalResistances()

std::array< double, BHECommonCoaxial::number_of_unknowns > ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::calcThermalResistances	(	double const	Nu_inner_pipe,
		double const	Nu_annulus_pipe ) const

Definition at line 106 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.diameter);
 
    // 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);
}

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 146 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 135 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::array< double, number_of_unknowns > 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 ) const

Definition at line 90 of file BHECommonCoaxial.cpp.

{
    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::RefrigerantProperties::density, ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant, ProcessLib::HeatTransportBHE::BHE::RefrigerantProperties::specific_heat_capacity, and velocities().

◆ pipeHeatCapacities()

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

Definition at line 44 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::RefrigerantProperties::density, ProcessLib::HeatTransportBHE::BHE::BHECommon::grout, ProcessLib::HeatTransportBHE::BHE::GroutParameters::heat_cap_g, ProcessLib::HeatTransportBHE::BHE::GroutParameters::porosity_g, ProcessLib::HeatTransportBHE::BHE::BHECommon::refrigerant, ProcessLib::HeatTransportBHE::BHE::GroutParameters::rho_g, and ProcessLib::HeatTransportBHE::BHE::RefrigerantProperties::specific_heat_capacity.

◆ pipeHeatConductions()

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

Definition at line 67 of file BHECommonCoaxial.cpp.

{
    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}};
}

◆ thermalResistance()

double ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::thermalResistance ( int const unknown_index ) const

inline

Definition at line 38 of file BHECommonCoaxial.h.

    {
        return _thermal_resistances[unknown_index];
    }

References _thermal_resistances.

◆ updateFlowRateAndTemperature()

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

Definition at line 57 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 154 of file BHECommonCoaxial.cpp.

{
    auto const tm_flow_properties_annulus =
        calculateThermoMechanicalFlowPropertiesAnnulus(_pipes.inner_pipe,
                                                       _pipes.outer_pipe,
                                                       borehole_geometry.length,
                                                       refrigerant,
                                                       flow_rate);
 
    _flow_velocity_annulus = tm_flow_properties_annulus.velocity;
 
    auto const tm_flow_properties = calculateThermoMechanicalFlowPropertiesPipe(
        _pipes.inner_pipe, borehole_geometry.length, refrigerant, flow_rate);
 
    _flow_velocity_inner = tm_flow_properties.velocity;
 
    _thermal_resistances =
        calcThermalResistances(tm_flow_properties.nusselt_number,
                               tm_flow_properties_annulus.nusselt_number);
}

References _flow_velocity_annulus, _flow_velocity_inner, _pipes, _thermal_resistances, ProcessLib::HeatTransportBHE::BHE::BHECommon::borehole_geometry, calcThermalResistances(), ProcessLib::HeatTransportBHE::BHE::calculateThermoMechanicalFlowPropertiesAnnulus(), ProcessLib::HeatTransportBHE::BHE::calculateThermoMechanicalFlowPropertiesPipe(), ProcessLib::HeatTransportBHE::BHE::PipeConfigurationCoaxial::inner_pipe, ProcessLib::HeatTransportBHE::BHE::BoreholeGeometry::length, ProcessLib::HeatTransportBHE::BHE::PipeConfigurationCoaxial::outer_pipe, 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()

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

protectedpure virtual

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

Referenced by pipeAdvectionVectors(), and pipeHeatConductions().

Member Data Documentation

◆ _flow_velocity_annulus

double ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::_flow_velocity_annulus = std::numeric_limits<double>::quiet_NaN()

protected

Definition at line 92 of file BHECommonCoaxial.h.

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

◆ _flow_velocity_inner

double ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::_flow_velocity_inner = std::numeric_limits<double>::quiet_NaN()

protected

Flow velocity inside the pipes and annulus. Depends on the flow_rate.

Definition at line 91 of file BHECommonCoaxial.h.

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

◆ _pipes

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

protected

Definition at line 76 of file BHECommonCoaxial.h.

Referenced by BHECommonCoaxial(), calcThermalResistances(), pipeHeatConductions(), and updateHeatTransferCoefficients().

◆ _thermal_resistances

std::array<double, number_of_unknowns> ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::_thermal_resistances

protected

Here we store the thermal resistances needed for computation of the heat exchange coefficients in the governing equations of BHE. These governing equations can be found in 1) Diersch (2013) FEFLOW book on page 958, M.3, or 2) Diersch (2011) Comp & Geosci 37:1122-1135, Eq. 90-97.

Definition at line 88 of file BHECommonCoaxial.h.

Referenced by BHECommonCoaxial(), thermalResistance(), and updateHeatTransferCoefficients().

◆ cross_section_area_annulus

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

Definition at line 70 of file BHECommonCoaxial.h.

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

◆ cross_section_area_grout

double ProcessLib::HeatTransportBHE::BHE::BHECommonCoaxial::cross_section_area_grout

Definition at line 71 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

Definition at line 70 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

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

staticconstexpr

Initial value:

= {

{0, 1}}

Definition at line 50 of file BHECommonCoaxial.h.

50 {

51 {0, 1}};

◆ number_of_grout_zones

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

staticconstexpr

Definition at line 36 of file BHECommonCoaxial.h.

◆ number_of_unknowns

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

staticconstexpr

Definition at line 35 of file BHECommonCoaxial.h.

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

ProcessLib/HeatTransportBHE/BHE/BHECommonCoaxial.h
ProcessLib/HeatTransportBHE/BHE/BHECommonCoaxial.cpp

Detailed Description

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Protected Attributes

Constructor & Destructor Documentation

◆ BHECommonCoaxial()

Member Function Documentation

◆ calcThermalResistances()

◆ getBHEBottomDirichletBCNodesAndComponents()

◆ getBHEInflowDirichletBCNodesAndComponents()

◆ getThermalResistances()

◆ pipeAdvectionVectors()

◆ pipeHeatCapacities()

◆ pipeHeatConductions()

◆ thermalResistance()

◆ updateFlowRateAndTemperature()

◆ updateHeatTransferCoefficients()

◆ velocities()

Member Data Documentation

◆ _flow_velocity_annulus

◆ _flow_velocity_inner

◆ _pipes

◆ _thermal_resistances

◆ cross_section_area_annulus

◆ cross_section_area_grout

◆ cross_section_area_inner_pipe

◆ inflow_outflow_bc_component_ids

◆ number_of_grout_zones

◆ number_of_unknowns