ThermalResistanceHelpers.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 <cmath>
#include <numbers>
 
#include "BaseLib/Error.h"
#include "Pipe.h"
 
namespace ProcessLib
{
namespace HeatTransportBHE
{
namespace BHE
{
inline double pipeWallThermalResistance(Pipe const& pipe,
                                        double const wall_thermal_conductivity)
{
    return std::log(pipe.outsideDiameter() / pipe.diameter) /
           (2.0 * std::numbers::pi * wall_thermal_conductivity);
}

inline double computeRgs(double const chi, double const R_g)
{
    return (1 - chi) * R_g;
}

inline double computeRgg(double const chi, double const R_gs, double const R_ar,
                         double const R_g)
{
    double const R_gg = 2.0 * R_gs * (R_ar - 2.0 * chi * R_g) /
                        (2.0 * R_gs - R_ar + 2.0 * chi * R_g);
    if (!std::isfinite(R_gg))
    {
        OGS_FATAL(
            "Error!!! Grout Thermal Resistance is an infinite number! The "
            "simulation will be stopped!");
    }
    return R_gg;
}
}  // namespace BHE
}  // namespace HeatTransportBHE
}  // namespace ProcessLib