RichardsSeffPrimePrimeAux Class Reference

Calculates derivative of effective saturation wrt specified porepressures. More...

#include <RichardsSeffPrimePrimeAux.h>

Inheritance diagram for RichardsSeffPrimePrimeAux:

Public Member Functions
	RichardsSeffPrimePrimeAux (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 dP_wrt2 so wrt1 is the porepressure number to differentiate wrt to. More...
int	_wrt2
	AuxKernel calculates d^2(Seff)/dP_wrt1 dP_wrt2 so wrt2 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< std::vector< Real > >	_mat
	matrix of 2nd derivtives: This auxkernel returns _mat[_wrt1][_wrt2]; More...

Detailed Description

Calculates derivative of effective saturation wrt specified porepressures.

Definition at line 26 of file RichardsSeffPrimePrimeAux.h.

Constructor & Destructor Documentation

RichardsSeffPrimePrimeAux()

RichardsSeffPrimePrimeAux::RichardsSeffPrimePrimeAux

(

const InputParameters &

parameters

)

Definition at line 34 of file RichardsSeffPrimePrimeAux.C.

  : AuxKernel(parameters),
    _seff_UO(getUserObject<RichardsSeff>("seff_UO")),
    _wrt1(getParam<int>("wrtnum1")),
    _wrt2(getParam<int>("wrtnum2"))
{
  int n = coupledComponents("pressure_vars");
  if (_wrt1 < 0 || _wrt1 >= n)
    mooseError("Your wrtnum1 is ", _wrt1, " but it must obey 0 <= wrtnum1 < ", n, ".");
  if (_wrt2 < 0 || _wrt2 >= n)
    mooseError("Your wrtnum2 is ", _wrt2, " but it must obey 0 <= wrtnum2 < ", n, ".");
  _pressure_vals.resize(n);

  for (int i = 0; i < n; ++i)
    _pressure_vals[i] = &coupledValue("pressure_vars", i);

  _mat.resize(n);
  for (int i = 0; i < n; ++i)
    _mat[i].resize(n);
}

Member Function Documentation

computeValue()

Real RichardsSeffPrimePrimeAux::computeValue ( )

protectedvirtual

Definition at line 56 of file RichardsSeffPrimePrimeAux.C.

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

Member Data Documentation

_mat

std::vector<std::vector<Real> > RichardsSeffPrimePrimeAux::_mat

protected

matrix of 2nd derivtives: This auxkernel returns _mat[_wrt1][_wrt2];

Definition at line 61 of file RichardsSeffPrimePrimeAux.h.

Referenced by computeValue(), and RichardsSeffPrimePrimeAux().

_pressure_vals

std::vector<const VariableValue *> RichardsSeffPrimePrimeAux::_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 58 of file RichardsSeffPrimePrimeAux.h.

Referenced by computeValue(), and RichardsSeffPrimePrimeAux().

_seff_UO

const RichardsSeff& RichardsSeffPrimePrimeAux::_seff_UO

protected

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

Definition at line 39 of file RichardsSeffPrimePrimeAux.h.

Referenced by computeValue().

_wrt1

int RichardsSeffPrimePrimeAux::_wrt1

protected

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

Definition at line 45 of file RichardsSeffPrimePrimeAux.h.

Referenced by computeValue(), and RichardsSeffPrimePrimeAux().

_wrt2

int RichardsSeffPrimePrimeAux::_wrt2

protected

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

Definition at line 51 of file RichardsSeffPrimePrimeAux.h.

Referenced by computeValue(), and RichardsSeffPrimePrimeAux().

---

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

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