MFC
Exascale flow solver
Loading...
Searching...
No Matches
m_bubbles Module Reference

Shared bubble-dynamics procedures (radial acceleration, wall pressure, sound speed) for ensemble- and volume-averaged models. More...

Functions/Subroutines

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.
elemental real(wp) function f_cpbw (fr0, fr, fv, fpb)
 Bubble wall pressure: stiffened gas with Laplace pressure and viscous stress.
elemental real(wp) function f_h (fcpbw, fcpinf, fntait, fbtait)
 Compute the bubble enthalpy.
elemental real(wp) function f_cgas (fcpinf, fntait, fbtait, fh)
 Compute the sound speed for the bubble.
elemental real(wp) function f_cpinfdot (frho, fp, falf, fntait, fbtait, advsrc, divu)
 Compute the time derivative of the driving pressure.
elemental real(wp) function f_hdot (fcpbw, fcpinf, fcpinf_dot, fntait, fbtait, fr, fv, fr0, fpbdot)
 Enthalpy derivative for Gilmore bubble model, Gilmore (1952).
elemental real(wp) function f_rddot_rp (fcp, frho, fr, fv, fcpbw)
 Rayleigh-Plesset bubble radial acceleration.
elemental real(wp) function f_rddot_g (fcpbw, fr, fv, fh, fhdot, fcgas, fntait, fbtait)
 Compute the Gilmore bubble radial acceleration.
elemental real(wp) function f_cpbw_km (fr0, fr, fv, fpb)
 Keller-Miksis bubble wall pressure.
elemental real(wp) function f_rddot_km (fpbdot, fcp, fcpbw, frho, fr, fv, fr0, fc)
 Keller-Miksis bubble radial acceleration.
elemental subroutine s_bwproperty (pb_in, ir0, chi_vw_out, k_mw_out, rho_mw_out)
 Compute bubble wall properties for vapor bubbles.
elemental subroutine s_vflux (fr, fv, fpb, fmass_v, ir0, vflux, fmass_g, fbeta_c, fr_m, fgamma_m)
 Compute the vapour flux.
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.
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. Wanner, Solving Ordinary Differential Equations I, Chapter II.4).
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, E. Hairer S.P.Norsett G. Wanner, Solving Ordinary Differential Equations I, Chapter II.4).
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 Runge-Kutta scheme.
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.

Variables

real(wp) chi_vw
 Bubble wall properties (Ando 2010).
real(wp) k_mw
 Bubble wall properties (Ando 2010).
real(wp) rho_mw
 Bubble wall properties (Ando 2010).

Detailed Description

Shared bubble-dynamics procedures (radial acceleration, wall pressure, sound speed) for ensemble- and volume-averaged models.

Function/Subroutine Documentation

◆ f_bpres_dot()

elemental real(wp) function m_bubbles::f_bpres_dot ( real(wp), intent(in) fvflux,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fpb,
real(wp), intent(in) fmass_v,
integer, intent(in) ir0,
real(wp), intent(in), optional fbeta_t,
real(wp), intent(in), optional fr_m,
real(wp), intent(in), optional fgamma_m )

Compute the time derivative of the internal bubble pressure.

Definition at line 788 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cgas()

elemental real(wp) function m_bubbles::f_cgas ( real(wp), intent(in) fcpinf,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait,
real(wp), intent(in) fh )

Compute the sound speed for the bubble.

Definition at line 459 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cpbw()

elemental real(wp) function m_bubbles::f_cpbw ( real(wp), intent(in) fr0,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fpb )

Bubble wall pressure: stiffened gas with Laplace pressure and viscous stress.

Definition at line 400 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cpbw_km()

elemental real(wp) function m_bubbles::f_cpbw_km ( real(wp), intent(in) fr0,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fpb )

Keller-Miksis bubble wall pressure.

Definition at line 615 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cpinfdot()

elemental real(wp) function m_bubbles::f_cpinfdot ( real(wp), intent(in) frho,
real(wp), intent(in) fp,
real(wp), intent(in) falf,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait,
real(wp), intent(in) advsrc,
real(wp), intent(in) divu )

Compute the time derivative of the driving pressure.

Definition at line 489 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_h()

elemental real(wp) function m_bubbles::f_h ( real(wp), intent(in) fcpbw,
real(wp), intent(in) fcpinf,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait )

Compute the bubble enthalpy.

Definition at line 429 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_hdot()

elemental real(wp) function m_bubbles::f_hdot ( real(wp), intent(in) fcpbw,
real(wp), intent(in) fcpinf,
real(wp), intent(in) fcpinf_dot,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fr0,
real(wp), intent(in) fpbdot )

Enthalpy derivative for Gilmore bubble model, Gilmore (1952).

Definition at line 523 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot()

elemental real(wp) function m_bubbles::f_rddot ( real(wp), intent(in) frho,
real(wp), intent(in) fp,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fr0,
real(wp), intent(in) fpb,
real(wp), intent(in) fpbdot,
real(wp), intent(in) alf,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait,
real(wp), intent(in) f_bub_adv_src,
real(wp), intent(in) f_divu,
real(wp), intent(in) fcson )

Compute the bubble radial acceleration based on the selected bubble model.

Definition at line 346 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot_g()

elemental real(wp) function m_bubbles::f_rddot_g ( real(wp), intent(in) fcpbw,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fh,
real(wp), intent(in) fhdot,
real(wp), intent(in) fcgas,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait )

Compute the Gilmore bubble radial acceleration.

Definition at line 584 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot_km()

elemental real(wp) function m_bubbles::f_rddot_km ( real(wp), intent(in) fpbdot,
real(wp), intent(in) fcp,
real(wp), intent(in) fcpbw,
real(wp), intent(in) frho,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fr0,
real(wp), intent(in) fc )

Keller-Miksis bubble radial acceleration.

Definition at line 648 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot_rp()

elemental real(wp) function m_bubbles::f_rddot_rp ( real(wp), intent(in) fcp,
real(wp), intent(in) frho,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fcpbw )

Rayleigh-Plesset bubble radial acceleration.

Definition at line 559 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_advance_el()

elemental subroutine m_bubbles::s_advance_el ( real(wp), intent(in) fr_tmp,
real(wp), intent(in) fv_tmp,
real(wp), intent(in) fpb_tmp,
real(wp), intent(in) fmv_tmp,
integer, intent(in) bub_id,
real(wp), intent(in) fmass_g,
real(wp), intent(in) fbeta_c,
real(wp), intent(in) fbeta_t,
real(wp), intent(inout) fdpbdt_tmp,
real(wp), intent(out) advance_el )

Changes of pressure and vapor mass in the lagrange bubbles.

Definition at line 1183 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_advance_step()

subroutine m_bubbles::s_advance_step ( real(wp), intent(in) frho,
real(wp), intent(in) fp,
real(wp), intent(inout) fr,
real(wp), intent(inout) fv,
real(wp), intent(in) fr0,
real(wp), intent(inout) fpb,
real(wp), intent(in) fpbdot,
real(wp), intent(in) alf,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait,
real(wp), intent(in) f_bub_adv_src,
real(wp), intent(in) f_divu,
integer, intent(in) bub_id,
real(wp), intent(inout) fmass_v,
real(wp), intent(in) fmass_g,
real(wp), intent(in) fbeta_c,
real(wp), intent(in) fbeta_t,
real(wp), intent(in) fcson,
integer, intent(inout) adap_dt_stop )

Adaptive time stepping routine for subgrid bubbles (See Heirer, E. Hairer S.P.Norsett G. Wanner, Solving Ordinary Differential Equations I, Chapter II.4).

Bubble radius, radial velocity, and radial acceleration for the inner loop

Definition at line 836 of file m_bubbles.fpp.f90.

Here is the call graph for this function:

◆ s_advance_substep()

subroutine m_bubbles::s_advance_substep ( real(wp), intent(out) err,
real(wp), intent(in) frho,
real(wp), intent(in) fp,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fr0,
real(wp), intent(in) fpb,
real(wp), intent(in) fpbdot,
real(wp), intent(in) alf,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait,
real(wp), intent(in) f_bub_adv_src,
real(wp), intent(in) f_divu,
integer, intent(in) bub_id,
real(wp), intent(in) fmass_v,
real(wp), intent(in) fmass_g,
real(wp), intent(in) fbeta_c,
real(wp), intent(in) fbeta_t,
real(wp), intent(in) fcson,
real(wp), intent(in) h,
real(wp), dimension(4), intent(out) myr_tmp,
real(wp), dimension(4), intent(out) myv_tmp,
real(wp), dimension(4), intent(out) mypb_tmp,
real(wp), dimension(4), intent(out) mymv_tmp )

Integrate bubble variables over the given time step size, h, using a third-order accurate embedded Runge-Kutta scheme.

Definition at line 1059 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_bwproperty()

elemental subroutine m_bubbles::s_bwproperty ( real(wp), intent(in) pb_in,
integer, intent(in) ir0,
real(wp), intent(out) chi_vw_out,
real(wp), intent(out) k_mw_out,
real(wp), intent(out) rho_mw_out )

Compute bubble wall properties for vapor bubbles.

Parameters
[out]chi_vw_outBubble wall properties (Ando 2010)
[out]k_mw_outBubble wall properties (Ando 2010)
[out]rho_mw_outBubble wall properties (Ando 2010)

Definition at line 691 of file m_bubbles.fpp.f90.

◆ s_initial_substep_h()

subroutine m_bubbles::s_initial_substep_h ( real(wp), intent(in) frho,
real(wp), intent(in) fp,
real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fr0,
real(wp), intent(in) fpb,
real(wp), intent(in) fpbdot,
real(wp), intent(in) alf,
real(wp), intent(in) fntait,
real(wp), intent(in) fbtait,
real(wp), intent(in) f_bub_adv_src,
real(wp), intent(in) f_divu,
real(wp), intent(in) fcson,
real(wp), intent(out) h )

Choose the initial time step size for the adaptive time stepping routine (See Heirer, E. Hairer S.P.Norsett G. Wanner, Solving Ordinary Differential Equations I, Chapter II.4).

Definition at line 980 of file m_bubbles.fpp.f90.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_vflux()

elemental subroutine m_bubbles::s_vflux ( real(wp), intent(in) fr,
real(wp), intent(in) fv,
real(wp), intent(in) fpb,
real(wp), intent(in) fmass_v,
integer, intent(in) ir0,
real(wp), intent(out) vflux,
real(wp), intent(in), optional fmass_g,
real(wp), intent(in), optional fbeta_c,
real(wp), intent(out), optional fr_m,
real(wp), intent(out), optional fgamma_m )

Compute the vapour flux.

Definition at line 726 of file m_bubbles.fpp.f90.

Here is the caller graph for this function:

Variable Documentation

◆ chi_vw

real(wp) m_bubbles::chi_vw

Bubble wall properties (Ando 2010).

Definition at line 328 of file m_bubbles.fpp.f90.

◆ k_mw

real(wp) m_bubbles::k_mw

Bubble wall properties (Ando 2010).

Definition at line 329 of file m_bubbles.fpp.f90.

◆ rho_mw

real(wp) m_bubbles::rho_mw

Bubble wall properties (Ando 2010).

Definition at line 330 of file m_bubbles.fpp.f90.