MFC
Exascale flow solver
Loading...
Searching...
No Matches
m_bubbles.fpp.f90
Go to the documentation of this file.
1# 1 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
2!>
3!! @file
4!! @brief Contains module m_bubbles
5
6# 1 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 1
7# 1 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 1
8# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
9# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
10# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
11# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
12# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
13# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
14
15# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
16# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
17# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
18
19# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
20
21# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
22
23# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
24
25# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
26
27# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
28
29# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
30
31# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
32
33# 145 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
34! New line at end of file is required for FYPP
35# 2 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
36# 1 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 1
37# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
38# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
39# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
40# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
41# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
42# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
43
44# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
45# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
46# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
47
48# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
49
50# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
51
52# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
53
54# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
55
56# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
57
58# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
59
60# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
61
62# 145 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
63! New line at end of file is required for FYPP
64# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
65
66# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
67# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
68# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
69# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
70# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
71
72# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
73
74# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
75
76# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
77
78# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
79
80# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
81
82# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
83
84# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
85
86# 76 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
87
88# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
89
90# 86 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
91
92# 91 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
93
94# 96 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
95
96# 101 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
97
98# 106 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
99
100# 111 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
101
102# 116 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
103
104# 121 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
105
106# 151 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
107
108# 192 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
109
110# 206 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
111
112# 231 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
113
114# 242 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
115
116# 244 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
117# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
118
119# 284 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
120
121# 294 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
122
123# 304 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
124
125# 313 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
126
127# 330 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
128
129# 340 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
130
131# 347 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
132
133# 353 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
134
135# 359 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
136
137# 365 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
138
139# 371 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
140
141# 377 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
142! New line at end of file is required for FYPP
143# 3 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
144# 1 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 1
145# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
146# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
147# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
148# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
149# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
150# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
151
152# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
153# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
154# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
155
156# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
157
158# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
159
160# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
161
162# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
163
164# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
165
166# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
167
168# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
169
170# 145 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
171! New line at end of file is required for FYPP
172# 2 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 2
173
174# 7 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
175
176# 17 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
177
178# 22 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
179
180# 27 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
181
182# 32 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
183
184# 37 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
185
186# 42 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
187
188# 47 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
189
190# 52 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
191
192# 57 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
193
194# 62 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
195
196# 73 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
197
198# 78 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
199
200# 83 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
201
202# 88 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
203
204# 103 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
205
206# 131 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
207
208# 160 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
209
210# 175 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
211
212# 193 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
213
214# 215 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
215
216# 244 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
217
218# 259 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
219
220# 269 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
221
222# 278 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
223
224# 294 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
225
226# 304 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
227
228# 311 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
229! New line at end of file is required for FYPP
230# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
231
232! GPU parallel region (scalar reductions, maxval/minval)
233# 23 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
234
235! GPU parallel loop over threads (most common GPU macro)
236# 43 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
237
238! Required closing for GPU_PARALLEL_LOOP
239# 55 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
240
241! Mark routine for device compilation
242# 112 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
243
244! Declare device-resident data
245# 130 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
246
247! Inner loop within a GPU parallel region
248# 145 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
249
250! Scoped GPU data region
251# 164 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
252
253! Host code with device pointers (for MPI with GPU buffers)
254# 193 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
255
256! Allocate device memory (unscoped)
257# 207 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
258
259! Free device memory
260# 219 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
261
262! Atomic operation on device
263# 231 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
264
265! End atomic capture block
266# 242 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
267
268! Copy data between host and device
269# 254 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
270
271! Synchronization barrier
272# 266 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
273
274! Import GPU library module (openacc or omp_lib)
275# 275 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
276
277! Emit code only for AMD compiler
278# 282 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
279
280! Emit code for non-Cray compilers
281# 289 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
282
283! Emit code only for Cray compiler
284# 296 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
285
286! Emit code for non-NVIDIA compilers
287# 303 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
288
289# 305 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
290# 306 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
291! New line at end of file is required for FYPP
292# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
293
294# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
295
296! Caution: This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI
297! rank. That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0. For an
298! example see misc/nvidia_uvm/bind.sh.
299# 57 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
300
301! Allocate and create GPU device memory
302# 77 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
303
304! Free GPU device memory and deallocate
305# 85 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
306
307! Cray-specific GPU pointer setup for vector fields
308# 109 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
309
310! Cray-specific GPU pointer setup for scalar fields
311# 125 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
312
313! Cray-specific GPU pointer setup for acoustic source spatials
314# 150 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
315
316# 156 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
317
318# 163 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
319! New line at end of file is required for FYPP
320# 6 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp" 2
321
322!> @brief Bubble-dynamics procedures for ensemble- and volume-averaged model
323module m_bubbles
324
325 use m_derived_types
326 use m_global_parameters
327 use m_mpi_proxy
328 use m_variables_conversion
329 use m_helper_basic
330 use m_bubbles_el_kernels
331 use m_constants, only: bubble_model_gilmore, bubble_model_keller_miksis, bubble_model_rayleigh_plesset
332
333 implicit none
334
335contains
336
337 !> Compute the bubble radial acceleration based on the selected bubble model
338 function f_rddot(fRho, fP, fR, fV, fR0, fpb, fpbdot, alf, fntait, fBtait, f_bub_adv_src, f_divu, fCson)
339
340
341# 25 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
342#if MFC_OpenACC
343# 25 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
344!$acc routine seq
345# 25 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
346#elif MFC_OpenMP
347# 25 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
348
349# 25 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
350
351# 25 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
352!$omp declare target device_type(any)
353# 25 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
354#endif
355 real(wp), intent(in) :: frho, fp, fr, fv, fr0, fpb, fpbdot, alf
356 real(wp), intent(in) :: fntait, fbtait, f_bub_adv_src, f_divu
357 real(wp), intent(in) :: fcson
358 real(wp) :: fcpbw, fcpinf, fcpinf_dot, fh, fhdot, c_gas, c_liquid
359 real(wp) :: f_rddot
360
361 if (bubble_model == bubble_model_gilmore) then
362 ! Gilmore bubbles
363 fcpinf = fp - eu
364 fcpbw = f_cpbw(fr0, fr, fv, fpb)
365 fh = f_h(fcpbw, fcpinf, fntait, fbtait)
366 c_gas = f_cgas(fcpinf, fntait, fbtait, fh)
367 fcpinf_dot = f_cpinfdot(frho, fp, alf, fntait, fbtait, f_bub_adv_src, f_divu)
368 fhdot = f_hdot(fcpbw, fcpinf, fcpinf_dot, fntait, fbtait, fr, fv, fr0, fpbdot)
369 f_rddot = f_rddot_g(fcpbw, fr, fv, fh, fhdot, c_gas, fntait, fbtait)
370 else if (bubble_model == bubble_model_keller_miksis) then
371 ! Keller-Miksis bubbles
372 fcpinf = fp
373 fcpbw = f_cpbw_km(fr0, fr, fv, fpb)
374 if (bubbles_euler) then
375 c_liquid = sqrt(fntait*(fp + fbtait)/(frho*(1._wp - alf)))
376 else
377 c_liquid = fcson
378 end if
379 f_rddot = f_rddot_km(fpbdot, fcpinf, fcpbw, frho, fr, fv, fr0, c_liquid)
380 else if (bubble_model == bubble_model_rayleigh_plesset) then
381 ! Rayleigh-Plesset bubbles
382 fcpbw = f_cpbw_km(fr0, fr, fv, fpb)
383 f_rddot = f_rddot_rp(fp, frho, fr, fv, fcpbw)
384 else
385 ! Default: No bubble dynamics
386 f_rddot = 0._wp
387 end if
388
389 end function f_rddot
390
391 !> Bubble wall pressure: stiffened gas with Laplace pressure and viscous stress
392 function f_cpbw(fR0, fR, fV, fpb)
393
394
395# 65 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
396#if MFC_OpenACC
397# 65 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
398!$acc routine seq
399# 65 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
400#elif MFC_OpenMP
401# 65 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
402
403# 65 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
404
405# 65 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
406!$omp declare target device_type(any)
407# 65 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
408#endif
409 real(wp), intent(in) :: fr0, fr, fv, fpb
410 real(wp) :: f_cpbw
411
412 if (polytropic) then
413 f_cpbw = (ca + 2._wp/web/fr0)*((fr0/fr)**(3._wp*gam)) - ca - 4._wp*re_inv*fv/fr - 2._wp/(fr*web)
414 else
415 f_cpbw = fpb - 1._wp - 4._wp*re_inv*fv/fr - 2._wp/(fr*web)
416 end if
417
418 end function f_cpbw
419
420 !> Compute the bubble enthalpy
421 function f_h(fCpbw, fCpinf, fntait, fBtait)
422
423
424# 80 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
425#if MFC_OpenACC
426# 80 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
427!$acc routine seq
428# 80 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
429#elif MFC_OpenMP
430# 80 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
431
432# 80 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
433
434# 80 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
435!$omp declare target device_type(any)
436# 80 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
437#endif
438 real(wp), intent(in) :: fcpbw, fcpinf, fntait, fbtait
439 real(wp) :: tmp1, tmp2, tmp3
440 real(wp) :: f_h
441
442 tmp1 = (fntait - 1._wp)/fntait
443 tmp2 = (fcpbw/(1._wp + fbtait) + 1._wp)**tmp1
444 tmp3 = (fcpinf/(1._wp + fbtait) + 1._wp)**tmp1
445
446 f_h = (tmp2 - tmp3)*fntait*(1._wp + fbtait)/(fntait - 1._wp)
447
448 end function f_h
449
450 !> Compute the sound speed for the bubble
451 function f_cgas(fCpinf, fntait, fBtait, fH)
452
453
454# 96 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
455#if MFC_OpenACC
456# 96 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
457!$acc routine seq
458# 96 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
459#elif MFC_OpenMP
460# 96 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
461
462# 96 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
463
464# 96 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
465!$omp declare target device_type(any)
466# 96 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
467#endif
468 real(wp), intent(in) :: fcpinf, fntait, fbtait, fh
469 real(wp) :: tmp
470 real(wp) :: f_cgas
471
472 ! get sound speed for Gilmore equations "C" -> c_gas
473 tmp = (fcpinf/(1._wp + fbtait) + 1._wp)**((fntait - 1._wp)/fntait)
474 tmp = fntait*(1._wp + fbtait)*tmp
475
476 f_cgas = sqrt(tmp + (fntait - 1._wp)*fh)
477
478 end function f_cgas
479
480 !> Compute the time derivative of the driving pressure
481 function f_cpinfdot(fRho, fP, falf, fntait, fBtait, advsrc, divu)
482
483
484# 112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
485#if MFC_OpenACC
486# 112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
487!$acc routine seq
488# 112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
489#elif MFC_OpenMP
490# 112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
491
492# 112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
493
494# 112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
495!$omp declare target device_type(any)
496# 112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
497#endif
498 real(wp), intent(in) :: frho, fp, falf, fntait, fbtait, advsrc, divu
499 real(wp) :: c2_liquid
500 real(wp) :: f_cpinfdot
501
502 ! get sound speed squared for liquid (only needed for pbdot) c_l^2 = gam (p+B) / (rho*(1-alf))
503 if (mpp_lim) then
504 c2_liquid = fntait*(fp + fbtait)/frho
505 else
506 c2_liquid = fntait*(fp + fbtait)/(frho*(1._wp - falf))
507 end if
508
509 ! \dot{Cp_inf} = rho sound^2 (alf_src - divu)
510 f_cpinfdot = frho*c2_liquid*(advsrc - divu)
511
512 end function f_cpinfdot
513
514 !> Enthalpy derivative for Gilmore bubble model, Gilmore (1952)
515 function f_hdot(fCpbw, fCpinf, fCpinf_dot, fntait, fBtait, fR, fV, fR0, fpbdot)
516
517
518# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
519#if MFC_OpenACC
520# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
521!$acc routine seq
522# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
523#elif MFC_OpenMP
524# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
525
526# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
527
528# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
529!$omp declare target device_type(any)
530# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
531#endif
532 real(wp), intent(in) :: fcpbw, fcpinf, fcpinf_dot, fntait, fbtait
533 real(wp), intent(in) :: fr, fv, fr0, fpbdot
534 real(wp) :: tmp1, tmp2
535 real(wp) :: f_hdot
536
537 if (polytropic) then
538 tmp1 = (fr0/fr)**(3._wp*gam)
539 tmp1 = -3._wp*gam*(ca + 2._wp/web/fr0)*tmp1*fv/fr
540 else
541 tmp1 = fpbdot
542 end if
543 tmp2 = (2._wp/web + 4._wp*re_inv*fv)*fv/(fr**2._wp)
544
545 f_hdot = (fcpbw/(1._wp + fbtait) + 1._wp)**(-1._wp/fntait)*(tmp1 + tmp2) - (fcpinf/(1._wp + fbtait) + 1._wp) &
546 & **(-1._wp/fntait)*fcpinf_dot
547
548 end function f_hdot
549
550 !> Rayleigh-Plesset bubble radial acceleration
551 function f_rddot_rp(fCp, fRho, fR, fV, fCpbw)
552
553
554# 154 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
555#if MFC_OpenACC
556# 154 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
557!$acc routine seq
558# 154 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
559#elif MFC_OpenMP
560# 154 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
561
562# 154 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
563
564# 154 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
565!$omp declare target device_type(any)
566# 154 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
567#endif
568 real(wp), intent(in) :: fcp, frho, fr, fv, fcpbw
569 real(wp) :: f_rddot_rp
570
571 f_rddot_rp = (-1.5_wp*(fv**2._wp) + (fcpbw - fcp)/frho)/fr
572
573 end function f_rddot_rp
574
575 !> Compute the Gilmore bubble radial acceleration
576 function f_rddot_g(fCpbw, fR, fV, fH, fHdot, fcgas, fntait, fBtait)
577
578
579# 165 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
580#if MFC_OpenACC
581# 165 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
582!$acc routine seq
583# 165 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
584#elif MFC_OpenMP
585# 165 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
586
587# 165 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
588
589# 165 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
590!$omp declare target device_type(any)
591# 165 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
592#endif
593 real(wp), intent(in) :: fcpbw, fr, fv, fh, fhdot
594 real(wp), intent(in) :: fcgas, fntait, fbtait
595 real(wp) :: tmp1, tmp2, tmp3
596 real(wp) :: f_rddot_g
597
598 tmp1 = fv/fcgas
599 tmp2 = 1._wp + 4._wp*re_inv/fcgas/fr*(fcpbw/(1._wp + fbtait) + 1._wp)**(-1._wp/fntait)
600 tmp3 = 1.5_wp*fv**2._wp*(tmp1/3._wp - 1._wp) + fh*(1._wp + tmp1) + fr*fhdot*(1._wp - tmp1)/fcgas
601
602 f_rddot_g = tmp3/(fr*(1._wp - tmp1)*tmp2)
603
604 end function f_rddot_g
605
606 !> Keller-Miksis bubble wall pressure
607 function f_cpbw_km(fR0, fR, fV, fpb)
608
609
610# 182 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
611#if MFC_OpenACC
612# 182 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
613!$acc routine seq
614# 182 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
615#elif MFC_OpenMP
616# 182 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
617
618# 182 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
619
620# 182 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
621!$omp declare target device_type(any)
622# 182 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
623#endif
624 real(wp), intent(in) :: fr0, fr, fv, fpb
625 real(wp) :: f_cpbw_km
626
627 if (polytropic) then
628 f_cpbw_km = ca*((fr0/fr)**(3._wp*gam)) - ca + eu
629 if (.not. f_is_default(web)) f_cpbw_km = f_cpbw_km + (2._wp/(web*fr0))*((fr0/fr)**(3._wp*gam))
630 else
631 f_cpbw_km = fpb
632 end if
633
634 if (.not. f_is_default(web)) f_cpbw_km = f_cpbw_km - 2._wp/(fr*web)
635 if (.not. f_is_default(re_inv)) f_cpbw_km = f_cpbw_km - 4._wp*re_inv*fv/fr
636
637 end function f_cpbw_km
638
639 !> Keller-Miksis bubble radial acceleration
640 function f_rddot_km(fpbdot, fCp, fCpbw, fRho, fR, fV, fR0, fC)
641
642
643# 201 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
644#if MFC_OpenACC
645# 201 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
646!$acc routine seq
647# 201 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
648#elif MFC_OpenMP
649# 201 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
650
651# 201 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
652
653# 201 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
654!$omp declare target device_type(any)
655# 201 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
656#endif
657 real(wp), intent(in) :: fpbdot, fcp, fcpbw
658 real(wp), intent(in) :: frho, fr, fv, fr0, fc
659 real(wp) :: tmp1, tmp2, cdot_star
660 real(wp) :: f_rddot_km
661 if (polytropic) then
662 cdot_star = -3._wp*gam*ca*((fr0/fr)**(3._wp*gam))*fv/fr
663 if (.not. f_is_default(web)) cdot_star = cdot_star - 3._wp*gam*(2._wp/(web*fr0))*((fr0/fr)**(3._wp*gam))*fv/fr
664 else
665 cdot_star = fpbdot
666 end if
667
668 if (.not. f_is_default(web)) cdot_star = cdot_star + (2._wp/web)*fv/(fr**2._wp)
669 if (.not. f_is_default(re_inv)) cdot_star = cdot_star + 4._wp*re_inv*((fv/fr)**2._wp)
670
671 tmp1 = fv/fc
672 tmp2 = 1.5_wp*(fv**2._wp)*(tmp1/3._wp - 1._wp) + (1._wp + tmp1)*(fcpbw - fcp)/frho + cdot_star*fr/(frho*fc)
673
674 if (f_is_default(re_inv)) then
675 f_rddot_km = tmp2/(fr*(1._wp - tmp1))
676 else
677 f_rddot_km = tmp2/(fr*(1._wp - tmp1) + 4._wp*re_inv/(frho*fc))
678 end if
679
680 end function f_rddot_km
681
682 !> Compute bubble wall properties for vapor bubbles
683 subroutine s_bwproperty(pb_in, iR0, chi_vw_out, k_mw_out, rho_mw_out)
684
685
686# 230 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
687#if MFC_OpenACC
688# 230 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
689!$acc routine seq
690# 230 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
691#elif MFC_OpenMP
692# 230 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
693
694# 230 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
695
696# 230 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
697!$omp declare target device_type(any)
698# 230 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
699#endif
700 real(wp), intent(in) :: pb_in
701 integer, intent(in) :: iR0
702 real(wp), intent(out) :: chi_vw_out !< Bubble wall properties (Ando 2010)
703 real(wp), intent(out) :: k_mw_out !< Bubble wall properties (Ando 2010)
704 real(wp), intent(out) :: rho_mw_out !< Bubble wall properties (Ando 2010)
705 real(wp) :: x_vw
706
707 ! mass fraction of vapor
708 chi_vw_out = 1._wp/(1._wp + r_v/r_g*(pb_in/pv - 1._wp))
709 ! mole fraction of vapor & thermal conductivity of gas mixture
710 x_vw = m_g*chi_vw_out/(m_v + (m_g - m_v)*chi_vw_out)
711 k_mw_out = x_vw*k_v(ir0)/(x_vw + (1._wp - x_vw)*phi_vg) + (1._wp - x_vw)*k_g(ir0)/(x_vw*phi_gv + 1._wp - x_vw)
712 ! gas mixture density
713 rho_mw_out = pv/(chi_vw_out*r_v*tw)
714
715 end subroutine s_bwproperty
716
717 !> Compute the vapour flux
718 subroutine s_vflux(fR, fV, fpb, fmass_v, iR0, vflux, fmass_g, fbeta_c, fR_m, fgamma_m, fchi_vw, frho_mw)
719
720
721# 251 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
722#if MFC_OpenACC
723# 251 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
724!$acc routine seq
725# 251 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
726#elif MFC_OpenMP
727# 251 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
728
729# 251 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
730
731# 251 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
732!$omp declare target device_type(any)
733# 251 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
734#endif
735 real(wp), intent(in) :: fR
736 real(wp), intent(in) :: fV
737 real(wp), intent(in) :: fpb
738 real(wp), intent(in) :: fmass_v
739 integer, intent(in) :: iR0
740 real(wp), intent(out) :: vflux
741 real(wp), intent(in), optional :: fmass_g, fbeta_c
742 real(wp), intent(out), optional :: fR_m, fgamma_m
743 real(wp), intent(in), optional :: fchi_vw, frho_mw !< Per-thread bubble-wall scratch (Euler-Euler path)
744 real(wp) :: chi_bar
745 real(wp) :: rho_mw_lag
746 real(wp) :: grad_chi
747 real(wp) :: conc_v
748
749 if (thermal == 3) then ! transfer
750 ! constant transfer model
751 if (bubbles_lagrange) then
752 ! Mixture properties (gas+vapor) in the bubble
753 conc_v = fmass_v/(fmass_v + fmass_g)
754 if (lag_params%massTransfer_model) then
755 conc_v = 1._wp/(1._wp + (r_v/r_g)*(fpb/pv - 1._wp))
756 end if
757 fr_m = (fmass_g*r_g + fmass_v*r_v)
758 fgamma_m = conc_v*gam_v + (1._wp - conc_v)*gam_g
759
760 ! Vapor flux
761 chi_bar = fmass_v/(fmass_v + fmass_g)
762 grad_chi = (chi_bar - conc_v)
763 rho_mw_lag = (fmass_g + fmass_v)/(4._wp/3._wp*pi*fr**3._wp)
764 vflux = 0._wp
765 if (lag_params%massTransfer_model) then
766 vflux = -fbeta_c*rho_mw_lag*grad_chi/(1._wp - conc_v)/fr
767 end if
768 else
769 chi_bar = fmass_v/(fmass_v + mass_g0(ir0))
770 grad_chi = -re_trans_c(ir0)*(chi_bar - fchi_vw)
771 vflux = frho_mw*grad_chi/pe_c/(1._wp - fchi_vw)/fr
772 end if
773 else
774 ! polytropic
775 vflux = pv*fv/(r_v*tw)
776 end if
777
778 end subroutine s_vflux
779
780 !> Compute the time derivative of the internal bubble pressure
781 function f_bpres_dot(fvflux, fR, fV, fpb, fmass_v, iR0, fbeta_t, fR_m, fgamma_m, fk_mw)
782
783
784# 300 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
785#if MFC_OpenACC
786# 300 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
787!$acc routine seq
788# 300 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
789#elif MFC_OpenMP
790# 300 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
791
792# 300 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
793
794# 300 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
795!$omp declare target device_type(any)
796# 300 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
797#endif
798 real(wp), intent(in) :: fvflux
799 real(wp), intent(in) :: fr
800 real(wp), intent(in) :: fv
801 real(wp), intent(in) :: fpb
802 real(wp), intent(in) :: fmass_v
803 integer, intent(in) :: ir0
804 real(wp), intent(in), optional :: fbeta_t, fr_m, fgamma_m
805 real(wp), intent(in), optional :: fk_mw !< Per-thread bubble-wall conductivity scratch (Euler-Euler path)
806 real(wp) :: t_bar
807 real(wp) :: grad_t
808 real(wp) :: f_bpres_dot
809 real(wp) :: heatflux
810
811 if (thermal == 3) then
812 if (bubbles_lagrange) then
813 t_bar = fpb*(4._wp/3._wp*pi*fr**3._wp)/fr_m
814 grad_t = -fbeta_t*(t_bar - tw)
815 heatflux = (fgamma_m - 1._wp)/fgamma_m*grad_t/fr
816 f_bpres_dot = 3._wp*fgamma_m*(-fv*fpb + fvflux*r_v*tw + heatflux)/fr
817 return
818 end if
819 grad_t = -re_trans_t(ir0)*((fpb/pb0(ir0))*(fr/r0(ir0))**3*(mass_g0(ir0) + mass_v0(ir0))/(mass_g0(ir0) + fmass_v) &
820 & - 1._wp)
821 f_bpres_dot = 3._wp*gam_m*(-fv*fpb + fvflux*r_v*tw + pb0(ir0)*fk_mw*grad_t/pe_t(ir0)/fr)/fr
822 else
823 f_bpres_dot = -3._wp*gam_m*fv/fr*(fpb - pv)
824 end if
825
826 end function f_bpres_dot
827
828 !> Adaptive time stepping routine for subgrid bubbles (See Heirer, E. Hairer S.P.Norsett G. Wanner, Solving Ordinary
829 !! Differential Equations I, Chapter II.4)
830 function f_advance_step(fRho, fP, fR, fV, fR0, fpb, fpbdot, alf, fntait, fBtait, f_bub_adv_src, f_divu, bub_id, fmass_v, &
831 & fmass_g, fbeta_c, fbeta_t, fCson, fRe, fPos, fVel, cell, q_prim_vf) result(adap_dt_stop)
832
833# 335 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
834#if MFC_OpenACC
835# 335 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
836!$acc routine seq
837# 335 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
838#elif MFC_OpenMP
839# 335 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
840
841# 335 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
842
843# 335 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
844!$omp declare target device_type(any)
845# 335 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
846#endif
847
848 real(wp), intent(inout) :: fr, fv, fpb, fmass_v
849 real(wp), intent(in) :: frho, fp, fr0, fpbdot, alf
850 real(wp), intent(in) :: fntait, fbtait, f_bub_adv_src, f_divu
851 integer, intent(in) :: bub_id
852 real(wp), intent(in) :: fmass_g, fbeta_c, fbeta_t, fcson
853 real(wp), intent(inout), dimension(3), optional :: fpos, fvel
854 real(wp), intent(in), optional :: fre
855 integer, intent(in), dimension(3), optional :: cell
856 type(scalar_field), intent(in), dimension(sys_size), optional :: q_prim_vf
857 real(wp), dimension(5) :: err !< Error estimates for adaptive time stepping
858 real(wp) :: t_new !< Updated time step size
859 real(wp) :: h0, h !< Time step size
860 !> Bubble radius, radial velocity, and radial acceleration for the inner loop
861 real(wp), dimension(4) :: myr_tmp1, myv_tmp1, myr_tmp2, myv_tmp2
862 real(wp), dimension(4) :: mypb_tmp1, mymv_tmp1, mypb_tmp2, mymv_tmp2 !< Gas pressure and vapor mass for the inner loop (EL)
863 real(wp) :: fr2, fv2, fpb2, fmass_v2, f_btemp
864 real(wp), dimension(3) :: vtemp, atemp
865 integer :: adap_dt_stop
866 integer :: l, iter_count
867
868 call s_initial_substep_h(frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, fcson, h0)
869 h = h0
870 ! Advancing one step
871 t_new = 0._wp
872 iter_count = 0
873 adap_dt_stop = 0
874
875 do
876 if (t_new + h > 0.5_wp*dt) then
877 h = 0.5_wp*dt - t_new
878 end if
879
880 ! Advancing one sub-step
881 do while (iter_count < adap_dt_max_iters)
882 iter_count = iter_count + 1
883
884 ! Advance one sub-step
885 call s_advance_substep(err(1), frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, &
886 & bub_id, fmass_v, fmass_g, fbeta_c, fbeta_t, fcson, h, myr_tmp1, myv_tmp1, mypb_tmp1, &
887 & mymv_tmp1)
888 if (err(1) > adap_dt_tol) then
889 h = 0.25_wp*h
890 cycle
891 end if
892
893 ! Advance one sub-step by advancing two half steps
894 call s_advance_substep(err(2), frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, &
895 & bub_id, fmass_v, fmass_g, fbeta_c, fbeta_t, fcson, 0.5_wp*h, myr_tmp2, myv_tmp2, &
896 & mypb_tmp2, mymv_tmp2)
897 if (err(2) > adap_dt_tol) then
898 h = 0.25_wp*h
899 cycle
900 end if
901
902 fr2 = myr_tmp2(4); fv2 = myv_tmp2(4)
903 fpb2 = mypb_tmp2(4); fmass_v2 = mymv_tmp2(4)
904
905 call s_advance_substep(err(3), frho, fp, fr2, fv2, fr0, fpb2, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, &
906 & bub_id, fmass_v2, fmass_g, fbeta_c, fbeta_t, fcson, 0.5_wp*h, myr_tmp2, myv_tmp2, &
907 & mypb_tmp2, mymv_tmp2)
908 if (err(3) > adap_dt_tol) then
909 h = 0.5_wp*h
910 cycle
911 end if
912
913 err(4) = abs((myr_tmp1(4) - myr_tmp2(4))/myr_tmp1(4))
914 err(5) = abs((myv_tmp1(4) - myv_tmp2(4))/myv_tmp1(4))
915 if (abs(myv_tmp1(4)) < verysmall) err(5) = 0._wp
916
917 ! Determine acceptance/rejection and update step size Rule 1: err1, err2, err3 < tol Rule 2: myR_tmp1(4) > 0._wp
918 ! Rule 3: abs((myR_tmp1(4) - myR_tmp2(4))/fR) < tol Rule 4: abs((myV_tmp1(4) - myV_tmp2(4))/fV) < tol
919 if ((err(1) <= adap_dt_tol) .and. (err(2) <= adap_dt_tol) .and. (err(3) <= adap_dt_tol) .and. (err(4) &
920 & <= adap_dt_tol) .and. (err(5) <= adap_dt_tol) .and. myr_tmp1(4) > 0._wp) then
921 ! Accepted. Finalize the sub-step
922 t_new = t_new + h
923
924 ! Update R and V
925 fr = myr_tmp1(4)
926 fv = myv_tmp1(4)
927
928 if (bubbles_lagrange) then
929 ! Update pb and mass_v
930 fpb = mypb_tmp1(4)
931 fmass_v = mymv_tmp1(4)
932
933 select case (lag_vel_model)
934 case (1)
935 do l = 1, num_dims
936 if (fd_order > 1) then
937 vtemp(l) = f_interpolate_velocity(fpos(l), cell, l, q_prim_vf)
938 else
939 vtemp(l) = q_prim_vf(eqn_idx%mom%beg + l - 1)%sf(cell(1), cell(2), cell(3))
940 end if
941 end do
942 do l = 1, num_dims
943 fvel(l) = vtemp(l)
944 fpos(l) = fpos(l) + h*vtemp(l)
945 end do
946 case (2)
947 do l = 1, num_dims
948 f_btemp = f_get_bubble_force(fpos(l), fr, fv, fvel(l), fmass_g, fmass_v, fre, frho, cell, l, &
949 & q_prim_vf)
950 atemp(l) = f_btemp/(fmass_g + fmass_v)
951 end do
952 do l = 1, num_dims
953 fvel(l) = fvel(l) + h*atemp(l)
954 fpos(l) = fpos(l) + h*fvel(l)
955 end do
956 end select
957 end if
958
959 ! Update step size for the next sub-step
960 h = h*min(2._wp, max(0.5_wp, (adap_dt_tol/err(1))**(1._wp/3._wp)))
961
962 exit
963 else
964 ! Rejected. Update step size for the next try on sub-step
965 if (err(2) <= adap_dt_tol) then
966 h = 0.5_wp*h
967 else
968 h = 0.25_wp*h
969 end if
970 end if
971 end do
972
973 ! Exit the loop if the final time reached dt
974 if (f_approx_equal(t_new, 0.5_wp*dt) .or. iter_count >= adap_dt_max_iters) exit
975 end do
976
977 if (iter_count >= adap_dt_max_iters) adap_dt_stop = 1
978
979 end function f_advance_step
980
981 !> Choose the initial time step size for the adaptive time stepping routine (See Heirer, E. Hairer S.P.Norsett G. Wanner,
982 !! Solving Ordinary Differential Equations I, Chapter II.4)
983 subroutine s_initial_substep_h(fRho, fP, fR, fV, fR0, fpb, fpbdot, alf, fntait, fBtait, f_bub_adv_src, f_divu, fCson, h)
984
985
986# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
987#ifdef _CRAYFTN
988# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
989#if MFC_OpenACC
990# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
991!$acc routine seq
992# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
993#elif MFC_OpenMP
994# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
995
996# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
997
998# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
999!$omp declare target device_type(any)
1000# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1001#else
1002# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1003!DIR$ INLINEALWAYS s_initial_substep_h
1004# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1005#endif
1006# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1007#elif MFC_OpenACC
1008# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1009!$acc routine seq
1010# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1011#elif MFC_OpenMP
1012# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1013
1014# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1015
1016# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1017!$omp declare target device_type(any)
1018# 474 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1019#endif
1020
1021 real(wp), intent(in) :: fRho, fP, fR, fV, fR0, fpb, fpbdot, alf
1022 real(wp), intent(in) :: fntait, fBtait, f_bub_adv_src, f_divu
1023 real(wp), intent(in) :: fCson
1024 real(wp), intent(out) :: h
1025 real(wp), dimension(2) :: h_size !< Time step size (h0, h1)
1026 real(wp), dimension(3) :: d_norms !< norms (d_0, d_1, d_2)
1027 real(wp), dimension(2) :: myR_tmp, myV_tmp, myA_tmp !< Bubble radius, radial velocity, and radial acceleration
1028 ! Determine the starting time step Evaluate f(x0,y0)
1029 myr_tmp(1) = fr
1030 myv_tmp(1) = fv
1031 mya_tmp(1) = f_rddot(frho, fp, myr_tmp(1), myv_tmp(1), fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, fcson)
1032
1033 ! Compute d_0 = ||y0|| and d_1 = ||f(x0,y0)||
1034 d_norms(1) = sqrt((myr_tmp(1)**2._wp + myv_tmp(1)**2._wp)/2._wp)
1035 d_norms(2) = sqrt((myv_tmp(1)**2._wp + mya_tmp(1)**2._wp)/2._wp)
1036 if (d_norms(1) < threshold_first_guess .or. d_norms(2) < threshold_first_guess) then
1037 h_size(1) = small_guess
1038 else
1039 h_size(1) = scale_guess*(d_norms(1)/d_norms(2))
1040 end if
1041
1042 ! Evaluate f(x0+h0,y0+h0*f(x0,y0))
1043 myr_tmp(2) = myr_tmp(1) + h_size(1)*myv_tmp(1)
1044 myv_tmp(2) = myv_tmp(1) + h_size(1)*mya_tmp(1)
1045 mya_tmp(2) = f_rddot(frho, fp, myr_tmp(2), myv_tmp(2), fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, fcson)
1046
1047 ! Compute d_2 = ||f(x0+h0,y0+h0*f(x0,y0))-f(x0,y0)||/h0
1048 d_norms(3) = sqrt(((myv_tmp(2) - myv_tmp(1))**2._wp + (mya_tmp(2) - mya_tmp(1))**2._wp)/2._wp)/h_size(1)
1049
1050 ! Set h1 = (0.01/max(d_1,d_2))^{1/(p+1)} if max(d_1,d_2) < 1.e-15_wp, h_size(2) = max(1.e-6_wp, h0*1.e-3_wp)
1051 if (max(d_norms(2), d_norms(3)) < threshold_second_guess) then
1052 h_size(2) = max(small_guess, h_size(1)*scale_first_guess)
1053 else
1054 h_size(2) = (scale_guess/max(d_norms(2), d_norms(3)))**(1._wp/3._wp)
1055 end if
1056
1057 h = min(h_size(1)/scale_guess, h_size(2))
1058
1059 end subroutine s_initial_substep_h
1060
1061 !> Integrate bubble variables over the given time step size, h, using a third-order accurate embedded Runge-Kutta scheme.
1062 subroutine s_advance_substep(err, fRho, fP, fR, fV, fR0, fpb, fpbdot, alf, fntait, fBtait, f_bub_adv_src, f_divu, bub_id, &
1063 & fmass_v, fmass_g, fbeta_c, fbeta_t, fCson, h, myR_tmp, myV_tmp, myPb_tmp, myMv_tmp)
1064
1065# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1066#ifdef _CRAYFTN
1067# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1068#if MFC_OpenACC
1069# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1070!$acc routine seq
1071# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1072#elif MFC_OpenMP
1073# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1074
1075# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1076
1077# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1078!$omp declare target device_type(any)
1079# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1080#else
1081# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1082!DIR$ INLINEALWAYS s_advance_substep
1083# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1084#endif
1085# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1086#elif MFC_OpenACC
1087# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1088!$acc routine seq
1089# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1090#elif MFC_OpenMP
1091# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1092
1093# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1094
1095# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1096!$omp declare target device_type(any)
1097# 519 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1098#endif
1099
1100 real(wp), intent(out) :: err
1101 real(wp), intent(in) :: fRho, fP, fR, fV, fR0, fpb, fpbdot, alf
1102 real(wp), intent(in) :: fntait, fBtait, f_bub_adv_src, f_divu, h
1103 integer, intent(in) :: bub_id
1104 real(wp), intent(in) :: fmass_v, fmass_g, fbeta_c, fbeta_t, fCson
1105 real(wp), dimension(4), intent(out) :: myR_tmp, myV_tmp, myPb_tmp, myMv_tmp
1106 real(wp), dimension(4) :: myA_tmp, mydPbdt_tmp, mydMvdt_tmp
1107 real(wp) :: err_R, err_V
1108
1109 mypb_tmp(1:4) = fpb
1110 mydpbdt_tmp(1:4) = fpbdot
1111
1112 ! Stage 0
1113 myr_tmp(1) = fr
1114 myv_tmp(1) = fv
1115 if (bubbles_lagrange) then
1116 mypb_tmp(1) = fpb
1117 mymv_tmp(1) = fmass_v
1118 call s_advance_el(myr_tmp(1), myv_tmp(1), mypb_tmp(1), mymv_tmp(1), bub_id, fmass_g, fbeta_c, fbeta_t, &
1119 & mydpbdt_tmp(1), mydmvdt_tmp(1))
1120 end if
1121 mya_tmp(1) = f_rddot(frho, fp, myr_tmp(1), myv_tmp(1), fr0, mypb_tmp(1), mydpbdt_tmp(1), alf, fntait, fbtait, &
1122 & f_bub_adv_src, f_divu, fcson)
1123
1124 ! Stage 1
1125 myr_tmp(2) = myr_tmp(1) + h*myv_tmp(1)
1126 if (myr_tmp(2) < 0._wp) then
1127 err = adap_dt_tol + 1._wp; return
1128 end if
1129 myv_tmp(2) = myv_tmp(1) + h*mya_tmp(1)
1130 if (bubbles_lagrange) then
1131 mypb_tmp(2) = mypb_tmp(1) + h*mydpbdt_tmp(1)
1132 mymv_tmp(2) = mymv_tmp(1) + h*mydmvdt_tmp(1)
1133 call s_advance_el(myr_tmp(2), myv_tmp(2), mypb_tmp(2), mymv_tmp(2), bub_id, fmass_g, fbeta_c, fbeta_t, &
1134 & mydpbdt_tmp(2), mydmvdt_tmp(2))
1135 end if
1136 mya_tmp(2) = f_rddot(frho, fp, myr_tmp(2), myv_tmp(2), fr0, mypb_tmp(2), mydpbdt_tmp(2), alf, fntait, fbtait, &
1137 & f_bub_adv_src, f_divu, fcson)
1138
1139 ! Stage 2
1140 myr_tmp(3) = myr_tmp(1) + (h/4._wp)*(myv_tmp(1) + myv_tmp(2))
1141 if (myr_tmp(3) < 0._wp) then
1142 err = adap_dt_tol + 1._wp; return
1143 end if
1144 myv_tmp(3) = myv_tmp(1) + (h/4._wp)*(mya_tmp(1) + mya_tmp(2))
1145 if (bubbles_lagrange) then
1146 mypb_tmp(3) = mypb_tmp(1) + (h/4._wp)*(mydpbdt_tmp(1) + mydpbdt_tmp(2))
1147 mymv_tmp(3) = mymv_tmp(1) + (h/4._wp)*(mydmvdt_tmp(1) + mydmvdt_tmp(2))
1148 call s_advance_el(myr_tmp(3), myv_tmp(3), mypb_tmp(3), mymv_tmp(3), bub_id, fmass_g, fbeta_c, fbeta_t, &
1149 & mydpbdt_tmp(3), mydmvdt_tmp(3))
1150 end if
1151 mya_tmp(3) = f_rddot(frho, fp, myr_tmp(3), myv_tmp(3), fr0, mypb_tmp(3), mydpbdt_tmp(3), alf, fntait, fbtait, &
1152 & f_bub_adv_src, f_divu, fcson)
1153
1154 ! Stage 3
1155 myr_tmp(4) = myr_tmp(1) + (h/6._wp)*(myv_tmp(1) + myv_tmp(2) + 4._wp*myv_tmp(3))
1156 if (myr_tmp(4) < 0._wp) then
1157 err = adap_dt_tol + 1._wp; return
1158 end if
1159 myv_tmp(4) = myv_tmp(1) + (h/6._wp)*(mya_tmp(1) + mya_tmp(2) + 4._wp*mya_tmp(3))
1160 if (bubbles_lagrange) then
1161 mypb_tmp(4) = mypb_tmp(1) + (h/6._wp)*(mydpbdt_tmp(1) + mydpbdt_tmp(2) + 4._wp*mydpbdt_tmp(3))
1162 mymv_tmp(4) = mymv_tmp(1) + (h/6._wp)*(mydmvdt_tmp(1) + mydmvdt_tmp(2) + 4._wp*mydmvdt_tmp(3))
1163 call s_advance_el(myr_tmp(4), myv_tmp(4), mypb_tmp(4), mymv_tmp(4), bub_id, fmass_g, fbeta_c, fbeta_t, &
1164 & mydpbdt_tmp(4), mydmvdt_tmp(4))
1165 end if
1166 mya_tmp(4) = f_rddot(frho, fp, myr_tmp(4), myv_tmp(4), fr0, mypb_tmp(4), mydpbdt_tmp(4), alf, fntait, fbtait, &
1167 & f_bub_adv_src, f_divu, fcson)
1168
1169 ! Estimate error
1170 err_r = (-5._wp*h/24._wp)*(myv_tmp(2) + myv_tmp(3) - 2._wp*myv_tmp(4))/max(abs(myr_tmp(1)), abs(myr_tmp(4)))
1171 err_v = (-5._wp*h/24._wp)*(mya_tmp(2) + mya_tmp(3) - 2._wp*mya_tmp(4))/max(abs(myv_tmp(1)), abs(myv_tmp(4)))
1172 ! Error correction for non-oscillating bubbles
1173 if (max(abs(myv_tmp(1)), abs(myv_tmp(4))) < 1.e-12_wp) then
1174 err_v = 0._wp
1175 end if
1176 if (bubbles_lagrange .and. f_approx_equal(mya_tmp(1), 0._wp) .and. f_approx_equal(mya_tmp(2), &
1177 & 0._wp) .and. f_approx_equal(mya_tmp(3), 0._wp) .and. f_approx_equal(mya_tmp(4), 0._wp)) then
1178 err_v = 0._wp
1179 end if
1180 err = sqrt((err_r**2._wp + err_v**2._wp)/2._wp)
1181
1182 end subroutine s_advance_substep
1183
1184 !> Changes of pressure and vapor mass in the lagrange bubbles.
1185 subroutine s_advance_el(fR_tmp, fV_tmp, fPb_tmp, fMv_tmp, bub_id, fmass_g, fbeta_c, fbeta_t, fdPbdt_tmp, advance_EL)
1186
1187
1188# 608 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1189#if MFC_OpenACC
1190# 608 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1191!$acc routine seq
1192# 608 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1193#elif MFC_OpenMP
1194# 608 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1195
1196# 608 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1197
1198# 608 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1199!$omp declare target device_type(any)
1200# 608 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles.fpp"
1201#endif
1202 real(wp), intent(in) :: fR_tmp, fV_tmp, fPb_tmp, fMv_tmp
1203 real(wp), intent(in) :: fmass_g, fbeta_c, fbeta_t
1204 integer, intent(in) :: bub_id
1205 real(wp), intent(inout) :: fdPbdt_tmp
1206 real(wp), intent(out) :: advance_EL
1207 real(wp) :: fVapFlux, myR_m, mygamma_m
1208
1209 call s_vflux(fr_tmp, fv_tmp, fpb_tmp, fmv_tmp, bub_id, fvapflux, fmass_g, fbeta_c, myr_m, mygamma_m)
1210 fdpbdt_tmp = f_bpres_dot(fvapflux, fr_tmp, fv_tmp, fpb_tmp, fmv_tmp, bub_id, fbeta_t, myr_m, mygamma_m)
1211 advance_el = 4._wp*pi*fr_tmp**2._wp*fvapflux
1212
1213 end subroutine s_advance_el
1214
1215end module m_bubbles
l
integer, intent(in) l
Definition m_phase_change.fpp.f90:990
m_bubbles_el_kernels
Kernel functions (Gaussian, delta) that smear Lagrangian bubble effects onto the Eulerian grid.
Definition m_bubbles_EL_kernels.fpp.f90:323
m_bubbles_el_kernels::f_get_bubble_force
real(wp) function f_get_bubble_force(pos, rad, rdot, vel, mg, mv, re, rho, cell, i, q_prim_vf)
Definition m_bubbles_EL_kernels.fpp.f90:1484
m_bubbles_el_kernels::f_interpolate_velocity
real(wp) function f_interpolate_velocity(pos, cell, i, q_prim_vf)
Definition m_bubbles_EL_kernels.fpp.f90:1353
m_bubbles
Bubble-dynamics procedures for ensemble- and volume-averaged model.
Definition m_bubbles.fpp.f90:323
m_bubbles::f_rddot_rp
real(wp) function f_rddot_rp(fcp, frho, fr, fv, fcpbw)
Rayleigh-Plesset bubble radial acceleration.
Definition m_bubbles.fpp.f90:552
m_bubbles::s_bwproperty
subroutine s_bwproperty(pb_in, ir0, chi_vw_out, k_mw_out, rho_mw_out)
Compute bubble wall properties for vapor bubbles.
Definition m_bubbles.fpp.f90:684
m_bubbles::s_initial_substep_h
subroutine s_initial_substep_h(frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, fcson, h)
Choose the initial time step size for the adaptive time stepping routine (See Heirer,...
Definition m_bubbles.fpp.f90:984
m_bubbles::f_hdot
real(wp) function f_hdot(fcpbw, fcpinf, fcpinf_dot, fntait, fbtait, fr, fv, fr0, fpbdot)
Enthalpy derivative for Gilmore bubble model, Gilmore (1952).
Definition m_bubbles.fpp.f90:516
m_bubbles::f_advance_step
integer function f_advance_step(frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, bub_id, fmass_v, fmass_g, fbeta_c, fbeta_t, fcson, fre, fpos, fvel, cell, q_prim_vf)
Adaptive time stepping routine for subgrid bubbles (See Heirer, E. Hairer S.P.Norsett G....
Definition m_bubbles.fpp.f90:832
m_bubbles::f_cpbw_km
real(wp) function f_cpbw_km(fr0, fr, fv, fpb)
Keller-Miksis bubble wall pressure.
Definition m_bubbles.fpp.f90:608
m_bubbles::f_rddot_g
real(wp) function f_rddot_g(fcpbw, fr, fv, fh, fhdot, fcgas, fntait, fbtait)
Compute the Gilmore bubble radial acceleration.
Definition m_bubbles.fpp.f90:577
m_bubbles::s_advance_el
subroutine s_advance_el(fr_tmp, fv_tmp, fpb_tmp, fmv_tmp, bub_id, fmass_g, fbeta_c, fbeta_t, fdpbdt_tmp, advance_el)
Changes of pressure and vapor mass in the lagrange bubbles.
Definition m_bubbles.fpp.f90:1186
m_bubbles::f_rddot
real(wp) function f_rddot(frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, fcson)
Compute the bubble radial acceleration based on the selected bubble model.
Definition m_bubbles.fpp.f90:339
m_bubbles::f_cgas
real(wp) function f_cgas(fcpinf, fntait, fbtait, fh)
Compute the sound speed for the bubble.
Definition m_bubbles.fpp.f90:452
m_bubbles::s_vflux
subroutine s_vflux(fr, fv, fpb, fmass_v, ir0, vflux, fmass_g, fbeta_c, fr_m, fgamma_m, fchi_vw, frho_mw)
Compute the vapour flux.
Definition m_bubbles.fpp.f90:719
m_bubbles::f_cpbw
real(wp) function f_cpbw(fr0, fr, fv, fpb)
Bubble wall pressure: stiffened gas with Laplace pressure and viscous stress.
Definition m_bubbles.fpp.f90:393
m_bubbles::f_bpres_dot
real(wp) function f_bpres_dot(fvflux, fr, fv, fpb, fmass_v, ir0, fbeta_t, fr_m, fgamma_m, fk_mw)
Compute the time derivative of the internal bubble pressure.
Definition m_bubbles.fpp.f90:782
m_bubbles::s_advance_substep
subroutine s_advance_substep(err, frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu, bub_id, fmass_v, fmass_g, fbeta_c, fbeta_t, fcson, h, myr_tmp, myv_tmp, mypb_tmp, mymv_tmp)
Integrate bubble variables over the given time step size, h, using a third-order accurate embedded Ru...
Definition m_bubbles.fpp.f90:1064
m_bubbles::f_rddot_km
real(wp) function f_rddot_km(fpbdot, fcp, fcpbw, frho, fr, fv, fr0, fc)
Keller-Miksis bubble radial acceleration.
Definition m_bubbles.fpp.f90:641
m_bubbles::f_h
real(wp) function f_h(fcpbw, fcpinf, fntait, fbtait)
Compute the bubble enthalpy.
Definition m_bubbles.fpp.f90:422
m_bubbles::f_cpinfdot
real(wp) function f_cpinfdot(frho, fp, falf, fntait, fbtait, advsrc, divu)
Compute the time derivative of the driving pressure.
Definition m_bubbles.fpp.f90:482
m_constants
Compile-time constant parameters: default values, tolerances, and physical constants.
Definition m_constants.fpp.f90:7
m_constants::bubble_model_gilmore
integer, parameter bubble_model_gilmore
Definition m_constants.fpp.f90:130
m_constants::scale_first_guess
real(wp), parameter scale_first_guess
Scale factor for initial step size.
Definition m_constants.fpp.f90:80
m_constants::threshold_second_guess
real(wp), parameter threshold_second_guess
Threshold for refined step size estimate.
Definition m_constants.fpp.f90:79
m_constants::pi
real(wp), parameter pi
Pi.
Definition m_constants.fpp.f90:16
m_constants::small_guess
real(wp), parameter small_guess
Minimum initial step size.
Definition m_constants.fpp.f90:82
m_constants::bubble_model_rayleigh_plesset
integer, parameter bubble_model_rayleigh_plesset
Definition m_constants.fpp.f90:132
m_constants::threshold_first_guess
real(wp), parameter threshold_first_guess
Threshold for initial step size estimate.
Definition m_constants.fpp.f90:78
m_constants::bubble_model_keller_miksis
integer, parameter bubble_model_keller_miksis
Definition m_constants.fpp.f90:131
m_constants::verysmall
real(wp), parameter verysmall
Very small number.
Definition m_constants.fpp.f90:17
m_constants::scale_guess
real(wp), parameter scale_guess
Scale factor for step size adjustment.
Definition m_constants.fpp.f90:81
m_derived_types
Shared derived types for field data, patch geometry, bubble dynamics, and MPI I/O structures.
Definition m_derived_types.fpp.f90:323
m_global_parameters
Global parameters for the computational domain, fluid properties, and simulation algorithm configurat...
Definition m_global_parameters.fpp.f90:334
m_global_parameters::gam
real(wp) gam
Definition m_global_parameters.fpp.f90:973
m_global_parameters::r_v
real(wp) r_v
Definition m_global_parameters.fpp.f90:986
m_global_parameters::pb0
real(wp), dimension(:), allocatable pb0
Definition m_global_parameters.fpp.f90:948
m_global_parameters::re_trans_t
real(wp), dimension(:), allocatable re_trans_t
Definition m_global_parameters.fpp.f90:949
m_global_parameters::lag_vel_model
integer lag_vel_model
Definition m_global_parameters.fpp.f90:1045
m_global_parameters::gam_m
real(wp) gam_m
Definition m_global_parameters.fpp.f90:973
m_global_parameters::m_v
real(wp) m_v
Definition m_global_parameters.fpp.f90:986
m_global_parameters::k_v
real(wp), dimension(:), allocatable k_v
Definition m_global_parameters.fpp.f90:948
m_global_parameters::r_g
real(wp) r_g
Definition m_global_parameters.fpp.f90:986
m_global_parameters::phi_vg
real(wp) phi_vg
Definition m_global_parameters.fpp.f90:935
m_global_parameters::r0
real(wp), dimension(:), allocatable r0
Bubble sizes.
Definition m_global_parameters.fpp.f90:861
m_global_parameters::k_g
real(wp), dimension(:), allocatable k_g
Definition m_global_parameters.fpp.f90:948
m_global_parameters::gam_v
real(wp) gam_v
Definition m_global_parameters.fpp.f90:986
m_global_parameters::phi_gv
real(wp) phi_gv
Definition m_global_parameters.fpp.f90:935
m_global_parameters::pe_c
real(wp) pe_c
Definition m_global_parameters.fpp.f90:935
m_global_parameters::m_g
real(wp) m_g
Definition m_global_parameters.fpp.f90:986
m_global_parameters::gam_g
real(wp) gam_g
Definition m_global_parameters.fpp.f90:986
m_global_parameters::mass_g0
real(wp), dimension(:), allocatable mass_g0
Definition m_global_parameters.fpp.f90:948
m_global_parameters::tw
real(wp) tw
Definition m_global_parameters.fpp.f90:935
m_global_parameters::re_trans_c
real(wp), dimension(:), allocatable re_trans_c
Definition m_global_parameters.fpp.f90:949
m_global_parameters::pv
real(wp) pv
Definition m_global_parameters.fpp.f90:986
m_global_parameters::pe_t
real(wp), dimension(:), allocatable pe_t
Definition m_global_parameters.fpp.f90:948
m_global_parameters::mass_v0
real(wp), dimension(:), allocatable mass_v0
Definition m_global_parameters.fpp.f90:948
m_global_parameters::eu
real(wp) eu
Euler number.
Definition m_global_parameters.fpp.f90:847
m_helper_basic
Basic floating-point utilities: approximate equality, default detection, and coordinate bounds.
Definition m_helper_basic.fpp.f90:323
m_helper_basic::f_approx_equal
logical elemental function, public f_approx_equal(a, b, tol_input)
Check if two floating point numbers of wp are within tolerance.
Definition m_helper_basic.fpp.f90:340
m_helper_basic::f_is_default
logical elemental function, public f_is_default(var)
Check if a real(wp) variable is of default value.
Definition m_helper_basic.fpp.f90:428
m_mpi_proxy
MPI halo exchange, domain decomposition, and buffer packing/unpacking for the simulation solver.
Definition m_mpi_proxy.fpp.f90:334
m_variables_conversion
Conservative-to-primitive variable conversion, mixture property evaluation, and pressure computation.
Definition m_variables_conversion.fpp.f90:334
m_derived_types::scalar_field
Derived type annexing a scalar field (SF).
Definition m_derived_types.fpp.f90:337