RichardsSeffPrimeAux Class Reference

Calculates derivative of effective saturation wrt a specified porepressure. More...

#include <RichardsSeffPrimeAux.h>

Inheritance diagram for RichardsSeffPrimeAux:

Public Member Functions
	RichardsSeffPrimeAux (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeValue ()

Protected Attributes
const RichardsSeff &	_seff_UO
	The user object that defines effective saturation as function of porepressure (or porepressures in the multiphase situation) More...
int	_wrt1
	AuxKernel calculates d^2(Seff)/dP_wrt1 so wrt1 is the porepressure number to differentiate wrt to. More...
std::vector< const VariableValue * >	_pressure_vals
	the porepressure values (this will be length N where N is the number of arguments that the _seff_UO requires) Eg, for twophase simulations N=2 More...
std::vector< Real >	_mat
	array of derivtives: This auxkernel returns _mat[_wrt1] More...

Detailed Description

Calculates derivative of effective saturation wrt a specified porepressure.

Definition at line 25 of file RichardsSeffPrimeAux.h.

Constructor & Destructor Documentation

RichardsSeffPrimeAux()

RichardsSeffPrimeAux::RichardsSeffPrimeAux

(

const InputParameters &

parameters

)

Definition at line 31 of file RichardsSeffPrimeAux.C.

  : AuxKernel(parameters),
    _seff_UO(getUserObject<RichardsSeff>("seff_UO")),
    _wrt1(getParam<int>("wrtnum"))
{
  int n = coupledComponents("pressure_vars");
  if (_wrt1 < 0 || _wrt1 >= n)
    mooseError("Your wrtnum is ", _wrt1, " but it must obey 0 <= wrtnum < ", n, ".");
  _pressure_vals.resize(n);
 
  for (int i = 0; i < n; ++i)
    _pressure_vals[i] = &coupledValue("pressure_vars", i);
 
  _mat.resize(n);
}

Member Function Documentation

computeValue()

Real RichardsSeffPrimeAux::computeValue ( )

protectedvirtual

Definition at line 48 of file RichardsSeffPrimeAux.C.

{
  _seff_UO.dseff(_pressure_vals, _qp, _mat);
  return _mat[_wrt1];
}

Member Data Documentation

_mat

std::vector<Real> RichardsSeffPrimeAux::_mat

protected

array of derivtives: This auxkernel returns _mat[_wrt1]

Definition at line 54 of file RichardsSeffPrimeAux.h.

Referenced by computeValue(), and RichardsSeffPrimeAux().

_pressure_vals

std::vector<const VariableValue *> RichardsSeffPrimeAux::_pressure_vals

protected

the porepressure values (this will be length N where N is the number of arguments that the _seff_UO requires) Eg, for twophase simulations N=2

Definition at line 51 of file RichardsSeffPrimeAux.h.

Referenced by computeValue(), and RichardsSeffPrimeAux().

_seff_UO

const RichardsSeff& RichardsSeffPrimeAux::_seff_UO

protected

The user object that defines effective saturation as function of porepressure (or porepressures in the multiphase situation)

Definition at line 38 of file RichardsSeffPrimeAux.h.

Referenced by computeValue().

_wrt1

int RichardsSeffPrimeAux::_wrt1

protected

AuxKernel calculates d^2(Seff)/dP_wrt1 so wrt1 is the porepressure number to differentiate wrt to.

Definition at line 44 of file RichardsSeffPrimeAux.h.

Referenced by computeValue(), and RichardsSeffPrimeAux().

---

The documentation for this class was generated from the following files:

• richards/include/auxkernels/RichardsSeffPrimeAux.h
• richards/src/auxkernels/RichardsSeffPrimeAux.C