#include "cs_base.h"
#include "cs_domain.h"

Include dependency graph for cs_vof.h:

Data Structures
struct	cs_vof_parameters_t
	VOF model parameters. Void fraction variable tracks fluid 2. More...

struct	cs_cavitation_parameters_t
	Cavitation model parameters. More...

Macros
#define	CS_VOF_ENABLED (1 << 0)

#define	CS_VOF_FREE_SURFACE (1 << 1)

#define	CS_VOF_MERKLE_MASS_TRANSFER (1 << 2)

Functions
cs_vof_parameters_t *	cs_get_glob_vof_parameters (void)

void	cs_vof_compute_linear_rho_mu (const cs_domain_t *domain)
	Compute the mixture density, mixture dynamic viscosity given fluid volume fractions and the reference density and dynamic viscosity $\rho_l, \mu_l$ (liquid), $\rho_v, \mu_v$ (gas).

void	cs_vof_update_phys_prop (const cs_domain_t *domain)
	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).

void	cs_vof_log_mass_budget (const cs_domain_t *domain)
	Write in main log the global mixture mass budget:

cs_cavitation_parameters_t *	cs_get_glob_cavitation_parameters (void)

Variables
const cs_vof_parameters_t *	cs_glob_vof_parameters

Function Documentation

◆ cs_get_glob_cavitation_parameters()

cs_cavitation_parameters_t * cs_get_glob_cavitation_parameters ( void )

◆ cs_get_glob_vof_parameters()

cs_vof_parameters_t * cs_get_glob_vof_parameters ( void )

◆ cs_vof_compute_linear_rho_mu()

void cs_vof_compute_linear_rho_mu ( const cs_domain_t * domain )

Compute the mixture density, mixture dynamic viscosity given fluid volume fractions and the reference density and dynamic viscosity $\rho_l, \mu_l$ (liquid), $\rho_v, \mu_v$ (gas).

Computation is done as follows on cells:

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

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

A similar linear formula is followed on boundary using fluid volume fraction value on the boundary.

Computation is done as follows on cells:

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

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

A similar linear formula is followed on boundary using fluid volume fraction value on the boundary.

◆ cs_vof_log_mass_budget()

void cs_vof_log_mass_budget ( const cs_domain_t * domain )

Write in main log 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).$

◆ cs_vof_update_phys_prop()

void cs_vof_update_phys_prop ( const cs_domain_t * domain )

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).

For the computation of mixture density, mixture dynamic viscosity, see cs_vof_compute_linear_rho_mu.

Computation of mass flux is as follows:

$\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.$

For the computation of mixture density, mixture dynamic viscosity, see cs_vof_compute_linear_rho_mu.

Computation of mass flux is as follows:

$\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.$

Variable Documentation

◆ cs_glob_vof_parameters

const cs_vof_parameters_t* cs_glob_vof_parameters

extern

Data Structures

Macros

Functions

Variables