#include <lsfv.h>

Inheritance diagram for CouplingLsFv< nDim_, SysEqn >:

Collaboration diagram for CouplingLsFv< nDim_, SysEqn >:

Public Types
using	LsSolver = LsCartesianSolver< nDim >

using	FvCartesianSolver = FvCartesianSolverXD< nDim, SysEqn >

using	Cell = typename maia::grid::tree::Tree< nDim >::Cell

Public Types inherited from CouplingLS< nDim_ >
using	solverType = LsCartesianSolver< nDim >

Public Member Functions
	CouplingLsFv (const MInt couplingId, LsSolver ls, FvCartesianSolver fv)

void	init () override

void	preCouple (MInt) override

void	subCouple (MInt, MInt, std::vector< MBool > &) override

void	postCouple (MInt recipeStep=0) override

void	finalizeCouplerInit () override

void	finalizeSubCoupleInit (MInt) override

void	postAdaptation () override
	finalizeAdaptation More...

void	cleanUp () override

Public Member Functions inherited from CouplingLS< nDim_ >
	CouplingLS (const MInt couplingId, solverType *b)

solverType &	lsSolver () const

MInt	a_noLsCells () const

MFloat	a_outsideGValue () const

MInt	a_noG0Cells (MInt set) const

MInt	a_noBandCells (MInt set) const

MInt	a_maxGCellLevel (const MInt setId) const

MFloat &	a_levelSetFunctionG (const MInt cellId, const MInt setId)

MInt	a_bodyIdG (const MInt cellId, const MInt set) const

MInt &	a_bodyIdG (const MInt cellId, const MInt set)

MFloat	a_coordinateG (const MInt gCellId, const MInt dim) const

MInt	a_G0CellId (const MInt id, const MInt set) const

MFloat	a_normalVectorG (const MInt gCellId, const MInt dim, const MInt set) const

MFloat	a_curvatureG (const MInt gCellId, const MInt set) const

MBool	a_inBandG (MInt gcellId, MInt set) const

MInt	a_potentialGapCellClose (MInt gcellId) const

MBool	a_nearGapG (const MInt gcellId) const

MInt	a_bodyToSet (const MInt bodyId) const

MInt	a_noEmbeddedBodies () const

MInt	a_noSets () const

MInt	a_maxnoSets () const

MInt	a_startSet () const

MFloat &	a_extensionVelocityG (const MInt cellId, const MInt dim, const MInt setId)

Public Member Functions inherited from Coupling
	Coupling (const MInt couplingId)

virtual	~Coupling ()=default

	Coupling (const Coupling &)=delete

Coupling &	operator= (const Coupling &)=delete

MInt	couplerId () const

virtual void	init ()=0

virtual void	finalizeSubCoupleInit (MInt solverId)=0

virtual void	finalizeCouplerInit ()=0

virtual void	preCouple (MInt recipeStep)=0

virtual void	subCouple (MInt recipeStep, MInt solverId, std::vector< MBool > &solverCompleted)=0

virtual void	postCouple (MInt recipeStep)=0

virtual void	cleanUp ()=0

virtual void	balancePre ()
	Load balancing. More...

virtual void	balancePost ()

virtual void	reinitAfterBalance ()

virtual void	prepareAdaptation ()

virtual void	postAdaptation ()

virtual void	finalizeAdaptation (const MInt)

virtual void	writeRestartFile (const MInt)

virtual MInt	noCellDataDlb () const
	Methods to inquire coupler data during balancing. More...

virtual MInt	cellDataTypeDlb (const MInt NotUsed(dataId)) const

virtual MInt	cellDataSizeDlb (const MInt NotUsed(dataId), const MInt NotUsed(cellId))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MInt const NotUsed(data))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MLong const NotUsed(data))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MFloat const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MInt *const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MLong *const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MFloat *const NotUsed(data))

virtual void	finalizeBalance (const MInt)

virtual MInt	noCouplingTimers (const MBool NotUsed(allTimings)) const
	Number of coupling timers. More...

virtual void	getCouplingTimings (std::vector< std::pair< MString, MFloat > > &NotUsed(timings), const MBool NotUsed(allTimings))
	Return coupling timings. More...

virtual void	getDomainDecompositionInformation (std::vector< std::pair< MString, MInt > > &NotUsed(domainInfo))
	Return information on current domain decomposition (e.g. number of coupled cells/elements/...) More...

void	setDlbTimer (const MInt timerId)

void	startLoadTimer (const MString &name) const
	Start the load timer of the coupler. More...

void	stopLoadTimer (const MString &name) const
	Stop the load timer of the coupler. More...

Static Public Attributes
static constexpr MInt	nDim = nDim_

Private Member Functions
void	initData ()
	Initialize coupling-class-specific Data. More...

void	checkProperties () override
	Checks property-data which is read in by both ls-and Fv-Solver. More...

void	readProperties () override

void	returnStep_semiLagrange ()
	mimics the behaviour of the rungeKuttaStep() methods with respect to increasing time More...

void	transferGapCellProperty ()
	Sets the gapcell-property. More...

void	testGapProperty ()
	transfers the LevelSetValues from the levelset to the moving boundary Part More...

void	computeGCellTimeStep ()
	computes the gcell time step More...

void	computeBodyProperties (MInt returnMode, MFloat *bodyData, MInt body, MFloat time)

void	testCoupling ()
	transfers the LevelSetValues from the levelset to the moving boundary Part More...

void	setLsInList (MIntScratchSpace &)

void	transferLevelSetValues ()
	Sets the Levelset-Values in fvSolver. More...

MInt	ls2fvId (const MInt lsId)

MInt	ls2fvIdParent (const MInt lsId)

MInt	fv2lsId (const MInt fvId)

MInt	fv2lsIdParent (const MInt fvId)

MInt	noLevelSetFieldData ()

MFloat	interpolateLevelSet (MInt cellId, MFloat *point, MInt set)

MBool	returnStep ()
	mimics the behaviour of the rungeKuttaStep() methods with respect to increasing time More...

MFloat	lsTimeStep () const

MFloat	crankAngle (MFloat)
	help-function for engine calculations which returns the crank-angle for a given time More...

MFloat	a_meanCoord (MInt dir) const

MFloat	a_UInfinity () const

MFloat	a_TInfinity () const

MFloat	a_Ma () const

MFloat	a_time () const

MFloat	a_timeRef () const

MFloat	a_physicalTime () const

MInt	a_noFvCells () const

MInt	a_noFvGridCells () const

MInt	a_G0CellId (MInt id, MInt set) const

MInt	a_RKStep () const

MInt	a_noRKSteps () const

Private Member Functions inherited from CouplingFv< nDim_, SysEqn >
	CouplingFv (const MInt couplingId, std::vector< FvCartesianSolverXD< nDim, SysEqn > * > fvSolvers, const MInt noSolvers)

	CouplingFv (const MInt couplingId, Solver *solvers)

	CouplingFv (const CouplingFv &)=delete

	~CouplingFv () override=default

CouplingFv &	operator= (const CouplingFv &)=delete

MInt	noSolvers () const

solverType &	fvSolver (const MInt solverId=0) const

MInt	a_noFvCells () const

MInt	a_noFvGridCells () const

Private Member Functions inherited from Coupling
	Coupling (const MInt couplingId)

virtual	~Coupling ()=default

	Coupling (const Coupling &)=delete

Coupling &	operator= (const Coupling &)=delete

MInt	couplerId () const

virtual void	init ()=0

virtual void	finalizeSubCoupleInit (MInt solverId)=0

virtual void	finalizeCouplerInit ()=0

virtual void	preCouple (MInt recipeStep)=0

virtual void	subCouple (MInt recipeStep, MInt solverId, std::vector< MBool > &solverCompleted)=0

virtual void	postCouple (MInt recipeStep)=0

virtual void	cleanUp ()=0

virtual void	balancePre ()
	Load balancing. More...

virtual void	balancePost ()

virtual void	reinitAfterBalance ()

virtual void	prepareAdaptation ()

virtual void	postAdaptation ()

virtual void	finalizeAdaptation (const MInt)

virtual void	writeRestartFile (const MInt)

virtual MInt	noCellDataDlb () const
	Methods to inquire coupler data during balancing. More...

virtual MInt	cellDataTypeDlb (const MInt NotUsed(dataId)) const

virtual MInt	cellDataSizeDlb (const MInt NotUsed(dataId), const MInt NotUsed(cellId))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MInt const NotUsed(data))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MLong const NotUsed(data))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MFloat const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MInt *const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MLong *const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MFloat *const NotUsed(data))

virtual void	finalizeBalance (const MInt)

virtual MInt	noCouplingTimers (const MBool NotUsed(allTimings)) const
	Number of coupling timers. More...

virtual void	getCouplingTimings (std::vector< std::pair< MString, MFloat > > &NotUsed(timings), const MBool NotUsed(allTimings))
	Return coupling timings. More...

virtual void	getDomainDecompositionInformation (std::vector< std::pair< MString, MInt > > &NotUsed(domainInfo))
	Return information on current domain decomposition (e.g. number of coupled cells/elements/...) More...

void	setDlbTimer (const MInt timerId)

void	startLoadTimer (const MString &name) const
	Start the load timer of the coupler. More...

void	stopLoadTimer (const MString &name) const
	Stop the load timer of the coupler. More...

MFloat	returnLoadRecord () const

MFloat	returnIdleRecord () const

Private Attributes
MInt	m_bandWidthRef

MInt	m_bandWidthRefMax

MInt	m_fvSolverId

MInt	m_lsSolverId

MInt	m_G0regionId

MFloat	m_initialCrankAngle

MInt *	m_hadGapCells {}

MInt *	m_hasGapCells {}

MFloat	m_cfl

MFloat	m_maxVelocity

MInt	m_timeStepMethod

MString	m_solverMethod

MBool	m_static_computeBodyProperties_first = true

MFloat	m_static_computeBodyProperties_amplitude [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_freqFactor [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_initialBodyCenter [m_maxNoEmbeddedBodies *3] {}

MFloat	m_static_computeBodyProperties_Strouhal {}

MFloat	m_static_computeBodyProperties_mu [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_mu2 [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_liftStartAngle1 [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_liftEndAngle1 [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_liftStartAngle2 [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_liftEndAngle2 [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_circleStartAngle [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_normal [m_maxNoEmbeddedBodies *3] {}

MInt	m_static_computeBodyProperties_bodyToFunction [m_maxNoEmbeddedBodies] {}

MFloat	m_static_computeBodyProperties_omega {}

MFloat	m_static_computeBodyProperties_rotAngle {}

MBool	m_static_crankAngle_first = true

MFloat	m_static_crankAngle_Strouhal

Private Attributes inherited from CouplingFv< nDim_, SysEqn >
std::vector< solverType * >	m_fvSolvers

Static Private Attributes
static constexpr const MInt	m_maxNoEmbeddedBodies = 20

Friends
class	FvCartesianSolverXD< nDim, SysEqn >

class	LsCartesianSolver< nDim >

Additional Inherited Members
Protected Member Functions inherited from Coupling
MFloat	returnLoadRecord () const

MFloat	returnIdleRecord () const

Private Types inherited from CouplingFv< nDim_, SysEqn >
using	solverType = FvCartesianSolverXD< nDim, SysEqn >

Detailed Description

template<MInt nDim_, class SysEqn>
class CouplingLsFv< nDim_, SysEqn >

Definition at line 27 of file lsfv.h.

Member Typedef Documentation

◆ Cell

template<MInt nDim_, class SysEqn >

using CouplingLsFv< nDim_, SysEqn >::Cell = typename maia::grid::tree::Tree<nDim>::Cell

Definition at line 43 of file lsfv.h.

◆ FvCartesianSolver

template<MInt nDim_, class SysEqn >

using CouplingLsFv< nDim_, SysEqn >::FvCartesianSolver = FvCartesianSolverXD<nDim, SysEqn>

Definition at line 39 of file lsfv.h.

◆ LsSolver

template<MInt nDim_, class SysEqn >

using CouplingLsFv< nDim_, SysEqn >::LsSolver = LsCartesianSolver<nDim>

Definition at line 37 of file lsfv.h.

Constructor & Destructor Documentation

◆ CouplingLsFv()

template<MInt nDim, class SysEqn >

CouplingLsFv< nDim, SysEqn >::CouplingLsFv	(	const MInt	couplingId,
		LsSolver *	ls,
		FvCartesianSolver *	fv
	)

Definition at line 26 of file lsfv.cpp.

  : Coupling(couplingId), CouplingLS<nDim>(couplingId, ls), CouplingFv<nDim, SysEqn>(couplingId, fv) {
  TRACE();
 
  initData();
  readProperties();
  checkProperties();
 
  for(MInt dir = 0; dir < nDim; dir++) {
    lsSolver().a_meanCoord(dir) = a_meanCoord(dir);
  }
}

Member Function Documentation

◆ a_G0CellId()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::a_G0CellId	(	MInt	id,
		MInt	set
	)		const

inlineprivate

Definition at line 167 of file lsfv.h.

167{ return lsSolver().a_G0CellId(id, set); }

◆ a_Ma()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::a_Ma ( ) const

inlineprivate

Definition at line 160 of file lsfv.h.

160{ return fvSolver().a_Ma(); }

CouplingFv< nDim_, SysEqn >::fvSolver

solverType & fvSolver(const MInt solverId=0) const

Definition: coupling.h:386

FvCartesianSolverXD::a_Ma

MFloat & a_Ma()

Returns the Mach number of the solver.

Definition: fvcartesiansolverxd.h:861

◆ a_meanCoord()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::a_meanCoord ( MInt dir ) const

inlineprivate

Definition at line 156 of file lsfv.h.

156{ return fvSolver().m_meanCoord[dir]; }

FvCartesianSolverXD::m_meanCoord

MFloat m_meanCoord[3]

Definition: fvcartesiansolverxd.h:2980

◆ a_noFvCells()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::a_noFvCells ( ) const

inlineprivate

Definition at line 165 of file lsfv.h.

165{ return fvSolver().a_noCells(); }

FvCartesianSolverXD::a_noCells

MInt a_noCells() const

Returns the number of cells.

Definition: fvcartesiansolverxd.h:336

◆ a_noFvGridCells()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::a_noFvGridCells ( ) const

inlineprivate

Definition at line 166 of file lsfv.h.

166{ return fvSolver().c_noCells(); }

maia::CartesianSolver::c_noCells

MInt c_noCells() const

Definition: cartesiansolver.h:164

◆ a_noRKSteps()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::a_noRKSteps ( ) const

inlineprivate

Definition at line 169 of file lsfv.h.

169{ return fvSolver().m_noRKSteps; }

FvCartesianSolverXD::m_noRKSteps

MInt m_noRKSteps

Definition: fvcartesiansolverxd.h:1710

◆ a_physicalTime()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::a_physicalTime ( ) const

inlineprivate

Definition at line 164 of file lsfv.h.

164{ return fvSolver().m_physicalTime; }

FvCartesianSolverXD::m_physicalTime

MFloat m_physicalTime

Definition: fvcartesiansolverxd.h:1993

◆ a_RKStep()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::a_RKStep ( ) const

inlineprivate

Definition at line 168 of file lsfv.h.

168{ return fvSolver().m_RKStep; }

FvCartesianSolverXD::m_RKStep

MInt m_RKStep

Definition: fvcartesiansolverxd.h:1717

◆ a_time()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::a_time ( ) const

inlineprivate

Definition at line 162 of file lsfv.h.

162{ return fvSolver().m_time; }

FvCartesianSolverXD::m_time

MFloat m_time

Definition: fvcartesiansolverxd.h:2029

◆ a_timeRef()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::a_timeRef ( ) const

inlineprivate

Definition at line 163 of file lsfv.h.

163{ return fvSolver().m_timeRef; }

FvCartesianSolverXD::m_timeRef

MFloat m_timeRef

Definition: fvcartesiansolverxd.h:2030

◆ a_TInfinity()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::a_TInfinity ( ) const

inlineprivate

Definition at line 159 of file lsfv.h.

159{ return fvSolver().m_TInfinity; }

FvCartesianSolverXD::m_TInfinity

MFloat m_TInfinity

Definition: fvcartesiansolverxd.h:2462

◆ a_UInfinity()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::a_UInfinity ( ) const

inlineprivate

Definition at line 158 of file lsfv.h.

158{ return fvSolver().m_UInfinity; }

FvCartesianSolverXD::m_UInfinity

MFloat m_UInfinity

Definition: fvcartesiansolverxd.h:2458

◆ checkProperties()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::checkProperties

overrideprivatevirtual

Author: Tim Wegmann

Implements Coupling.

Definition at line 137 of file lsfv.cpp.

                                                 {
  TRACE();
 
  ASSERT(lsSolver().a_maxGCellLevel() == fvSolver().maxRefinementLevel(), "");
}

◆ cleanUp()

template<MInt nDim_, class SysEqn >

void CouplingLsFv< nDim_, SysEqn >::cleanUp ( )

inlineoverridevirtual

Implements Coupling.

Definition at line 119 of file lsfv.h.

119{};

◆ computeBodyProperties()

template<MInt nDim_, class SysEqn >

void CouplingLsFv< nDim_, SysEqn >::computeBodyProperties	(	MInt	returnMode,
		MFloat *	bodyData,
		MInt	body,
		MFloat	time
	)

private

◆ computeGCellTimeStep()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::computeGCellTimeStep

private

Author: Daniel Hartmann, Tim Wegmann

Definition at line 214 of file lsfv.cpp.

                                                      {
  TRACE();
 
  if(m_timeStepMethod < 6 || m_timeStepMethod > 8) {
    if(string2enum(m_solverMethod) == MAIA_SEMI_LAGRANGE_LEVELSET
       || string2enum(m_solverMethod) == MAIA_RUNGE_KUTTA_LEVELSET) {
      mTerm(1, AT_, "Computation with pure LVS not possible with timeStepMethod other than 6 ,7 or 8! Please check!");
    }
  }
 
  ASSERT(m_cfl > 0, "Couldn't read cfl property for the ls-solver timestepping!");
 
  if(m_timeStepMethod == 8) {
    lsSolver().m_timeStep = m_cfl * lsSolver().m_gCellDistance;
 
  } else if(m_timeStepMethod == 6) {
    lsSolver().m_timeStep = m_cfl * lsSolver().m_gCellDistance;
    fvSolver().forceTimeStep(lsSolver().m_timeStep / a_timeRef());
 
  } else if(m_timeStepMethod == 7) {
    lsSolver().m_timeStep = m_cfl * lsSolver().m_gCellDistance / m_maxVelocity;
    fvSolver().forceTimeStep(lsSolver().m_timeStep / a_timeRef());
 
  } else {
    lsSolver().m_timeStep = lsTimeStep();
  }
}

◆ crankAngle()

template<MInt nDim, class SysEqn >

MFloat CouplingLsFv< nDim, SysEqn >::crankAngle ( MFloat elapsedTime )

private

Author: Tim Wegmann

Definition at line 1886 of file lsfv.cpp.

                                                                {
  TRACE();
 
  MFloat& Strouhal = m_static_crankAngle_Strouhal;
  MBool& first = m_static_crankAngle_first;
 
  if(first) {
    Strouhal = Context::getSolverProperty<MFloat>("Strouhal", m_lsSolverId, AT_, &Strouhal);
    first = false;
  }
 
  const MFloat mu2 = Strouhal * F2 * PI;
  MFloat cad = mu2 * elapsedTime;
  const MInt maxNoCycles = 20;
 
  for(MInt cycle = maxNoCycles; cycle > 0; cycle--) {
    if(cad >= 4 * PI * cycle) cad = cad - 4 * PI * cycle;
  }
  cad = cad * 180 / PI;
 
  // consider initial crank angle
  cad = cad + m_initialCrankAngle;
 
  return cad;
}

◆ finalizeCouplerInit()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::finalizeCouplerInit

overridevirtual

\initializes the coupler in context of the unified run loop

Author: Thomas Lürkens

Implements Coupling.

Definition at line 59 of file lsfv.cpp.

                                                     {
  for(MInt s = 0; s < lsSolver().m_noSets; s++) {
    fvSolver().m_levelSetValues[s].resize(fvSolver().a_noCells());
    fvSolver().m_curvatureG[s].resize(fvSolver().a_noCells());
  }
 
  transferGapCellProperty();
  transferLevelSetValues();
  computeGCellTimeStep();
}

◆ finalizeSubCoupleInit()

template<MInt nDim_, class SysEqn >

void CouplingLsFv< nDim_, SysEqn >::finalizeSubCoupleInit ( MInt )

inlineoverridevirtual

Implements Coupling.

Definition at line 117 of file lsfv.h.

117{};

◆ fv2lsId()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::fv2lsId ( const MInt fvId )

inlineprivate

Definition at line 134 of file lsfv.h.

                                {
    // TODO labels:COUPLER,FVMB Is this check needed? This functions should never be called for these modes i guess..
    if(fvSolver().m_levelSetMb && fvSolver().m_constructGField) {
      return -1;
    }
    return convertId(fvSolver(), lsSolver(), fvId);
  };

◆ fv2lsIdParent()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::fv2lsIdParent ( const MInt fvId )

inlineprivate

Definition at line 141 of file lsfv.h.

                                      {
    // TODO labels:COUPLER,FVMB Is this check needed? This functions should never be called for these modes i guess..
    if(fvSolver().m_levelSetMb && fvSolver().m_constructGField) {
      return -1;
    }
    return convertIdParent(fvSolver(), lsSolver(), fvId);
  };

◆ init()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::init

overridevirtual

\initializes the coupler in context of the unified run loop

Author: Thomas Lürkens

Implements Coupling.

Definition at line 43 of file lsfv.cpp.

                                      {
  mAlloc(fvSolver().m_levelSetValues, lsSolver().m_noSets, "fvSolver().m_levelSetValues", AT_);
  mAlloc(fvSolver().m_curvatureG, lsSolver().m_noSets, "fvSolver().m_curvatureG", AT_);
 
  // set the time in the lsSolver:
  lsSolver().m_time = -99;
  lsSolver().m_time = a_time();
  // transferGapCellProperty();
  // computeGCellTimeStep();
  return;
}

◆ initData()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::initData

private

Author: Tim Wegmann

Definition at line 120 of file lsfv.cpp.

                                          {
  TRACE();
 
  // solver-specific data:
  m_fvSolverId = fvSolver().m_solverId;
  m_lsSolverId = lsSolver().m_solverId;
 
  if(lsSolver().m_closeGaps) {
    mAlloc(m_hadGapCells, lsSolver().m_noGapRegions, "m_hadGapCells", 0, AT_);
    mAlloc(m_hasGapCells, lsSolver().m_noGapRegions, "m_hasGapCells", 0, AT_);
  }
}

◆ interpolateLevelSet()

template<MInt nDim, class SysEqn >

MFloat CouplingLsFv< nDim, SysEqn >::interpolateLevelSet	(	MInt	cellId,
		MFloat *	point,
		MInt	set
	)

private

Definition at line 268 of file lsfv.cpp.

                                                                                           {
  MInt interpolationCells[8] = {0, 0, 0, 0, 0, 0, 0, 0};
  MInt position = 0;
 
  position = lsSolver().setUpLevelSetInterpolationStencil(cellId, interpolationCells, point);
 
  // Interpolate level set
  if(position > -1) {
    return lsSolver().interpolateLevelSet(interpolationCells, point, set);
  } else {
    return lsSolver().a_levelSetFunctionG(cellId, set);
  }
}

◆ ls2fvId()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::ls2fvId ( const MInt lsId )

inlineprivate

Definition at line 132 of file lsfv.h.

132{ return convertId(lsSolver(), fvSolver(), lsId); };

◆ ls2fvIdParent()

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::ls2fvIdParent ( const MInt lsId )

inlineprivate

Definition at line 133 of file lsfv.h.

133{ return convertIdParent(lsSolver(), fvSolver(), lsId); };

◆ lsTimeStep()

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::lsTimeStep ( ) const

inlineprivate

Definition at line 153 of file lsfv.h.

153{ return fvSolver().m_timeStep; }

◆ noLevelSetFieldData()

template<MInt nDim, class SysEqn >

MInt CouplingLsFv< nDim, SysEqn >::noLevelSetFieldData

private

Definition at line 284 of file lsfv.cpp.

                                                     {
  MInt noLevelSetFieldData = 0;
  ASSERT(fvSolver().m_levelSet, "");
 
  if(lsSolver().m_writeOutAllLevelSetFunctions) {
    noLevelSetFieldData += noLevelSetFieldData + 2 * a_noSets();
  } else {
    noLevelSetFieldData += 1;
  }
  if(lsSolver().m_writeOutAllCurvatures) {
    noLevelSetFieldData += a_noSets();
  } else {
    noLevelSetFieldData += 1;
  }
 
  return noLevelSetFieldData;
}

◆ postAdaptation()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::postAdaptation

overridevirtual

Author: Jannik Borgelt

Reimplemented from Coupling.

Definition at line 102 of file lsfv.cpp.

                                                {
  TRACE();
 
  for(MInt s = 0; s < lsSolver().m_noSets; s++) {
    fvSolver().m_levelSetValues[s].clear();
    fvSolver().m_levelSetValues[s].resize(fvSolver().a_noCells());
    fvSolver().m_curvatureG[s].clear();
    fvSolver().m_curvatureG[s].resize(fvSolver().a_noCells());
  }
 
  transferLevelSetValues();
}

◆ postCouple()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::postCouple ( MInt recipeStep = 0 )

overridevirtual

\ postCouple routine in context of the unified run loop

Author: Thomas Lürkens

Implements Coupling.

Definition at line 86 of file lsfv.cpp.

                                                  {
  lsSolver().m_time = a_time();
 
  transferLevelSetValues();
  transferGapCellProperty();
  testCoupling();
 
  if(!lsSolver().m_combustion && !lsSolver().m_levelSetMb && !lsSolver().m_semiLagrange && !lsSolver().m_LSSolver) {
    returnStep();
  }
}

◆ preCouple()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::preCouple ( MInt )

overridevirtual

\ preCouple routine in context of the unified run loop

Author: Jannik Borgelt

Implements Coupling.

Definition at line 74 of file lsfv.cpp.

                                                 {
  if(!lsSolver().m_combustion && !lsSolver().m_levelSetMb) {
    if(lsSolver().m_semiLagrange) {
      returnStep_semiLagrange();
    }
  }
}

◆ readProperties()

template<MInt nDim_, class SysEqn >

void CouplingLsFv< nDim_, SysEqn >::readProperties ( )

overrideprivatevirtual

Implements Coupling.

◆ returnStep()

template<MInt nDim, class SysEqn >

MBool CouplingLsFv< nDim, SysEqn >::returnStep

private

Author: Tim Wegmann

Definition at line 246 of file lsfv.cpp.

                                             {
  TRACE();
 
  fvSolver().m_time += lsTimeStep();
  fvSolver().m_physicalTime += lsTimeStep() * a_timeRef();
 
  return true;
}

◆ returnStep_semiLagrange()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::returnStep_semiLagrange

private

Author: Tim Wegmann

Definition at line 259 of file lsfv.cpp.

                                                         {
  TRACE();
 
  fvSolver().m_time += lsTimeStep() * a_timeRef();
  fvSolver().m_physicalTime += lsTimeStep() * a_timeRef();
}

◆ setLsInList()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::setLsInList ( MIntScratchSpace & )

private

Definition at line 304 of file lsfv.cpp.

                                                              {
  TRACE();
 
  cerr0 << "Setting Sensors for Fv-Solver adaptation!" << endl;
 
  MInt bandWidthRef = m_bandWidthRef;
  MInt bandWidthRefmax = m_bandWidthRefMax;
 
  MIntScratchSpace sendBufferSize(fvSolver().grid().noNeighborDomains(), AT_, "sendBufferSize");
  MIntScratchSpace receiveBufferSize(fvSolver().grid().noNeighborDomains(), AT_, "receiveBufferSize");
  MInt listCount = 0;
  MIntScratchSpace inList(a_noFvGridCells(), AT_, "inList");
  for(MInt c = 0; c < a_noFvGridCells(); c++)
    inList[c] = 0;
 
  // 2) add the g0-Cells and all parents to the list:
  MInt endSet = a_noSets();
  if(lsSolver().m_buildCollectedLevelSetFunction) {
    endSet = 1;
  }
 
  for(MInt set = 0; set < endSet; set++) {
    for(MInt id = 0; id < lsSolver().a_noG0Cells(set); id++) {
      MInt fvCellId = ls2fvIdParent(a_G0CellId(id, set));
      if(fvCellId < 0) continue;
      if(fvSolver().a_isHalo(fvCellId)) continue;
      inList[fvCellId] = 1;
      listCount++;
      MInt parentId = fvSolver().c_parentId(fvCellId);
      while(parentId > -1) {
        if(parentId < a_noFvGridCells()) inList[parentId] = 1;
        parentId = fvSolver().c_parentId(parentId);
      }
    }
  }
 
  // Exchange the listCount on all Domains
  MPI_Allreduce(&listCount, &listCount, 1, MPI_INT, MPI_SUM, fvSolver().mpiComm(), AT_, "listCount", "listCount");
 
  if(listCount == 0) {
    if(fvSolver().domainId() == 0) cerr << "No G0-Cells found!" << endl;
  }
 
  fvSolver().exchangeData(&inList[0], 1);
 
  for(MInt level = fvSolver().minLevel(); level < fvSolver().maxRefinementLevel(); level++) {
    if(level == fvSolver().maxRefinementLevel() - 1) bandWidthRef = bandWidthRefmax;
 
    // 4) Loop over the number of revinement-grid-cells and add those to the list:
    for(MInt loopMarker = 1; loopMarker < bandWidthRef; loopMarker++) {
      for(MInt cellId = 0; cellId < a_noFvGridCells(); cellId++) {
        if(fvSolver().a_level(cellId) != level) continue;
        if(inList[cellId] != loopMarker) continue;
 
        // direct neighbor
        for(MInt n = 0; n < fvSolver().m_noDirs; n++) {
          MInt nghbrId = fvSolver().c_neighborId(cellId, n, false);
          if(nghbrId < 0) continue;
          if(inList[nghbrId] == 0) inList[nghbrId] = loopMarker + 1;
 
          // diagonal neighbors
          if(n == 0 || n == 1) {
            for(MInt nn = 2; nn < fvSolver().m_noDirs; nn++) {
              MInt nghbrId2 = fvSolver().c_neighborId(nghbrId, nn, false);
              if(nghbrId2 < 0) continue;
              if(inList[nghbrId2] == 0) inList[nghbrId2] = loopMarker + 1;
            }
          }
          if(n == 4 || n == 5) {
            for(MInt nn = 2; nn < 4; nn++) {
              MInt nghbrId2 = fvSolver().c_neighborId(nghbrId, nn, false);
              if(nghbrId2 < 0) continue;
              if(inList[nghbrId2] == 0) inList[nghbrId2] = loopMarker + 1;
 
              for(MInt nnn = 0; nnn < 1; nnn++) {
                MInt nghbrId3 = fvSolver().c_neighborId(nghbrId2, nnn, false);
                if(nghbrId3 < 0) continue;
                if(inList[nghbrId3] == 0) inList[nghbrId3] = loopMarker + 1;
              }
            }
          }
        }
      }
 
      fvSolver().exchangeData(&inList[0], 1);
    }
  }
}

◆ subCouple()

template<MInt nDim_, class SysEqn >

void CouplingLsFv< nDim_, SysEqn >::subCouple	(	MInt	,
		MInt	,
		std::vector< MBool > &
	)

inlineoverridevirtual

Implements Coupling.

Definition at line 114 of file lsfv.h.

114{};

◆ testCoupling()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::testCoupling

private

Author: Tim Wegmann

Definition at line 583 of file lsfv.cpp.

                                              {
  TRACE();
 
  if(lsSolver().m_constructGField) return;
 
  for(MInt fc = 0; fc < a_noFvGridCells(); fc++) {
    ASSERT(fvSolver().grid().tree().solver2grid(fc) >= 0, "");
    ASSERT(fvSolver().grid().solverFlag(fvSolver().grid().tree().solver2grid(fc), m_fvSolverId), "");
#ifdef COUPLING_DEBUG_
    if(lsSolver().grid().solverFlag(fvSolver().grid().tree().solver2grid(fc), m_lsSolverId)) {
      ASSERT(ls2fvId(fv2lsId(fc)) == fc,
             to_string(fc) + " " + to_string(fv2lsId(fc)) + " " + to_string(ls2fvId(fv2lsId(fc))));
    }
#endif
  }
  for(MInt cellId = 0; cellId < a_noLsCells(); cellId++) {
    ASSERT(lsSolver().grid().tree().solver2grid(cellId) >= 0, "");
    ASSERT(lsSolver().grid().solverFlag(lsSolver().grid().tree().solver2grid(cellId), m_lsSolverId), "");
#ifdef COUPLING_DEBUG_
    if(fvSolver().grid().solverFlag(lsSolver().grid().tree().solver2grid(cellId), m_fvSolverId)) {
      ASSERT(fv2lsId(ls2fvId(cellId)) == cellId,
             to_string(cellId) + " " + to_string(ls2fvId(cellId)) + " " + to_string(fv2lsId(ls2fvId(cellId))));
    }
#endif
  }
 
  if(!lsSolver().m_levelSetMb) return;
 
 
#ifdef COUPLING_DEBUG_
  for(MInt fc = 0; fc < a_noFvGridCells(); fc++) {
    for(MInt dir = 0; dir < nDim; dir++) {
      if(lsSolver().grid().solverFlag(fvSolver().grid().tree().solver2grid(fc), m_lsSolverId)) {
        ASSERT(abs(fvSolver().c_coordinate(fc, dir) - a_coordinateG(fv2lsId(fc), dir)) < 0.00000001, "");
      }
    }
  }
#endif
}

◆ testGapProperty()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::testGapProperty

private

Author: Tim Wegmann

Definition at line 560 of file lsfv.cpp.

                                                 {
  TRACE();
 
  if(!lsSolver().m_closeGaps) return;
  if((lsSolver().m_levelSetMb) && fvSolver().m_constructGField) return;
 
#ifdef COUPLING_DEBUG_
  for(MInt fc = 0; fc < a_noFvGridCells(); fc++) {
    if(fvSolver().a_isGapCell(fc) || fvSolver().a_wasGapCell(fc)) {
      MInt gc = lsSolver().grid().tree().grid2solver(fvSolver().grid().tree().solver2grid(fc));
      ASSERT(gc > -1, "");
      ASSERT(lsSolver().a_potentialGapCellClose(gc) > 0
                 && lsSolver().a_potentialGapCellClose(gc) <= fvSolver().m_noEmbeddedBodies,
             "");
    }
  }
#endif
}

◆ transferGapCellProperty()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::transferGapCellProperty

private

Author: Tim Wegmann

Definition at line 433 of file lsfv.cpp.

                                                         {
  TRACE();
 
  if(!lsSolver().m_closeGaps) return;
  if(lsSolver().m_noGapRegions <= 0) return;
  if((lsSolver().m_levelSetMb) && fvSolver().m_constructGField) return;
 
#ifdef COUPLING_DEBUG_
  testCoupling();
#endif
 
  // 1) reset has/had Gap-Cells
  for(MInt region = 0; region < lsSolver().m_noGapRegions; region++) {
    m_hadGapCells[region] = 0;
    m_hasGapCells[region] = 0;
  }
 
  // 2) set a_wasGapCell, m_hadGapCells and initialize a_isGapCell
  for(MInt fc = 0; fc < a_noFvGridCells(); fc++) {
    if(fvSolver().a_isGapCell(fc)) {
      MInt lsId = fv2lsId(fc);
      ASSERT(lsId > -1, "");
      // skip G0-regions
      if(m_G0regionId > -1 && a_potentialGapCellClose(lsId) == m_G0regionId) continue;
      MInt regionId = a_potentialGapCellClose(lsId) - 2;
      ASSERT(a_potentialGapCellClose(lsId) > 0 && a_potentialGapCellClose(lsId) <= lsSolver().m_noEmbeddedBodies, "");
      if(regionId > -1 && regionId < lsSolver().m_noGapRegions && globalTimeStep > 0) {
        // don't set hadGapCells during initialisation, so that initGapClosure is still called
        // at the first timeStep!
        m_hadGapCells[regionId]++;
      }
    }
    if(globalTimeStep > 0) {
      fvSolver().a_wasGapCell(fc) = fvSolver().a_isGapCell(fc);
    } else {
      fvSolver().a_wasGapCell(fc) = false;
    }
    fvSolver().a_isGapCell(fc) = false;
  }
 
  for(MInt fc = a_noFvGridCells(); fc < a_noFvCells(); fc++) {
    if(globalTimeStep > 0) {
      fvSolver().a_wasGapCell(fc) = fvSolver().a_isGapCell(fc);
    } else {
      fvSolver().a_wasGapCell(fc) = false;
    }
    fvSolver().a_isGapCell(fc) = false;
  }
 
  // 4) set a_isGapCell and hasGapCells
  for(MInt gCellId = 0; gCellId < lsSolver().a_noCells(); gCellId++) {
    if(a_nearGapG(gCellId)) {
      MInt fvId = ls2fvId(gCellId);
      if(fvId < 0) {
        ASSERT(abs(abs(a_levelSetFunctionG(gCellId, 0)) - a_outsideGValue()) < 0.000001,
               "ERROR, no fv-Cell found for relevant Gap-Cells!");
      }
      if(a_potentialGapCellClose(gCellId) > 0 && a_potentialGapCellClose(gCellId) <= lsSolver().m_noEmbeddedBodies) {
        fvSolver().a_isGapCell(fvId) = true;
        MInt regionId = a_potentialGapCellClose(gCellId) - 2;
        if(regionId > -1 && regionId < lsSolver().m_noGapRegions) {
          m_hasGapCells[regionId]++;
        }
      }
    }
  }
  testGapProperty();
 
  // set isGapCell for g0regionId:
  if(m_G0regionId > -1) {
    MInt noG0regionCells = 0;
    for(MInt gCellId = 0; gCellId < lsSolver().a_noCells(); gCellId++) {
      if(lsSolver().a_potentialGapCellClose(gCellId) == m_G0regionId && lsSolver().a_gapWidth(gCellId) > 0) {
        MInt fvId = ls2fvId(gCellId);
        if(fvId < 0) {
          ASSERT(abs(abs(a_levelSetFunctionG(gCellId, 0)) - a_outsideGValue()) < 0.000001,
                 "ERROR, no fv-Cell found for relevant Gap-Cells!");
        }
        noG0regionCells++;
        fvSolver().a_isGapCell(fvId) = true;
        ASSERT(!lsSolver().a_nearGapG(gCellId), "");
      }
    }
#if defined COUPLING_DEBUG_ || !defined NDEBUG
    MPI_Allreduce(MPI_IN_PLACE, &noG0regionCells, 1, MPI_INT, MPI_SUM, fvSolver().mpiComm(), AT_, "INPLACE",
                  "noG0regionCells");
    if(lsSolver().domainId() == 0) {
      cerr << " No of g0-region Cells " << noG0regionCells << endl;
    }
#endif
  }
 
 
  // at the beginning of the restart
  if(lsSolver().m_restart && globalTimeStep == fvSolver().m_restartTimeStep) {
    for(MInt fc = a_noFvGridCells(); fc < a_noFvCells(); fc++) {
      fvSolver().a_wasGapCell(fc) = fvSolver().a_isGapCell(fc);
    }
    for(MInt region = 0; region < lsSolver().m_noGapRegions; region++) {
      m_hadGapCells[region] = m_hasGapCells[region];
    }
  }
 
  // 3) exchange hadGapCells and hasGapCells for all regions
  MPI_Allreduce(MPI_IN_PLACE, &m_hadGapCells[0], lsSolver().m_noGapRegions, MPI_INT, MPI_SUM, fvSolver().mpiComm(), AT_,
                "INPLACE", "m_hadGapCells");
 
  MPI_Allreduce(MPI_IN_PLACE, &m_hasGapCells[0], lsSolver().m_noGapRegions, MPI_INT, MPI_SUM, fvSolver().mpiComm(), AT_,
                "INPLACE", "m_hasGapCells");
 
  fvSolver().exchangeGapInfo();
 
#if defined COUPLING_DEBUG_ || !defined NDEBUG
  for(MInt region = 0; region < lsSolver().m_noGapRegions; region++) {
    if(lsSolver().domainId() == 0) {
      cerr << globalTimeStep << " region " << region << " has " << m_hasGapCells[region] << " had "
           << m_hadGapCells[region] << endl;
    }
  }
#endif
}

◆ transferLevelSetValues()

template<MInt nDim, class SysEqn >

void CouplingLsFv< nDim, SysEqn >::transferLevelSetValues

private

Author: Jannik Borgelt

Definition at line 398 of file lsfv.cpp.

                                                        {
  TRACE();
 
  fvSolver().a_noSets() = a_noSets();
  fvSolver().a_noLevelSetFieldData() = noLevelSetFieldData();
 
  for(MInt s = 0; s < lsSolver().m_noSets; s++) {
    fvSolver().m_levelSetValues[s].resize(fvSolver().a_noCells());
  }
 
  for(MInt s = 0; s < lsSolver().m_noSets; s++) {
    for(MInt fc = 0; fc < a_noFvCells(); fc++) {
      MInt lsId = fv2lsIdParent(fc);
      fvSolver().a_levelSetFunction(fc, s) = a_outsideGValue();
      if(lsId > -1) {
        fvSolver().a_curvatureG(fc, s) = lsSolver().a_curvatureG(lsId, 0);
        if(fvSolver().a_isInterface(fc)) {
          // Interplate levelset for cut-cells
          MFloat point[nDim];
          for(MInt d = 0; d < nDim; d++) {
            point[d] = fvSolver().a_coordinate(fc, d);
          }
          fvSolver().a_levelSetFunction(fc, s) = interpolateLevelSet(fv2lsIdParent(fc), point, s);
        } else {
          fvSolver().a_levelSetFunction(fc, s) = lsSolver().a_levelSetFunctionG(fv2lsIdParent(fc), s);
        }
      }
    }
  }
}

Friends And Related Function Documentation

◆ FvCartesianSolverXD< nDim, SysEqn >

template<MInt nDim_, class SysEqn >

friend class FvCartesianSolverXD< nDim, SysEqn >

friend

Definition at line 29 of file lsfv.h.

◆ LsCartesianSolver< nDim >

template<MInt nDim_, class SysEqn >

friend class LsCartesianSolver< nDim >

friend

Definition at line 29 of file lsfv.h.

Member Data Documentation

◆ m_bandWidthRef

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::m_bandWidthRef

private

Definition at line 63 of file lsfv.h.

◆ m_bandWidthRefMax

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::m_bandWidthRefMax

private

Definition at line 64 of file lsfv.h.

◆ m_cfl

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_cfl

private

Definition at line 75 of file lsfv.h.

◆ m_fvSolverId

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::m_fvSolverId

private

Definition at line 66 of file lsfv.h.

◆ m_G0regionId

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::m_G0regionId

private

Definition at line 69 of file lsfv.h.

◆ m_hadGapCells

template<MInt nDim_, class SysEqn >

MInt* CouplingLsFv< nDim_, SysEqn >::m_hadGapCells {}

private

Definition at line 72 of file lsfv.h.

◆ m_hasGapCells

template<MInt nDim_, class SysEqn >

MInt* CouplingLsFv< nDim_, SysEqn >::m_hasGapCells {}

private

Definition at line 73 of file lsfv.h.

◆ m_initialCrankAngle

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_initialCrankAngle

private

Definition at line 70 of file lsfv.h.

◆ m_lsSolverId

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::m_lsSolverId

private

Definition at line 67 of file lsfv.h.

◆ m_maxNoEmbeddedBodies

template<MInt nDim_, class SysEqn >

constexpr const MInt CouplingLsFv< nDim_, SysEqn >::m_maxNoEmbeddedBodies = 20

staticconstexprprivate

Definition at line 80 of file lsfv.h.

◆ m_maxVelocity

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_maxVelocity

private

Definition at line 76 of file lsfv.h.

◆ m_solverMethod

template<MInt nDim_, class SysEqn >

MString CouplingLsFv< nDim_, SysEqn >::m_solverMethod

private

Definition at line 78 of file lsfv.h.

◆ m_static_computeBodyProperties_amplitude

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_amplitude[m_maxNoEmbeddedBodies] {}

private

Definition at line 84 of file lsfv.h.

◆ m_static_computeBodyProperties_bodyToFunction

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_bodyToFunction[m_maxNoEmbeddedBodies] {}

private

Definition at line 96 of file lsfv.h.

◆ m_static_computeBodyProperties_circleStartAngle

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_circleStartAngle[m_maxNoEmbeddedBodies] {}

private

Definition at line 94 of file lsfv.h.

◆ m_static_computeBodyProperties_first

template<MInt nDim_, class SysEqn >

MBool CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_first = true

private

Definition at line 83 of file lsfv.h.

◆ m_static_computeBodyProperties_freqFactor

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_freqFactor[m_maxNoEmbeddedBodies] {}

private

Definition at line 85 of file lsfv.h.

◆ m_static_computeBodyProperties_initialBodyCenter

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_initialBodyCenter[m_maxNoEmbeddedBodies *3] {}

private

Definition at line 86 of file lsfv.h.

◆ m_static_computeBodyProperties_liftEndAngle1

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_liftEndAngle1[m_maxNoEmbeddedBodies] {}

private

Definition at line 91 of file lsfv.h.

◆ m_static_computeBodyProperties_liftEndAngle2

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_liftEndAngle2[m_maxNoEmbeddedBodies] {}

private

Definition at line 93 of file lsfv.h.

◆ m_static_computeBodyProperties_liftStartAngle1

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_liftStartAngle1[m_maxNoEmbeddedBodies] {}

private

Definition at line 90 of file lsfv.h.

◆ m_static_computeBodyProperties_liftStartAngle2

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_liftStartAngle2[m_maxNoEmbeddedBodies] {}

private

Definition at line 92 of file lsfv.h.

◆ m_static_computeBodyProperties_mu

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_mu[m_maxNoEmbeddedBodies] {}

private

Definition at line 88 of file lsfv.h.

◆ m_static_computeBodyProperties_mu2

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_mu2[m_maxNoEmbeddedBodies] {}

private

Definition at line 89 of file lsfv.h.

◆ m_static_computeBodyProperties_normal

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_normal[m_maxNoEmbeddedBodies *3] {}

private

Definition at line 95 of file lsfv.h.

◆ m_static_computeBodyProperties_omega

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_omega {}

private

Definition at line 97 of file lsfv.h.

◆ m_static_computeBodyProperties_rotAngle

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_rotAngle {}

private

Definition at line 98 of file lsfv.h.

◆ m_static_computeBodyProperties_Strouhal

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_computeBodyProperties_Strouhal {}

private

Definition at line 87 of file lsfv.h.

◆ m_static_crankAngle_first

template<MInt nDim_, class SysEqn >

MBool CouplingLsFv< nDim_, SysEqn >::m_static_crankAngle_first = true

private

Definition at line 101 of file lsfv.h.

◆ m_static_crankAngle_Strouhal

template<MInt nDim_, class SysEqn >

MFloat CouplingLsFv< nDim_, SysEqn >::m_static_crankAngle_Strouhal

private

Definition at line 102 of file lsfv.h.

◆ m_timeStepMethod

template<MInt nDim_, class SysEqn >

MInt CouplingLsFv< nDim_, SysEqn >::m_timeStepMethod

private

Definition at line 77 of file lsfv.h.

◆ nDim

template<MInt nDim_, class SysEqn >

constexpr MInt CouplingLsFv< nDim_, SysEqn >::nDim = nDim_

staticconstexpr

Definition at line 29 of file lsfv.h.

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

/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/COUPLER/lsfv.h
/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/COUPLER/lsfv.cpp

Public Types

Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Friends

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ Cell

◆ FvCartesianSolver

◆ LsSolver

Constructor & Destructor Documentation

◆ CouplingLsFv()

Member Function Documentation

◆ a_G0CellId()

◆ a_Ma()

◆ a_meanCoord()

◆ a_noFvCells()

◆ a_noFvGridCells()

◆ a_noRKSteps()

◆ a_physicalTime()

◆ a_RKStep()

◆ a_time()

◆ a_timeRef()

◆ a_TInfinity()

◆ a_UInfinity()

◆ checkProperties()

◆ cleanUp()

◆ computeBodyProperties()

◆ computeGCellTimeStep()

◆ crankAngle()

◆ finalizeCouplerInit()

◆ finalizeSubCoupleInit()

◆ fv2lsId()

◆ fv2lsIdParent()

◆ init()

◆ initData()

◆ interpolateLevelSet()

◆ ls2fvId()

◆ ls2fvIdParent()

◆ lsTimeStep()

◆ noLevelSetFieldData()

◆ postAdaptation()

◆ postCouple()

◆ preCouple()

◆ readProperties()

◆ returnStep()

◆ returnStep_semiLagrange()

◆ setLsInList()

◆ subCouple()

◆ testCoupling()

◆ testGapProperty()

◆ transferGapCellProperty()

◆ transferLevelSetValues()

Friends And Related Function Documentation

◆ FvCartesianSolverXD< nDim, SysEqn >

◆ LsCartesianSolver< nDim >

Member Data Documentation

◆ m_bandWidthRef

◆ m_bandWidthRefMax

◆ m_cfl

◆ m_fvSolverId

◆ m_G0regionId

◆ m_hadGapCells

◆ m_hasGapCells

◆ m_initialCrankAngle

◆ m_lsSolverId

◆ m_maxNoEmbeddedBodies

◆ m_maxVelocity

◆ m_solverMethod

◆ m_static_computeBodyProperties_amplitude

◆ m_static_computeBodyProperties_bodyToFunction

◆ m_static_computeBodyProperties_circleStartAngle

◆ m_static_computeBodyProperties_first

◆ m_static_computeBodyProperties_freqFactor

◆ m_static_computeBodyProperties_initialBodyCenter

◆ m_static_computeBodyProperties_liftEndAngle1

◆ m_static_computeBodyProperties_liftEndAngle2