Introduction to computational fluid dynamics
1 / 23

Introduction to Computational Fluid Dynamics - PowerPoint PPT Presentation

  • Updated On :

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)

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

PowerPoint Slideshow about 'Introduction to Computational Fluid Dynamics' - marinel

An Image/Link below is provided (as is) to download presentation

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 - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Introduction to computational fluid dynamics l.jpg

Introduction to Computational Fluid Dynamics

Course Notes (CFD 4)

Karthik Duraisamy

Department of Aerospace Engineering

University of Glasgow

Contents l.jpg

  • 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 l.jpg
Numerical Simulations : Estimates of cost of fixed wing calculation

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

Turbulent flow and the ns equations l.jpg
Turbulent flow and the NS Equations calculation

Characteristics of turbulent flow l.jpg
Characteristics of turbulent flow calculation

  • 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 l.jpg
Scales of turbulence calculation

  • 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.

Scales l.jpg
Scales calculation

Movies l.jpg
MOVIES calculation

Time and length scales l.jpg
Time and length scales calculation

  • 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 l.jpg
RANS and LES calculation

  • 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 l.jpg
In Summary calculation

  • 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 l.jpg
RANS Equations – What are they? calculation



Finally (Substituting & averaging)


Slide13 l.jpg
RANS calculation

  • 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 l.jpg
Eddy viscosity models calculation

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 l.jpg
Zero equation models calculation

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 l.jpg
2 Equation models calculation

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 l.jpg
Reynolds Stress Models calculation

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 l.jpg
Discretization calculation

  • 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 l.jpg
Wall bounded flows calculation

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 l.jpg
Turbulent Boundary layer calculation

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

Turbulent boundary layer scaling l.jpg
Turbulent Boundary Layer Scaling calculation

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

Slide22 l.jpg
LES calculation

  • 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.

Slide23 l.jpg
RANS calculation

  • 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