moeng_flux Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for moeng_flux:

Public Member Functions
	moeng_flux ()
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_q_flux (const int &, const int &, const int &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_t_flux (const int &i, const int &j, const int &k, const amrex::Array4< const amrex::Real > &cons_arr, const amrex::Array4< const amrex::Real > &velx_arr, const amrex::Array4< const amrex::Real > &vely_arr, const amrex::Array4< const amrex::Real > &umm_arr, const amrex::Array4< const amrex::Real > &tm_arr, const amrex::Array4< const amrex::Real > &u_star_arr, const amrex::Array4< const amrex::Real > &t_star_arr, const amrex::Array4< const amrex::Real > &t_surf_arr) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_u_flux (const int &i, const int &j, const int &k, const amrex::Array4< const amrex::Real > &cons_arr, const amrex::Array4< const amrex::Real > &velx_arr, const amrex::Array4< const amrex::Real > &vely_arr, const amrex::Array4< const amrex::Real > &umm_arr, const amrex::Array4< const amrex::Real > &um_arr, const amrex::Array4< const amrex::Real > &u_star_arr) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_v_flux (const int &i, const int &j, const int &k, const amrex::Array4< const amrex::Real > &cons_arr, const amrex::Array4< const amrex::Real > &velx_arr, const amrex::Array4< const amrex::Real > &vely_arr, const amrex::Array4< const amrex::Real > &umm_arr, const amrex::Array4< const amrex::Real > &vm_arr, const amrex::Array4< const amrex::Real > &u_star_arr) const

Private Attributes
const amrex::Real	eps = 1e-15
const amrex::Real	WSMIN = 0.1

Detailed Description

Moeng flux formulation

Constructor & Destructor Documentation

moeng_flux()

moeng_flux::moeng_flux ( )

inline

{}

Member Function Documentation

compute_q_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux::compute_q_flux ( const int &  , const int &  , const int &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &    ) const

inline

    {
        // NOTE: this is rho*<q'w'> = -K dqdz
        amrex::Real moflux  = 0.0;
 
        return moflux;
    }

compute_t_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux::compute_t_flux ( const int &  i, const int &  j, const int &  k, const amrex::Array4< const amrex::Real > &  cons_arr, const amrex::Array4< const amrex::Real > &  velx_arr, const amrex::Array4< const amrex::Real > &  vely_arr, const amrex::Array4< const amrex::Real > &  umm_arr, const amrex::Array4< const amrex::Real > &  tm_arr, const amrex::Array4< const amrex::Real > &  u_star_arr, const amrex::Array4< const amrex::Real > &  t_star_arr, const amrex::Array4< const amrex::Real > &  t_surf_arr  ) const

inline

    {
        amrex::Real rho   = cons_arr(i,j,k,Rho_comp);
        amrex::Real theta = cons_arr(i,j,k,RhoTheta_comp) / rho;
        amrex::Real velx  = 0.5 * ( velx_arr(i,j,k) + velx_arr(i+1,j  ,k) );
        amrex::Real vely  = 0.5 * ( vely_arr(i,j,k) + vely_arr(i  ,j+1,k) );
 
        amrex::Real theta_mean  = tm_arr(i,j,k);
        amrex::Real ustar       = u_star_arr(i,j,k);
        amrex::Real tstar       = t_star_arr(i,j,k);
        amrex::Real theta_surf  = t_surf_arr(i,j,k);
        amrex::Real wsp_mean    = umm_arr(i,j,k);
        wsp_mean = std::max(wsp_mean, WSMIN);
 
        amrex::Real wsp     = sqrt(velx*velx+vely*vely);
        amrex::Real num1    = wsp * (theta_mean-theta_surf);
        amrex::Real num2    = wsp_mean * (theta-theta_mean);
 
        // NOTE: this is rho*<T'w'> = -K dTdz
        amrex::Real moflux  = (std::abs(tstar) > eps) ?
                              -rho*tstar*ustar*(num1+num2)/((theta_mean-theta_surf)*wsp_mean) : 0.0;
 
        return moflux;
    }

compute_u_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux::compute_u_flux ( const int &  i, const int &  j, const int &  k, const amrex::Array4< const amrex::Real > &  cons_arr, const amrex::Array4< const amrex::Real > &  velx_arr, const amrex::Array4< const amrex::Real > &  vely_arr, const amrex::Array4< const amrex::Real > &  umm_arr, const amrex::Array4< const amrex::Real > &  um_arr, const amrex::Array4< const amrex::Real > &  u_star_arr  ) const

inline

    {
        amrex::Real velx  = velx_arr(i,j,k);
        amrex::Real vely  = 0.25 * ( vely_arr(i  ,j,k) + vely_arr(i  ,j+1,k)
                                   + vely_arr(i-1,j,k) + vely_arr(i-1,j+1,k) );
        amrex::Real rho   = 0.5  * ( cons_arr(i-1,j,k,Rho_comp) + cons_arr(i,j,k,Rho_comp) );
 
        amrex::Real umean    = um_arr(i,j,k);
        amrex::Real ustar    = 0.5 * ( u_star_arr(i-1,j,k) + u_star_arr(i,j,k) );
        amrex::Real wsp_mean = 0.5 * (    umm_arr(i-1,j,k) +    umm_arr(i,j,k) );
        wsp_mean = std::max(wsp_mean, WSMIN);
 
        // Note: The surface mean shear stress is decomposed into tau_xz by
        //       multiplying the modeled shear stress (rho*ustar^2) with
        //       a factor of umean/wsp_mean for directionality; this factor
        //       modifies the denominator from what is in Moeng 1984.
        amrex::Real wsp     = sqrt(velx*velx+vely*vely);
        amrex::Real num1    = wsp * umean;
        amrex::Real num2    = wsp_mean * (velx-umean);
 
        // NOTE: this is rho*<u'w'> = -K dudz
        amrex::Real stressx = -rho*ustar*ustar * (num1+num2)/(wsp_mean*wsp_mean);
 
        return stressx;
    }

compute_v_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux::compute_v_flux ( const int &  i, const int &  j, const int &  k, const amrex::Array4< const amrex::Real > &  cons_arr, const amrex::Array4< const amrex::Real > &  velx_arr, const amrex::Array4< const amrex::Real > &  vely_arr, const amrex::Array4< const amrex::Real > &  umm_arr, const amrex::Array4< const amrex::Real > &  vm_arr, const amrex::Array4< const amrex::Real > &  u_star_arr  ) const

inline

    {
        amrex::Real velx  = 0.25 * ( velx_arr(i,j  ,k) + velx_arr(i+1,j  ,k)
                                   + velx_arr(i,j-1,k) + velx_arr(i+1,j-1,k) );
        amrex::Real vely  = vely_arr(i,j,k);
        amrex::Real rho   = 0.5  * ( cons_arr(i,j-1,k,Rho_comp) + cons_arr(i,j,k,Rho_comp) );
 
        amrex::Real vmean     = vm_arr(i,j,k);
        amrex::Real ustar     = 0.5 * ( u_star_arr(i,j-1,k) + u_star_arr(i,j,k) );
        amrex::Real wsp_mean  = 0.5 * (    umm_arr(i,j-1,k) +    umm_arr(i,j,k) );
        wsp_mean = std::max(wsp_mean, WSMIN);
 
        // Note: The surface mean shear stress is decomposed into tau_yz by
        //       multiplying the modeled shear stress (rho*ustar^2) with
        //       a factor of vmean/wsp_mean for directionality; this factor
        //       modifies the denominator from what is in Moeng 1984.
        amrex::Real wsp     = sqrt(velx*velx+vely*vely);
        amrex::Real num1    = wsp * vmean;
        amrex::Real num2    = wsp_mean * (vely-vmean);
 
        // NOTE: this is rho*<v'w'> = -K dvdz
        amrex::Real stressy = -rho*ustar*ustar * (num1+num2)/(wsp_mean*wsp_mean);
 
        return stressy;
    }

Member Data Documentation

eps

const amrex::Real moeng_flux::eps = 1e-15

private

Referenced by compute_t_flux().

WSMIN

const amrex::Real moeng_flux::WSMIN = 0.1

private

Referenced by compute_t_flux(), compute_u_flux(), and compute_v_flux().

---

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

• Source/BoundaryConditions/ERF_MOSTStress.H