RussanovForSystem Class Reference

#include <russanovforsystem.h>

Inheritance diagram for RussanovForSystem:

[legend]

Collaboration diagram for RussanovForSystem:

[legend]

List of all members.

Public Member Functions
	~RussanovForSystem ()
	RussanovForSystem (OrthoMesh &mesh, Function3D &fInitU, const VecDouble &cPor, Function3D &fPrescribedU, FaceFluxFunction &flux, FixedValueCondition &fixedC, double CFL)
virtual void	iterateN (unsigned nSteps, double dt)

---

Detailed Description

Definition at line 9 of file russanovforsystem.h.

---

Constructor & Destructor Documentation

RussanovForSystem::~RussanovForSystem

(

)

Definition at line 12 of file russanovforsystem.cpp.

{

}

RussanovForSystem::RussanovForSystem	(	OrthoMesh &	mesh,
		Function3D &	fInitU,
		const VecDouble &	cPor,
		Function3D &	fPrescribedU,
		FaceFluxFunction &	flux,
		FixedValueCondition &	fixedC,
		double	CFL
	)

Definition at line 4 of file russanovforsystem.cpp.

  :LaxFriedrichsForSystem(mesh,fInitU,cPor,fPrescribedU,flux,fixedC,CFL)
{
  
}

---

Member Function Documentation

void RussanovForSystem::iterateN	(	unsigned	nSteps,
		double	dt
	)			[virtual]

Reimplemented from LaxFriedrichsForSystem.

Definition at line 20 of file russanovforsystem.cpp.

{
  //To contain the left and right indices of the adjacent cells
  //for each face
  unsigned posCell,negCell;
  double a;
  //We iterate the method along faces, so 
  //alloc the vectors to contain the solutions in neg
  //and pos cells of the faces.  Remeber that the solution
  //in a cell is vector with a value for each component
  static VecDouble  vBC(n_components),vFlux(n_components),vA(n_components);

  VecDoubleRef  vQPos(n_components),vQNeg(n_components);

  //For each time step
  for (unsigned count = 0; count < nSteps; count++)
  {
    //Get the face iterator and the number of faces.
    //OrthoMesh::Cash_Face_It face = m_mesh.begin_cash_face();
    //OrthoMesh::Cash_Face_It endf = m_mesh.end_cash_face();
    OrthoMesh::Face_It face = m_mesh.begin_face();
    OrthoMesh::Face_It endf = m_mesh.end_face();
    bool hasbc;
    m_vbc.rewind();
    //m_solAnt receives the solution in time tn since
    //m_sol will contain the solution in tn+1. 
    m_solAnt = m_sol;

    //For each face
    for (;face!=endf; face++)
    {
      //Get the values of the previous solution in the right and left of the face
      face->getAdjCellIndices(negCell,posCell);
      hasbc=this->getValuesOfTheFaceCells(m_mesh,m_vbc,m_solAnt,face->index(),negCell,posCell,vQNeg,vQPos);

      FaceInfo fInfo(face);

      
      



      //Get the porsity
      double posPor = getPorosity()(fInfo.cell2);
      double negPor = getPorosity()(fInfo.cell1);
      //double mPor = (posPor + negPor)/2.0;

      //get Flux
      m_flux.fluxAtFace(vFlux,fInfo,vQNeg,vQPos);
      if (!hasbc)
      {
        a=fabs(m_flux.maxLocalCharVelocity(fInfo,vQNeg,vQPos));
      }

      //Now update the cell values component by compoenent
      for (unsigned cmp=0;cmp<n_components;cmp++)
      {
        //See if the face has prescribed boundary condition BC      
        double flux;
        if (!hasbc)    
        {
          /*
            The face doesnt have boundary condition BC
            proceed to the usual computation of the flux
            to calculate the amount of the mass passed through
            the face and divide it by the cell volume
          */
          assert(a*dt<=m_mesh.get_dx()/2.0);
          flux = dt*face->areaPerCellVol()*(vFlux(cmp) -a*(vQPos(cmp)-vQNeg(cmp))/2.0);
        }
        else //the face has BC, so dont add the stability term
        {
          flux = dt*face->areaPerCellVol() * vFlux(cmp);
        }

      
      
        //printf("face %d = %g,Pos: %d,Neg %d, SwPos = %g,SwNeg = %g \n",faceIndex,flux,face->getPosCellIndex(),face->getNegCellIndex(),SwPos,SwNeg);
        if (face->hasPosCell())
          m_sol(posCell,cmp) +=  + flux/posPor;

        if (face->hasNegCell())
          m_sol(negCell,cmp) +=  - flux/negPor;
        
      } //end for each component
    } //end for each face
  } //end for each time step
}

---

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

• russanovforsystem.h
• russanovforsystem.cpp