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

Class to handle sub-grid scale model (Smagorinsky) More...

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

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

Public Member Functions
	LbSrcTerm_smagorinsky (LbSolverDxQy< nDim, nDist, SysEqn > *p_solver)
void	init () override
void	apply_preCollision () override
void	apply_postCollision () override
void	apply_postPropagation () override
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

Protected Member Functions
void	readProperties () override
virtual void	readProperties ()=0

Protected Attributes
LbSolverDxQy< nDim, nDist, SysEqn > *	m_solver
MFloat	m_Cs = 0.1
	Smagorinsky constant. More...
const MFloat	m_deltaX = 1.0
	Filter width. More...

Additional Inherited Members
Public Types inherited from maia::lb::LbSrcTerm< nDim, nDist, SysEqn >
using	SysEqn = SysEqn
Static Public Attributes inherited from maia::lb::LbSrcTerm< nDim, nDist, SysEqn >
static constexpr MInt	nDim
static constexpr MInt	nDist

Detailed Description

template<MInt nDim, MInt nDist, class SysEqn, MBool bubble = false>
class maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >

Author

Miro Gondrum

Date

01.02.2022

Definition at line 23 of file lbsrctermsgs.cpp.

Constructor & Destructor Documentation

LbSrcTerm_smagorinsky()

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble = false>

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

inline

Definition at line 34 of file lbsrctermsgs.cpp.

: m_solver(p_solver) { readProperties(); };

Member Function Documentation

apply_postCollision()

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble = false>

void maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >::apply_postCollision ( )

inlineoverridevirtual

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

Definition at line 38 of file lbsrctermsgs.cpp.

{/*NOT USED*/};

apply_postPropagation()

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble = false>

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

inlineoverridevirtual

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

Definition at line 39 of file lbsrctermsgs.cpp.

{/*NOT USED*/};

apply_preCollision()

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble>

void maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >::apply_preCollision

overridevirtual

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

Definition at line 53 of file lbsrctermsgs.cpp.

                                                                            {
  TRACE();
  LbSolverDxQy<nDim, nDist, SysEqn>* const s = m_solver; // alias for readability
 
  s->m_nu = s->m_Ma * LBCS / s->m_Re * s->m_referenceLength;
  maia::parallelFor<false>(0, s->m_currentMaxNoCells, [=](MInt index) {
    const MInt pCellId = s->m_activeCellList[index];
    const MInt lvlDiff = s->maxLevel() - s->a_level(pCellId);
    if((globalTimeStep - 1) % IPOW2(lvlDiff) != 0) return;
    // TODO labels:LB daniell: This is critical for setting a_oldNu ?! But without
    // test case 3D_LB-FV_Euler-Euler_bubble_sphere fails
    s->a_nu(pCellId) = s->m_nu;
    // Ref: Hou, Sterling, Chen, and Doolen 1994
    // Also worth reading: Siebler, Krafczyk, Freudiger, and Geier 2011
    constexpr MInt nDimSqr = nDim * nDim;
    MFloat c[nDimSqr]{};
    const MFloat rho = s->a_variable(pCellId, s->PV->RHO);
#ifndef WAR_NVHPC_PSTL
    const MFloat* const u = &s->a_variable(pCellId, s->PV->U);
    const MFloat* const dist = &s->a_oldDistribution(pCellId, 0);
    s->sysEqn().calcMomentumFlux(rho, u, dist, c);
#else
    MFloat u[nDim] = {F0};
    MFloat dist[nDist] = {F0};
    for(MInt d = 0; d < nDim; d++) {
      u[d] = s->a_variable(pCellId, d);
    }
    for(MInt d = 0; d < nDist; d++) {
      dist[d] = s->a_oldDistribution(pCellId, d);
    }
    s->sysEqn().calcMomentumFlux(rho, u, dist, c);
#endif
    const MFloat Q = sqrt(2.0 * std::inner_product(&c[0], &c[nDimSqr], &c[0], .0));
    MFloat nu_BIT = 0.0;
    if constexpr(bubble) {
      if(s->m_isEELiquid) {
        // BIT (bubble induced turbulence) see Mohammadi 19 / Sato 81 / Milelli 01
        const MFloat u_d = sqrt(POW2(s->a_uOtherPhase(pCellId, 0) - u[0]) + POW2(s->a_uOtherPhase(pCellId, 1) - u[1])
                                + POW2(s->a_uOtherPhase(pCellId, 2) - u[2]));
        nu_BIT = m_Cs * s->a_alphaGasLim(pCellId) * u_d;
      }
    }
    // The stress tensor S_mean depends on tau and tau depends on S_mean.
    // This leads to a quadratic equation, that is solved by this.
    MFloat tauBase = F1B2 + 3.0 * (s->m_nu + nu_BIT) * FFPOW2(s->maxLevel() - s->a_level(pCellId));
    const MFloat tau = F1B2
                       * (tauBase
                          + sqrt(tauBase * tauBase
                                 + 2.0 * m_Cs * m_Cs * m_deltaX * m_deltaX * (F1BCSsq * F1BCSsq) * Q
                                       * FPOW2(s->maxLevel() - s->a_level(pCellId))));
 
    // eddy viscosity
    if constexpr(bubble) {
      if(s->m_cells.saveOldNu()) {
        s->a_oldNu(pCellId) = s->a_nu(pCellId);
      }
    }
    s->a_nu(pCellId) = (2.0 * tau - 1.0) * CSsq * FPOW2(s->maxLevel() - s->a_level(pCellId)) / 2.0;
    if constexpr(bubble) {
      if(s->m_isEELiquid) {
        if(s->m_cells.saveOldNu()) {
          s->a_oldNuT(pCellId) = s->a_nuT(pCellId);
        }
        s->a_nuT(pCellId) = s->a_nu(pCellId) - s->m_nu - nu_BIT;
      }
    }
  });
}

init()

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble = false>

void maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >::init ( )

inlineoverridevirtual

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

Definition at line 36 of file lbsrctermsgs.cpp.

{/*NOT USED*/};

readProperties()

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble>

void maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >::readProperties

overrideprotectedvirtual

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

Definition at line 44 of file lbsrctermsgs.cpp.

                                                                        {
  TRACE();
  const MInt solverId = m_solver->m_solverId;
  if(Context::propertyExists("smagorinskyConstant")) {
    m_Cs = Context::getSolverProperty<MFloat>("smagorinskyConstant", solverId, AT_);
  }
}

Member Data Documentation

m_Cs

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble = false>

MFloat maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >::m_Cs = 0.1

protected

Definition at line 26 of file lbsrctermsgs.cpp.

m_deltaX

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble = false>

const MFloat maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >::m_deltaX = 1.0

protected

Definition at line 27 of file lbsrctermsgs.cpp.

m_solver

template<MInt nDim, MInt nDist, class SysEqn , MBool bubble = false>

LbSolverDxQy<nDim, nDist, SysEqn>* maia::lb::LbSrcTerm_smagorinsky< nDim, nDist, SysEqn, bubble >::m_solver

protected

Definition at line 25 of file lbsrctermsgs.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/lbsrctermsgs.cpp