#include "tutorials/cor-0to1/src/MaterialUnsaturatedFlow.hpp"

Inheritance diagram for MixTransport::CommonMaterialData:

Collaboration diagram for MixTransport::CommonMaterialData:

Public Types
typedef boost::function< boost::shared_ptr< CommonMaterialData >(const CommonMaterialData &data)>	RegisterHook

Public Member Functions
	CommonMaterialData ()

double	initialPcEval () const
	Initialize head.

void	addOptions (po::options_description &o, const std::string &prefix)

void	printMatParameters (const int id, const std::string &prefix) const

Public Member Functions inherited from MixTransport::GenericMaterial
virtual	~GenericMaterial ()

virtual MoFEMErrorCode	calK ()

virtual MoFEMErrorCode	calDiffK ()

virtual MoFEMErrorCode	calC ()

virtual MoFEMErrorCode	calDiffC ()

virtual MoFEMErrorCode	calTheta ()

virtual MoFEMErrorCode	calSe ()

Public Attributes
int	blockId
	Block Id.

std::string	matName
	material name

double	Ks
	Saturated hydraulic conductivity [m/day].

double	hS
	minimum capillary height [m]

double	thetaS
	saturated water content

double	thetaR
	residual water contents

double	thetaM
	model parameter

double	alpha
	model parameter

double	n
	model parameter

double	Ah
	Initial hydraulic head coefficient.

double	AhZ
	Initial hydraulic head coefficient.

double	AhZZ
	Initial hydraulic head coefficient.

Public Attributes inherited from MixTransport::GenericMaterial
double	sCale
	Scale time dependent eq.

double	h
	hydraulic head

double	h_t
	rate of hydraulic head

double	K
	Hydraulic conductivity [L/s].

double	diffK
	Derivative of hydraulic conductivity [L/s * L^2/F].

double	C
	Capacity [S^2/L^2].

double	diffC
	Derivative of capacity [S^2/L^2 * L^2/F ].

double	tHeta
	Water content.

double	Se
	Effective saturation.

Range	tEts
	Elements with this material.

double	x

double	y

double	z
	in meters (L)

Additional Inherited Members
Static Public Attributes inherited from MixTransport::GenericMaterial
static double	ePsilon0 = 0
	Regularization parameter.

static double	ePsilon1 = 0
	Regularization parameter.

static double	scaleZ = 1
	Scale z direction.

Detailed Description

Examples: mofem/tutorials/cor-0to1/unsaturated_transport.cpp.

Definition at line 13 of file MaterialUnsaturatedFlow.hpp.

Member Typedef Documentation

◆ RegisterHook

typedef boost::function<boost::shared_ptr<CommonMaterialData>( const CommonMaterialData &data)> MixTransport::CommonMaterialData::RegisterHook

Definition at line 99 of file MaterialUnsaturatedFlow.hpp.

Constructor & Destructor Documentation

◆ CommonMaterialData()

MixTransport::CommonMaterialData::CommonMaterialData ( )

inline

Definition at line 15 of file MaterialUnsaturatedFlow.hpp.

                       {
    blockId = 0;
    matName = "SimpleDarcy";
    sCale = 1;
    thetaS = 0.38;
    thetaM = 0.38;
    thetaR = 0.068;
    alpha = 0.8;
    n = 1.09;
    Ks = 0.048;
    hS = 0;
    Ah = 0.;
    AhZ = 0.;
    AhZZ = 0;
  }

Member Function Documentation

◆ addOptions()

void MixTransport::CommonMaterialData::addOptions	(	po::options_description &	o,
		const std::string &	prefix
	)

inline

Definition at line 49 of file MaterialUnsaturatedFlow.hpp.

                                                                     {
    o.add_options()((prefix + ".material_name").c_str(),
                    po::value<std::string>(&matName)->default_value(matName))(
        (prefix + ".ePsilon0").c_str(),
        po::value<double>(&ePsilon0)->default_value(ePsilon0))(
        (prefix + ".ePsilon1").c_str(),
        po::value<double>(&ePsilon1)->default_value(ePsilon1))(
        (prefix + ".sCale").c_str(),
        po::value<double>(&sCale)->default_value(sCale))(
        (prefix + ".thetaS").c_str(),
        po::value<double>(&thetaS)->default_value(thetaS))(
        (prefix + ".thetaR").c_str(),
        po::value<double>(&thetaR)->default_value(thetaR))(
        (prefix + ".thetaM").c_str(),
        po::value<double>(&thetaM)->default_value(thetaM))(
        (prefix + ".alpha").c_str(),
        po::value<double>(&alpha)->default_value(alpha))(
        (prefix + ".n").c_str(), po::value<double>(&n)->default_value(n))(
        (prefix + ".hS").c_str(), po::value<double>(&hS)->default_value(hS))(
        (prefix + ".Ks").c_str(), po::value<double>(&Ks)->default_value(Ks))(
        (prefix + ".Ah").c_str(), po::value<double>(&Ah)->default_value(Ah))(
        (prefix + ".AhZ").c_str(), po::value<double>(&AhZ)->default_value(AhZ))(
        (prefix + ".AhZZ").c_str(),
        po::value<double>(&AhZZ)->default_value(AhZZ))(
        (prefix + ".scaleZ").c_str(),
        po::value<double>(&scaleZ)->default_value(scaleZ));
  }

◆ initialPcEval()

double MixTransport::CommonMaterialData::initialPcEval ( ) const

inlinevirtual

Initialize head.

Returns: value of head

Implements MixTransport::GenericMaterial.

Definition at line 47 of file MaterialUnsaturatedFlow.hpp.

47{ return Ah + AhZ * z + AhZZ * z * z; }

MixTransport::GenericMaterial::z

double z

in meters (L)

Definition UnsaturatedFlow.hpp:44

◆ printMatParameters()

void MixTransport::CommonMaterialData::printMatParameters	(	const int	id,
		const std::string &	prefix
	)		const

inlinevirtual

Implements MixTransport::GenericMaterial.

Examples: mofem/tutorials/cor-0to1/unsaturated_transport.cpp.

Definition at line 77 of file MaterialUnsaturatedFlow.hpp.

                                                                       {
    PetscPrintf(PETSC_COMM_WORLD, "Mat name %s-%s block id %d\n",
                prefix.c_str(), matName.c_str(), id);
    PetscPrintf(PETSC_COMM_WORLD, "Material name: %s\n", matName.c_str());
    PetscPrintf(PETSC_COMM_WORLD, "thetaS=%6.4g\n", thetaS);
    PetscPrintf(PETSC_COMM_WORLD, "thetaR=%6.4g\n", thetaR);
    PetscPrintf(PETSC_COMM_WORLD, "thetaM=%6.4g\n", thetaM);
    PetscPrintf(PETSC_COMM_WORLD, "alpha=%6.4g\n", alpha);
    PetscPrintf(PETSC_COMM_WORLD, "n=%6.4g\n", n);
    PetscPrintf(PETSC_COMM_WORLD, "hS=%6.4g\n", hS);
    PetscPrintf(PETSC_COMM_WORLD, "Ks=%6.4g\n", Ks);
    PetscPrintf(PETSC_COMM_WORLD, "Ah=%6.4g\n", Ah);
    PetscPrintf(PETSC_COMM_WORLD, "AhZ=%6.4g\n", AhZ);
    PetscPrintf(PETSC_COMM_WORLD, "AhZZ=%6.4g\n", AhZZ);
    PetscPrintf(PETSC_COMM_WORLD, "ePsilon0=%6.4g\n", ePsilon0);
    PetscPrintf(PETSC_COMM_WORLD, "ePsilon1=%6.4g\n", ePsilon1);
    PetscPrintf(PETSC_COMM_WORLD, "sCale=%6.4g\n", sCale);
    PetscPrintf(PETSC_COMM_WORLD, "scaleZ=%6.4g\n", scaleZ);
  }