laxfriedrichsforsystem.h

Go to the documentation of this file.

#ifndef  _MY_LaxFriedrichsForSystem_
#define _MY_LaxFriedrichsForSystem_
#include "conservativemethodforsystem.h"
#include "facefluxfunction.h"
#include "diffusivestep.h"
#include "vecincbc.h"
#include "arrayofvecdouble.h"
#include "fixedvaluecondition.h"


class LaxFriedrichsForSystem : public ConservativeMethodForSystem
{
 protected:
  VecIncBC m_vbc;
  VecIndex m_cellsOfBCFaces;
  OrthoMesh  &m_mesh; 
  const VecDouble& m_cPor; 
  FaceFluxFunction &m_flux;
  const unsigned n_components;
  ArrayOfVecDouble m_sol,m_solAnt;
  FixedValueCondition m_fixedC;
  //std::vector<VecDouble> m_fBC; /*!<Array of vectors with boundary conditions of the solution in the boundary faces*/
  //  ArrayOfVecDouble m_vSol;
  //ArrayOfVecDouble m_vSolAnt;
  
  double m_CFL;


 protected:

 public:
  LaxFriedrichsForSystem(OrthoMesh &mesh,Function3D &fInitU,const VecDouble &cPor,Function3D &fPrescribedU,FaceFluxFunction &flux, FixedValueCondition &fixedC, double CFL);
  virtual ArrayOfVecDouble& getSolutionAtCells();
 
  virtual void updateVelocities(DynamicBase &dynMod);
  const VecDouble& getPorosity(){return m_cPor;} 
  virtual void iterateN(unsigned nSteps, double dt);
  virtual double getDt(double t,double tEnd);
  virtual ~LaxFriedrichsForSystem(){}
  virtual void printOutput();
  virtual void updateFixedPressureBC(Function3D *fDP,FlashCompositional &flash);
  void getAdjCellOfBCFaces(OrthoMesh &mesh,VecIncBC &vbc,VecIndex &vCells);
  void updateBCForCompressibility(VecIncBC &vbc,DynamicBase &dyn,FlashCompositional &flash); 


  //  virtual void updateValuesAtBoundary(Function3D *fDP,Function3D *fFlux,FlashCompositional &flash);
  

  
};


#endif

---