CreateBHEUType.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateBHEUType.h"
 
#include <cmath>
 
#include "BHE_1U.h"
#include "BHE_2U.h"
#include "BaseLib/ConfigTree.h"
#include "BaseLib/Error.h"
#include "CreateFlowAndTemperatureControl.h"
namespace ProcessLib
{
namespace HeatTransportBHE
{
namespace BHE
{
static std::tuple<BoreholeGeometry, RefrigerantProperties, GroutParameters,
                  FlowAndTemperatureControl, PipeConfigurationUType, bool>
parseBHEUTypeConfig(
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>>& parameters,
    std::map<std::string,
             std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
        curves)
{
    // if the BHE is using python boundary condition
    auto const bhe_if_use_python_bc_conf =
        config.getConfigParameter<bool>("use_bhe_pipe_network", false);
    DBUG("If using python boundary condition : {:s}",
         (bhe_if_use_python_bc_conf) ? "true" : "false");
 
    auto const borehole_geometry =
        createBoreholeGeometry(config.getConfigSubtree("borehole"), parameters);

    auto const& pipes_config = config.getConfigSubtree("pipes");
    Pipe const inlet_pipe =
        createPipe(pipes_config.getConfigSubtree("inlet"), parameters);
    Pipe const outlet_pipe =
        createPipe(pipes_config.getConfigSubtree("outlet"), parameters);
    const auto pipe_distance =
        pipes_config.getConfigParameter<double>("distance_between_pipes");
    const auto pipe_longitudinal_dispersion_length =
        pipes_config.getConfigParameter<double>(
            "longitudinal_dispersion_length");
    PipeConfigurationUType const pipes{inlet_pipe, outlet_pipe, pipe_distance,
                                       pipe_longitudinal_dispersion_length};
 
    if (pipe_distance <= 0)
    {
        OGS_FATAL(
            "distance_between_pipes must be positive for U-type BHEs, got "
            "{:g}.",
            pipe_distance);
    }
 
    // The resistance formulas use a single pipe outside diameter; the BHE_1U /
    // BHE_2U constructor enforces that inlet and outlet share it, so the inlet
    // value is representative here.
    double const d0 = pipes.inlet.outsideDiameter();
    if (pipe_distance <= d0)
    {
        OGS_FATAL(
            "distance_between_pipes ({:g}) must be > pipe outside diameter "
            "({:g}) for valid U-type thermal resistance formulas. A distance "
            "equal to the pipe outside diameter makes the inter-grout acosh "
            "argument exactly 1, i.e. a zero resistance and an infinite "
            "assembly coefficient.",
            pipe_distance, d0);
    }

    auto const grout = createGroutParameters(config.getConfigSubtree("grout"));
 
    auto const refrigerant =
        createRefrigerantProperties(config.getConfigSubtree("refrigerant"));
 
    auto const flowAndTemperatureControl = createFlowAndTemperatureControl(
        config.getConfigSubtree("flow_and_temperature_control"),
        parameters,
        curves,
        refrigerant);
 
    return {borehole_geometry,         refrigerant, grout,
            flowAndTemperatureControl, pipes,       bhe_if_use_python_bc_conf};
}
 
template <typename T_BHE>
T_BHE createBHEUType(
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>>& parameters,
    std::map<std::string,
             std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
        curves)
{
    auto UType = parseBHEUTypeConfig(config, parameters, curves);
    return {std::get<0>(UType), std::get<1>(UType), std::get<2>(UType),
            std::get<3>(UType), std::get<4>(UType), std::get<5>(UType)};
}
 
template BHE_1U createBHEUType<BHE_1U>(
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>>& parameters,
    std::map<std::string,
             std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
        curves);
 
template BHE_2U createBHEUType<BHE_2U>(
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>>& parameters,
    std::map<std::string,
             std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
        curves);
}  // namespace BHE
}  // namespace HeatTransportBHE
}  // namespace ProcessLib