bulk_coeff_flux Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for bulk_coeff_flux:

Public Member Functions
	bulk_coeff_flux (amrex::Real m_Cd, amrex::Real m_Ch, amrex::Real m_Cq)
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_q_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 > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &q_surf_arr) 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 > &, const amrex::Array4< const amrex::Real > &, 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 > &, 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 > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &) const

Private Attributes
most_data	mdata

Detailed Description

Bulk coefficient flux formulation

Constructor & Destructor Documentation

bulk_coeff_flux()

bulk_coeff_flux::bulk_coeff_flux ( amrex::Real  m_Cd, amrex::Real  m_Ch, amrex::Real  m_Cq  )

inline

Construct the bulk-coefficient flux calculator.

Parameters

[in]	m_Cd	drag coefficient
[in]	m_Ch	heat-transfer coefficient
[in]	m_Cq	moisture-transfer coefficient

    {
        mdata.Cd = m_Cd;
        mdata.Ch = m_Ch;
        mdata.Cq = m_Cq;
    }

Member Function Documentation

compute_q_flux()

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

inline

Compute moisture flux for one surface point.

Calling sequence

i, j, k, cons_arr, velx_arr, vely_arr, umm_arr, qvm_arr, u_star_arr, q_star_arr, q_surf_arr.

    {
        amrex::Real rho    = cons_arr(i,j,k,Rho_comp);
        amrex::Real qv     = cons_arr(i,j,k,RhoQ1_comp)/cons_arr(i,j,k,Rho_comp);
        amrex::Real qvsurf = q_surf_arr(i,j,0);
        amrex::Real velx   = myhalf * ( velx_arr(i,j,k) + velx_arr(i+1,j  ,k) );
        amrex::Real vely   = myhalf * ( vely_arr(i,j,k) + vely_arr(i  ,j+1,k) );
        amrex::Real wsp    = std::sqrt(velx*velx+vely*vely);
 
        // NOTE: this is rho*<q'w'> = -K dqdz
        amrex::Real moflux  = rho * mdata.Cq * wsp * (qvsurf - qv);
 
        return moflux;
    }

compute_t_flux()

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

inline

Compute heat flux for one surface point.

Calling sequence

i, j, k, cons_arr, velx_arr, vely_arr, umm_arr, tm_arr, u_star_arr, t_star_arr, t_surf_arr.

    {
        amrex::Real rho    = cons_arr(i,j,k,Rho_comp);
        amrex::Real th     = cons_arr(i,j,k,RhoTheta_comp)/cons_arr(i,j,k,Rho_comp);
        amrex::Real thsurf = t_surf_arr(i,j,0);
        amrex::Real velx   = myhalf * ( velx_arr(i,j,k) + velx_arr(i+1,j  ,k) );
        amrex::Real vely   = myhalf * ( vely_arr(i,j,k) + vely_arr(i  ,j+1,k) );
        amrex::Real wsp    = std::sqrt(velx*velx+vely*vely);
 
        // NOTE: this is rho*<T'w'> = -K dTdz
        amrex::Real moflux  = rho * mdata.Ch * wsp * (thsurf - th);
 
        return moflux;
    }

compute_u_flux()

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

inline

Compute x-momentum flux for one surface face.

Calling sequence

i, j, k, cons_arr, velx_arr, vely_arr, umm_arr, um_arr, u_star_arr.

    {
        amrex::Real rho   = myhalf * ( cons_arr(i-1,j,k,Rho_comp) + cons_arr(i,j,k,Rho_comp) );
        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 wsp   = std::sqrt(velx*velx+vely*vely);
 
        // NOTE: this is rho*<u'w'> = -K dudz
        // NOTE: tau_tot = rho * Cd * wsp^2, multiply by u/wsp to get tau_13
        amrex::Real stressx = -rho * mdata.Cd * wsp * velx;
 
        return stressx;
    }

compute_v_flux()

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

inline

Compute y-momentum flux for one surface face.

Calling sequence

i, j, k, cons_arr, velx_arr, vely_arr, umm_arr, vm_arr, u_star_arr.

    {
        amrex::Real rho   = myhalf * ( cons_arr(i,j-1,k,Rho_comp) + cons_arr(i,j,k,Rho_comp) );
        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 wsp   = std::sqrt(velx*velx+vely*vely);
 
        // NOTE: this is rho*<v'w'> = -K dvdz
        // NOTE: tau_tot = rho * Cd * wsp^2, multiply by v/wsp to get tau_13
        amrex::Real stressy = -rho * mdata.Cd * wsp * vely;
 
        return stressy;
    }

Member Data Documentation

mdata

most_data bulk_coeff_flux::mdata

private

Referenced by bulk_coeff_flux(), compute_q_flux(), 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