TrioCFD 1.9.8
TrioCFD documentation
Loading...
Searching...
No Matches
Keywords derived from mor_eqn

Keywords introduced by TrioCFD that extend the Mor_eqn family. See Keywords derived from mor_eqn for the base-class documentation and the canonical keyword list.


convection_diffusion_concentration_ft_disc

not_set

Parameters:

  • [equation_interface] (type: string) his is the name of the interface tracking equation to watch. The scalar will not diffuse through the interface of this equation.
  • phase (type: int into [0, 1]) tells whether the scalar must be confined in phase 0 or in phase 1
  • [option] (type: string) Experimental features used to prevent the concentration to leak through the interface between phases due to numerical diffusion. RIEN: do nothing RAMASSE_MIETTES_SIMPLE: at each timestep, this algorithm takes all the mass located in the opposite phase and spreads it uniformly in the given phase.
  • [nom_inconnue] (type: string) Keyword Nom_inconnue will rename the unknown of this equation with the given name. In the postprocessing part, the concentration field will be accessible with this name. This is usefull if you want to track more than one concentration (otherwise, only the concentration field in the first concentration equation can be accessed).
  • [alias] (type: string) not_set
  • [masse_molaire] (type: float) not_set
  • [is_multi_scalar | is_multi_scalar_diffusion] (type: flag) Flag to activate the multi_scalar diffusion operator
  • [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) 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.

convection_diffusion_concentration_turbulent_ft_disc

equation_non_resolue

Parameters:

  • [equation_interface] (type: string) his is the name of the interface tracking equation to watch. The scalar will not diffuse through the interface of this equation.
  • phase (type: int into [0, 1]) tells whether the scalar must be confined in phase 0 or in phase 1
  • [option] (type: string) Experimental features used to prevent the concentration to leak through the interface between phases due to numerical diffusion. RIEN: do nothing RAMASSE_MIETTES_SIMPLE: at each timestep, this algorithm takes all the mass located in the opposite phase and spreads it uniformly in the given phase.
  • [equations_source_chimie] (type: list of str) This term specifies the name of the concentration equation of the reagents. It should be specified only in the bloc that concerns the convection/diffusion equation of the product.
  • [modele_cinetique] (type: int) This is the keyword that the user defines for the reaction model that he wants to use. Four reaction models are currently offered (1 to 4). Model 1 is the default one and is based on the laminar rate formulation. Model 2 employs an LES diffusive EDC formulation. Model 3 defines an LES variance formulation. Model 4 is a mix between models 2 and 3.
  • [equation_nu_t] (type: string) This specifies the name of the hydraulic equation used which defines the turbulent (basically SGS) viscosity.
  • [constante_cinetique] (type: float) This is the constant kinetic rate of the reaction and is used for the laminar model 1 only.
  • [modele_turbulence] (type: modele_turbulence_scal_base) Turbulence model to be used in the constituent transport equations. The only model currently available is Schmidt.
  • [nom_inconnue] (type: string) Keyword Nom_inconnue will rename the unknown of this equation with the given name. In the postprocessing part, the concentration field will be accessible with this name. This is usefull if you want to track more than one concentration (otherwise, only the concentration field in the first concentration equation can be accessed).
  • [alias] (type: string) not_set
  • [masse_molaire] (type: float) not_set
  • [is_multi_scalar | is_multi_scalar_diffusion] (type: flag) Flag to activate the multi_scalar diffusion operator
  • [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) 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.

convection_diffusion_espece_multi_qc

Species conservation equation for a multi-species quasi-compressible fluid.

Parameters:

  • [espece] (type: espece) Assosciate a species (with its properties) to the equation
  • [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) 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.

convection_diffusion_phase_field

Cahn-Hilliard equation of the Phase Field problem. The unknown of this equation is the concentration C.

Parameters:

  • [mu_1] (type: float) Dynamic viscosity of the first phase.
  • [mu_2] (type: float) Dynamic viscosity of the second phase.
  • [rho_1] (type: float) Density of the first phase.
  • [rho_2] (type: float) Density of the second phase.
  • potentiel_chimique_generalise (type: string into ['avec_energie_cinetique', 'sans_energie_cinetique']) To define (chaine set to avec_energie_cinetique) or not (chaine set to sans_energie_cinetique) if the Cahn-Hilliard equation contains the cinetic energy term.
  • [nom_inconnue] (type: string) Keyword Nom_inconnue will rename the unknown of this equation with the given name. In the postprocessing part, the concentration field will be accessible with this name. This is usefull if you want to track more than one concentration (otherwise, only the concentration field in the first concentration equation can be accessed).
  • [alias] (type: string) not_set
  • [masse_molaire] (type: float) not_set
  • [is_multi_scalar | is_multi_scalar_diffusion] (type: flag) Flag to activate the multi_scalar diffusion operator
  • [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) 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.

convection_diffusion_temperature_ft_disc

not_set

Parameters:

  • [equation_interface] (type: string) The name of the interface equation should be given.
  • phase (type: int into [0, 1]) Phase in which the temperature equation will be solved. The temperature, which may be postprocessed with the keyword temperature_EquationName, in the orther phase may be negative: the code only computes the temperature field in the specified phase. The other phase is supposed to physically stay at saturation temperature. The code uses a ghost fluid numerical method to work on a smooth temperature field at the interface. In the opposite phase (1-X) the temperature will therefore be extrapolated in the vicinity of the interface and have the opposite sign, saturation temperature is zero by convention).
  • [equation_navier_stokes] (type: string) The name of the Navier Stokes equation of the problem should be given.
  • [stencil_width] (type: int) distance in mesh elements over which the temperature field should be extrapolated in the opposite phase.
  • [maintien_temperature] (type: objet_lecture_maintien_temperature) maintien_temperature SOUS_ZONE_NAME VALUE : experimental, this acts as a dynamic source term that heats or cools the fluid to maintain the average temperature to VALUE within the specified region. At this time, this is done by multiplying the temperature within the SOUS_ZONE by an appropriate uniform value at each timestep. This feature might be implemented in a separate source term in the future.
  • [prescribed_mpoint] (type: float) User defined value of the phase-change rate (override the value computed based on the temperature field)
  • [correction_mpoint_diff_conv_energy] (type: list of float) not_set
  • [correction_gradt_ordre] (type: int) not_set
  • [divergence_free_velocity_extension] (type: flag) not_set
  • [solveur_pression_fictive] (type: solveur_sys_base) not_set
  • [bc_opening_pressure] (type: list of str) not_set
  • [penalisation_l2_ftd] (type: bloc_lecture) to activate or not (the default is Direct Forcing method) the Penalized Direct Forcing method to impose the specified temperature on the solid-fluid interface.
  • [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) 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.

convection_diffusion_temperature_sensibility

Energy sensitivity equation (temperature diffusion convection)

Parameters:

  • [convection_sensibility | sensibility] (type: convection_deriv) Choice between: amont and muscl Example: convection { Sensibility { amont } }
  • velocity_state (type: bloc_lecture) Block to indicate the state problem. Between the braces, you must specify the key word 'pb_champ_evaluateur' then the name of the state problem and the velocity unknown Example: velocity_state { pb_champ_evaluateur pb_state velocity }
  • temperature_state (type: bloc_lecture) Block to indicate the state problem. Between the braces, you must specify the key word 'pb_champ_evaluateur' then the name of the state problem and the temperature unknown Example: velocity_state { pb_champ_evaluateur pb_state temperature }
  • uncertain_variable (type: bloc_lecture) Block to indicate the name of the uncertain variable. Between the braces, you must specify the name of the unknown variable (choice between: temperature, beta_th, boussinesq_temperature, Cp and lambda . Example: uncertain_variable { temperature }
  • [polynomial_chaos] (type: float) It is the method that we will use to study the sensitivity of the
  • [penalisation_l2_ftd] (type: bloc_lecture) to activate or not (the default is Direct Forcing method) the Penalized Direct Forcing method to impose the specified temperature on the solid-fluid interface.
  • [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) 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.

ijk_interfaces

Interface tracking equation for Front-Tracking problem in the discontinuous version.

Parameters:

  • [fichier_reprise_interface] (type: string) not_set
  • [timestep_reprise_interface] (type: float) not_set
  • [lata_meshname] (type: string) not_set
  • [remaillage_ft_ijk] (type: bloc_lecture) not_set
  • [parcours_interface] (type: bloc_lecture) not_set
  • [maillage_ft_ijk] (type: bloc_lecture) not_set
  • [terme_gravite] (type: string) not_set
  • [use_tryggvason_interfacial_source] (type: flag) not_set
  • [compute_distance_autres_interfaces] (type: flag) not_set
  • [avoid_duplicata] (type: flag) not_set
  • [portee_force_repulsion] (type: float) not_set
  • [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) 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_aposteriori

Modification of the Navier_Stokes_standard class in order to accept the estimateur_aposteriori post-processing. To post-process estimateur_aposteriori, add this keyword into the list of fields to be post-processed. This estimator whill generate a map of aposteriori error estimators; it is defined on each mesh cell and is a measure of the local discretisation error. This will serve for adaptive mesh refinement

Parameters:

  • [solveur_pression] (type: solveur_sys_base) 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) value factor : 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 evualated 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) 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.
  • **[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) 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_ft_disc

Two-phase momentum balance equation.

Parameters:

  • [equation_interfaces_proprietes_fluide] (type: string) This keyword is used for liquid-gas, liquid-vapor and fluid-fluid deformable interface, which transported at the Eulerian velocity. When this case is selected, the keyword sequence Methode_transport vitesse_interpolee is used in the block Transport_Interfaces_FT_Disc to define the velocity field for the displacement of the interface.
  • [equation_interfaces_vitesse_imposee] (type: string) This keyword is used to specify the velocity field to be used when using an interface that mimics a solid interface moving with a given solid speed of displacement. When this case is selected, the keyword sequence Methode_transport vitesse_imposee in the Transport_Interfaces_FT_Disc block will define the velocity field for the displacement of the interface.
  • [equations_interfaces_vitesse_imposee] (type: list of str) This keyword is used to specify the velocity field to be used when using an interface that mimics a solid interface moving with a given solid speed of displacement. When this case is selected, the keyword sequence Methode_transport vitesse_imposee in the Transport_Interfaces_FT_Disc block will define the velocity field for the displacement of the interface. If two or more solid interfaces are defined, then the keyword equations_interfaces_vitesse_imposee should be used.
  • [clipping_courbure_interface] (type: int) This keyword is used to numerically limit the values of curvature used in the momentum balance equation. Curvature is computed as usual, but values exceeding the clipping value are replaced by this threshold, before using the clipped curvature in the momentum balance. Each time a curvature value is clipped, a counter is increased by one unity and the value of the counter is written in the .err file at the end of the time step. This clipping allows not reducing drastically the time stepping when a geometrical singularity occurs in the interface mesh. However, physical phenomena may be concealed with the use of such a clipping.
  • [terme_gravite] (type: string into ['rho_g', 'grad_i']) The Terme_gravite keyword changes the numerical scheme used for the gravity source term. The default is grad_i, which is designed to remove spurious currents around the interface. In this case, the pressure field does not contain the hydrostatic part but only a jump across the interface. This scheme seems not to work very well in vef. The rho_g option uses the more traditional source term, equal to rho*g in the volume. In this case, the hydrostatic pressure is visible in the pressure field and the boundary conditions in pressure must be set accordingly. This model produces spurious currents in the vicinity of the fluid-fluid interfaces and with the immersed boundary conditions.
  • [equation_temperature_mpoint] (type: string) The equation_temperature_mpoint should be used in the case of liquid-vapor flow with phase-change (see the TRUST_ROOT/doc/TRUST/ft_chgt_phase.pdf written in French for more information about the model). The name of the temperature equation, defined with the convection_diffusion_temperature_ft_disc keyword, should be given.
  • [matrice_pression_invariante] (type: flag) This keyword is a shortcut to be used only when the flow is a single-phase one, with interface tracking only used for solid-fluid interfaces. In this peculiar case, the density of the fluid does not evolve during the computation and the pressure matrix does not need to be actuated at each time step.
  • [penalisation_forcage] (type: penalisation_forcage) This keyword is used to specify a strong formulation (value set to 0) or a weak formulation (value set to 1) for an imposed pressure boundary condition. The first formulation converges quicker and is stable in general cases except some rare cases (see Ecoulement_Neumann test case for example) where the second one should be used despite of its slow convergence.
  • [equation_temperature_mpoint_vapeur] (type: string) not_set
  • [mpoint_inactif_sur_qdm] (type: flag) not_set
  • [mpoint_vapeur_inactif_sur_qdm] (type: flag) not_set
  • [new_mass_source] (type: flag) Flag for localised computation of velocity jump based on interfacial area AI (advanced option)
  • [shift_secmem2] (type: flag) not_set
  • [interpol_indic_pour_di_dt] (type: string into ['interp_ai_based', 'interp_standard', 'interp_modifiee']) Specific interpolation of phase indicator function in VoF mass-preserving method (advanced option)
  • [outletcorrection_pour_di_dt] (type: string into ['correction_ghost_indic']) not_set
  • [boussinesq_approximation] (type: flag) not_set
  • [modele_turbulence] (type: modele_turbulence_hyd_deriv) Turbulence model for Navier-Stokes equations.
  • [solveur_pression] (type: solveur_sys_base) 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) value factor : 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 evualated 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) 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.
  • **[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) 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_ftd_ijk

Navier-Stokes equations.

Parameters:

  • multigrid_solver (type: multigrid_solver) not_set
  • [vitesse_entree] (type: float) Velocity to prescribe at inlet
  • [expression_vx_init] (type: string) initial field for x-velocity component (parser of x,y,z)
  • [expression_vy_init] (type: string) initial field for y-velocity component (parser of x,y,z)
  • [expression_vz_init] (type: string) initial field for z-velocity component (parser of x,y,z)
  • [expression_p_init] (type: string) initial pressure field (optional)
  • [velocity_diffusion_op] (type: string) Type of velocity diffusion scheme
  • [velocity_convection_op] (type: string) Type of velocity convection scheme
  • [fichier_reprise_vitesse] (type: string) not_set
  • [timestep_reprise_vitesse] (type: string) not_set
  • [disable_solveur_poisson] (type: flag) Disable pressure poisson solver
  • [disable_diffusion_qdm] (type: flag) Disable diffusion operator in momentum
  • [disable_convection_qdm] (type: flag) Disable convection operator in momentum
  • [frozen_velocity] (type: string) not_set
  • [velocity_reset] (type: string) not_set
  • [resolution_fluctuations] (type: flag) Disable pressure poisson solver
  • [harmonic_nu_in_diff_operator] (type: flag) Disable pressure poisson solver
  • [use_inv_rho_for_mass_solver_and_calculer_rho_v] (type: string) not_set
  • [use_inv_rho_in_poisson_solver] (type: flag) not_set
  • [diffusion_alternative] (type: string) not_set
  • [test_etapes_et_bilan] (type: flag) not_set
  • [ajout_init_a_reprise] (type: string) not_set
  • [improved_initial_pressure_guess] (type: string) not_set
  • [include_pressure_gradient_in_ustar] (type: string) not_set
  • [upstream_dir] (type: int) Direction to prescribe the velocity
  • [vitesse_upstream] (type: float) Velocity to prescribe at 'nb_diam_upstream_' before bubble 0.
  • [expression_vitesse_upstream] (type: string) Analytical expression to set the upstream velocity
  • [upstream_stencil] (type: int) Width on which the velocity is set
  • [nb_diam_upstream] (type: float) Number of bubble diameters upstream of bubble 0 to prescribe the velocity.
  • [nb_diam_ortho_shear_perio] (type: string) not_set
  • [vol_bulle_monodisperse] (type: string) not_set
  • [diam_bulle_monodisperse] (type: string) not_set
  • [coeff_evol_volume] (type: string) not_set
  • [vol_bulles] (type: string) not_set
  • [reprise_vap_velocity_tmoy] (type: string) not_set
  • [reprise_liq_velocity_tmoy] (type: string) not_set
  • [disable_source_interf] (type: flag) Disable computation of the interfacial source term
  • [harmonic_nu_in_calc_with_indicatrice] (type: flag) Disable pressure poisson solver
  • [refuse_patch_conservation_qdm_rk3_source_interf] (type: flag) experimental Keyword, not for use
  • [suppression_rejetons] (type: string) not_set
  • [p_seuil_max] (type: float) not_set, default 10000000
  • [p_seuil_min] (type: float) not_set, default -10000000
  • [coef_ammortissement] (type: float) not_set
  • [coef_immobilisation] (type: float) not_set
  • [expression_derivee_force] (type: string) expression of the time-derivative of the X-component of a source-term (see terme_force_ini for the initial value). terme_force_ini : initial value of the X-component of the source term (see expression_derivee_force for time evolution)
  • [terme_force_init] (type: string) not_set
  • [terme_force_init_z] (type: string) not_set
  • [correction_force] (type: string) not_set
  • [compute_force_init] (type: flag) not_set
  • [expression_variable_source_x] (type: string) not_set
  • [expression_variable_source_y] (type: string) not_set
  • [expression_variable_source_z] (type: string) not_set
  • [facteur_variable_source_init] (type: string) not_set
  • [expression_derivee_facteur_variable_source] (type: string) not_set
  • [expression_potential_phi] (type: string) parser to define phi and make a momentum source Nabla phi.
  • [forcage] (type: bloc_lecture) not_set
  • [corrections_qdm] (type: bloc_lecture) not_set
  • [coef_mean_force] (type: float) not_set
  • [coef_force_time_n] (type: float) not_set
  • [coef_rayon_force_rappel] (type: float) not_set
  • [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) 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_phase_field

Navier Stokes equation for the Phase Field problem.

Parameters:

  • [solveur_pression] (type: solveur_sys_base) 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) value factor : 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 evualated 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) 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.
  • **[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) 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_standard_sensibility

Resolution of Navier-Stokes sensitivity problem

Parameters:

  • state (type: bloc_lecture) Block to indicate the state problem. Between the braces, you must specify the key word 'pb_champ_evaluateur' then the name of the state problem and the velocity unknown Example: state { pb_champ_evaluateur pb_state velocity }
  • [uncertain_variable] (type: bloc_lecture) Block to indicate the name of the uncertain variable. Between the braces, you must specify the name of the unknown variable. Choice between velocity and mu. Example: uncertain_variable { velocity }
  • [polynomial_chaos] (type: float) It is the method that we will use to study the sensitivity of the Navier Stokes equation: if poly_chaos=0, the sensitivity will be treated by the standard sentivity method. If different than 0, it will be treated by the polynomial chaos method
  • [adjoint] (type: flag) A keyword to indicate that the adjoint Navier-Stokes equations will be solved
  • [solveur_pression] (type: solveur_sys_base) 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) value factor : 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 evualated 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) 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.
  • **[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) 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.

transport_2eq_base

Base class of the RANS transport equations in triocfd

Parameters:

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

transport_epsilon

The eps transport equation in bicephale (standard or realisable) k-eps model.

Parameters:

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

transport_interfaces_ft_disc

Interface tracking equation for Front-Tracking problem in the discontinuous version.

Parameters:

  • [conditions_initiales | initial_conditions] (type: bloc_lecture) The keyword conditions_initiales is used to define the shape of the initial interfaces through the zero level-set of a function, or through a mesh fichier_geom. Indicator function is set to 0, that is fluide0, where the function is negative; indicator function is set to 1, that is fluide1, where the function is positive; the interfaces are the level-set 0 of that function: conditions_initiales { fonction $(-((x-0.002)^2+(y-0.002)^2+z^2-(0.00125)^2))*((x-0.005)^2+(y-0.007)^2+z^2 (0.00150)^2))*(0.020-z))$ } In the above example, there are three interfaces: two bubbles in a liquid with a free surface. One bubble has a radius of 0.00125, i.e. 1.25 mm, and its center is {0.002, 0.002, 0.000}. The other bubble has a radius of 0.00150, i.e. 1.5 mm, and its center is {0.005, 0.007, 0.000}. The free surface is above the two bubble, at a level z=0.02. Additional feature in this block concerns the keywords ajout_phase0 and ajout_phase1. They can be used to simplify the composition of different interfaces. When using these keywords, the initial function defines the indicator function; ajout_phase0 and ajout_phase1 are used to modify this initial field. Each time ajout_phase0 is used, the field is untouched where the function is positive whereas the indicator field is set to 0 where the function is negative. The keyword ajout_phase1 has the symmetrical use, keeping the field value where the function is negative and setting the indicator field to 1 where the function is positive. The previous example can also be written: conditions_initiales { fonction z-0.020 , NL fonction ajout_phase1 $(x-0.002)^2+(y-0.002)^2+z^2-(0.00125)^2$ , fonction ajout_phase1 $(x-0.005)^2+(y-0.007)^2+z^2-(0.00150)^2$ }
  • [methode_transport] (type: methode_transport_deriv) Method of transport of interface.
  • [n_iterations_distance] (type: int) Keyword to specify the number or iterations requested for the smoothing process of computing the field corresponding to the signed distance to the interfaces and located at the center of the Eulerian elements. This smoothing is necessary when there are more Lagrangian nodes than Eulerian two-phase cells.
  • [maillage] (type: bloc_lecture) This optional block is used to specify that we want a Gnuplot drawing of the initial mesh. There is only one keyword, niveau_plot, that is used only to define if a Gnuplot drawing is active (value 1) or not active (value -1). By default, skipping the block will produce non Gnuplot drawing. This option is to be used only in a debug process.
  • [remaillage] (type: bloc_lecture_remaillage) This block is used to specify the operations that are used to keep the solid interfaces in a proper condition. The remaillage block only contains parameter\'s values.
  • [collisions] (type: bloc_lecture) This block is used to specify the operations that are used when a collision occurs between two parts of interfaces. When this occurs, it is necessary to build a new mesh that has locally a clear definition of what is inside and what is outside of the mesh. The collisions can either be active or inactive. If the collisions are active (highly recommended), a Juric level-set reconstruction method will be used to re-create the new mesh after each coalescence or breakup. An option Juric_local phase_continue N can be used to force the remeshing to impact only a local portion of the mesh, near the collision. The next line (type_remaillage) is used to state whose field will be used for the level-set computation. Main option is Juric, a remeshing that is compatible with parallel computing. When using Juric level-set remeshing, the source field (source_isovaleur) that is used to compute the level-sets is then defined. It can be either the indicator function (indicatrice), a choice which is the default one and the most robust, or a geometrical distance computed from the mesh at the beginning of the time step (fonction_distance), a choice that may be more accurate in specific situations. Type_remaillage can be either Juric or Thomas. When Thomas is used, it is an enhancement of the Juric remeshing algorithm designed to compensate for mass loss during remeshing. The mesh is always reconstructed with the indicator function (not with the distance function). After having reconstructed the mesh with the Juric algorithm, the difference between the old indicator function (before remeshing) and the new indicator function is computed. The differences occuring at a distance below or equal to N elements from the interface are summed up and used to move the interface in the normal direction. The displacement of the interface is such that the volume of each phase after displacement is equal to the volume of the phase before remeshing. N (default value 1) must be smaller than n_iterations_distance (suggested value: 2).
  • [methode_interpolation_v] (type: string into ['valeur_a_elem', 'vdf_lineaire', 'mean_volumic_velocity']) In this block, two keywords are possible for method to select the way the interpolation is performed. With the choice valeur_a_elem the speed of displacement of the nodes of the interfaces is the velocity at the center of the Eulerian element in which each node is located at the beginning of the time step. This choice is the default interpolation method. The choice VDF_lineaire is only available with a VDF discretization (VDF). In this case, the speed of displacement of the nodes of the interfaces is linearly interpolated on the 4 (in 2D) or the 6 (in 3D) Eulerian velocities closest the location of each node at the beginning of the time step. In peculiar situation, this choice may provide a better interpolated value. Of course, this choice is not available with a VEF discretization (VEFPreP1B).
  • [volume_impose_phase_1] (type: float) this keyword is used to specify the volume of one phase to keep the volume of the phases constant during the remeshing process. It is an alternate solution to trouble in mass conservation. This option is mainly realistic when only one inclusion of phase 1 is present in the domain. In most other situations, the iterations_correction_volume keyword seems easier to justify. The volume to be keep is in m3 and should agree with initial condition.
  • [parcours_interface] (type: parcours_interface) Parcours_interface allows you to configure the algorithm that computes the surface mesh to volume mesh intersection. This algorithm has some serious trouble when the surface mesh points coincide with some faces of the volume mesh. Effects are visible on the indicator function, in VDF when a plane interface coincides with a volume mesh surface. To overcome these problems, the keyword correction_parcours_thomas keyword can be used: it allows the algorithm to slightly move some mesh points. This algorithm is experimental and is NOT activated by default.
  • [interpolation_repere_local] (type: flag) Triggers a new transport algorithm for the interface: the velocity vector of lagrangian nodes is computed in the moving frame of reference of the center of each connex component, in such a way that relative displacements of nodes within a connex component of the lagrangian mesh are minimized, hence reducing the necessity of barycentering, smooting and local remeshing. Very efficient for bubbly flows.
  • [interpolation_champ_face] (type: interpolation_champ_face_deriv) It is possible to compute the imposed velocity for the solid-fluid interface by direct affectation (interpolation_scheme would be set to base) or by multi-linear interpolation (interpolation_scheme would be set to lineaire). The default value is base.
  • [n_iterations_interpolation_ibc] (type: int) Useful only with interpolation_champ_face positioned to lineaire. Set the value concerning the width of the region of the linear interpolation. For the Penalized Direct Forcing model, a value equals to 1 is enough.
  • [type_vitesse_imposee] (type: string into ['uniforme', 'analytique']) Useful only with interpolation_champ_face positioned to lineaire. Value of the keyword is uniforme (for an uniform solid-fluide interface\'s velocity, i.e. zero for instance) or analytique (for an analytic expression of the solid-fluide interface\'s velocity depending on the spatial coordinates). The default value is uniforme.
  • [nombre_facettes_retenues_par_cellule] (type: int) Keyword to specify the default number (3) of facets per cell used to describe the geometry of the solid-solid interface. This number should be increased if the geometry of the solid-solid interface is complex in each cell (eulerian mesh too coarse for example).
  • [seuil_convergence_uzawa] (type: float) Optional option to change the default value (10-8) of the threshold convergence for the Uzawa algorithm if used in the Penalized Direct Forcing model. Sometime, the value should be decreased to insure a better convergence to force equality between sequential and parallel results.
  • [nb_iteration_max_uzawa] (type: int) Optional option to change the default value (10-8) of the threshold convergence for the Uzawa algorithm if used in the Penalized Direct Forcing model. Sometime, the value should be decreased to insure a better convergence to force equality between sequential and parallel results.
  • [injecteur_interfaces] (type: string) not_set
  • [vitesse_imposee_regularisee] (type: int) not_set
  • [indic_faces_modifiee] (type: bloc_lecture) not_set
  • [distance_projete_faces] (type: string into ['simplifiee', 'initiale', 'modifiee']) not_set
  • [voflike_correction_volume] (type: int) flag (0 or 1). If activated, code will perfom the requested number of iterations for the correction process that can be used to keep the volume of the phases constant during the transport process.
  • [nb_lissage_correction_volume] (type: int) Select 0 or N, the number of smothing to apply to avoid potential spikes due to volume correction.
  • [nb_iterations_correction_volume] (type: int) to iterate on the volume correction until seuil is reached.
  • [type_indic_faces] (type: type_indic_faces_deriv) kind of interpolation to compute the face value of the phase indicator function (advanced option). Could be STANDARD, MODIFIEE or AI_BASED
  • [compute_particles_rms] (type: flag) not_set
  • [post_process_hydrodynamic_forces] (type: bloc_lecture) not_set
  • [collision_model_fpi] (type: collision_model_fpi_deriv) not_set
  • [suppression_sous_domaine] (type: string) not_set
  • [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.
  • [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) 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.

transport_k

The k transport equation in bicephale (standard or realisable) k-eps model.

Parameters:

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

transport_k_eps_base

Base equation for RANS k-eps model. Should not be used directly

Parameters:

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

transport_k_eps_non_std

Base Transport equation for non standard k-epsilon models

Parameters:

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

transport_k_keps

Transport equation for the k-epsilon model using the two-layer wall function.

Parameters:

  • [nb_couches] (type: int) Number of cells layer where the model is applied
  • [ystar_switch] (type: int) Value of y* to delimit the two layers
  • [conv_forcee] (type: string) Activate the forced convection model
  • [conv_nat] (type: string) Activate the natural convection model
  • [nut_switch] (type: int) Value of nut/nu which delimites the two layers
  • [impr] (type: flag) To print or not the 2 layers zone in the out file (default = false)
  • [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) 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.

transport_k_ou_eps

Transport equation for the uncoupled resolution version of the k-epsilon model.

Parameters:

  • [with_nu] (type: int) Add the kinematic viscosity in the diffusion operator
  • [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) 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.

transport_k_ou_eps_base

Base class of the transport equations for the uncoupled resolution version RANS models

Parameters:

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

transport_k_ou_eps_realisable

Transport equation for the uncoupled resolution version of the realisable k-epsilon model.

Parameters:

  • [with_nu] (type: int) Add the kinematic viscosity in the diffusion operator
  • [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) 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.

transport_marqueur_ft

not_set

Parameters:

  • [conditions_initiales | initial_conditions] (type: bloc_lecture) ne semble pas standard
  • [injection] (type: injection_marqueur) The keyword injection can be used to inject periodically during the calculation some other particles. The syntax for ensemble_points and proprietes_particles is the same than the initial conditions for the particles. The keyword t_debut_injection give the injection initial time (by default, given by t_debut_integration) and dt_injection gives the injection time period (by default given by dt_min).
  • [transformation_bulles] (type: bloc_lecture) This keyword will activate the transformation of an inclusion (small bubbles) into a particle. localisation gives the sub-zones (N number of sub-zones and their names) where the transformation may happen. The diameter size for the inclusion transformation is given by either diameter_min option, in this case the inclusion will be suppressed for a diameter less than diameter_size, either by the beta_transfo option, in this case the inclusion will be suppressed for a diameter less than diameter_size*cell_volume (cell_volume is the volume of the cell containing the inclusion). interface specifies the name of the inclusion interface and t_debut_transfo is the beginning time for the inclusion transformation operation (by default, it is t_debut_integr value) and dt_transfo is the period transformation (by default, it is dt_min value). In a two phase flow calculation, the particles will be suppressed when entring into the non marked phase
  • [phase_marquee] (type: int) Phase number giving the marked phase, where the particles are located (when they leave this phase, they are suppressed). By default, for a the two phase fluide, the particles are supposed to be into the phase 0 (liquid).
  • [methode_transport] (type: string into ['vitesse_interpolee', 'vitesse_particules']) Kind of transport method for the particles. With vitesse_interpolee, the velocity of the particles is the velocity a fluid interpolation velocity (option by default). With vitesse_particules, the velocity of the particules is governed by the resolution of a momentum equation for the particles.
  • [methode_couplage] (type: string into ['suivi', 'one_way_coupling', 'two_way_coupling']) Way of coupling between the fluid and the particles. By default, (keyword suivi), there is no interaction between both. With one_way_coupling keyword, the fluid act on the particles. With two_way_coupling keyword, besides, particles act on the fluid.
  • [nb_iterations] (type: int) Number of sub-timesteps to solve the momentum equation for the particles (1 per default).
  • [contribution_one_way] (type: int into [0, 1]) Activate (1, default) or not (0) the fluid forces on the particles when one_way_coupling or two_way_coupling coupling method is used.
  • [implicite] (type: int into [0, 1]) Impliciting (1) or not (0) the time scheme when weight added source term is used in the momentum equation
  • [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.
  • [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) 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.