While many simulations involving chemical species may require no special procedures during the solution process, you may find that one or more of the solution techniques noted in this section helps to accelerate the convergence or improve the stability of more complex simulations. The techniques outlined below may be of particular importance if your problem involves many species and/or chemical reactions, especially when modeling combusting flows.
Stability and Convergence in Reacting Flows
Obtaining a converged solution in a reacting flow can be difficult for a number of reasons. First, the impact of the chemical reaction on the basic flow pattern may be strong, leading to a model in which there is strong coupling between the mass/momentum balances and the species transport equations. This is especially true in combustion , where the reactions lead to a large heat release and subsequent density changes and large accelerations in the flow. All reacting systems have some degree of coupling, however, when the flow properties depend on the species concentrations. These coupling issues are best addressed by the use of a two-step solution process, as described below, and by the use of under-relaxation as described in Section 25.9.2.
A second convergence issue in reacting flows involves the magnitude of the reaction source term. When the FLUENT model involves very rapid reaction rates (reaction time scales are much faster than convection and diffusion time scales), the solution of the species transport equations becomes numerically difficult. Such systems are termed "stiff'' systems. Stiff systems with laminar chemistry can be solved using either the pressure-based solver with the Stiff Chemistry Solver option enabled, or the density-based solver (see Section 14.1.7). The laminar chemistry model may also be used for turbulent flames, where turbulence-chemistry interactions are neglected. However, for such flames, the EDC or PDF Transport models, which account for turbulence-chemistry interactions, may be a better choice.
Two-Step Solution Procedure (Cold Flow Simulation)
Solving a reacting flow as a two-step process can be a practical method for reaching a stable converged solution to your FLUENT problem. In this process, you begin by solving the flow, energy, and species equations with reactions disabled (the "cold-flow'', or unreacting flow). When the basic flow pattern has thus been established, you can re-enable the reactions and continue the calculation. The cold-flow solution provides a good starting solution for the calculation of the combusting system. This two-step approach to combustion modeling can be accomplished using the following procedure:
Define Models Species Transport & Reaction...
Solve Controls Solution...
One of the main reasons a combustion calculation can have difficulty converging is that large changes in temperature cause large changes in density, which can, in turn, cause instabilities in the flow solution. When you use the pressure-based solver, FLUENT allows you to under-relax the change in density to alleviate this difficulty. The default value for density under-relaxation is 1, but if you encounter convergence trouble you may wish to reduce this to a value between 0.5 and 1 (in the Solution Controls panel).
Ignition in Combustion Simulations
If you introduce fuel to an oxidant, spontaneous ignition does not occur unless the temperature of the mixture exceeds the activation energy threshold required to maintain combustion. This physical issue manifests itself in a FLUENT simulation as well. If you are using the laminar finite-rate, finite-rate/eddy-dissipation, EDC or PDF Transport model for turbulence-chemistry interaction, you have to supply an ignition source to initiate combustion. This ignition source may be a heated surface or inlet mass flow that heats the gas mixture above the required ignition temperature. Often, however, it is the equivalent of a spark: an initial solution state that causes combustion to proceed. You can supply this initial spark by patching a hot temperature into a region of the FLUENT model that contains a sufficient fuel/air mixture for ignition to occur.
Solve Initialize Patch...
Depending on the model, you may need to patch both the temperature and the fuel/ oxidant/product concentrations to produce ignition in your model. The initial patch has no impact on the final steady-state solution--no more than the location of a match determines the final flow pattern of the torch that it lights. See Section 25.14.2 for details about patching initial values.
Solution of Stiff Laminar Chemistry Systems
When modeling stiff laminar flames with the laminar finite-rate model, you can either use the pressure-based solver with the Stiff Chemistry Solver option enabled as seen in the Species Model panel (Figure 14.1.1), or the density-based solver.
When using the pressure-based solver for unsteady simulations, the Stiff Chemistry Solver option applies a fractional step algorithm. In the first fractional step, the chemistry in each cell is reacted at constant pressure for the flow time-step, using the ISAT integrator. In the second fractional step, the convection and diffusion terms are treated just as in a non-reacting simulation.
For steady simulations using the pressure-based solver, the Stiff Chemistry Solver option approximates the reaction rate in the species transport equation (see Equation 14.1-5) as,
where is some time-step that requires careful selection. Note that as tends to zero the approximation becomes exact but the stiff numerics will cause the pressure-based solver to diverge. On the other hand, as tends to infinity, the approximated reaction rate tends to zero and, while the numerical stiffness is alleviated, there is no reaction. In FLUENT, the default value for is set to one-tenth of the minimum convective or diffusive time-scale in the cell. This value was found to be sufficiently accurate and robust, although it can be modified via the solve/set/stiff-chemistry text command. ISAT is employed to integrate the stiff chemistry in Equation 14.1-34.
Details about the ISAT algorithm may be found in Sections 18.2.4 and Section 18.3.9. For efficient and accurate use of ISAT, a review of this section is highly recommended.
Choosing the density-based implicit solver can provide further solution stability by enabling the Stiff Chemistry Solver option. This option allows a larger stable Courant (CFL) number specification, although additional calculations are required to calculate the eigenvalues of the chemical Jacobian [ 393]. When enabling the stiff-chemistry solver, the following must be specified:
If the density-based explicit solver is used, then the stiff-chemistry solver has to be enabled via the text command:
solve set stiff-chemistry
You will be prompted to specify the following:
The default values of these parameters are applicable in most cases.
EDC Model Solution Procedure
Due to the high computational expense of the EDC model, it is recommended that you use the following procedure to obtain a solution using the pressure-based solver:
Define Models Species...
Define Boundary Conditions...
The Acceleration Factor can be set from 0 (slow but stable) to 1 (fast but unstable).