#include <EBLevelAdvect.H>

Public Member Functions
	EBLevelAdvect ()
	Default constructor. More...

	~EBLevelAdvect ()
	Destructor. More...

	EBLevelAdvect (const DisjointBoxLayout &a_thisDBL, const DisjointBoxLayout &a_coarDBL, const EBISLayout &a_thisEBISL, const EBISLayout &a_coarEBISL, const ProblemDomain &a_DProblem, const int &a_nRefine, const RealVect &a_dx, const bool &a_hasCoarser, const bool &a_hasFiner, const EBPatchGodunovFactory const a_patchGodunov, const bool &a_forceNoEBCF=false, const EBIndexSpace const a_eb=Chombo_EBIS::instance())

void	define (const DisjointBoxLayout &a_thisDBL, const DisjointBoxLayout &a_coarDBL, const EBISLayout &a_thisEBISL, const EBISLayout &a_coarEBISL, const ProblemDomain &a_DProblem, const int &a_nRefine, const RealVect &a_dx, const bool &a_hasCoarser, const bool &a_hasFiner, const EBPatchGodunovFactory const a_patchGodunov, const bool &a_forceNoEBCF=false, const EBIndexSpace const a_eb=Chombo_EBIS::instance())
	Actual constructor. More...

void	advectToFacesCol (LevelData< EBFluxFAB > &a_extrapState, LayoutData< Vector< BaseIVFAB< Real > * > > &a_coveredPrimLo, LayoutData< Vector< BaseIVFAB< Real > * > > &a_coveredPrimHi, const LayoutData< Vector< Vector< VolIndex > > > &a_coveredFaceLo, const LayoutData< Vector< Vector< VolIndex > > > &a_coveredFaceHi, const LayoutData< Vector< IntVectSet > > &a_coveredSetsLo, const LayoutData< Vector< IntVectSet > > &a_coveredSetsHi, const LevelData< EBCellFAB > &a_consState, const LevelData< EBCellFAB > &a_normalVel, const LevelData< EBFluxFAB > &a_advectionVel, const LevelData< EBCellFAB > a_consStateCoarseOld, const LevelData< EBCellFAB > a_consStateCoarseNew, const LevelData< EBCellFAB > a_normalVelCoarseOld, const LevelData< EBCellFAB > a_normalVelCoarseNew, const Real &a_timeCoarseOld, const Real &a_timeCoarseNew, const Real &a_timeFine, const Real &a_dt, const LevelData< EBCellFAB > const a_source=NULL, const LevelData< EBCellFAB > const a_sourceCoarOld=NULL, const LevelData< EBCellFAB > *const a_sourceCoarNew=NULL)

void	advectToFaces (EBFluxFAB &a_extrapState, BaseIVFAB< Real > &a_boundaryPrim, Vector< BaseIVFAB< Real > * > &a_coveredPrimLo, Vector< BaseIVFAB< Real > * > &a_coveredPrimHi, const Vector< Vector< VolIndex > > &a_coveredFaceLo, const Vector< Vector< VolIndex > > &a_coveredFaceHi, const Vector< IntVectSet > &a_coveredSetsLo, const Vector< IntVectSet > &a_coveredSetsHi, const EBCellFAB &a_consState, const EBCellFAB &a_normalVel, const EBFluxFAB &a_advectionVel, const Box &a_cellBox, const EBISBox &a_ebisBox, const Real &a_dt, const Real &a_time, const EBCellFAB &a_source, const DataIndex &a_dit, bool a_doBoundaryPrim)

void	advectToFacesBCG (LevelData< EBFluxFAB > &a_extrapState, const LevelData< EBCellFAB > &a_consState, const LevelData< EBCellFAB > &a_normalVel, const LevelData< EBFluxFAB > &a_advectionVel, const LevelData< EBCellFAB > a_consStateCoarseOld, const LevelData< EBCellFAB > a_consStateCoarseNew, const LevelData< EBCellFAB > a_normalVelCoarseOld, const LevelData< EBCellFAB > a_normalVelCoarseNew, const Real &a_timeCoarseOld, const Real &a_timeCoarseNew, const Real &a_timeFine, const Real &a_dt, const LevelData< EBCellFAB > const a_source=NULL, const LevelData< EBCellFAB > const a_sourceCoarOld=NULL, const LevelData< EBCellFAB > *const a_sourceCoarNew=NULL)

void	advectToFacesBCG (EBFluxFAB &a_extrapState, BaseIVFAB< Real > &a_boundaryPrim, const EBCellFAB &a_consState, const EBCellFAB &a_normalVel, const EBFluxFAB &a_advectionVel, const Box &a_cellBox, const EBISBox &a_ebisBox, const Real &a_dt, const Real &a_time, const EBCellFAB &a_source, const DataIndex &a_dit, bool a_doBoundaryPrim)

void	computeNormalVel (LevelData< EBCellFAB > &a_normalVel, const LevelData< EBFluxFAB > &a_advectionVel, const LayoutData< Vector< BaseIVFAB< Real > * > > &a_coveredVeloLo, const LayoutData< Vector< BaseIVFAB< Real > * > > &a_coveredVeloHi, const LayoutData< Vector< Vector< VolIndex > > > &a_coveredFaceLo, const LayoutData< Vector< Vector< VolIndex > > > &a_coveredFaceHi) const

bool	isDefined () const

EBPatchAdvect &	getPatchAdvect (const DataIndex &a_dit)

void	resetBCs (const RefCountedPtr< EBPhysIBCFactory > &a_advectBC)

Protected Attributes
bool	m_hasCoarser

bool	m_hasFiner

bool	m_forceNoEBCF

bool	m_isDefined

RealVect	m_dx

int	m_refRatCrse

int	m_nVar

int	m_nGhost

DisjointBoxLayout	m_thisGrids

DisjointBoxLayout	m_coarGrids

EBISLayout	m_thisEBISL

EBISLayout	m_coarEBISL

AggEBPWLFillPatch	m_fillPatch

AggEBPWLFillPatch	m_fillPatchVel

ProblemDomain	m_domain

LayoutData< EBPatchAdvect * >	m_ebPatchAdvect

Private Member Functions
void	operator= (const EBLevelAdvect &a_input)

	EBLevelAdvect (const EBLevelAdvect &a_input)

Detailed Description

Advection and diffusion on a level using

Constructor & Destructor Documentation

◆ EBLevelAdvect() [1/3]

EBLevelAdvect::EBLevelAdvect ( )

Default constructor.

Object requires define(..) to be called before all other functions.

◆ ~EBLevelAdvect()

EBLevelAdvect::~EBLevelAdvect ( )

Destructor.

destroys all objects created by define(..). Passed in data references of define(..) are left alone.

◆ EBLevelAdvect() [2/3]

EBLevelAdvect::EBLevelAdvect	(	const DisjointBoxLayout &	a_thisDBL,
		const DisjointBoxLayout &	a_coarDBL,
		const EBISLayout &	a_thisEBISL,
		const EBISLayout &	a_coarEBISL,
		const ProblemDomain &	a_DProblem,
		const int &	a_nRefine,
		const RealVect &	a_dx,
		const bool &	a_hasCoarser,
		const bool &	a_hasFiner,
		const EBPatchGodunovFactory *const	a_patchGodunov,
		const bool &	a_forceNoEBCF = `false`,
		const EBIndexSpace *const	a_eb = `Chombo_EBIS::instance()`
	)

forceNoEBCF = true promises that the embedded boundary will not cross the coarse-fine interface. This will save you a lot of compute time if you are certain you are tagging all irregulkar cells.

◆ EBLevelAdvect() [3/3]

EBLevelAdvect::EBLevelAdvect ( const EBLevelAdvect & a_input )

inlineprivate

References MayDay::Error().

Member Function Documentation

◆ define()

void EBLevelAdvect::define	(	const DisjointBoxLayout &	a_thisDBL,
		const DisjointBoxLayout &	a_coarDBL,
		const EBISLayout &	a_thisEBISL,
		const EBISLayout &	a_coarEBISL,
		const ProblemDomain &	a_DProblem,
		const int &	a_nRefine,
		const RealVect &	a_dx,
		const bool &	a_hasCoarser,
		const bool &	a_hasFiner,
		const EBPatchGodunovFactory *const	a_patchGodunov,
		const bool &	a_forceNoEBCF = `false`,
		const EBIndexSpace *const	a_eb = `Chombo_EBIS::instance()`
	)

Actual constructor.

For the coarsest level, an empty DisjointBoxLayout is passed in for coaserDisjointBoxLayout. Inside the routine, we cast away const-ness on the data members for the assignment. The arguments passed in are maintained const. (coding standards). a_nRefine is the refinement ratio between this level and the next coarser level. a_numGhosts is the number of ghost cells in each direction. forceNoEBCF = true promises that the embedded boundary will not cross the coarse-fine interface. This will save you a lot of compute time if you are certain you are tagging all irregulkar cells.

◆ advectToFacesCol()

void EBLevelAdvect::advectToFacesCol	(	LevelData< EBFluxFAB > &	a_extrapState,
		LayoutData< Vector< BaseIVFAB< Real > * > > &	a_coveredPrimLo,
		LayoutData< Vector< BaseIVFAB< Real > * > > &	a_coveredPrimHi,
		const LayoutData< Vector< Vector< VolIndex > > > &	a_coveredFaceLo,
		const LayoutData< Vector< Vector< VolIndex > > > &	a_coveredFaceHi,
		const LayoutData< Vector< IntVectSet > > &	a_coveredSetsLo,
		const LayoutData< Vector< IntVectSet > > &	a_coveredSetsHi,
		const LevelData< EBCellFAB > &	a_consState,
		const LevelData< EBCellFAB > &	a_normalVel,
		const LevelData< EBFluxFAB > &	a_advectionVel,
		const LevelData< EBCellFAB > *	a_consStateCoarseOld,
		const LevelData< EBCellFAB > *	a_consStateCoarseNew,
		const LevelData< EBCellFAB > *	a_normalVelCoarseOld,
		const LevelData< EBCellFAB > *	a_normalVelCoarseNew,
		const Real &	a_timeCoarseOld,
		const Real &	a_timeCoarseNew,
		const Real &	a_timeFine,
		const Real &	a_dt,
		const LevelData< EBCellFAB > *const	a_source = `NULL`,
		const LevelData< EBCellFAB > *const	a_sourceCoarOld = `NULL`,
		const LevelData< EBCellFAB > *const	a_sourceCoarNew = `NULL`
	)

Average advection velocity to cell centers to get unorm.
Do linear C/F interpolation on a_cellState and compute the diffusion source term if necessary. Quadratic is unnecessary because the laplacian does its own CF interpolation in the viscous case.
Do linear C/F interpolation in space and time
Extrapolate cellState to to faces using patchgodunov->extrapolatePrim

Notes: The advection velocity that is sent in is not used. It will be taken out in later designs

◆ advectToFaces()

void EBLevelAdvect::advectToFaces	(	EBFluxFAB &	a_extrapState,
		BaseIVFAB< Real > &	a_boundaryPrim,
		Vector< BaseIVFAB< Real > * > &	a_coveredPrimLo,
		Vector< BaseIVFAB< Real > * > &	a_coveredPrimHi,
		const Vector< Vector< VolIndex > > &	a_coveredFaceLo,
		const Vector< Vector< VolIndex > > &	a_coveredFaceHi,
		const Vector< IntVectSet > &	a_coveredSetsLo,
		const Vector< IntVectSet > &	a_coveredSetsHi,
		const EBCellFAB &	a_consState,
		const EBCellFAB &	a_normalVel,
		const EBFluxFAB &	a_advectionVel,
		const Box &	a_cellBox,
		const EBISBox &	a_ebisBox,
		const Real &	a_dt,
		const Real &	a_time,
		const EBCellFAB &	a_source,
		const DataIndex &	a_dit,
		bool	a_doBoundaryPrim
	)

◆ advectToFacesBCG() [1/2]

void EBLevelAdvect::advectToFacesBCG	(	LevelData< EBFluxFAB > &	a_extrapState,
		const LevelData< EBCellFAB > &	a_consState,
		const LevelData< EBCellFAB > &	a_normalVel,
		const LevelData< EBFluxFAB > &	a_advectionVel,
		const LevelData< EBCellFAB > *	a_consStateCoarseOld,
		const LevelData< EBCellFAB > *	a_consStateCoarseNew,
		const LevelData< EBCellFAB > *	a_normalVelCoarseOld,
		const LevelData< EBCellFAB > *	a_normalVelCoarseNew,
		const Real &	a_timeCoarseOld,
		const Real &	a_timeCoarseNew,
		const Real &	a_timeFine,
		const Real &	a_dt,
		const LevelData< EBCellFAB > *const	a_source = `NULL`,
		const LevelData< EBCellFAB > *const	a_sourceCoarOld = `NULL`,
		const LevelData< EBCellFAB > *const	a_sourceCoarNew = `NULL`
	)

◆ advectToFacesBCG() [2/2]

void EBLevelAdvect::advectToFacesBCG	(	EBFluxFAB &	a_extrapState,
		BaseIVFAB< Real > &	a_boundaryPrim,
		const EBCellFAB &	a_consState,
		const EBCellFAB &	a_normalVel,
		const EBFluxFAB &	a_advectionVel,
		const Box &	a_cellBox,
		const EBISBox &	a_ebisBox,
		const Real &	a_dt,
		const Real &	a_time,
		const EBCellFAB &	a_source,
		const DataIndex &	a_dit,
		bool	a_doBoundaryPrim
	)

◆ computeNormalVel()

void EBLevelAdvect::computeNormalVel	(	LevelData< EBCellFAB > &	a_normalVel,
		const LevelData< EBFluxFAB > &	a_advectionVel,
		const LayoutData< Vector< BaseIVFAB< Real > * > > &	a_coveredVeloLo,
		const LayoutData< Vector< BaseIVFAB< Real > * > > &	a_coveredVeloHi,
		const LayoutData< Vector< Vector< VolIndex > > > &	a_coveredFaceLo,
		const LayoutData< Vector< Vector< VolIndex > > > &	a_coveredFaceHi
	)		const

◆ isDefined()

bool EBLevelAdvect::isDefined ( ) const

◆ getPatchAdvect()

EBPatchAdvect& EBLevelAdvect::getPatchAdvect ( const DataIndex & a_dit )

inline

References m_ebPatchAdvect, and resetBCs().

◆ resetBCs()

void EBLevelAdvect::resetBCs ( const RefCountedPtr< EBPhysIBCFactory > & a_advectBC )

Referenced by getPatchAdvect().

◆ operator=()

void EBLevelAdvect::operator= ( const EBLevelAdvect & a_input )

inlineprivate

References MayDay::Error().

Member Data Documentation

◆ m_hasCoarser

bool EBLevelAdvect::m_hasCoarser

protected

◆ m_hasFiner

bool EBLevelAdvect::m_hasFiner

protected

◆ m_forceNoEBCF

bool EBLevelAdvect::m_forceNoEBCF

protected

◆ m_isDefined

bool EBLevelAdvect::m_isDefined

protected

◆ m_dx

RealVect EBLevelAdvect::m_dx

protected

◆ m_refRatCrse

int EBLevelAdvect::m_refRatCrse

protected

◆ m_nVar

int EBLevelAdvect::m_nVar

protected

◆ m_nGhost

int EBLevelAdvect::m_nGhost

protected

◆ m_thisGrids

DisjointBoxLayout EBLevelAdvect::m_thisGrids

protected

◆ m_coarGrids

DisjointBoxLayout EBLevelAdvect::m_coarGrids

protected

◆ m_thisEBISL

EBISLayout EBLevelAdvect::m_thisEBISL

protected

◆ m_coarEBISL

EBISLayout EBLevelAdvect::m_coarEBISL

protected

◆ m_fillPatch

AggEBPWLFillPatch EBLevelAdvect::m_fillPatch

protected

◆ m_fillPatchVel

AggEBPWLFillPatch EBLevelAdvect::m_fillPatchVel

protected

◆ m_domain

ProblemDomain EBLevelAdvect::m_domain

protected

◆ m_ebPatchAdvect

LayoutData<EBPatchAdvect*> EBLevelAdvect::m_ebPatchAdvect

protected

Referenced by getPatchAdvect().

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

EBLevelAdvect.H

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ EBLevelAdvect() [1/3]

◆ ~EBLevelAdvect()

◆ EBLevelAdvect() [2/3]

◆ EBLevelAdvect() [3/3]

Member Function Documentation

◆ define()

◆ advectToFacesCol()

◆ advectToFaces()

◆ advectToFacesBCG() [1/2]

◆ advectToFacesBCG() [2/2]

◆ computeNormalVel()

◆ isDefined()

◆ getPatchAdvect()

◆ resetBCs()

◆ operator=()

Member Data Documentation

◆ m_hasCoarser

◆ m_hasFiner

◆ m_forceNoEBCF

◆ m_isDefined

◆ m_dx

◆ m_refRatCrse

◆ m_nVar

◆ m_nGhost

◆ m_thisGrids

◆ m_coarGrids

◆ m_thisEBISL

◆ m_coarEBISL

◆ m_fillPatch

◆ m_fillPatchVel

◆ m_domain

◆ m_ebPatchAdvect