Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Solver Settings
5-1 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Outline
u Using the Solver
l Setting Solver Parameters
l Convergence
n Definition
n Monitoring
n Stability
n Accelerating Convergence
l Accuracy
n Grid Independence
n Adaption
u Appendix: Background
l Finite Volume Method
l Explicit vs. Implicit
l Segregated vs. Coupled
l Transient Solutions
5-2 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Solution Procedure Overview
u Solution Parameters Set the solution parameters
l Choosing the Solver
l Discretization Schemes Initialize the solution
u Initialization
u Convergence Enable the solution monitors of interest
l Monitoring Convergence
Calculate a solution
l Stability Modify solution
parameters or grid
n Setting Under-relaxation
n Setting Courant number Check for convergence
l Accelerating Convergence Yes
No
u Accuracy Check for accuracy
l Grid Independence No
Yes
l Adaption Stop
5-3 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Choosing a Solver
u Choices are Coupled-Implicit, Coupled-Explicit, or Segregated (Implicit)
u The Coupled solvers are recommended if a strong inter-dependence exists
between density, energy, momentum, and/or species.
l e.g., high speed compressible flow or finite-rate reaction modeled flows.
l In general, the Coupled-Implicit solver is recommended over the coupled-explicit
solver.
n Time required: Implicit solver runs roughly twice as fast.
n Memory required: Implicit solver requires roughly twice as much memory as coupled-
explicit or segregated-implicit solvers!
l The Coupled-Explicit solver should only be used for unsteady flows when the
characteristic time scale of problem is on same order as that of the acoustics.
n e.g., tracking transient shock wave
u The Segregated (implicit) solver is preferred in all other cases.
l Lower memory requirements than coupled-implicit solver.
l Segregated approach provides flexibility in solution procedure.
5-4 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Initialization
u Iterative procedure requires that all solution variables be initialized
before calculating a solution.
Solve → Initialize → Initialize...
l Realistic ‘guesses’ improves solution stability and accelerates convergence.
l In some cases, correct initial guess is required:
n Example: high temperature region to initiate chemical reaction.
u “Patch” values for individual
variables in certain regions.
Solve → Initialize → Patch...
l Free jet flows
(patch high velocity for jet)
l Combustion problems
(patch high temperature
for ignition)
5-8 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Convergence
u At convergence:
l All discrete conservation equations (momentum, energy, etc.) are
obeyed in all cells to a specified tolerance.
l Solution no longer changes with more iterations.
l Overall mass, momentum, energy, and scalar balances are obtained.
u Monitoring convergence with residuals:
l Generally, a decrease in residuals by 3 orders of magnitude indicates at
least qualitative convergence.
n Major flow features established.
l Scaled energy residual must decrease to 10-6 for segregated solver.
l Scaled species residual may need to decrease to 10-5 to achieve species
balance.
u Monitoring quantitative convergence:
l Monitor other variables for changes.
l Ensure that property conservation is satisfied.
5-10 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Convergence Monitors: Residuals
u Residual plots show when the residual values have reached the
specified tolerance.
Solve → Monitors → Residual...
All equations converged.
10-3
10-6
5-11 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Convergence Monitors: Forces/Surfaces
u In addition to residuals, you can also monitor:
l Lift, drag, or moment
Solve → Monitors → Force...
l Variables or functions (e.g., surface integrals)
at a boundary or any defined surface:
Solve → Monitors → Surface...
5-12 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Checking for Property Conservation
u In addition to monitoring residual and variable histories, you should
also check for overall heat and mass balances.
l At a minimum, the net imbalance should be less than 1% of smallest flux
through domain boundary.
Report → Fluxes...
5-13 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Decreasing the Convergence Tolerance
u If your monitors indicate that the solution is converged, but the
solution is still changing or has a large mass/heat imbalance:
l Reduce Convergence Criterion
or disable Check Convergence.
l Then calculate until solution
converges to the new tolerance.
5-14 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Convergence Difficulties
u Numerical instabilities can arise with an ill-posed problem, poor
quality mesh, and/or inappropriate solver settings.
l Exhibited as increasing (diverging) or “stuck” residuals.
l Diverging residuals imply increasing imbalance in conservation equations.
l Unconverged results can be misleading!
u Troubleshooting: Continuity equation convergence
l Ensure problem is well posed. trouble affects convergence of
all equations.
l Compute an initial solution with
a first-order discretization scheme.
l Decrease under-relaxation for
equations having convergence
trouble (segregated).
l Reduce Courant number (coupled).
l Re-mesh or refine grid with high
aspect ratio or highly skewed cells.
5-15 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Modifying Under-relaxation Factors
u Under-relaxation factor, α, is φ p = φ p ,old + α∆φ p
included to stabilize the iterative
process for the segregated solver.
u Use default under-relaxation factors
to start a calculation.
Solve → Controls → Solution...
u Decreasing under-relaxation for
momentum often aids convergence.
l Default settings are aggressive but
suitable for wide range of problems.
l ‘Appropriate’ settings best learned
from experience.
u For coupled solvers, under-relaxation factors for equations outside coupled
set are modified as in segregated solver.
5-16 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Modifying the Courant Number
u Courant number defines a ‘time
step’ size for steady-state problems.
l A transient term is included in the
coupled solver even for steady state
problems.
u For coupled-explicit solver:
l Stability constraints impose a
maximum limit on Courant number.
n Cannot be greater than 2.
s Default value is 1.
n Reduce Courant number when
having difficulty converging.
(CFL )∆x
u For coupled-implicit solver: ∆t =
u
l Courant number is not limited by stability constraints.
n Default is set to 5.
5-17 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Accelerating Convergence
u Convergence can be accelerated by:
l Supplying good initial conditions
n Starting from a previous solution.
l Increasing under-relaxation factors or Courant number
n Excessively high values can lead to instabilities.
n Recommend saving case and data files before continuing iterations.
l Controlling multigrid solver settings.
n Default settings define robust Multigrid solver and typically do not need
to be changed.
5-18 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Accuracy
u A converged solution is not necessarily an accurate one.
l Solve using 2nd order discretization.
l Ensure that solution is grid-independent.
n Use adaption to modify grid.
u If flow features do not seem reasonable:
l Reconsider physical models and boundary conditions.
l Examine grid and re-mesh.
5-21 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Mesh Quality and Solution Accuracy
u Numerical errors are associated with calculation of cell gradients and cell
face interpolations.
u These errors can be contained:
l Use higher order discretization schemes.
l Attempt to align grid with flow.
l Refine the mesh.
n Sufficient mesh density is necessary to resolve salient features of flow.
s Interpolation errors decrease with decreasing cell size.
n Minimize variations in cell size.
s Truncation error is minimized in a uniform mesh.
s Fluent provides capability to adapt mesh based on cell size variation.
n Minimize cell skewness and aspect ratio.
s In general, avoid aspect ratios higher than 5:1 (higher ratios allowed in b.l.).
s Optimal quad/hex cells have bounded angles of 90 degrees
s Optimal tri/tet cells are equilateral.
5-22 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Determining Grid Independence
u When solution no longer changes with further grid refinement, you
have a “grid-independent” solution.
u Procedure:
l Obtain new grid:
n Adapt
s Save original mesh before adapting.
– If you know where large gradients are expected, concentrate the
original grid in that region, e.g., boundary layer.
s Adapt grid.
– Data from original grid is automatically interpolated to finer grid.
n file → write-bc and file → read-bc facilitates set up of new problem
n file → reread-grid and File → Interpolate...
l Continue calculation to convergence.
l Compare results obtained w/different grids.
l Repeat procedure if necessary.
5-23 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Unsteady Flow Problems
u Transient solutions are possible with both segregated and coupled solvers.
l Solver iterates to convergence at each time level,
then advances automatically.
l Solution Initialization defines initial condition
and must be realistic.
u For segregated solver:
l Time step size, ∆t, is input in Iterate panel.
n ∆t must be small enough to resolve time
dependent features and to ensure convergence
within 20 iterations.
n May need to start solution with small ∆t.
l Number of time steps, N, is also required.
n N*∆t = total simulated time.
l To iterate without advancing time step, use ‘0’ time steps.
l PISO may aid in accelerating convergence for each time step.
5-24 © Fluent Inc. 1/29/02
Fluent User Services Center Introductory FLUENT Notes
FLUENT v6.0 Jan 2002
[Link] .com
Summary
u Solution procedure for the segregated and coupled solvers is the same:
l Calculate until you get a converged solution.
l Obtain second-order solution (recommended).
l Refine grid and recalculate until grid-independent solution is obtained.
u All solvers provide tools for judging and improving convergence and
ensuring stability.
u All solvers provide tools for checking and improving accuracy.
u Solution accuracy will depend on the appropriateness of the physical
models that you choose and the boundary conditions that you specify.
5-26 © Fluent Inc. 1/29/02