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

Bubble-dynamics procedures 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

Bubble-dynamics procedures 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 787 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 458 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 399 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 614 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 488 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 428 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 522 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 345 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 583 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 647 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 558 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 1182 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 835 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 1058 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 690 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 979 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 725 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 327 of file m_bubbles.fpp.f90.

◆ k_mw

real(wp) m_bubbles::k_mw

Bubble wall properties (Ando 2010).

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

◆ rho_mw

real(wp) m_bubbles::rho_mw

Bubble wall properties (Ando 2010).

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