Understanding CFD Simulation Process with Examples


Having already explained the background & evolution of Computational Fluid Dynamics (CFD) in the earlier blog Introduction to CFD, let us now try to understand the CFD simulation process with a few examples."CFD is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena by solving the mathematical equations which govern these processes using numerical methods"


It provides a qualitative and quantitative prediction of fluid flows by means of :

  • Mathematical modeling (partial differential equations)
  • Numerical methods (discretization and solution techniques)
  • Software tools (solvers, pre- and post-processing utilities)

Simulation Process :

To understand the simulation process and the steps involved in it let us consider an example of a flow through a pipe bend. The figure below gives series of the steps that would be involved in its analysis. For a fluid flow through a pipe bend we have the geometry built up, segregated into smaller fragments/ segments, called a mesh. With this mesh we actually define our probe-points where we want the analysis to be done. We then define the boundary conditions to get a unique solution solving it with a computer. The results obtained gives us a lot of data along these probe points that are then post-processed with visualization tools to analyse the results.

Thus CFD process in overall is a 3 step procedure :

  1. Pre Processing : This step consist of defining a geometry to define our domain of interest. The domain of interest is then divided into segments, called as mesh generation step and the problem is set-up defining the boundary conditions. Gridgen, CFD-GEOM, or ANSYS Workbench Environment & Modules, ANSYS ICEM CFD, TGrid etc., are some of the popular pre-processing softwares.
  2. Solver : Once the problem is set-up defining the boundary conditions we solve it with the software on the computer, (can also be done by hand-calculations, but would take long time). We have different popular commercial softwares available for this like Star-CD and Star CCM+ (CD-Adapco), FLUENT and CFX (ANSYS, Inc), GASP (Aerosoft, Inc), CFD++ (Metacomp Technologies) etc. Also there are free to use softwares like OpenFOAM, CFL3D, Typhon, OVERFLOW, Wind-US etc, all with different capabilities. These softwares are capable of solving the equations at every probe-point defined during the mesh generation step and also we can include additional models as required by the physics. The numerical methods are also defined at the this stage and we solve the whole problem.
  3. Post-processing : Once we get the results as values at our probe points we analyse them by means of color plots, contour plots, appropriate graphical representations & can generate reports. Tecplot 360, EnSight, FieldView, ParaView, ANSYS CFD-Post etc., are some of the popular post processing softwares.

How CFD Works ?

Now let us try to analyze a real life problem, with 2 examples discussed below.

1. Test Case: Fin-Tube Heat Exchanger

A Fin-Tube Heat Exchanger
 
The above image is of a fin-tube heat exchanger typically used for transferring heat in radiators in automobiles or in household applications (in cold countries where room heating is essential). In this we have cold/ hot fluid being poured through these tubes and other fluid (air/water) flowing over the tubes.
  Domain of Interest: Area between two fins
 
Now looking at the geometry, it can be seen that the fin tube heat exchanger is a cascade of a large number of fins attached to the tubes and thus seems to be a complex problem. However in CFD analysis of the fintube heat exchanger, the problem can be simplified to a great extent by identifying our domain of interest and considering only a small section of the it (as above). Simplifying assumptions are made in order to make the problem tractable (e.g., in this case: steady-state, incompressible, inviscid, two-dimensional). Also for solving the problem we consider the conservation of mass, momentum and energy as is required in the study. In pre-processing step we take the geometry and divide it into smaller fragments as in figure below, called meshing or the grid generation step.
{tex} \oint \rho \phi \vec v \cdot d\vec A = \oint \Gamma_{\phi} \nabla \phi \cdot d\vec A + \int_V S_{\phi} dV {/tex}
 
{tex}\sum_{f}^{N_{faces}} \rho_f \vec V_f \phi_f \cdot \vec{A_f} = \sum_{f} ^{N_{faces}} \Gamma_\phi (\nabla \phi)_n \cdot \vec{A_f} + S_{\phi} V \mbox{~} {/tex}
 
Schematic representation of domain discretization
 
Thus as the geometry is discretized so are the equations (as above) at each cell which on solving gives us the values that is obtained in the form of colorful contour plots using the visualization techniques that can give us a very good insight to locate the hot-spots, recirculation, & dead zones. So its not only the qualitative depiction of values that we generate but also the quantitative data (from figure we can see temperature varying from 464 K to 361 K) that can help us analyze the overall flow phenomena. Also if you see the regions between the fins as in figure below, it is clear that there is rapid fall in temperature from the base (hot region) to the top as is typically encountered with flue gases.
 

Temperature contours

2. Water flow over a tube bank

Consider the analysis of a cascade of tubes placed inside a domain across which a fluid (water) is flowing. The objective here is to compute average pressure drop and heat transfer per tube row. So it is a physical system in which we have a complicated setup of several cascade of tubes. To start with, we shall not proceed to solve the real problem but would try to simplify it taking a section for analysis, set-up the problem, first gain confidence over it and then if the computational resources are available and if time permits, proceed to implement on the real problem.

Geometry simplification

Assumptions

  • flow is two-dimensional, laminar, incompressible
  • flow approaching tube bank is steady with a known velocity
  • body forces due to gravity are negligible
  • flow is translationally periodic (i.e. geometry repeats itself)

As can be seen from the discretized geometry above, we have a fine mesh size at near the wall of the tubes and coarser at the other regions, so as to resolve the boundary layer flow at these regions. We have a whole system being modeled as we set up the problem on a software for example, Fluent as in figure below. In the solver we import the mesh, select the appropriate solver methodology, define operating conditions (no-slip, Qw or Tw at walls), initialize and iterate to get a converged solution.

 
ANSYS FLUENT setup
 
Once the problem is solved we can have results as in figure below. It shows contours of temperature around the tubes. As seen, the regions near the tube wall have re-circulation zones as a result of which there is a heat built-up (red color) also shown by color variations. So we see that we can have a very nice depiction of real life situation through a simple 2-D analysis.
 
 
Temperature contours around tubes (ANSYS FLUENT)
 
Also you can analyse a case wherein you have flow over a single tube with a obstruction-free region ahead, as in figure below. The figure shows the analysis for a typical phenomena called vortex shedding typically encountered in cases of flow around tubes because of which a familiar process called the von Karman vortices are generated. The tube encounters a lot of lift force that is sinusoidal in nature. The plot of time versus lift force clearly shows a sinusoidal nature in an analysis performed in ANSYS CFX.
Animation showing vortex shedding phenomena
 
 
 
Variation of lift with time (ANSYS CFX)

 

The Author

{module [317]}