1 #ifndef ERF_DIFFUSION_H_
2 #define ERF_DIFFUSION_H_
5 #include "AMReX_MultiFab.H"
6 #include "AMReX_BCRec.H"
14 const amrex::Array4< amrex::Real>& rho_u_rhs,
15 const amrex::Array4< amrex::Real>& rho_v_rhs,
16 const amrex::Array4< amrex::Real>& rho_w_rhs,
17 const amrex::Array4<const amrex::Real>&
tau11 ,
18 const amrex::Array4<const amrex::Real>&
tau22 ,
19 const amrex::Array4<const amrex::Real>&
tau33 ,
20 const amrex::Array4<const amrex::Real>&
tau12 ,
const amrex::Array4<const amrex::Real>&
tau21,
21 const amrex::Array4<const amrex::Real>&
tau13 ,
const amrex::Array4<const amrex::Real>&
tau31,
22 const amrex::Array4<const amrex::Real>&
tau23 ,
const amrex::Array4<const amrex::Real>&
tau32,
23 const amrex::Array4<const amrex::Real>& detJ_arr ,
24 const amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
25 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
26 const amrex::Array4<const amrex::Real>& mf_mx ,
27 const amrex::Array4<const amrex::Real>& mf_ux ,
28 const amrex::Array4<const amrex::Real>& mf_vx ,
29 const amrex::Array4<const amrex::Real>& mf_my ,
30 const amrex::Array4<const amrex::Real>& mf_uy ,
31 const amrex::Array4<const amrex::Real>& mf_vy ,
32 const bool use_stretched_dz,
const bool use_variable_dz);
35 [[maybe_unused]]
const amrex::Box& domain,
36 const amrex::Box& bxx,
const amrex::Box& bxy,
const amrex::Box& bxz,
37 const amrex::Array4< amrex::Real>& rho_u_rhs,
38 const amrex::Array4< amrex::Real>& rho_v_rhs,
39 const amrex::Array4< amrex::Real>& rho_w_rhs,
40 const amrex::Array4<const amrex::Real>& u_arr ,
41 const amrex::Array4<const amrex::Real>& v_arr ,
42 const amrex::Array4<const amrex::Real>& w_arr ,
43 const amrex::Array4<const amrex::Real>&
tau11 ,
44 const amrex::Array4<const amrex::Real>&
tau22 ,
45 const amrex::Array4<const amrex::Real>&
tau33 ,
46 const amrex::Array4<const amrex::Real>&
tau12 ,
47 const amrex::Array4<const amrex::Real>&
tau13 ,
48 const amrex::Array4<const amrex::Real>&
tau23 ,
50 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
51 const amrex::Array4<const amrex::Real>& mf_mx ,
52 const amrex::Array4<const amrex::Real>& mf_ux ,
53 const amrex::Array4<const amrex::Real>& mf_vx ,
54 const amrex::Array4<const amrex::Real>& mf_my ,
55 const amrex::Array4<const amrex::Real>& mf_uy ,
56 const amrex::Array4<const amrex::Real>& mf_vy ,
59 [[maybe_unused]]
const amrex::BCRec* bc_ptr);
62 int start_comp,
int num_comp,
63 const amrex::Array4<const amrex::Real>& u,
64 const amrex::Array4<const amrex::Real>& v,
65 const amrex::Array4<const amrex::Real>& cell_data,
66 const amrex::Array4<const amrex::Real>& cell_prim,
67 const amrex::Array4<amrex::Real>& cell_rhs,
68 const amrex::Array4<amrex::Real>& xflux,
69 const amrex::Array4<amrex::Real>& yflux,
70 const amrex::Array4<amrex::Real>& zflux,
71 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
72 const amrex::Array4<const amrex::Real>& SmnSmn_a,
73 const amrex::Array4<const amrex::Real>& mf_mx,
74 const amrex::Array4<const amrex::Real>& mf_ux,
75 const amrex::Array4<const amrex::Real>& mf_vx,
76 const amrex::Array4<const amrex::Real>& mf_my,
77 const amrex::Array4<const amrex::Real>& mf_uy,
78 const amrex::Array4<const amrex::Real>& mf_vy,
79 amrex::Array4< amrex::Real>& hfx_z,
80 amrex::Array4< amrex::Real>& qfx1_z,
81 amrex::Array4< amrex::Real>& qfx2_z,
82 amrex::Array4< amrex::Real>& diss,
83 const amrex::Array4<const amrex::Real>& mu_turb,
86 const amrex::Array4<const amrex::Real>& tm_arr,
87 const amrex::GpuArray<amrex::Real,AMREX_SPACEDIM> grav_gpu,
88 const amrex::BCRec* bc_ptr,
89 const bool use_SurfLayer,
93 int start_comp,
int num_comp,
94 const amrex::Array4<const amrex::Real>& u,
95 const amrex::Array4<const amrex::Real>& v,
96 const amrex::Array4<const amrex::Real>& cell_data,
97 const amrex::Array4<const amrex::Real>& cell_prim,
98 const amrex::Array4<amrex::Real>& cell_rhs,
99 const amrex::Array4<amrex::Real>& xflux,
100 const amrex::Array4<amrex::Real>& yflux,
101 const amrex::Array4<amrex::Real>& zflux,
102 const amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
103 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
104 const amrex::Array4<const amrex::Real>& SmnSmn_a,
105 const amrex::Array4<const amrex::Real>& mf_mx,
106 const amrex::Array4<const amrex::Real>& mf_ux,
107 const amrex::Array4<const amrex::Real>& mf_vx,
108 const amrex::Array4<const amrex::Real>& mf_my,
109 const amrex::Array4<const amrex::Real>& mf_uy,
110 const amrex::Array4<const amrex::Real>& mf_vy,
111 amrex::Array4< amrex::Real>& hfx_z,
112 amrex::Array4< amrex::Real>& qfx1_z,
113 amrex::Array4< amrex::Real>& qfx2_z,
114 amrex::Array4< amrex::Real>& diss,
115 const amrex::Array4<const amrex::Real>& mu_turb,
118 const amrex::Array4<const amrex::Real>& tm_arr,
119 const amrex::GpuArray<amrex::Real,AMREX_SPACEDIM> grav_gpu,
120 const amrex::BCRec* bc_ptr,
121 const bool use_SurfLayer,
125 int start_comp,
int num_comp,
127 const amrex::Array4<const amrex::Real>& u,
128 const amrex::Array4<const amrex::Real>& v,
129 const amrex::Array4<const amrex::Real>& cell_data,
130 const amrex::Array4<const amrex::Real>& cell_prim,
131 const amrex::Array4<amrex::Real>& cell_rhs,
132 const amrex::Array4<amrex::Real>& xflux,
133 const amrex::Array4<amrex::Real>& yflux,
134 const amrex::Array4<amrex::Real>& zflux,
135 const amrex::Array4<const amrex::Real>& z_nd,
136 const amrex::Array4<const amrex::Real>& z_cc,
137 const amrex::Array4<const amrex::Real>& ax,
138 const amrex::Array4<const amrex::Real>& ay,
139 const amrex::Array4<const amrex::Real>& az,
140 const amrex::Array4<const amrex::Real>& detJ,
141 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
142 const amrex::Array4<const amrex::Real>& SmnSmn_a,
143 const amrex::Array4<const amrex::Real>& mf_mx,
144 const amrex::Array4<const amrex::Real>& mf_ux,
145 const amrex::Array4<const amrex::Real>& mf_vx,
146 const amrex::Array4<const amrex::Real>& mf_my,
147 const amrex::Array4<const amrex::Real>& mf_uy,
148 const amrex::Array4<const amrex::Real>& mf_vy,
149 amrex::Array4< amrex::Real>& hfx_x,
150 amrex::Array4< amrex::Real>& hfx_y,
151 amrex::Array4< amrex::Real>& hfx_z,
152 amrex::Array4< amrex::Real>& qfx1_x,
153 amrex::Array4< amrex::Real>& qfx1_y,
154 amrex::Array4< amrex::Real>& qfx1_z,
155 amrex::Array4< amrex::Real>& qfx2_z,
156 amrex::Array4< amrex::Real>& diss,
157 const amrex::Array4<const amrex::Real>& mu_turb,
160 const amrex::Array4<const amrex::Real>& tm_arr,
161 const amrex::GpuArray<amrex::Real,AMREX_SPACEDIM> grav_gpu,
162 const amrex::BCRec* bc_ptr,
163 const bool use_SurfLayer,
172 const amrex::Array4<const amrex::Real>& cell_data,
173 amrex::Array4<amrex::Real>&
tau11,
174 amrex::Array4<amrex::Real>&
tau22,
175 amrex::Array4<amrex::Real>&
tau33,
176 amrex::Array4<amrex::Real>&
tau12,
177 amrex::Array4<amrex::Real>&
tau13,
178 amrex::Array4<amrex::Real>&
tau23,
179 const amrex::Array4<const amrex::Real>& er_arr,
180 amrex::Array4<amrex::Real>& tau13i,
181 amrex::Array4<amrex::Real>& tau23i,
182 amrex::Array4<amrex::Real>& tau33i);
189 const amrex::Array4<const amrex::Real>& cell_data,
190 amrex::Array4<amrex::Real>&
tau11,
191 amrex::Array4<amrex::Real>&
tau22,
192 amrex::Array4<amrex::Real>&
tau33,
193 amrex::Array4<amrex::Real>&
tau12, amrex::Array4<amrex::Real>&
tau21,
194 amrex::Array4<amrex::Real>&
tau13, amrex::Array4<amrex::Real>&
tau31,
195 amrex::Array4<amrex::Real>&
tau23, amrex::Array4<amrex::Real>&
tau32,
196 const amrex::Array4<const amrex::Real>& er_arr,
197 const amrex::Array4<const amrex::Real>& mf_mx,
198 const amrex::Array4<const amrex::Real>& mf_ux,
199 const amrex::Array4<const amrex::Real>& mf_vx,
200 const amrex::Array4<const amrex::Real>& mf_my,
201 const amrex::Array4<const amrex::Real>& mf_uy,
202 const amrex::Array4<const amrex::Real>& mf_vy,
203 amrex::Array4<amrex::Real>& tau13i,
204 amrex::Array4<amrex::Real>& tau23i,
205 amrex::Array4<amrex::Real>& tau33i);
212 const amrex::Array4<const amrex::Real>& cell_data,
213 amrex::Array4<amrex::Real>&
tau11,
214 amrex::Array4<amrex::Real>&
tau22,
215 amrex::Array4<amrex::Real>&
tau33,
216 amrex::Array4<amrex::Real>&
tau12, amrex::Array4<amrex::Real>&
tau21,
217 amrex::Array4<amrex::Real>&
tau13, amrex::Array4<amrex::Real>&
tau31,
218 amrex::Array4<amrex::Real>&
tau23, amrex::Array4<amrex::Real>&
tau32,
219 const amrex::Array4<const amrex::Real>& er_arr,
220 const amrex::Array4<const amrex::Real>& z_nd,
221 const amrex::Array4<const amrex::Real>& detJ,
222 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
223 const amrex::Array4<const amrex::Real>& mf_mx,
224 const amrex::Array4<const amrex::Real>& mf_ux,
225 const amrex::Array4<const amrex::Real>& mf_vx,
226 const amrex::Array4<const amrex::Real>& mf_my,
227 const amrex::Array4<const amrex::Real>& mf_uy,
228 const amrex::Array4<const amrex::Real>& mf_vy,
229 amrex::Array4<amrex::Real>& tau13i,
230 amrex::Array4<amrex::Real>& tau23i,
231 amrex::Array4<amrex::Real>& tau33i);
238 const amrex::Array4<const amrex::Real>& mu_turb,
239 const amrex::Array4<const amrex::Real>& cell_data,
240 amrex::Array4<amrex::Real>&
tau11,
241 amrex::Array4<amrex::Real>&
tau22,
242 amrex::Array4<amrex::Real>&
tau33,
243 amrex::Array4<amrex::Real>&
tau12,
244 amrex::Array4<amrex::Real>&
tau13,
245 amrex::Array4<amrex::Real>&
tau23,
246 const amrex::Array4<const amrex::Real>& er_arr,
247 amrex::Array4<amrex::Real>& tau13i,
248 amrex::Array4<amrex::Real>& tau23i,
249 amrex::Array4<amrex::Real>& tau33i);
256 const amrex::Array4<const amrex::Real>& mu_turb,
257 const amrex::Array4<const amrex::Real>& cell_data,
258 amrex::Array4<amrex::Real>&
tau11,
259 amrex::Array4<amrex::Real>&
tau22,
260 amrex::Array4<amrex::Real>&
tau33,
261 amrex::Array4<amrex::Real>&
tau12, amrex::Array4<amrex::Real>&
tau21,
262 amrex::Array4<amrex::Real>&
tau13, amrex::Array4<amrex::Real>&
tau31,
263 amrex::Array4<amrex::Real>&
tau23, amrex::Array4<amrex::Real>&
tau32,
264 const amrex::Array4<const amrex::Real>& er_arr,
265 const amrex::Array4<const amrex::Real>& mf_mx,
266 const amrex::Array4<const amrex::Real>& mf_ux,
267 const amrex::Array4<const amrex::Real>& mf_vx,
268 const amrex::Array4<const amrex::Real>& mf_my,
269 const amrex::Array4<const amrex::Real>& mf_uy,
270 const amrex::Array4<const amrex::Real>& mf_vy,
271 amrex::Array4<amrex::Real>& tau13i,
272 amrex::Array4<amrex::Real>& tau23i,
273 amrex::Array4<amrex::Real>& tau33i);
280 const amrex::Array4<const amrex::Real>& mu_turb,
281 const amrex::Array4<const amrex::Real>& cell_data,
282 amrex::Array4<amrex::Real>&
tau11,
283 amrex::Array4<amrex::Real>&
tau22,
284 amrex::Array4<amrex::Real>&
tau33,
285 amrex::Array4<amrex::Real>&
tau12, amrex::Array4<amrex::Real>&
tau21,
286 amrex::Array4<amrex::Real>&
tau13, amrex::Array4<amrex::Real>&
tau31,
287 amrex::Array4<amrex::Real>&
tau23, amrex::Array4<amrex::Real>&
tau32,
288 const amrex::Array4<const amrex::Real>& er_arr,
289 const amrex::Array4<const amrex::Real>& z_nd,
290 const amrex::Array4<const amrex::Real>& detJ,
291 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
292 const amrex::Array4<const amrex::Real>& mf_mx,
293 const amrex::Array4<const amrex::Real>& mf_ux,
294 const amrex::Array4<const amrex::Real>& mf_vx,
295 const amrex::Array4<const amrex::Real>& mf_my,
296 const amrex::Array4<const amrex::Real>& mf_uy,
297 const amrex::Array4<const amrex::Real>& mf_vy,
298 amrex::Array4<amrex::Real>& tau13i,
299 amrex::Array4<amrex::Real>& tau23i,
300 amrex::Array4<amrex::Real>& tau33i);
305 amrex::Box tbxxz, amrex::Box tbxyz, amrex::Box domain,
306 const amrex::Array4<const amrex::Real>& u,
307 const amrex::Array4<const amrex::Real>& v,
308 const amrex::Array4<const amrex::Real>& w,
309 amrex::Array4<amrex::Real>&
tau11,
310 amrex::Array4<amrex::Real>&
tau22,
311 amrex::Array4<amrex::Real>&
tau33,
312 amrex::Array4<amrex::Real>&
tau12,
313 amrex::Array4<amrex::Real>&
tau13,
314 amrex::Array4<amrex::Real>&
tau23,
315 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
316 const amrex::Array4<const amrex::Real>& mf_mx,
317 const amrex::Array4<const amrex::Real>& mf_ux,
318 const amrex::Array4<const amrex::Real>& mf_vx,
319 const amrex::Array4<const amrex::Real>& mf_my,
320 const amrex::Array4<const amrex::Real>& mf_uy,
321 const amrex::Array4<const amrex::Real>& mf_vy,
322 const amrex::BCRec* bc_ptr,
323 amrex::Array4<amrex::Real>& tau13i,
324 amrex::Array4<amrex::Real>& tau23i);
327 amrex::Box tbxxz, amrex::Box tbxyz, amrex::Box domain,
328 const amrex::Array4<const amrex::Real>& u,
329 const amrex::Array4<const amrex::Real>& v,
330 const amrex::Array4<const amrex::Real>& w,
331 amrex::Array4<amrex::Real>&
tau11,
332 amrex::Array4<amrex::Real>&
tau22,
333 amrex::Array4<amrex::Real>&
tau33,
334 amrex::Array4<amrex::Real>&
tau12, amrex::Array4<amrex::Real>&
tau21,
335 amrex::Array4<amrex::Real>&
tau13, amrex::Array4<amrex::Real>&
tau31,
336 amrex::Array4<amrex::Real>&
tau23, amrex::Array4<amrex::Real>&
tau32,
337 const amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
338 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
339 const amrex::Array4<const amrex::Real>& mf_mx,
340 const amrex::Array4<const amrex::Real>& mf_ux,
341 const amrex::Array4<const amrex::Real>& mf_vx,
342 const amrex::Array4<const amrex::Real>& mf_my,
343 const amrex::Array4<const amrex::Real>& mf_uy,
344 const amrex::Array4<const amrex::Real>& mf_vy,
345 const amrex::BCRec* bc_ptr,
346 amrex::Array4<amrex::Real>& tau13i,
347 amrex::Array4<amrex::Real>& tau23i);
350 amrex::Box tbxxz, amrex::Box tbxyz, amrex::Box domain,
351 const amrex::Array4<const amrex::Real>& u,
352 const amrex::Array4<const amrex::Real>& v,
353 const amrex::Array4<const amrex::Real>& w,
354 amrex::Array4<amrex::Real>&
tau11,
355 amrex::Array4<amrex::Real>&
tau22,
356 amrex::Array4<amrex::Real>&
tau33,
357 amrex::Array4<amrex::Real>&
tau12, amrex::Array4<amrex::Real>&
tau21,
358 amrex::Array4<amrex::Real>&
tau13, amrex::Array4<amrex::Real>&
tau31,
359 amrex::Array4<amrex::Real>&
tau23, amrex::Array4<amrex::Real>&
tau32,
360 const amrex::Array4<const amrex::Real>& z_nd,
361 const amrex::Array4<const amrex::Real>& detJ,
362 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
363 const amrex::Array4<const amrex::Real>& mf_mx,
364 const amrex::Array4<const amrex::Real>& mf_ux,
365 const amrex::Array4<const amrex::Real>& mf_vx,
366 const amrex::Array4<const amrex::Real>& mf_my,
367 const amrex::Array4<const amrex::Real>& mf_uy,
368 const amrex::Array4<const amrex::Real>& mf_vy,
369 const amrex::BCRec* bc_ptr,
370 amrex::Array4<amrex::Real>& tau13i,
371 amrex::Array4<amrex::Real>& tau23i);
376 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM*2>& bc_neumann_vals,
377 const amrex::Array4< amrex::Real>& cell_data,
378 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
379 const amrex::Array4<const amrex::Real>& hfx_z,
380 const amrex::Array4<const amrex::Real>& mu_turb,
382 const amrex::BCRec* bc_ptr,
383 const bool use_SurfLayer,
389 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM*2>& bc_neumann_vals,
390 const amrex::Array4< amrex::Real>& cell_data,
391 const amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
392 const amrex::Array4<const amrex::Real>& hfx_z,
393 const amrex::Array4<const amrex::Real>& mu_turb,
395 const amrex::BCRec* bc_ptr,
396 const bool use_SurfLayer,
402 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM*2>& bc_neumann_vals,
403 const amrex::Array4< amrex::Real>& cell_data,
404 const amrex::Array4<const amrex::Real>& z_nd,
405 const amrex::Array4<const amrex::Real>& detJ,
406 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
407 const amrex::Array4<const amrex::Real>& hfx_z,
408 const amrex::Array4<const amrex::Real>& mu_turb,
410 const amrex::BCRec* bc_ptr,
411 const bool use_SurfLayer,
414 template <
int stagdir>
417 const amrex::Array4<const amrex::Real>& cell_data,
418 const amrex::Array4< amrex::Real>& face_data,
419 const amrex::Array4<const amrex::Real>& tau_corr,
420 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
421 const amrex::Array4<const amrex::Real>& mu_turb,
423 const amrex::BCRec* bc_ptr,
424 const bool use_SurfLayer,
427 template <
int stagdir>
430 const amrex::Array4<const amrex::Real>& cell_data,
431 const amrex::Array4< amrex::Real>& face_data,
432 const amrex::Array4<const amrex::Real>& tau_corr,
433 const amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
434 const amrex::Array4<const amrex::Real>& mu_turb,
436 const amrex::BCRec* bc_ptr,
437 const bool use_SurfLayer,
440 template <
int stagdir>
443 const amrex::Array4<const amrex::Real>& cell_data,
444 const amrex::Array4< amrex::Real>& face_data,
445 const amrex::Array4<const amrex::Real>& tau_corr,
446 const amrex::Array4<const amrex::Real>& z_nd,
447 const amrex::Array4<const amrex::Real>& detJ,
448 const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
449 const amrex::Array4<const amrex::Real>& mu_turb,
451 const amrex::BCRec* bc_ptr,
452 const bool use_SurfLayer,
@ tau12
Definition: ERF_DataStruct.H:31
@ tau23
Definition: ERF_DataStruct.H:31
@ tau33
Definition: ERF_DataStruct.H:31
@ tau22
Definition: ERF_DataStruct.H:31
@ tau11
Definition: ERF_DataStruct.H:31
@ tau32
Definition: ERF_DataStruct.H:31
@ tau31
Definition: ERF_DataStruct.H:31
@ tau21
Definition: ERF_DataStruct.H:31
@ tau13
Definition: ERF_DataStruct.H:31
void ComputeStressVarVisc_T(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Real mu_eff, const amrex::Array4< const amrex::Real > &mu_turb, const amrex::Array4< const amrex::Real > &cell_data, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau21, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau31, amrex::Array4< amrex::Real > &tau23, amrex::Array4< amrex::Real > &tau32, const amrex::Array4< const amrex::Real > &er_arr, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i, amrex::Array4< amrex::Real > &tau33i)
void ImplicitDiffForMom_T(const amrex::Box &bx, const amrex::Box &domain, const int level, const amrex::Real dt, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< amrex::Real > &face_data, const amrex::Array4< const amrex::Real > &tau_corr, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void DiffusionSrcForMom(const amrex::Box &bxx, const amrex::Box &bxy, const amrex::Box &bxz, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &tau11, const amrex::Array4< const amrex::Real > &tau22, const amrex::Array4< const amrex::Real > &tau33, const amrex::Array4< const amrex::Real > &tau12, const amrex::Array4< const amrex::Real > &tau21, const amrex::Array4< const amrex::Real > &tau13, const amrex::Array4< const amrex::Real > &tau31, const amrex::Array4< const amrex::Real > &tau23, const amrex::Array4< const amrex::Real > &tau32, const amrex::Array4< const amrex::Real > &detJ_arr, const amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const bool use_stretched_dz, const bool use_variable_dz)
void ComputeStressConsVisc_N(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Real mu_eff, const amrex::Array4< const amrex::Real > &cell_data, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau23, const amrex::Array4< const amrex::Real > &er_arr, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i, amrex::Array4< amrex::Real > &tau33i)
void DiffusionSrcForMom_EB(const amrex::MFIter &mfi, [[maybe_unused]] const amrex::Box &domain, const amrex::Box &bxx, const amrex::Box &bxy, const amrex::Box &bxz, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &u_arr, const amrex::Array4< const amrex::Real > &v_arr, const amrex::Array4< const amrex::Real > &w_arr, const amrex::Array4< const amrex::Real > &tau11, const amrex::Array4< const amrex::Real > &tau22, const amrex::Array4< const amrex::Real > &tau33, const amrex::Array4< const amrex::Real > &tau12, const amrex::Array4< const amrex::Real > &tau13, const amrex::Array4< const amrex::Real > &tau23, const amrex::Real *dx_arr, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const SolverChoice &solverChoice, const eb_ &ebfact, [[maybe_unused]] const amrex::BCRec *bc_ptr)
void DiffusionSrcForState_S(const amrex::Box &bx, const amrex::Box &domain, int start_comp, int num_comp, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &cell_prim, const amrex::Array4< amrex::Real > &cell_rhs, const amrex::Array4< amrex::Real > &xflux, const amrex::Array4< amrex::Real > &yflux, const amrex::Array4< amrex::Real > &zflux, const amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &SmnSmn_a, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, amrex::Array4< amrex::Real > &hfx_z, amrex::Array4< amrex::Real > &qfx1_z, amrex::Array4< amrex::Real > &qfx2_z, amrex::Array4< amrex::Real > &diss, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const int level, const amrex::Array4< const amrex::Real > &tm_arr, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > grav_gpu, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void DiffusionSrcForState_T(const amrex::Box &bx, const amrex::Box &domain, int start_comp, int num_comp, const bool &rotate, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &cell_prim, const amrex::Array4< amrex::Real > &cell_rhs, const amrex::Array4< amrex::Real > &xflux, const amrex::Array4< amrex::Real > &yflux, const amrex::Array4< amrex::Real > &zflux, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &z_cc, const amrex::Array4< const amrex::Real > &ax, const amrex::Array4< const amrex::Real > &ay, const amrex::Array4< const amrex::Real > &az, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &SmnSmn_a, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, amrex::Array4< amrex::Real > &hfx_x, amrex::Array4< amrex::Real > &hfx_y, amrex::Array4< amrex::Real > &hfx_z, amrex::Array4< amrex::Real > &qfx1_x, amrex::Array4< amrex::Real > &qfx1_y, amrex::Array4< amrex::Real > &qfx1_z, amrex::Array4< amrex::Real > &qfx2_z, amrex::Array4< amrex::Real > &diss, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const int level, const amrex::Array4< const amrex::Real > &tm_arr, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > grav_gpu, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void ComputeStrain_T(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Box domain, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau21, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau31, amrex::Array4< amrex::Real > &tau23, amrex::Array4< amrex::Real > &tau32, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const amrex::BCRec *bc_ptr, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i)
void ComputeStressConsVisc_T(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Real mu_eff, const amrex::Array4< const amrex::Real > &cell_data, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau21, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau31, amrex::Array4< amrex::Real > &tau23, amrex::Array4< amrex::Real > &tau32, const amrex::Array4< const amrex::Real > &er_arr, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i, amrex::Array4< amrex::Real > &tau33i)
void ImplicitDiffForMom_N(const amrex::Box &bx, const amrex::Box &domain, const int level, const amrex::Real dt, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< amrex::Real > &face_data, const amrex::Array4< const amrex::Real > &tau_corr, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void ImplicitDiffForMom_S(const amrex::Box &bx, const amrex::Box &domain, const int level, const amrex::Real dt, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< amrex::Real > &face_data, const amrex::Array4< const amrex::Real > &tau_corr, const amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void DiffusionSrcForState_N(const amrex::Box &bx, const amrex::Box &domain, int start_comp, int num_comp, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &cell_prim, const amrex::Array4< amrex::Real > &cell_rhs, const amrex::Array4< amrex::Real > &xflux, const amrex::Array4< amrex::Real > &yflux, const amrex::Array4< amrex::Real > &zflux, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &SmnSmn_a, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, amrex::Array4< amrex::Real > &hfx_z, amrex::Array4< amrex::Real > &qfx1_z, amrex::Array4< amrex::Real > &qfx2_z, amrex::Array4< amrex::Real > &diss, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const int level, const amrex::Array4< const amrex::Real > &tm_arr, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > grav_gpu, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void ComputeStressConsVisc_S(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Real mu_eff, const amrex::Array4< const amrex::Real > &cell_data, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau21, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau31, amrex::Array4< amrex::Real > &tau23, amrex::Array4< amrex::Real > &tau32, const amrex::Array4< const amrex::Real > &er_arr, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i, amrex::Array4< amrex::Real > &tau33i)
void ComputeStrain_N(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Box domain, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau23, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const amrex::BCRec *bc_ptr, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i)
void ComputeStressVarVisc_S(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Real mu_eff, const amrex::Array4< const amrex::Real > &mu_turb, const amrex::Array4< const amrex::Real > &cell_data, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau21, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau31, amrex::Array4< amrex::Real > &tau23, amrex::Array4< amrex::Real > &tau32, const amrex::Array4< const amrex::Real > &er_arr, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i, amrex::Array4< amrex::Real > &tau33i)
void ImplicitDiffForState_N(const amrex::Box &bx, const amrex::Box &domain, const int level, const amrex::Real dt, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM *2 > &bc_neumann_vals, const amrex::Array4< amrex::Real > &cell_data, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &hfx_z, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void ImplicitDiffForState_S(const amrex::Box &bx, const amrex::Box &domain, const int level, const amrex::Real dt, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM *2 > &bc_neumann_vals, const amrex::Array4< amrex::Real > &cell_data, const amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::Array4< const amrex::Real > &hfx_z, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void ImplicitDiffForState_T(const amrex::Box &bx, const amrex::Box &domain, const int level, const amrex::Real dt, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM *2 > &bc_neumann_vals, const amrex::Array4< amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &hfx_z, const amrex::Array4< const amrex::Real > &mu_turb, const SolverChoice &solverChoice, const amrex::BCRec *bc_ptr, const bool use_SurfLayer, const amrex::Real implicit_fac)
void ComputeStressVarVisc_N(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Real mu_eff, const amrex::Array4< const amrex::Real > &mu_turb, const amrex::Array4< const amrex::Real > &cell_data, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau23, const amrex::Array4< const amrex::Real > &er_arr, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i, amrex::Array4< amrex::Real > &tau33i)
void ComputeStrain_S(amrex::Box bxcc, amrex::Box tbxxy, amrex::Box tbxxz, amrex::Box tbxyz, amrex::Box domain, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, amrex::Array4< amrex::Real > &tau11, amrex::Array4< amrex::Real > &tau22, amrex::Array4< amrex::Real > &tau33, amrex::Array4< amrex::Real > &tau12, amrex::Array4< amrex::Real > &tau21, amrex::Array4< amrex::Real > &tau13, amrex::Array4< amrex::Real > &tau31, amrex::Array4< amrex::Real > &tau23, amrex::Array4< amrex::Real > &tau32, const amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const amrex::BCRec *bc_ptr, amrex::Array4< amrex::Real > &tau13i, amrex::Array4< amrex::Real > &tau23i)
amrex::Real Real
Definition: ERF_ShocInterface.H:19
Definition: ERF_DataStruct.H:128