TraceFaceFluxWithGravity Class Reference

#include <tracefacefluxwithgravity.h>

Inheritance diagram for TraceFaceFluxWithGravity:

[legend]

Collaboration diagram for TraceFaceFluxWithGravity:

[legend]

List of all members.

Public Member Functions
	TraceFaceFluxWithGravity (unsigned numFaces, Function1D &fMobW, Function1D &fMobGrav, VecDouble &cK, double gc)
	~TraceFaceFluxWithGravity ()
virtual double	fluxAtFace (OrthoMesh::Face_It &faceIt, int face, int cell1, int cell2, double &u1, double &u2)
virtual void	setVelocity (DynamicBase &dynMod)
virtual void	setSw (VecDouble &cSw)
virtual VecDouble &	getSw ()
virtual void	maxGlobalCharVelocity (double vel[3])
virtual double	maxLocalCharVelocity (OrthoMesh::Face_It &faceIt, int face, int cell1, int cell2, double &u1, double &u2)
virtual	~FaceFluxFunction ()
Private Attributes
Function1D &	m_fMobW
Function1D &	m_fMobGrav
VecDouble &	m_cK
VecDouble *	m_pcSw
double	m_gc
double	m_maxK
Matrix	m_Vel

---

Detailed Description

Implement the flux function of the dissolved CO2 equation given by the expression: ((1-mobW(Sw))*Vdt + mobGrav(Sw)*gc*Grad(Z))/(1-Sw)

IMPORTANT: This flux function assumers that the domain of the functions mobW e mobGrav is beetween [0,1] and that the mobW is a non decreasing function ( derivative is always greater than zero). Remeber that grad (Z) == <0, 0, 1>.

Definition at line 16 of file tracefacefluxwithgravity.h.

---

Constructor & Destructor Documentation

TraceFaceFluxWithGravity::TraceFaceFluxWithGravity	(	unsigned	numFaces,
		Function1D &	fMobW,
		Function1D &	fMobGrav,
		VecDouble &	cK,
		double	gc
	)

Definition at line 4 of file tracefacefluxwithgravity.cpp.

  :m_fMobW(fMobW),m_fMobGrav(fMobGrav),m_cK(cK),m_Vel(numFaces,3)
{
  m_gc = gc;
  m_maxK = cK.linfty_norm();
}

TraceFaceFluxWithGravity::~TraceFaceFluxWithGravity

(

)

virtual TraceFaceFluxWithGravity::~FaceFluxFunction

(

)

[inline, virtual]

Reimplemented from FaceFluxFunction.

Definition at line 37 of file tracefacefluxwithgravity.h.

{}

---

Member Function Documentation

double TraceFaceFluxWithGravity::fluxAtFace	(	OrthoMesh::Face_It &	faceIt,
		int	face,
		int	cell1,
		int	cell2,
		double &	u1,
		double &	u2
	)			[virtual]

This method compute the flux at a face using the values of u1 and u2 located in the two cells that share the face.

Parameters:

	faceIt	The face iteratator of the OrthoMesh.
	face	Index of the face
	cell1	Index of the cell that contains the face
	cell2	Index of the second cell that contains the face
	u1	Value of the numeric solution for cell1
	u2	Value of the numeric solution for cell2

Returns:

The flux at face

Definition at line 21 of file tracefacefluxwithgravity.cpp.

{
  assert( (unsigned) face < m_Vel.m());
  assert( (unsigned) cell1 < m_cK.size());
  assert( (unsigned) cell2 < m_cK.size());
  assert(m_cK.size() == getSw().size());
  VecDouble &cVSw = getSw();
  double S1 = cVSw(cell1);
  double S2 = cVSw(cell2);
  
  double dMobW1 = 1-m_fMobW(S1);
  double dMobW2 = 1-m_fMobW(S2);
  double *vel = &m_Vel(face,0);
  double flux = (dMobW1+dMobW2)/2.0*faceIt->normal_multiply(vel);
  
  //If the face has normal in the Z direction, then add the gravity term.
  if (faceIt->getNormalNonZeroComponent() == Z && !faceIt->at_boundary() )
  {
    VecDouble &cK = getPermeability();
    double K1 = cK(cell1);
    double K2 = cK(cell2);
    flux += (m_fMobGrav(S1) + m_fMobGrav(S2)) * K1*K2/(K1+K2) * m_gc; 
  }
  return (u1+u2)/(S1+S2)*flux/(S;
}

virtual VecDouble& TraceFaceFluxWithGravity::getSw

(

)

[inline, virtual]

Definition at line 34 of file tracefacefluxwithgravity.h.

{assert(m_pcS);return *m_pcSw;}

void TraceFaceFluxWithGravity::maxGlobalCharVelocity

(

double

vel[3]

)

[virtual]

Reimplemented from FaceFluxFunction.

Definition at line 49 of file tracefacefluxwithgravity.cpp.

{
  VecDouble minVel(3);
  VecDouble maxVel(3);

  NumericMethods::getMinMaxValueByColumn(m_Vel,minVel,maxVel);
  vel[X]= NumericMethods::maxMod(minVel(X),maxVel(X))*m_fMobW.getMaxNormGrad(0,1);
  vel[Y]= NumericMethods::maxMod(minVel(Y),maxVel(Y))*m_fMobW.getMaxNormGrad(0,1);
  vel[Z]= NumericMethods::maxMod(minVel(Z),maxVel(Z))*m_fMobW.getMaxNormGrad(0,1) + m_maxK*fabs(m_fMobGrav.getMaxNormGrad(0,1)*m_gc);
}

double TraceFaceFluxWithGravity::maxLocalCharVelocity	(	OrthoMesh::Face_It &	faceIt,
		int	face,
		int	cell1,
		int	cell2,
		double &	u1,
		double &	u2
	)			[virtual]

Definition at line 62 of file tracefacefluxwithgravity.cpp.

{
  assert(u1<=u2);
  double *vel = &(m_Vel(face,0));
  
  if (faceIt->getNormalNonZeroComponent() != Z)
    return  faceIt->normal_multiply(vel)*m_fMobW.getMaxNormGrad(u1,u2);
  else
  {
    double term1 = faceIt->normal_multiply(vel)*m_fMobW.getMaxNormGrad(u1,u2);

    //Get the max permeability
    double K = fmax(getPermeability()(cell1), getPermeability()(cell2));

    //The m_fMobGrav function can have negative derivative.
    //An the gravity term of the flux can be negative or positive.
    double dMax,dMin;
    m_fMobGrav.getMinMaxGradValues(u1,u2,dMin,dMax); 

    double result1 = term1 + K*dMin*m_gc;
    double result2 = term1 + K*dMax*m_gc;

    if (fabs(result1) < fabs(result2))
      return result2;
    else
      return result1;
  }
}

virtual void TraceFaceFluxWithGravity::setSw

(

VecDouble &

cSw

)

[inline, virtual]

Definition at line 33 of file tracefacefluxwithgravity.h.

{m_pcSw = &cSw;}

void TraceFaceFluxWithGravity::setVelocity

(

DynamicBase &

dynMod

)

[virtual]

Definition at line 92 of file tracefacefluxwithgravity.cpp.

{
  dynMod.getVelocitiesAtFaces(m_Vel);
}

---

Member Data Documentation

VecDouble& TraceFaceFluxWithGravity::m_cK [private]

Definition at line 20 of file tracefacefluxwithgravity.h.

Function1D & TraceFaceFluxWithGravity::m_fMobGrav [private]

Definition at line 19 of file tracefacefluxwithgravity.h.

Function1D& TraceFaceFluxWithGravity::m_fMobW [private]

Definition at line 19 of file tracefacefluxwithgravity.h.

double TraceFaceFluxWithGravity::m_gc [private]

Definition at line 22 of file tracefacefluxwithgravity.h.

double TraceFaceFluxWithGravity::m_maxK [private]

Definition at line 23 of file tracefacefluxwithgravity.h.

VecDouble* TraceFaceFluxWithGravity::m_pcSw [private]

Definition at line 21 of file tracefacefluxwithgravity.h.

Matrix TraceFaceFluxWithGravity::m_Vel [private]

Definition at line 24 of file tracefacefluxwithgravity.h.

---

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

• tracefacefluxwithgravity.h
• tracefacefluxwithgravity.cpp