Différences

Ci-dessous, les différences entre deux révisions de la page.

--- sunfluidh:new_numerical_methods_setup_namelist [2017/09/22 17:44] – créée yann
+++ sunfluidh:new_numerical_methods_setup_namelist [2018/05/29 15:06] (Version actuelle) – [Poisson_NumericalMethod] witko
@@ Ligne 1: / Ligne 1: @@
-===== Namelist "Numerical_Methods" =====
+===== Namelist "Numerical_Methods" (new version) =====
+__** Not for the release SUNFLUIDH_EDU**__
 <WRAP info>
-This data setup is devoted to the selection of the numerical methods and schemes used in order to solve the conservation equations for velocity components, temperature, species mass fractions and density (in particuliar cases) and Poisson equation for the pressure. Some parameter setting could be also considered in respect with the numerical method selected.
+This new data setup is also devoted to the selection of the numerical methods and schemes used in order to solve the conservation equations for velocity components, temperature, species mass fractions and density (in particuliar cases) and Poisson equation for the pressure. Some parameter setting could be also considered in respect with the numerical method selected.
+The data are divided in three groups in order to define:
+    * the numerical method applied for solving the conservation equations (for velocity, temperature, species, ...)
+    * the choice of advective or convective flux discretization (for 2nd order schemes only). The viscous, conductive or diffusive fluxes are always discretized with a centered 2nd order or 4th order scheme according to the previous choice.
+    * the numerical methods (and associated parameters) for solving the Poisson's equation according to the type of problem considered. Several methods are presented based on direct or iterative approaches.
 </WRAP>
 <note important>
-Two versions of the data setup exits. This is related to the newer release of the code.
+   * This version allows the user to select the numerical methods by means of characters strings instead of option numbers as in the old version [[Numerical_Methods_Setup_Namelist |Numerical_Methods (old version)]] .
+   * For the sake of clarity, the iterative method parameters used for solving the Poisson's equation are set in individual namelists :
+     * [[HomeData_PoissonSolver_Namelist |HomeData_PoissonSolver]] for methods directly implemented in the code ("homemade" development)
+     * [[HypreData_PoissonSolver_Namelist |HypreData_PoissonSolver]] for methods provided by the HYPRE library
+   * [[Numerical_Methods_Examples |Click here]] to get __some examples__ illustrating how to use the namelists "Numerical_Methods", "HomeData_PoissonSolver" or "HypreData_PoissonSolver"
 </note>
 ===== Full data set of the namelist =====
@@ Ligne 19: / Ligne 27: @@
 <note important>
-   * Following the problem treated, some numerical methods are better adapted than others, this point is precised for each available option.
+   * Depending on the problem, some numerical methods are better adapted than others, this point is precised for each available option.
-   * The partial diagonalization method (Poisson_NumericalMethod="Home-PartialDiagonalization") used for solving the Poisson equation is only permitted for separable problem.
+   * The partial diagonalization method (Poisson_NumericalMethod = "Home-PartialDiagonalization") used for solving the Poisson equation is only permitted for separable problems.
    * The HYPRE library solvers for solving the Poisson's equation are only available if the HYPRE library has been installed.
 </note>
@@ Ligne 30: / Ligne 39: @@
 ===== Solving the conservation equations =====
 -----
-====  NS_NumericalMethod (strictly equivalent to the variable "Numerical_Scheme" in the previous release) ====
+====  NS_NumericalMethod  ====
+(equivalent to "Numerical_Scheme" in the previous release, see [[Numerical_Methods_Setup_Namelist |Numerical_Methods (old version)]] )
    * Type: character string
    * Selection of the numerical scheme for solving the conservation equations :
-     * "BDF2-SchemeO2" : 2nd order Backward Differentiation Formula for the time discretization. Semi-implicit scheme on the viscous or diffusion terms. 2nd order spatial discretization (centered for the viscous/diffusion terms, different schemes are available for the convection/advection terms - see further). \\ Preferentially used for incompressible or low Mach-number flows without multi-species components (except two phase incompressible flows).
+     * __"BDF2-SchemeO2"__ : 2nd order Backward Differentiation Formula for the time discretization. Semi-implicit scheme on the viscous or diffusion terms. 2nd order spatial discretization (centered for the viscous/diffusion terms, different schemes are available for the convection/advection terms - see further). \\ Preferentially used for incompressible or low Mach-number flows without multi-species components (except two phase incompressible flows).
-     * "CN-SchemeO2" : 2nd order Crank-Nicolson scheme (semi-implicit scheme on the viscous or diffusion terms. 2nd order spatial discretization (centered for the viscous/diffusion terms, different schemes are available for the convection/advection terms - see further).\\ Preferentially used for incompressible or low Mach-number flows without multi-species components (except two phase incompressible flows).\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"CN-SchemeO2"__ : 2nd order Crank-Nicolson scheme (semi-implicit scheme on the viscous or diffusion terms. 2nd order spatial discretization (centered for the viscous/diffusion terms, different schemes are available for the convection/advection terms - see further).\\ Preferentially used for incompressible or low Mach-number flows without multi-species components (except two phase incompressible flows). __** Not for the release SUNFLUIDH_EDU**__ .
-     * "BDF1-SchemeO2" : 1st order Backward Differentiation Formula. Semi-implicit scheme on the viscous or diffusion terms. 2nd order spatial discretization (centered for the viscous/diffusion terms, different schemes are available for the convection/advection terms - see further).\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"BDF1-SchemeO2"__ : 1st order Backward Differentiation Formula. Semi-implicit scheme on the viscous or diffusion terms. 2nd order spatial discretization (centered for the viscous/diffusion terms, different schemes are available for the convection/advection terms - see further). __** Not for the release SUNFLUIDH_EDU**__ .
-     * "ExplicitPredCorO2-SchemeO2-Knio" : Explicit predictor-corrector scheme : 2nd order Adams Bashforth - Runge-Kutta scheme (Knio et al., JCP 1998). \\ This scheme is mainly used for low Mach Number flows with multi-species components.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"ExplicitPredCorO2-SchemeO2-Knio"__ : Explicit predictor-corrector scheme : 2nd order Adams Bashforth - Runge-Kutta scheme (Knio et al., JCP 1998). \\ This scheme is mainly used for low Mach Number flows with multi-species components. __** Not for the release SUNFLUIDH_EDU**__ .
-     * "ExplicitPredCorO2-SchemeO2-Njam" : 2nd order Explicit scheme based on Strang's operator-splitting with a multi-timestep procedure on the viscous/diffusion terms (2nd order Adams Bashforth - Runge-Kutta schemes) The density is explicitly solved from the differential formulation of the equation of state and the Temperature is deduced from the equation of state. The projection method is based on the Chorin's formulation by considering Poisson's Operator with constant coefficients (Njam et al., JCP 1999). \\ This scheme is mainly used for reactive flows under low Mach number hypothesis.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"ExplicitPredCorO2-SchemeO2-Njam"__ : 2nd order Explicit scheme based on Strang's operator-splitting with a multi-timestep procedure on the viscous/diffusion terms (2nd order Adams Bashforth - Runge-Kutta schemes) The density is explicitly solved from the differential formulation of the equation of state and the Temperature is deduced from the equation of state. The projection method is based on the Chorin's formulation by considering Poisson's Operator with constant coefficients (Njam et al., JCP 1999). \\ This scheme is mainly used for reactive flows under low Mach number hypothesis. __** Not for the release SUNFLUIDH_EDU**__ .
-     * "ExplicitPredCorO2-SchemeO2-Rixen" : 2nd order Explicit scheme based on Strang's operator-splitting with multistep  a multi-timestep procedure on the viscous/diffusion terms (2nd order Adams Bashforth - Runge-Kutta schemes). The temperature is explicitly solved from the enthalpy equation and the density is deduced from the equation of state. The projection method is based on the Goda's formulation by considering Poisson's Operator with variable coefficients (Rixin Yu et al., JCP 2012).\\__** Not for the release SUNFLUIDH_EDU**__ .
+     * __"ExplicitPredCorO2-SchemeO2-Rixen"__ : 2nd order Explicit scheme based on Strang's operator-splitting with multistep  a multi-timestep procedure on the viscous/diffusion terms (2nd order Adams Bashforth - Runge-Kutta schemes). The temperature is explicitly solved from the enthalpy equation and the density is deduced from the equation of state. The projection method is based on the Goda's formulation by considering Poisson's Operator with variable coefficients (Rixin Yu et al., JCP 2012). __** Not for the release SUNFLUIDH_EDU**__ .
-     * "ExplRKO3-CompactSchemeO4" : 4th order hermittian scheme in space + 3rd order explicit Runge-Kutta scheme in time. The projection method is based on the Goda's formulation by considering Poisson's Operator with constant coefficients (Knikker, ijnmf 2008, 2009). \\ This numerical scheme is used for incompressible or low Mach number flows with or without heat transfer.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"ExplRKO3-CompactSchemeO4"__ : 4th order hermittian scheme in space + 3rd order explicit Runge-Kutta scheme in time. The projection method is based on the Goda's formulation by considering Poisson's Operator with constant coefficients (Knikker, ijnmf 2008, 2009). \\ This numerical scheme is used for incompressible or low Mach number flows with or without heat transfer. __** Not for the release SUNFLUIDH_EDU**__ .
-     * "CN-CompactSchemeO4" : 4th order hermittian scheme in space + semi-implicit Crank-Nicolson scheme in time. The projection method is based on the Goda's formulation by considering Poisson's Operator with constant coefficients (Knikker, ijnmf 2008, 2009). \\ This numerical scheme is used for incompressible or low Mach number flows with or without heat transfer.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"CN-CompactSchemeO4"__ : 4th order hermittian scheme in space + semi-implicit Crank-Nicolson scheme in time. The projection method is based on the Goda's formulation by considering Poisson's Operator with constant coefficients (Knikker, ijnmf 2008, 2009). \\ This numerical scheme is used for incompressible or low Mach number flows with or without heat transfer. __** Not for the release SUNFLUIDH_EDU**__ .
-     * "RK3CN-CompactSchemeO4" : 4th order hermittian scheme in space + hybrid RK3/Crank-Nicolson scheme in time. The projection method is based on the Goda's formulation by considering Poisson's Operator with constant coefficients (Knikker, ijnmf 2008, 2009). \\ This numerical scheme is used for incompressible or low Mach number flows with or without heat transfer.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"RK3CN-CompactSchemeO4"__ : 4th order hermittian scheme in space + hybrid RK3/Crank-Nicolson scheme in time. The projection method is based on the Goda's formulation by considering Poisson's Operator with constant coefficients (Knikker, ijnmf 2008, 2009). \\ This numerical scheme is used for incompressible or low Mach number flows with or without heat transfer. __** Not for the release SUNFLUIDH_EDU**__ .
-     * "CN-SchemeO2-SpecialLowMachFlow": Crank-Nicolson semi-implicit scheme (as the option 2) with a predictor-corrector procedure for solving species mass fractions and temperature. This numerical scheme is used for low Mach number flows.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"CN-SchemeO2-SpecialLowMachFlow"__: Crank-Nicolson semi-implicit scheme (as the option 2) with a predictor-corrector procedure for solving species mass fractions and temperature. This numerical scheme is used for low Mach number flows. __** Not for the release SUNFLUIDH_EDU**__ .
-==== Convective_Flux_Discretization_Type ====
+==== MomentumConvection_Scheme ====
-   * Type : integer value
+(equivalent to "Convective_Flux_Discretization_Type" in the previous release, see [[Numerical_Methods_Setup_Namelist |Numerical_Methods (old version)]] )
+   * Type : character string
    * Selection of the 2nd order spatial discretization for the convection flux in the momentum equations. The options are :
-     * 1 : 2nd order centered scheme in the conservative formulation.
+     * __"Centered-O2-Conservative"__ : 2nd order centered scheme in the conservative formulation.
-     * 2 : 2nd order centered scheme in the advective formulation.
+     * __"Centered-O2-Convective"__   : 2nd order centered scheme in the advective formulation.
-     * 3 : 2nd order centered scheme in the skew-symemetric formulation (developer's test).
+     * __"Centered-O2-Skewsymmetric"__: 2nd order centered scheme in the skew-symemetric formulation (developer's test).
-     * 4 : not commented (develper's test)
+     * __"Upwind-O1-Conservative"__   : not commented (developer's test)
-     * 5 : QUICK scheme in the conservative formulation.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"Quick-O2-Conservative"__ : QUICK scheme in the conservative formulation. __** Not for the release SUNFLUIDH_EDU**__ .
-     * 6 : 2nd order Upwind scheme in the conservative formulation.\\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"Upwind-O2-Conservative"__ : 2nd order Upwind scheme in the conservative formulation. __** Not for the release SUNFLUIDH_EDU**__ .
-     * 7 : 2nd order ENO scheme in the conservative formulation (in progress). \\ __** Not for the release SUNFLUIDH_EDU**__ .
+     * __"Eno-O2-Conservative"__ : 2nd order ENO scheme in the conservative formulation (in progress). __** Not for the release SUNFLUIDH_EDU**__ .
-==== Temperature_Advective_Flux_Discretization_Type ====
+==== TemperatureAdvection_Scheme ====
-   * Type : integer value
+(equivalent to "Temperature_Advective_Flux_Discretization_Type" in the previous release)
+   * Type : character string
    * Selection of the 2nd order spatial discretization for the advection flux in the temperature/enthalpy equation. The options are the same as previously.
-==== species_Advective_Flux_Discretization_Type ====
+==== SpeciesAdvection_Scheme ====
-   * Type : integer value.  __** Not for the release SUNFLUIDH_EDU**__ .
+(equivalent to "Species_Advective_Flux_Discretization_Type" in the previous release)
+   * Type : character string.  __** Not for the release SUNFLUIDH_EDU**__ .
    * Selection of the 2nd order spatial discretization for the advection flux in the species equations. The options are the same as previously.
 ==== Explicit_Solving_of_Density ====
@@ Ligne 77: / Ligne 90: @@
 -----
-==== Numerical_Method_Poisson_Equation ====
+====  Poisson_NumericalMethod ====
-   * Type : Integer value
+(equivalent to "Numerical_Method_Poisson_Equation " in the previous release, see [[Numerical_Methods_Setup_Namelist |Numerical_Methods (old version)]] )
+   * Type : Character string
    * Selection of the numerical method for solving the Poisson equation in accordance to the projection method. The solution is the pressure time increment ($\Phi= P^{n+1}-P^{n}$, Goda's method) used to update the velocity field according to the principle of the projection method (When the numerical method of Njam et al. is used to solve the Navier-Stokes equations, the pressure is solved in place of its time increment (Chorin's method). The options are :
-     * 1 : Successive Over-Relaxed method (SOR) coupled with a multigrid method in order to accelerate the convergence. The matrix elements depend on the cell size only and the source term depends on the density, it is defined from the divergence of the momentum.
+     * __"Home-SORMultigrid-ConstantMatrixCoef"__ : Successive Over-Relaxed method (SOR) coupled with a nV-cycle multigrid method in order to accelerate the convergence. The matrix elements depend on the cell size only (constant elements). This method is directly implemented in the code ("homemade" development). No external library is required to use it. The associated parameters are set in the namelist [[HomeData_PoissonSolver_Namelist |HomeData_PoissonSolver]] .
-     * 2 : Poisson's operator is approximated by a Helmholtz's operator (experimental method proposed by J.L. Guermond) - For incompressible flow only.\\ __** Not for the release SUNFLUIDH_EDU**__
+     * __"Home-HelmholtzApproximation"__ : Poisson's operator (constant matrix elements) is approximated by a Helmholtz's operator (experimental method proposed by J.L. Guermond) - For incompressible flow only. This method is directly implemented in the code ("homemade" development). No external library is required to use it. No parameter setting is needed.
-     * 3: Partial diagonalisation of the Laplacian operator. BEWARE : the problem must be separable.\\ __** Not for the release SUNFLUIDH_EDU**__
+     * __"Home-PartialDiagonalization"__ : Partial diagonalisation of the Laplacian operator (constant elements). BEWARE : the problem must be separable. This method is directly implemented in the code ("homemade" development). External library Lapack and Blas is required to use it. The link must be specified in the makefile. No parameter setting is needed.
-     * 4 : SOR iterative method coupled with a multigrid method in order to accelerate the convergence. The matrix elements depend on the cell size and the density, the source term is defined from the divergence of the velocity.\\ __** Not for the release SUNFLUIDH_EDU**__
+     * __"Home-SORMultigrid-VariableMatrixCoef"__ : SOR iterative method coupled with a multigrid method in order to accelerate the convergence. The matrix elements depend on the cell size and the density, the source term is defined from the divergence of the velocity. This method is directly implemented in the code ("homemade" development). No external library is required to use it. The associated parameters are set in the namelist [[HomeData_PoissonSolver_Namelist |HomeData_PoissonSolver]] .
-     * 5 : Poisson's equation is solved by the HYPRE library tools. The matrix elements are constants. The parameters of the HYPRE library are set in the namelist [[ | ]]
+     * __"Hypre-ConstantMatrixCoef"__ : Poisson's equation is solved by the HYPRE library tools. The matrix elements are constants. The parameters of the HYPRE library are set in the namelist [[HypreData_PoissonSolver_Namelist |HypreData_PoissonSolver]] .
-==== Iterative_Method_Selection (only if the SOR method is used) ====
+     * __"Hypre-VariableMatrixCoef"__: Poisson's equation is solved by the HYPRE library tools. The matrix elements are non-constants. The parameters of the HYPRE library are set in the namelist [[HypreData_PoissonSolver_Namelist |HypreData_PoissonSolver]] .
-   * Type : Integer value
-   * Select the algorithm used for solving the Poisson equation from an SOR method. These options are :
-     * 1 : The SOR method is based on the Jacobi's algorithm
-     * 2 : The SOR method is based on the Gauss-Seidel algorithm
-     * 3 : The SOR method is based on the Gauss-Seidel algorithm with a two-color approach (Red-Black type). This method is used in the context of parallel simulations (MPI domain decomposition, Multithreading) in order to ensure an independence of the result in respect with the distribution of the calculation over the MPI processes or the threads. \\ __** Not for the release SUNFLUIDH_EDU**__
-==== Number_max_Grid (multigrid method only) ====
-   * Type : Integer value
-   * Maximum number of grid levels used in the multigrid method.
-==== Number_max_Cycle (multigrid method only) ====
-   * Type : Integer value
-   * Maximum number of V-cycles permitted in the multigrid method.
-==== Number_Iteration (only if the SOR method is used) ====
-   * Type : integer value
-   * Number of iterations performed by the SOR method in order to approach the solution on each grid level.
-==== Relaxation_Coefficient (only if the SOR method is used) ====
-   * Type : real value between one and two
-   * over-relaxation coefficient used in the SOR method.
-==== Convergence_Criterion (only if the SOR method is used) ====
-   * Type : real value
-   * The nV-cycle procedure of the multigrid method is stopped when the L2-norm of the Poisson equation's residu is lower than the criterion chosen.
-==== Off_Set_Poisson_Source_term ====
-   * Type : Boolean value. __** Not for the release SUNFLUIDH_EDU**__.
-   * The source term is shifted from its averaged value such as its integral value over the domain is zero. This could facilitate the convergence in some particuliar cases. this option must be used with caution. For expert users only.