BHECommonCoaxial.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <Eigen/Core>
#include <optional>
 
#include "BHECommon.h"
#include "FlowAndTemperatureControl.h"
#include "PipeConfigurationCoaxial.h"
 
namespace ProcessLib
{
namespace HeatTransportBHE
{
namespace BHE
{
class BHECommonCoaxial : public BHECommon
{
public:
    BHECommonCoaxial(BoreholeGeometry const& borehole,
                     RefrigerantProperties const& refrigerant,
                     GroutParameters const& grout,
                     FlowAndTemperatureControl const& flowAndTemperatureControl,
                     PipeConfigurationCoaxial const& pipes,
                     bool const use_python_bcs);
 
    static constexpr int number_of_unknowns = 3;
    static constexpr int number_of_grout_zones = 1;
 
    double thermalResistance(int const unknown_index) const
    {
        return _thermal_resistances[unknown_index];
    }
 
    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;
 
    static constexpr std::pair<int, int> inflow_outflow_bc_component_ids[] = {
        {0, 1}};
 
    static std::array<std::pair<std::size_t /*node_id*/, int /*component*/>, 2>
    getBHEInflowDirichletBCNodesAndComponents(
        std::size_t const top_node_id,
        std::size_t const /*bottom_node_id*/,
        int const in_component_id);
 
    static std::optional<
        std::array<std::pair<std::size_t /*node_id*/, int /*component*/>, 2>>
    getBHEBottomDirichletBCNodesAndComponents(std::size_t const bottom_node_id,
                                              int const in_component_id,
                                              int const out_component_id);
 
    std::array<double, number_of_unknowns> pipeHeatConductions() const;
 
    std::array<Eigen::Vector3d, number_of_unknowns> pipeAdvectionVectors(
        Eigen::Vector3d const& /*elem_direction*/) const;
 
    double cross_section_area_inner_pipe, cross_section_area_annulus,
        cross_section_area_grout;
 
    void updateHeatTransferCoefficients(double const flow_rate);
 
protected:
    PipeConfigurationCoaxial const _pipes;
 
    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;

    std::array<double, number_of_unknowns> _thermal_resistances;

    double _flow_velocity_inner = std::numeric_limits<double>::quiet_NaN(),
           _flow_velocity_annulus = std::numeric_limits<double>::quiet_NaN();
};
}  // namespace BHE
}  // namespace HeatTransportBHE
}  // end of namespace ProcessLib