ProcessLib::HT::anonymous_namespace{HTFEM.cpp} Namespace Reference

Functions
void	checkBiotStorageRelationValues (ParameterLib::SpatialPosition const &pos, double const biot_coefficient, double const specific_storage)

Function Documentation

checkBiotStorageRelationValues()

void ProcessLib::HT::anonymous_namespace{HTFEM.cpp}::checkBiotStorageRelationValues	(	ParameterLib::SpatialPosition const &	pos,
		double const	biot_coefficient,
		double const	specific_storage )

Checks the requirement \(\alpha_B=1 \Rightarrow S_s=0\) on already evaluated values and ends the run with OGS_FATAL if it is violated.

A specific storage that is not a number is not a violation. Where the check runs before the primary variables are known, a storage property depending on them evaluates to NaN, and NaN compares unequal to zero; such a setting is left to the checks during assembly, where the storage has a value. The same reasoning covers a Biot coefficient of NaN, which compares unequal to one.

Definition at line 25 of file HTFEM.cpp.

{
    if (biot_coefficient != 1.0 || specific_storage == 0.0 ||
        std::isnan(specific_storage))
    {
        return;
    }
 
    OGS_FATAL(
        "At {} the Biot coefficient evaluates to 1.0, which requires the "
        "specific storage of the solid phase to be 0.0, but it evaluates to "
        "{:g}.",
        pos, specific_storage);
}

References checkBiotStorageRelationValues(), and OGS_FATAL.

Referenced by checkBiotStorageRelationValues().