BHE_1P.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 <vector>
 
#include "BHECommon.h"
#include "BaseLib/Error.h"
#include "FlowAndTemperatureControl.h"
#include "PipeConfiguration1PType.h"
 
namespace ProcessLib
{
namespace HeatTransportBHE
{
namespace BHE
{
class BHE_1P final : public BHECommon
{
public:
    BHE_1P(BoreholeGeometry const& borehole,
           RefrigerantProperties const& refrigerant,
           GroutParameters const& grout,
           FlowAndTemperatureControl const& flowAndTemperatureControl,
           PipeConfiguration1PType const& pipes,
           bool const use_python_bcs);

    BHE_1P withGeometry(BoreholeGeometry const& g) const
    {
        return {g,     refrigerant,   grout, flowAndTemperatureControl,
                _pipe, use_python_bcs};
    }
 
    static constexpr int number_of_unknowns = 2;
    static constexpr int number_of_grout_zones = 1;
 
    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& elem_direction,
        int const section_index = 0) const;
 
    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)
    {
        // Here we are looping over two resistance terms
        // First PHI_fg is the resistance between pipe and grout
        // Second PHI_gs is the resistance between grout and soil
        switch (idx_bhe_unknowns)
        {
            case 0:  // PHI_fg
                R_matrix.block(0, NPoints, NPoints, NPoints) +=
                    -1.0 * matBHE_loc_R;
                R_matrix.block(NPoints, 0, NPoints, NPoints) +=
                    -1.0 * matBHE_loc_R;
 
                R_matrix.block(0, 0, NPoints, NPoints) +=
                    matBHE_loc_R;  // K_i/o
                R_matrix.block(NPoints, NPoints, NPoints, NPoints) +=
                    matBHE_loc_R;  // K_fg
                return;
            case 1:  // PHI_gs
                R_s_matrix += matBHE_loc_R;
 
                R_pi_s_matrix.block(NPoints, 0, NPoints, NPoints) +=
                    -1.0 * matBHE_loc_R;
 
                R_matrix.block(NPoints, NPoints, NPoints, NPoints) +=
                    matBHE_loc_R;  // K_fg
                return;
            default:
                OGS_FATAL(
                    "BHE_1P::assembleRMatrices: unknown index {:d} "
                    "out of range.",
                    idx_bhe_unknowns);
        }
    }

    double updateFlowRateAndTemperature(double T_out, double current_time);
 
    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*/);
 
public:
    std::array<double, number_of_unknowns> crossSectionAreas(
        int const section_index = 0) const;
 
    void updateHeatTransferCoefficients(double const flow_rate);
 
protected:
    PipeConfiguration1PType const _pipe;
 
private:
    std::vector<double> calcThermalResistances(
        double const Nu, int const section_index = 0) const;
};
}  // namespace BHE
}  // namespace HeatTransportBHE
}  // namespace ProcessLib