Skip to content
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

Merged

Conversation

Scinopode
Copy link
Contributor

Brooks&Corey relative permeability property for the generalized material property library.

  1. Feature description was added to the changelog
  2. Tests covering your feature were added?
  3. Any new feature or behavior change was documented?

@Scinopode Scinopode force-pushed the mpl_brooks_corey_relative_permeability branch 2 times, most recently from bfafbc7 to e7ad48a Compare September 11, 2019 09:30
@Scinopode Scinopode force-pushed the mpl_brooks_corey_relative_permeability branch from e7ad48a to 27dafff Compare September 12, 2019 05:08
auto const exponent =
//! \ogs_file_param{prj__media__medium__properties__property__RelPermBrooksCorey__lambda}
config.getConfigParameter<double>("lambda");
if (exponent == 0.)
Copy link
Member

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.

Suggested change
if (exponent == 0.)
if (exponent <= 0.)

const double min_relative_permeability_liquid,
const double min_relative_permeability_gas,
const double exponent)

Copy link
Member

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) &&
Copy link
Member

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?

Copy link
Member

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!

Copy link
Member

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
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
/// 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

Copy link
Member

@endJunction endJunction left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks ok.

Copy link
Member

@TomFischer TomFischer left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

👍

@TomFischer TomFischer merged commit 5b00fe3 into ufz:master Sep 12, 2019
@endJunction
Copy link
Member

@Scinopode Please add this and the saturation to the Changelog!!

@Scinopode Scinopode deleted the mpl_brooks_corey_relative_permeability branch September 18, 2019 05:56
@ogsbot
Copy link
Member

ogsbot commented Jun 19, 2020

OpenGeoSys development has been moved to GitLab.

See this pull request on GitLab.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants