donelan_flux Struct Reference

#include <ERF_MOSTStress.H>

Public Member Functions
	donelan_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 > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &umm_arr, const amrex::Array4< const amrex::Real > &tm_arr, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, 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 > &, const amrex::Array4< const amrex::Real > &) 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 > &, const amrex::Array4< const amrex::Real > &) const

Detailed Description

Donelan flux formulation

Constructor & Destructor Documentation

donelan_flux()

donelan_flux::donelan_flux ( )

inline

{}

Member Function Documentation

compute_q_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real donelan_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  = zero;
 
        return moflux;
    }

compute_t_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real donelan_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 > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  umm_arr, const amrex::Array4< const amrex::Real > &  tm_arr, const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  t_surf_arr  ) const

inline

    {
        amrex::Real Ch          = amrex::Real(0.0012);
        amrex::Real rho         = cons_arr(i,j,k,Rho_comp);
 
        amrex::Real theta_surf  = t_surf_arr(i,j,0);
        amrex::Real theta_mean  =     tm_arr(i,j,0);
        amrex::Real wsp_mean    =    umm_arr(i,j,0);
 
        // NOTE: this is rho*<T'w'> = -K dTdz
        amrex::Real moflux = -rho * Ch * wsp_mean * (theta_mean - theta_surf);
 
        return moflux;
    }

compute_u_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real donelan_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 > &  , const amrex::Array4< const amrex::Real > &    ) const

inline

    {
        amrex::Real velx  = velx_arr(i,j,k);
        amrex::Real vely  = fourth * ( 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   = myhalf  * ( cons_arr(i-1,j,k,Rho_comp) + cons_arr(i,j,k,Rho_comp) );
 
        amrex::Real Cd       = amrex::Real(0.001);
        const amrex::Real c  = amrex::Real(7e-5);
        amrex::Real wsp      = std::sqrt(velx*velx+vely*vely);
        amrex::Real wsp_mean = myhalf * ( umm_arr(i-1,j,0) + umm_arr(i,j,0) );
        if (wsp_mean <= amrex::Real(5.0)) {
            Cd = amrex::Real(0.001);
        } else if (wsp_mean < amrex::Real(25.0) && wsp_mean > amrex::Real(5.0)) {
            Cd = amrex::Real(0.001) + c * (wsp_mean - amrex::Real(5.0));
        } else {
            Cd = amrex::Real(0.0024);
        }
 
        // NOTE: this is rho*<u'w'> = -K dudz
        amrex::Real stressx = -rho * Cd * velx * wsp;
 
        return stressx;
    }

compute_v_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real donelan_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 > &  , const amrex::Array4< const amrex::Real > &    ) const

inline

    {
        amrex::Real velx  = fourth * ( 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   = myhalf  * ( cons_arr(i,j-1,k,Rho_comp) + cons_arr(i,j,k,Rho_comp) );
 
        amrex::Real Cd       = amrex::Real(0.001);
        const amrex::Real c  = amrex::Real(7e-5);
        amrex::Real wsp      = std::sqrt(velx*velx+vely*vely);
        amrex::Real wsp_mean = myhalf * ( umm_arr(i,j-1,0) + umm_arr(i,j,0) );
        if (wsp_mean <= amrex::Real(5.0)) {
            Cd = amrex::Real(0.001);
        } else if (wsp_mean < amrex::Real(25.0) && wsp_mean > amrex::Real(5.0)) {
            Cd = amrex::Real(0.001) + c * (wsp_mean - amrex::Real(5.0));
        } else {
            Cd = amrex::Real(0.0024);
        }
 
        // NOTE: this is rho*<v'w'> = -K dvdz
        amrex::Real stressy = -rho * Cd * vely * wsp;
 
        return stressy;
    }

---

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

• Source/BoundaryConditions/ERF_MOSTStress.H