Keywords introduced by TrioCFD that extend the Equation_base family. See Keywords derived from equation_base for the base-class documentation and the canonical keyword list.
convection_diffusion_concentration_ft_disc
not_set
Parameters:
- [constante_cinetique] (type: float) Constant kinetic rate of the reaction (laminar model 1 only).
- [equations_source_chimie] (type: list of str) Names of the reagent concentration equations.
- [modele_cinetique] (type: int) Reaction kinetic model index (1 to 4).
- [constante_cinetique_nu_t] (type: float) Kinetic constant scaling with the turbulent viscosity.
- [equation_nu_t] (type: string) Name of the hydraulic equation defining the turbulent viscosity.
- [domaine_sortie] (type: string) Name of the outlet sub-domain for the balance output.
- [phase] (type: int into [0, 1]) tells whether the scalar must be confined in phase 0 or in phase 1
- 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.
[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.
- Inherited from convection_diffusion_concentration:
- [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
- 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.
convection_diffusion_concentration_turbulent_ft_disc
equation_non_resolue
Parameters:
- [constante_cinetique] (type: float) This is the constant kinetic rate of the reaction and is used for the laminar model 1 only.
- [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.
- [constante_cinetique_nu_t] (type: float) Kinetic constant scaling with the turbulent viscosity.
- [equation_nu_t] (type: string) This specifies the name of the hydraulic equation used which defines the turbulent (basically SGS) viscosity.
- [domaine_sortie] (type: string) Name of the outlet sub-domain for the balance output.
- [phase] (type: int into [0, 1]) tells whether the scalar must be confined in phase 0 or in phase 1
- 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.
[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.
- Inherited from convection_diffusion_concentration_turbulent:
- [modele_turbulence] (type: modele_turbulence_scal_base — one of: null, prandtl, schmidt … (+1 more)) Turbulence model to be used in the constituent transport equations. The only model currently available is Schmidt.
- Inherited from convection_diffusion_concentration:
- [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
- 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.
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
- 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.
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.
- Inherited from convection_diffusion_concentration:
- [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
- 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.
convection_diffusion_temperature_ft_disc
not_set
Parameters:
- [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).
- [stencil_width] (type: int) distance in mesh elements over which the temperature field should be extrapolated in the opposite phase.
- [correction_courbure_ordre] (type: int) Order of the curvature correction for the temperature gradient.
- [correction_gradt_ordre] (type: int) not_set
- equation_interface (type: string) The name of the interface equation should be given.
- [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.
- equation_navier_stokes (type: string) The name of the Navier Stokes equation of the problem should be given.
- [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
- [divergence_free_velocity_extension] (type: flag) not_set
- [solveur_pression_fictive] (type: solveur_sys_base — one of: amg, amgx, cholesky … (+7 more)) not_set
- [bc_opening_pressure] (type: list of str) not_set
- [reinjection_tcl] (type: flag) Enable triple-contact-line (TCL) reinjection.
- [tempc] (type: float) Contact temperature for the TCL reinjection model.
- [rc_inject] (type: float) Injection radius for the TCL reinjection model.
- [thetac] (type: float) Contact angle for the TCL reinjection model.
- [rc_gridn] (type: int) Number of grid cells used by the TCL reinjection model.
- Inherited from convection_diffusion_temperature:
- [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.
- 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.
convection_diffusion_temperature_sensibility
Energy sensitivity equation (temperature diffusion convection)
Parameters:
[convection_sensibility | sensibility] (type: convection_deriv — one of: ale, amont, amont_old … (+20 more)) 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
- Inherited from convection_diffusion_temperature:
- [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.
- 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.
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
- 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_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:
- 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_ft_disc
Two-phase momentum balance equation.
Parameters:
- [boussinesq_approximation] (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
- [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.
- [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.
- [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.
- [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.
- [equation_temperature_mpoint] (type: string) The equation_temperature_mpoint should be used in the case of liquid-vapor flow with phase-change (see Front-Tracking VDF/VEF, section "Heat and mass transfer at the interfaces", for more information about the model). The name of the temperature equation, defined with the convection_diffusion_temperature_ft_disc keyword, should be given.
- [equation_temperature_mpoint_vapeur] (type: string) not_set
- [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.
- [equations_concentration_source_vortex] (type: list of str) Names of the concentration equations feeding the vortex source.
- [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.
- [mpoint_inactif_sur_qdm] (type: flag) not_set
- [mpoint_vapeur_inactif_sur_qdm] (type: flag) not_set
- [correction_courbure_ordre] (type: int) Order of the curvature correction in the interfacial velocity jump.
- [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
- [correction_diffusion_pch] (type: flag) Enable the diffusion correction for phase change.
- Inherited from navier_stokes_turbulent:
- [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_ftd_ijk
Navier-Stokes equations.
Parameters:
- multigrid_solver (type: multigrid_solver) not_set
- [vitesse_entree_dir] (type: int) Direction index of the imposed inlet velocity.
- [vitesse_entree_compo_to_force] (type: int) Velocity component forced at the inlet.
- [stencil_vitesse_entree] (type: int) Stencil width over which the inlet velocity is imposed.
- [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
- [projection_initiale] (type: flag) Flag to perform an initial projection of the velocity field.
- [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
- [use_harmonic_viscosity] (type: flag) Use a harmonic mean of the viscosity at the faces.
- [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_velocity_bubble_factor] (type: float) Bubble factor used by the upstream-velocity forcing.
- [upstream_velocity_bubble_factor_deriv] (type: float) Time derivative of the upstream-velocity bubble factor.
- [upstream_velocity_bubble_factor_integral] (type: float) Time integral of the upstream-velocity bubble factor.
- [velocity_bubble_scope] (type: float) Spatial extent over which the bubble velocity is evaluated.
- [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
- [upstream_velocity_measured] (type: flag) Flag to use the measured (rather than imposed) upstream velocity.
- [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
- [kp] (type: float) Proportional gain of the force/flow-rate controller.
- [kd] (type: float) Derivative gain of the force/flow-rate controller.
- [epaisseur_maille] (type: int) Cell-layer thickness used by the control algorithm.
- [no_sym] (type: int) Symmetry treatment option.
- [delta] (type: float) Characteristic length delta of the control box.
- [l] (type: float) Length L of the control box.
- [l_vol_controle] (type: float) Length of the control-volume box.
- [coef_mean_force] (type: float) not_set
- [coef_force_time_n] (type: float) not_set
- [coef_rayon_force_rappel] (type: float) not_set
- [correction_semi_locale_volume_bulle] (type: flag) Flag to enable the semi-local bubble-volume correction.
- 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_phase_field
Navier Stokes equation for the Phase Field problem.
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_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
- 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.
transport_2eq_base
Base class of the RANS transport equations in triocfd
Usable keywords (pick one): transport_k_epsilon, transport_k_epsilon_realisable, transport_k_keps, transport_k_omega, transport_k_ou_eps, transport_k_ou_eps_realisable
Parameters:
- 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.
transport_epsilon
- Defined in: src/TRAD2_turbulence.xd
The eps transport equation in bicephale (standard or realisable) k-eps model.
Parameters:
- 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.
transport_interfaces_ft_disc
Interface tracking equation for Front-Tracking problem in the discontinuous version.
Parameters:
- [indic_faces_modifiee] (type: bloc_lecture) not_set
- [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_transport (type: methode_transport_deriv — one of: loi_horaire, vitesse_imposee, vitesse_interpolee) 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.
- [iterations_correction_volume] (type: int) Number of volume-correction iterations (obsolete keyword).
- [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.
- [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.
- [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).
[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.
- [injecteur_interfaces] (type: string) not_set
- [suppression_sous_domaine] (type: string) not_set
- [sous_domaine_volume_impose] (type: string) Name of the sub-domain over which the phase volume is imposed.
- [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 — one of: lineaire, base) 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.
- [distance_ibc_faces] (type: string into ['initiale', 'modifiee']) Method to compute the IBC face distance.
- [distance_projete_faces] (type: string into ['simplifiee', 'initiale', 'modifiee']) not_set
- [temps_debut] (type: float) Physical time at which the interface transport starts.
- [vitesse_imposee_regularisee] (type: int) not_set
- [type_indic_faces] (type: type_indic_faces_deriv — one of: ai_based, modifiee, standard) kind of interpolation to compute the face value of the phase indicator function (advanced option). Could be STANDARD, MODIFIEE or AI_BASED
- [collision_model_fpi] (type: collision_model_ft_base — one of: collision_model_ft_arbitrary, collision_model_ft_sphere, collision_model_ft_spheroid) Collision model to use for the Fluid-Particle Interaction feature
- [compute_particles_rms] (type: flag) not_set
- [post_process_hydrodynamic_forces] (type: bloc_lecture) not_set
- Inherited from eqn_base:
[conditions_initiales | initial_conditions] (type: list of Cond_init_ft_base) Definition of the initial interface(s): a braced, comma-separated list of entries, each one of:
- fonction [ajout_phase0|ajout_phase1] EXPR
- fonction_ignorer_collision [ajout_phase0|ajout_phase1] EXPR
- fichier_geom { ... } (build from a geometry/mesh file)
- reprise FILE.xyz (resume the interfaces from a save file).
EXPR is a single-token analytical function of (x,y,z): each fonction entry adds an interface as its zero level-set, with the indicator set to 0 (fluide0) where EXPR<0 and to 1 (fluide1) where EXPR>0.
The first entry sets the indicator field; ajout_phase0/ajout_phase1 (only allowed from the second entry on) compose further interfaces onto it: ajout_phase0 keeps the field where its function is positive and sets the indicator to 0 where it is negative; ajout_phase1 is symmetric (keeps it where negative, sets it to 1 where positive).
Example - two bubbles below a free surface:
conditions_initiales {
fonction z-0.020 ,
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
}
- [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 — 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.
transport_k
- Defined in: src/TRAD2_turbulence.xd
The k transport equation in bicephale (standard or realisable) k-eps model.
Parameters:
- 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.
transport_k_eps_base
Base equation for RANS k-eps model. Should not be used directly
Usable keywords (pick one): transport_k_epsilon, transport_k_epsilon_realisable, transport_k_keps
Parameters:
- 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.
transport_k_eps_non_std_base
Base Transport equation for non standard k-epsilon models
Usable keywords (pick one): transport_k_keps
Parameters:
- 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.
transport_k_epsilon
The (k-eps) transport equation. To resume from a previous mixing length calculation, an external MED-format file containing reconstructed K and Epsilon quantities can be read (see fichier_ecriture_k_eps) thanks to the Champ_fonc_MED keyword.
Warning, When used with the Quasi-compressible model, k and eps should be viewed as rho k and rho epsilon when defining initial and boundary conditions or when visualizing values for k and eps. This bug will be fixed in a future version.
Parameters:
- [exit_on_negative_k_eps] (type: flag) Flag to exit (with postprocessing of fields) if a negative value is found for k or epsilon
- [exit_on_big_eps] (type: flag) Flag to exit (with postprocessing of fields) if an excessively big values of epsilon is found
- [disable_vef_mean_value_corrections] (type: flag) not_set
- [with_nu] (type: string into ['yes', 'no']) yes/no
- 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.
transport_k_epsilon_realisable
Realizable K-Epsilon Turbulence Model Transport Equations for K and Epsilon.
Parameters:
- [with_nu] (type: int) Add the kinematic viscosity in the diffusion operator
- 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.
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)
- [loi_2couches] (type: bloc_lecture) Two-layer near-wall model block.
- 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.
transport_k_omega
The (k-omega) transport equation.
Parameters:
- [exit_on_negative_k_omega] (type: flag) Flag to exit (with postprocessing of fields) if a negative value is found for k or omega
- [exit_on_big_omega] (type: flag) Flag to exit (with postprocessing of fields) if an excessively big values of omega are found
- [with_nu] (type: string into ['yes', 'no']) yes/no (default no)
- 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.
transport_k_omega_base
Base equation for RANS k-omega model. Should not be used directly
Usable keywords (pick one): transport_k_omega
Parameters:
- 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.
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
- 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.
transport_k_ou_eps_base
Base class of the transport equations for the uncoupled resolution version RANS models
Usable keywords (pick one): transport_k_ou_eps, transport_k_ou_eps_realisable
Parameters:
- 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.
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
- 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.
transport_marqueur_ft
not_set
Parameters:
- [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.
- [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).
- [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
- [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.
- [nb_iterations] (type: int) Number of sub-timesteps to solve the momentum equation for the particles (1 per default).
- [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
- Inherited from eqn_base:
- [conditions_initiales | initial_conditions] (type: bloc_lecture) ne semble pas standard
- [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 — 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.