#include "cs_base.h"#include "cs_field.h"#include "cs_halo.h"#include "cs_gradient.h"#include "cs_mesh.h"#include "cs_mesh_quantities.h"Go to the source code of this file.
Functions | |
| void | cs_mass_flux (const cs_mesh_t *m, const cs_mesh_quantities_t *fvq, int f_id, int itypfl, int iflmb0, int init, int inc, int imrgra, int nswrgu, cs_gradient_limit_t imligu, int iwarnu, double epsrgu, double climgu, const cs_real_t rom[], const cs_real_t romb[], const cs_real_3_t vel[], cs_field_bc_coeffs_t *bc_coeffs_v, cs_real_t *i_massflux, cs_real_t *b_massflux) |
| Add | |
| void | cs_divergence (const cs_mesh_t *m, int init, const cs_real_t i_massflux[], const cs_real_t b_massflux[], cs_real_t *diverg) |
| Add the integrated mass flux on the cells. | |
| void | cs_tensor_divergence (const cs_mesh_t *m, int init, const cs_real_3_t i_massflux[], const cs_real_3_t b_massflux[], cs_real_3_t *diverg) |
| Add the integrated mass flux on the cells for a tensor variable. | |
| void | cs_ext_force_flux (const cs_mesh_t *m, cs_mesh_quantities_t *fvq, int init, int nswrgu, const cs_real_3_t frcxt[], const cs_real_t cofbfp[], cs_real_t *i_massflux, cs_real_t *b_massflux, const cs_real_t i_visc[], const cs_real_t b_visc[], const cs_real_t viselx[], const cs_real_t visely[], const cs_real_t viselz[]) |
| Project the external source terms to the faces in coherence with cs_face_diffusion_scalar for the improved hydrostatic pressure algorithm (iphydr=1). | |
| void | cs_ext_force_anisotropic_flux (const cs_mesh_t *m, cs_mesh_quantities_t *fvq, int init, int nswrgp, int ircflp, const cs_real_3_t frcxt[], const cs_real_t cofbfp[], const cs_real_t i_visc[], const cs_real_t b_visc[], cs_real_6_t viscel[], const cs_real_2_t weighf[], cs_real_t *i_massflux, cs_real_t *b_massflux) |
| Project the external source terms to the faces in coherence with cs_face_anisotropic_diffusion_scalar for the improved hydrostatic pressure algorithm (iphydr=1). | |
| void | cs_tensor_face_flux (const cs_mesh_t *m, const cs_mesh_quantities_t *fvq, int f_id, int itypfl, int iflmb0, int init, int inc, int imrgra, int nswrgu, cs_gradient_limit_t imligu, int iwarnu, double epsrgu, double climgu, const cs_real_t c_rho[], const cs_real_t b_rho[], const cs_real_6_t c_var[], const cs_field_bc_coeffs_t *bc_coeffs_ts, cs_real_3_t *i_massflux, cs_real_3_t *b_massflux) |
| Add | |
| void cs_divergence | ( | const cs_mesh_t * | m, |
| int | init, | ||
| const cs_real_t | i_massflux[], | ||
| const cs_real_t | b_massflux[], | ||
| cs_real_t * | diverg ) |
Add the integrated mass flux on the cells.
![\[\dot{m}_i = \dot{m}_i + \sum_{\fij \in \Facei{\celli}} \dot{m}_\ij
\]](form_42.png)
| [in] | m | pointer to mesh |
| [in] | init | indicator
|
| [in] | i_massflux | mass flux at interior faces |
| [in] | b_massflux | mass flux at boundary faces |
| [in,out] | diverg | mass flux divergence |
| void cs_ext_force_anisotropic_flux | ( | const cs_mesh_t * | m, |
| cs_mesh_quantities_t * | fvq, | ||
| int | init, | ||
| int | nswrgp, | ||
| int | ircflp, | ||
| const cs_real_3_t | frcxt[], | ||
| const cs_real_t | cofbfp[], | ||
| const cs_real_t | i_visc[], | ||
| const cs_real_t | b_visc[], | ||
| cs_real_6_t | viscel[], | ||
| const cs_real_2_t | weighf[], | ||
| cs_real_t * | i_massflux, | ||
| cs_real_t * | b_massflux ) |
Project the external source terms to the faces in coherence with cs_face_anisotropic_diffusion_scalar for the improved hydrostatic pressure algorithm (iphydr=1).
| [in] | m | pointer to mesh |
| [in] | fvq | pointer to finite volume quantities |
| [in] | init | indicator
|
| [in] | nswrgp | number of reconstruction sweeps for the gradients |
| [in] | ircflp | indicator
|
| [in] | frcxt | body force creating the hydrostatic pressure |
| [in] | cofbfp | boundary condition array for the diffusion of the variable (implicit part) |
| [in] | i_visc | |
| [in] | b_visc | |
| [in] | viscel | symmetric cell tensor |
| [in] | weighf | internal face weight between cells i j in case of tensor diffusion |
| [in,out] | i_massflux | mass flux at interior faces |
| [in,out] | b_massflux | mass flux at boundary faces |
| void cs_ext_force_flux | ( | const cs_mesh_t * | m, |
| cs_mesh_quantities_t * | fvq, | ||
| int | init, | ||
| int | nswrgu, | ||
| const cs_real_3_t | frcxt[], | ||
| const cs_real_t | cofbfp[], | ||
| cs_real_t * | i_massflux, | ||
| cs_real_t * | b_massflux, | ||
| const cs_real_t | i_visc[], | ||
| const cs_real_t | b_visc[], | ||
| const cs_real_t | viselx[], | ||
| const cs_real_t | visely[], | ||
| const cs_real_t | viselz[] ) |
Project the external source terms to the faces in coherence with cs_face_diffusion_scalar for the improved hydrostatic pressure algorithm (iphydr=1).
| [in] | m | pointer to mesh |
| [in] | fvq | pointer to finite volume quantities |
| [in] | init | indicator
|
| [in] | nswrgu | number of reconstruction sweeps for the gradients |
| [in] | frcxt | body force creating the hydrostatic pressure |
| [in] | cofbfp | boundary condition array for the diffusion of the variable (implicit part) |
| [in,out] | i_massflux | mass flux at interior faces |
| [in,out] | b_massflux | mass flux at boundary faces |
| [in] | i_visc | |
| [in] | b_visc | |
| [in] | viselx | viscosity by cell, dir x |
| [in] | visely | viscosity by cell, dir y |
| [in] | viselz | viscosity by cell, dir z |
| void cs_mass_flux | ( | const cs_mesh_t * | m, |
| const cs_mesh_quantities_t * | fvq, | ||
| int | f_id, | ||
| int | itypfl, | ||
| int | iflmb0, | ||
| int | init, | ||
| int | inc, | ||
| int | imrgra, | ||
| int | nswrgu, | ||
| cs_gradient_limit_t | imligu, | ||
| int | iwarnu, | ||
| double | epsrgu, | ||
| double | climgu, | ||
| const cs_real_t | rom[], | ||
| const cs_real_t | romb[], | ||
| const cs_real_3_t | vel[], | ||
| cs_field_bc_coeffs_t * | bc_coeffs_v, | ||
| cs_real_t * | i_massflux, | ||
| cs_real_t * | b_massflux ) |
Add 

For the reconstruction, 


For the mass flux at the boundary we have:
![\[\dot{m}_\ib = \left[ \rho_\fib \vect{a}_u + \rho_\fib \tens{b}_u \vect{u}
+ \tens{b}_u \left(\gradt \vect{u} \cdot \vect{\centi \centip}\right)\right]
\cdot \vect{s}_\ij
\]](form_37.png)
The last equation uses some approximations detailed in the theory guide.
| [in] | m | pointer to mesh |
| [in] | fvq | pointer to finite volume quantities |
| [in] | f_id | field id (or -1) |
| [in] | itypfl | indicator (take rho into account or not)
|
| [in] | iflmb0 | the mass flux is set to 0 on walls and symmetries if = 1 |
| [in] | init | the mass flux is initialized to 0 if > 0 |
| [in] | inc | indicator
|
| [in] | imrgra | indicator
|
| [in] | nswrgu | number of sweeps for the reconstruction of the gradients |
| [in] | imligu | clipping gradient method
|
| [in] | iwarnu | verbosity |
| [in] | epsrgu | relative precision for the gradient reconstruction |
| [in] | climgu | clipping coefficient for the computation of the gradient |
| [in] | rom | cell density |
| [in] | romb | density at boundary faces |
| [in] | vel | vector variable |
| [in] | bc_coeff_v | BC structure for the vector variable |
| [in,out] | i_massflux | mass flux at interior faces |
| [in,out] | b_massflux | mass flux at boundary faces |
| void cs_tensor_divergence | ( | const cs_mesh_t * | m, |
| int | init, | ||
| const cs_real_3_t | i_massflux[], | ||
| const cs_real_3_t | b_massflux[], | ||
| cs_real_3_t * | diverg ) |
Add the integrated mass flux on the cells for a tensor variable.
![\[\dot{m}_i = \dot{m}_i + \sum_{\fij \in \Facei{\celli}} \dot{m}_\ij
\]](form_42.png)
| [in] | m | pointer to mesh |
| [in] | init | indicator
|
| [in] | i_massflux | mass flux vector at interior faces |
| [in] | b_massflux | mass flux vector at boundary faces |
| [in,out] | diverg | mass flux divergence vector |
| void cs_tensor_face_flux | ( | const cs_mesh_t * | m, |
| const cs_mesh_quantities_t * | fvq, | ||
| int | f_id, | ||
| int | itypfl, | ||
| int | iflmb0, | ||
| int | init, | ||
| int | inc, | ||
| int | imrgra, | ||
| int | nswrgu, | ||
| cs_gradient_limit_t | imligu, | ||
| int | iwarnu, | ||
| double | epsrgu, | ||
| double | climgu, | ||
| const cs_real_t | c_rho[], | ||
| const cs_real_t | b_rho[], | ||
| const cs_real_6_t | c_var[], | ||
| const cs_field_bc_coeffs_t * | bc_coeffs_ts, | ||
| cs_real_3_t * | i_massflux, | ||
| cs_real_3_t * | b_massflux ) |
Add 
| [in] | m | pointer to mesh |
| [in] | fvq | pointer to finite volume quantities |
| [in] | f_id | field id (or -1) |
| [in] | itypfl | indicator (take rho into account or not)
|
| [in] | iflmb0 | the mass flux is set to 0 on walls and symmetries if = 1 |
| [in] | init | the mass flux is initialized to 0 if > 0 |
| [in] | inc | indicator
|
| [in] | imrgra | indicator
|
| [in] | nswrgu | number of sweeps for the reconstruction of the gradients |
| [in] | imligu | clipping gradient method
|
| [in] | iwarnu | verbosity |
| [in] | epsrgu | relative precision for the gradient reconstruction |
| [in] | climgu | clipping coefficient for the computation of the gradient |
| [in] | c_rho | cell density |
| [in] | b_rho | density at boundary faces |
| [in] | c_var | variable |
| [in] | bc_coeffs_ts | boundary condition structure for the variable |
| [in,out] | i_massflux | mass flux at interior faces |
| [in,out] | b_massflux | mass flux at boundary faces |