Coefficients and stability functions for the Level 2 MYNN closure. More...

#include <ERF_MYNNStruct.H>

Collaboration diagram for MYNNLevel2:

Public Member Functions
void	init_coeffs (amrex::Real A1_lvl25, amrex::Real A2_lvl25, amrex::Real B1, amrex::Real B2, amrex::Real C1, amrex::Real C2, amrex::Real C3, amrex::Real, amrex::Real C5)
	Initialize Level 2 closure coefficients from Level 2.5 coefficients. More...

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_Rf (const amrex::Real GM, const amrex::Real GH) const
	Calculate the reduced Richardson number. More...

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_SH (const amrex::Real Rf) const
	Calculate the heat stability function. More...

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_SM (const amrex::Real Rf) const
	Calculate the momentum stability function. More...

Public Attributes
bool	initialized {false}

amrex::Real	A1

amrex::Real	A2

amrex::Real	F1

amrex::Real	F2

const amrex::Real	gam1 = amrex::Real(0.235)

amrex::Real	gam2

amrex::Real	Rf1

amrex::Real	Rf2

amrex::Real	Rfc

amrex::Real	Ri1

amrex::Real	Ri2

amrex::Real	Ri3

const amrex::Real	eps = std::numeric_limits<amrex::Real>::epsilon()

Detailed Description

Coefficients and stability functions for the Level 2 MYNN closure.

Member Function Documentation

◆ calc_Rf()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real MYNNLevel2::calc_Rf	(	const amrex::Real	GM,
		const amrex::Real	GH
	)		const

inline

Calculate the reduced Richardson number.

Parameters

[in]	GM	Momentum gradient.
[in]	GH	Heat gradient.

Returns: Reduced Richardson number.

     {
         AMREX_ASSERT(initialized);
         amrex::Real lGM = std::copysign(std::max(std::fabs(GM),eps),GM);
         amrex::Real Ri  = -GH/lGM;
         return Ri1*(Ri + Ri2 - std::sqrt(Ri*Ri - Ri3*Ri + Ri2*Ri2));
     }

◆ calc_SH()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real MYNNLevel2::calc_SH ( const amrex::Real Rf ) const

inline

Calculate the heat stability function.

Parameters

[in] Rf Reduced Richardson number.

Returns: Heat stability function.

     {
         AMREX_ASSERT(initialized);
         return three*A2*(gam1+gam2)*(Rfc-Rf)/(one-Rf);
     }

Referenced by calc_SM().

Here is the caller graph for this function:

◆ calc_SM()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real MYNNLevel2::calc_SM ( const amrex::Real Rf ) const

inline

Calculate the momentum stability function.

Parameters

[in] Rf Reduced Richardson number.

Returns: Momentum stability function.

     {
         AMREX_ASSERT(initialized);
         return A1*F1/(A2*F2) * (Rf1-Rf) / (Rf2-Rf) * calc_SH(Rf);
     }

Here is the call graph for this function:

◆ init_coeffs()

void MYNNLevel2::init_coeffs	(	amrex::Real	A1_lvl25,
		amrex::Real	A2_lvl25,
		amrex::Real	B1,
		amrex::Real	B2,
		amrex::Real	C1,
		amrex::Real	C2,
		amrex::Real	C3,
		amrex::Real	,
		amrex::Real	C5
	)

inline

Initialize Level 2 closure coefficients from Level 2.5 coefficients.

Parameters

[in]	A1_lvl25	Level 2.5 coefficient A1.
[in]	A2_lvl25	Level 2.5 coefficient A2.
[in]	B1	Closure coefficient B1.
[in]	B2	Closure coefficient B2.
[in]	C1	Closure coefficient C1.
[in]	C2	Closure coefficient C2.
[in]	C3	Closure coefficient C3.
[in]	C5	Closure coefficient C5.

     {
         A1 = A1_lvl25;
         A2 = A2_lvl25;
         gam2 = (two*A1*(three-two*C2) + B2*(one-C3)) / B1;
         F1  = B1*(gam1-C1) + two*A1*(three-two*C2) + three*A2*(one-C2)*(one-C5);
         F2  = B1*(gam1+gam2) - three*A1*(one-C2);
         Rf1 = B1*(gam1-C1) / F1;
         Rf2 = B1*gam1 / F2;
         Rfc = gam1 / (gam1 + gam2);
         Ri1 = myhalf*A2*F2/(A1*F1);
         Ri2 = myhalf*Rf1/Ri1;
         Ri3 = (two*Rf2-Rf1)/Ri1;
  
         initialized = true;
     }