65 c_dyn, test_filter, mij, lij, num, den)
66 real(kind=
rp),
intent(in) :: t
67 integer,
intent(in) :: tstep
68 type(
coef_t),
intent(in) :: coef
69 type(
field_t),
intent(inout) :: nut
70 type(
field_t),
intent(in) :: delta
71 type(
field_t),
intent(inout) :: c_dyn
73 type(
field_t),
intent(inout) :: mij(6), lij(6)
74 type(
field_t),
intent(inout) :: num, den
76 type(
field_t),
pointer :: u, v, w
79 type(
field_t),
pointer :: s11, s22, s33, s12, s13, s23, s_abs
80 real(kind=
rp) :: alpha
81 integer :: temp_indices(7)
84 if (tstep .eq. 1)
then
103 call strain_rate(s11%x, s22%x, s33%x, s12%x, s13%x, s23%x, u, v, w, coef)
105 call coef%gs_h%op(s11%x, s11%dof%size(),
gs_op_add)
106 call coef%gs_h%op(s22%x, s11%dof%size(),
gs_op_add)
107 call coef%gs_h%op(s33%x, s11%dof%size(),
gs_op_add)
108 call coef%gs_h%op(s12%x, s11%dof%size(),
gs_op_add)
109 call coef%gs_h%op(s13%x, s11%dof%size(),
gs_op_add)
110 call coef%gs_h%op(s23%x, s11%dof%size(),
gs_op_add)
112 do concurrent(i = 1:s11%dof%size())
113 s11%x(i,1,1,1) = s11%x(i,1,1,1) * coef%mult(i,1,1,1)
114 s22%x(i,1,1,1) = s22%x(i,1,1,1) * coef%mult(i,1,1,1)
115 s33%x(i,1,1,1) = s33%x(i,1,1,1) * coef%mult(i,1,1,1)
116 s12%x(i,1,1,1) = s12%x(i,1,1,1) * coef%mult(i,1,1,1)
117 s13%x(i,1,1,1) = s13%x(i,1,1,1) * coef%mult(i,1,1,1)
118 s23%x(i,1,1,1) = s23%x(i,1,1,1) * coef%mult(i,1,1,1)
121 do concurrent(i = 1:u%dof%size())
122 s_abs%x(i,1,1,1) = sqrt(2.0_rp * (s11%x(i,1,1,1)*s11%x(i,1,1,1) + &
123 s22%x(i,1,1,1)*s22%x(i,1,1,1) + &
124 s33%x(i,1,1,1)*s33%x(i,1,1,1)) + &
125 4.0_rp * (s12%x(i,1,1,1)*s12%x(i,1,1,1) + &
126 s13%x(i,1,1,1)*s13%x(i,1,1,1) + &
127 s23%x(i,1,1,1)*s23%x(i,1,1,1)))
132 s_abs, test_filter, delta, u%dof%size())
135 do concurrent(i =1:u%dof%size())
136 if (den%x(i,1,1,1) .gt. 0.0_rp)
then
137 c_dyn%x(i,1,1,1) = 0.5_rp * (num%x(i,1,1,1)/den%x(i,1,1,1))
139 c_dyn%x(i,1,1,1) = 0.0_rp
141 c_dyn%x(i,1,1,1) =
max(c_dyn%x(i,1,1,1),0.0_rp)
142 nut%x(i,1,1,1) = c_dyn%x(i,1,1,1) * delta%x(i,1,1,1)**2 * s_abs%x(i,1,1,1)
146 call col2(nut%x, coef%mult, nut%dof%size())
160 type(
field_t),
intent(inout) :: lij(6)
161 type(
field_t),
pointer,
intent(in) :: u, v, w
163 integer,
intent(in) :: n
166 integer :: temp_indices(3)
167 type(
field_t),
pointer :: fu, fv, fw
173 call test_filter%apply(fu, u)
174 call test_filter%apply(fv, v)
175 call test_filter%apply(fw, w)
178 do concurrent(i = 1:n)
179 lij(1)%x(i,1,1,1) = fu%x(i,1,1,1) * fu%x(i,1,1,1)
180 lij(2)%x(i,1,1,1) = fv%x(i,1,1,1) * fv%x(i,1,1,1)
181 lij(3)%x(i,1,1,1) = fw%x(i,1,1,1) * fw%x(i,1,1,1)
182 lij(4)%x(i,1,1,1) = fu%x(i,1,1,1) * fv%x(i,1,1,1)
183 lij(5)%x(i,1,1,1) = fu%x(i,1,1,1) * fw%x(i,1,1,1)
184 lij(6)%x(i,1,1,1) = fv%x(i,1,1,1) * fw%x(i,1,1,1)
190 call col3(fu%x, u%x, u%x, n)
191 call test_filter%apply(fv, fu)
192 call sub2(lij(1)%x, fv%x, n)
194 call col3(fu%x, v%x, v%x, n)
195 call test_filter%apply(fv, fu)
196 call sub2(lij(2)%x, fv%x, n)
198 call col3(fu%x, w%x, w%x, n)
199 call test_filter%apply(fv, fu)
200 call sub2(lij(3)%x, fv%x, n)
202 call col3(fu%x, u%x, v%x, n)
203 call test_filter%apply(fv, fu)
204 call sub2(lij(4)%x, fv%x, n)
206 call col3(fu%x, u%x, w%x, n)
207 call test_filter%apply(fv, fu)
208 call sub2(lij(5)%x, fv%x, n)
210 call col3(fu%x, v%x, w%x, n)
211 call test_filter%apply(fv, fu)
212 call sub2(lij(6)%x, fv%x, n)
225 s_abs, test_filter, delta, n)
226 type(
field_t),
intent(inout) :: mij(6)
227 type(
field_t),
intent(inout) :: s11, s22, s33, s12, s13, s23, s_abs
229 type(
field_t),
intent(in) :: delta
230 integer,
intent(in) :: n
232 integer :: temp_indices(7)
233 type(
field_t),
pointer :: fs11, fs22, fs33, fs12, fs13, fs23, fs_abs
234 real(kind=
rp) :: delta_ratio2
236 real(kind=
rp) :: delta2
238 delta_ratio2 = ((test_filter%nx-1.0_rp)/(test_filter%nt-1.0_rp))**2
250 call test_filter%apply(fs_abs, s_abs)
252 call test_filter%apply(fs11, s11)
253 call col3(mij(1)%x, fs_abs%x, fs11%x, n)
254 call cmult(mij(1)%x, delta_ratio2, n)
256 call test_filter%apply(fs22, s22)
257 call col3(mij(2)%x, fs_abs%x, fs22%x, n)
258 call cmult(mij(2)%x, delta_ratio2, n)
260 call test_filter%apply(fs33, s33)
261 call col3(mij(3)%x, fs_abs%x, fs33%x, n)
262 call cmult(mij(3)%x, delta_ratio2, n)
264 call test_filter%apply(fs12, s12)
265 call col3(mij(4)%x, fs_abs%x, fs12%x, n)
266 call cmult(mij(4)%x, delta_ratio2, n)
268 call test_filter%apply(fs13, s13)
269 call col3(mij(5)%x, fs_abs%x, fs13%x, n)
270 call cmult(mij(5)%x, delta_ratio2, n)
272 call test_filter%apply(fs23, s23)
273 call col3(mij(6)%x, fs_abs%x, fs23%x, n)
274 call cmult(mij(6)%x, delta_ratio2, n)
280 call col3(fs11%x, s_abs%x, s11%x, n)
281 call test_filter%apply(fs22, fs11)
282 call sub2(mij(1)%x, fs22%x, n)
284 call col3(fs11%x, s_abs%x, s22%x, n)
285 call test_filter%apply(fs22, fs11)
286 call sub2(mij(2)%x, fs22%x, n)
288 call col3(fs11%x, s_abs%x, s33%x, n)
289 call test_filter%apply(fs22, fs11)
290 call sub2(mij(3)%x, fs22%x, n)
292 call col3(fs11%x, s_abs%x, s12%x, n)
293 call test_filter%apply(fs22, fs11)
294 call sub2(mij(4)%x, fs22%x, n)
296 call col3(fs11%x, s_abs%x, s13%x, n)
297 call test_filter%apply(fs22, fs11)
298 call sub2(mij(5)%x, fs22%x, n)
300 call col3(fs11%x, s_abs%x, s23%x, n)
301 call test_filter%apply(fs22, fs11)
302 call sub2(mij(6)%x, fs22%x, n)
305 do concurrent(i = 1:n)
306 delta2 = delta%x(i,1,1,1)**2
307 mij(1)%x(i,1,1,1) = mij(1)%x(i,1,1,1) * delta2
308 mij(2)%x(i,1,1,1) = mij(2)%x(i,1,1,1) * delta2
309 mij(3)%x(i,1,1,1) = mij(3)%x(i,1,1,1) * delta2
310 mij(4)%x(i,1,1,1) = mij(4)%x(i,1,1,1) * delta2
311 mij(5)%x(i,1,1,1) = mij(5)%x(i,1,1,1) * delta2
312 mij(6)%x(i,1,1,1) = mij(6)%x(i,1,1,1) * delta2
324 type(
field_t),
intent(inout) :: num, den
325 type(
field_t),
intent(in) :: lij(6), mij(6)
326 real(kind=
rp),
intent(in) :: alpha
327 integer,
intent(in) :: n
329 real(kind=
rp),
dimension(n) :: num_curr, den_curr
332 do concurrent(i = 1:n)
333 num_curr(i) = mij(1)%x(i,1,1,1)*lij(1)%x(i,1,1,1) + &
334 mij(2)%x(i,1,1,1)*lij(2)%x(i,1,1,1) + &
335 mij(3)%x(i,1,1,1)*lij(3)%x(i,1,1,1) + &
336 2.0_rp*(mij(4)%x(i,1,1,1)*lij(4)%x(i,1,1,1) + &
337 mij(5)%x(i,1,1,1)*lij(5)%x(i,1,1,1) + &
338 mij(6)%x(i,1,1,1)*lij(6)%x(i,1,1,1))
339 den_curr(i) = mij(1)%x(i,1,1,1)*mij(1)%x(i,1,1,1) + &
340 mij(2)%x(i,1,1,1)*mij(2)%x(i,1,1,1) + &
341 mij(3)%x(i,1,1,1)*mij(3)%x(i,1,1,1) + &
342 2.0_rp*(mij(4)%x(i,1,1,1)*mij(4)%x(i,1,1,1) + &
343 mij(5)%x(i,1,1,1)*mij(5)%x(i,1,1,1) + &
344 mij(6)%x(i,1,1,1)*mij(6)%x(i,1,1,1))
348 do concurrent(i = 1:n)
349 num%x(i,1,1,1) = alpha * num%x(i,1,1,1) + (1.0_rp - alpha) * num_curr(i)
350 den%x(i,1,1,1) = alpha * den%x(i,1,1,1) + (1.0_rp - alpha) * den_curr(i)