FluidThermalExpansion.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "FluidThermalExpansion.h"
 
#include "MaterialLib/MPL/Utils/GetLiquidThermalExpansivity.h"
 
namespace ProcessLib::ThermoRichardsMechanics
{
template <int DisplacementDim>
void FluidThermalExpansionModel<DisplacementDim>::eval(
    SpaceTimeData const& x_t, MediaData const& media_data,
    CapillaryPressureData<DisplacementDim> const& p_cap_data,
    TemperatureData<DisplacementDim> const& T_data,
    SolidThermalExpansionData<DisplacementDim> const& s_therm_exp_data,
    PorosityData const& poro_data, LiquidDensityData const& rho_L_data,
    BiotData const& biot_data, FluidThermalExpansionData& out) const
{
    static constexpr int kelvin_vector_size =
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
    using Invariants = MathLib::KelvinVector::Invariants<kelvin_vector_size>;
 
    namespace MPL = MaterialPropertyLib;
    MPL::VariableArray variables;
    variables.liquid_phase_pressure = -p_cap_data.p_cap;
    variables.temperature = T_data.T;
 
    double const phi = poro_data.phi;
 
    double const fluid_volumetric_thermal_expansion =
        phi * MPL::getLiquidThermalExpansivity(media_data.liquid, variables,
                                               rho_L_data.rho_LR, x_t.x, x_t.t,
                                               x_t.dt);
 
    out.eff_thermal_expansion =
        (biot_data() - phi) *
            Invariants::trace(
                s_therm_exp_data.solid_linear_thermal_expansivity_vector) +
        fluid_volumetric_thermal_expansion;
}
 
template struct FluidThermalExpansionModel<2>;
template struct FluidThermalExpansionModel<3>;
}  // namespace ProcessLib::ThermoRichardsMechanics