ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
ERF_EddyViscosity.H
Go to the documentation of this file.
1 /** \file ERF_EddyViscosity.H */
2 
3 #ifndef ERF_EDDY_VISCOSITY_H_
4 #define ERF_EDDY_VISCOSITY_H_
5 
6 #include "AMReX_BCRec.H"
7 #include "ERF_SurfaceLayer.H"
8 #include "ERF_DataStruct.H"
9 #include "ERF_IndexDefines.H"
10 #include "ERF_Constants.H"
11 #include "ERF_EOS.H"
12 
13 void
15  const amrex::MultiFab& xvel ,
16  const amrex::MultiFab& yvel ,
17  amrex::Vector<std::unique_ptr<amrex::MultiFab>>& Tau_lev,
18  amrex::MultiFab& cons_in,
19  const amrex::MultiFab& wdist,
20  amrex::MultiFab& eddyViscosity,
21  amrex::MultiFab& Hfx1,
22  amrex::MultiFab& Hfx2,
23  amrex::MultiFab& Hfx3,
24  amrex::MultiFab& Diss,
25  const amrex::Geometry& geom,
26  amrex::Vector<std::unique_ptr<amrex::MultiFab>>& mapfac,
27  const std::unique_ptr<amrex::MultiFab>& z_phys_nd,
28  const SolverChoice& solverChoice,
29  std::unique_ptr<SurfaceLayer>& SurfLayer,
30  const amrex::MultiFab* z_0,
31  const bool& use_terrain_fitted_coords,
32  const bool& use_moisture,
33  int level,
34  const amrex::BCRec* bc_ptr,
35  bool vert_only = false);
36 
37 AMREX_GPU_DEVICE
38 AMREX_FORCE_INLINE
40 ComputeSmnSmn (int& i, int& j, int& k,
41  const amrex::Array4<amrex::Real const>& tau11,
42  const amrex::Array4<amrex::Real const>& tau22,
43  const amrex::Array4<amrex::Real const>& tau33,
44  const amrex::Array4<amrex::Real const>& tau12,
45  const amrex::Array4<amrex::Real const>& tau13,
46  const amrex::Array4<amrex::Real const>& tau23)
47 {
48  amrex::Real s11bar = tau11(i,j,k);
49  amrex::Real s22bar = tau22(i,j,k);
50  amrex::Real s33bar = tau33(i,j,k);
51  amrex::Real s12bar = 0.25 * ( tau12(i , j , k ) + tau12(i , j+1, k )
52  + tau12(i+1, j , k ) + tau12(i+1, j+1, k ) );
53  amrex::Real s13bar = 0.25 * ( tau13(i , j , k ) + tau13(i , j , k+1)
54  + tau13(i+1, j , k ) + tau13(i+1, j , k+1) );
55  amrex::Real s23bar = 0.25 * ( tau23(i , j , k ) + tau23(i , j , k+1)
56  + tau23(i , j+1, k ) + tau23(i , j+1, k+1) );
57 
58  amrex::Real SmnSmn = s11bar*s11bar + s22bar*s22bar + s33bar*s33bar
59  + 2.0*s12bar*s12bar + 2.0*s13bar*s13bar + 2.0*s23bar*s23bar;
60 
61  return SmnSmn;
62 }
63 
64 AMREX_GPU_DEVICE
65 AMREX_FORCE_INLINE
67 ComputeSmnSmn2D (int& i, int& j, int& k,
68  const amrex::Array4<amrex::Real const>& tau11,
69  const amrex::Array4<amrex::Real const>& tau22,
70  const amrex::Array4<amrex::Real const>& tau12)
71 {
72  amrex::Real sdiff = tau11(i,j,k) - tau22(i,j,k);
73  amrex::Real s12bar = 0.25 * ( tau12(i , j , k ) + tau12(i , j+1, k )
74  + tau12(i+1, j , k ) + tau12(i+1, j+1, k ) );
75  return 0.5 * sdiff*sdiff + 2.0*s12bar*s12bar;
76 }
77 #endif
@ tau12
Definition: ERF_DataStruct.H:30
@ tau23
Definition: ERF_DataStruct.H:30
@ tau33
Definition: ERF_DataStruct.H:30
@ tau22
Definition: ERF_DataStruct.H:30
@ tau11
Definition: ERF_DataStruct.H:30
@ tau13
Definition: ERF_DataStruct.H:30
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)
Definition: ERF_EddyViscosity.H:67
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, bool vert_only=false)
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)
Definition: ERF_EddyViscosity.H:40
amrex::Real Real
Definition: ERF_ShocInterface.H:19
@ xvel
Definition: ERF_IndexDefines.H:141
@ yvel
Definition: ERF_IndexDefines.H:142
Definition: ERF_DataStruct.H:123