DFChaventMobility Class Reference

#include <dfchaventmobility.h>

Inheritance diagram for DFChaventMobility:

[legend]

Collaboration diagram for DFChaventMobility:

[legend]

List of all members.

Public Member Functions
	DFChaventMobility (double vw, double vo, double Srw, double Sro)
virtual	~DFChaventMobility ()
virtual double	operator() (double x, unsigned cmp=0) const
virtual void	getMinMaxValues (double a, double b, double &min, double &max) const
virtual double	getMaxNorm (double a, double b) const
void	updateInflectionPoint ()
double	getInflectionDerivativePoint (double M)
virtual void	setParameters (const VecDouble &v)
Protected Attributes
double	m_M
double	MinSo
double	MinSw
double	MaxSo
double	MaxSw
double	_c
double	_d
double	_a1
double	_a2
double	_a3
double	m_inflectPoint
double	m_inflectPointDeriv

---

Detailed Description

Definition at line 9 of file dfchaventmobility.h.

---

Constructor & Destructor Documentation

DFChaventMobility::DFChaventMobility	(	double	vw,
		double	vo,
		double	Srw,
		double	Sro
	)

Definition at line 4 of file dfchaventmobility.cpp.

:Function1D(1)
{
  m_M=vw/vo;

  MaxSo=1-Srw;
  MaxSw=1-Sro;
  MinSw=Srw;
  MinSo=Sro;


  //1_Srw2=1_Srw*1_Srw;
  //1_Sro=1-Sro;

  _c = m_M*MaxSo*MaxSo/(MaxSw*MaxSw);
  _d = 2*_c*(MaxSw-MinSw);
  _a1=_c+1;
  _a2=-2*(MinSw + _c*MaxSw);
  _a3=_c*MaxSw*MaxSw+MinSw*MinSw;

  m_inflectPoint = getInflectionDerivativePoint(m_M);
  m_inflectPointDeriv=(*this)(m_inflectPoint);
}

DFChaventMobility::~DFChaventMobility

(

)

[virtual]

Definition at line 125 of file dfchaventmobility.cpp.

{

}

---

Member Function Documentation

double DFChaventMobility::getInflectionDerivativePoint

(

double

M

)

Definition at line 94 of file dfchaventmobility.cpp.

{
  static double sqrt3 = sqrt(3.0);
  

  double dd=MaxSw-MinSw;
  double dd2=dd*dd;
  double dd3=dd*dd2;

  double b = -sqrt(_c)*dd3/(4*(_c+1));
  double a = (1-_c)*dd3/(8*(_c+1));
  double r=pow(hypot(a,b),1.0/3.0);
  double theta=atan2(b,a)/3.0;
  a=r*cos(theta);
  b=r*sin(theta);
  double aux1=(-sqrt3*b-a)/2.0;
    
  return aux1*dd2/(4*(b*b + a*a)) + aux1 + (MaxSw+MinSw)/2.0;
  
  
  
}

double DFChaventMobility::getMaxNorm	(	double	a,
		double	b
	)			const [virtual]

Reimplemented from Function1D.

Definition at line 65 of file dfchaventmobility.cpp.

{
  assert(a<=b);
  if (a<0.0) a=0.0;
  if (b>1.0) b=1.0;
  if (b< m_inflectPoint)
  {
    return (*this)(b,0);
  }
  else if (a>m_inflectPoint)
  {
    return (*this)(a,0);
  }
  else
  {
    return m_inflectPointDeriv;
  }

}

void DFChaventMobility::getMinMaxValues	(	double	a,
		double	b,
		double &	min,
		double &	max
	)			const [virtual]

Reimplemented from Function1D.

Definition at line 41 of file dfchaventmobility.cpp.

{
  
  assert(a<=b);
  if (a<0.0) a=0.0;
  if (b>1.0) b=1.0;
  if (b< m_inflectPoint)
  {
    min=(*this)(a,0);
    max=(*this)(b,0);
  }
  else if (a>m_inflectPoint)
  {
    max=(*this)(a,0);
    min=(*this)(b,0);
  }
  else
  {
    max=m_inflectPointDeriv;
    min=std::min( (*this)(a,0), (*this)(b,0));
  }
}

double DFChaventMobility::operator()	(	double	x,
		unsigned	cmp = 0
	)			const [virtual]

Implements Function1D.

Definition at line 31 of file dfchaventmobility.cpp.

{
   if (x >= MaxSw || x <= MinSw ) 
     return 0.0;
   register double quot=(_a1*x + _a2)*x + _a3;
   quot*=quot;
   return _d*(MaxSw-x)*(x-MinSw)/quot;
}

void DFChaventMobility::setParameters

(

const VecDouble &

v

)

[virtual]

Reimplemented from GeneralFunctionInterface.

Definition at line 118 of file dfchaventmobility.cpp.

{
  assert(v.size() ==2);
  m_M=v(0)/v(1);
}

void DFChaventMobility::updateInflectionPoint

(

)

Definition at line 85 of file dfchaventmobility.cpp.

{
  m_inflectPoint = getInflectionDerivativePoint(m_M);
  m_inflectPointDeriv=(*this)(m_inflectPoint);
}

---

Member Data Documentation

double DFChaventMobility::_a1 [protected]

Definition at line 16 of file dfchaventmobility.h.

double DFChaventMobility::_a2 [protected]

Definition at line 16 of file dfchaventmobility.h.

double DFChaventMobility::_a3 [protected]

Definition at line 16 of file dfchaventmobility.h.

double DFChaventMobility::_c [protected]

Definition at line 16 of file dfchaventmobility.h.

double DFChaventMobility::_d [protected]

Definition at line 16 of file dfchaventmobility.h.

double DFChaventMobility::m_inflectPoint [protected]

Definition at line 18 of file dfchaventmobility.h.

double DFChaventMobility::m_inflectPointDeriv [protected]

Definition at line 19 of file dfchaventmobility.h.

double DFChaventMobility::m_M [protected]

Definition at line 14 of file dfchaventmobility.h.

double DFChaventMobility::MaxSo [protected]

Definition at line 15 of file dfchaventmobility.h.

double DFChaventMobility::MaxSw [protected]

Definition at line 15 of file dfchaventmobility.h.

double DFChaventMobility::MinSo [protected]

Definition at line 15 of file dfchaventmobility.h.

double DFChaventMobility::MinSw [protected]

Definition at line 15 of file dfchaventmobility.h.

---

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

• dfchaventmobility.h
• dfchaventmobility.cpp