programmer's documentation
Functions/Subroutines | Variables
vof Module Reference

Functions/Subroutines

subroutine init_vof
 Default initialization of the module variables. More...
 
subroutine vof_update_phys_prop (voidf, coavoi, cobvoi, ivoifl, bvoifl, crom, brom, imasfl, bmasfl)
 Compute the mixture density, mixture dynamic viscosity and mixture mass flux given the volumetric flux, the volume fraction and the reference density and dynamic viscosity $ \rho_l, \mu_l $ (liquid), $ \rho_v, \mu_v $ (gas) as follows:

\[ \rho_\celli = \alpha_\celli \rho_v + (1-\alpha_\celli) \rho_l, \]

\[ \mu_\celli = \alpha_\celli \mu_v + (1-\alpha_\celli) \mu_l, \]

\[ \left( \rho\vect{u}\cdot\vect{S} \right)_\ij = \\ \left\lbrace \begin{array}{ll} \rho_\celli (\vect{u}\cdot\vect{S})_\ij &\text{ if } (\vect{u}\cdot\vect{S})_\ij>0, \\ \rho_\cellj (\vect{u}\cdot\vect{S})_\ij &\text{ otherwise }, \end{array} \right. \]

\[ \left( \rho\vect{u}\cdot\vect{S} \right)_\ib = \\ \left\lbrace \begin{array}{ll} \rho_\celli (\vect{u}\cdot\vect{S})_\ib &\text{ if } (\vect{u}\cdot\vect{S})_\ib>0, \\ \rho_b (\vect{u}\cdot\vect{S})_\ib &\text{ otherwise }. \end{array} \right. \]

. More...

 
subroutine vof_print_mass_budget (crom, croma, brom, dt, imasfl_rel, bmasfl_rel)
 Print the global mixture mass budget:

\[ \sum_i\left( |\Omega_i|\dfrac{\alpha_i^n - \alpha_i^{n-1}}{\Delta t} + \sum_{j\in\Face{\celli}}\left(\rho\vect{u}\vect{S}\right)_{ij}^n \right). \]

. More...

 

Variables

double precision, save rho1
 reference density of fluid 1 (kg/m3). By convention, liquid phase for cavitation model. More...
 
double precision, save rho2
 reference density of fluid 2 (kg/m3). By convention, gas phase for cavitation model. More...
 
double precision, save mu1
 reference molecular viscosity of fluid 1 (kg/(m s)) More...
 
double precision, save mu2
 reference molecular viscosity of fluid 2 (kg/(m s)) More...
 
double precision, save clvfmn
 clipping min. for the volume fraction More...
 
double precision, save clvfmx
 clipping max. for the volume fraction More...
 

Function/Subroutine Documentation

◆ init_vof()

subroutine vof::init_vof ( )

Default initialization of the module variables.

◆ vof_print_mass_budget()

subroutine vof::vof_print_mass_budget ( double precision, dimension(ncelet)  crom,
double precision, dimension(ncelet)  croma,
double precision, dimension(nfabor)  brom,
double precision, dimension(ncelet)  dt,
double precision, dimension(:), target  imasfl_rel,
double precision, dimension(:), target  bmasfl_rel 
)

Print the global mixture mass budget:

\[ \sum_i\left( |\Omega_i|\dfrac{\alpha_i^n - \alpha_i^{n-1}}{\Delta t} + \sum_{j\in\Face{\celli}}\left(\rho\vect{u}\vect{S}\right)_{ij}^n \right). \]

.

Parameters
[in]cromDensity at cell centers at current time step
[in]cromaDensity at cell centers at previous time step
[in]bromDensity at boundary faces at current time step
[in]dtTime step
[in]imasfl_relMass flux at internal faces array
[in]bmasfl_relMass flux at internal boundary faces array

◆ vof_update_phys_prop()

subroutine vof::vof_update_phys_prop ( double precision, dimension(ncelet)  voidf,
double precision, dimension(nfabor)  coavoi,
double precision, dimension(nfabor)  cobvoi,
double precision, dimension(nfac)  ivoifl,
double precision, dimension(nfabor)  bvoifl,
double precision, dimension(ncelet)  crom,
double precision, dimension(nfabor)  brom,
double precision, dimension(nfac)  imasfl,
double precision, dimension(nfabor)  bmasfl 
)

Compute the mixture density, mixture dynamic viscosity and mixture mass flux given the volumetric flux, the volume fraction and the reference density and dynamic viscosity $ \rho_l, \mu_l $ (liquid), $ \rho_v, \mu_v $ (gas) as follows:

\[ \rho_\celli = \alpha_\celli \rho_v + (1-\alpha_\celli) \rho_l, \]

\[ \mu_\celli = \alpha_\celli \mu_v + (1-\alpha_\celli) \mu_l, \]

\[ \left( \rho\vect{u}\cdot\vect{S} \right)_\ij = \\ \left\lbrace \begin{array}{ll} \rho_\celli (\vect{u}\cdot\vect{S})_\ij &\text{ if } (\vect{u}\cdot\vect{S})_\ij>0, \\ \rho_\cellj (\vect{u}\cdot\vect{S})_\ij &\text{ otherwise }, \end{array} \right. \]

\[ \left( \rho\vect{u}\cdot\vect{S} \right)_\ib = \\ \left\lbrace \begin{array}{ll} \rho_\celli (\vect{u}\cdot\vect{S})_\ib &\text{ if } (\vect{u}\cdot\vect{S})_\ib>0, \\ \rho_b (\vect{u}\cdot\vect{S})_\ib &\text{ otherwise }. \end{array} \right. \]

.

Parameters
[in]voidfVoid fraction array
[in]coavoiVoid fraction boundary coefficient array
[in]cobvoiVoid fraction boundary coefficient array
[in]ivoiflVolumetric flux at internal faces array
[in]bvoiflVolumetric flux at boundary faces array
[out]cromDensity at cell center array
[out]bromDensity at boudary faces array
[out]imasflMass flux at internal faces array
[out]bmasflMass flux at internal boundary faces array