Fluid density assuming constant bulk modulus, for p>cut_limit Then following a cubic for zero_point <= p <= cut_limit Then zero for p<zero_point. More...

#include <RichardsDensityConstBulkCut.h>

Inheritance diagram for RichardsDensityConstBulkCut:

[legend]

Public Member Functions
	RichardsDensityConstBulkCut (const InputParameters &parameters)

Real	density (Real p) const
	fluid density as a function of porepressure More...

Real	ddensity (Real p) const
	derivative of fluid density wrt porepressure More...

Real	d2density (Real p) const
	second derivative of fluid density wrt porepressure More...

void	initialize ()

void	execute ()

void	finalize ()

Protected Attributes
Real	_dens0
	density = _dens0*exp(p/_bulk), modified by cubic More...

Real	_bulk
	density = _dens0*exp(p/_bulk), modified by cubic More...

Real	_cut_limit
	where the cubic starts More...

Real	_zero_point
	where the density is zero More...

Real	_c3
	(cut_limit-zero_point)^3 More...

Detailed Description

Fluid density assuming constant bulk modulus, for p>cut_limit Then following a cubic for zero_point <= p <= cut_limit Then zero for p<zero_point.

The cubic is chosen so the function and its first derivative is continuous throughout

Definition at line 25 of file RichardsDensityConstBulkCut.h.

Constructor & Destructor Documentation

◆ RichardsDensityConstBulkCut()

RichardsDensityConstBulkCut::RichardsDensityConstBulkCut ( const InputParameters & parameters )

Definition at line 38 of file RichardsDensityConstBulkCut.C.

   : RichardsDensity(parameters),
     _dens0(getParam<Real>("dens0")),
     _bulk(getParam<Real>("bulk_mod")),
     _cut_limit(getParam<Real>("cut_limit")),
     _zero_point(getParam<Real>("zero_point")),
     _c3(std::pow(_cut_limit - _zero_point, 3))
 {
   if (_zero_point >= _cut_limit)
     mooseError("RichardsDensityConstantbulkCut: zero_point must be less than cut_limit");
 }

Member Function Documentation

◆ d2density()

Real RichardsDensityConstBulkCut::d2density ( Real p ) const

virtual

second derivative of fluid density wrt porepressure

Parameters

p	porepressure

Implements RichardsDensity.

Definition at line 82 of file RichardsDensityConstBulkCut.C.

 {
   if (p <= _zero_point)
     return 0;
  
   Real unmodified = _dens0 * std::exp(p / _bulk);
   if (p >= _cut_limit)
     return unmodified / _bulk / _bulk;
  
   Real modifier =
       (3 * _cut_limit - 2 * p - _zero_point) * (p - _zero_point) * (p - _zero_point) / _c3;
   Real dmodifier = 6 * (_cut_limit - p) * (p - _zero_point) / _c3;
   Real d2modifier = (6 * _cut_limit + 6 * _zero_point - 12 * p) / _c3;
   return (modifier / _bulk / _bulk + 2 * dmodifier / _bulk + d2modifier) * unmodified;
 }

◆ ddensity()

Real RichardsDensityConstBulkCut::ddensity ( Real p ) const

virtual

derivative of fluid density wrt porepressure

Parameters

p	porepressure

Implements RichardsDensity.

Definition at line 66 of file RichardsDensityConstBulkCut.C.

 {
   if (p <= _zero_point)
     return 0;
  
   Real unmodified = _dens0 * std::exp(p / _bulk);
   if (p >= _cut_limit)
     return unmodified / _bulk;
  
   Real modifier =
       (3 * _cut_limit - 2 * p - _zero_point) * (p - _zero_point) * (p - _zero_point) / _c3;
   Real dmodifier = 6 * (_cut_limit - p) * (p - _zero_point) / _c3;
   return (modifier / _bulk + dmodifier) * unmodified;
 }

◆ density()

Real RichardsDensityConstBulkCut::density ( Real p ) const

virtual

fluid density as a function of porepressure

Parameters

p	porepressure

Implements RichardsDensity.

Definition at line 51 of file RichardsDensityConstBulkCut.C.

 {
   if (p <= _zero_point)
     return 0;
  
   Real unmodified = _dens0 * std::exp(p / _bulk);
   if (p >= _cut_limit)
     return unmodified;
  
   Real modifier =
       (3 * _cut_limit - 2 * p - _zero_point) * (p - _zero_point) * (p - _zero_point) / _c3;
   return modifier * unmodified;
 }