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]	crom	Density at cell centers at current time step
[in]	croma	Density at cell centers at previous time step
[in]	brom	Density at boundary faces at current time step
[in]	dt	Time step
[in]	imasfl_rel	Mass flux at internal faces array
[in]	bmasfl_rel	Mass 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]	voidf	Void fraction array
[in]	coavoi	Void fraction boundary coefficient array
[in]	cobvoi	Void fraction boundary coefficient array
[in]	ivoifl	Volumetric flux at internal faces array
[in]	bvoifl	Volumetric flux at boundary faces array
[out]	crom	Density at cell center array
[out]	brom	Density at boudary faces array
[out]	imasfl	Mass flux at internal faces array
[out]	bmasfl	Mass flux at internal boundary faces array

Functions/Subroutines

Variables

Function/Subroutine Documentation

◆ init_vof()

◆ vof_print_mass_budget()

◆ vof_update_phys_prop()