3 #ifndef ERF_EDDY_VISCOSITY_H_
4 #define ERF_EDDY_VISCOSITY_H_
6 #include "AMReX_BCRec.H"
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,
29 std::unique_ptr<SurfaceLayer>& SurfLayer,
30 const amrex::MultiFab* z_0,
31 const bool& use_terrain_fitted_coords,
32 const bool& use_moisture,
34 const amrex::BCRec* bc_ptr,
35 bool vert_only =
false);
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)
58 amrex::Real SmnSmn = s11bar*s11bar + s22bar*s22bar + s33bar*s33bar
59 + 2.0*s12bar*s12bar + 2.0*s13bar*s13bar + 2.0*s23bar*s23bar;
68 const amrex::Array4<amrex::Real const>&
tau11,
69 const amrex::Array4<amrex::Real const>&
tau22,
70 const amrex::Array4<amrex::Real const>&
tau12)
75 return 0.5 * sdiff*sdiff + 2.0*s12bar*s12bar;
@ 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