UFOperators2D::BlockData Struct Reference

#include <users_modules/softmech/unsaturated_2Dflow/src/Unsaturated2DFlowUDOP.hpp>

Collaboration diagram for UFOperators2D::BlockData:

Public Member Functions
	BlockData ()
MoFEMErrorCode	setOptions ()
template<typename T >
T	get_waterContent (T head)
template<typename T >
T	get_conductivity (T head, double head_norm)
template<typename T >
T	get_relativeConductivity (T head)
template<typename T >
T	get_Fe (T eff_satu)
template<typename T >
T	get_eff_satuStar (T eff_satu)
template<typename T >
T	get_effSaturation (T head)
template<typename T >
T	get_capacity (T head, double head_norm)
double	get_diffCapacity (double head)
double	get_diffConductivity (double head)

Public Attributes
int	block_id
double	K_s
double	ll
double	theta_r
double	theta_s
double	theta_m
double	nn
double	alpha
double	h_s
Range	block_ents

Detailed Description

Definition at line 26 of file Unsaturated2DFlowUDOP.hpp.

Constructor & Destructor Documentation

BlockData()

UFOperators2D::BlockData::BlockData ( )

inline

Definition at line 39 of file Unsaturated2DFlowUDOP.hpp.

      : K_s(1.00), ll(0.5000), theta_r(0.0450), theta_s(0.4300),
        theta_m(0.4300000), nn(5.3800000), alpha(1.812500), h_s(0.0000) {}

Member Function Documentation

get_capacity()

template<typename T >

T UFOperators2D::BlockData::get_capacity ( T  head, double  head_norm  )

inline

Definition at line 102 of file Unsaturated2DFlowUDOP.hpp.

                                                                 {
    const T m = 1 - 1.0 / nn;
    T ret_val;
    if (head_norm < h_s) {
      ret_val = alpha * m * nn * (theta_m - theta_r) *
                pow(-alpha * head, nn - 1) /
                pow(1 + pow(-alpha * head, nn), m + 1);
    } else {
      ret_val = 0;
    }
    return ret_val;
  }

get_conductivity()

template<typename T >

T UFOperators2D::BlockData::get_conductivity ( T  head, double  head_norm  )

inline

Definition at line 71 of file Unsaturated2DFlowUDOP.hpp.

                                                                     {
    T ret_val;
    if (head_norm < h_s) {
      ret_val = K_s * get_relativeConductivity(head);
    } else {
      ret_val = K_s;
    }
    return ret_val;
  }

get_diffCapacity()

double UFOperators2D::BlockData::get_diffCapacity ( double  head )

inline

Definition at line 115 of file Unsaturated2DFlowUDOP.hpp.

                                       {
    std::complex<double> cx_head = head + eps * 1i;
    auto capacity = get_capacity(cx_head, head);
    return capacity.imag() / eps;
  }

get_diffConductivity()

double UFOperators2D::BlockData::get_diffConductivity ( double  head )

inline

Definition at line 121 of file Unsaturated2DFlowUDOP.hpp.

                                           {
    std::complex<double> cx_head = head + eps * 1i;
    auto conductivity = get_conductivity(cx_head, head);
    return conductivity.imag() / eps;
  }

get_eff_satuStar()

template<typename T >

T UFOperators2D::BlockData::get_eff_satuStar ( T  eff_satu )

inline

Definition at line 93 of file Unsaturated2DFlowUDOP.hpp.

                                                       {
    return (theta_s - theta_r) / (theta_m - theta_r) * eff_satu;
  }

get_effSaturation()

template<typename T >

T UFOperators2D::BlockData::get_effSaturation ( T  head )

inline

Definition at line 97 of file Unsaturated2DFlowUDOP.hpp.

                                                    {
    const T theta = get_waterContent(head);
    return (theta - theta_r) / (theta_s - theta_r);
  }

get_Fe()

template<typename T >

T UFOperators2D::BlockData::get_Fe ( T  eff_satu )

inline

Definition at line 88 of file Unsaturated2DFlowUDOP.hpp.

                                             {
    const T m = 1 - 1.0 / nn;
    const T S_eStar = get_eff_satuStar(eff_satu);
    return pow(1 - pow(S_eStar, 1.0 / m), m);
  }

get_relativeConductivity()

template<typename T >

T UFOperators2D::BlockData::get_relativeConductivity ( T  head )

inline

Definition at line 81 of file Unsaturated2DFlowUDOP.hpp.

                                                           {
    const T S_e = get_effSaturation(head);
    const T F_e = get_Fe(S_e);
    const T F_1 = get_Fe(1.0);
    return pow(S_e, ll) * pow(((1 - F_e) / (1 - F_1)), 2);
  }

get_waterContent()

template<typename T >

T UFOperators2D::BlockData::get_waterContent ( T  head )

inline

Definition at line 64 of file Unsaturated2DFlowUDOP.hpp.

                                                   {
    T ret_val;
    T m = 1 - 1.0 / nn;
    return ret_val = theta_r +
                     (theta_m - theta_r) / pow(1 + pow(-alpha * head, nn), m);
  }

setOptions()

MoFEMErrorCode UFOperators2D::BlockData::setOptions ( )

inline

Definition at line 43 of file Unsaturated2DFlowUDOP.hpp.

                              {
    MoFEMFunctionBegin;
    ierr = PetscOptionsBegin(PETSC_COMM_WORLD, "", "Material options", "none");
    CHKERR PetscOptionsScalar("-K_s", "K_s", "", K_s, &K_s, PETSC_NULL);
    CHKERR PetscOptionsScalar("-saturated_conductivity", "K_s", "", K_s, &K_s,
                              PETSC_NULL);
    CHKERR PetscOptionsScalar("-l", "l", "", ll, &ll, PETSC_NULL);
    CHKERR PetscOptionsScalar("-theta_r", "theta_r", "", theta_r, &theta_r,
                              PETSC_NULL);
    CHKERR PetscOptionsScalar("-theta_s", "theta_s", "", theta_s, &theta_s,
                              PETSC_NULL);
    CHKERR PetscOptionsScalar("-theta_m", "theta_m", "", theta_m, &theta_m,
                              PETSC_NULL);
    CHKERR PetscOptionsScalar("-n", "n", "", nn, &nn, PETSC_NULL);
    CHKERR PetscOptionsScalar("-alpha", "alpha", "", alpha, &alpha, PETSC_NULL);
    CHKERR PetscOptionsScalar("-h_s", "h_s", "", h_s, &h_s, PETSC_NULL);
    ierr = PetscOptionsEnd();
    CHKERRG(ierr);
    MoFEMFunctionReturn(0);
  };

Member Data Documentation

alpha

double UFOperators2D::BlockData::alpha

Definition at line 34 of file Unsaturated2DFlowUDOP.hpp.

block_ents

Range UFOperators2D::BlockData::block_ents

Definition at line 37 of file Unsaturated2DFlowUDOP.hpp.

block_id

int UFOperators2D::BlockData::block_id

Definition at line 27 of file Unsaturated2DFlowUDOP.hpp.

h_s

double UFOperators2D::BlockData::h_s

Definition at line 35 of file Unsaturated2DFlowUDOP.hpp.

K_s

double UFOperators2D::BlockData::K_s

Definition at line 28 of file Unsaturated2DFlowUDOP.hpp.

ll

double UFOperators2D::BlockData::ll

Definition at line 29 of file Unsaturated2DFlowUDOP.hpp.

nn

double UFOperators2D::BlockData::nn

Definition at line 33 of file Unsaturated2DFlowUDOP.hpp.

theta_m

double UFOperators2D::BlockData::theta_m

Definition at line 32 of file Unsaturated2DFlowUDOP.hpp.

theta_r

double UFOperators2D::BlockData::theta_r

Definition at line 30 of file Unsaturated2DFlowUDOP.hpp.

theta_s

double UFOperators2D::BlockData::theta_s

Definition at line 31 of file Unsaturated2DFlowUDOP.hpp.

---

The documentation for this struct was generated from the following file:

• users_modules/softmech/unsaturated_2Dflow/src/Unsaturated2DFlowUDOP.hpp