#include <lbrb.h>

Inheritance diagram for LbRb< nDim, nDist, SysEqn >:

Collaboration diagram for LbRb< nDim, nDist, SysEqn >:

Classes
struct	ConversionFactors

struct	Timers

Public Types
using	CouplingLb = CouplingLB< nDim, nDist, SysEqn >

using	CouplingRb = CouplingRigidBodies< nDim >

using	LbSolver = LbSolverDxQy< nDim, nDist, SysEqn >

using	RBodies = typename CouplingRb::RBodies

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 Types inherited from CouplingRigidBodies< nDim >
using	RBodies = RigidBodies< nDim >

Public Member Functions
	LbRb (MInt couplingId, LbSolver lb, RBodies rb)
	C'tor for the lattice Boltzmann rigid bodies coupler. More...

	~LbRb ()
	D'tor for the lattice Boltzmann rigid bodies coupler. More...

void	init () override
	Initialize solver data needed for this coupling. More...

void	finalizeSubCoupleInit (MInt) override

void	finalizeCouplerInit () override

void	preCouple (MInt) override

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

void	postCouple (MInt) override

void	postAdaptation () override

void	finalizeAdaptation (const MInt solverId) override
	Coupling between solver substeps. More...

void	cleanUp () override

void	reinitAfterBalance ()

void	getDomainDecompositionInformation (std::vector< std::pair< MString, MInt > > &) override

void	initializeSolidDomain ()
	Initialize cells which are inisde the solid domain or just entered the fluid domain. More...

void	getBodyVelocity (const MInt body, MFloat *const velocity)
	Get body velocity converted to LB units. More...

void	getBodyAngularVelocity (const MInt body, MFloat *const angularVelocity)
	Get angular body velocity converted to LB units. More...

void	createBodyTree ()

void	bc3060 (MInt)

void	transferLevelSetFieldValues (MBool)

void	computeGCellTimeStep ()

MInt	returnNoActiveCorners (MInt)

void	returnLevelSetSignForFluidCellCorners (MInt, MIntScratchSpace *levelSetCornerSigns, MInt set=0)

MInt	noLevelSetFieldData ()

MInt	returnLevelSetSignForFluidFaceCentroid (MInt, MInt, MInt set=0)

void	constructGField ()
	Dispatch function for the body specific construction of the level set. More...

template<MInt bodyType>
void	constructGField_ ()
	Constructs the level-set field after each time step. More...

template<MInt bodyType>
void	descendLevelSetValue (const MInt cellId, const MInt *bodyId, const MInt bodyCnt)
	Descend the level set value to the cells children. More...

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

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

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

MFloat	a_initTemperatureKelvin (const MInt id=0) const

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

MInt	a_isThermal (const MInt id=0) const

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

MFloat	a_Ma (const MInt id=0) const

MbCellCollector &	a_mbCell (const MInt id=0)

MInt	a_noCells (const MInt id=0) const

MInt	a_noDistributions (const MInt id=0) const

MInt	a_noEmbeddedBodiesLB (const MInt id=0) const

MInt	a_noLbCells (const MInt id=0) const

MInt	a_noLevelSetsMb (const MInt id=0) const

MInt	a_noVariables (const MInt id=0) const

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

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

MInt	a_pvrho (const MInt id=0) const

MInt	a_pvt (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

MFloat	a_Re (const MInt id=0) const

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

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

LbBndCnd &	lbBndCnd (const MInt id=0)

solverType &	lbSolver (const MInt solverId=0) const

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

MInt	noMinCells (const MInt id=0) const

MInt	a_noCollectorBodies () const

RBodies &	bodies () const

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 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 CouplingRigidBodies< nDim >
	CouplingRigidBodies (const MInt couplingId, RBodies *solver)

MInt	a_noEmbeddedBodies () const

MInt	a_noCollectorBodies () const

Public Attributes
MInt	m_lbSolverId {}

MString	m_timerType

maia::coupling::Mapping	bndryToBodyMapping

maia::coupling::Mapping	bodyToBndryMapping

MFloat **	forces = nullptr

MBool	m_static_updateLevelSetFlowSolver_firstRun = true

Static Public Attributes
static constexpr MInt	m_noCorners = (nDim == 2) ? 4 : 8

static constexpr MInt	nRot = (nDim == 3) ? 3 : 1

Private Member Functions
void	initData ()
	Initializes coupler specific data. More...

void	initTimers ()
	Creates all timers and subtimers. More...

void	averageTimer ()

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

void	readProperties ()
	Read all relevant properties. More...

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

Private Attributes
ConversionFactors	conversionRbLb

ConversionFactors	conversionLbRb

std::array< MInt, Timers::_count >	m_timers

Additional Inherited Members
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)

Protected Member Functions inherited from Coupling
MFloat	returnLoadRecord () const

MFloat	returnIdleRecord () const

Protected Member Functions inherited from CouplingRigidBodies< nDim >
RBodies &	bodies () const

Detailed Description

template<MInt nDim, MInt nDist, class SysEqn>
class LbRb< nDim, nDist, SysEqn >

Definition at line 28 of file lbrb.h.

Member Typedef Documentation

◆ CouplingLb

template<MInt nDim, MInt nDist, class SysEqn >

using LbRb< nDim, nDist, SysEqn >::CouplingLb = CouplingLB<nDim, nDist, SysEqn>

Definition at line 31 of file lbrb.h.

◆ CouplingRb

template<MInt nDim, MInt nDist, class SysEqn >

using LbRb< nDim, nDist, SysEqn >::CouplingRb = CouplingRigidBodies<nDim>

Definition at line 32 of file lbrb.h.

◆ LbSolver

template<MInt nDim, MInt nDist, class SysEqn >

using LbRb< nDim, nDist, SysEqn >::LbSolver = LbSolverDxQy<nDim, nDist, SysEqn>

Definition at line 34 of file lbrb.h.

◆ RBodies

template<MInt nDim, MInt nDist, class SysEqn >

using LbRb< nDim, nDist, SysEqn >::RBodies = typename CouplingRb::RBodies

Definition at line 35 of file lbrb.h.

Constructor & Destructor Documentation

◆ LbRb()

template<MInt nDim, MInt nDist, class SysEqn >

LbRb< nDim, nDist, SysEqn >::LbRb	(	MInt	couplingId,
		LbSolver *	lb,
		RBodies *	rb
	)

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Parameters

couplingId[in]	Unique id to identify the couplre
lb[in]	Pointer to the lattice Botzmann solver
rb[in]	Pointer to the rigid bodies solver

Definition at line 30 of file lbrb.cpp.

  : Coupling(couplingId), CouplingLb(couplingId, lb), CouplingRb(couplingId, rb) {
  TRACE();
 
  // Init timers as the first action
  initTimers();
 
  initData();
  readProperties();
  checkProperties();
  updateGeometry();
 
  RECORD_TIMER_STOP(m_timers[Timers::Constructor]);
}

◆ ~LbRb()

template<MInt nDim, MInt nDist, class SysEqn >

LbRb< nDim, nDist, SysEqn >::~LbRb

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Definition at line 51 of file lbrb.cpp.

                                 {
  TRACE();
 
  RECORD_TIMER_STOP(m_timers[Timers::Class]);
 
  averageTimer();
}

Member Function Documentation

◆ a_associatedBodyIdsMb() [1/2]

template<MInt nDim, MInt nDist, class SysEqn >

MInt & CouplingLB< nDim, nDist, SysEqn >::a_associatedBodyIdsMb	(	const MInt	cellId,
		const MInt	set,
		const MInt	id = `0`
	)

inline

Definition at line 498 of file coupling.h.

                                                                                    {
    return lbSolver(id).a_associatedBodyIds(cellId, set);
  }

◆ a_associatedBodyIdsMb() [2/2]

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_associatedBodyIdsMb	(	const MInt	cellId,
		const MInt	set,
		const MInt	id = `0`
	)		const

inline

Definition at line 502 of file coupling.h.

                                                                                         {
    return lbSolver(id).a_associatedBodyIds(cellId, set);
  }

◆ a_cellLengthAtLevel()

template<MInt nDim, MInt nDist, class SysEqn >

MFloat CouplingLB< nDim, nDist, SysEqn >::a_cellLengthAtLevel	(	MInt	level,
		const MInt	id = `0`
	)

inline

Definition at line 518 of file coupling.h.

518{ return lbSolver(id).grid().cellLengthAtLevel(level); }

◆ a_childId()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_childId	(	const MInt	cellId,
		const MInt	child,
		const MInt	id = `0`
	)

inline

Definition at line 508 of file coupling.h.

                                                                         {
    return lbSolver(id).c_childId(cellId, child);
  }

◆ a_initTemperatureKelvin()

template<MInt nDim, MInt nDist, class SysEqn >

MFloat CouplingLB< nDim, nDist, SysEqn >::a_initTemperatureKelvin ( const MInt id = 0 ) const

inline

Definition at line 482 of file coupling.h.

482{ return lbSolver(id).m_initTemperatureKelvin; }

◆ a_isActive()

template<MInt nDim, MInt nDist, class SysEqn >

MBool CouplingLB< nDim, nDist, SysEqn >::a_isActive	(	const MInt	cellId,
		const MInt	id = `0`
	)		const

inline

Definition at line 522 of file coupling.h.

522{ return lbSolver(id).a_isActive(cellId); }

◆ a_isThermal()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_isThermal ( const MInt id = 0 ) const

inline

Definition at line 479 of file coupling.h.

479{ return lbSolver(id).m_isThermal; }

◆ a_levelSetFunctionMb() [1/2]

template<MInt nDim, MInt nDist, class SysEqn >

MFloat & CouplingLB< nDim, nDist, SysEqn >::a_levelSetFunctionMb	(	const MInt	cellId,
		const MInt	set,
		const MInt	id = `0`
	)

inline

Definition at line 490 of file coupling.h.

                                                                                     {
    return lbSolver(id).a_levelSetFunctionMB(cellId, set);
  }

◆ a_levelSetFunctionMb() [2/2]

template<MInt nDim, MInt nDist, class SysEqn >

MFloat CouplingLB< nDim, nDist, SysEqn >::a_levelSetFunctionMb	(	const MInt	cellId,
		const MInt	set,
		const MInt	id = `0`
	)		const

inline

Definition at line 494 of file coupling.h.

                                                                                          {
    return lbSolver(id).a_levelSetFunctionMB(cellId, set);
  }

◆ a_Ma()

template<MInt nDim, MInt nDist, class SysEqn >

MFloat CouplingLB< nDim, nDist, SysEqn >::a_Ma ( const MInt id = 0 ) const

inline

Definition at line 472 of file coupling.h.

472{ return lbSolver(id).m_Ma; }

◆ a_mbCell()

template<MInt nDim, MInt nDist, class SysEqn >

MbCellCollector & CouplingLB< nDim, nDist, SysEqn >::a_mbCell ( const MInt id = 0 )

inline

Definition at line 486 of file coupling.h.

486{ return lbBndCnd(id).m_boundaryCellsMb; }

LbBndCnd::m_boundaryCellsMb

MbCellCollector m_boundaryCellsMb

Definition: lbbndcnd.h:57

LbRb::lbBndCnd

LbBndCnd & lbBndCnd(const MInt id=0)

Definition: coupling.h:461

◆ a_noCells()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_noCells ( const MInt id = 0 ) const

inline

Definition at line 516 of file coupling.h.

516{ return lbSolver(id).grid().noCells(); }

◆ a_noCollectorBodies()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingRigidBodies< nDim >::a_noCollectorBodies ( ) const

inline

Definition at line 591 of file coupling.h.

591{ return bodies().noCollectorBodies(); }

LbRb::bodies

RBodies & bodies() const

Definition: coupling.h:594

◆ a_noDistributions()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_noDistributions ( const MInt id = 0 ) const

inline

Definition at line 481 of file coupling.h.

481{ return lbSolver(id).m_noDistributions; }

LbSolverDxQy::m_noDistributions

static constexpr MInt m_noDistributions

Definition: lbsolverdxqy.h:180

◆ a_noEmbeddedBodiesLB()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_noEmbeddedBodiesLB ( const MInt id = 0 ) const

inline

Definition at line 520 of file coupling.h.

520{ return lbSolver(id).m_noEmbeddedBodies; }

◆ a_noLbCells()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_noLbCells ( const MInt id = 0 ) const

inline

Definition at line 470 of file coupling.h.

470{ return lbSolver(id).a_noCells(); }

◆ a_noLevelSetsMb()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_noLevelSetsMb ( const MInt id = 0 ) const

inline

Definition at line 471 of file coupling.h.

471{ return lbSolver(id).m_noLevelSetsUsedForMb; }

◆ a_noVariables()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_noVariables ( const MInt id = 0 ) const

inline

Definition at line 526 of file coupling.h.

526{ return lbSolver(id).noVariables(); }

◆ a_oldVariable()

template<MInt nDim, MInt nDist, class SysEqn >

MFloat & CouplingLB< nDim, nDist, SysEqn >::a_oldVariable	(	const MInt	cellId,
		const MInt	varId,
		const MInt	id = `0`
	)

inline

Definition at line 533 of file coupling.h.

                                                                                {
    return lbSolver(id).a_oldVariable(cellId, varId);
  }

◆ a_parentId()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_parentId	(	const MInt	cellId,
		const MInt	id = `0`
	)

inline

Definition at line 506 of file coupling.h.

506{ return lbSolver(id).c_parentId(cellId); }

◆ a_pvrho()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_pvrho ( const MInt id = 0 ) const

inline

Definition at line 477 of file coupling.h.

477{ return lbSolver(id).PV->RHO; }

◆ a_pvt()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_pvt ( const MInt id = 0 ) const

inline

Definition at line 478 of file coupling.h.

478{ return lbSolver(id).PV->T; }

◆ a_pvu()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_pvu ( const MInt id = 0 ) const

inline

Definition at line 474 of file coupling.h.

474{ return lbSolver(id).PV->U; }

◆ a_pvv()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_pvv ( const MInt id = 0 ) const

inline

Definition at line 475 of file coupling.h.

475{ return lbSolver(id).PV->V; }

◆ a_pvw()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::a_pvw ( const MInt id = 0 ) const

inline

Definition at line 476 of file coupling.h.

476{ return lbSolver(id).PV->W; }

◆ a_Re()

template<MInt nDim, MInt nDist, class SysEqn >

MFloat CouplingLB< nDim, nDist, SysEqn >::a_Re ( const MInt id = 0 ) const

inline

Definition at line 473 of file coupling.h.

473{ return lbSolver(id).m_Re; }

◆ a_variable()

template<MInt nDim, MInt nDist, class SysEqn >

MFloat & CouplingLB< nDim, nDist, SysEqn >::a_variable	(	const MInt	cellId,
		const MInt	varId,
		const MInt	id = `0`
	)

inline

Definition at line 528 of file coupling.h.

                                                                             {
    return lbSolver(id).a_variable(cellId, varId);
  }

◆ a_wasActive()

template<MInt nDim, MInt nDist, class SysEqn >

MBool CouplingLB< nDim, nDist, SysEqn >::a_wasActive	(	const MInt	cellId,
		const MInt	id = `0`
	)		const

inline

Definition at line 524 of file coupling.h.

524{ return lbSolver(id).a_wasActive(cellId); }

◆ averageTimer()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::averageTimer

private

Definition at line 580 of file lbrb.cpp.

                                             {
  TRACE();
  if(!lbSolver().grid().isActive()) return;
 
  // Get timer operation
  m_timerType = "max";
  m_timerType = Context::getSolverProperty<MString>("timerType", lbSolver().solverId(), AT_, &m_timerType);
 
  // 0) map timer ids for safety
  const MInt noTimers = m_timers.size();
 
  // 1) fill buffer with local timer values
  std::vector<MFloat> timerValues_;
  timerValues_.reserve(noTimers);
  for(MInt i = 0; i < noTimers; i++) {
    timerValues_.emplace_back(RETURN_TIMER_TIME(m_timers[i]));
  }
 
  // 2) collect values from all ranks
  if(m_timerType == "average") {
    MPI_Allreduce(MPI_IN_PLACE, timerValues_.data(), noTimers, maia::type_traits<MFloat>::mpiType(), MPI_SUM,
                  lbSolver().mpiComm(), AT_, "MPI_IN_PLACE", "timerValues_");
  } else {
    MPI_Allreduce(MPI_IN_PLACE, timerValues_.data(), noTimers, maia::type_traits<MFloat>::mpiType(), MPI_MAX,
                  lbSolver().mpiComm(), AT_, "MPI_IN_PLACE", "timerValues_");
  }
 
  // 3) perform averaging on timer and4) set new timer values
  if(m_timerType == "average") {
    const MInt noDomains_ = lbSolver().noDomains();
    for(MInt i = 0; i < noTimers; i++) {
      const MFloat meanValue = timerValues_[i] / noDomains_;
      SET_RECORD(m_timers[i], meanValue);
    }
  } else {
    for(MInt i = 0; i < noTimers; i++) {
      SET_RECORD(m_timers[i], timerValues_[i]);
    }
  }
}

◆ bc3060()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::bc3060 ( MInt )

◆ bodies()

template<MInt nDim, MInt nDist, class SysEqn >

RBodies & CouplingRigidBodies< nDim >::bodies ( ) const

inline

Definition at line 594 of file coupling.h.

594{ return *m_rigidBodies; }

CouplingRigidBodies::m_rigidBodies

RBodies * m_rigidBodies

Definition: coupling.h:597

◆ checkProperties()

template<MInt nDim, MInt nDist, class SysEqn >

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

privatevirtual

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Implements Coupling.

Definition at line 152 of file lbrb.cpp.

                                                {
  TRACE();
 
  lbSolver().m_noLevelSetsUsedForMb = 1;
  lbSolver().m_maxNoSets = 1;
}

◆ cleanUp()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::cleanUp ( )

inlineoverridevirtual

Implements Coupling.

Definition at line 152 of file lbrb.h.

152{};

◆ computeGCellTimeStep()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::computeGCellTimeStep ( )

◆ constructGField()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::constructGField

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Definition at line 442 of file lbrb.cpp.

                                                {
  if(bodies().a_bodyType() == 1)
    constructGField_<1>();
  else if(bodies().a_bodyType() == 2)
    constructGField_<2>();
  else if(bodies().a_bodyType() == 3)
    constructGField_<3>();
  else if(bodies().a_bodyType() == 4)
    constructGField_<4>();
  else if(bodies().a_bodyType() == 7)
    constructGField_<7>();
  else
    mTerm(1, AT_, "Body type not implemented!");
}

◆ constructGField_()

template<MInt nDim, MInt nDist, class SysEqn >

template<MInt bodyType>

void LbRb< nDim, nDist, SysEqn >::constructGField_

Author: Lennart Schneiders, Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Template Parameters

bodyType Body type for which the level set is constructed

Definition at line 466 of file lbrb.cpp.

                                                 {
  TRACE();
 
  RECORD_TIMER_START(m_timers[Timers::Preparation]);
 
  /* bandwidth at "level" has the no. of cells at "level",
  whose total width equates to all summed up bands up to "level" */
  MInt level = lbSolver().maxUniformRefinementLevel() - 1;
  MFloat m_bodyDistThreshold = 0.0;
  if(lbSolver().m_adaptation) {
    m_bodyDistThreshold = lbSolver().m_bandWidth[level] * a_cellLengthAtLevel(level);
    // m_bodyDistThreshold =  4 * a_cellLengthAtLevel(lbSolver().maxLevel());
  } else {
    m_bodyDistThreshold = 2 * a_cellLengthAtLevel(lbSolver().maxLevel());
  }
 
  // To be changed later ~jv
  const MInt m_noLevelSetsUsedForMb = 1;
 
  // Reset all sets
  for(MInt cellId = 0; cellId < a_noCells(); cellId++) {
    for(MInt set = 0; set < m_noLevelSetsUsedForMb; set++) {
      a_associatedBodyIdsMb(cellId, set) = -1;
      a_levelSetFunctionMb(cellId, set) = m_bodyDistThreshold + 1e-14;
    }
  }
 
  const MInt noRelevantBodies = a_noCollectorBodies();
 
  if(noRelevantBodies == 0) {
    RECORD_TIMER_STOP(m_timers[Timers::Preparation]);
    return;
  }
 
  RECORD_TIMER_STOP(m_timers[Timers::Preparation]);
 
  // construct vector with connecting Bodies
  std::vector<MInt> collectorBodyIds(a_noCollectorBodies());
  for(MInt i = 0; i < a_noCollectorBodies(); i++) {
    collectorBodyIds[i] = i;
  }
 
  for(MInt i = 0; i < noMinCells(); i++) {
    const MInt cellId = minCell(i);
    const MFloat minDist = m_bodyDistThreshold + 1e-14;
 
    for(MInt set = 0; set < m_noLevelSetsUsedForMb; set++) {
      a_associatedBodyIdsMb(cellId, set) = -1;
      a_levelSetFunctionMb(cellId, set) = minDist;
    }
 
    if(a_noCollectorBodies() > 0) {
      descendLevelSetValue<bodyType>(cellId, collectorBodyIds.data(), noRelevantBodies);
    }
 
    MInt parentId = a_parentId(cellId);
    while(parentId > -1) {
      for(MInt set = 0; set < m_noLevelSetsUsedForMb; set++) {
        a_levelSetFunctionMb(parentId, set) = a_levelSetFunctionMb(cellId, set);
        a_associatedBodyIdsMb(parentId, set) = a_associatedBodyIdsMb(cellId, set);
      }
      parentId = a_parentId(parentId);
    }
  }
}

◆ createBodyTree()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::createBodyTree ( )

◆ descendLevelSetValue()

template<MInt nDim, MInt nDist, class SysEqn >

template<MInt bodyType>

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

This function is used recursively starting from the min level cells. The descend is skipped if the distance on the lower level is greater than a given threshold.

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Parameters

[in]	cellId	Cell id of the coarse cell
[in]	bodyIds	List of all relevant body ids
[in]	bodyCnt	Number of bodies

Definition at line 546 of file lbrb.cpp.

                                                                                                                     {
  const MFloat minLevelThreshold = a_cellLengthAtLevel(lbSolver().a_level(cellId));
 
  MBool skipDescend = true;
 
  for(MInt b = 0; b < bodyCnt; b++) {
    const MInt k = bodyIds[b];
    const MInt set = 0;
 
    MFloat dist = bodies().template getDistance<bodyType>(&lbSolver().a_coordinate(cellId, 0), k);
 
    if(dist < minLevelThreshold) {
      skipDescend = false;
    }
 
    if(fabs(dist) < fabs(a_levelSetFunctionMb(cellId, set)) || (dist < F0)) {
      a_levelSetFunctionMb(cellId, set) = dist;
      a_associatedBodyIdsMb(cellId, set) = k;
    }
  }
 
  if(skipDescend || !lbSolver().grid().tree().hasChildren(cellId)) {
    return;
  }
 
  // Recursively descend via children
  constexpr MInt maxNoChildren = nDim == 3 ? 8 : 4;
  for(MInt child = 0; child < maxNoChildren; child++) {
    if(a_childId(cellId, child) < 0) continue;
    descendLevelSetValue<bodyType>(a_childId(cellId, child), bodyIds, bodyCnt);
  }
}

◆ finalizeAdaptation()

template<MInt nDim, MInt nDist, class SysEqn >

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

overridevirtual

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Parameters

[in]	rs	Current recipe step
[in]	sid	Current solver id after which this function is called
[in]	solverCompleted	Current completion status of all solvers

Reimplemented from Coupling.

Definition at line 213 of file lbrb.cpp.

                                                                      {
  // LB bnd cnd is restarted after adaptation and needs to be initialized again
  if(solverId == lbSolver().solverId()) {
    lbSolver().initializeMovingBoundaries();
    lbBndCnd().initializeBndMovingBoundaries();
  }
}

◆ finalizeCouplerInit()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::finalizeCouplerInit ( )

inlineoverridevirtual

Implements Coupling.

Definition at line 143 of file lbrb.h.

143{};

◆ finalizeSubCoupleInit()

template<MInt nDim, MInt nDist, class SysEqn >

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

inlineoverridevirtual

Implements Coupling.

Definition at line 142 of file lbrb.h.

142{};

◆ getBodyAngularVelocity()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::getBodyAngularVelocity	(	const MInt	body,
		MFloat *const	angularVelocity
	)

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Parameters

[in]	body	Body id
[in]	bodyVelocity	Angular body velocity

Definition at line 331 of file lbrb.cpp.

                                                                                                     {
  for(MInt n = 0; n < nRot; n++) {
    angularVelocity[n] = bodies().a_angularVelocity(body, n) / conversionRbLb.time;
  }
}

◆ getBodyVelocity()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::getBodyVelocity	(	const MInt	body,
		MFloat *const	bodyVelocity
	)

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Parameters

[in]	body	Body id
[in]	bodyVelocity	Body velocity

Definition at line 316 of file lbrb.cpp.

                                                                                           {
  for(MInt n = 0; n < nDim; n++) {
    bodyVelocity[n] = bodies().a_bodyVelocity(body, n) * conversionRbLb.velocity;
  }
}

◆ getDomainDecompositionInformation()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::getDomainDecompositionInformation ( std::vector< std::pair< MString, MInt > > & )

inlineoverride

Definition at line 154 of file lbrb.h.

154{};

◆ init()

template<MInt nDim, MInt nDist, class SysEqn >

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

overridevirtual

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Implements Coupling.

Definition at line 173 of file lbrb.cpp.

                                     {
  TRACE();
 
  // LB solver
  lbSolver().initializeMovingBoundaries();
  lbBndCnd().initializeBndMovingBoundaries();
 
  // RB solver
  bodies().setTimestep(conversionLbRb.time);
 
  // Construct level set for initial adaptation
  RECORD_TIMER_START(m_timers[Timers::CouplePostRb]);
  RECORD_TIMER_START(m_timers[Timers::ConstructGField]);
  constructGField();
  RECORD_TIMER_STOP(m_timers[Timers::ConstructGField]);
  RECORD_TIMER_STOP(m_timers[Timers::CouplePostRb]);
}

◆ initData()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::initData

private

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Definition at line 105 of file lbrb.cpp.

                                         {
  TRACE();
 
  // update RB grid meta data
  bodies().updateMaxLevel(lbSolver().maxLevel());
 
  // Conversion factors
  const MFloat dx = a_cellLengthAtLevel(lbSolver().maxLevel());
 
  conversionLbRb.length = dx;
  conversionRbLb.length = 1 / conversionLbRb.length;
 
  conversionLbRb.velocity = sqrt(3) / a_Ma();
  conversionRbLb.velocity = 1 / conversionLbRb.velocity;
 
  conversionLbRb.time = conversionLbRb.length / conversionLbRb.velocity;
  conversionRbLb.time = 1 / conversionLbRb.time;
 
  conversionLbRb.force = POW2(conversionLbRb.velocity) * POW2(conversionLbRb.length);
  conversionRbLb.force = 1 / conversionLbRb.force;
 
  conversionLbRb.torque = conversionLbRb.force * conversionLbRb.length;
  conversionRbLb.torque = 1 / conversionLbRb.torque;
 
  bodies().setTimestep(conversionLbRb.time);
 
#ifndef NDEBUG
  if(!lbSolver().domainId()) {
    std::cout << "CONVERSION" << std::endl
              << "LENGTH LB RB " << conversionLbRb.length << std::endl
              << "LENGTH RB LB " << conversionRbLb.length << std::endl
              << "VEL LB RB " << conversionLbRb.velocity << std::endl
              << "VEL RB LB " << conversionRbLb.velocity << std::endl
              << "TIME LB RB " << conversionLbRb.time << std::endl
              << "TIME RB LB " << conversionRbLb.time << std::endl
              << "FORCE LB RB " << conversionLbRb.force << std::endl
              << "FORCE RB LB " << conversionRbLb.force << std::endl;
  }
#endif
}

◆ initializeSolidDomain()

template<MInt nDim, MInt nDist, class SysEqn >

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

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Definition at line 357 of file lbrb.cpp.

                                                      {
  TRACE();
 
  if(lbSolver().noNeighborDomains() > 0) {
    lbSolver().exchangeOldDistributions();
  }
 
  MFloatScratchSpace bodyVelocities(a_noCollectorBodies(), nDim, AT_, "bodyVelocities");
 
  // Get the body Velocity for each embedded body
  for(MInt body = 0; body < a_noCollectorBodies(); body++) {
    getBodyVelocity(body, &bodyVelocities(body, 0));
  }
 
  for(MInt i = 0; i < a_noCells(); i++) {
    // Regular fluid cell
    if(a_isActive(i) && a_wasActive(i)) {
      continue;
    }
 
    // Regular Solid cell
    if(!a_isActive(i)) {
      if(!lbSolver().c_isLeafCell(i)) {
        continue;
      }
 
      // determine the 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;
        }
      }
 
      // if body was deleted, the halo cells that were occupied by the body in the previous timeStep, must be refilled
      // now
      if(bodyId == -1 && lbSolver().a_isHalo(i) && bodies().m_bodyWasDeleted) {
        lbBndCnd().refillEmergedCell(i);
        lbSolver().a_isActive(i) = 1;
        continue;
      }
 
      ASSERT(bodyId > -1, "No valid bodyId for solid cell! (bodyId=" << bodyId << ") (" << lbSolver().a_coordinate(i, 0)
                                                                     << " " << lbSolver().a_coordinate(i, 1) << ") ("
                                                                     << i << " " << lbSolver().c_globalId(i) << ")"
                                                                     << " isHalo " << lbSolver().a_isHalo(i));
 
      // the Velocity of the deactivated cell is set to the body velocity
      // the Density is set to 1.0
      if((bodyId >= 0) && (bodyId < a_noCollectorBodies()) && (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);
        }
      } else {
        for(MInt j = 0; j < a_noVariables(); j++) {
          a_variable(i, j) = F0;
          a_oldVariable(i, j) = F0;
        }
      }
      a_variable(i, a_pvrho()) = 1.0;
      a_oldVariable(i, a_pvrho()) = 1.0;
 
      if(a_isThermal()) {
        a_variable(i, a_pvt()) = a_initTemperatureKelvin();
      }
 
      // Distributions are not set since they are not used for solid nodes
    } else {
      // New fluid cell
      lbBndCnd().refillEmergedCell(i);
    }
  }
  bodies().m_bodyWasDeleted = false;
}

◆ initTimers()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::initTimers

private

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Definition at line 65 of file lbrb.cpp.

                                           {
  TRACE();
 
  // Create timer group and coupler timer, and start the timer
  NEW_TIMER_GROUP_NOCREATE(m_timers[Timers::TimerGroup], "LbRb");
  NEW_TIMER_NOCREATE(m_timers[Timers::Class], "Total object lifetime", m_timers[Timers::TimerGroup]);
  RECORD_TIMER_START(m_timers[Timers::Class]);
 
  // Create and start constructor timer
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::Constructor], "Constructor", m_timers[Timers::Class]);
  RECORD_TIMER_START(m_timers[Timers::Constructor]);
 
  // Create Couple timers
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::CouplePostRb], "Couple post RB", m_timers[Timers::Class]);
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::ConstructGField], "ConstructGField", m_timers[Timers::CouplePostRb]);
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::Preparation], "Prepare construct", m_timers[Timers::ConstructGField]);
 
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::FindBoundaryCells], "FindBoundaryCells", m_timers[Timers::CouplePostRb]);
 
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::PreCouple], "Pre Couple", m_timers[Timers::FindBoundaryCells]);
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::FindBoundaryMapping], "Find Boundary Mapping",
                         m_timers[Timers::FindBoundaryCells]);
 
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::InitSolidDomain], "InitSolidDomain", m_timers[Timers::CouplePostRb]);
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::SetBoundaryVelocity], "SetBoundaryVelocity", m_timers[Timers::CouplePostRb]);
 
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::CouplePostLb], "Couple post LB", m_timers[Timers::Class]);
 
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::CreateComm], "CreateComm", m_timers[Timers::CouplePostLb]);
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::ApplyBC], "ApplyBC", m_timers[Timers::CouplePostLb]);
  NEW_SUB_TIMER_NOCREATE(m_timers[Timers::ComputeBodyForces], "ComputeBodyForces", m_timers[Timers::CouplePostLb]);
}

◆ lbBndCnd()

template<MInt nDim, MInt nDist, class SysEqn >

LbBndCnd & CouplingLB< nDim, nDist, SysEqn >::lbBndCnd ( const MInt id = 0 )

inline

Definition at line 461 of file coupling.h.

461{ return *lbSolver(id).m_bndCnd; }

LbSolverDxQy::m_bndCnd

LbBndCnd * m_bndCnd

Pointers for the Boundary Conditions, for flow solving.

Definition: lbsolverdxqy.h:332

◆ lbSolver()

template<MInt nDim, MInt nDist, class SysEqn >

solverType & CouplingLB< nDim, nDist, SysEqn >::lbSolver ( const MInt solverId = 0 ) const

inline

Definition at line 454 of file coupling.h.

                                                      {
    ASSERT(solverId < noSolvers(), "Invalid solverId " + std::to_string(solverId) + "/" + std::to_string(noSolvers()));
    return *m_lbSolvers[solverId];
  }

◆ minCell()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::minCell	(	const MInt	index,
		const MInt	id = `0`
	)		const

inline

Definition at line 512 of file coupling.h.

512{ return lbSolver(id).grid().minCell(index); }

◆ noLevelSetFieldData()

template<MInt nDim, MInt nDist, class SysEqn >

MInt LbRb< nDim, nDist, SysEqn >::noLevelSetFieldData ( )

◆ noMinCells()

template<MInt nDim, MInt nDist, class SysEqn >

MInt CouplingLB< nDim, nDist, SysEqn >::noMinCells ( const MInt id = 0 ) const

inline

Definition at line 514 of file coupling.h.

514{ return lbSolver(id).grid().noMinCells(); }

◆ postAdaptation()

template<MInt nDim, MInt nDist, class SysEqn >

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

overridevirtual

Reimplemented from Coupling.

Definition at line 192 of file lbrb.cpp.

                                               {
  // Update conversion factors
  initData();
  // Construct level set to initialize newly refined cells
  RECORD_TIMER_START(m_timers[Timers::CouplePostRb]);
  RECORD_TIMER_START(m_timers[Timers::ConstructGField]);
  constructGField();
  RECORD_TIMER_STOP(m_timers[Timers::ConstructGField]);
  RECORD_TIMER_STOP(m_timers[Timers::CouplePostRb]);
}

◆ postCouple()

template<MInt nDim, MInt nDist, class SysEqn >

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

inlineoverridevirtual

Implements Coupling.

Definition at line 147 of file lbrb.h.

147{};

◆ preCouple()

template<MInt nDim, MInt nDist, class SysEqn >

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

inlineoverridevirtual

Implements Coupling.

Definition at line 145 of file lbrb.h.

145{};

◆ readProperties()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::readProperties

privatevirtual

Author: Julian Vorspohl j.vor.nosp@m.spoh.nosp@m.l@aia.nosp@m..rwt.nosp@m.h-aac.nosp@m.hen..nosp@m.de

Implements Coupling.

Definition at line 165 of file lbrb.cpp.

165{}

◆ reinitAfterBalance()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::reinitAfterBalance ( )

inlinevirtual

Reimplemented from Coupling.

Definition at line 153 of file lbrb.h.

153{};

◆ returnLevelSetSignForFluidCellCorners()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::returnLevelSetSignForFluidCellCorners	(	MInt	,
		MIntScratchSpace *	levelSetCornerSigns,
		MInt	set = `0`
	)

◆ returnLevelSetSignForFluidFaceCentroid()

template<MInt nDim, MInt nDist, class SysEqn >

MInt LbRb< nDim, nDist, SysEqn >::returnLevelSetSignForFluidFaceCentroid	(	MInt	,
		MInt	,
		MInt	set = `0`
	)

◆ returnNoActiveCorners()

template<MInt nDim, MInt nDist, class SysEqn >

MInt LbRb< nDim, nDist, SysEqn >::returnNoActiveCorners ( MInt )

◆ subCouple()

template<MInt nDim, MInt nDist, class SysEqn >

void LbRb< nDim, nDist, SysEqn >::subCouple	(	MInt	,
		MInt	sid,
		std::vector< MBool > &	solverCompleted
	)

overridevirtual

Implements Coupling.

Definition at line 222 of file lbrb.cpp.

                                                                                             {
  // Do nothing if both solver already finished
  if(solverCompleted[lbSolver().solverId()] && solverCompleted[bodies().solverId()]) {
    return;
  }
 
  if(sid == bodies().solverId()) {
    // POST RB
    RECORD_TIMER_START(m_timers[Timers::CouplePostRb]);
    RECORD_TIMER_START(m_timers[Timers::ConstructGField]);
    constructGField();
    RECORD_TIMER_STOP(m_timers[Timers::ConstructGField]);
 
    std::vector<MInt> maxGCellLevels(lbSolver().m_maxNoSets);
    for(MInt set = 0; set < lbSolver().m_maxNoSets; set++) {
      maxGCellLevels[set] = lbSolver().maxLevel();
    }
 
    RECORD_TIMER_START(m_timers[Timers::FindBoundaryCells]);
    RECORD_TIMER_START(m_timers[Timers::PreCouple]);
    lbSolver().preCoupleLs(maxGCellLevels);
    RECORD_TIMER_STOP(m_timers[Timers::PreCouple]);
 
    RECORD_TIMER_START(m_timers[Timers::FindBoundaryMapping]);
    lbSolver().createBndryToBodyMapping(bndryToBodyMapping, bodyToBndryMapping);
    RECORD_TIMER_STOP(m_timers[Timers::FindBoundaryMapping]);
    RECORD_TIMER_STOP(m_timers[Timers::FindBoundaryCells]);
 
    RECORD_TIMER_START(m_timers[Timers::InitSolidDomain]);
    initializeSolidDomain();
    RECORD_TIMER_STOP(m_timers[Timers::InitSolidDomain]);
 
    RECORD_TIMER_START(m_timers[Timers::SetBoundaryVelocity]);
 
    struct {
      MFloat velocity;
      MFloat angularVelocity;
      MFloat length;
      MFloat temperature;
    } conversion{conversionRbLb.velocity, conversionLbRb.time, 1.0 / a_cellLengthAtLevel(lbSolver().maxLevel()), 1.0};
 
    maia::coupling::setBoundaryVelocity<nDim>(bodies(), a_mbCell(), bodyToBndryMapping, conversion);
 
    RECORD_TIMER_STOP(m_timers[Timers::SetBoundaryVelocity]);
    RECORD_TIMER_STOP(m_timers[Timers::CouplePostRb]);
 
  } else if(sid == lbSolver().solverId()) {
    // POST LB
    RECORD_TIMER_START(m_timers[Timers::CouplePostLb]);
    RECORD_TIMER_START(m_timers[Timers::CreateComm]);
    lbBndCnd().createMBComm();
    RECORD_TIMER_STOP(m_timers[Timers::CreateComm]);
 
    RECORD_TIMER_START(m_timers[Timers::ApplyBC]);
    lbBndCnd().postCouple();
    RECORD_TIMER_STOP(m_timers[Timers::ApplyBC]);
 
    RECORD_TIMER_START(m_timers[Timers::ComputeBodyForces]);
 
    struct {
      MFloat force;
      MFloat length;
    } const conversion{conversionLbRb.force, 1.0};
 
    if(a_noCollectorBodies() > 0) {
      maia::coupling::setBoundaryForceAndTorque<nDim>(a_mbCell(), bodies(), bndryToBodyMapping, conversion);
    }
 
#ifndef NDEBUG
    for(MInt b = 0; b < a_noCollectorBodies(); b++) {
      for(MInt n = 0; n < nDim; n++) {
        if(std::isnan(bodies().a_bodyForce(b, n))) {
          std::cout << "SUM FORCE IS NAN AFTER SETBOUNDANDFORCE AT TIMESETEP: " << lbSolver().getCurrentTimeStep()
                    << std::endl;
          // TERMM(1, "SUM FORCE AFTER IS NAN");
        }
      }
    }
#endif
 
    RECORD_TIMER_STOP(m_timers[Timers::ComputeBodyForces]);
    RECORD_TIMER_STOP(m_timers[Timers::CouplePostLb]);
  }
}

◆ transferLevelSetFieldValues()

template<MInt nDim, MInt nDist, class SysEqn >

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

◆ updateGeometry()

template<MInt nDim, MInt nDist, class SysEqn >

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

private

Author: Tim Wegmann

Definition at line 343 of file lbrb.cpp.

                                               {
  TRACE();
 
  lbSolver().m_noEmbeddedBodies = 1;
  lbSolver().m_geometryIntersection->m_noLevelSetsUsedForMb = 1;
  lbSolver().m_geometryIntersection->m_noEmbeddedBodies = bodies().size();
}

Member Data Documentation

◆ bndryToBodyMapping

template<MInt nDim, MInt nDist, class SysEqn >

maia::coupling::Mapping LbRb< nDim, nDist, SysEqn >::bndryToBodyMapping

Definition at line 44 of file lbrb.h.

◆ bodyToBndryMapping

template<MInt nDim, MInt nDist, class SysEqn >

maia::coupling::Mapping LbRb< nDim, nDist, SysEqn >::bodyToBndryMapping

Definition at line 45 of file lbrb.h.

◆ conversionLbRb

template<MInt nDim, MInt nDist, class SysEqn >

ConversionFactors LbRb< nDim, nDist, SysEqn >::conversionLbRb

private

Definition at line 98 of file lbrb.h.

◆ conversionRbLb

template<MInt nDim, MInt nDist, class SysEqn >

ConversionFactors LbRb< nDim, nDist, SysEqn >::conversionRbLb

private

Definition at line 97 of file lbrb.h.

◆ forces

template<MInt nDim, MInt nDist, class SysEqn >

MFloat** LbRb< nDim, nDist, SysEqn >::forces = nullptr

Definition at line 46 of file lbrb.h.

◆ m_lbSolverId

template<MInt nDim, MInt nDist, class SysEqn >

MInt LbRb< nDim, nDist, SysEqn >::m_lbSolverId {}

Definition at line 41 of file lbrb.h.

◆ m_noCorners

template<MInt nDim, MInt nDist, class SysEqn >

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

staticconstexpr

Definition at line 48 of file lbrb.h.

◆ m_static_updateLevelSetFlowSolver_firstRun

template<MInt nDim, MInt nDist, class SysEqn >

MBool LbRb< nDim, nDist, SysEqn >::m_static_updateLevelSetFlowSolver_firstRun = true

Definition at line 51 of file lbrb.h.

◆ m_timers

template<MInt nDim, MInt nDist, class SysEqn >

std::array<MInt, Timers::_count> LbRb< nDim, nDist, SysEqn >::m_timers

private

Definition at line 126 of file lbrb.h.

◆ m_timerType

template<MInt nDim, MInt nDist, class SysEqn >

MString LbRb< nDim, nDist, SysEqn >::m_timerType

Definition at line 42 of file lbrb.h.

◆ nRot

template<MInt nDim, MInt nDist, class SysEqn >

constexpr MInt LbRb< nDim, nDist, SysEqn >::nRot = (nDim == 3) ? 3 : 1

staticconstexpr

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

Classes

Public Types

Public Member Functions

Public Attributes

Static Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ CouplingLb

◆ CouplingRb

◆ LbSolver

◆ RBodies

Constructor & Destructor Documentation

◆ LbRb()

◆ ~LbRb()

Member Function Documentation

◆ a_associatedBodyIdsMb() [1/2]

◆ a_associatedBodyIdsMb() [2/2]

◆ a_cellLengthAtLevel()

◆ a_childId()

◆ a_initTemperatureKelvin()

◆ a_isActive()

◆ a_isThermal()

◆ a_levelSetFunctionMb() [1/2]

◆ a_levelSetFunctionMb() [2/2]

◆ a_Ma()

◆ a_mbCell()

◆ a_noCells()

◆ a_noCollectorBodies()

◆ a_noDistributions()

◆ a_noEmbeddedBodiesLB()

◆ a_noLbCells()

◆ a_noLevelSetsMb()

◆ a_noVariables()

◆ a_oldVariable()

◆ a_parentId()

◆ a_pvrho()

◆ a_pvt()

◆ a_pvu()

◆ a_pvv()

◆ a_pvw()

◆ a_Re()

◆ a_variable()

◆ a_wasActive()

◆ averageTimer()

◆ bc3060()

◆ bodies()

◆ checkProperties()

◆ cleanUp()

◆ computeGCellTimeStep()

◆ constructGField()

◆ constructGField_()

◆ createBodyTree()

◆ descendLevelSetValue()

◆ finalizeAdaptation()

◆ finalizeCouplerInit()

◆ finalizeSubCoupleInit()

◆ getBodyAngularVelocity()

◆ getBodyVelocity()

◆ getDomainDecompositionInformation()

◆ init()

◆ initData()

◆ initializeSolidDomain()

◆ initTimers()

◆ lbBndCnd()

◆ lbSolver()

◆ minCell()

◆ noLevelSetFieldData()

◆ noMinCells()

◆ postAdaptation()

◆ postCouple()

◆ preCouple()

◆ readProperties()

◆ reinitAfterBalance()

◆ returnLevelSetSignForFluidCellCorners()

◆ returnLevelSetSignForFluidFaceCentroid()

◆ returnNoActiveCorners()

◆ subCouple()