Numerical Solution of Partial Differential Equations governing compressible flows

Size: px
Start display at page:

Download "Numerical Solution of Partial Differential Equations governing compressible flows"

Transcription

1 Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore January 2009 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

2 Description of flow Gas treated as a continuum At every point of space occupied by gas, define density ρ(x, y, z, t) = lim volume 0 mass volume Velocity vector u(x, y, z, t) v(x, y, z, t) w(x, y, z, t) Pressure p(x, y, z, t) = Force lim area 0 area Praveen. C (TIFR-CAM) CFD Jan 31, / 45

3 Description of flow Temperature: T (x, y, z, t) Thermodynamics: relates ρ, p, T Five quantities: ρ, u, v, w, p completely specify the state of the fluid Require five equations to describe the evolution of the state Praveen. C (TIFR-CAM) CFD Jan 31, / 45

4 Divergence theorem In Cartesian coordinates, q = (u, v, w) q = u x + v y + w z Divergence theorem dv n V V qdv = V q ˆndS Praveen. C (TIFR-CAM) CFD Jan 31, / 45

5 Conservation laws Basic laws of physics are conservation laws: mass, momentum, energy, charge Control volume n Streamlines Rate of change of φ in V = -(net flux across V ) φdv = F ˆndS t V V Praveen. C (TIFR-CAM) CFD Jan 31, / 45

6 Conservation laws If F is differentiable, use divergence theorem φdv = t F dv or V [ ] φ t + F dv = 0, V V φ t + F = 0 In Cartesian coordinates, F = (f, g, h) φ t + f x + g y + h z = 0 for every control volume V Praveen. C (TIFR-CAM) CFD Jan 31, / 45

7 Mass conservation equation φ = ρ, F = ρ q F ˆndS = amount of mass flowing across ds per unit time Using divergence theorem ρdv + t or ρdv + ρ q ˆndS = 0 t V V V V V (ρ q)ds = 0 [ ] ρ t + (ρ q) dv = 0 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

8 Mass conservation equation Differential form ρ t + (ρ q) = 0 In Cartesian coordinates, q = (u, v, w) ρ t + x (ρu) + y (ρv) + z (ρw) = 0 Non-conservative form ρ t + u ρ x + v ρ y + w ρ ( u z + ρ x + v y + w ) = 0 z or in vector notation ρ + q ρ + ρ q = 0 t Praveen. C (TIFR-CAM) CFD Jan 31, / 45

9 Momentum conservation Newton s law of motion Rate of change of = Total body force + Total surface force momentum Body forces like gravity Surface forces: pressure and viscous shearing x, y, z momentum equations t (ρu)+ x (p+ρu2 )+ y (ρuv)+ z (ρuw) = x τ xx+ y τ xy+ z τ xz t (ρv)+ x (ρvu)+ y (p+ρv2 )+ z (ρvw) = x τ yx+ y τ yy+ z τ yz t (ρw)+ x (ρwu)+ y (ρwv)+ z (p+ρw2 ) = x τ zx+ y τ zy+ z τ zz Praveen. C (TIFR-CAM) CFD Jan 31, / 45

10 Energy conservation Energy E = internal + kinetic + etc. Rate of Flux of Work done by Transfer of change = energy due + surface and + heat by of energy to flow body forces conduction Energy equation E t + x [(E + p)u] + y [(E + p)v] + [(E + p)w] = etc. z Praveen. C (TIFR-CAM) CFD Jan 31, / 45

11 Navier-Stokes equations Mass: Momentum: Energy: t ρ t + (ρ q) = 0 (ρ q) + (ρ q q) + p = τ E t + (E + p) q = ( q τ) + Q Praveen. C (TIFR-CAM) CFD Jan 31, / 45

12 Convection and diffusion Mass: Momentum: Energy: t ρ t + (ρ q) = 0 (ρ q) + (ρ q q) + p = τ E t + (E + p) q = ( q τ) + Q Praveen. C (TIFR-CAM) CFD Jan 31, / 45

13 Convection and diffusion Continuity equation in non-conservative form ρ t + u ρ x + v ρ y + w ρ ( u z + ρ x + v y + w ) z = Linear convection equation x momentum equation u t + u u x +... = µ 2 u x = Non-linear convection and diffusion phenomena = 0 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

14 Differential and integral form Differential form φ t + F = 0 This is the basis of Finite Difference Methods (conserved quantity) + divergence(flux) = 0 t Integral form φdv = F ˆndS t V V This is the basis of Finite Volume Methods Praveen. C (TIFR-CAM) CFD Jan 31, / 45

15 Flow with shocks Praveen. C (TIFR-CAM) CFD Jan 31, / 45

16 Flow with shocks X-15 model at mach = 3.5 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

17 Computational domain (Josy P. P., CTFD, NAL) Praveen. C (TIFR-CAM) CFD Jan 31, / 45

18 Grids and finite volumes Praveen. C (TIFR-CAM) CFD Jan 31, / 45

19 Structured grid Praveen. C (TIFR-CAM) CFD Jan 31, / 45

20 Unstructured grid Praveen. C (TIFR-CAM) CFD Jan 31, / 45

21 Unstructured hybrid grid Praveen. C (TIFR-CAM) CFD Jan 31, / 45

22 Unstructured hybrid grid Praveen. C (TIFR-CAM) CFD Jan 31, / 45

23 Cartesian grid Saras (Josy) Praveen. C (TIFR-CAM) CFD Jan 31, / 45

24 Block-structured grid Fighter aircraft (Nair) Praveen. C (TIFR-CAM) CFD Jan 31, / 45

25 Forces on an airplane Lif t = W eight Drag = T hrust Praveen. C (TIFR-CAM) CFD Jan 31, / 45

26 Wing and Airfoil Praveen. C (TIFR-CAM) CFD Jan 31, / 45

27 Effect of Mach number: M = fluid speed/sound speed Praveen. C (TIFR-CAM) CFD Jan 31, / 45

28 CFD process Given PDE, initial condition, boundary condition Define computational domain Divide computational domain into a grid Solve using a numerical scheme Compute quantities of interest lift force and drag force Praveen. C (TIFR-CAM) CFD Jan 31, / 45

29 Flow over NACA0012 airfoil: M = 0.63, α = 2 Pressure Praveen. C (TIFR-CAM) CFD Jan 31, / 45

30 Flow over NACA0012 airfoil: M = 0.85, α = 1 Pressure 5 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

31 Flow over Suddhoo-Hall airfoil Pressure contours Praveen. C (TIFR-CAM) CFD Jan 31, / 45

32 Flow over Suddhoo-Hall airfoil Pressure on airfoil 5 KMM Exact Cp x Praveen. C (TIFR-CAM) CFD Jan 31, / 45

33 Flow through scramjet intake Density Praveen. C (TIFR-CAM) CFD Jan 31, / 45

34 CFD example: Falcon Mach = 0.85 Grid: nodes, Solver: Num3sis Praveen. C (TIFR-CAM) CFD Jan 31, / 45

35 Applications of numerical simulation Forward problem Given initial and boundary conditions, shape of computational domain solve for the flow Inverse problem Given the flow solution compute initial and/or boundary condition, and/or shape of computational domain Optimization problem Given some measure of system performance compute initial and/or boundary condition, and/or shape of computational domain to maximize the performance measure = Optimization of systems governed by PDE Praveen. C (TIFR-CAM) CFD Jan 31, / 45

36 Example of optimization: RAE2822 Initial shape Solver: euler2d Praveen. C (TIFR-CAM) CFD Jan 31, / 45

37 Example of optimization: RAE2822 Optimized shape Optimizer: Torczon Simplex, 20 Hicks-Henne parameters Praveen. C (TIFR-CAM) CFD Jan 31, / 45

38 Example of optimization: RAE2822 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

39 Shape optimization Very small shape modifications: cannot be obtained using trial and error methods Transonic wings: possible to obtain shock-free flow with small shape changes 1% - 3% reduction in drag = significant fuel savings Many other applications Turbomachinery Ship hull design Automobiles Chemical industry Praveen. C (TIFR-CAM) CFD Jan 31, / 45

40 Wing optimization Initial shape Praveen. C (TIFR-CAM) CFD Jan 31, / 45

41 Wing optimization Optimized shape Praveen. C (TIFR-CAM) CFD Jan 31, / 45

42 Shape optimization problem Cost function J = J(S, u) u satisfies the flow equations (Euler or Navier-Stokes) N(S, u) = 0 u depends on shape S Minimization problem min S S ad J(S, u), s.t. N(S, u) = 0 S ad = Set of admissible shapes Praveen. C (TIFR-CAM) CFD Jan 31, / 45

43 Classical method Let the shape S be given by a function F (x). Variation in cost function: F F + δf δj = G(x)δF (x)dx G: Shape gradient Steepest descent: choose δf = λg, λ > 0 δj = λ G 2 (x)dx 0 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

44 Steepest descent method 1 Choose an initial shape S o = S(F o ) 2 Set n = 0 3 Solve the flow equations N(S n, u n ) = 0 4 Compute shape gradient G n 5 Change the shape S n+1 = S(F n + δf n ), δf n = λg n 6 If G n < TOL, then STOP, else n = n + 1, go to step 3 Praveen. C (TIFR-CAM) CFD Jan 31, / 45

45 The End Theory + Experiment + Simulation John D. Anderson: Introduction to Computational Fluid Dynamics Thank You You can download these slides from Praveen. C (TIFR-CAM) CFD Jan 31, / 45

Finite Volume Method

Finite Volume Method Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline

More information

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

ENGR Heat Transfer II

ENGR Heat Transfer II ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value

More information

2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations

2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations . Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations We need to review the governing equations of fluid mechanics before examining the methods of computational aerodynamics in detail.

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis

More information

1 PART1: Bratu problem

1 PART1: Bratu problem D9: Advanced Numerical Analysis: 3 Computer Assignment 3: Recursive Projection Method(RPM) Joakim Möller PART: Bratu problem. Background The Bratu problem in D is given by xx u + yy u + λe u =, u Γ = ()

More information

Module 2: Governing Equations and Hypersonic Relations

Module 2: Governing Equations and Hypersonic Relations Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume

More information

Governing Equations of Fluid Dynamics

Governing Equations of Fluid Dynamics Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present

More information

Chapter 9: Differential Analysis

Chapter 9: Differential Analysis 9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

A Study on Numerical Solution to the Incompressible Navier-Stokes Equation

A Study on Numerical Solution to the Incompressible Navier-Stokes Equation A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field

More information

Convection Heat Transfer

Convection Heat Transfer Convection Heat Transfer Department of Chemical Eng., Isfahan University of Technology, Isfahan, Iran Seyed Gholamreza Etemad Winter 2013 Heat convection: Introduction Difference between the temperature

More information

Chapter 9: Differential Analysis of Fluid Flow

Chapter 9: Differential Analysis of Fluid Flow of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known

More information

Airfoil shape optimization using adjoint method and automatic differentiation

Airfoil shape optimization using adjoint method and automatic differentiation Airfoil shape optimization using adjoint method and automatic differentiation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 5665

More information

Computation of NACA0012 Airfoil Transonic Buffet Phenomenon with Unsteady Navier-Stokes Equations

Computation of NACA0012 Airfoil Transonic Buffet Phenomenon with Unsteady Navier-Stokes Equations 5th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 9-2 January 22, Nashville, Tennessee AIAA 22-699 Computation of NACA2 Airfoil Transonic Buffet Phenomenon with

More information

A Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007

A Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007 A Study of Transonic Flow and Airfoils Presented by: Huiliang Lui 3 th April 7 Contents Background Aims Theory Conservation Laws Irrotational Flow Self-Similarity Characteristics Numerical Modeling Conclusion

More information

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1 AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01

More information

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation

More information

Investigation of Co-Flow Jet Airfoil Mixing Mechanism Using Large Eddy Simulation

Investigation of Co-Flow Jet Airfoil Mixing Mechanism Using Large Eddy Simulation 4st AIAA Fluid Dynamics Conference and Exhibit 27-3 June 2, Honolulu, Hawaii AIAA 2-398 Investigation of Co-Flow Jet Airfoil Mixing Mechanism Using Large Eddy Simulation Hong-Sik IM, Ge-Cheng Zha, Bertrand

More information

AE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Computational Fluid Dynamics (CFD) 1...in the phrase computational fluid dynamics the word computational is simply an adjective to fluid dynamics.... -John D. Anderson 2 1 Equations of Fluid

More information

Follow this and additional works at:

Follow this and additional works at: Washington University in St. Louis Washington University Open Scholarship Mechanical Engineering and Materials Science Independent Study Mechanical Engineering & Materials Science 12-19-2016 The Effects

More information

Challenges and Complexity of Aerodynamic Wing Design

Challenges and Complexity of Aerodynamic Wing Design Chapter 1 Challenges and Complexity of Aerodynamic Wing Design Kasidit Leoviriyakit and Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA kasidit@stanford.edu and

More information

PEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru

PEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru Governing Equations of Fluid Flow Session delivered by: M. Sivapragasam 1 Session Objectives -- At the end of this session the delegate would have understood The principle of conservation laws Different

More information

Aerothermodynamics of High Speed Flows

Aerothermodynamics of High Speed Flows Aerothermodynamics of High Speed Flows Lecture 1: Introduction G. Dimitriadis 1 The sound barrier Supersonic aerodynamics and aircraft design go hand in hand Aspects of supersonic flow theory were developed

More information

Adjoint code development and optimization using automatic differentiation (AD)

Adjoint code development and optimization using automatic differentiation (AD) Adjoint code development and optimization using automatic differentiation (AD) Praveen. C Computational and Theoretical Fluid Dynamics Division National Aerospace Laboratories Bangalore - 560 037 CTFD

More information

Aalto University School of Science and Technology CFD-group/ Department of Applied Mechanics. MEMO No CFD/MECHA DATE: March 15, 2012

Aalto University School of Science and Technology CFD-group/ Department of Applied Mechanics. MEMO No CFD/MECHA DATE: March 15, 2012 Aalto University School of Science and Technology CFD-group/ Department of Applied Mechanics MEMO No CFD/MECHA-2-212 DATE: March 15, 212 TITLE The Effect of Free-Stream Turbulence Parameters on the SST

More information

Entropy stable schemes for compressible flows on unstructured meshes

Entropy stable schemes for compressible flows on unstructured meshes Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/

More information

The Shallow Water Equations

The Shallow Water Equations The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow

More information

Notes 4: Differential Form of the Conservation Equations

Notes 4: Differential Form of the Conservation Equations Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.

More information

Is My CFD Mesh Adequate? A Quantitative Answer

Is My CFD Mesh Adequate? A Quantitative Answer Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011

More information

Adjoint approach to optimization

Adjoint approach to optimization Adjoint approach to optimization Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in Health, Safety

More information

Chapter 2: Fluid Dynamics Review

Chapter 2: Fluid Dynamics Review 7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Chapter 2: Basic Governing Equations

Chapter 2: Basic Governing Equations -1 Reynolds Transport Theorem (RTT) - Continuity Equation -3 The Linear Momentum Equation -4 The First Law of Thermodynamics -5 General Equation in Conservative Form -6 General Equation in Non-Conservative

More information

FLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes

FLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam

More information

Application of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations. Xiaoyue Shen

Application of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations. Xiaoyue Shen Application of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations Xiaoyue Shen Supervisor: Prof. David W. Zingg A thesis submitted in conformity with the requirements

More information

Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations

Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University

More information

An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond

An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics and Astronautics AIAA Aerospace Sciences Meeting, Reno, NV AIAA Paper 2007-0037

More information

High-Order Hyperbolic Navier-Stokes Reconstructed Discontinuous Galerkin Method

High-Order Hyperbolic Navier-Stokes Reconstructed Discontinuous Galerkin Method AIAA SciTech Forum 7- January 29 San Diego California AIAA Scitech 29 Forum.254/6.29-5 High-Order Hyperbolic Navier-Stokes Reconstructed Discontinuous Galerkin Method Lingquan Li Jialin Lou and Hong Luo

More information

Introduction to Physical Acoustics

Introduction to Physical Acoustics Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with

More information

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous

More information

Study on Numerical Simulation Method of Gust Response in Time Domain Jun-Li WANG

Study on Numerical Simulation Method of Gust Response in Time Domain Jun-Li WANG International Conference on Mechanics and Civil Engineering (ICMCE 4) Study on Numerical Simulation Method of Gust Response in Time Domain Jun-Li WANG School of Mechanical Engineering, Shaanxi University

More information

RANS Equations in Curvilinear Coordinates

RANS Equations in Curvilinear Coordinates Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the

More information

Section 26.1: Reporting Conventions. Section 26.2: Fluxes Through Boundaries. Section 26.3: Forces on Boundaries

Section 26.1: Reporting Conventions. Section 26.2: Fluxes Through Boundaries. Section 26.3: Forces on Boundaries Chapter 26. Alphanumeric Reporting FLUENT provides tools for computing and reporting integral quantities at surfaces and boundaries. These tools enable you to find the mass flow rate and heat transfer

More information

for what specific application did Henri Pitot develop the Pitot tube? what was the name of NACA s (now NASA) first research laboratory?

for what specific application did Henri Pitot develop the Pitot tube? what was the name of NACA s (now NASA) first research laboratory? 1. 5% short answers for what specific application did Henri Pitot develop the Pitot tube? what was the name of NACA s (now NASA) first research laboratory? in what country (per Anderson) was the first

More information

Cost target result Temperature

Cost target result Temperature Adjoint Techniques for Aerodynamic Shape Optimization Problems: New Mathematical Formulations and Software Assessment for Well Posed Problems with Acceptable Solutions I. Kampolis, Student D. Papadimitriou,

More information

CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki

CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing Lakhdar Remaki Outline What we doing in CFD? CFD in Industry Shock-filter model for mesh adaptation

More information

Introduction to Fluid Mechanics

Introduction to Fluid Mechanics Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the

More information

APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS

APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK

More information

Aerodynamic study of a small horizontal-axis wind turbine

Aerodynamic study of a small horizontal-axis wind turbine Aerodynamic study of a small horizontal-axis wind turbine Mihai Leonida NICULESCU*,1, Marius Gabriel COJOCARU 1, Cornelia NITA 1 *Corresponding author *,1 INCAS - National Institute for Aerospace Research

More information

University of Bristol - Explore Bristol Research. Peer reviewed version. Link to published version (if available): /6.

University of Bristol - Explore Bristol Research. Peer reviewed version. Link to published version (if available): /6. Huntley, S., Jones, D., & Gaitonde, A. (217). Aeroelastic Gust Response of an Aircraft Using a Prescribed Velocity Method in Viscous Flows. In 23rd AIAA Computational Fluid Dynamics Conference. [AIAA 217-3616]

More information

In this section, mathematical description of the motion of fluid elements moving in a flow field is

In this section, mathematical description of the motion of fluid elements moving in a flow field is Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small

More information

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»

More information

Chapter 4: Fluid Kinematics

Chapter 4: Fluid Kinematics Overview Fluid kinematics deals with the motion of fluids without considering the forces and moments which create the motion. Items discussed in this Chapter. Material derivative and its relationship to

More information

Investigations of Double Surface Co-Flow Jet Transonic Airfoil

Investigations of Double Surface Co-Flow Jet Transonic Airfoil University of Miami Scholarly Repository Open Access Theses Electronic Theses and Dissertations 2019-02-09 Investigations of Double Surface Co-Flow Jet Transonic Airfoil Luchen Wang University of Miami,

More information

EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009)

EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009) EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009) Dr. Mohamad Hekarl Uzir-chhekarl@eng.usm.my School of Chemical Engineering Engineering Campus, Universiti

More information

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit

More information

COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL

COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL Kengo Asada 1 and Kozo Fujii 2 ABSTRACT The effects of body force distribution on the

More information

FLUID MECHANICS. ! Atmosphere, Ocean. ! Aerodynamics. ! Energy conversion. ! Transport of heat/other. ! Numerous industrial processes

FLUID MECHANICS. ! Atmosphere, Ocean. ! Aerodynamics. ! Energy conversion. ! Transport of heat/other. ! Numerous industrial processes SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.)! INL 1 (3 cr.)! 3 sets of home work problems (for 10 p. on written exam)! 1 laboration! TEN1 (4.5 cr.)! 1 written

More information

Lecture: Wave-induced Momentum Fluxes: Radiation Stresses

Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Chapter 4 Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Here we derive the wave-induced depth-integrated momentum fluxes, otherwise known as the radiation stress tensor S. These are the 2nd-order

More information

Performance. 5. More Aerodynamic Considerations

Performance. 5. More Aerodynamic Considerations Performance 5. More Aerodynamic Considerations There is an alternative way of looking at aerodynamic flow problems that is useful for understanding certain phenomena. Rather than tracking a particle of

More information

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil

More information

Introduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)

Introduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD) Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The

More information

Viscous Fluids. Amanda Meier. December 14th, 2011

Viscous Fluids. Amanda Meier. December 14th, 2011 Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport

More information

Chapter 6: Momentum Analysis

Chapter 6: Momentum Analysis 6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of

More information

The Phenomena of Fluid Flow

The Phenomena of Fluid Flow The Phenomena of Fluid Flow Nicholas S. Vlachos Lab. Fluid Mechanics & Turbomachines Department of Mechanical Engineering University of Thessaly Program of Graduate Studies Academic Year 2009-2010 2010

More information

Chapter 9 Flow over Immersed Bodies

Chapter 9 Flow over Immersed Bodies 57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,

More information

Basic Aspects of Discretization

Basic Aspects of Discretization Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)

More information

Module 2 : Lecture 1 GOVERNING EQUATIONS OF FLUID MOTION (Fundamental Aspects)

Module 2 : Lecture 1 GOVERNING EQUATIONS OF FLUID MOTION (Fundamental Aspects) Module : Lecture 1 GOVERNING EQUATIONS OF FLUID MOTION (Fundamental Aspects) Descriptions of Fluid Motion A fluid is composed of different particles for which the properties may change with respect to

More information

An Introduction to the Discontinuous Galerkin Method

An Introduction to the Discontinuous Galerkin Method An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar

More information

Feedback Control of Aerodynamic Flows

Feedback Control of Aerodynamic Flows 44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-843 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Feedback Control of Aerodynamic

More information

Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme

Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of

More information

Computational Astrophysics

Computational Astrophysics Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture

More information

Cavitation modeling using compressible Navier Stokes and Korteweg equations

Cavitation modeling using compressible Navier Stokes and Korteweg equations Computational Methods in Multiphase Flow VIII 425 Cavitation modeling using compressible Navier Stokes and Korteweg equations T. Ait Ali 1, S. Khelladi 1, L. Ramirez 2 & X. Nogueira 2 1 DynFluid Laboratory,

More information

A Multi-Dimensional Limiter for Hybrid Grid

A Multi-Dimensional Limiter for Hybrid Grid APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy

More information

On limiting for higher order discontinuous Galerkin method for 2D Euler equations

On limiting for higher order discontinuous Galerkin method for 2D Euler equations On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation

More information

Reduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis

Reduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis Journal of Mechanical Science and Technology 24 (2) (2010) 601~608 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-010-0106-0 Reduced order model of three-dimensional Euler equations using proper

More information

Flux - definition: (same format for all types of transport, momentum, energy, mass)

Flux - definition: (same format for all types of transport, momentum, energy, mass) Fundamentals of Transport Flu - definition: (same format for all types of transport, momentum, energy, mass) flu in a given direction Quantity of property being transferred ( time)( area) More can be transported

More information

Continuity Equation for Compressible Flow

Continuity Equation for Compressible Flow Continuity Equation for Compressible Flow Velocity potential irrotational steady compressible Momentum (Euler) Equation for Compressible Flow Euler's equation isentropic velocity potential equation for

More information

Validation of the Chemistry Module for the Euler Solver in Unified Flow Solver

Validation of the Chemistry Module for the Euler Solver in Unified Flow Solver Air Force Institute of Technology AFIT Scholar Theses and Dissertations 3-22-2012 Validation of the Chemistry Module for the Euler Solver in Unified Flow Solver William C. Humphrey Jr. Follow this and

More information

Compressible Potential Flow: The Full Potential Equation. Copyright 2009 Narayanan Komerath

Compressible Potential Flow: The Full Potential Equation. Copyright 2009 Narayanan Komerath Compressible Potential Flow: The Full Potential Equation 1 Introduction Recall that for incompressible flow conditions, velocity is not large enough to cause density changes, so density is known. Thus

More information

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16. CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo

More information

PDE Solvers for Fluid Flow

PDE Solvers for Fluid Flow PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic

More information

Overview of Convection Heat Transfer

Overview of Convection Heat Transfer Overview of Convection Heat Transfer Lecture Notes for ME 448/548 Gerald Recktenwald March 16, 2006 1 Introduction CFD is a powerful tool for engineering design involving fluid flow and heat transfer.

More information

Y. Abe, N. Iizuka, T. Nonomura, K. Fujii Corresponding author:

Y. Abe, N. Iizuka, T. Nonomura, K. Fujii Corresponding author: Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2012-2801 ICCFD7-2801 Symmetric-conservative metric evaluations for higher-order finite

More information

n i,j+1/2 q i,j * qi+1,j * S i+1/2,j

n i,j+1/2 q i,j * qi+1,j * S i+1/2,j Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer

More information

Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics

Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/

More information

Introduction to numerical simulation of fluid flows

Introduction to numerical simulation of fluid flows Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural

More information

Dynamics of Glaciers

Dynamics of Glaciers Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers

More information

5. FVM discretization and Solution Procedure

5. FVM discretization and Solution Procedure 5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized

More information

A parallel Newton-Krylov-Schur flow solver for the Reynolds-Averaged Navier-Stokes equations

A parallel Newton-Krylov-Schur flow solver for the Reynolds-Averaged Navier-Stokes equations A parallel Newton-Krylov-Schur flow solver for the Reynolds-Averaged Navier-Stokes equations Michal Osusky and David W. Zingg Institute for Aerospace Studies, University of Toronto, Toronto, Ontario, M3H

More information

LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort

LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort - 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [

More information

CHAPTER 2 INVISCID FLOW

CHAPTER 2 INVISCID FLOW CHAPTER 2 INVISCID FLOW Changes due to motion through a field; Newton s second law (f = ma) applied to a fluid: Euler s equation; Euler s equation integrated along a streamline: Bernoulli s equation; Bernoulli

More information

Differential relations for fluid flow

Differential relations for fluid flow Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow

More information

Mathematical Theory of Non-Newtonian Fluid

Mathematical Theory of Non-Newtonian Fluid Mathematical Theory of Non-Newtonian Fluid 1. Derivation of the Incompressible Fluid Dynamics 2. Existence of Non-Newtonian Flow and its Dynamics 3. Existence in the Domain with Boundary Hyeong Ohk Bae

More information

Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2

Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2 Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2 1 M. Tech Scholar, 2 Associate Professor Department of Mechanical Engineering, Bipin Tripathi

More information

Applied Computational Fluid Dynamics. in Marine Engineering

Applied Computational Fluid Dynamics. in Marine Engineering Applied Computational Fluid Dynamics in Marine Engineering Objectives Understand basic CFD theory Learn how to set up and run simulations in Star CCM+ and interpret results Learn about limitations and

More information

Finite volume method on unstructured grids

Finite volume method on unstructured grids Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen

More information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information