Direct Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
|
|
- Shonda Hudson
- 5 years ago
- Views:
Transcription
1 19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme Yves Allaneau and Antony Jameson Stanford University, Stanford, California, 94305, USA This paper presents the results of flow computations done using Jameson s finite volume kinetic energy preserving scheme. Direct numerical simulations were performed in a twodimensional viscous shocktube. This testcase proves that the method is stable, extremely low dissipative and computationally fast. Finally unexpected patterns observed in the flow are described. Re Reynolds Number P r Prandtl Number x Spatial Coordinate t Time variable v i i th component of flow velocity ρ Density p Pressure E Total Energy H Specific Total Enthalpy k Specific Kinetic Energy µ Dynamic Viscosity λ Volume Viscosity κ Coefficient of heat conduction Specific Heat C p Nomenclature Ω Ω Fixed spatial domain Boundary of the domain Ω Superscript i, j Coordinate directions Subscript o, p Grid Locations I. Introduction The massive use of Direct Numerical Simulation in today s flow computations is confronted by many obstacles. The complexity of the fluid dynamics equations and their discretized counterparts make the process extremely expensive in terms of computational costs. Additionally, some problems are still encountered when trying to describe complex geometries as automatic mesh generation remains an open question. Graduate Student, Department of Aeronautics and Astronautics, Stanford University, AIAA Member. Professor, Department of Aeronautics and Astronautics, Stanford University, AIAA Member. 1 of 11 Copyright 009 by the American Institute of Aeronautics and American Astronautics, Institute Inc. All of rights Aeronautics reserved. and Astronautics
2 Furthermore, even if one manages to set up and run such a computation, results can be disappointing and unrepresentative of a real phenomenon. One example is when excessive dissipation is introduced by the numerical scheme. Very often, dissipation is added in the computation to help capture shockwaves or enhance stability. If the addition of artificial viscosity is not completely controlled and strictly limited to where it is needed, dampening could occur over the entire computational domain. In that case, some features of the flow, such as eddies in turbulent flow, can be completely erased by the simulation. Jameson recently developed a scheme that overcomes some of these problems. His finite volume Kinetic Energy Preserving scheme respects the global balance of kinetic energy over the computational domain, along with the balance of mass, momentum and total energy. The respect of the kinetic energy law for smooth flow provides enough stability to perform viscous computations of more complex flows, even involving the computation of shockwaves. As this method does not add any artificial viscosity, we can expect to observe fine features in the flow that would not necessarily be observed with the use of other schemes. The purpose of this paper is to prove that the Kinetic Energy Preserving scheme can be efficiently implemented and provides good quality results on direct numerical simulations of viscous flows. Section of this paper briefly describes the Kinetic Energy Preserving scheme developed by Jameson and the formalism used to solve the equations. Section 3 presents a computation of the flow in a shocktube using the Kinetic Energy Preserving scheme and the Roe scheme. A comparison is made based on the results obtained. In Section 4, we study the flow more carefully in the pseudosteady area. We define pseudosteady area as the zone between the expansion and the shockwave containing the contact discontinuity. In the inviscid approximation, the flow is uniform and steady in this area. When viscosity is present, some interesting patterns can be observed. II. Kinetic Energy Preserving scheme for viscous flow In a recent paper, Jameson 1, presented a new semi-discrete finite volume scheme to solve the Navier- Stokes equations. This scheme has the property to satisfy the global conservation law for kinetic energy. We shall briefly describe this scheme in the present section. A. Continuous Model First, consider the three-dimensional Navier-Stokes equations in their conservative form: where ρ ρv 1 u = ρv ρv 3 ρe u t + x i f i u) = 0 1) and f i = ρv i ρv i v 1 + pδ i1 σ i1 ρv i v + pδ i σ i ρv i v 3 + pδ i3 σ i3 ρv i H v j σ ij q j The viscous stress tensor σ ij is given for a Newtonian fluid by σ ij ij vk v = λδ + µ i x k x + vj j x ). Often i in aerodynamics, λ is taken to be equal to 3 µ. The heat flux is proportional to the temperature gradient Fourier s Law) q j = κ T x. j An equation for the kinetic energy k = 1 ρvi momentum equations. Indeed, It follows by substituting t ρv i ) and ρ t k t = ) 1 t ρvi = v i t ) can be derived by combining the continuity and the ρv i ) vi ρ t by their corresponding fluxes that: of 11
3 k t + ) ] [v j x j p + ρ vi v i σ ij = p vj vi σij xj x j 3) We assume that we are interested in a domain Ω fixed in space. Ω denotes the boundary of Ω. By integrating 3) over the domain Ω, we get a global conservation law for kinetic energy. t Ω kdv = Ω [v j p + ρ vi ) v i σ ij ] n j ds + Ω ) p vj vi σij xj x j dv 4) Definition 1 A numerical scheme to solve the viscous Navier-Stokes equations is said to be Kinetic Energy Preserving if it satisfies a discrete analog of 4). Here we have assumed that the domain contains no discontinuity. If a shockwave is present in Ω, the relation 4) does not hold anymore. B. Semi-discrete approach Now, we consider a finite volume discretization of the governing equations in the domain Ω. The generic cell is a polyhedral control volume o. Each cell has one or more neighbors. The face separating cell o and cell p has an area A op, and we define n i op to be the unit normal to this face, directed from o to p. Evidently n i op = n i po. We also define S i op = A op n i op. S i op can be interpreted as the projected face area in the coordinate direction i. Boundary control volumes are closed by an outer face of directed area S i o = p Si op a control volume is delimited by a closed surface). In this framework, the semi-discrete finite volume approximation of the governing equations takes the form: u o vol o t + fop i n i opa op = 0 5) or vol o u o t + p neighbor p neighbor f i op S i op = 0 6) For a boundary control volume b, another contribution to the fluxes f i b S b comes from the outer face. Now we assume that u o and f i op take the form: ρ o ρ o v 1 o u o = ρ o vo ρ o vo 3 ρ o E o and fop i = ρv i ) o ρv i v 1 ) op + pδ i1 σ i1 ) op ρv i v ) op + pδ i σ i ) op ρv i v 3 ) op + pδ i3 σ i3 ) op ρv i H) op v j σ ij + q j ) op 7) Jameson has exhibited a set of sufficient conditions on the elements of f i op that lead to a Kinetic Energy Preserving KEP) scheme. Proposition 1 If the elements of f i op defined in 7) satisfy the following conditions: 3 of 11
4 a - ρv i v j ) op = 1 ρvi ) op v j p + v j o) b - pδ ij σ ij ) op = 1 pδij σ ij ) o + 1 pδij σ ij ) p and if the fluxes at the boundaries are evaluated such that: c - fb i = f i u b ) where b is a boundary control volume then the semi discrete finite volume scheme 6) satifies the discrete global variation law for kinetic energy. Indeed in that case, the discrete kinetic energy k o satisfies the following relation: d dt vol o k o = o b + o S j b p o p v j b vb i p b + ρ b ) v i bσ ij b vo i + vp i Sop i σo ij p ) vo i + vp i Sop i ) 8) which is indeed a discretization of 4). Condition a of the previous proposition is not very restrictive and allows some degrees of freedom in the construction of the fluxes defined in 7). Let s denote by g op the arithmetic average of the quantity g between cell o and cell p : g op = g o + g p )/. We can rewrite condition a ρv i v j) op = ρv i) op vj op 9) We can evaluate the average ρv i) op by any means ρv i) op = ρ opv i op or ρv i op for example) and then deduce ρv i v j) by using 9) to satisfy condition a. op This degree of freedom could be used to design schemes that also satisfy other properties than Kinetic Energy Preservation. III. Direct Numerical Simulation of a two-dimensional flow in a shocktube The scheme developed by Jameson has been applied on two-dimensional computations. We showed experimentally that the conditions stated in section II, which are sufficient to conserve numerical global kinetic energy in the case of a smooth flow are actually providing enough stability to perform computations in a flow containing discontinuities such as shockwaves. This improved stability removes the need of adding artificial viscosity to capture shockwaves, the main drawback of artificial viscosity being that it tends to add dissipation everywhere on the computational domain. The chosen test case is a two-dimensional viscous Sod shocktube 3, 4. It allows us to study the behavior of the scheme in boundary layers, shockwaves and expansions. Although the scheme is designed for general unstructured mesh, we implemented it on rectilinear and curvilinear grids, as the shape of the shocktube is quite simple. The shocktube is characterized by its left and right initial states, both at rest at the beginning of the simulation. The left state is characterized by its pressure p l, its density ρ l and its temperature T l. The right state is also characterized by its pressure p r, its density ρ r and its temperature T r. Pressure, density and temperature can be related by the perfect gas law. The length of the shocktube is L, the height is h. The aspect ratio of the shocktube is defined by α = h/l. We assume that the walls of the shocktube are rigid and adiabatic. Viscosity is evaluated using Sutherland s formula µt ) = C T 3/ T + S For air, at reasonable temperatures, C = kg/ ms K) and S = 11K. 4 of 11
5 We define the velocity V l = p l /ρ l proportional to the speed of sound in the left region and the Reynolds number Re = ρ llv l µ l 10) where µ l = µt l ). The Prandtl number is given by it was taken to be equal to 5. P r = µc p κ 11) Numerical computations were done for the Kinetic Energy Preserving scheme using this averaging formula for the convective terms: ρv i ) op = ρ opv i op ρv i v j) op = ρ opv i opv j op ρv i H ) 1) op = ρ opv i oph op Note that condition a of proposition 1) does not require a specific form for ρv i H ). We just chose it to op be consistent with the continuity and momentum fluxes. Viscous stress was evaluated in each cell by introducing a complementary mesh, for which cell vertices are the centers of the original control volumes. Time integration was performed using the Total Variation Diminishing Runge Kutta 3 like scheme proposed by Shu 5. For a semi discrete scheme in the form this 3-stages scheme advances from time n to time n + 1 by u + Ru) = 0 13) t u 1 = u n tru n ) u = 3 4 un u1 1 4 tru1 ) u n+1 = 1 3 un + 3 u 3 tru ) This time discretization does not guarantee the preservation of kinetic energy in time. One could use a Crank-Nicholson semi implicit scheme as suggested by Jameson 1 to ensure conservation in time, but the computational costs would increase in important proportions. Eventually, results were checked using a classic Roe scheme 6, advanced explicitly in time by a Runge Kutta scheme. A. Simple shocktube and comparison with Roe Scheme First, we ran the simulation for the case ρ r = 5 ρ l 14) p r = p l 15) v l = v r = 0 16) T l = T r = 300 K 17) The Reynolds number was Re = 5000 and the aspect ratio α =. This odd shaped shocktube guarantees that the boundary layers will not affect the flow at the centerline too much. The grid size was 4096 cells in the x-direction and 56 cells in the y-direction only half of the domain is computed using this 5 of 11
6 mesh, the second half is obtained by symmetry). The grid is uniform in the x-direction but stretched in the y-direction such that y min /L = α/4000. Figure 1 shows at time t = 136L/V l the variations of nondimensional pressure, density, velocity and energy at the centerline. Figure presents a comparison of the pressure profile along the centerline for the Roe scheme and the KEP scheme. We can notice that the KEP scheme provides sharper results than the Roe scheme which prooves that it introduces less dissipation. Results are especially convincing at the shockwave. 1 Pressure 1 Density 0 1 a) Pressure 0 1 b) Density 1 x-velocity 3 Energy c) x-velocity d) Energy Figure 1. Variation of state variables along the centerline. Re = 5000, α = Figure 3 is a global picture of the x-velocity in the shocktube. The shape of the boundary layer and the curved aspect of the shockwave near the walls are in agreement with the usual results observed in a viscous shocktube 8, 9. Lighter colors coincide with faster flow. B. Study of nonclassical effects in the pseudosteady flow area In the previous part, α was chosen to be equal to for a Reynolds number of 5000 in order to observe phenomena similar to the inviscid case on the centerline. Actually, if we had plot the pressure on a line closer to the wall, the picture would have been quite different. The quasi complete absence of dissipation introduced by the KEP scheme allowed us to capture some unexpected features of the flow in the pseudosteady region, where the contact discontinuity is see Figure 4). Figure 5 is a sketch of the pattern observed in the flow for the case described above. A +) indicates a pressure wave or pressure point, where the pressure is larger than in the inviscid case. A -) indicates a depression, where the pressure is lower than in the inviscid computation. 6 of 11
7 1 Pressure KEP Pressure Roe Pressure KEP Pressure Roe a) At the start of the expansion b) Through the shockwave Figure. Comparison of pressures on the centerline for the KEP scheme and the Roe scheme at two locations. y/l 0 0 Figure 3. Distribution of nondimensional x -velocity in the shocktube The first obvious pattern that can be observed is the pressure waves developing at the base of the expansion, near the walls. These waves are starting in the boundary layer and are curved towards the direction of the flow. Figure 6 shows the shape of pressure waves in the pseudosteady area for 3 different values of α. When α is varied, the shape of the waves remains the same. As a consequence, when α is reduced too much, the waves end up by crossing each other case α = ). If α is further reduced case α = on the figure), waves will reflect on the walls. On the other side of the pseudosteady area, near the shockwave this time, we can observe depression waves. These are visible on figures 4 and 6 again. They seem to start in the boundary layer near the shockwave and extend upstream in the pseudosteady area. These wave are much more smooth than the one previously described. Their s shape is particularly obvious in the case α =, but the way these waves interfere for smaller values of α is not very clear. Eventually, a third pressure pattern can be observed in the flow, this time not visible on the previous figures. Just after the shockwave after the shockwave is in the pseudosteady region), in the boundary layer, a high pressure point can be observed, as shown on figure 7. Figure 7 represents the variation of pressure along the wall for the case α =, Re = Interesting fact is that this high pressure point is located at the root of the depression wave described above. The small bump observed in figure 7-a for coincides with the reflection of the pressure waves decribed earlier. Note on 7-b the overshoot in pressure located near the shockwave. Once again, it seems that the KEP scheme is less dissipative than the Roe scheme. The interpretation of these various patterns have not yet been established and more simulations are 7 of 11
8 a) x-velocity 1 b) y-velocity c) Pressure Figure 4. Distributions of nondimensional velocities and pressure in the pseudosteady area contact discontinuity area) in the case Re = 5000, α =, t = 136L/V l expansion shockwave + + centerline + + Figure 5. Pressure pattern observed in the pseudosteady area of the flow. A +) is a surpressure compared to the inviscid case while a ) corresponds to a depression. Re = 5000, α =. required to understand how they develop in time. 8 of 11
9 a) α =, y min /L = α/000 b) α =, y min /L = α/400 c) α =, y min /L = α/4000 Figure 6. Pressure Waves pattern in the pseudosteady flow area for α =, α =, α = 06. Re = of 11
10 1 Wall Pressure KEP 5 Wall pressure KEP Wall pressure Roe a) Wall pressure KEP b) Shockwave wall pressure Figure 7. Pressure distribution along the walls of the shocktube. α =, Re = 5000 IV. Conclusion This paper proves that the Kinetic Energy Scheme can be applied efficiently on more complex geometries. It as been seen that it allows a correct description of shockwaves, expansion fans and boundary layers. Experimentally, it has been confirmed that the number of grid points required to ensure numerical stability is proportional to the local cell Reynolds number Re o = ρv 3 vol o µ. For a global Reynolds number of 5000, stability was ensured using grids as large as 4096 mesh cell in the x-direction by 56 mesh cells in the y-direction. This might seem to be a lot, but the KEP fluxes are extremly simple and cheap to evaluate on each cell face. Furthermore, the use of an explicit code makes it easily parallelizable with a large scalability, by simple domain decomposition. For example, when computing the flow in the shocktube for a global Reynolds number of 5000 and α =, using a cells mesh and y min /L = α/000, the computation took only 9 hours on a small 16 CPUs cluster. Acknowledgments Yves Allaneau is supported by a Stanford Graduate Fellowship. This work is also supported by the AFOSR grant #FA from the Computational Math Program under the direction of Dr. Fariba Fahroo. References 1 A. Jameson, The Construction of Discretely Conservative Finite Volume Schemes that Also Globally Converve Energy or Entropy, Journal of Scientific Computing, Vol. 34-, p , 008 A. Jameson, Formulation of Kinetic Energy Preserving Conservative Schemes for Gas Dynamics and Direct Numerical Simulation of One-Dimensional Viscous Compressible Flow in a Shock Tube Using Entropy and Kinetic Energy Preserving Schemes, Journal of Scientific Computing, Vol. 34-, p , G. A. Sod, Survey of Several Finite Difference Methods for Systems of Nonlinear Hyperbolique Conservation Laws, Journal of Computational Physics, Vol. 6, p. 1-31, P. Jorgenson and E. Turkel, Central Difference TVD and TVB Schemes for Time Dependant and Steady State Problems, 30th Aerospace Sciences Meeting and Exhibit sponsored by the AIAA, Reno, Nevada, January C. W. Shu, Total Variation Diminishing Time Discretizations, SIAM Journal on Scientific and Statistical Computing, Vol. 9, p P. L. Roe, Approximate Riemann Solvers, Parameter Vectors and Difference Scheme, Journal of Computational Physics, Vol. 14, p , A. Jameson, Analysis and Design of Numerical Schemes for Gas Dynamics 1 Artificial Diffusion, Upwind Biasing, Limiters and Their Effects on Accuracy and Multigrid Convergence, International Journal of Computational Fluid Dynamics, Vol. 5, p. 1-38, H. Mirels, Test Time in Low-Pressure Shock Tubes, The physics of Fluids, vol. 6-9, p , of 11
11 9 H. Mirels, Flow Nonuniformities in Shock Tubes Operating at Maximimum Test Times, The Physics of Fluids, Vol. 9-10, p , of 11
Direct Numerical Simulations of Plunging Airfoils
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 010, Orlando, Florida AIAA 010-78 Direct Numerical Simulations of Plunging Airfoils Yves Allaneau
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationAntony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
More informationA 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 informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationModeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique
Modeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique Torsten Palenschat 317220 30.04.2014 Master seminar CES WS 2013/2014 Center for Computational Engineering Sciences RWTH Aachen
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More information41st AIAA Aerospace Sciences Meeting and Exhibit January 6 9, 2003/Reno, NV
AIAA 3 267 Calculation of Unsteady Transonic Flow by an Euler Method with Small Perturbation Boundary Conditions C. Gao, S. Luo, F. Liu Department of Mechanical and Aerospace Engineering University of
More informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationFeedback 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 informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationApplication 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 informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationTurbomachinery Applications with the Time Spectral Method
35th AIAA Fluid Dynamics Conference and Exhibit, June 6-9,005, Toronto,Ontario Turbomachinery Applications with the Time Spectral Method Edwin van der Weide, Arathi K Gopinath and Antony Jameson Stanford
More informationAntony Jameson. Stanford University Aerospace Computing Laboratory Report ACL
Formulation of Kinetic Energy Preserving Conservative Schemes for Gas Dynamics and Direct Numerical Simulation of One-dimensional Viscous Comressible Flow in a Shock Tube Using Entroy and Kinetic Energy
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
More informationImprovements 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 informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationIntroduction 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 informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationSub-Cell Shock Capturing for Discontinuous Galerkin Methods
Sub-Cell Shock Capturing for Discontinuous Galerkin Methods Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 39, U.S.A. A shock capturing strategy for higher order
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationNew DNS Results of Shockwave/Turbulent Boundary Layer Interaction
36th AIAA Fluid Dynamics Conference, June 5 8, 26/San Francisco,California New DNS Results of Shockwave/Turbulent Boundary Layer Interaction M. Wu and M.P. Martin Mechanical and Aerospace Engineering Department
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationEffect of Applied Magnetic Field on Shock Boundary Layer Interaction
University of Kentucky UKnowledge Mechanical Engineering Faculty Publications Mechanical Engineering 1-2012 Effect of Applied Magnetic Field on Shock Boundary Layer Interaction Ovais U. Khan University
More informationA KINETIC SCHEME FOR THE NAVIER-STOKES EQUATIONS AND HIGH-ORDER METHODS FOR HYPERBOLIC CONSERVATION LAWS
A KINETIC SCHEME FOR THE NAVIER-STOKES EQUATIONS AND HIGH-ORDER METHODS FOR HYPERBOLIC CONSERVATION LAWS A DISSERTATION SUBMITTED TO THE DEPARTMENT OF AERONAUTICS AND ASTRONAUTICS AND THE COMMITTEE ON
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationOptimum aerodynamic design using the Navier-Stokes equations
Copyright 1997, American Institute of Aeronautics and Astronautics, Inc. AIAA Meeting Papers on Disc, January 1997 A9715189, AF-AFOSR-91-0391, AIAA Paper 97-0101 Optimum aerodynamic design using the Navier-Stokes
More informationLarge Eddy Simulation of Three-Stream Jets
Large Eddy Simulation of Three-Stream Jets J. Xiong 1, F. Liu 2, and D. Papamoschou 3 University of California, Irvine, Irvine, CA, 92697 We present a computational study of three-stream jets simulating
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationNumerical Simulation of Flow Separation Control using Multiple DBD Plasma Actuators
Journal of Applied Fluid Mechanics, Vol. 9, No. 4, pp. 1865-1875, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. DOI: 10.18869/acadpub.jafm.68.235.25325 Numerical Simulation
More informationn 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 informationAssessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow
Assessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering - University
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationA Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems
A Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems Dia Zeidan,a) and Eric Goncalves 2 School of Basic Sciences and Humanities, German Jordanian University, Amman, Jordan
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationURANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver
URANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver *Yohan Choe 1), Hyeongjun Kim 1), Chongam Kim 2) 1), 2) Department of Aerospace Engineering, Seoul National
More informationAnalysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data
45th AIAA Aerospace Science Meeting and Exhibit, January 8 11, 25/Reno, Nevada Analysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data M. Wu and M.P. Martin Mechanical
More informationNonlinear iterative solvers for unsteady Navier-Stokes equations
Proceedings of Symposia in Applied Mathematics Nonlinear iterative solvers for unsteady Navier-Stokes equations Philipp Birken and Antony Jameson This paper is dedicated to Gene Golub. Abstract. The application
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationCOMPUTATIONAL 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 informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationEntropy 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 informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationRANS Solutions Using High Order Discontinuous Galerkin Methods
RANS Solutions Using High Order Discontinuous Galerkin Methods Ngoc Cuong Nguyen, Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 2139, U.S.A. We present a practical
More informationActive Flutter Control using an Adjoint Method
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-844 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Active Flutter Control using an
More informationConservation 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 informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationValidation of an Entropy-Viscosity Model for Large Eddy Simulation
Validation of an Entropy-Viscosity Model for Large Eddy Simulation J.-L. Guermond, A. Larios and T. Thompson 1 Introduction A primary mainstay of difficulty when working with problems of very high Reynolds
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
More informationComputational issues and algorithm assessment for shock/turbulence interaction problems
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln NASA Publications National Aeronautics and Space Administration 2007 Computational issues and algorithm assessment for shock/turbulence
More informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationJun 22-25, 2009/San Antonio, TX
19th AIAA Computational Fluid Dynamics 22-25 June 2009, San Antonio, Texas AIAA 2009-4273 AIAA 2009 4273 An Assessment of Dual-Time Stepping, Time Spectral and Artificial Compressibility based Numerical
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationSistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas
Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Pedro González-Casanova Henríquez Unidad de Investigación en Cómputo Aplicado DGSCA, UNAM Elvira Torondel y Ricard Garrido
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationNonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50
Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations Matthew McMullen Advisor: Antony Jameson Co-advisor: Juan Alonso Sponsor: Accelerated Strategic Computing Initiative
More informationAPPLICATION 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 informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationLES Study of Shock Wave and Turbulent Boundary Layer Interaction
51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 07-10 January 2013, Grapevine (Dallas/Ft. Worth Region), Texas AIAA 2013-0984 LES Study of Shock Wave and
More informationA CFD Approach to Modeling Spacecraft Fuel Slosh
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 2009, Orlando, Florida AIAA 2009-366 A CFD Approach to Modeling Spacecraft Fuel Slosh Introduction:
More informationGT UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE USING AN EFFICIENT NON-LINEAR HARMONIC BALANCE METHOD
Proceedings of ASME Turbo Expo 213: Turbine Technical Conference and Exposition GT213 June 3-7, 213, San Antonio, Texas, USA GT213-94574 UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationNumerical investigation of swirl flow inside a supersonic nozzle
Advances in Fluid Mechanics IX 131 Numerical investigation of swirl flow inside a supersonic nozzle E. Eslamian, H. Shirvani & A. Shirvani Faculty of Science and Technology, Anglia Ruskin University, UK
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
More informationFEDSM COMPUTATIONAL AEROACOUSTIC ANALYSIS OF OVEREXPANDED SUPERSONIC JET IMPINGEMENT ON A FLAT PLATE WITH/WITHOUT HOLE
Proceedings of FEDSM2007: 5 th Joint ASME/JSME Fluids Engineering Conference July 30-August 2, 2007, San Diego, CA, USA FEDSM2007-37563 COMPUTATIONAL AEROACOUSTIC ANALYSIS OF OVEREXPANDED SUPERSONIC JET
More informationMultiscale Hydrodynamic Phenomena
M2, Fluid mechanics 2014/2015 Friday, December 5th, 2014 Multiscale Hydrodynamic Phenomena Part I. : 90 minutes, NO documents 1. Quick Questions In few words : 1.1 What is dominant balance? 1.2 What is
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationNon-Equilibrium Flow Modeling Using High- Order Schemes for the Boltzmann Model Equations
Purdue University Purdue e-pubs School of Aeronautics and Astronautics Faculty Publications School of Aeronautics and Astronautics 28 Non-Equilibrium Flow Modeling Using High- Order Schemes for the Boltzmann
More informationIntroduction of compressible turbulence
Introduction of compressible turbulence 1 Main topics Derive averaged equations for compressible turbulence Introduce a math. technique to perform averaging in presence of density variation Favre average
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More information