1. Computational Fluid Dynamics¶
1.2. What is a solution?¶
- The velocity field (vector)
- The associated pressure field (scalar)
1.3. Why do we Need CFD?¶
- There are very few known analytical solutions to the Navier Stokes Equations, e.g. when the convective term goes to zero, when there is flow in a pipe, flow between parallel plates and flow between concentric circles.
- Systems may be difficult to test through experimentation - e.g. experiment doesn’t allow us to see inside or the instrumentation is limited, or it’s dangerous to do the experiment.
- Faster and easier than experiment - CFD allows us to ask “what if” questions about a situation. Don’t have to build prototype.
- CFD can be used in animation for films
1.4. Components of a CFD Model¶
1.4.1. Mathematical Model¶
- A set of partial differential equations or integral-differential equations
- Associated boundary conditions
- The non linearity is a source of turbulence, vorticity, shock waves, combustion, multi-phase, bubble dynamics, evaporation, condensation. Therefore the model is associated with a target application - e.g. incompressible, inviscid, turbulent, 2D or 3D model
1.4.2. Discretization Method - the major part of CFD¶
- Method for approximating the partial differential equations or the integral-differential equations by a system of algebraic equations, i.e. we have a PDE \(\mathcal{L} [u(\underline{x})] = f(\underline{x})\) and we need to convert it to arithmetic \(A\underline{x} = \underline{b}\)
- The most important methods for doing this are:
- Finite Difference (1950s)
- Finite Element (1960s)
- Spectral Methods (1970s)
- Finite Volume (1980s)
- Boundary Element Method
- Particle Methods
- Discretization has two aspects:
- The Geometry (grid or mesh or particles) - gives us a vessel for the solution
- The Model - all mathematical operators converted into arithmetic operations on grid
1.4.3. Analyse Numerical Scheme¶
- All numerical schemes must satisfy certain conditions to be accepted:
- Consistency
- Stability
- Convergence
- Accuracy
1.4.4. Solve¶
- Obtain grid/point values of all flow variables
- Two situations:
- Time dependent \(\Rightarrow\) ODEs
- Steady \(\Rightarrow\) algebraic system of equations
- To solve these equations we require:
- Time integrators
- Linear solvers
1.4.5. Post Processing¶
- Visualization