ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Functions
ERF_EddyViscosity.H File Reference
#include "AMReX_BCRec.H"
#include "ERF_SurfaceLayer.H"
#include "ERF_DataStruct.H"
#include "ERF_IndexDefines.H"
#include "ERF_Constants.H"
#include "ERF_EOS.H"
#include "ERF_EB.H"
Include dependency graph for ERF_EddyViscosity.H:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions

void ComputeTurbulentViscosity (amrex::Real dt, const amrex::MultiFab &xvel, const amrex::MultiFab &yvel, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Tau_lev, amrex::MultiFab &cons_in, const amrex::MultiFab &wdist, amrex::MultiFab &eddyViscosity, amrex::MultiFab &Hfx1, amrex::MultiFab &Hfx2, amrex::MultiFab &Hfx3, amrex::MultiFab &Diss, const amrex::Geometry &geom, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, const std::unique_ptr< amrex::MultiFab > &z_phys_nd, const SolverChoice &solverChoice, std::unique_ptr< SurfaceLayer > &SurfLayer, const amrex::MultiFab *z_0, const bool &use_terrain_fitted_coords, const bool &use_moisture, int level, const amrex::BCRec *bc_ptr, const eb_ &ebfact, bool vert_only=false)
 
void ComputeTurbulentViscosityLES (amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Tau_lev, const amrex::MultiFab &cons_in, amrex::MultiFab &eddyViscosity, amrex::MultiFab &Hfx1, amrex::MultiFab &Hfx2, amrex::MultiFab &Hfx3, amrex::MultiFab &Diss, const amrex::Geometry &geom, bool use_terrain_fitted_coords, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, const std::unique_ptr< amrex::MultiFab > &z_phys_nd, const TurbChoice &turbChoice, const amrex::Real const_grav, std::unique_ptr< SurfaceLayer > &SurfLayer, const MoistureComponentIndices &moisture_indices, const amrex::MultiFab *xvel=nullptr, const amrex::MultiFab *yvel=nullptr)
 
void ComputeTurbulentViscosityLES_EB (amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Tau_lev, const amrex::MultiFab &cons_in, amrex::MultiFab &eddyViscosity, amrex::MultiFab &Hfx1, amrex::MultiFab &Hfx2, amrex::MultiFab &Hfx3, const amrex::Geometry &geom, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, const TurbChoice &turbChoice, const amrex::Real const_grav, const SolverChoice &solverChoice, std::unique_ptr< SurfaceLayer > &SurfLayer, const MoistureComponentIndices &moisture_indices, const eb_ &ebfact, const amrex::MultiFab *xvel=nullptr, const amrex::MultiFab *yvel=nullptr)
 
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real ComputeSmnSmn (int &i, int &j, int &k, const amrex::Array4< amrex::Real const > &tau11, const amrex::Array4< amrex::Real const > &tau22, const amrex::Array4< amrex::Real const > &tau33, const amrex::Array4< amrex::Real const > &tau12, const amrex::Array4< amrex::Real const > &tau13, const amrex::Array4< amrex::Real const > &tau23)
 
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real ComputeSmnSmn_EB (int &i, int &j, int &k, const amrex::Array4< amrex::Real const > &tau11, const amrex::Array4< amrex::Real const > &tau22, const amrex::Array4< amrex::Real const > &tau33, const amrex::Array4< amrex::Real const > &tau12, const amrex::Array4< amrex::Real const > &tau13, const amrex::Array4< amrex::Real const > &tau23, const amrex::Array4< const amrex::EBCellFlag > &c_cflag, const amrex::Array4< const amrex::EBCellFlag > &u_cflag, const amrex::Array4< const amrex::EBCellFlag > &v_cflag, const amrex::Array4< const amrex::EBCellFlag > &w_cflag)
 
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real ComputeSmnSmn2D (int &i, int &j, int &k, const amrex::Array4< amrex::Real const > &tau11, const amrex::Array4< amrex::Real const > &tau22, const amrex::Array4< amrex::Real const > &tau12)
 

Function Documentation

◆ ComputeSmnSmn()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real ComputeSmnSmn ( int &  i,
int &  j,
int &  k,
const amrex::Array4< amrex::Real const > &  tau11,
const amrex::Array4< amrex::Real const > &  tau22,
const amrex::Array4< amrex::Real const > &  tau33,
const amrex::Array4< amrex::Real const > &  tau12,
const amrex::Array4< amrex::Real const > &  tau13,
const amrex::Array4< amrex::Real const > &  tau23 
)
90 {
91  amrex::Real s11bar = tau11(i,j,k);
92  amrex::Real s22bar = tau22(i,j,k);
93  amrex::Real s33bar = tau33(i,j,k);
94  amrex::Real s12bar = 0.25 * ( tau12(i , j , k ) + tau12(i , j+1, k )
95  + tau12(i+1, j , k ) + tau12(i+1, j+1, k ) );
96  amrex::Real s13bar = 0.25 * ( tau13(i , j , k ) + tau13(i , j , k+1)
97  + tau13(i+1, j , k ) + tau13(i+1, j , k+1) );
98  amrex::Real s23bar = 0.25 * ( tau23(i , j , k ) + tau23(i , j , k+1)
99  + tau23(i , j+1, k ) + tau23(i , j+1, k+1) );
100 
101  amrex::Real SmnSmn = s11bar*s11bar + s22bar*s22bar + s33bar*s33bar
102  + 2.0*s12bar*s12bar + 2.0*s13bar*s13bar + 2.0*s23bar*s23bar;
103 
104  return SmnSmn;
105 }
tau12
@ tau12
Definition: ERF_DataStruct.H:32
tau23
@ tau23
Definition: ERF_DataStruct.H:32
tau33
@ tau33
Definition: ERF_DataStruct.H:32
tau22
@ tau22
Definition: ERF_DataStruct.H:32
tau11
@ tau11
Definition: ERF_DataStruct.H:32
tau13
@ tau13
Definition: ERF_DataStruct.H:32
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19

Referenced by ComputeTurbulentViscosityLES(), ComputeTurbulentViscosityLES_EB(), and erf_make_tau_terms().

Here is the caller graph for this function:

◆ ComputeSmnSmn2D()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real ComputeSmnSmn2D ( int &  i,
int &  j,
int &  k,
const amrex::Array4< amrex::Real const > &  tau11,
const amrex::Array4< amrex::Real const > &  tau22,
const amrex::Array4< amrex::Real const > &  tau12 
)
208 {
209  amrex::Real sdiff = tau11(i,j,k) - tau22(i,j,k);
210  amrex::Real s12bar = 0.25 * ( tau12(i , j , k ) + tau12(i , j+1, k )
211  + tau12(i+1, j , k ) + tau12(i+1, j+1, k ) );
212  return 0.5 * sdiff*sdiff + 2.0*s12bar*s12bar;
213 }

Referenced by ComputeTurbulentViscosityLES().

Here is the caller graph for this function:

◆ ComputeSmnSmn_EB()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real ComputeSmnSmn_EB ( int &  i,
int &  j,
int &  k,
const amrex::Array4< amrex::Real const > &  tau11,
const amrex::Array4< amrex::Real const > &  tau22,
const amrex::Array4< amrex::Real const > &  tau33,
const amrex::Array4< amrex::Real const > &  tau12,
const amrex::Array4< amrex::Real const > &  tau13,
const amrex::Array4< amrex::Real const > &  tau23,
const amrex::Array4< const amrex::EBCellFlag > &  c_cflag,
const amrex::Array4< const amrex::EBCellFlag > &  u_cflag,
const amrex::Array4< const amrex::EBCellFlag > &  v_cflag,
const amrex::Array4< const amrex::EBCellFlag > &  w_cflag 
)
121 {
122  amrex::Real SmnSmn = 0.0;
123  if (!c_cflag(i,j,k).isCovered()) {
124  amrex::Real s11bar = tau11(i,j,k);
125  amrex::Real s22bar = tau22(i,j,k);
126  amrex::Real s33bar = tau33(i,j,k);
127  amrex::Real s12bar = 0.0;
128  amrex::Real s13bar = 0.0;
129  amrex::Real s23bar = 0.0;
130 
131  amrex::Real count_s12 = 0.0;
132  if (!u_cflag(i,j,k).isCovered() || !v_cflag(i,j,k).isCovered()) {
133  s12bar += tau12(i,j,k);
134  count_s12 += 1.0;
135  }
136  if (!u_cflag(i+1,j,k).isCovered() || !v_cflag(i,j,k).isCovered()) {
137  s12bar += tau12(i+1,j,k);
138  count_s12 += 1.0;
139  }
140  if (!u_cflag(i,j,k).isCovered() || !v_cflag(i,j+1,k).isCovered()) {
141  s12bar += tau12(i,j+1,k);
142  count_s12 += 1.0;
143  }
144  if (!u_cflag(i+1,j,k).isCovered() || !v_cflag(i,j+1,k).isCovered()) {
145  s12bar += tau12(i+1,j+1,k);
146  count_s12 += 1.0;
147  }
148  if (count_s12 > 0.0) {
149  s12bar /= count_s12;
150  }
151 
152  amrex::Real count_s13 = 0.0;
153  if (!u_cflag(i,j,k).isCovered() || !w_cflag(i,j,k).isCovered()) {
154  s13bar += tau13(i, j, k);
155  count_s13 += 1.0;
156  }
157  if (!u_cflag(i,j,k).isCovered() || !w_cflag(i,j,k+1).isCovered()) {
158  s13bar += tau13(i, j, k+1);
159  count_s13 += 1.0;
160  }
161  if (!u_cflag(i+1,j,k).isCovered() || !w_cflag(i,j,k).isCovered()) {
162  s13bar += tau13(i+1, j, k);
163  count_s13 += 1.0;
164  }
165  if (!u_cflag(i+1,j,k).isCovered() || !w_cflag(i,j,k+1).isCovered()) {
166  s13bar += tau13(i+1, j, k+1);
167  count_s13 += 1.0;
168  }
169  if (count_s13 > 0.0) {
170  s13bar /= count_s13;
171  }
172 
173  amrex::Real count_s23 = 0.0;
174  if (!v_cflag(i,j,k).isCovered() || !w_cflag(i,j,k).isCovered()) {
175  s23bar += tau23(i, j, k);
176  count_s23 += 1.0;
177  }
178  if (!v_cflag(i,j,k).isCovered() || !w_cflag(i,j,k+1).isCovered()) {
179  s23bar += tau23(i, j, k+1);
180  count_s23 += 1.0;
181  }
182  if (!v_cflag(i,j+1,k).isCovered() || !w_cflag(i,j,k).isCovered()) {
183  s23bar += tau23(i, j+1, k);
184  count_s23 += 1.0;
185  }
186  if (!v_cflag(i,j+1,k).isCovered() || !w_cflag(i,j,k+1).isCovered()) {
187  s23bar += tau23(i, j+1, k+1);
188  count_s23 += 1.0;
189  }
190  if (count_s23 > 0.0) {
191  s23bar /= count_s23;
192  }
193 
194  SmnSmn = s11bar*s11bar + s22bar*s22bar + s33bar*s33bar
195  + 2.0*s12bar*s12bar + 2.0*s13bar*s13bar + 2.0*s23bar*s23bar;
196  }
197 
198  return SmnSmn;
199 }

Referenced by ComputeTurbulentViscosityLES_EB().

Here is the caller graph for this function:

◆ ComputeTurbulentViscosity()

void ComputeTurbulentViscosity ( amrex::Real  dt,
const amrex::MultiFab &  xvel,
const amrex::MultiFab &  yvel,
amrex::Vector< std::unique_ptr< amrex::MultiFab >> &  Tau_lev,
amrex::MultiFab &  cons_in,
const amrex::MultiFab &  wdist,
amrex::MultiFab &  eddyViscosity,
amrex::MultiFab &  Hfx1,
amrex::MultiFab &  Hfx2,
amrex::MultiFab &  Hfx3,
amrex::MultiFab &  Diss,
const amrex::Geometry &  geom,
amrex::Vector< std::unique_ptr< amrex::MultiFab >> &  mapfac,
const std::unique_ptr< amrex::MultiFab > &  z_phys_nd,
const SolverChoice &  solverChoice,
std::unique_ptr< SurfaceLayer > &  SurfLayer,
const amrex::MultiFab *  z_0,
const bool &  use_terrain_fitted_coords,
const bool &  use_moisture,
int  level,
const amrex::BCRec *  bc_ptr,
const eb_ &  ebfact,
bool  vert_only = false 
)

◆ ComputeTurbulentViscosityLES()

void ComputeTurbulentViscosityLES ( amrex::Vector< std::unique_ptr< amrex::MultiFab >> &  Tau_lev,
const amrex::MultiFab &  cons_in,
amrex::MultiFab &  eddyViscosity,
amrex::MultiFab &  Hfx1,
amrex::MultiFab &  Hfx2,
amrex::MultiFab &  Hfx3,
amrex::MultiFab &  Diss,
const amrex::Geometry &  geom,
bool  use_terrain_fitted_coords,
amrex::Vector< std::unique_ptr< amrex::MultiFab >> &  mapfac,
const std::unique_ptr< amrex::MultiFab > &  z_phys_nd,
const TurbChoice &  turbChoice,
const amrex::Real  const_grav,
std::unique_ptr< SurfaceLayer > &  SurfLayer,
const MoistureComponentIndices &  moisture_indices,
const amrex::MultiFab *  xvel = nullptr,
const amrex::MultiFab *  yvel = nullptr 
)

◆ ComputeTurbulentViscosityLES_EB()

void ComputeTurbulentViscosityLES_EB ( amrex::Vector< std::unique_ptr< amrex::MultiFab >> &  Tau_lev,
const amrex::MultiFab &  cons_in,
amrex::MultiFab &  eddyViscosity,
amrex::MultiFab &  Hfx1,
amrex::MultiFab &  Hfx2,
amrex::MultiFab &  Hfx3,
const amrex::Geometry &  geom,
amrex::Vector< std::unique_ptr< amrex::MultiFab >> &  mapfac,
const TurbChoice &  turbChoice,
const amrex::Real  const_grav,
const SolverChoice &  solverChoice,
std::unique_ptr< SurfaceLayer > &  SurfLayer,
const MoistureComponentIndices &  moisture_indices,
const eb_ &  ebfact,
const amrex::MultiFab *  xvel = nullptr,
const amrex::MultiFab *  yvel = nullptr 
)