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