Go to the source code of this file.
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_mesh_t *m) |
Compute the mixture density, mixture dynamic viscosity given fluid volume fractions and the reference density and dynamic viscosity ![]() ![]() | |
void | cs_vof_update_phys_prop (const cs_mesh_t *m) |
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 ![]() ![]() | |
void | cs_vof_surface_tension (const cs_mesh_t *m, const cs_mesh_quantities_t *mq, cs_real_3_t stf[]) |
Compute the surface tension momentum source term following the CSF model of Brackbill et al. (1992). More... | |
void | cs_vof_log_mass_budget (const cs_mesh_t *m, const cs_mesh_quantities_t *mq) |
Write in main log the global mixture mass budget: More... | |
void | cs_vof_deshpande_drift_flux (const cs_mesh_t *m, const cs_mesh_quantities_t *mq) |
Compute the flux of the drift velocity ![]() ![]() | |
void | cs_vof_drift_term (int imrgra, int nswrgp, int imligp, int iwarnp, cs_real_t epsrgp, cs_real_t climgp, cs_real_t *restrict pvar, const cs_real_t *restrict pvara, cs_real_t *restrict rhs) |
Add the explicit part of the convection/diffusion terms of a standard transport equation of a scalar field ![]() | |
cs_cavitation_parameters_t * | cs_get_glob_cavitation_parameters (void) |
void | cs_vof_solve_void_fraction (cs_real_t dt[], int iterns) |
Solve the void fraction ![]() | |
Variables | |
const cs_vof_parameters_t * | cs_glob_vof_parameters |
cs_cavitation_parameters_t* cs_get_glob_cavitation_parameters | ( | void | ) |
cs_vof_parameters_t* cs_get_glob_vof_parameters | ( | void | ) |
void cs_vof_compute_linear_rho_mu | ( | const cs_mesh_t * | m | ) |
Compute the mixture density, mixture dynamic viscosity given fluid volume fractions and the reference density and dynamic viscosity (liquid),
(gas).
Computation is done as follows on cells:
A similar linear formula is followed on boundary using fluid volume fraction value on the boundary.
[in] | m | pointer to mesh structure |
void cs_vof_deshpande_drift_flux | ( | const cs_mesh_t * | m, |
const cs_mesh_quantities_t * | mq | ||
) |
Compute the flux of the drift velocity , by using the flux of the standard velocity
following the approach described by Suraj S Deshpande et al 2012 Comput. Sci. Disc. 5 014016. It is activated with the option idrift = 1.
Using the notation:
The drift flux is computed as:
Where is the drift flux factor defined with the variable cdrift,
the normal vector to the interface. The gradient is computed using a centered scheme:
[in] | m | pointer to mesh structure |
[in] | mq | pointer to mesh quantities structure |
Compute the flux of the drift velocity , by using the flux of the standard velocity
following the approach described by Suraj S Deshpande et al 2012 Comput. Sci. Disc. 5 014016. It is activated with the option idrift = 1.
Using the notation:
The drift flux is computed as:
Where is the drift flux factor defined with the variable cdrift,
the normal vector to the interface. The gradient is computed using a centered scheme:
[in] | m | pointer to mesh structure |
[in] | mq | pointer to mesh quantities structure |
void cs_vof_drift_term | ( | int | imrgra, |
int | nswrgp, | ||
int | imligp, | ||
int | iwarnp, | ||
cs_real_t | epsrgp, | ||
cs_real_t | climgp, | ||
cs_real_t *restrict | pvar, | ||
const cs_real_t *restrict | pvara, | ||
cs_real_t *restrict | rhs | ||
) |
Add the explicit part of the convection/diffusion terms of a standard transport equation of a scalar field .
More precisely, the right hand side is updated as follows:
[in] | imrgra | indicator
|
[in] | nswrgp | number of reconstruction sweeps for the gradients |
[in] | imligp | clipping gradient method
|
[in] | iwarnp | verbosity |
[in] | epsrgp | relative precision for the gradient reconstruction |
[in] | climgp | clipping coefficient for the computation of the gradient |
[in] | pvar | solved variable (current time step) |
[in] | pvara | solved variable (previous time step) |
[in,out] | rhs | right hand side ![]() |
Add the explicit part of the convection/diffusion terms of a standard transport equation of a scalar field .
More precisely, the right hand side is updated as follows:
[in] | imrgra | indicator
|
[in] | nswrgp | number of reconstruction sweeps for the gradients |
[in] | imligp | clipping gradient method
|
[in] | iwarnp | verbosity |
[in] | epsrgp | relative precision for the gradient reconstruction |
[in] | climgp | clipping coefficient for the computation of the gradient |
[in] | pvar | solved variable (current time step) |
[in] | pvara | solved variable (previous time step) |
[in,out] | rhs | right hand side ![]() |
void cs_vof_log_mass_budget | ( | const cs_mesh_t * | m, |
const cs_mesh_quantities_t * | mq | ||
) |
Write in main log the global mixture mass budget:
[in] | m | pointer to mesh structure |
[in] | mq | pointer to mesh quantities structure |
[in] | m | pointer to mesh structure |
[in] | mq | pointer to mesh quantities structure |
void cs_vof_solve_void_fraction | ( | cs_real_t | dt[], |
int | iterns | ||
) |
Solve the void fraction for the Volume of Fluid method (and hence for cavitating flows).
This function solves:
with the eventual vaporization source term (Merkle model) in case the cavitation model is enabled,
the reference gas density and
the drift velocity for the compressed interface.
void cs_vof_surface_tension | ( | const cs_mesh_t * | m, |
const cs_mesh_quantities_t * | mq, | ||
cs_real_3_t | stf[] | ||
) |
Compute the surface tension momentum source term following the CSF model of Brackbill et al. (1992).
[in] | m | pointer to mesh structure |
[in] | mq | pointer to mesh quantities structure |
[out] | stf | surface tension momentum source term |
void cs_vof_update_phys_prop | ( | const cs_mesh_t * | m | ) |
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 (liquid),
(gas).
For the computation of mixture density, mixture dynamic viscosity, see cs_vof_compute_linear_rho_mu.
Computation of mass flux is as follows:
[in] | m | pointer to mesh structure |
|
extern |