#include <lslb.h>

Inheritance diagram for LsLb< nDim_, nDist, SysEqn >:

Collaboration diagram for LsLb< nDim_, nDist, SysEqn >:

Public Types
using	LsSolver = LsCartesianSolver< nDim >

using	LbSolver = LbSolverDxQy< nDim, nDist, SysEqn >

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

Public Types inherited from CouplingLB< nDim_, nDist, SysEqn >
using	solverType = LbSolverDxQy< nDim, nDist, SysEqn >

using	LbBndCnd = LbBndCndDxQy< nDim, nDist, SysEqn >

using	MbCellCollector = maia::lb::collector::LbMbCellCollector< nDim >

Public Member Functions
	LsLb (MInt couplingId, LsSolver ls, LbSolver lb)

void	init () override

void	preCouple (MInt step) override

void	postCouple (MInt step) override

void	finalizeCouplerInit () override

void	finalizeSubCoupleInit (MInt couplingStep) override

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

void	cleanUp () override

void	postAdaptation () override
	finalizeAdaptation More...

void	finalizeAdaptation (const MInt solverId) override

MInt	ls2lbId (const MInt lsId)

MInt	ls2lbIdParent (const MInt lsId)

MInt	lb2lsId (const MInt lbId)

MInt	lb2lsIdParent (const MInt lbId)

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

void	updateBoundaryCellsFromGField ()

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

void	updateLevelSetFlowSolver ()
	Updates the fv-mb-solver flow solver (after a completed levelSet TimeStep and finalizeLevelSet() ) More...

void	updateFlowSolverLevelSet ()

void	transferLevelSetFieldValues (MBool)
	transfers the LevelSetValues for all cells from the levelset to the moving boundary Part More...

void	buildCollectedLevelSet (const MInt cellId)
	build the combined levelSet for the given cellId More...

MFloat	interpolateLevelSet (MInt interpolationCells, MFloat point, const MInt set)

void	interpolateLsLb (const MInt from, const MInt to)
	interpolate levelset values on the lb-grid More...

MInt	noLevelSetFieldData ()

void	initializeSolidDomain ()

void	constructGField ()

template<MInt bodyType>
void	constructGField_ ()

template<MInt bodyType>
void	descendLevelSetValue (const MInt cellId, const MInt *bodyId, const MInt bodyCnt)

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...

Public Member Functions inherited from CouplingLB< nDim_, nDist, SysEqn >
	CouplingLB (const MInt couplingId, Solver *solvers, const MInt noSolvers=1)

	CouplingLB (const MInt couplingId, std::vector< solverType * > solvers)

MFloat	a_physicalTime () const

MFloat	lsTimeStep () const

MInt	a_RKStep () const

MInt	a_noLbCells (const MInt id=0) const

MInt	a_noLevelSetsMb (const MInt id=0) const

MFloat	a_Ma (const MInt id=0) const

MFloat	a_Re (const MInt id=0) const

MInt	a_pvu (const MInt id=0) const

MInt	a_pvv (const MInt id=0) const

MInt	a_pvw (const MInt id=0) const

MInt	a_pvrho (const MInt id=0) const

MInt	a_pvt (const MInt id=0) const

MInt	a_isThermal (const MInt id=0) const

MInt	a_noDistributions (const MInt id=0) const

MFloat	a_initTemperatureKelvin (const MInt id=0) const

MFloat	a_time () const

MbCellCollector &	a_mbCell (const MInt id=0)

MInt	a_boundaryCellMb (const MInt cellId, const MInt id=0)

MFloat &	a_levelSetFunctionMb (const MInt cellId, const MInt set, const MInt id=0)

MFloat	a_levelSetFunctionMb (const MInt cellId, const MInt set, const MInt id=0) const

MInt &	a_associatedBodyIdsMb (const MInt cellId, const MInt set, const MInt id=0)

MInt	a_associatedBodyIdsMb (const MInt cellId, const MInt set, const MInt id=0) const

MInt	a_parentId (const MInt cellId, const MInt id=0)

MInt	a_childId (const MInt cellId, const MInt child, const MInt id=0)

MInt	minCell (const MInt index, const MInt id=0) const

MInt	noMinCells (const MInt id=0) const

MInt	a_noCells (const MInt id=0) const

MFloat	a_cellLengthAtLevel (MInt level, const MInt id=0)

MInt	a_noEmbeddedBodiesLB (const MInt id=0) const

MBool	a_isActive (const MInt cellId, const MInt id=0) const

MBool	a_wasActive (const MInt cellId, const MInt id=0) const

MInt	a_noVariables (const MInt id=0) const

MFloat &	a_variable (const MInt cellId, const MInt varId, const MInt id=0)

MFloat &	a_oldVariable (const MInt cellId, const MInt varId, const MInt id=0)

MInt	a_bndCellId (const MInt bndCell, const MInt id=0)

MInt	a_noBndCells (const MInt id=0)

Public Attributes
MFloat *	m_transferBoundingBox {}

Static Public Attributes
static constexpr MInt	nDim = nDim_

static constexpr const MBool	m_constructGField = false

Private Member Functions
void	initData ()

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

void	readProperties ()

void	updateGeometry ()
	Updates the member-variables in the geometry-intersection class. More...

Private Attributes
MInt	m_lbSolverId {}

MInt	m_lsSolverId {}

MBool	m_outsideDefault {}

maia::coupling::Mapping	bndryToBodyMapping

maia::coupling::Mapping	bodyToBndryMapping

MBool	m_static_computeBodyProperties_first = true

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_amplitude {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_freqFactor {}

std::array< MFloat, m_maxNoEmbeddedBodies *3 >	m_static_computeBodyProperties_initialBodyCenter {}

MFloat	m_static_computeBodyProperties_Strouhal {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_mu {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_mu2 {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_liftStartAngle1 {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_liftEndAngle1 {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_liftStartAngle2 {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_liftEndAngle2 {}

std::array< MFloat, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_circleStartAngle {}

std::array< MFloat, m_maxNoEmbeddedBodies *3 >	m_static_computeBodyProperties_normal {}

std::array< MInt, m_maxNoEmbeddedBodies >	m_static_computeBodyProperties_bodyToFunction {}

MFloat	m_static_computeBodyProperties_omega {}

MFloat	m_static_computeBodyProperties_rotAngle {}

MBool	m_static_updateLevelSetFlowSolver_firstRun = true

Static Private Attributes
static constexpr const MInt	m_maxNoEmbeddedBodies = 1

static constexpr const MInt	m_noCorners = (nDim == 2) ? 4 : 8

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

MFloat	returnIdleRecord () const

Protected Member Functions inherited from CouplingLB< nDim_, nDist, SysEqn >
MInt	noSolvers () const

solverType &	lbSolver (const MInt solverId=0) const

LbBndCnd &	lbBndCnd (const MInt id=0)

Detailed Description

template<MInt nDim_, MInt nDist, class SysEqn>
class LsLb< nDim_, nDist, SysEqn >

Definition at line 25 of file lslb.h.

Member Typedef Documentation

◆ LbSolver

template<MInt nDim_, MInt nDist, class SysEqn >

using LsLb< nDim_, nDist, SysEqn >::LbSolver = LbSolverDxQy<nDim, nDist, SysEqn>

Definition at line 32 of file lslb.h.

◆ LsSolver

template<MInt nDim_, MInt nDist, class SysEqn >

using LsLb< nDim_, nDist, SysEqn >::LsSolver = LsCartesianSolver<nDim>

Definition at line 31 of file lslb.h.

Constructor & Destructor Documentation

◆ LsLb()

template<MInt nDim, MInt nDist, class SysEqn >

LsLb< nDim, nDist, SysEqn >::LsLb	(	MInt	couplingId,
		LsSolver *	ls,
		LbSolver *	lb
	)

Definition at line 23 of file lslb.cpp.

  : Coupling(couplingId), CouplingLS<nDim>(couplingId, ls), CouplingLB<nDim, nDist, SysEqn>(couplingId, lb) {
  TRACE();
 
  initData();
  readProperties();
  checkProperties();
  updateGeometry();
}

Member Function Documentation

◆ buildCollectedLevelSet()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::buildCollectedLevelSet ( const MInt cellId )

Author: Tim Wegmann

Definition at line 598 of file lslb.cpp.

                                                                        {
  TRACE();
 
  lbSolver().a_levelSetFunctionMB(cellId, 0) = lbSolver().a_levelSetFunctionMB(cellId, 1);
  lbSolver().a_associatedBodyIds(cellId, 0) = lbSolver().a_associatedBodyIds(cellId, 1);
  for(MInt set = 2; set < a_noLevelSetsMb(); set++) {
    MFloat phi0 = lbSolver().a_levelSetFunctionMB(cellId, 0);
    MFloat phi1 = lbSolver().a_levelSetFunctionMB(cellId, set);
    MInt body0 = lbSolver().a_associatedBodyIds(cellId, 0);
    MInt body1 = lbSolver().a_associatedBodyIds(cellId, set);
 
    //
    if(phi0 >= F0 && phi1 >= F0) {
      if(phi1 < phi0) {
        body0 = body1;
      }
      phi0 = mMin(phi0, phi1);
    } else if(phi0 <= F0 && phi1 <= F0) {
      if(abs(phi1) > abs(phi0)) {
        body0 = body1;
      }
      phi0 = -sqrt(phi0 * phi0 + phi1 * phi1);
    } else if(phi0 * phi1 <= F0) {
      if(phi0 < F0) {
        // phi0 = phi0;
      } else {
        phi0 = phi1;
        body0 = body1;
      }
    }
 
    lbSolver().a_levelSetFunctionMB(cellId, 0) = phi0;
    lbSolver().a_associatedBodyIds(cellId, 0) = body0;
  }
}

◆ checkProperties()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::checkProperties

privatevirtual

Implements Coupling.

Definition at line 70 of file lslb.cpp.

                                                {
  TRACE();
 
  lbSolver().m_noEmbeddedBodies = lsSolver().m_noEmbeddedBodies;
 
  if(lsSolver().m_maxNoSets > 0) lbSolver().m_noLevelSetsUsedForMb = lsSolver().m_maxNoSets;
  if(lbSolver().m_useOnlyCollectedLS) lbSolver().m_noLevelSetsUsedForMb = 1;
 
  lbSolver().m_maxNoSets = lbSolver().m_noLevelSetsUsedForMb;
 
  lbSolver().m_constructGField = lsSolver().m_constructGField;
}

◆ cleanUp()

template<MInt nDim_, MInt nDist, class SysEqn >

void LsLb< nDim_, nDist, SysEqn >::cleanUp ( )

inlineoverridevirtual

Implements Coupling.

Definition at line 128 of file lslb.h.

128{};

◆ constructGField()

template<MInt nDim_, MInt nDist, class SysEqn >

void LsLb< nDim_, nDist, SysEqn >::constructGField ( )

◆ constructGField_()

template<MInt nDim_, MInt nDist, class SysEqn >

template<MInt bodyType>

void LsLb< nDim_, nDist, SysEqn >::constructGField_ ( )

◆ createBodyTree()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::createBodyTree

Definition at line 344 of file lslb.cpp.

                                               {
  TRACE();
  IF_CONSTEXPR(nDim == 3) { updateGeometry(); }
}

◆ descendLevelSetValue()

template<MInt nDim_, MInt nDist, class SysEqn >

template<MInt bodyType>

void LsLb< nDim_, nDist, SysEqn >::descendLevelSetValue	(	const MInt	cellId,
		const MInt *	bodyId,
		const MInt	bodyCnt
	)

◆ finalizeAdaptation()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::finalizeAdaptation ( const MInt solverId )

overridevirtual

Reimplemented from Coupling.

Definition at line 505 of file lslb.cpp.

                                                                      {
  TRACE();
 
  if(solverId == 1) {
    if(lbSolver().isActive()) {
      lbSolver().initializeMovingBoundaries();
      lbBndCnd().initializeBndMovingBoundaries();
    }
    updateLevelSetFlowSolver();
  }
}

◆ finalizeCouplerInit()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::finalizeCouplerInit

overridevirtual

Implements Coupling.

Definition at line 422 of file lslb.cpp.

                                                    {
  TRACE();
  lsSolver().m_timeStep = F1;
}

◆ finalizeSubCoupleInit()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::finalizeSubCoupleInit ( MInt couplingStep )

overridevirtual

Implements Coupling.

Definition at line 428 of file lslb.cpp.

                                                                       {
  TRACE();
  std::ignore = couplingStep;
}

◆ init()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::init

overridevirtual

Implements Coupling.

Definition at line 411 of file lslb.cpp.

                                     {
  TRACE();
 
  if(lbSolver().isActive()) {
    lbSolver().initializeMovingBoundaries();
    lbBndCnd().initializeBndMovingBoundaries();
  }
  updateLevelSetFlowSolver();
}

◆ initData()

template<MInt nDim_, MInt nDist, class SysEqn >

void LsLb< nDim_, nDist, SysEqn >::initData ( )

private

◆ initializeSolidDomain()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::initializeSolidDomain

Definition at line 519 of file lslb.cpp.

                                                      {
  TRACE();
 
  if(!lbSolver().isActive()) {
    return;
  }
 
  MFloatScratchSpace bodyVelocities(a_noEmbeddedBodies(), nDim, AT_, "bodyVelocities");
 
  // Get the body Velocity for each embedded body
  for(MInt body = 0; body < a_noEmbeddedBodies(); body++) {
    lsSolver().computeBodyPropertiesForced(2, &bodyVelocities(body, 0), body, globalTimeStep);
  }
 
  for(MInt i = 0; i < lbSolver().a_noCells(); i++) {
    // Regular fluid cell
    if(a_isActive(i) && a_wasActive(i)) {
      continue;
    }
 
    // Regular Solid cell
    if(!a_isActive(i)) {
      // TODO labels:COUPLER,toremove Skip non-leaf cells ...
      /*if(!lbSolver().c_isLeafCell(i)){
        continue;
      }*/
 
      // determine the LS-Body to which the cell belongs, to set the right body velocity
      MInt bodyId = -1;
      MInt setOfBody = 0;
      for(MInt set = lbSolver().m_levelSetId; set < lbSolver().m_maxNoSets; set++) {
        if(a_associatedBodyIdsMb(i, set) >= 0) {
          bodyId = a_associatedBodyIdsMb(i, set);
          setOfBody = set;
          break;
        }
      }
 
      // the Velocity of the deactivated cell is set to the body velocity
      // the Density is set to 1.0
      if((bodyId >= 0) && (bodyId < a_noEmbeddedBodies()) && (a_levelSetFunctionMb(i, setOfBody) < 0)) {
        for(MInt j = 0; j < nDim; j++) {
          a_variable(i, j) = bodyVelocities(bodyId, j);
          a_oldVariable(i, j) = bodyVelocities(bodyId, j);
        }
        MFloat squaredVelocity = 0.0;
        for(MInt n = 0; n < nDim; n++) {
          squaredVelocity += bodyVelocities(bodyId, n) * bodyVelocities(bodyId, n);
        }
        lbSolver().setEqDists(i, 1.0, &bodyVelocities(bodyId, 0));
        if(a_isThermal()) {
          a_variable(i, a_pvt()) = a_initTemperatureKelvin();
          lbSolver().setEqDistsThermal(i, a_initTemperatureKelvin(), F1, squaredVelocity, &bodyVelocities(bodyId, 0));
        }
      } else {
        for(MInt j = 0; j < a_noVariables(); j++) {
          a_variable(i, j) = F0;
          a_oldVariable(i, j) = F0;
        }
        MFloat squaredVelocity = 0.0;
        MFloat bodyVelocitiesOutsideLs[nDim] = {F0};
        lbSolver().setEqDists(i, 1.0, bodyVelocitiesOutsideLs);
        if(a_isThermal()) {
          a_variable(i, a_pvt()) = a_initTemperatureKelvin();
          lbSolver().setEqDistsThermal(i, a_initTemperatureKelvin(), F1, squaredVelocity, bodyVelocitiesOutsideLs);
        }
      }
      a_variable(i, a_pvrho()) = 1.0;
      a_oldVariable(i, a_pvrho()) = 1.0;
    } else {
      lbBndCnd().refillEmergedCell(i);
    }
  }
}

◆ interpolateLevelSet()

template<MInt nDim, MInt nDist, class SysEqn >

MFloat LsLb< nDim, nDist, SysEqn >::interpolateLevelSet	(	MInt *	interpolationCells,
		MFloat *	point,
		const MInt	set
	)

interpolates the levelSet value further on fv-cells

Author: Tim Wegmann

Date: 2020-04-01

Definition at line 663 of file lslb.cpp.

                                                                                                             {
  TRACE();
 
  std::function<MFloat(const MInt, const MInt)> scalarField = [&](const MInt cellId, const MInt refSet) {
    return static_cast<MFloat>(lbSolver().a_levelSetFunctionMB(cellId, refSet));
  };
 
  std::function<MFloat(const MInt, const MInt)> coordinate = [&](const MInt cellId, const MInt id) {
    return static_cast<MFloat>(lbSolver().a_coordinate(cellId, id));
  };
 
  return lbSolver().interpolateFieldDataLb(&interpolationCells[0], &point[0], set, scalarField, coordinate);
}

◆ interpolateLsLb()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::interpolateLsLb	(	const MInt	from,
		const MInt	to
	)

Author: Tim Wegmann

Date: March 2020

Definition at line 640 of file lslb.cpp.

                                                                              {
  TRACE();
 
  MInt interpolationCells[8] = {0, 0, 0, 0, 0, 0, 0, 0};
  MFloat point[3] = {lbSolver().a_coordinate(to, 0), lbSolver().a_coordinate(to, 1), lbSolver().a_coordinate(to, 2)};
 
  const MInt position = lbSolver().setUpLbInterpolationStencil(from, interpolationCells, point);
 
  for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
    if(position < 0) {
      lbSolver().a_levelSetFunctionMB(to, set) = lbSolver().a_levelSetFunctionMB(from, set);
    } else {
      MFloat phi = interpolateLevelSet(interpolationCells, point, set);
      lbSolver().a_levelSetFunctionMB(to, set) = phi;
    }
  }
  buildCollectedLevelSet(to);
}

◆ lb2lsId()

template<MInt nDim_, MInt nDist, class SysEqn >

MInt LsLb< nDim_, nDist, SysEqn >::lb2lsId ( const MInt lbId )

inline

Definition at line 137 of file lslb.h.

137{ return convertId(lbSolver(), lsSolver(), lbId); };

convertId

MInt convertId(SolverA &solverA, SolverB &solverB, const MInt solverAId)

Conversion from solverA id to the solverB id on the same-level only!

Definition: couplingutils.h:21

◆ lb2lsIdParent()

template<MInt nDim_, MInt nDist, class SysEqn >

MInt LsLb< nDim_, nDist, SysEqn >::lb2lsIdParent ( const MInt lbId )

inline

Definition at line 138 of file lslb.h.

138{ return convertIdParent(lbSolver(), lsSolver(), lbId); };

convertIdParent

MInt convertIdParent(SolverA &solverA, SolverB &solverB, const MInt solverAId)

Conversion from solverA id to the solverB id If no cell on the same level is found,...

Definition: couplingutils.h:46

◆ ls2lbId()

template<MInt nDim_, MInt nDist, class SysEqn >

MInt LsLb< nDim_, nDist, SysEqn >::ls2lbId ( const MInt lsId )

inline

Definition at line 135 of file lslb.h.

135{ return convertId(lsSolver(), lbSolver(), lsId); };

◆ ls2lbIdParent()

template<MInt nDim_, MInt nDist, class SysEqn >

MInt LsLb< nDim_, nDist, SysEqn >::ls2lbIdParent ( const MInt lsId )

inline

Definition at line 136 of file lslb.h.

136{ return convertIdParent(lsSolver(), lbSolver(), lsId); };

◆ noLevelSetFieldData()

template<MInt nDim, MInt nDist, class SysEqn >

MInt LsLb< nDim, nDist, SysEqn >::noLevelSetFieldData

Definition at line 352 of file lslb.cpp.

                                                    {
  MInt noLevelSetFieldData = 0;
 
  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, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::postAdaptation

overridevirtual

Author: Moritz Waldmann

Reimplemented from Coupling.

Definition at line 500 of file lslb.cpp.

                                               {
  TRACE();
}

◆ postCouple()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::postCouple ( MInt step )

overridevirtual

Implements Coupling.

Definition at line 478 of file lslb.cpp.

                                                    {
  TRACE();
 
  switch(step) {
    case 0: {
      break;
    }
    case 1: {
      lbBndCnd().postCouple();
      lsSolver().m_timeStep = F1;
      break;
    }
    default: {
      break;
    }
  }
}

◆ preCouple()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::preCouple ( MInt step )

overridevirtual

Implements Coupling.

Definition at line 434 of file lslb.cpp.

                                                   {
  TRACE();
 
  switch(step) {
    case 0: {
      break;
    }
    case 1: {
      updateLevelSetFlowSolver();
      std::vector<MInt> maxGCellLevels(lsSolver().m_maxNoSets);
      for(MInt set = 0; set < lsSolver().m_maxNoSets; set++) {
        maxGCellLevels[set] = lsSolver().a_maxGCellLevel(set);
      }
 
      lbSolver().preCoupleLs(maxGCellLevels);
 
      lbSolver().createBndryToBodyMapping(bndryToBodyMapping, bodyToBndryMapping);
 
      lbBndCnd().createMBComm();
 
      initializeSolidDomain();
 
      MFloatScratchSpace bodyVelocities(m_maxNoEmbeddedBodies, nDim, AT_, "bodyVelocities");
      for(MInt body = 0; body < m_maxNoEmbeddedBodies; body++) {
        lsSolver().computeBodyPropertiesForced(2, &bodyVelocities(body, 0), body, globalTimeStep);
      }
 
      for(MInt mbCell = 0; mbCell < a_mbCell().size(); mbCell++) {
        for(auto body : bndryToBodyMapping[mbCell]) {
          for(MInt n = 0; n < nDim; n++) {
            a_mbCell().velocity(mbCell, n) = bodyVelocities(body, n);
          }
        }
      }
 
      break;
    }
    default: {
      break;
    }
  }
}

◆ readProperties()

template<MInt nDim_, MInt nDist, class SysEqn >

void LsLb< nDim_, nDist, SysEqn >::readProperties ( )

privatevirtual

Implements Coupling.

◆ subCouple()

template<MInt nDim_, MInt nDist, class SysEqn >

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

inlineoverridevirtual

Implements Coupling.

Definition at line 127 of file lslb.h.

127{};

◆ testCoupling()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::testCoupling

Definition at line 114 of file lslb.cpp.

                                             {
  TRACE();
 
  if(lsSolver().m_constructGField) return;
 
  for(MInt lbCellId = 0; lbCellId < a_noLbCells(); lbCellId++) {
    ASSERT(lbSolver().grid().tree().solver2grid(lbCellId) >= 0, "");
    ASSERT(lbSolver().grid().solverFlag(lbSolver().grid().tree().solver2grid(lbCellId), m_lbSolverId), "");
#ifdef COUPLING_DEBUG_
    if(lsSolver().grid().solverFlag(lbSolver().grid().tree().solver2grid(lbCellId), m_lsSolverId)) {
      ASSERT(ls2lbId(lb2lsId(lbCellId)) == lbCellId,
             to_string(lbCellId) + " " + to_string(lb2lsId(lbCellId)) + " " + to_string(ls2lbId(lb2lsId(lbCellId))));
    }
#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(lbSolver().grid().solverFlag(lsSolver().grid().tree().solver2grid(cellId), m_lbSolverId)) {
      ASSERT(lb2lsId(ls2lbId(cellId)) == cellId,
             to_string(cellId) + " " + to_string(ls2lbId(cellId)) + " " + to_string(lb2lsId(ls2lbId(cellId))));
    }
#endif
  }
 
  if(!lsSolver().m_levelSetLb) return;
 
#ifdef COUPLING_DEBUG_
  for(MInt lbCellId = 0; lbCellId < a_noLbCells(); lbCellId++) {
    for(MInt dir = 0; dir < nDim; dir++) {
      if(lsSolver().grid().solverFlag(lbSolver().grid().tree().solver2grid(lbCellId), m_lsSolverId)) {
        ASSERT(abs(lbSolver().c_coordinate(lbCellId, dir) - a_coordinateG(lb2lsId(lbCellId), dir)) < 0.00000001, "");
      }
    }
  }
#endif
}

◆ transferLevelSetFieldValues()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::transferLevelSetFieldValues ( MBool exchangeLVS )

Author: Christoph Siewert, Tim Wegmann

Definition at line 159 of file lslb.cpp.

                                                                             {
  TRACE();
 
  const MInt noRfJumps = mMax(0, lbSolver().maxLevel() - lsSolver().maxRefinementLevel());
 
  // Nothing to transfer if lbSolver is inactive!
  if(!lbSolver().isActive()) return;
 
  if(exchangeLVS) {
    if(lsSolver().isActive()) {
      lsSolver().exchangeAllLevelSetData();
      lsSolver().checkHaloCells();
    }
  }
 
  if(lsSolver().isActive()) {
    std::list<MInt> interpolationParents;
    interpolationParents.clear();
 
    for(MInt cellId = 0; cellId < a_noLbCells(); cellId++) {
      // reset values with invalid/defaults
      for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
        if(m_outsideDefault) {
          lbSolver().a_levelSetFunctionMB(cellId, set) = a_outsideGValue();
        } else {
          lbSolver().a_levelSetFunctionMB(cellId, set) = -a_outsideGValue();
        }
        lbSolver().a_associatedBodyIds(cellId, set) = -2;
      }
      if(m_transferBoundingBox != nullptr) {
        if(m_transferBoundingBox[0] > lbSolver().c_coordinate(cellId, 0)
           || m_transferBoundingBox[0 + nDim] < lbSolver().c_coordinate(cellId, 0)) {
          continue;
        }
        if(m_transferBoundingBox[1] > lbSolver().c_coordinate(cellId, 1)
           || m_transferBoundingBox[1 + nDim] < lbSolver().c_coordinate(cellId, 1)) {
          continue;
        }
        IF_CONSTEXPR(nDim == 3) {
          if(m_transferBoundingBox[2] > lbSolver().c_coordinate(cellId, 2)
             || m_transferBoundingBox[2 + nDim] < lbSolver().c_coordinate(cellId, 2)) {
            continue;
          }
        }
      }
      const MInt gCellId = lb2lsId(cellId);
      // direct transfer of matching cells
      if(gCellId > -1) {
        for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
          lbSolver().a_levelSetFunctionMB(cellId, set) = a_levelSetFunctionG(gCellId, set);
          lbSolver().a_associatedBodyIds(cellId, set) = a_bodyIdG(gCellId, set);
        }
      } else {
        if(!g_multiSolverGrid) {
          ASSERT(cellId >= a_noLbCells(), "");
        } else {
          // multiSolverGrid, where the level and the grid-extension might differ!
          if(cellId <= a_noLbCells()) { // grid-cell!
            const MInt gCellParent = lb2lsIdParent(cellId);
 
            if(gCellParent > -1) { // just a different level
 
              if(noRfJumps == 0) {
                // same maxRefinementLevel, meaning just a difference due to different refinement-widths
                for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
                  lbSolver().a_associatedBodyIds(cellId, set) = a_bodyIdG(gCellParent, set);
                  // ASSERT(!lsSolver().a_inBandG( gCellParent ,  set ), "");
                  lbSolver().a_levelSetFunctionMB(cellId, set) = a_levelSetFunctionG(gCellParent, set);
                }
              } else {
                // ls-interpolation towards higher lb-level
                if(lbSolver().a_isHalo(cellId) || lsSolver().a_isHalo(gCellParent)) continue;
                // simple transfer if not a band-cell
                if(lsSolver().m_buildCollectedLevelSetFunction && !lsSolver().a_inBandG(gCellParent, 0)) {
                  for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
                    lbSolver().a_levelSetFunctionMB(cellId, set) = a_levelSetFunctionG(gCellParent, set);
                    lbSolver().a_associatedBodyIds(cellId, set) = a_bodyIdG(gCellParent, set);
                  }
                  continue;
                }
                const MInt levelDifference = lbSolver().c_level(cellId) - lsSolver().a_level(gCellParent);
                if(levelDifference == 1) {
                  // direct interpolation from ls to fv for only 1 level difference
                  MFloat point[3] = {lbSolver().a_coordinate(cellId, 0), lbSolver().a_coordinate(cellId, 1),
                                     lbSolver().a_coordinate(cellId, 2)};
                  MInt interpolationCells[8] = {0, 0, 0, 0, 0, 0, 0, 0};
                  MInt position = lsSolver().setUpLevelSetInterpolationStencil(gCellParent, interpolationCells, point);
 
                  for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
                    lbSolver().a_associatedBodyIds(cellId, set) = a_bodyIdG(gCellParent, set);
                    if(position < 0) {
                      // no valid interpolation stencil found
                      lbSolver().a_levelSetFunctionMB(cellId, set) = a_levelSetFunctionG(gCellParent, set);
                    } else { // interpolation for all sets
                      const MFloat phi = lsSolver().interpolateLevelSet(interpolationCells, point, set);
 
                      lbSolver().a_levelSetFunctionMB(cellId, set) = phi;
                    }
                  }
                  // build collected levelset data
                  buildCollectedLevelSet(cellId);
 
                } else {
                  // add cells based on which the fv-interpolation needs to start
                  // that is cells for which the level difference is above 1
                  for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
                    lbSolver().a_associatedBodyIds(cellId, set) = a_bodyIdG(gCellParent, set);
                  }
                  MInt parent = cellId;
                  for(MInt i = 0; i < levelDifference - 1; i++) {
                    parent = lbSolver().c_parentId(parent);
                  }
                  ASSERT(parent > -1, "");
                  // parent of this parent has a direct link, meaning that the values for the parent
                  // are set above!
                  ASSERT(lb2lsId(lbSolver().c_parentId(parent)) == gCellParent, "");
                  interpolationParents.push_back(parent);
                }
              }
            }
          }
        }
      }
    }
 
    // now interpolate through all childs
    // iterate downwards level by level and exchange values in between!
    if(!interpolationParents.empty()) {
      cerr << "This part is not tested yet!!!" << endl;
      std::list<MInt> newInterpolationParents;
 
      for(MInt lvl = lsSolver().maxUniformRefinementLevel(); lvl < lbSolver().maxRefinementLevel(); lvl++) {
        lbSolver().exchangeData(&(lbSolver().a_levelSetFunctionMB(0, 0)), lbSolver().m_maxNoSets);
        lbSolver().exchangeData(&(lbSolver().a_associatedBodyIds(0, 0)), lbSolver().m_maxNoSets);
 
        MInt noInterpolationCells = (MInt)interpolationParents.size();
        MPI_Allreduce(MPI_IN_PLACE, &noInterpolationCells, 1, MPI_INT, MPI_MAX, lbSolver().mpiComm(), AT_, "INPLACE",
                      "noInterpolationCells");
        if(noInterpolationCells == 0) break;
 
        if(interpolationParents.empty()) continue;
 
        ASSERT(noRfJumps > 1, "");
 
        interpolationParents.sort();
        interpolationParents.unique();
        newInterpolationParents.clear();
 
        for(auto it = interpolationParents.begin(); it != interpolationParents.end(); it++) {
          const MInt parent = (*it);
          for(MInt c = 0; c < lbSolver().grid().m_maxNoChilds; c++) {
            const MInt childId = lbSolver().c_childId(parent, c);
            if(childId < 0) continue;
            // interpolateLsLb(parent, childId);
            if(!lbSolver().c_isLeafCell(childId)) {
              newInterpolationParents.push_back(childId);
            }
          }
        }
        interpolationParents.clear();
        for(auto it = newInterpolationParents.begin(); it != newInterpolationParents.end(); it++) {
          const MInt cellId = (*it);
          interpolationParents.push_back(cellId);
        }
      }
    }
  } else {
    for(MInt cellId = 0; cellId < a_noLbCells(); cellId++) {
      for(MInt set = 0; set < a_noLevelSetsMb(); set++) {
        if(m_outsideDefault) {
          lbSolver().a_levelSetFunctionMB(cellId, set) = a_outsideGValue();
        } else {
          lbSolver().a_levelSetFunctionMB(cellId, set) = -a_outsideGValue();
        }
        lbSolver().a_associatedBodyIds(cellId, set) = -2;
      }
    }
  }
}

◆ updateBoundaryCellsFromGField()

template<MInt nDim_, MInt nDist, class SysEqn >

void LsLb< nDim_, nDist, SysEqn >::updateBoundaryCellsFromGField ( )

◆ updateFlowSolverLevelSet()

template<MInt nDim_, MInt nDist, class SysEqn >

void LsLb< nDim_, nDist, SysEqn >::updateFlowSolverLevelSet ( )

◆ updateGeometry()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::updateGeometry

private

Author: Tim Wegmann

Definition at line 374 of file lslb.cpp.

                                               {
  TRACE();
 
  lbSolver().m_geometryIntersection->m_noEmbeddedBodies = lbSolver().m_noEmbeddedBodies;
  lbSolver().m_geometryIntersection->m_noLevelSetsUsedForMb = lbSolver().m_noLevelSetsUsedForMb;
  lbSolver().m_geometryIntersection->m_bodyToSetTable = lsSolver().m_bodyToSetTable;
  lbSolver().m_geometryIntersection->m_setToBodiesTable = lsSolver().m_setToBodiesTable;
  lbSolver().m_geometryIntersection->m_noBodiesInSet = lsSolver().m_noBodiesInSet;
}

◆ updateLevelSetFlowSolver()

template<MInt nDim, MInt nDist, class SysEqn >

void LsLb< nDim, nDist, SysEqn >::updateLevelSetFlowSolver

Author: Tim Wegmann

Definition at line 389 of file lslb.cpp.

                                                         {
  TRACE();
 
  if(lbSolver().m_constructGField) return;
 
  MBool& firstRun = m_static_updateLevelSetFlowSolver_firstRun;
 
  if(firstRun
     || (lbSolver().m_trackMovingBndry && globalTimeStep >= lbSolver().m_trackMbStart
         && globalTimeStep < lbSolver().m_trackMbEnd)) {
    testCoupling();
 
    // 1) transfer the levelSet Data to the flow solver
    transferLevelSetFieldValues(true);
  }
 
  firstRun = false;
}

Member Data Documentation

◆ bndryToBodyMapping

template<MInt nDim_, MInt nDist, class SysEqn >

maia::coupling::Mapping LsLb< nDim_, nDist, SysEqn >::bndryToBodyMapping

private

Definition at line 43 of file lslb.h.

◆ bodyToBndryMapping

template<MInt nDim_, MInt nDist, class SysEqn >

maia::coupling::Mapping LsLb< nDim_, nDist, SysEqn >::bodyToBndryMapping

private

Definition at line 44 of file lslb.h.

◆ m_constructGField

template<MInt nDim_, MInt nDist, class SysEqn >

constexpr const MBool LsLb< nDim_, nDist, SysEqn >::m_constructGField = false

staticconstexpr

Definition at line 110 of file lslb.h.

◆ m_lbSolverId

template<MInt nDim_, MInt nDist, class SysEqn >

MInt LsLb< nDim_, nDist, SysEqn >::m_lbSolverId {}

private

Definition at line 38 of file lslb.h.

◆ m_lsSolverId

template<MInt nDim_, MInt nDist, class SysEqn >

MInt LsLb< nDim_, nDist, SysEqn >::m_lsSolverId {}

private

Definition at line 39 of file lslb.h.

◆ m_maxNoEmbeddedBodies

template<MInt nDim_, MInt nDist, class SysEqn >

constexpr const MInt LsLb< nDim_, nDist, SysEqn >::m_maxNoEmbeddedBodies = 1

staticconstexprprivate

Definition at line 46 of file lslb.h.

◆ m_noCorners

template<MInt nDim_, MInt nDist, class SysEqn >

constexpr const MInt LsLb< nDim_, nDist, SysEqn >::m_noCorners = (nDim == 2) ? 4 : 8

staticconstexprprivate

Definition at line 47 of file lslb.h.

◆ m_outsideDefault

template<MInt nDim_, MInt nDist, class SysEqn >

MBool LsLb< nDim_, nDist, SysEqn >::m_outsideDefault {}

private

Definition at line 41 of file lslb.h.

◆ m_static_computeBodyProperties_amplitude

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_amplitude {}

private

Definition at line 51 of file lslb.h.

◆ m_static_computeBodyProperties_bodyToFunction

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MInt, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_bodyToFunction {}

private

Definition at line 63 of file lslb.h.

◆ m_static_computeBodyProperties_circleStartAngle

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_circleStartAngle {}

private

Definition at line 61 of file lslb.h.

◆ m_static_computeBodyProperties_first

template<MInt nDim_, MInt nDist, class SysEqn >

MBool LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_first = true

private

Definition at line 50 of file lslb.h.

◆ m_static_computeBodyProperties_freqFactor

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_freqFactor {}

private

Definition at line 52 of file lslb.h.

◆ m_static_computeBodyProperties_initialBodyCenter

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies * 3> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_initialBodyCenter {}

private

Definition at line 53 of file lslb.h.

◆ m_static_computeBodyProperties_liftEndAngle1

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_liftEndAngle1 {}

private

Definition at line 58 of file lslb.h.

◆ m_static_computeBodyProperties_liftEndAngle2

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_liftEndAngle2 {}

private

Definition at line 60 of file lslb.h.

◆ m_static_computeBodyProperties_liftStartAngle1

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_liftStartAngle1 {}

private

Definition at line 57 of file lslb.h.

◆ m_static_computeBodyProperties_liftStartAngle2

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_liftStartAngle2 {}

private

Definition at line 59 of file lslb.h.

◆ m_static_computeBodyProperties_mu

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_mu {}

private

Definition at line 55 of file lslb.h.

◆ m_static_computeBodyProperties_mu2

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_mu2 {}

private

Definition at line 56 of file lslb.h.

◆ m_static_computeBodyProperties_normal

template<MInt nDim_, MInt nDist, class SysEqn >

std::array<MFloat, m_maxNoEmbeddedBodies * 3> LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_normal {}

private

Definition at line 62 of file lslb.h.

◆ m_static_computeBodyProperties_omega

template<MInt nDim_, MInt nDist, class SysEqn >

MFloat LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_omega {}

private

Definition at line 64 of file lslb.h.

◆ m_static_computeBodyProperties_rotAngle

template<MInt nDim_, MInt nDist, class SysEqn >

MFloat LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_rotAngle {}

private

Definition at line 65 of file lslb.h.

◆ m_static_computeBodyProperties_Strouhal

template<MInt nDim_, MInt nDist, class SysEqn >

MFloat LsLb< nDim_, nDist, SysEqn >::m_static_computeBodyProperties_Strouhal {}

private

Definition at line 54 of file lslb.h.

◆ m_static_updateLevelSetFlowSolver_firstRun

template<MInt nDim_, MInt nDist, class SysEqn >

MBool LsLb< nDim_, nDist, SysEqn >::m_static_updateLevelSetFlowSolver_firstRun = true

private

Definition at line 68 of file lslb.h.

◆ m_transferBoundingBox

template<MInt nDim_, MInt nDist, class SysEqn >

MFloat* LsLb< nDim_, nDist, SysEqn >::m_transferBoundingBox {}

Definition at line 111 of file lslb.h.

◆ nDim

template<MInt nDim_, MInt nDist, class SysEqn >

constexpr MInt LsLb< nDim_, nDist, SysEqn >::nDim = nDim_

staticconstexpr

Definition at line 27 of file lslb.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/lslb.h
/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/COUPLER/lslb.cpp

Public Types

Public Member Functions

Public Attributes

Static Public Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ LbSolver

◆ LsSolver

Constructor & Destructor Documentation

◆ LsLb()

Member Function Documentation

◆ buildCollectedLevelSet()

◆ checkProperties()

◆ cleanUp()

◆ constructGField()

◆ constructGField_()

◆ createBodyTree()

◆ descendLevelSetValue()

◆ finalizeAdaptation()

◆ finalizeCouplerInit()

◆ finalizeSubCoupleInit()

◆ init()

◆ initData()

◆ initializeSolidDomain()

◆ interpolateLevelSet()

◆ interpolateLsLb()

◆ lb2lsId()

◆ lb2lsIdParent()

◆ ls2lbId()

◆ ls2lbIdParent()

◆ noLevelSetFieldData()

◆ postAdaptation()

◆ postCouple()

◆ preCouple()

◆ readProperties()

◆ subCouple()

◆ testCoupling()

◆ transferLevelSetFieldValues()

◆ updateBoundaryCellsFromGField()

◆ updateFlowSolverLevelSet()

◆ updateGeometry()

◆ updateLevelSetFlowSolver()

Member Data Documentation

◆ bndryToBodyMapping

◆ bodyToBndryMapping

◆ m_constructGField

◆ m_lbSolverId

◆ m_lsSolverId

◆ m_maxNoEmbeddedBodies

◆ m_noCorners

◆ m_outsideDefault

◆ 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

◆ m_static_computeBodyProperties_liftStartAngle1

◆ m_static_computeBodyProperties_liftStartAngle2

◆ m_static_computeBodyProperties_mu

◆ m_static_computeBodyProperties_mu2

◆ m_static_computeBodyProperties_normal

◆ m_static_computeBodyProperties_omega

◆ m_static_computeBodyProperties_rotAngle

◆ m_static_computeBodyProperties_Strouhal

◆ m_static_updateLevelSetFlowSolver_firstRun

◆ m_transferBoundingBox

◆ nDim