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 = 1e-15
 
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
299 {}

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

◆ 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
370  {
371  amrex::Real rho = cons_arr(i,j,k,Rho_comp);
372  amrex::Real theta = cons_arr(i,j,k,RhoTheta_comp) / rho;
373 
374  // Volume-weighted average of x-face velocities to cell center
375  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k);
376  amrex::Real velx = (u_vfrac_sum > eps) ?
377  (velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k))
378  / u_vfrac_sum : zero;
379 
380  // Volume-weighted average of y-face velocities to cell center
381  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k);
382  amrex::Real vely = (v_vfrac_sum > eps) ?
383  (vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k))
384  / v_vfrac_sum : zero;
385 
386  // Volume-weighted average of z-face velocities to cell center
387  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1);
388  amrex::Real velz = (w_vfrac_sum > eps) ?
389  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1))
390  / w_vfrac_sum : zero;
391 
392  // Get boundary normal components
393  amrex::Real nx = bnorm_arr(i,j,k,0);
394  amrex::Real ny = bnorm_arr(i,j,k,1);
395  amrex::Real nz = bnorm_arr(i,j,k,2);
396 
397  // Project velocity onto tangent plane (remove normal component)
398  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
399  amrex::Real velx_tangent = velx - v_dot_n * nx;
400  amrex::Real vely_tangent = vely - v_dot_n * ny;
401 
402  amrex::Real theta_mean = tm_arr(i,j,0);
403  amrex::Real ustar = u_star_arr(i,j,k);
404  amrex::Real tstar = t_star_arr(i,j,k);
405  amrex::Real theta_surf = t_surf_arr(i,j,k);
406  amrex::Real wsp_mean = umm_arr(i,j,0);
407  wsp_mean = std::max(wsp_mean, WSMIN);
408 
409  // Use tangential velocity magnitude instead of Cartesian
410  amrex::Real wsp = sqrt(velx_tangent*velx_tangent+vely_tangent*vely_tangent);
411  amrex::Real num1 = wsp * (theta_mean-theta_surf);
412  amrex::Real num2 = wsp_mean * (theta-theta_mean);
413 
414  // NOTE: this is rho*<T'w'> = -K dTdz
415  amrex::Real moflux = (std::abs(tstar) > eps) ?
416  -rho*tstar*ustar*(num1+num2)/((theta_mean-theta_surf)*wsp_mean) : zero;
417 
418  return moflux;
419  }
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
441  {
442  amrex::Real velx, vely, rho, ustar, wsp_mean;
443  amrex::Real velx_tangent, vely_tangent;
444 
445  if (idir == 0) {
446  // x-face: average to x-face
447  velx = velx_arr(i,j,k);
448 
449  // Volume-weighted average of y-face velocities to x-face
450  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k) +
451  v_vfrac_arr(i-1,j,k) + v_vfrac_arr(i-1,j+1,k);
452  vely = (v_vfrac_sum > eps) ?
453  (vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k) +
454  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))
455  / v_vfrac_sum : zero;
456 
457  // Volume-weighted average of z-face velocities to x-face
458  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1) +
459  w_vfrac_arr(i-1,j,k) + w_vfrac_arr(i-1,j,k+1);
460  amrex::Real velz = (w_vfrac_sum > eps) ?
461  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1) +
462  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))
463  / w_vfrac_sum : zero;
464 
465  // Get boundary normal at x-face (already at correct staggered location)
466  amrex::Real nx = bnorm_arr(i,j,k,0);
467  amrex::Real ny = bnorm_arr(i,j,k,1);
468  amrex::Real nz = bnorm_arr(i,j,k,2);
469 
470  // Project velocity onto tangent plane
471  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
472  velx_tangent = velx - v_dot_n * nx;
473  vely_tangent = vely - v_dot_n * ny;
474 
475  // Volume-weighted average of cell-centered density to x-face
476  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i-1,j,k) + cc_vfrac_arr(i,j,k);
477  rho = (cc_vfrac_sum > eps) ?
478  (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))
479  / cc_vfrac_sum : zero;
480 
481  // Average cell-centered u_star and wsp_mean to x-face, using only valid (SingleValued) cells
482  bool low_valid = cc_flag_arr(i-1,j,k).isSingleValued();
483  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
484 
485  if (low_valid && high_valid) {
486  ustar = myhalf * (u_star_arr(i-1,j,k) + u_star_arr(i,j,k));
487  wsp_mean = myhalf * (umm_arr(i-1,j,0) + umm_arr(i,j,0));
488  } else if (low_valid) {
489  ustar = u_star_arr(i-1,j,k);
490  wsp_mean = umm_arr(i-1,j,0);
491  } else if (high_valid) {
492  ustar = u_star_arr(i,j,k);
493  wsp_mean = umm_arr(i,j,0);
494  } else {
495  ustar = zero;
496  wsp_mean = WSMIN;
497  }
498 
499  } else if (idir == 1) {
500  // y-face: average to y-face
501  vely = vely_arr(i,j,k);
502 
503  // Volume-weighted average of x-face velocities to y-face
504  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k) +
505  u_vfrac_arr(i,j-1,k) + u_vfrac_arr(i+1,j-1,k);
506  velx = (u_vfrac_sum > eps) ?
507  (velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k) +
508  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))
509  / u_vfrac_sum : zero;
510 
511  // Volume-weighted average of z-face velocities to y-face
512  amrex::Real w_vfrac_sum = w_vfrac_arr(i,j,k) + w_vfrac_arr(i,j,k+1) +
513  w_vfrac_arr(i,j-1,k) + w_vfrac_arr(i,j-1,k+1);
514  amrex::Real velz = (w_vfrac_sum > eps) ?
515  (velz_arr(i,j,k) * w_vfrac_arr(i,j,k) + velz_arr(i,j,k+1) * w_vfrac_arr(i,j,k+1) +
516  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))
517  / w_vfrac_sum : zero;
518 
519  // Get boundary normal at y-face (already at correct staggered location)
520  amrex::Real nx = bnorm_arr(i,j,k,0);
521  amrex::Real ny = bnorm_arr(i,j,k,1);
522  amrex::Real nz = bnorm_arr(i,j,k,2);
523 
524  // Project velocity onto tangent plane
525  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
526  velx_tangent = velx - v_dot_n * nx;
527  vely_tangent = vely - v_dot_n * ny;
528 
529  // Volume-weighted average of cell-centered density to y-face
530  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i,j-1,k) + cc_vfrac_arr(i,j,k);
531  rho = (cc_vfrac_sum > eps) ?
532  (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))
533  / cc_vfrac_sum : zero;
534 
535  // Average cell-centered u_star and wsp_mean to y-face, using only valid (SingleValued) cells
536  bool low_valid = cc_flag_arr(i,j-1,k).isSingleValued();
537  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
538 
539  if (low_valid && high_valid) {
540  ustar = myhalf * (u_star_arr(i,j-1,k) + u_star_arr(i,j,k));
541  wsp_mean = myhalf * (umm_arr(i,j-1,0) + umm_arr(i,j,0));
542  } else if (low_valid) {
543  ustar = u_star_arr(i,j-1,k);
544  wsp_mean = umm_arr(i,j-1,0);
545  } else if (high_valid) {
546  ustar = u_star_arr(i,j,k);
547  wsp_mean = umm_arr(i,j,0);
548  } else {
549  ustar = zero;
550  wsp_mean = WSMIN;
551  }
552 
553  } else {
554  // z-face: average to z-face
555  // Volume-weighted average of x-face velocities to z-face
556  amrex::Real u_vfrac_sum = u_vfrac_arr(i,j,k-1) + u_vfrac_arr(i+1,j,k-1) +
557  u_vfrac_arr(i,j,k) + u_vfrac_arr(i+1,j,k);
558  velx = (u_vfrac_sum > eps) ?
559  (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) +
560  velx_arr(i,j,k) * u_vfrac_arr(i,j,k) + velx_arr(i+1,j,k) * u_vfrac_arr(i+1,j,k))
561  / u_vfrac_sum : zero;
562 
563  // Volume-weighted average of y-face velocities to z-face
564  amrex::Real v_vfrac_sum = v_vfrac_arr(i,j,k-1) + v_vfrac_arr(i,j+1,k-1) +
565  v_vfrac_arr(i,j,k) + v_vfrac_arr(i,j+1,k);
566  vely = (v_vfrac_sum > eps) ?
567  (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) +
568  vely_arr(i,j,k) * v_vfrac_arr(i,j,k) + vely_arr(i,j+1,k) * v_vfrac_arr(i,j+1,k))
569  / v_vfrac_sum : zero;
570 
571  // z-velocity is already at z-face
572  amrex::Real velz = velz_arr(i,j,k);
573 
574  // Get boundary normal at z-face (already at correct staggered location)
575  amrex::Real nx = bnorm_arr(i,j,k,0);
576  amrex::Real ny = bnorm_arr(i,j,k,1);
577  amrex::Real nz = bnorm_arr(i,j,k,2);
578 
579  // Project velocity onto tangent plane
580  amrex::Real v_dot_n = velx*nx + vely*ny + velz*nz;
581  velx_tangent = velx - v_dot_n * nx;
582  vely_tangent = vely - v_dot_n * ny;
583 
584  // Volume-weighted average of cell-centered density to z-face
585  amrex::Real cc_vfrac_sum = cc_vfrac_arr(i,j,k-1) + cc_vfrac_arr(i,j,k);
586  rho = (cc_vfrac_sum > eps) ?
587  (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))
588  / cc_vfrac_sum : zero;
589 
590  // Average cell-centered u_star and wsp_mean to z-face, using only valid (SingleValued) cells
591  bool low_valid = cc_flag_arr(i,j,k-1).isSingleValued();
592  bool high_valid = cc_flag_arr(i,j,k).isSingleValued();
593 
594  if (low_valid && high_valid) {
595  ustar = myhalf * (u_star_arr(i,j,k-1) + u_star_arr(i,j,k));
596  wsp_mean = myhalf * (umm_arr(i,j,0) + umm_arr(i,j,0));
597  } else if (low_valid) {
598  ustar = u_star_arr(i,j,k-1);
599  wsp_mean = umm_arr(i,j,0);
600  } else if (high_valid) {
601  ustar = u_star_arr(i,j,k);
602  wsp_mean = umm_arr(i,j,0);
603  } else {
604  ustar = zero;
605  wsp_mean = WSMIN;
606  }
607  }
608 
609  wsp_mean = std::max(wsp_mean, WSMIN);
610  amrex::Real umean = um_arr(i,j,0);
611 
612  // Note: The surface mean shear stress is decomposed into tau_xz by
613  // multiplying the modeled shear stress (rho*ustar^2) with
614  // a factor of umean/wsp_mean for directionality; this factor
615  // modifies the denominator from what is in Moeng amrex::Real(1984.)
616  amrex::Real wsp = sqrt(velx_tangent*velx_tangent+vely_tangent*vely_tangent);
617  amrex::Real num1 = wsp * umean;
618  amrex::Real num2 = wsp_mean * (velx_tangent-umean);
619 
620  // NOTE: this is rho*<u'w'> = -K dudz
621  amrex::Real stressx = -rho*ustar*ustar * (num1+num2)/(wsp_mean*wsp_mean);
622 
623  return stressx;
624  }
myhalf
constexpr amrex::Real myhalf
Definition: ERF_Constants.H:11

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

Member Data Documentation

◆ eps

const amrex::Real moeng_flux_eb::eps = 1e-15
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: