maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn > Class Template Reference

LB sponge source term : sponge towards const target f_eq(rho_trg, u) More...

Inheritance diagram for maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >:

Collaboration diagram for maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >:

Public Types
using	Base = LbSrcTerm_sponge< nDim, nDist, SysEqn >
Public Types inherited from maia::lb::LbSrcTerm< nDim, nDist, SysEqn >
using	SysEqn = SysEqn

Public Member Functions
	LbSrcTerm_spongeRhoConst (LbSolverDxQy< nDim, nDist, SysEqn > *p_solver)
void	apply_preCollision () override
void	apply_postCollision () override
void	apply_postPropagation () override
Public Member Functions inherited from maia::lb::LbSrcTerm_sponge< nDim, nDist, SysEqn >
	LbSrcTerm_sponge (LbSolverDxQy< nDim, nDist, SysEqn > *p_solver)
void	init () override
	Initialize properties common by all lb sponge terms. More...
Public Member Functions inherited from maia::lb::LbSrcTerm< nDim, nDist, SysEqn >
virtual void	init ()=0
virtual void	apply_preCollision ()=0
virtual void	apply_postCollision ()=0
virtual void	apply_postPropagation ()=0
virtual	~LbSrcTerm ()=default

Additional Inherited Members
Static Public Attributes inherited from maia::lb::LbSrcTerm< nDim, nDist, SysEqn >
static constexpr MInt	nDim
static constexpr MInt	nDist
Protected Member Functions inherited from maia::lb::LbSrcTerm_sponge< nDim, nDist, SysEqn >
MInt	a_spongeCellId (const MInt cellId)
	Get index of cellId in m_spongeCells, returns -1 if unavailable. More...
void	readProperties () override
	Reading basics properties common for all lb sponge terms. More...
virtual void	readProperties ()=0
Protected Attributes inherited from maia::lb::LbSrcTerm_sponge< nDim, nDist, SysEqn >
LbSolverDxQy< nDim, nDist, SysEqn > *	m_solver
MBool	m_isActive = false
std::vector< SpongeCell >	m_spongeCells
std::unordered_map< MInt, MInt >	m_mapCellId2SpongeCellId
MFloat	m_spongeLayerThickness = 0.0
std::array< MFloat, 2 *nDim >	m_spongeThicknessFactor {}
MFloat	m_spongeSigma = 0.0
MFloat	m_spongeBeta = 2.0
MFloat	m_trgRho = 1.0
std::array< MFloat, nDim >	m_trgU {}

Detailed Description

template<MInt nDim, MInt nDist, class SysEqn>
class maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >

Author

Miro Gondrum

Date

04.01.2024 The PPDF f is damped towards an equilibrium state f_a following: f -= sigma *(f_eq - f_a) , with f_a = f_eq(rho_trg, u) Here, rho_trg is a user defined input and for this sponge are constant in space and time and u is the local velocity

Definition at line 25 of file lbsrctermsponge.cpp.

Member Typedef Documentation

Base

template<MInt nDim, MInt nDist, class SysEqn >

using maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >::Base = LbSrcTerm_sponge<nDim, nDist, SysEqn>

Definition at line 27 of file lbsrctermsponge.cpp.

Constructor & Destructor Documentation

LbSrcTerm_spongeRhoConst()

template<MInt nDim, MInt nDist, class SysEqn >

maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >::LbSrcTerm_spongeRhoConst ( LbSolverDxQy< nDim, nDist, SysEqn > *  p_solver )

inline

Definition at line 28 of file lbsrctermsponge.cpp.

: Base(p_solver){};

Member Function Documentation

apply_postCollision()

template<MInt nDim, MInt nDist, class SysEqn >

void maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >::apply_postCollision

overridevirtual

Implements maia::lb::LbSrcTerm< nDim, nDist, SysEqn >.

Definition at line 37 of file lbsrctermsponge.cpp.

                                                                        {
  TRACE();
  LbSolverDxQy<nDim, nDist, SysEqn>* const s = this->m_solver; // alias for readability
  using Ld = LbLatticeDescriptor<nDim, nDist>;
 
  const MInt timeStep = s->getCurrentTimeStep();
  const MUint noSpongeCells = this->m_spongeCells.size();
  maia::parallelFor<true>(0, noSpongeCells, [=](MInt spongeCellId) {
    auto& cell = this->m_spongeCells[spongeCellId];
    const MInt cellId = cell.cellId;
    if((timeStep - 1) % IPOW2(s->maxLevel() - s->a_level(cellId)) != 0) return;
    // compute the forcing term for the density
    const MFloat deltaRho = s->a_variable(cellId, s->PV->RHO) - this->m_trgRho;
    s->a_variable(cellId, s->PV->RHO) -= cell.spongeFactor * deltaRho;
 
    // prescribe new density
    s->a_distribution(cellId, Ld::lastId()) = s->a_variable(cellId, s->PV->RHO);
    for(MInt dist = 0; dist < nDist - 1; dist++) {
      s->a_distribution(cellId, nDist - 1) -= s->a_distribution(cellId, dist);
    }
  });
}

apply_postPropagation()

template<MInt nDim, MInt nDist, class SysEqn >

void maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >::apply_postPropagation ( )

inlineoverridevirtual

Implements maia::lb::LbSrcTerm< nDim, nDist, SysEqn >.

Definition at line 31 of file lbsrctermsponge.cpp.

{};

apply_preCollision()

template<MInt nDim, MInt nDist, class SysEqn >

void maia::lb::LbSrcTerm_spongeRhoConst< nDim, nDist, SysEqn >::apply_preCollision ( )

inlineoverridevirtual

Implements maia::lb::LbSrcTerm< nDim, nDist, SysEqn >.

Definition at line 29 of file lbsrctermsponge.cpp.

{};

---

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

• /home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/LB/lbsrctermsponge.cpp