Introduction to Computational Fluid Dynamics - PowerPoint PPT Presentation

1 / 23

Introduction to Computational Fluid Dynamics. Course Notes (CFD 4). Karthik Duraisamy Department of Aerospace Engineering University of Glasgow. Contents. Introduction (1.5) Classification of PDE, Model equations (1.5) Finite difference methods:  Spatial discretization (3)

Related searches for Introduction to Computational Fluid Dynamics

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

Introduction to Computational Fluid Dynamics

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Introduction to Computational Fluid Dynamics

Course Notes (CFD 4)

Karthik Duraisamy

Department of Aerospace Engineering

University of Glasgow

Contents

• Introduction (1.5)

• Classification of PDE, Model equations (1.5)

• Finite difference methods:

 Spatial discretization (3)

 Temporal discretization (2)

 Convergence, Consistency, Stability (2)

• Grids/Boundary conditions (1)

• RANS Equations and Turbulence modeling (2)

• DNS/LES (1)

• Best practices in CFD (1)

• Case studies/Demonstrations (3)

(.) – Approximate number of lectures

Numerical Simulations : Estimates of cost of fixed wing calculation

Sample fixed wing of AR=10, Re=5e6

Turbulent flow and the NS Equations

http://web.mit.edu/fluids/www/Shapiro/ncfmf.html

Characteristics of turbulent flow

• Irregularity (Random and chaotic nature of flow)

• Increased exchange of momentum (Diffusivity - spreading rate of jets, boundary layers etc.)

• Large Reynolds numbers

• Dissipation of kinetic energy to internal energy

• Wide range of time and length scales

• Almost all practical flows are turbulent.

Scales of turbulence

• Turbulent flows are characterized by a wide range of length scales

• Think of turbulent flow as a collection of eddies of different sizes

• Eddies  time/length/velocity scales

• The largest “energy containing” eddies are of the order of the length scale of the object that generated turbulence in the first place (vortex shedding from a cylinder, boundary layer thickness)

• The smallest eddies are the ones where the energy is dissipated

• Roughly speaking, there is a “cascade” of energy from the large scale to the smallest scales  This happens because the large eddies interact with each other and breakdown into smaller eddies

• The smallest scales are called the “kolmogorov” scales (η). At these scales, the Reynolds number of the eddies is small enough that viscous effects become dominant and the energy is dissipated.

Time and length scales

• Kolmogorov (1941) showed that

L/η ~ (Re)3/4

T/τ ~ (Re)1/2

The problem with DNS of turbulent flows is that you have to simulate (or resolve) all these scales and in 3 dimensions

Therefore, number of points in each direction ~ (Re)3/4

Therefore, total number of points ~ (Re)9/4

Therefore, total number of time steps ~ (Re)1/2

Therefore, total number of operations ~ (Re)11/4

Remember, on top of this, you have to do hundreds of operations per point

RANS and LES

• DNS is obviously not feasible for high Reynolds number flows, therefore, LES and RANS to the rescue

• LES: Simulate (resolve) the large scales of turbulence and model the effect of the smaller scales (smaller scales are universal) by adding a subgrid viscosity

• RANS: Model all scales of turbulence

• Problems: LES still impractical for many flows, RANS is inaccurate in many flows

• But: One can get very useful and sometimes accurate, affordable engineering solutions with a good knowledge of the flow and the pros and cons of RANS and LES.

In Summary

• DNS : All scales of turbulence are resolved. Therefore, smallest grid size is of the order of the Kolmogorov scale η ~ L/(Re)3/4

• LES: Only the “Large” turbulent scales are

resolved. The “smaller” scales are modeled

• RANS : All the turbulent scales are modeled

In DNS, you just solve the Navier Stokes Equations

In LES, you solve a filtered version of the Navier-Stokes Equations along with another equation to represent the turbulent small scales

In RANS, you solve the averaged version of the Navier-Stokes equation along with another equation to represent all the turbulent scales [The extra equations are called turbulence model equations. Typical turbulence models are k-ε, k-ω, Spalart-Allmaras, etc]

RANS Equations – What are they?

NSE

Decompose

Finally (Substituting & averaging)

Where,

RANS

• Therefore, the objective in RANS turbulence modeling is to represent the unknown (Reynolds stress) in terms of the knowns (mean flow)

• Establishing such a relationship is called the closure model

• Types of closure models:

• - Algebraic models (Zero equation) [Mixing length model, Baldwin Lomax,Cebecci Smith]

• - One equation [K-model, Spalart Allmaras]

• - Two equation [K-e model, k-w model, SST]

• - Algebraic Reynolds Stress Models (ASBM)

• - Full Reynolds stress closures

Eddy viscosity models

Therefore, we determine eddy viscosity at each point and add it to the laminar viscosity

Closure problem is therefore, to determine eddy viscosity

Zero equation models

Lmix is dependant on the problem. Near a wall, it will be a function of distance to the wall and Reynolds number.

Models are reasonably accurate in attached flows

Very easy to code up

Poor correlation in off design conditions

No time history effects

Fully local.

2 Equation models

5 free constants

Simple, can work well in a variety of flows, history effects, more non-local

Reasonable results in many flows, but separation is a problem

Very diffusive, poor convergence near walls

Reynolds Stress Models

Solve directly for Reynolds stresses

Therefore, in theory more accurate than other methods – Production term is very important.

But still have to model some terms (other equations required for epsilon)

Very expensive, convergence is poor near the wall

Discretization

• Lets now look in little more detail about the mesh spacing required in various types of flows

• Turbulent flow

 Wall bounded flow

 Free-shear flow

Wall bounded flows

In reality, there exist eddies of a wide range of length scales in a boundary layer.

In DNS, we simulate (resolve) all these eddies

In LES we simulate some of the “significant” eddies and model the rest  Problem is that some of the “significant” eddies are still small

In RANS we model all the eddies

Turbulent Boundary layer

• Knowledge of boundary layer is essential to understand resolution requirements.

Turbulent Boundary Layer Scaling

u+ = u / U*, y+ = U* y / ν

U* =

Inner layer: u+ = y+

Outer layer u+ = 2.5 ln(y+)+5.5

Only for zero pressure gradient boundary layer on a flat plate

Inner layer

Outer layer

LES

• Problem with LES is that the near wall streaks are of size

L+ = 1000, W+=20, H+=30. Therefore to resolve this, we need

Δx+ = 100 (streamwise) , Δy+ = 1 (wall normal), Δz+ = 5 (spanwise)

• When you calculate requirements you find out that it is nearly as expensive as DNS!!

• Note, this is only for wall bounded flows

• Remember, in LES, you look to simulate those eddies that are significant. It so happens that in near-wall flows, eddies of this size are significant.

• In free shear flows, the “significant” eddies are not as small – loosely speaking.

RANS

• In RANS, you are only concerned about the wall-normal direction. In this direction, you need Δy+ = 1 – In the other two directions, grid sizes can be thousands of times larger and hence the savings.

• Spacing in other two directions is driven by accuracy considerations

• Remember, in DNS, only numerical error.

• In LES and RANS, both numerical and turbulence modeling errors