navier_stokes_standard

Navier-Stokes equations.

Parameters:

[solveur_pression] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) Linear pressure system resolution method.
[dt_projection] (type: deuxmots) nb value : This keyword checks every nb time-steps the equality of velocity divergence to zero. value is the criteria convergency for the solver used.
[traitement_particulier] (type: traitement_particulier) Keyword to post-process particular values.
[seuil_divu] (type: floatfloat) this keyword is intended to minimise the number of iterations during the pressure system resolution. The convergence criteria during this step (seuil in solveur_pression) is dynamically adapted according to the mass conservation. At tn , the linear system Ax=B is considered as solved if the residual ||Ax-B||<seuil(tn). For tn+1, the threshold value seuil(tn+1) will be evaluated as:

If ( |max(DivU)*dt| < value )

Seuil(tn+1)= Seuil(tn)*factor

Else

Seuil(tn+1)= Seuil(tn)*factor

Endif

The first parameter (value) is the mass evolution the user is ready to accept per timestep, and the second one (factor) is the factor of evolution for seuil (for example 1.1, so 10% per timestep). Investigations has to be lead to know more about the effects of these two last parameters on the behaviour of the simulations.
[solveur_bar] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) This keyword is used to define when filtering operation is called (typically for EF convective scheme, standard diffusion operator and Source_Qdm_lambdaup ). A file (solveur.bar) is then created and used for inversion procedure. Syntax is the same then for pressure solver (GCP is required for multi-processor calculations and, in a general way, for big meshes).
[projection_initiale] (type: int) Keyword to suppress, if boolean equals 0, the initial projection which checks DivU=0. By default, boolean equals 1.
[postraiter_gradient_pression_sans_masse] (type: flag) Avoid mass matrix multiplication for the gradient postprocessing
[methode_calcul_pression_initiale] (type: string into ['avec_les_cl', 'avec_sources', 'avec_sources_et_operateurs', 'sans_rien']) Keyword to select an option for the pressure calculation before the fist time step. Options are : avec_les_cl (default option lapP=0 is solved with Neuman boundary conditions on pressure if any), avec_sources (lapP=f is solved with Neuman boundaries conditions and f integrating the source terms of the Navier-Stokes equations) and avec_sources_et_operateurs (lapP=f is solved as with the previous option avec_sources but f integrating also some operators of the Navier-Stokes equations). The two last options are useful and sometime necessary when source terms are implicited when using an implicit time scheme to solve the Navier-Stokes equations.
Inherited from eqn_base:
[disable_equation_residual] (type: int) The equation residual will not be used for the problem residual used when checking time convergence or computing dynamic time-step
[convection] (type: bloc_convection) Keyword to alter the convection scheme.
[diffusion] (type: bloc_diffusion) Keyword to specify the diffusion operator.
[conditions_limites | boundary_conditions] (type: list of Condlimlu) Boundary conditions.
[conditions_initiales | initial_conditions] (type: list of Condinit) Initial conditions.
[sources] (type: list of Source_base) The sources.
[ecrire_fichier_xyz_valeur] (type: ecrire_fichier_xyz_valeur) This keyword is used to write the values of a field only for some boundaries in a text file
[parametre_equation] (type: parametre_equation_base — one of: parametre_diffusion_implicite, parametre_implicite) Keyword used to specify additional parameters for the equation
[equation_non_resolue] (type: string) The equation will not be solved while condition(t) is verified if equation_non_resolue keyword is used. Exemple: The Navier-Stokes equations are not solved between time t0 and t1.

Navier_Sokes_Standard

{ equation_non_resolue (t>t0)*(t<t1) }
[renommer_equation | rename_equation] (type: string) Rename the equation with a specific name.

Keywords derived from navier_stokes_turbulent

navier_stokes_ibm

Inherits from: navier_stokes_standard
Defined in: src/ThHyd/IBM/Equations/Navier_Stokes_IBM.cpp
Used by: pb_hydraulique_ibm, pb_thermohydraulique_ibm

IBM Navier-Stokes equations.

Parameters:

[correction_matrice_projection_initiale] (type: int) (IBM advanced) fix matrix of initial projection for PDF
[correction_calcul_pression_initiale] (type: int) (IBM advanced) fix initial pressure computation for PDF
[correction_vitesse_projection_initiale] (type: int) (IBM advanced) fix initial velocity computation for PDF
[correction_matrice_pression] (type: int) (IBM advanced) fix pressure matrix for PDF
[matrice_pression_penalisee_h1] (type: int) (IBM advanced) fix pressure matrix for PDF
[correction_vitesse_modifie] (type: int) (IBM advanced) fix velocity for PDF
[correction_pression_modifie] (type: int) (IBM advanced) fix pressure for PDF
[gradient_pression_qdm_modifie] (type: int) (IBM advanced) fix pressure gradient
[correction_variable_initiale] (type: int) Modify initial variable
Inherited from navier_stokes_standard:
[solveur_pression] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) Linear pressure system resolution method.
[dt_projection] (type: deuxmots) nb value : This keyword checks every nb time-steps the equality of velocity divergence to zero. value is the criteria convergency for the solver used.
[traitement_particulier] (type: traitement_particulier) Keyword to post-process particular values.
[seuil_divu] (type: floatfloat) this keyword is intended to minimise the number of iterations during the pressure system resolution. The convergence criteria during this step (seuil in solveur_pression) is dynamically adapted according to the mass conservation. At tn , the linear system Ax=B is considered as solved if the residual ||Ax-B||<seuil(tn). For tn+1, the threshold value seuil(tn+1) will be evaluated as:

If ( |max(DivU)*dt| < value )

Seuil(tn+1)= Seuil(tn)*factor

Else

Seuil(tn+1)= Seuil(tn)*factor

Endif

The first parameter (value) is the mass evolution the user is ready to accept per timestep, and the second one (factor) is the factor of evolution for seuil (for example 1.1, so 10% per timestep). Investigations has to be lead to know more about the effects of these two last parameters on the behaviour of the simulations.
[solveur_bar] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) This keyword is used to define when filtering operation is called (typically for EF convective scheme, standard diffusion operator and Source_Qdm_lambdaup ). A file (solveur.bar) is then created and used for inversion procedure. Syntax is the same then for pressure solver (GCP is required for multi-processor calculations and, in a general way, for big meshes).
[projection_initiale] (type: int) Keyword to suppress, if boolean equals 0, the initial projection which checks DivU=0. By default, boolean equals 1.
[postraiter_gradient_pression_sans_masse] (type: flag) Avoid mass matrix multiplication for the gradient postprocessing
[methode_calcul_pression_initiale] (type: string into ['avec_les_cl', 'avec_sources', 'avec_sources_et_operateurs', 'sans_rien']) Keyword to select an option for the pressure calculation before the fist time step. Options are : avec_les_cl (default option lapP=0 is solved with Neuman boundary conditions on pressure if any), avec_sources (lapP=f is solved with Neuman boundaries conditions and f integrating the source terms of the Navier-Stokes equations) and avec_sources_et_operateurs (lapP=f is solved as with the previous option avec_sources but f integrating also some operators of the Navier-Stokes equations). The two last options are useful and sometime necessary when source terms are implicited when using an implicit time scheme to solve the Navier-Stokes equations.
Inherited from eqn_base:
[disable_equation_residual] (type: int) The equation residual will not be used for the problem residual used when checking time convergence or computing dynamic time-step
[convection] (type: bloc_convection) Keyword to alter the convection scheme.
[diffusion] (type: bloc_diffusion) Keyword to specify the diffusion operator.
[conditions_limites | boundary_conditions] (type: list of Condlimlu) Boundary conditions.
[conditions_initiales | initial_conditions] (type: list of Condinit) Initial conditions.
[sources] (type: list of Source_base) The sources.
[ecrire_fichier_xyz_valeur] (type: ecrire_fichier_xyz_valeur) This keyword is used to write the values of a field only for some boundaries in a text file
[parametre_equation] (type: parametre_equation_base — one of: parametre_diffusion_implicite, parametre_implicite) Keyword used to specify additional parameters for the equation
[equation_non_resolue] (type: string) The equation will not be solved while condition(t) is verified if equation_non_resolue keyword is used. Exemple: The Navier-Stokes equations are not solved between time t0 and t1.

Navier_Sokes_Standard

{ equation_non_resolue (t>t0)*(t<t1) }
[renommer_equation | rename_equation] (type: string) Rename the equation with a specific name.

navier_stokes_ibm_turbulent

Inherits from: navier_stokes_standard
Defined in: src/ThHyd/IBM/Equations/Navier_Stokes_IBM_Turbulent.cpp
Used by: pb_hydraulique_ibm_turbulent, pb_thermohydraulique_ibm_turbulent

IBM Navier-Stokes equations as well as the associated turbulence model equations.

Parameters:

[modele_turbulence] (type: modele_turbulence_hyd_deriv — one of: combinaison, k_epsilon_realisable, k_epsilon … (+21 more)) Turbulence model for Navier-Stokes equations.
Inherited from navier_stokes_standard:
[solveur_pression] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) Linear pressure system resolution method.
[dt_projection] (type: deuxmots) nb value : This keyword checks every nb time-steps the equality of velocity divergence to zero. value is the criteria convergency for the solver used.
[traitement_particulier] (type: traitement_particulier) Keyword to post-process particular values.
[seuil_divu] (type: floatfloat) this keyword is intended to minimise the number of iterations during the pressure system resolution. The convergence criteria during this step (seuil in solveur_pression) is dynamically adapted according to the mass conservation. At tn , the linear system Ax=B is considered as solved if the residual ||Ax-B||<seuil(tn). For tn+1, the threshold value seuil(tn+1) will be evaluated as:

If ( |max(DivU)*dt| < value )

Seuil(tn+1)= Seuil(tn)*factor

Else

Seuil(tn+1)= Seuil(tn)*factor

Endif

The first parameter (value) is the mass evolution the user is ready to accept per timestep, and the second one (factor) is the factor of evolution for seuil (for example 1.1, so 10% per timestep). Investigations has to be lead to know more about the effects of these two last parameters on the behaviour of the simulations.
[solveur_bar] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) This keyword is used to define when filtering operation is called (typically for EF convective scheme, standard diffusion operator and Source_Qdm_lambdaup ). A file (solveur.bar) is then created and used for inversion procedure. Syntax is the same then for pressure solver (GCP is required for multi-processor calculations and, in a general way, for big meshes).
[projection_initiale] (type: int) Keyword to suppress, if boolean equals 0, the initial projection which checks DivU=0. By default, boolean equals 1.
[postraiter_gradient_pression_sans_masse] (type: flag) Avoid mass matrix multiplication for the gradient postprocessing
[methode_calcul_pression_initiale] (type: string into ['avec_les_cl', 'avec_sources', 'avec_sources_et_operateurs', 'sans_rien']) Keyword to select an option for the pressure calculation before the fist time step. Options are : avec_les_cl (default option lapP=0 is solved with Neuman boundary conditions on pressure if any), avec_sources (lapP=f is solved with Neuman boundaries conditions and f integrating the source terms of the Navier-Stokes equations) and avec_sources_et_operateurs (lapP=f is solved as with the previous option avec_sources but f integrating also some operators of the Navier-Stokes equations). The two last options are useful and sometime necessary when source terms are implicited when using an implicit time scheme to solve the Navier-Stokes equations.
Inherited from eqn_base:
[disable_equation_residual] (type: int) The equation residual will not be used for the problem residual used when checking time convergence or computing dynamic time-step
[convection] (type: bloc_convection) Keyword to alter the convection scheme.
[diffusion] (type: bloc_diffusion) Keyword to specify the diffusion operator.
[conditions_limites | boundary_conditions] (type: list of Condlimlu) Boundary conditions.
[conditions_initiales | initial_conditions] (type: list of Condinit) Initial conditions.
[sources] (type: list of Source_base) The sources.
[ecrire_fichier_xyz_valeur] (type: ecrire_fichier_xyz_valeur) This keyword is used to write the values of a field only for some boundaries in a text file
[parametre_equation] (type: parametre_equation_base — one of: parametre_diffusion_implicite, parametre_implicite) Keyword used to specify additional parameters for the equation
[equation_non_resolue] (type: string) The equation will not be solved while condition(t) is verified if equation_non_resolue keyword is used. Exemple: The Navier-Stokes equations are not solved between time t0 and t1.

Navier_Sokes_Standard

{ equation_non_resolue (t>t0)*(t<t1) }
[renommer_equation | rename_equation] (type: string) Rename the equation with a specific name.

navier_stokes_qc

Inherits from: navier_stokes_standard
Defined in: src/ThHyd/Dilatable/Quasi_Compressible/Equations/Navier_Stokes_QC.cpp
Used by: pb_hydraulique_melange_binaire_qc, pb_thermohydraulique_especes_qc, pb_thermohydraulique_qc

Navier-Stokes equation for a quasi-compressible fluid.

Parameters:

Inherited from navier_stokes_standard:
[solveur_pression] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) Linear pressure system resolution method.
[dt_projection] (type: deuxmots) nb value : This keyword checks every nb time-steps the equality of velocity divergence to zero. value is the criteria convergency for the solver used.
[traitement_particulier] (type: traitement_particulier) Keyword to post-process particular values.
[seuil_divu] (type: floatfloat) this keyword is intended to minimise the number of iterations during the pressure system resolution. The convergence criteria during this step (seuil in solveur_pression) is dynamically adapted according to the mass conservation. At tn , the linear system Ax=B is considered as solved if the residual ||Ax-B||<seuil(tn). For tn+1, the threshold value seuil(tn+1) will be evaluated as:

If ( |max(DivU)*dt| < value )

Seuil(tn+1)= Seuil(tn)*factor

Else

Seuil(tn+1)= Seuil(tn)*factor

Endif

The first parameter (value) is the mass evolution the user is ready to accept per timestep, and the second one (factor) is the factor of evolution for seuil (for example 1.1, so 10% per timestep). Investigations has to be lead to know more about the effects of these two last parameters on the behaviour of the simulations.
[solveur_bar] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) This keyword is used to define when filtering operation is called (typically for EF convective scheme, standard diffusion operator and Source_Qdm_lambdaup ). A file (solveur.bar) is then created and used for inversion procedure. Syntax is the same then for pressure solver (GCP is required for multi-processor calculations and, in a general way, for big meshes).
[projection_initiale] (type: int) Keyword to suppress, if boolean equals 0, the initial projection which checks DivU=0. By default, boolean equals 1.
[postraiter_gradient_pression_sans_masse] (type: flag) Avoid mass matrix multiplication for the gradient postprocessing
[methode_calcul_pression_initiale] (type: string into ['avec_les_cl', 'avec_sources', 'avec_sources_et_operateurs', 'sans_rien']) Keyword to select an option for the pressure calculation before the fist time step. Options are : avec_les_cl (default option lapP=0 is solved with Neuman boundary conditions on pressure if any), avec_sources (lapP=f is solved with Neuman boundaries conditions and f integrating the source terms of the Navier-Stokes equations) and avec_sources_et_operateurs (lapP=f is solved as with the previous option avec_sources but f integrating also some operators of the Navier-Stokes equations). The two last options are useful and sometime necessary when source terms are implicited when using an implicit time scheme to solve the Navier-Stokes equations.
Inherited from eqn_base:
[disable_equation_residual] (type: int) The equation residual will not be used for the problem residual used when checking time convergence or computing dynamic time-step
[convection] (type: bloc_convection) Keyword to alter the convection scheme.
[diffusion] (type: bloc_diffusion) Keyword to specify the diffusion operator.
[conditions_limites | boundary_conditions] (type: list of Condlimlu) Boundary conditions.
[conditions_initiales | initial_conditions] (type: list of Condinit) Initial conditions.
[sources] (type: list of Source_base) The sources.
[ecrire_fichier_xyz_valeur] (type: ecrire_fichier_xyz_valeur) This keyword is used to write the values of a field only for some boundaries in a text file
[parametre_equation] (type: parametre_equation_base — one of: parametre_diffusion_implicite, parametre_implicite) Keyword used to specify additional parameters for the equation
[equation_non_resolue] (type: string) The equation will not be solved while condition(t) is verified if equation_non_resolue keyword is used. Exemple: The Navier-Stokes equations are not solved between time t0 and t1.

Navier_Sokes_Standard

{ equation_non_resolue (t>t0)*(t<t1) }
[renommer_equation | rename_equation] (type: string) Rename the equation with a specific name.

navier_stokes_wc

Inherits from: navier_stokes_standard
Defined in: src/ThHyd/Dilatable/Weakly_Compressible/Equations/Navier_Stokes_WC.cpp
Used by: pb_hydraulique_melange_binaire_wc, pb_thermohydraulique_especes_wc, pb_thermohydraulique_wc

Navier-Stokes equation for a weakly-compressible fluid.

Parameters:

[mass_source] (type: mass_source) Mass source used in a dilatable simulation to add/reduce a mass at the boundary (volumetric source in the first cell of a given boundary).
Inherited from navier_stokes_standard:
[solveur_pression] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) Linear pressure system resolution method.
[dt_projection] (type: deuxmots) nb value : This keyword checks every nb time-steps the equality of velocity divergence to zero. value is the criteria convergency for the solver used.
[traitement_particulier] (type: traitement_particulier) Keyword to post-process particular values.
[seuil_divu] (type: floatfloat) this keyword is intended to minimise the number of iterations during the pressure system resolution. The convergence criteria during this step (seuil in solveur_pression) is dynamically adapted according to the mass conservation. At tn , the linear system Ax=B is considered as solved if the residual ||Ax-B||<seuil(tn). For tn+1, the threshold value seuil(tn+1) will be evaluated as:

If ( |max(DivU)*dt| < value )

Seuil(tn+1)= Seuil(tn)*factor

Else

Seuil(tn+1)= Seuil(tn)*factor

Endif

The first parameter (value) is the mass evolution the user is ready to accept per timestep, and the second one (factor) is the factor of evolution for seuil (for example 1.1, so 10% per timestep). Investigations has to be lead to know more about the effects of these two last parameters on the behaviour of the simulations.
[solveur_bar] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) This keyword is used to define when filtering operation is called (typically for EF convective scheme, standard diffusion operator and Source_Qdm_lambdaup ). A file (solveur.bar) is then created and used for inversion procedure. Syntax is the same then for pressure solver (GCP is required for multi-processor calculations and, in a general way, for big meshes).
[projection_initiale] (type: int) Keyword to suppress, if boolean equals 0, the initial projection which checks DivU=0. By default, boolean equals 1.
[postraiter_gradient_pression_sans_masse] (type: flag) Avoid mass matrix multiplication for the gradient postprocessing
[methode_calcul_pression_initiale] (type: string into ['avec_les_cl', 'avec_sources', 'avec_sources_et_operateurs', 'sans_rien']) Keyword to select an option for the pressure calculation before the fist time step. Options are : avec_les_cl (default option lapP=0 is solved with Neuman boundary conditions on pressure if any), avec_sources (lapP=f is solved with Neuman boundaries conditions and f integrating the source terms of the Navier-Stokes equations) and avec_sources_et_operateurs (lapP=f is solved as with the previous option avec_sources but f integrating also some operators of the Navier-Stokes equations). The two last options are useful and sometime necessary when source terms are implicited when using an implicit time scheme to solve the Navier-Stokes equations.
Inherited from eqn_base:
[disable_equation_residual] (type: int) The equation residual will not be used for the problem residual used when checking time convergence or computing dynamic time-step
[convection] (type: bloc_convection) Keyword to alter the convection scheme.
[diffusion] (type: bloc_diffusion) Keyword to specify the diffusion operator.
[conditions_limites | boundary_conditions] (type: list of Condlimlu) Boundary conditions.
[conditions_initiales | initial_conditions] (type: list of Condinit) Initial conditions.
[sources] (type: list of Source_base) The sources.
[ecrire_fichier_xyz_valeur] (type: ecrire_fichier_xyz_valeur) This keyword is used to write the values of a field only for some boundaries in a text file
[parametre_equation] (type: parametre_equation_base — one of: parametre_diffusion_implicite, parametre_implicite) Keyword used to specify additional parameters for the equation
[equation_non_resolue] (type: string) The equation will not be solved while condition(t) is verified if equation_non_resolue keyword is used. Exemple: The Navier-Stokes equations are not solved between time t0 and t1.

Navier_Sokes_Standard

{ equation_non_resolue (t>t0)*(t<t1) }
[renommer_equation | rename_equation] (type: string) Rename the equation with a specific name.