rotate_flux Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for rotate_flux:

Public Member Functions
	rotate_flux (int l_zlo)
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_q_flux (const int &i, const int &j, const int &k, const int &n, const int &icomp, const amrex::Real &dz, const amrex::Real &dz1, const bool &, const amrex::Array4< const amrex::Real > &, 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 > &, 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< amrex::Real > &dest_arr) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_t_flux (const int &i, const int &j, const int &k, const int &n, const int &icomp, const amrex::Real &dz, const amrex::Real &dz1, const bool &, const amrex::Array4< const amrex::Real > &, 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 > &, const amrex::Array4< const amrex::Real > &tm_arr, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &t_star_arr, const amrex::Array4< const amrex::Real > &t_surf_arr, const amrex::Array4< amrex::Real > &dest_arr) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_u_flux (const int &i, const int &j, const int &k, const int &icomp, const amrex::Real &dz, const amrex::Real &dz1, const bool &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &cons_arr, const amrex::Array4< const amrex::Real > &velx_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 > &u_star_arr, const amrex::Array4< amrex::Real > &dest_arr) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	compute_v_flux (const int &i, const int &j, const int &k, const int &icomp, const amrex::Real &dz, const amrex::Real &dz1, const bool &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, 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< amrex::Real > &dest_arr) const

Private Attributes
int	zlo
const amrex::Real	eps = 1e-15

Detailed Description

Rotate flux formulation

Constructor & Destructor Documentation

rotate_flux()

rotate_flux::rotate_flux ( int  l_zlo )

inline

      :  zlo(l_zlo) {}

Member Function Documentation

compute_q_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real rotate_flux::compute_q_flux ( const int &  i, const int &  j, const int &  k, const int &  n, const int &  icomp, const amrex::Real &  dz, const amrex::Real &  dz1, const bool &  , const amrex::Array4< const amrex::Real > &  , 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 > &  , 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< amrex::Real > &  dest_arr  ) const

inline

    {
        int ic, jc;
        ic = i  < lbound(cons_arr).x ? lbound(cons_arr).x : i;
        jc = j  < lbound(cons_arr).y ? lbound(cons_arr).y : j;
        ic = ic > ubound(cons_arr).x ? ubound(cons_arr).x : ic;
        jc = jc > ubound(cons_arr).y ? ubound(cons_arr).y : jc;
 
        // surface gradient equal to gradient at first zface
        amrex::Real deltaz      = dz * (zlo - k);
        amrex::Real rqvgrad  = ( cons_arr(ic,jc,zlo+1,RhoQ1_comp)
                               - cons_arr(ic,jc,zlo  ,RhoQ1_comp) ) / (0.5*(dz+dz1));
        dest_arr(i,j,k,icomp+n) = cons_arr(ic,jc,zlo,RhoQ1_comp) - rqvgrad * deltaz;
 
        // NOTE: this is the total stress
        amrex::Real rho    = cons_arr(ic,jc,zlo,Rho_comp);
        amrex::Real qstar  = 0.0; //q_star_arr(ic,jc,zlo);
        amrex::Real moflux = (std::abs(qstar) > eps) ? -rho * qstar : 0.0;
 
        return moflux;
    }

compute_t_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real rotate_flux::compute_t_flux ( const int &  i, const int &  j, const int &  k, const int &  n, const int &  icomp, const amrex::Real &  dz, const amrex::Real &  dz1, const bool &  , const amrex::Array4< const amrex::Real > &  , 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 > &  , const amrex::Array4< const amrex::Real > &  tm_arr, const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  t_star_arr, const amrex::Array4< const amrex::Real > &  t_surf_arr, const amrex::Array4< amrex::Real > &  dest_arr  ) const

inline

    {
        int ic, jc;
        ic = i  < lbound(cons_arr).x ? lbound(cons_arr).x : i;
        jc = j  < lbound(cons_arr).y ? lbound(cons_arr).y : j;
        ic = ic > ubound(cons_arr).x ? ubound(cons_arr).x : ic;
        jc = jc > ubound(cons_arr).y ? ubound(cons_arr).y : jc;
 
        // surface gradient equal to gradient at first zface
        amrex::Real deltaz      = dz * (zlo - k);
        amrex::Real rthetagrad  = ( cons_arr(ic,jc,zlo+1,RhoTheta_comp)
                                  - cons_arr(ic,jc,zlo  ,RhoTheta_comp) ) / (0.5*(dz+dz1));
        dest_arr(i,j,k,icomp+n) = cons_arr(ic,jc,zlo,RhoTheta_comp) - rthetagrad * deltaz;
 
        // NOTE: this is the total stress
        amrex::Real theta_mean = tm_arr(ic,jc,zlo);
        amrex::Real theta_surf = t_surf_arr(ic,jc,zlo);
        amrex::Real rho    = cons_arr(ic,jc,zlo,Rho_comp);
        amrex::Real tstar  = t_star_arr(ic,jc,zlo);
        amrex::Real moflux = (std::abs(tstar) > eps) ? -rho * tstar * (theta_mean-theta_surf) : 0.0;
 
        return moflux;
    }

compute_u_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real rotate_flux::compute_u_flux ( const int &  i, const int &  j, const int &  k, const int &  icomp, const amrex::Real &  dz, const amrex::Real &  dz1, const bool &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  cons_arr, const amrex::Array4< const amrex::Real > &  velx_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 > &  u_star_arr, const amrex::Array4< amrex::Real > &  dest_arr  ) const

inline

    {
        int ic, jc;
        ic = i  < lbound(cons_arr).x+1 ? lbound(cons_arr).x+1 : i;
        jc = j  < lbound(cons_arr).y   ? lbound(cons_arr).y   : j;
        ic = ic > ubound(cons_arr).x   ? ubound(cons_arr).x   : ic;
        jc = jc > ubound(cons_arr).y   ? ubound(cons_arr).y   : jc;
 
        // surface gradient equal to gradient at first zface
        amrex::Real deltaz    = dz * (zlo - k);
        amrex::Real ugrad     = (velx_arr(i,j,zlo+1) - velx_arr(i,j,zlo)) / (0.5*(dz+dz1));
        dest_arr(i,j,k,icomp) = velx_arr(i,j,zlo) - ugrad * deltaz;
 
        // NOTE: this is the total stress
        amrex::Real rho     = 0.5 *( cons_arr(ic-1,jc,zlo,Rho_comp)
                                   + cons_arr(ic  ,jc,zlo,Rho_comp) );
        amrex::Real ustar   = 0.5 * ( u_star_arr(ic-1,jc,zlo) + u_star_arr(ic,jc,zlo) );
        amrex::Real stressx = -rho * ustar * ustar;
 
        return stressx;
    }

compute_v_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real rotate_flux::compute_v_flux ( const int &  i, const int &  j, const int &  k, const int &  icomp, const amrex::Real &  dz, const amrex::Real &  dz1, const bool &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , 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< amrex::Real > &  dest_arr  ) const

inline

    {
        // surface gradient equal to gradient at first zface
        amrex::Real deltaz    = dz * (zlo - k);
        amrex::Real vgrad     = (vely_arr(i,j,zlo+1) - vely_arr(i,j,zlo)) / (0.5*(dz+dz1));
        dest_arr(i,j,k,icomp) = vely_arr(i,j,zlo) - vgrad * deltaz;
 
        // NOTE: this is the total stress
        amrex::Real stressy = 0.0;
 
        return stressy;
    }

Member Data Documentation

eps

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

private

Referenced by compute_q_flux(), and compute_t_flux().

zlo

int rotate_flux::zlo

private

Referenced by 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