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

Inheritance diagram for MixTransport::MaterialVanGenuchten:

Collaboration diagram for MixTransport::MaterialVanGenuchten:

Public Member Functions
	MaterialVanGenuchten (const CommonMaterialData &data)

template<typename TYPE >
TYPE	funTheta (TYPE &h, const double m)

template<typename TYPE >
TYPE	funFunSeStar (TYPE &SeStar, const double m)

adouble	funKr (adouble &ah)

virtual void	recordTheta ()

virtual void	recordKr ()

void	printTheta (const double b, const double e, double s, const std::string &prefix)

void	printKappa (const double b, const double e, double s, const std::string &prefix)

void	printC (const double b, const double e, double s, const std::string &prefix)

Public Member Functions inherited from MixTransport::MaterialWithAutomaticDifferentiation
	MaterialWithAutomaticDifferentiation (const CommonMaterialData &data)

template<typename TYPE >
TYPE	funSe (TYPE &theta)

MoFEMErrorCode	calK ()

MoFEMErrorCode	calDiffK ()

MoFEMErrorCode	calC ()

MoFEMErrorCode	calDiffC ()

MoFEMErrorCode	calTheta ()

MoFEMErrorCode	calSe ()

Public Member Functions inherited from MixTransport::CommonMaterialData
	CommonMaterialData ()

double	initialPcEval () const
	Initialize head. More...

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 ()

Static Public Member Functions
static boost::shared_ptr< CommonMaterialData >	createMatPtr (const CommonMaterialData &data)

Public Attributes
adouble	ah

adouble	aTheta

adouble	aKr

adouble	aSe

adouble	aSeStar

Public Attributes inherited from MixTransport::MaterialWithAutomaticDifferentiation
double	Kr

double	diffKr

Public Attributes inherited from MixTransport::CommonMaterialData
int	blockId
	Block Id. More...

std::string	matName
	material name More...

double	Ks
	Saturated hydraulic conductivity [m/day]. More...

double	hS
	minimum capillary height [m] More...

double	thetaS
	saturated water content More...

double	thetaR
	residual water contents More...

double	thetaM
	model parameter More...

double	alpha
	model parameter More...

double	n
	model parameter More...

double	Ah
	Initial hydraulic head coefficient. More...

double	AhZ
	Initial hydraulic head coefficient. More...

double	AhZZ
	Initial hydraulic head coefficient. More...

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

double	h
	hydraulic head More...

double	h_t
	rate of hydraulic head More...

double	K
	Hydraulic conductivity [L/s]. More...

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

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

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

double	tHeta
	Water content. More...

double	Se
	Effective saturation. More...

Range	tEts
	Elements with this material. More...

double	x

double	y

double	z
	in meters (L) More...

Additional Inherited Members
Public Types inherited from MixTransport::CommonMaterialData
typedef boost::function< boost::shared_ptr< CommonMaterialData > const CommonMaterialData &data)>	RegisterHook

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

static double	ePsilon1 = 0
	Regularization parameter. More...

static double	scaleZ = 1
	Scale z direction. More...

Detailed Description

Examples: unsaturated_transport.cpp.

Definition at line 254 of file MaterialUnsaturatedFlow.hpp.

Constructor & Destructor Documentation

◆ MaterialVanGenuchten()

MixTransport::MaterialVanGenuchten::MaterialVanGenuchten ( const CommonMaterialData & data )

inline

Definition at line 262 of file MaterialUnsaturatedFlow.hpp.

       : MaterialWithAutomaticDifferentiation(data) {
     recordTheta();
     recordKr();
   }

Member Function Documentation

◆ createMatPtr()

static boost::shared_ptr<CommonMaterialData> MixTransport::MaterialVanGenuchten::createMatPtr ( const CommonMaterialData & data )

inlinestatic

Definition at line 257 of file MaterialUnsaturatedFlow.hpp.

                                                {
     return boost::shared_ptr<CommonMaterialData>(
         new MaterialVanGenuchten(data));
   }

◆ funFunSeStar()

template<typename TYPE >

TYPE MixTransport::MaterialVanGenuchten::funFunSeStar	(	TYPE &	SeStar,
		const double	m
	)

inline

Definition at line 279 of file MaterialUnsaturatedFlow.hpp.

                                                          {
     return pow(1 - pow(SeStar, 1 / m), m);
   }

◆ funKr()

adouble MixTransport::MaterialVanGenuchten::funKr ( adouble & ah )

inline

Definition at line 283 of file MaterialUnsaturatedFlow.hpp.

                                     {
     const double m = 1 - 1 / n;
     aTheta = funTheta(ah, m);
     aSe = funSe(aTheta);
     aSeStar = aSe * (thetaS - thetaR) / (thetaM - thetaR);
     double one = 1;
     const double c = funFunSeStar<double>(one, m);
     return sqrt(aSe) *
            pow((1 - funFunSeStar<adouble>(aSeStar, m)) / (1 - c), 2);
   }

◆ funTheta()

template<typename TYPE >

TYPE MixTransport::MaterialVanGenuchten::funTheta	(	TYPE &	h,
		const double	m
	)

inline

Definition at line 274 of file MaterialUnsaturatedFlow.hpp.

                                                                          {
     return thetaR + (thetaM - thetaR) / pow(1 + pow(-alpha * h, n), m);
   }

◆ printC()

void MixTransport::MaterialVanGenuchten::printC	(	const double	b,
		const double	e,
		double	s,
		const std::string &	prefix
	)

inline

Examples: unsaturated_transport.cpp.

Definition at line 339 of file MaterialUnsaturatedFlow.hpp.

                                        {
     h = b;
     for (; h >= e; h += s) {
       s = -pow(-s, 0.9);
       calC();
       PetscPrintf(PETSC_COMM_SELF, "%s %6.4e %6.4e\n", prefix.c_str(), h, C);
     }
   }

◆ printKappa()

void MixTransport::MaterialVanGenuchten::printKappa	(	const double	b,
		const double	e,
		double	s,
		const std::string &	prefix
	)

inline

Examples: unsaturated_transport.cpp.

Definition at line 328 of file MaterialUnsaturatedFlow.hpp.

                                            {
     h = b;
     for (; h >= e; h += s) {
       s = -pow(-s, 0.9);
       calK();
       PetscPrintf(PETSC_COMM_SELF, "%s %6.4e %6.4e %6.4e\n", prefix.c_str(), h,
                   Kr, K);
     }
   }

◆ printTheta()

void MixTransport::MaterialVanGenuchten::printTheta	(	const double	b,
		const double	e,
		double	s,
		const std::string &	prefix
	)

inline

Examples: unsaturated_transport.cpp.

Definition at line 315 of file MaterialUnsaturatedFlow.hpp.

                                            {
     const double m = 1 - 1 / n;
     h = b;
     for (; h >= e; h += s) {
       s = -pow(-s, 0.9);
       double theta = funTheta(h, m);
       double Se = (theta - thetaR) / (thetaS - thetaR);
       PetscPrintf(PETSC_COMM_SELF, "%s %6.4e %6.4e %6.4e\n", prefix.c_str(), h,
                   theta, Se);
     }
   }

◆ recordKr()

virtual void MixTransport::MaterialVanGenuchten::recordKr ( )

inlinevirtual

Implements MixTransport::MaterialWithAutomaticDifferentiation.

Definition at line 305 of file MaterialUnsaturatedFlow.hpp.

                           {
     h = -1 - hS;
     trace_on(2 * blockId + 1, true);
     ah <<= h;
     aKr = funKr(ah);
     double r_Kr;
     aKr >>= r_Kr;
     trace_off();
   }

◆ recordTheta()

virtual void MixTransport::MaterialVanGenuchten::recordTheta ( )

inlinevirtual

Implements MixTransport::MaterialWithAutomaticDifferentiation.

Definition at line 294 of file MaterialUnsaturatedFlow.hpp.

                              {
     h = -1 - hS;
     trace_on(2 * blockId + 0, true);
     ah <<= h;
     const double m = 1 - 1 / n;
     aTheta = funTheta(ah, m);
     double r_theta;
     aTheta >>= r_theta;
     trace_off();
   }