Skip to content
Snippets Groups Projects
Commit a53bcb11 authored by wenqing's avatar wenqing
Browse files

[Unsat] Check m<0 for van Genuchten relative permeability model

parent 02610383
No related branches found
No related tags found
No related merge requests found
...@@ -52,7 +52,7 @@ std::unique_ptr<RelativePermeability> createWettingPhaseVanGenuchten( ...@@ -52,7 +52,7 @@ std::unique_ptr<RelativePermeability> createWettingPhaseVanGenuchten(
//! \ogs_file_param{material_property__porous_medium__porous_medium__relative_permeability__WettingPhaseVanGenuchten__m} //! \ogs_file_param{material_property__porous_medium__porous_medium__relative_permeability__WettingPhaseVanGenuchten__m}
const double m = config.getConfigParameter<double>("m"); const double m = config.getConfigParameter<double>("m");
if (m > 1.0) // m <= 1 if (m < 0. || m > 1.0) // m <= 1
{ {
OGS_FATAL( OGS_FATAL(
"The exponent parameter of WettingPhaseVanGenuchten relative\n" "The exponent parameter of WettingPhaseVanGenuchten relative\n"
...@@ -84,7 +84,7 @@ std::unique_ptr<RelativePermeability> createNonWettingPhaseVanGenuchten( ...@@ -84,7 +84,7 @@ std::unique_ptr<RelativePermeability> createNonWettingPhaseVanGenuchten(
//! \ogs_file_param{material_property__porous_medium__porous_medium__relative_permeability__NonWettingPhaseVanGenuchten__m} //! \ogs_file_param{material_property__porous_medium__porous_medium__relative_permeability__NonWettingPhaseVanGenuchten__m}
const double m = config.getConfigParameter<double>("m"); const double m = config.getConfigParameter<double>("m");
if (m > 1.0) // m <= 1 if (m < 0.|| m > 1.0) // m <= 1
{ {
OGS_FATAL( OGS_FATAL(
"The exponent parameter of NonWettingPhaseVanGenuchten relative\n" "The exponent parameter of NonWettingPhaseVanGenuchten relative\n"
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment