#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.

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.

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
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.

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 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: mofem/tutorials/cor-0to1/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: mofem/tutorials/cor-0to1/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: mofem/tutorials/cor-0to1/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();
  }