MFC: Simulation
High-fidelity multiphase flow simulation
Loading...
Searching...
No Matches
m_bubbles.fpp.f90 File Reference

Functions/Subroutines

program __m_bubbles_fpp_f90__
 
Bubble dynamic source terms
subroutine s_initialize_bubbles_module ()
 
subroutine s_comp_alpha_from_n (q_cons_vf)
 
subroutine s_compute_bubbles_rhs (idir, q_prim_vf)
 
subroutine s_compute_bubble_source (q_cons_vf, q_prim_vf, t_step, rhs_vf)
 The purpose of this procedure is to compute the source terms that are needed for the bubble modeling.
 
real(kind(0d0)) function f_cpbw (fr0, fr, fv, fpb)
 Function that computes that bubble wall pressure for Gilmore bubbles.
 
real(kind(0d0)) function f_h (fcpbw, fcpinf, fntait, fbtait)
 Function that computes the bubble enthalpy.
 
real(kind(0d0)) function f_cgas (fcpinf, fntait, fbtait, fh)
 Function that computes the sound speed for the bubble.
 
real(kind(0d0)) function f_cpinfdot (frho, fp, falf, fntait, fbtait, advsrc, divu)
 Function that computes the time derivative of the driving pressure.
 
real(kind(0d0)) function f_hdot (fcpbw, fcpinf, fcpinf_dot, fntait, fbtait, fr, fv, fr0, fpbdot)
 Function that computes the time derivative of the enthalpy.
 
real(kind(0d0)) function f_rddot (frho, fp, fr, fv, fr0, fpb, fpbdot, alf, fntait, fbtait, f_bub_adv_src, f_divu)
 Function that computes the bubble radial acceleration based on bubble models.
 
real(kind(0d0)) function f_rddot_rp (fcp, frho, fr, fv, fr0, fcpbw)
 Function that computes the bubble radial acceleration for Rayleigh-Plesset bubbles.
 
real(kind(0d0)) function f_rddot_g (fcpbw, fr, fv, fh, fhdot, fcgas, fntait, fbtait)
 Function that computes the bubble radial acceleration.
 
real(kind(0d0)) function f_cpbw_km (fr0, fr, fv, fpb)
 Function that computes the bubble wall pressure for Keller–Miksis bubbles.
 
real(kind(0d0)) function f_rddot_km (fpbdot, fcp, fcpbw, frho, fr, fv, fr0, fc)
 Function that computes the bubble radial acceleration for Keller–Miksis bubbles.
 
subroutine s_bwproperty (pb, ir0)
 Subroutine that computes bubble wall properties for vapor bubbles.
 
real(kind(0.d0)) function f_vflux (fr, fv, fmass_v, ir0)
 Function that computes the vapour flux.
 
real(kind(0.d0)) function f_bpres_dot (fvflux, fr, fv, fpb, fmass_v, ir0)
 Function that computes the time derivative of the internal bubble pressure.
 

Function/Subroutine Documentation

◆ __m_bubbles_fpp_f90__()

program __m_bubbles_fpp_f90__

◆ f_bpres_dot()

real(kind(0.d0)) function __m_bubbles_fpp_f90__::f_bpres_dot ( real(kind(0.d0)), intent(in) fvflux,
real(kind(0.d0)), intent(in) fr,
real(kind(0.d0)), intent(in) fv,
real(kind(0.d0)), intent(in) fpb,
real(kind(0.d0)), intent(in) fmass_v,
integer, intent(in) ir0 )
private

Function that computes the time derivative of the internal bubble pressure.

Parameters
fvfluxVapour flux
fRCurrent bubble radius
fVCurrent bubble velocity
fpbCurrent internal bubble pressure
fmass_vCurrent mass of vapour
iR0Bubble size index
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cgas()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_cgas ( real(kind(0d0)), intent(in) fcpinf,
real(kind(0d0)), intent(in) fntait,
real(kind(0d0)), intent(in) fbtait,
real(kind(0d0)), intent(in) fh )
private

Function that computes the sound speed for the bubble.

Parameters
fCpinfDriving bubble pressure
fntaitTait EOS parameter
fBtaitTait EOS parameter
fHBubble enthalpy
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cpbw()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_cpbw ( real(kind(0d0)), intent(in) fr0,
real(kind(0d0)), intent(in) fr,
real(kind(0d0)), intent(in) fv,
real(kind(0d0)), intent(in) fpb )
private

Function that computes that bubble wall pressure for Gilmore bubbles.

Parameters
fR0Equilibrium bubble radius
fRCurrent bubble radius
fVCurrent bubble velocity
fpbInternal bubble pressure
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cpbw_km()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_cpbw_km ( real(kind(0d0)), intent(in) fr0,
real(kind(0d0)), intent(in) fr,
real(kind(0d0)), intent(in) fv,
real(kind(0d0)), intent(in) fpb )
private

Function that computes the bubble wall pressure for Keller–Miksis bubbles.

Parameters
fR0Equilibrium bubble radius
fRCurrent bubble radius
fVCurrent bubble velocity
fpbInternal bubble pressure
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_cpinfdot()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_cpinfdot ( real(kind(0d0)), intent(in) frho,
real(kind(0d0)), intent(in) fp,
real(kind(0d0)), intent(in) falf,
real(kind(0d0)), intent(in) fntait,
real(kind(0d0)), intent(in) fbtait,
real(kind(0d0)), intent(in) advsrc,
real(kind(0d0)), intent(in) divu )
private

Function that computes the time derivative of the driving pressure.

Parameters
fRhoLocal liquid density
fPLocal pressure
falfLocal void fraction
fntaitTait EOS parameter
fBtaitTait EOS parameter
advsrcAdvection equation source term
divuDivergence of velocity
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_h()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_h ( real(kind(0d0)), intent(in) fcpbw,
real(kind(0d0)), intent(in) fcpinf,
real(kind(0d0)), intent(in) fntait,
real(kind(0d0)), intent(in) fbtait )
private

Function that computes the bubble enthalpy.

Parameters
fCpbwBubble wall pressure
fCpinfDriving bubble pressure
fntaitTait EOS parameter
fBtaitTait EOS parameter
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_hdot()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_hdot ( real(kind(0d0)), intent(in) fcpbw,
real(kind(0d0)), intent(in) fcpinf,
real(kind(0d0)), intent(in) fcpinf_dot,
real(kind(0d0)), intent(in) fntait,
real(kind(0d0)), intent(in) fbtait,
real(kind(0d0)), intent(in) fr,
real(kind(0d0)), intent(in) fv,
real(kind(0d0)), intent(in) fr0,
real(kind(0d0)), intent(in) fpbdot )
private

Function that computes the time derivative of the enthalpy.

Parameters
fCpbwBubble wall pressure
fCpinfDriving bubble pressure
fCpinf_dotTime derivative of the driving pressure
fntaitTait EOS parameter
fBtaitTait EOS parameter
fR0Equilibrium bubble radius
fRCurrent bubble radius
fVCurrent bubble velocity
fpbdotTime derivative of the internal bubble pressure
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_rddot ( real(kind(0d0)), intent(in) frho,
real(kind(0d0)), intent(in) fp,
real(kind(0d0)), intent(in) fr,
real(kind(0d0)), intent(in) fv,
real(kind(0d0)), intent(in) fr0,
real(kind(0d0)), intent(in) fpb,
real(kind(0d0)), intent(in) fpbdot,
real(kind(0d0)), intent(in) alf,
real(kind(0d0)), intent(in) fntait,
real(kind(0d0)), intent(in) fbtait,
real(kind(0d0)), intent(in) f_bub_adv_src,
real(kind(0d0)), intent(in) f_divu )
private

Function that computes the bubble radial acceleration based on bubble models.

Parameters
fRhoCurrent density
fPCurrent driving pressure
fRCurrent bubble radius
fVCurrent bubble velocity
fR0Equilibrium bubble radius
fpbInternal bubble pressure
fpbdotTime-derivative of internal bubble pressure
alfbubble volume fraction
fntaitTait EOS parameter
fBtaitTait EOS parameter
f_bub_adv_srcSource for bubble volume fraction
f_divuDivergence of velocity
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot_g()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_rddot_g ( real(kind(0d0)), intent(in) fcpbw,
real(kind(0d0)), intent(in) fr,
real(kind(0d0)), intent(in) fv,
real(kind(0d0)), intent(in) fh,
real(kind(0d0)), intent(in) fhdot,
real(kind(0d0)), intent(in) fcgas,
real(kind(0d0)), intent(in) fntait,
real(kind(0d0)), intent(in) fbtait )
private

Function that computes the bubble radial acceleration.

Parameters
fCpbwBubble wall pressure
fRCurrent bubble radius
fVCurrent bubble velocity
fHCurrent enthalpy
fHdotCurrent time derivative of the enthalpy
fcgasCurrent gas sound speed
fntaitTait EOS parameter
fBtaitTait EOS parameter
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot_km()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_rddot_km ( real(kind(0d0)), intent(in) fpbdot,
real(kind(0d0)), intent(in) fcp,
real(kind(0d0)), intent(in) fcpbw,
real(kind(0d0)), intent(in) frho,
real(kind(0d0)), intent(in) fr,
real(kind(0d0)), intent(in) fv,
real(kind(0d0)), intent(in) fr0,
real(kind(0d0)), intent(in) fc )
private

Function that computes the bubble radial acceleration for Keller–Miksis bubbles.

Parameters
fpbdotTime-derivative of internal bubble pressure
fCpDriving pressure
fCpbwBubble wall pressure
fRhoCurrent density
fRCurrent bubble radius
fVCurrent bubble velocity
fR0Equilibrium bubble radius
fCCurrent sound speed
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_rddot_rp()

real(kind(0d0)) function __m_bubbles_fpp_f90__::f_rddot_rp ( real(kind(0d0)), intent(in) fcp,
real(kind(0d0)), intent(in) frho,
real(kind(0d0)), intent(in) fr,
real(kind(0d0)), intent(in) fv,
real(kind(0d0)), intent(in) fr0,
real(kind(0d0)), intent(in) fcpbw )
private

Function that computes the bubble radial acceleration for Rayleigh-Plesset bubbles.

Parameters
fCpDriving pressure
fRhoCurrent density
fRCurrent bubble radius
fVCurrent bubble velocity
fR0Equilibrium bubble radius
fCpbwBoundary wall pressure
Here is the call graph for this function:
Here is the caller graph for this function:

◆ f_vflux()

real(kind(0.d0)) function __m_bubbles_fpp_f90__::f_vflux ( real(kind(0.d0)), intent(in) fr,
real(kind(0.d0)), intent(in) fv,
real(kind(0.d0)), intent(in) fmass_v,
integer, intent(in) ir0 )
private

Function that computes the vapour flux.

Parameters
fRCurrent bubble radius
fVCurrent bubble velocity
fmass_vCurrent mass of vapour
iR0Bubble size index
Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_bwproperty()

subroutine __m_bubbles_fpp_f90__::s_bwproperty ( real(kind(0.d0)), intent(in) pb,
integer, intent(in) ir0 )
private

Subroutine that computes bubble wall properties for vapor bubbles.

Parameters
pbInternal bubble pressure
iR0Current bubble size index
Here is the caller graph for this function:

◆ s_comp_alpha_from_n()

subroutine __m_bubbles_fpp_f90__::s_comp_alpha_from_n ( type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf)
private
Here is the caller graph for this function:

◆ s_compute_bubble_source()

subroutine __m_bubbles_fpp_f90__::s_compute_bubble_source ( type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf,
type(scalar_field), dimension(sys_size), intent(in) q_prim_vf,
integer, intent(in) t_step,
type(scalar_field), dimension(sys_size), intent(inout) rhs_vf )
private

The purpose of this procedure is to compute the source terms that are needed for the bubble modeling.

Parameters
q_prim_vfPrimitive variables
q_cons_vfConservative variables
divuDivergence of velocity
bub_adv_srcAdvection equation source due to bubble compression/expansion
bub_r_srcBubble radius equation source
bub_v_srcBubble velocity equation source
bub_p_srcBubble pressure equation source
bub_m_srcBubble mass equation source
Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_compute_bubbles_rhs()

subroutine __m_bubbles_fpp_f90__::s_compute_bubbles_rhs ( integer idir,
type(scalar_field), dimension(sys_size), intent(in) q_prim_vf )
private
Here is the caller graph for this function:

◆ s_initialize_bubbles_module()

subroutine __m_bubbles_fpp_f90__::s_initialize_bubbles_module
private
Here is the caller graph for this function: