Fluid Mechanics Numerical Methods Algorithms Boundary Conditions Supercomputers Turbulent Wind Tunnel Flight Tests
Fluid Mechanics Numerical Methods Algorithms Boundary Conditions Supercomputers Turbulent Wind Tunnel Flight Tests
- what is the name of the fundamental fluid dynamics equation describe the term involved in the
equation
is it possible to solve the equation directly by analytical methods
if yes, how?
If not, how one can solve the equation?
- What is laminar and turbulent flow, does the fundamental equation describe this two kinds of flow? If
not, what is the additional equation required to describe that?
- How can you describe turbulent behaviour in the sense of mathematical modelling?
What is turbulent viscosity and why we need the parameter?
What is RANS and LES modelling? What is the difference between both?
What is DNS? Give advantages and disadvantages
Computational fluid dynamics, usually abbreviated as CFD, is a branch of fluid mechanics that
uses numerical methods and algorithms to solve and analyze problems that involve fluid flows.
Computers are used to perform the calculations required to simulate the interaction of liquids and
gases with surfaces defined by boundary conditions. With high-speed supercomputers, better
solutions can be achieved. Ongoing research, however, yields software that improves the accuracy
and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of
such software is performed using a wind tunnel with the final validation coming in flight tests.
With advancements in micro-processor technology and the wide availability of powerful personal
computers, complex flow problems can now be solved on the desktop.
In addition to the computation of fluid flow, CFD techniques can also be used to compute the dynamic
behaviour of related processes including heat-transfer, solidification, sedimentation, and dispersion.
During pre-processing
The geometry (physical bounds) of the problem is defined.
The volume occupied by the fluid is divided into discrete cells (the mesh). The mesh may be
uniform or non uniform.
The physical modelling is defined – for example, the equations of motions + enthalpy +
radiation + species conservation
Boundary conditions are defined. This involves specifying the fluid behaviour and properties
at the boundaries of the problem. For transient problems, the initial conditions are also
defined.
The simulation is started and the equations are solved iteratively as a steady-state or transient.
Finally a postprocessor is used for the analysis and visualization of the resulting solution.
Laminar flow
Re < 2000
'low' velocity
Dye does not mix with water
Fluid particles move in straight lines
Simple mathematical analysis possible
Rare in practice in water systems.
Turbulent flow
Re > 4000
'high' velocity
Dye mixes rapidly and completely
Particle paths completely irregular
Average motion is in the direction of the flow
Cannot be seen by the naked eye
Changes/fluctuations are very difficult to detect. Must use laser.
Mathematical analysis very difficult - so experimental measures are used
Most common type of flow.
No so for that the new equation has been developed named as follows:
Reynolds-averaged Navier–Stokes
Main article: Reynolds-averaged Navier–Stokes equations
Large eddy simulation (LES) is a technique in which the smallest scales of the flow are
removed through a filtering operation, and their effect modelled using sub grid scale models.
This allows the largest and most important scales of the turbulence to be resolved, while
greatly reducing the computational cost incurred by the smallest scales. This method
requires greater computational resources than RANS methods, but is far cheaper than DNS.
Q1) Computational fluid dynamics, usually abbreviated as CFD, is a branch of fluid mechanics
that uses numerical methods and algorithms to solve and analyse problems that involve fluid
flows. Computers are used to perform the calculations required to simulate the interaction of
liquids and gases with surfaces defined by boundary conditions. With high-speed
supercomputers, better solutions can be achieved. Ongoing research, however, yields software
that improves the accuracy and speed of complex simulation scenarios such as transonic or
turbulent flows.
With advancements in micro-processor technology and the wide availability of powerful personal
computers, complex flow problems can now be solved on the desktop.
In addition to the computation of fluid flow, CFD techniques can also be used to compute the
dynamic behaviour of related processes including heat-transfer, solidification, sedimentation, and
dispersion.
Terms Involves
Analytical Soln of Navier stokes Pls Find
Q2
Laminar flow, sometimes known as streamline flow, occurs when a fluid flows in parallel layers, with
no disruption between the layers.[1] At low velocities the fluid tends to flow without lateral mixing, and
adjacent layers slide past one another like playing cards. There are no cross currents perpendicular
to the direction of flow, nor eddies or swirls of fluids.[2] In laminar flow the motion of the particles of
fluid is very orderly with all particles moving in straight lines parallel to the pipe walls.
• Re < 2000
• 'low' velocity
• Dye does not mix with water
• Fluid particles move in straight lines Simple mathematical analysis possible
Rare in practice in water systems.
In fluid dynamics, turbulence or turbulent flow is a flow regime characterized by chaotic and
stochastic property changes. This includes lowmomentum diffusion, high momentum convection,
and rapid variation of pressure and velocity in space and time.
• Re > 4000
• 'high' velocity
• Dye mixes rapidly and completely Particle paths completely irregular
• Average motion is in the direction of the flow
• Cannot be seen by the naked eye
• Changes/fluctuations are very difficult to detect. Must use laser. Mathematical
analysis very difficult - so experimental measures are used most common type of
flow.
• The numerical solution of the Navier–Stokes equations for turbulent flow is extremely
difficult, and due to the significantly different mixing-length scales that are involved in
turbulent flow, the stable solution of this requires such a fine mesh resolution that the
computational time becomes significantly infeasible for calculation (see Direct numerical
simulation). Attempts to solve turbulent flow using a laminar solver typically result in a
time-unsteady solution, which fails to converge appropriately. To counter this, time-
averaged equations such as the Reynolds-averaged Navier– Stokes equations (RANS),
supplemented with turbulence models, are used in practical computational fluid
dynamics (CFD) applications when modelling turbulent flows.
Q3
Turbulence models
• A turbulence model is a computational procedure to close the system of mean flow equations.
• For most engineering applications it is unnecessary to resolve the details of the turbulent
fluctuations.
• Turbulence models allow the calculation of the mean flow without first calculating the full time-
dependent flow field.
• We only need to know how turbulence affected the mean flow.
• In particular we need expressions for the Reynolds stresses.
• For a turbulence model to be useful it: must have wide applicability, be accurate, simple, and
economical to run.
A new quantity appears: the turbulent viscosity
• Its unit is the same as that of the molecular viscosity: Pa.s.
• It is also called the eddy viscosity.
• We can also define a kinematic turbulent viscosity . Its unit is m2/s.
• The turbulent viscosity is not homogeneous, i.e. it varies in space.
• It is, however, assumed to be isotropic. It is the same in all directions. This assumption is valid
for many flows, but not for all (e.g. flows with strong separation or swirl).
Eddy viscosity/ Turbulent Viscosity
In the study of turbulence in fluids, a common practical strategy for calculation is to ignore the
small-scale vortices (or eddies) in the motion and to calculate a large-scale motion with an eddy
viscosity that characterizes the transport and dissipation of energy in the smaller-scale flow (see
large eddy simulation).
RANS Modelling
Reynolds-averaged Navier-Stokes (RANS) equations are the oldest approach to turbulence
modelling. An ensemble version of the governing equations is solved, which introduces new
apparent stresses known as Reynolds stresses
RANS models can be divided into two broad approaches:
Boussinesq hypothesis
This method involves using an algebraic equation for the Reynolds stresses which include
and depending on the level of sophistication of the
determining the turbulent viscosity, model,
Solving transport equations for determining the turbulent kinetic energy and dissipation.
Models include k-ε (Launder and Spalding),[6] Mixing Length Model (Prandtl),[7] and Zero
Equation Model (CebeciandSmith).[7] The models available in this approach are often
referred to by the number of transport equations associated with the method. For example,
the Mixing Length model is a "Zero Equation" model because no transport equations are
solved; the k − ε is a "Two Equation" model because two transport equations (one for k and
one for ε) are solved.
Reynolds stress model (RSM)
This approach attempts to actually solve transport equations for the Reynolds stresses. This means
introduction of several transport equations for all the Reynolds stresses and hence this approach is
much more costly in CPU effort.
Large eddy simulation (LES) is a technique in which the smallest scales of the flow are
removed through a filtering operation, and their effect modelled using sub grid scale
models. This allows the largest and most important scales of the turbulence to be resolved,
while greatly reducing the computational cost incurred by the smallest scales. This method
requires greater computational resources than RANS methods, but is far cheaper than
DNS.
A direct numerical simulation (DNS) is a simulation in computational fluid dynamics in which
the Navier-Stokes equations are numerically solved without any turbulence model. This means
that the whole range of spatial and temporal scales of the turbulence must be resolved. All the
spatial scales of the turbulence must be resolved in the computational mesh, from the smallest
dissipative scales (Kolmogorov microscales), up to the integral scale L, associated with the
motions containing most of the kinetic energy.
One can estimate that the number of floating-point operations required to complete the simulation is
proportional to the number of mesh points and the number of time steps, and in conclusion, the
number of operations grows as Re3.
Therefore, the computational cost of DNS is very high, even at low Reynolds numbers. For
the Reynolds numbers encountered in most industrial applications, the computational
resources required by a DNS would exceed the capacity of the most powerful computers
currently available. However, direct numerical simulation is a useful tool in fundamental
research in turbulence. Using DNS it is possible to perform "numerical experiments", and
extract from them information difficult or impossible
The contributions of DNS to turbulence research in the last decade have been impressive
and the future seems bright. The greatest advantage of DNS is the stringent control it
provides over the flow being studied. It is expected that as flow geometries become more
complex, the numerical methods used in DNS will evolve. However, the significantly
higher numerical fidelity required by DNS will have to be kept in mind. It is expected that
use of non-conventional methodologies (e.g. multigrid) will lead to DNS solutions at an
affordable cost, and that development of nonlinear methods of analysis are likely to
prove very productive.
However, the main technical challenges of DNS remains the memory and computa-
tional speed requirements.
This makes it difficult to perform the DNS of turbulence with a higher
Reynolds number.
Q4
Q5
Consider a typical chemical reaction:
aA + bB → pP + qQ
The lowercase letters (a, b, p, and q) represent stoichiometric coefficients, while the capital
letters represent the reactants (A and B) and the products(P and Q).
According to IUPAC's Gold Book definition[1] the reaction rate r for a chemical reaction
occurring in a closed system under constant-volume conditions, without a build-up of reaction
intermediates, is defined as:
Where [X] denotes the concentration (Molarity, mol/L) of the substance X. (NOTE: Rate of
a reaction is always positive. '-' sign is present in the reactant involving terms because the reactant
concentration is decreasing.)
Q6
A grey body is defined as a body with constant emissivity over all wavelengths and
temperatures. Such an ideal body does not exist in practice but the assumption is a good
approximation for many objects used in engineering.
A black body is an idealized physical body that absorbs all incident electromagnetic
radiation. Because of this perfect absorptivity at all wavelengths, a black body is also the best
possible emitter of thermal radiation, which it radiates incandescently in a characteristic,
continuous spectrum that depends on the body's temperature. At Earth-ambient temperatures
this emission is in the infrared region of the electromagnetic spectrum and is not visible. The
object appears black, since it does not reflect or emit any visible light.
Emissivity
The emissivity ε indicates the ratio between the emitted energy of a real (not
black, grey) body and a black one.
Most earthly materials behave nearly as black bodies,i.e. are grey bodies (Figure 30) and emissivities
range from 0.9 to 1.0
"Grey body" emissivity ε ≤ 1 (P = εσT4)
Explicit and implicit methods are approaches used in numerical analysis for obtaining numerical
solutions of time-dependent ordinary and partial differential equations, as is required in computer
simulations of physical processes.
Explicit methods calculate the state of a system at a later time from the state of the system at the
current time, while implicit methods find a solution by solving an equation involving both the current
state of the system and the later one. Mathematically, if is the current system state and
is the state at the later time ( is a small time step), then, for an explicit method
to find
It is clear that implicit methods require an extra computation (solving the above equation),
and they can be much harder to implement. Implicit methods are used because many
problems arising in practice are stiff, for which the use of an explicit method requires
impractically small time steps to keep the error in the result bounded (see numerical
stability). For such problems, to achieve given accuracy, it takes much less computational
time to use an implicit method with larger time steps, even taking into account that one
needs to solve an equation of the form (1) at each time step.
Discretization Methods
Questionnaire on CFD calculation methods (Goerner, WS 2010)
Chapter 1
Not Important
Chapter 2
1.2. What are two main components, a simulation model of a technical system composed of?
Mathematical model and Physical model (based on probability density models, empirical models and /
or models physicochemical
2.2. What is the difference between system analysis and system identification? Which Method
is applied? What tools are available?
- System Analysis: output values unknown, input values and system characteristics are known which
are used to simulate the results of a process, such as Effect of a new combustion chamber
geometry on the exhaust gas composition; Tools: Flow Sheeting, CFD simulation - System ID:
System Properties unknown, input and output values are known. System is a black box, functional
relationships to be determined; Tools: fuzzy logic, artificial neural networks (ANN)
4.2. Which components and functions has an ANN and what is the analogy to the human
brain?
Human Brain: axon, neuron, synapses, dendrites (Function & order)
ANN: input variables, weighting factors, accounting (clearing) units, output variables
5.2. What is the typical procedure when using an ANN system identification? 1 Data
recording, 2 Data preparation, 3 Modelling, 4 Model training, 5 Model analysis, the model is
improved and recovery after 3, and extension of the data recovery after 1, 6 Optimization
9.2. What is the basic procedure for the CFD simulation of a system?
Problem Formulation → geometry creation and discretization (meshing) → select the applicable
models and algorithms → Numerical solution of the equation system → Display of the results or
simplified (( Pre-processing → Calculation → post-processing))
Chapter 3
1.3. Which criteria will be classified the different calculation models?
Type equation (Algebraic Eq., Differential Eq., Difference Eq., Integral Eq.), differential equation type
(Simple, Partial, Difference), single or multi-dimensional, homogeneous / inhomogeneous, stationary /
non-stationary (steady / unsteady)
Chapter 4
Not Important
Chapter 5 (conservation equations; Balance Equations)
1.5. What are the basic calculation approaches (frameworks) are available for the flow
simulation?
Eulerian framework (volume-based calculation approach), Lagrangian framework (streamline theory-
based calculation approach)
2.5. What is a conserved quantity? Give examples!
-Physical size, the quantity does not change within the limits of the system under consideration, i.e.,
remains constant (d / dx, y, z = 0, d / dt = 0)
Examples: energy, mass, amount of substance, pulse
3.5. Which major terms, the conservation equation is used to calculate the volume-based
model to together?
Time-dependent behaviour (transitory changes, volume integral of the time derivative of the
conserved quantity)
Convective-transport term (energy, momentum, and fabric) is the volume integral of the product of the
conserved quantity and velocity perpendicular to the surface, derived by a spatial direction, e.g. d / dx
Diffusion term: volume integral of the product of the transport coefficients and conserved, doubly
derived for the considered spatial coordinate
-Sources / Sink term: Volume integral over the sources / sink size
4.5. What is the flux density of a conserved quantity? Equations which can be used?
Flux = negative product of a transport coefficient and the condition of the size of the discharge space
coordinate, in which takes place the transport
Used: Newtonian flow law (flow), Fick's law (diffusion), Fourier law (heat conduction)
5.5. From which distribution functions and reactions can be set as a combustion model
together? Which sub-models are used to calculate?
Velocity, temperature, concentration and pressure field, gas volume and temperature, radiation,
convection and particle motion
Sub-models: flow model, heat transfer model reaction model
6.5. What can be made up a CFD model? (Honeycomb model)
Flow model (with turbulence model, particle turbulence interaction, two-phase model)
Heat transfer model (model with radiation, suspension temperature, radiation quantities) Reaction
model (heterogeneous and homogeneous gas phase reactions, combustion and gasification model)
Chapter 6
1.6. What coordinate systems are available for modelling? What are the respective coordinates
are defined?
Cartesian KS: Only axial coordination in 3 directions
Cylinder KS: One axial, radial (r to rmax) and azimuthal (0 to 2π)
Polar (spherical) KS: Arranged radial (r to rmax), azimuth and polar coordinate (0 to 2π), disposed
about 90 ° to each other
2.6. How should the conservation equation can be transformed for numerical computation?
The condition equation is a partial differential equation, and must be transformed into a differential
quotients.
3.6. What types of grid are available for the crosslinking of the geometric data set available?
Orthogonal grid (equidistant and variable) with 4 (2d) and 8 (3d) nodes per FE, unstructured (trigonal,
tetrahedral) network with three (2d) or 4 (3d) nodes per FE
Chapter 7
1.7. What are the basic types of differential equations? What are their characteristics?
Elliptic, parabolic and hyperbolic differential equation, characteristics = coefficient a j, bj, etc.
Need to be compared, where j stands for the respective iteration Elliptic
D.E.: All coefficients aj are non-zero and equal in magnitude
Hyperbolic D.E.: All coefficients for a j are equal to zero and a coefficient equal in magnitude
Parabolic D.E.: aj and bj equal to zero equal to zero
2.7. What are numerical methods to solve partial differential equation available? Finite
Difference Method: broken Taylor series development solution for all nodes, advantage: best
developed
Finite element method: minimizing errors between exact and approximate solution, weighted solution
for the entire space between the nodes, advantage: Networking is problem-adjusted Interpolation:
polynomials are used to approximate the dependent variable, (What is a dependent variable), FAM
advantage: Numerically stable (which causes the shorter computation times)
Spectral method: Fourier transform of the dependent variable in the frequency bands
3.7. What types of discretization are there?
-Temporal: explicit and implicit
-Spatial: forward, backward and central difference
4.7. What is the difference between an explicit and an implicit algorithm? Draw the iteration in
the space-time plane!
An explicit algorithm, the calculation of the value is i in a position for the next time step on the basis of
the local reference and change function for the current time step
-In an implicit algorithm, the value at location i for the next time step based on the values for i-1 and i
+1 are computed for the next time step (interpolated), i.e. it needs only the values for i-1 and i +1 are
calculated
Explicit method Implicit Method
5.7. Which method for the spatial discretization (spatial discretization) exist and how do they
work?
Central-Difference: The difference in values for i +1 and i-1, divided by twice the increment Forward-
Difference: The difference of the value of the next node (i +1) and the local value, divided by the
simple increment
Backward-Difference: The difference of the local value and the value at the previous node (i-1),
divided by the simple increment
6.7. What can you calculate by the Navier-Stokes equations and how are used?
Chapter 8
1.8. Which of the following equations are sub-tasks of a numerical simulation are used:
- Computational fluid dynamics: Navier-Stokes equation (momentum Eq., velocity and density change
per unit time = spatial variation of the stress plus volume forces acting on the fluid stress tensor
And
- balance of momentum:
- Volume forces: gravity and flow resistance
Chapter 9
1.9. What is the particularity of the velocity vector in the NS equation?
Speed transported other properties and itself is transported
Chapter 10
1.10. What is the general formulation of the conservation equation for a homogeneous reaction
2 order?
aA + bB → cC Σ iAi = 0 (can be expressed in matrix or vector form)
Matrix form:
Vector form:
10.10. Sketch the temperature and concentration gradients at the interface of an oil droplet
during the evaporation process and simultaneous evaporation and combustion!
11.10. Which pathways are important for the modelling of a heterogeneous gas-solid reaction
(e.g. coal combustion)?
Direct-particle reaction (Particle path way)
Pyrolysis via the intermediate tar (Gas liquid path way) with subsequent combustion of the gaseous
and solid reaction products
12.10. Which models and which parameters are to be considered for modelling the combustion
of coal?
Models: Pyrolysis model, gas and solid (char) reaction models, turbulence model, radiation model,
particle size distribution model, particle tracking
Parameters: Kinetic data for coal, balanced species, composition of the coal, physical properties of
coal (e.g. density, Cp), particle sizes, gas phase properties in the high temperature range, eddy
breakup parameters (vortex breakdown parameters), radiation parameters, turbulence parameters.
14.10. Draw a simplified reaction scheme for the combustion of coal! What time constants are
typical for mixing / reaction?
Time constants (in ms): Pyrolysis: 10(mixing) & 50 (reaction), Char combustion 200-400 (mixing)
/100-1000 (reaction)
15.10. Which model are exist for pyrolysis?
- Simple Reaction (Arrhenius or non-Arrhenius reaction)
- Various parallel reactions
- Various concurrent parallel reactions
- Complex reaction schemes (multiple consecutive parallel or competing parallel reactions)
- Reaction with secondary char formation
- Solomon two-stage model with primary and secondary pyrolysis
(Constant release rate, linear model, Arrhenius model, model of concurrent reactions e.g.
Kobayashi, DAE model, Component depending models e.g. Flashchain)
16.10. What are the reasons for the application of the simple reaction model?
-Heterogeneous composition of coal and associated uncertainty / variation in the indication of reaction
parameters
- Spatial and temporal discretization opportunities are not yet sufficient to resolve the pyrolysis
process at large boiler furnace
17.10. Which factors influence the release of volatiles (volatiles) in the pyrolysis?
-Heating time and heating rate, coal composition, particle size and size distribution, ambient
atmosphere
18.10. Draw the thermogravimetric curve of a pyrolysis reaction after Knill and after Kobayashi
model:
19.10. What are limited the rate constant k for the char combustion?
- Through (into CO on the grain surface and oxygen in the pores) transport processes
- Through a chemical reaction limited speed
- Processes in the boundary layer and in the pores running from diffusion controlled, processes
outside the boundary layer reaction controlled!
20.10. Name the most important reactions during coal gasification in order of their activation
energy!
C + 0.5 O2 → CO (lowest Ea, about 60 kJ / kmol)
C + CO2 → 2CO
C + H2O → CO + H2
CO2 + H2 → CO + H2O (highest Ea, about 327 MJ / kmol)
Chapter 11
No Important
Chapter 12
1.12. What proportions is composed of the total enthalpy change term?
- Enthalpy change by pressure change
- Change the sensible heat cp* T by change in temperature
- Enthalpy by changing the chemical composition
2.12. What proportions is composed of the species enthalpy hi?
- Chemically bound energy fh
- Sensible heat-cp* T
3.12. What proportions is composed of the transport term for the total enthalpy?
-Time change of the total enthalpy
- Convective transport net- - Net-diffusive transport
-Quellen-/Senkenterm
6.12. On what terms, the conservation equation for the total radiation intensity sets together?
- Transitory change (Temporal variation) of the radiation intensity
- Change of Iv along the radiation pathway S (Path-dependent change in the radiation intensity) - Sum
of reduction of radiation intensity by absorption and out-scattering (Radiation loss due to scattering
and absorption)
- Increase of radiation intensity by emission
- Change be in-scattering from all direction ( Change in the radiation intensity by interference from all
directions)
7.12. Which operations are to be considered in the calculation of the radiant flux (radiant
model as a "black-box")?
- Incoming radiation is absorbed, scattered or transmitted. Only the Transmitted and scattered
component is booked together with the interference from sources other than Escaping radiation.
9.12. What are the differences in the radiation and transmission behavior of real gases? - Band
radiation source: radiation emission, absorption and transmission only in specific wavelength ranges
(bands) Line radiation source: Radiation only certain wavelengths from linear (e.g. noble gases)
What are the types of discretization methods?
Time: Explicit and Implicit
Spatial: Forward, Backward and Central Difference
What is explicit and implicit? Which one of these methods is more stable and why?
Explicit methods calculate the state of a system at a later time from the state of the system at the
current time.
Implicit methods find a solution by solving an equation involving both the current state of the system
and the later one.
Implicit methods need extra computation and are harder to implement.
What is grey body, black body and define emissivity. Give value for emissivity of a grey body.
A grey body is a body with constant emissivity over all wavelengths and temperature.
A black body is an idealized physical body that absorbs all incident electromagnetic radiation.