-
Notifications
You must be signed in to change notification settings - Fork 239
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
New MPL property: Brooks&Corey relative permeability #2656
New MPL property: Brooks&Corey relative permeability #2656
Conversation
bfafbc7
to
e7ad48a
Compare
e7ad48a
to
27dafff
Compare
auto const exponent = | ||
//! \ogs_file_param{prj__media__medium__properties__property__RelPermBrooksCorey__lambda} | ||
config.getConfigParameter<double>("lambda"); | ||
if (exponent == 0.) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The error message below is different from the check.
if (exponent == 0.) | |
if (exponent <= 0.) |
const double min_relative_permeability_liquid, | ||
const double min_relative_permeability_gas, | ||
const double exponent) | ||
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The empty line not necessary.
ParameterLib::SpatialPosition const& pos, double const t) const | ||
{ | ||
(void)primary_variable; | ||
assert((primary_variable == Variable::liquid_saturation) && |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Maybe it is better to check this also in the release version?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
not sure here, let's measure it!
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Talked with @Scinopode : The primary variable isn't user input. It is 'hard' coded in the process and mistakes should be detected while the development of a new process.
const double /*_min_relative_permeability_gas*/, | ||
const double /*exponent*/ | ||
); | ||
/// This method assigns a pointer to the meterial object that is the owner |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
/// This method assigns a pointer to the meterial object that is the owner | |
/// This method assigns a pointer to the material object that is the owner |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks ok.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
👍
@Scinopode Please add this and the saturation to the Changelog!! |
OpenGeoSys development has been moved to GitLab. |
Brooks&Corey relative permeability property for the generalized material property library.