ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Public Member Functions | Private Attributes | List of all members
moeng_flux_eb Struct Reference

#include <ERF_EBMOSTStress.H>

Collaboration diagram for moeng_flux_eb:
Collaboration graph
[legend]

Public Member Functions

 moeng_flux_eb ()
 
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 > &umm_arr, const amrex::Array4< const amrex::Real > &qvm_arr, const amrex::Array4< const amrex::Real > &u_star_arr, const amrex::Array4< const amrex::Real > &q_star_arr, const amrex::Array4< const amrex::Real > &q_surf_arr, const amrex::Array4< const amrex::Real > &u_vfrac_arr, const amrex::Array4< const amrex::Real > &v_vfrac_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 > &velz_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::Array4< const amrex::Real > &u_vfrac_arr, const amrex::Array4< const amrex::Real > &v_vfrac_arr, const amrex::Array4< const amrex::Real > &w_vfrac_arr, const amrex::Array4< const amrex::Real > &bnorm_arr) const
 
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real compute_u_flux (int i, int j, 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 > &velz_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::Array4< const amrex::Real > &u_vfrac_arr, const amrex::Array4< const amrex::Real > &v_vfrac_arr, const amrex::Array4< const amrex::Real > &w_vfrac_arr, const amrex::Array4< const amrex::Real > &cc_vfrac_arr, const amrex::Array4< const amrex::EBCellFlag > &cc_flag_arr, const amrex::Array4< const amrex::Real > &bnorm_arr, int idir=0) const
 
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real compute_v_flux (int i, int j, 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 > &velz_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 amrex::Array4< const amrex::Real > &u_vfrac_arr, const amrex::Array4< const amrex::Real > &v_vfrac_arr, const amrex::Array4< const amrex::Real > &w_vfrac_arr, const amrex::Array4< const amrex::Real > &cc_vfrac_arr, const amrex::Array4< const amrex::EBCellFlag > &cc_flag_arr, const amrex::Array4< const amrex::Real > &bnorm_arr, int idir=0) const
 

Private Attributes

const amrex::Real eps = amrex::Real(1e-12)
 
const amrex::Real WSMIN = amrex::Real(0.1)
 

Detailed Description

Moeng flux formulation

Constructor & Destructor Documentation

◆ moeng_flux_eb()

moeng_flux_eb::moeng_flux_eb ( )
inline
301 {}

Member Function Documentation

◆ compute_q_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux_eb::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 > &  umm_arr,
const amrex::Array4< const amrex::Real > &  qvm_arr,
const amrex::Array4< const amrex::Real > &  u_star_arr,
const amrex::Array4< const amrex::Real > &  q_star_arr,
const amrex::Array4< const amrex::Real > &  q_surf_arr,
const amrex::Array4< const amrex::Real > &  u_vfrac_arr,
const amrex::Array4< const amrex::Real > &  v_vfrac_arr 
) const
inline
319  {
320  amrex::Real rho = cons_arr(i,j,k,Rho_comp);
321  amrex::Real qv = cons_arr(i,j,k,RhoQ1_comp) / rho;
322 
323  // Volume-weighted average of x-face velocities to cell center
324  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k);
325  amrex::Real velx = (u_vfrac_sum > eps) ?
326  (velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k))
327  / u_vfrac_sum : zero;
328 
329  // Volume-weighted average of y-face velocities to cell center
330  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k);
331  amrex::Real vely = (v_vfrac_sum > eps) ?
332  (vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k))
333  / v_vfrac_sum : zero;
334 
335  amrex::Real qv_mean = qvm_arr(i,j,0);
336  amrex::Real ustar = u_star_arr(i,j,k);
337  amrex::Real qstar = q_star_arr(i,j,k);
338  amrex::Real qv_surf = q_surf_arr(i,j,k);
339  amrex::Real wsp_mean = umm_arr(i,j,0);
340  wsp_mean = std::max(wsp_mean, WSMIN);
341 
342  amrex::Real wsp = std::sqrt(velx*velx+vely*vely);
343  amrex::Real num1 = wsp * (qv_mean-qv_surf);
344  amrex::Real num2 = wsp_mean * (qv-qv_mean);
345 
346  // NOTE: this is rho*<Qv'w'> = -K dQvdz
347  amrex::Real moflux = (std::abs(qstar) > eps) ?
348  -rho*qstar*ustar*(num1+num2)/((qv_mean-qv_surf)*wsp_mean) : zero;
349 
350  return moflux;
351  }
zero
constexpr amrex::Real zero
Definition: ERF_Constants.H:8
Rho_comp
#define Rho_comp
Definition: ERF_IndexDefines.H:36
RhoQ1_comp
#define RhoQ1_comp
Definition: ERF_IndexDefines.H:42
rho
rho
Definition: ERF_InitCustomPert_Bubble.H:106
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19
MicVar_Kess::qv
@ qv
Definition: ERF_Kessler.H:29
moeng_flux_eb::WSMIN
const amrex::Real WSMIN
Definition: ERF_EBMOSTStress.H:839
moeng_flux_eb::eps
const amrex::Real eps
Definition: ERF_EBMOSTStress.H:837

◆ compute_t_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux_eb::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 > &  velz_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::Array4< const amrex::Real > &  u_vfrac_arr,
const amrex::Array4< const amrex::Real > &  v_vfrac_arr,
const amrex::Array4< const amrex::Real > &  w_vfrac_arr,
const amrex::Array4< const amrex::Real > &  bnorm_arr 
) const
inline
372  {
373  amrex::Real rho = cons_arr(i,j,k,Rho_comp);
374  amrex::Real theta = cons_arr(i,j,k,RhoTheta_comp) / rho;
375 
376  // Volume-weighted average of x-face velocities to cell center
377  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k);
378  amrex::Real velx = (u_vfrac_sum > eps) ?
379  (velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k))
380  / u_vfrac_sum : zero;
381 
382  // Volume-weighted average of y-face velocities to cell center
383  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k);
384  amrex::Real vely = (v_vfrac_sum > eps) ?
385  (vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k))
386  / v_vfrac_sum : zero;
387 
388  // Volume-weighted average of z-face velocities to cell center
389  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1);
390  amrex::Real velz = (w_vfrac_sum > eps) ?
391  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1))
392  / w_vfrac_sum : zero;
393 
394  // Get boundary normal components
395  amrex::Real nx = bnorm_arr(i,j,k,0);
396  amrex::Real ny = bnorm_arr(i,j,k,1);
397  amrex::Real nz = bnorm_arr(i,j,k,2);
398 
399  // Project velocity onto tangent plane (remove normal component)
400  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
401  amrex::Real velx_tangent = velx - v_dot_n * nx;
402  amrex::Real vely_tangent = vely - v_dot_n * ny;
403 
404  amrex::Real theta_mean = tm_arr(i,j,0);
405  amrex::Real ustar = u_star_arr(i,j,k);
406  amrex::Real tstar = t_star_arr(i,j,k);
407  amrex::Real theta_surf = t_surf_arr(i,j,k);
408  amrex::Real wsp_mean = umm_arr(i,j,0);
409  wsp_mean = std::max(wsp_mean, WSMIN);
410 
411  // Use tangential velocity magnitude instead of Cartesian
412  amrex::Real wsp = std::sqrt(velx_tangent*velx_tangent+vely_tangent*vely_tangent);
413  amrex::Real num1 = wsp * (theta_mean-theta_surf);
414  amrex::Real num2 = wsp_mean * (theta-theta_mean);
415 
416  // NOTE: this is rho*<T'w'> = -K dTdz
417  amrex::Real moflux = (std::abs(tstar) > eps) ?
418  -rho*tstar*ustar*(num1+num2)/((theta_mean-theta_surf)*wsp_mean) : zero;
419 
420  return moflux;
421  }
RhoTheta_comp
#define RhoTheta_comp
Definition: ERF_IndexDefines.H:37
LsmVar_MM5::theta
@ theta
Definition: ERF_MM5.H:20

◆ compute_u_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux_eb::compute_u_flux ( int  i,
int  j,
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 > &  velz_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::Array4< const amrex::Real > &  u_vfrac_arr,
const amrex::Array4< const amrex::Real > &  v_vfrac_arr,
const amrex::Array4< const amrex::Real > &  w_vfrac_arr,
const amrex::Array4< const amrex::Real > &  cc_vfrac_arr,
const amrex::Array4< const amrex::EBCellFlag > &  cc_flag_arr,
const amrex::Array4< const amrex::Real > &  bnorm_arr,
int  idir = 0 
) const
inline
443  {
444  amrex::Real velx, vely, rho, ustar, wsp_mean;
445  amrex::Real velx_tangent, vely_tangent;
446 
447  if (idir == 0) {
448  // x-face: average to x-face
449  velx = velx_arr(i,j,k);
450 
451  // Volume-weighted average of y-face velocities to x-face
452  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k) +
453  v_vfrac_arr(i-1,j,k) + v_vfrac_arr(i-1,j+1,k);
454  vely = (v_vfrac_sum > eps) ?
455  (vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k) +
456  vely_arr(i-1,j,k) * v_vfrac_arr(i-1,j,k) + vely_arr(i-1,j+1,k) * v_vfrac_arr(i-1,j+1,k))
457  / v_vfrac_sum : zero;
458 
459  // Volume-weighted average of z-face velocities to x-face
460  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1) +
461  w_vfrac_arr(i-1,j,k) + w_vfrac_arr(i-1,j,k+1);
462  amrex::Real velz = (w_vfrac_sum > eps) ?
463  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1) +
464  velz_arr(i-1,j,k) * w_vfrac_arr(i-1,j,k) + velz_arr(i-1,j,k+1) * w_vfrac_arr(i-1,j,k+1))
465  / w_vfrac_sum : zero;
466 
467  // Get boundary normal at x-face (already at correct staggered location)
468  amrex::Real nx = bnorm_arr(i,j,k,0);
469  amrex::Real ny = bnorm_arr(i,j,k,1);
470  amrex::Real nz = bnorm_arr(i,j,k,2);
471 
472  // Project velocity onto tangent plane
473  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
474  velx_tangent = velx - v_dot_n * nx;
475  vely_tangent = vely - v_dot_n * ny;
476 
477  // Volume-weighted average of cell-centered density to x-face
478  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i-1,j,k) + cc_vfrac_arr(i,j,k);
479  rho = (cc_vfrac_sum > eps) ?
480  (cons_arr(i-1,j,k,Rho_comp) * cc_vfrac_arr(i-1,j,k) + cons_arr(i,j,k,Rho_comp) * cc_vfrac_arr(i,j,k))
481  / cc_vfrac_sum : zero;
482 
483  // Average cell-centered u_star and wsp_mean to x-face, using only valid (SingleValued) cells
484  bool low_valid = cc_flag_arr(i-1,j,k).isSingleValued();
485  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
486 
487  if (low_valid && high_valid) {
488  ustar = myhalf * (u_star_arr(i-1,j,k) + u_star_arr(i,j,k));
489  wsp_mean = myhalf * (umm_arr(i-1,j,0) + umm_arr(i,j,0));
490  } else if (low_valid) {
491  ustar = u_star_arr(i-1,j,k);
492  wsp_mean = umm_arr(i-1,j,0);
493  } else if (high_valid) {
494  ustar = u_star_arr(i,j,k);
495  wsp_mean = umm_arr(i,j,0);
496  } else {
497  ustar = zero;
498  wsp_mean = WSMIN;
499  }
500 
501  } else if (idir == 1) {
502  // y-face: average to y-face
503  vely = vely_arr(i,j,k);
504 
505  // Volume-weighted average of x-face velocities to y-face
506  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k) +
507  u_vfrac_arr(i,j-1,k) + u_vfrac_arr(i+1,j-1,k);
508  velx = (u_vfrac_sum > eps) ?
509  (velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k) +
510  velx_arr(i,j-1,k) * u_vfrac_arr(i,j-1,k) + velx_arr(i+1,j-1,k) * u_vfrac_arr(i+1,j-1,k))
511  / u_vfrac_sum : zero;
512 
513  // Volume-weighted average of z-face velocities to y-face
514  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1) +
515  w_vfrac_arr(i,j-1,k) + w_vfrac_arr(i,j-1,k+1);
516  amrex::Real velz = (w_vfrac_sum > eps) ?
517  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1) +
518  velz_arr(i,j-1,k) * w_vfrac_arr(i,j-1,k) + velz_arr(i,j-1,k+1) * w_vfrac_arr(i,j-1,k+1))
519  / w_vfrac_sum : zero;
520 
521  // Get boundary normal at y-face (already at correct staggered location)
522  amrex::Real nx = bnorm_arr(i,j,k,0);
523  amrex::Real ny = bnorm_arr(i,j,k,1);
524  amrex::Real nz = bnorm_arr(i,j,k,2);
525 
526  // Project velocity onto tangent plane
527  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
528  velx_tangent = velx - v_dot_n * nx;
529  vely_tangent = vely - v_dot_n * ny;
530 
531  // Volume-weighted average of cell-centered density to y-face
532  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i,j-1,k) + cc_vfrac_arr(i,j,k);
533  rho = (cc_vfrac_sum > eps) ?
534  (cons_arr(i,j-1,k,Rho_comp) * cc_vfrac_arr(i,j-1,k) + cons_arr(i,j,k,Rho_comp) * cc_vfrac_arr(i,j,k))
535  / cc_vfrac_sum : zero;
536 
537  // Average cell-centered u_star and wsp_mean to y-face, using only valid (SingleValued) cells
538  bool low_valid = cc_flag_arr(i,j-1,k).isSingleValued();
539  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
540 
541  if (low_valid && high_valid) {
542  ustar = myhalf * (u_star_arr(i,j-1,k) + u_star_arr(i,j,k));
543  wsp_mean = myhalf * (umm_arr(i,j-1,0) + umm_arr(i,j,0));
544  } else if (low_valid) {
545  ustar = u_star_arr(i,j-1,k);
546  wsp_mean = umm_arr(i,j-1,0);
547  } else if (high_valid) {
548  ustar = u_star_arr(i,j,k);
549  wsp_mean = umm_arr(i,j,0);
550  } else {
551  ustar = zero;
552  wsp_mean = WSMIN;
553  }
554 
555  } else {
556  // z-face: average to z-face
557  // Volume-weighted average of x-face velocities to z-face
558  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k-1) + u_vfrac_arr(i+1,j,k-1) +
559  u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k);
560  velx = (u_vfrac_sum > eps) ?
561  (velx_arr(i,j,k-1) * u_vfrac_arr(i,j,k-1) + velx_arr(i+1,j,k-1) * u_vfrac_arr(i+1,j,k-1) +
562  velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k))
563  / u_vfrac_sum : zero;
564 
565  // Volume-weighted average of y-face velocities to z-face
566  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k-1) + v_vfrac_arr(i,j+1,k-1) +
567  v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k);
568  vely = (v_vfrac_sum > eps) ?
569  (vely_arr(i,j,k-1) * v_vfrac_arr(i,j,k-1) + vely_arr(i,j+1,k-1) * v_vfrac_arr(i,j+1,k-1) +
570  vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k))
571  / v_vfrac_sum : zero;
572 
573  // z-velocity is already at z-face
574  amrex::Real velz = velz_arr(i,j,k);
575 
576  // Get boundary normal at z-face (already at correct staggered location)
577  amrex::Real nx = bnorm_arr(i,j,k,0);
578  amrex::Real ny = bnorm_arr(i,j,k,1);
579  amrex::Real nz = bnorm_arr(i,j,k,2);
580 
581  // Project velocity onto tangent plane
582  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
583  velx_tangent = velx - v_dot_n * nx;
584  vely_tangent = vely - v_dot_n * ny;
585 
586  // Volume-weighted average of cell-centered density to z-face
587  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i,j,k-1) + cc_vfrac_arr(i,j,k);
588  rho = (cc_vfrac_sum > eps) ?
589  (cons_arr(i,j,k-1,Rho_comp) * cc_vfrac_arr(i,j,k-1) + cons_arr(i,j,k,Rho_comp) * cc_vfrac_arr(i,j,k))
590  / cc_vfrac_sum : zero;
591 
592  // Average cell-centered u_star and wsp_mean to z-face, using only valid (SingleValued) cells
593  bool low_valid = cc_flag_arr(i,j,k-1).isSingleValued();
594  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
595 
596  if (low_valid && high_valid) {
597  ustar = myhalf * (u_star_arr(i,j,k-1) + u_star_arr(i,j,k));
598  wsp_mean = myhalf * (umm_arr(i,j,0) + umm_arr(i,j,0));
599  } else if (low_valid) {
600  ustar = u_star_arr(i,j,k-1);
601  wsp_mean = umm_arr(i,j,0);
602  } else if (high_valid) {
603  ustar = u_star_arr(i,j,k);
604  wsp_mean = umm_arr(i,j,0);
605  } else {
606  ustar = zero;
607  wsp_mean = WSMIN;
608  }
609  }
610 
611  wsp_mean = std::max(wsp_mean, WSMIN);
612  amrex::Real umean = um_arr(i,j,0);
613 
614  // Note: The surface mean shear stress is decomposed into tau_xz by
615  // multiplying the modeled shear stress (rho*ustar^2) with
616  // a factor of umean/wsp_mean for directionality; this factor
617  // modifies the denominator from what is in Moeng amrex::Real(1984.)
618  amrex::Real wsp = std::sqrt(velx_tangent*velx_tangent+vely_tangent*vely_tangent);
619  amrex::Real num1 = wsp * umean;
620  amrex::Real num2 = wsp_mean * (velx_tangent-umean);
621 
622  // NOTE: this is rho*<u'w'> = -K dudz
623  amrex::Real stressx = -rho*ustar*ustar * (num1+num2)/(wsp_mean*wsp_mean);
624 
625  return stressx;
626  }
myhalf
constexpr amrex::Real myhalf
Definition: ERF_Constants.H:13

◆ compute_v_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real moeng_flux_eb::compute_v_flux ( int  i,
int  j,
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 > &  velz_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 amrex::Array4< const amrex::Real > &  u_vfrac_arr,
const amrex::Array4< const amrex::Real > &  v_vfrac_arr,
const amrex::Array4< const amrex::Real > &  w_vfrac_arr,
const amrex::Array4< const amrex::Real > &  cc_vfrac_arr,
const amrex::Array4< const amrex::EBCellFlag > &  cc_flag_arr,
const amrex::Array4< const amrex::Real > &  bnorm_arr,
int  idir = 0 
) const
inline
648  {
649  amrex::Real velx, vely, rho, ustar, wsp_mean;
650  amrex::Real velx_tangent, vely_tangent;
651 
652  if (idir == 0) {
653  // x-face: average from cells (i-1) and (i)
654  // Volume-weighted average of y-face velocities to x-face
655  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k) +
656  v_vfrac_arr(i-1,j,k) + v_vfrac_arr(i-1,j+1,k);
657  vely = (v_vfrac_sum > eps) ?
658  (vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k) +
659  vely_arr(i-1,j,k) * v_vfrac_arr(i-1,j,k) + vely_arr(i-1,j+1,k) * v_vfrac_arr(i-1,j+1,k))
660  / v_vfrac_sum : zero;
661 
662  velx = velx_arr(i,j,k);
663 
664  // Volume-weighted average of z-face velocities to x-face
665  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1) +
666  w_vfrac_arr(i-1,j,k) + w_vfrac_arr(i-1,j,k+1);
667  amrex::Real velz = (w_vfrac_sum > eps) ?
668  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1) +
669  velz_arr(i-1,j,k) * w_vfrac_arr(i-1,j,k) + velz_arr(i-1,j,k+1) * w_vfrac_arr(i-1,j,k+1))
670  / w_vfrac_sum : zero;
671 
672  // Get boundary normal at x-face (already at correct staggered location)
673  amrex::Real nx = bnorm_arr(i,j,k,0);
674  amrex::Real ny = bnorm_arr(i,j,k,1);
675  amrex::Real nz = bnorm_arr(i,j,k,2);
676 
677  // Project velocity onto tangent plane
678  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
679  velx_tangent = velx - v_dot_n * nx;
680  vely_tangent = vely - v_dot_n * ny;
681 
682  // Volume-weighted average of cell-centered density to x-face
683  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i-1,j,k) + cc_vfrac_arr(i,j,k);
684  rho = (cc_vfrac_sum > eps) ?
685  (cons_arr(i-1,j,k,Rho_comp) * cc_vfrac_arr(i-1,j,k) + cons_arr(i,j,k,Rho_comp) * cc_vfrac_arr(i,j,k))
686  / cc_vfrac_sum : zero;
687 
688  // Average cell-centered u_star and wsp_mean to x-face, using only valid (SingleValued) cells
689  bool low_valid = cc_flag_arr(i-1,j,k).isSingleValued();
690  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
691 
692  if (low_valid && high_valid) {
693  ustar = myhalf * (u_star_arr(i-1,j,k) + u_star_arr(i,j,k));
694  wsp_mean = myhalf * (umm_arr(i-1,j,0) + umm_arr(i,j,0));
695  } else if (low_valid) {
696  ustar = u_star_arr(i-1,j,k);
697  wsp_mean = umm_arr(i-1,j,0);
698  } else if (high_valid) {
699  ustar = u_star_arr(i,j,k);
700  wsp_mean = umm_arr(i,j,0);
701  } else {
702  ustar = zero;
703  wsp_mean = WSMIN;
704  }
705 
706  } else if (idir == 1) {
707  // y-face: average from cells (j-1) and (j)
708  // Volume-weighted average of x-face velocities to y-face
709  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k) +
710  u_vfrac_arr(i,j-1,k) + u_vfrac_arr(i+1,j-1,k);
711  velx = (u_vfrac_sum > eps) ?
712  (velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k) +
713  velx_arr(i,j-1,k) * u_vfrac_arr(i,j-1,k) + velx_arr(i+1,j-1,k) * u_vfrac_arr(i+1,j-1,k))
714  / u_vfrac_sum : zero;
715 
716  vely = vely_arr(i,j,k);
717 
718  // Volume-weighted average of z-face velocities to y-face
719  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1) +
720  w_vfrac_arr(i,j-1,k) + w_vfrac_arr(i,j-1,k+1);
721  amrex::Real velz = (w_vfrac_sum > eps) ?
722  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1) +
723  velz_arr(i,j-1,k) * w_vfrac_arr(i,j-1,k) + velz_arr(i,j-1,k+1) * w_vfrac_arr(i,j-1,k+1))
724  / w_vfrac_sum : zero;
725 
726  // Get boundary normal at y-face (already at correct staggered location)
727  amrex::Real nx = bnorm_arr(i,j,k,0);
728  amrex::Real ny = bnorm_arr(i,j,k,1);
729  amrex::Real nz = bnorm_arr(i,j,k,2);
730 
731  // Project velocity onto tangent plane
732  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
733  velx_tangent = velx - v_dot_n * nx;
734  vely_tangent = vely - v_dot_n * ny;
735 
736  // Volume-weighted average of cell-centered density to y-face
737  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i,j-1,k) + cc_vfrac_arr(i,j,k);
738  rho = (cc_vfrac_sum > eps) ?
739  (cons_arr(i,j-1,k,Rho_comp) * cc_vfrac_arr(i,j-1,k) + cons_arr(i,j,k,Rho_comp) * cc_vfrac_arr(i,j,k))
740  / cc_vfrac_sum : zero;
741 
742  // Average cell-centered u_star and wsp_mean to y-face, using only valid (SingleValued) cells
743  bool low_valid = cc_flag_arr(i,j-1,k).isSingleValued();
744  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
745 
746  if (low_valid && high_valid) {
747  ustar = myhalf * (u_star_arr(i,j-1,k) + u_star_arr(i,j,k));
748  wsp_mean = myhalf * (umm_arr(i,j-1,0) + umm_arr(i,j,0));
749  } else if (low_valid) {
750  ustar = u_star_arr(i,j-1,k);
751  wsp_mean = umm_arr(i,j-1,0);
752  } else if (high_valid) {
753  ustar = u_star_arr(i,j,k);
754  wsp_mean = umm_arr(i,j,0);
755  } else {
756  ustar = zero;
757  wsp_mean = WSMIN;
758  }
759 
760  } else {
761  // z-face: average from cells (k-1) and (k)
762  // Volume-weighted average of x-face velocities to z-face
763  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k-1) + u_vfrac_arr(i+1,j,k-1) +
764  u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k);
765  velx = (u_vfrac_sum > eps) ?
766  (velx_arr(i,j,k-1) * u_vfrac_arr(i,j,k-1) + velx_arr(i+1,j,k-1) * u_vfrac_arr(i+1,j,k-1) +
767  velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k))
768  / u_vfrac_sum : zero;
769 
770  // Volume-weighted average of y-face velocities to z-face
771  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k-1) + v_vfrac_arr(i,j+1,k-1) +
772  v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k);
773  vely = (v_vfrac_sum > eps) ?
774  (vely_arr(i,j,k-1) * v_vfrac_arr(i,j,k-1) + vely_arr(i,j+1,k-1) * v_vfrac_arr(i,j+1,k-1) +
775  vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k))
776  / v_vfrac_sum : zero;
777 
778  // z-velocity is already at z-face
779  amrex::Real velz = velz_arr(i,j,k);
780 
781  // Get boundary normal at z-face (already at correct staggered location)
782  amrex::Real nx = bnorm_arr(i,j,k,0);
783  amrex::Real ny = bnorm_arr(i,j,k,1);
784  amrex::Real nz = bnorm_arr(i,j,k,2);
785 
786  // Project velocity onto tangent plane
787  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
788  velx_tangent = velx - v_dot_n * nx;
789  vely_tangent = vely - v_dot_n * ny;
790 
791  // Volume-weighted average of cell-centered density to z-face
792  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i,j,k-1) + cc_vfrac_arr(i,j,k);
793  rho = (cc_vfrac_sum > eps) ?
794  (cons_arr(i,j,k-1,Rho_comp) * cc_vfrac_arr(i,j,k-1) + cons_arr(i,j,k,Rho_comp) * cc_vfrac_arr(i,j,k))
795  / cc_vfrac_sum : zero;
796 
797  // Average cell-centered u_star and wsp_mean to z-face, using only valid (SingleValued) cells
798  bool low_valid = cc_flag_arr(i,j,k-1).isSingleValued();
799  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
800 
801  if (low_valid && high_valid) {
802  ustar = myhalf * (u_star_arr(i,j,k-1) + u_star_arr(i,j,k));
803  wsp_mean = myhalf * (umm_arr(i,j,0) + umm_arr(i,j,0));
804  } else if (low_valid) {
805  ustar = u_star_arr(i,j,k-1);
806  wsp_mean = umm_arr(i,j,0);
807  } else if (high_valid) {
808  ustar = u_star_arr(i,j,k);
809  wsp_mean = umm_arr(i,j,0);
810  } else {
811  ustar = zero;
812  wsp_mean = WSMIN;
813  }
814  }
815 
816  wsp_mean = std::max(wsp_mean, WSMIN);
817  amrex::Real vmean = vm_arr(i,j,0);
818 
819  // Note: The surface mean shear stress is decomposed into tau_yz by
820  // multiplying the modeled shear stress (rho*ustar^2) with
821  // a factor of vmean/wsp_mean for directionality; this factor
822  // modifies the denominator from what is in Moeng amrex::Real(1984.)
823  amrex::Real wsp = std::sqrt(velx_tangent*velx_tangent+vely_tangent*vely_tangent);
824  amrex::Real num1 = wsp * vmean;
825  amrex::Real num2 = wsp_mean * (vely_tangent-vmean);
826 
827  // NOTE: this is rho*<v'w'> = -K dvdz
828  amrex::Real stressy = -rho*ustar*ustar * (num1+num2)/(wsp_mean*wsp_mean);
829 
830  return stressy;
831  }

Member Data Documentation

◆ eps

const amrex::Real moeng_flux_eb::eps = amrex::Real(1e-12)
private

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

◆ WSMIN

const amrex::Real moeng_flux_eb::WSMIN = amrex::Real(0.1)
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: