Modified on
25 Nov 2022 06:26 pm
Skill-Lync
In CFD, most of the fluid flow problems are non-linear in nature and deal with inherently unstable phenomena such as turbulence. The solution techniques use an iterative process to successively improve a solution until ‘convergence’ is reached.
Convergence is the limiting behavior observed in fluid dynamics. The exact solution to the iterative problem is unknown, but you want to be sufficiently close to the solution for a particular required level of accuracy. Convergence, therefore, does need to be associated with a requirement for a particular level of accuracy. This requirement depends upon the purpose to which the solution will be applied.
It is also important to keep in mind that a converged solution is not necessarily an accurate one. Many people called it “budgetary convergence,” which means your solution is converged when you can’t afford any more computer time.
Since the point at which the analysis is deemed converged is defined by the judgment of the analyst, users should have a solid understanding of when the analysis has reached its final solution. Typically, when assessing the convergence of a steady-state CFD analysis, the following three criteria need to be monitored as the analysis progresses:
The residual is one of the most fundamental measures of an iterative solution’s convergence, as it directly quantifies the error in the solution of the system of equations.
In a CFD analysis, the residual measures the local imbalance of a conserved variable in each control volume. Therefore, every cell in your model will have its own residual value for each of the equations being solved. A typical CFD simulation includes equations for momentum, pressure, and turbulence, and the solver performs iterative solutions of these equations. The residuals plot shows the difference between successive solutions of these equations.
The solutions are aggregated and normalized so that each equation gets represented by a single number, with all numbers scaled to a common range. All residuals get plotted on one graph (with Y-axis as log scale) to allow the CFD engineer to study interactions between the equations.
In an iterative numerical solution, the residual will never be exactly zero. However, the lower the residual value is, the more numerically accurate the solution. Each CFD code will have its own procedure for normalizing the solution residuals. It is best to check your code’s documentation for guidance on an appropriate criterion when judging convergence.
A good residuals plot has several characteristics to identify. The key element is decreasing lines and downward slopes. Residuals should always decrease. The sawtooth pattern in Figure 2 happens in unsteady simulations. Each spike represents a new timestep. Successive spikes should show decreasing peaks, or at least show each peak at the same height.
Figure 1: Residuals for Steady-state simulation.
Figure 2: Residuals for Unsteady state simulation.
For CFD, RMS residual levels of 1E-4 are considered to be loosely converged, levels of 1E-5 are considered to be well converged, and levels of 1E-6 are considered to be tightly converged. For complicated problems, however, it's not always possible to achieve residual levels as low as 1E-6 or even 1E-5.
As the word conservation suggests, the quantities from conservation equations (mass, momentum & energy) should conserve till the last iteration.
The solution imbalances will never be exactly zero because we are representing the physical system into an algebraic(numerical) form. However, the imbalances should be sufficiently small before considering the solution converged. As a good practice, we can aim for solution imbalances of below 1 % as a starting point but, highly sensitive analysis can demand tighter convergence.
Figure 3: Solution Imbalance of temperature.
As the figure shows, after the initial start-up period, the solution imbalances gradually decrease as the solution progresses. In general, the small values of solution residuals lead to small solution imbalances.
There are a few exceptions to the above statement that sometimes small residuals can have large imbalances. For example, in CHT simulation conduction timescales are much larger than flow timescales.
In an iterative solution, the calculation process starts with initial estimations of quantities (pressure, velocity, temperature, etc.). In every iteration, quantities are updated with respect to the previous iteration's results. In a converged steady-state solution, one should expect to see those quantities come to an expected level and don’t change anymore.
Monitoring integrated quantities such as force, drag, or average temperature can help the user judge when his or her analysis has reached this point. There might be instances where one quantity may have reached its final point, yet the other still shows signs of fluctuations. Thus, considering this point as the final solution can be misleading. It is better to iterate to the point where values flatten out, in order to consider the solution converged.
Figure 4: Temperature over the number of iterations.
We can see in the above figure that the temperature after 25 iterations is within just a few percent of its final value. Once the monitor point values have "flattened out", we can safely assume the solution is converged.
The residuals are more than just a tool for judging convergence. They also provide feedback for troubleshooting bad simulations. The pattern of the residuals provides clues as to why a simulation fails.
As CFD engineers, we do not have the luxury of judging a binary system: pass/fail. Convergence is not an exact science, so we should rely on evidence from all the above-mentioned criteria. As discussed, we have three tools to judge convergence which are residuals, monitors, and flow patterns. We need to combine all three tools and build a composite picture to reliably judge convergence in any simulation.
Author
Navin Baskar
Author
Skill-Lync
Subscribe to Our Free Newsletter
Continue Reading
Related Blogs
Learn how to render a shock-tube-simulation and how to work on similar projects after enrolling into anyone of Skill-Lync's CAE courses.
10 May 2020
In this blog, read how to design the frontal BIW enclosure of a car (Bonnet) and learn how Skill-Lync Master's Program in Automotive Design using CATIA V5 will help you get employed as a design engineer.
10 May 2020
Tetrahedral is a four- nodded solid element that can be generated through the tria element by creating a volume and also through the existing volume of the geometry. These elements are used where the geometry has high thickness and complexity. The image attached below is a representation of a Tetra element. The Tetra element will have 4 triangular faces with four nodes joining them together
02 Aug 2022
A connector is a mechanism that specifies how an object (vertex, edge, or face) is connected to another object or the ground. By often simulating the desired behaviour without having to build the precise shape or specify contact circumstances, connectors make modeling simpler.
03 Aug 2022
One of the most crucial processes in carrying out an accurate simulation using FEA is meshing. A mesh is composed of elements that have nodes—coordinate positions in space that might change depending on the element type—that symbolise the geometry's shape.
04 Aug 2022
Author
Skill-Lync
Subscribe to Our Free Newsletter
Continue Reading
Related Blogs
Learn how to render a shock-tube-simulation and how to work on similar projects after enrolling into anyone of Skill-Lync's CAE courses.
10 May 2020
In this blog, read how to design the frontal BIW enclosure of a car (Bonnet) and learn how Skill-Lync Master's Program in Automotive Design using CATIA V5 will help you get employed as a design engineer.
10 May 2020
Tetrahedral is a four- nodded solid element that can be generated through the tria element by creating a volume and also through the existing volume of the geometry. These elements are used where the geometry has high thickness and complexity. The image attached below is a representation of a Tetra element. The Tetra element will have 4 triangular faces with four nodes joining them together
02 Aug 2022
A connector is a mechanism that specifies how an object (vertex, edge, or face) is connected to another object or the ground. By often simulating the desired behaviour without having to build the precise shape or specify contact circumstances, connectors make modeling simpler.
03 Aug 2022
One of the most crucial processes in carrying out an accurate simulation using FEA is meshing. A mesh is composed of elements that have nodes—coordinate positions in space that might change depending on the element type—that symbolise the geometry's shape.
04 Aug 2022
Related Courses