A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
|
|
- Shanna Wilcox
- 6 years ago
- Views:
Transcription
1 A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC Dipartimento di Ingegneria Meccanica - Università di Brescia Eccomas CFD 2006 Sept. 5 8, Egmond aan Zee, NL
2 Presentation outline Governing equations Compressible RANS + k-ω equations Realizability constraints Numerical solution DG space discretization Time integration Results Stanitz elbow ONERA M6 wing Conclusions
3 Compressible RANS + k-ω equations Compressible RANS + k-ω equations RANS + k-ω equations ( ω = log ω) ρ t + (ρu j ) = 0 x j t (ρu i) + (ρu j u i ) = p + τ ji x j x i x j t (ρe 0) + (ρu j h 0 ) = u i [u i τ ij q j ] τ ij + β ρke eωr x j x j x j t (ρk) + (ρu j k) = [ (µ + σ µ x t ) k ] u i + τ ij j x j t (ρ ω) + (ρu j ω) = x j x j x j [ (µ + σµ t ) ω x j + (µ + σµ t ) ω x k ω x k β ρke eωr x j ] + α k τ u i ij βρe eωr x j
4 Compressible RANS + k-ω equations Compressible RANS + k-ω equations RANS + k-ω equations ( ω = log ω) ρ t + (ρu j ) = 0 x j t (ρu i) + (ρu j u i ) = p + τ ji x j x i x j t (ρe 0) + (ρu j h 0 ) = u i [u i τ ij q j ] τ ij + β ρke eωr x j x j x j t (ρk) + (ρu j k) = [ (µ + σ µ x t ) k ] u i + τ ij j x j t (ρ ω) + (ρu j ω) = x j x j x j [ (µ + σµ t ) ω x j + (µ + σµ t ) ω x k ω x k β ρke eωr x j ] + α k τ u i ij βρe eωr x j
5 Compressible RANS + k-ω equations Compressible RANS + k-ω equations Pressure and heat flux vector p = (γ 1)ρ (e 0 u k u k /2), ( µ q j = Pr + µ ) t h Pr t x j Mean strain-rate, molecular, and turbulent stress tensors S ij = 1 ( ui + u ) j, Z ij = S ij 1 u k δ ij 2 x j x i 3 x k τ ij = 2µ t Z ij 2 3 ρkδ ij, τ ij = 2µZ ij + τ ij Turbulence kinetic energy and viscosity coefficient k = max (0, k), µ t = α ρke eωr
6 Realizability constraints Realizability constraints The variable ω = log ω has been used since the near wall distribution of ω is easier to compute than that of ω The limited value k is introduced to deal with possible negative values of k (the solution of the k equation could take negative values) The variable ω r = max( ω r0, ω) is introduced to preclude too small values of ω ω r0, the minimum allowed value of ω r, is determined by enforcing the following realizability conditions Positive normal turbulent stresses ρu 2 i 0 Schwartz inequality for the shear turbulent stresses 2 ρu i u 2 j ρu i ρu 2 j
7 Realizability constraints Realizability for the k-ω model Realizability constraints in terms of the modeled turbulent stresses ( 2µ t S ij ) ρk 2µ tz ii 0, i = 1, 2, 3 ( 2 ) ( 2 ) 3 ρk 2µ tz ii 3 ρk 2µ tz jj i, j = 1, 2, 3, i j Dividing by µ t and recalling that ρk/µ t = e eωr /α (because of the definition µ t = α ρke eωr ) we obtain ( e eω r ) 2 3 (Z ii + Z e eωr jj) α e eωr α 3Z ii 0, i = 1, 2, 3 α +9 ( Z ii Z jj S 2 ij ) 0, i, j = 1, 2, 3, i j
8 Realizability constraints ω r0 High-Reynolds k-ω model Let s a denote the minimum value of e eωr /α allowed by the preceding inequalities, i.e. the maximum root of ( e eω r ) 2 3 (Z ii + Z e eωr jj) α e eωr α 3Z ii = 0, i = 1, 2, 3 α +9 ( Z ii Z jj S 2 ij ) = 0 i, j = 1, 2, 3, i j Standard k-ωmodel: The minimum value ω r0 is trivially determined since α is a constant parameter, i.e. ω r0 = log (α a)
9 Realizability constraints ω r0 Low-Reynolds k-ω model Low-Reynolds k-ωmodel: in this case α = α t α 0 + Re t/r k 1 + Re t /R k where α t, α 0 and R k are constants parameters but Re t = k/(e eωr ν) The minimum value ω r0 is obtained from the solution of the quadratic equation [ e 2eωr0 αt α0a k ] e eωr0 αt R k ν k R k ν a = 0 where a is the maximum root as computed in the High-Reynolds model case
10 DG space discretization Numerical approximation Differential form u t + F c(u) + F v (u, u) + s(u, u) = 0 Weak form φ u Ω t dx φ F(u, u) dx + φf(u, u) n dσ Ω Ω + φs(u, u) dx = 0 Triangulation T h = {K} of an approximation Ω h of Ω Components of the approximate solution u h on T h u hi Φ h def = { φ h L 2 (Ω h ) : φ h K P k (K) K T h } Ω
11 DG space discretization Notation We introduce the symbols f and F to denote a generic face of the triangulation and the set of all the faces, respectively Let f K + K for K +, K T h. We define the following trace operators {q} def = q+ + q f 2 P [[q]] def = q + n + + q n n K + n + K These definitions can be extended to boundary faces by replacing the value on the exterior element by the boundary data
12 DG space discretization DG space discretization Discrete counterpart of weak form for an element K T h K u h φ h dx h φ h F(u h, h u h ) dx t K + φ h F(u h K, h u h K ) n dσ + φ h s(u h, u h ) dx = 0 K K
13 DG space discretization DG space discretization Discrete counterpart of weak form for an element K T h K u h φ h dx h φ h F(u h, h u h ) dx t K + φ h F(u h K, h u h K ) n dσ + φ h s(u h, u h ) dx = 0 K Flux function of boundary integral not uniquely defined K
14 DG space discretization DG space discretization Discrete counterpart of weak form for an element K T h K u h φ h dx h φ h F(u h, h u h ) dx t K + φ h F(u h K, h u h K ) n dσ + φ h s(u h, u h ) dx = 0 K Flux function of boundary integral not uniquely defined A consistent and stable discretization of viscous fluxes must account for the effect of interface jumps on h u h K
15 DG space discretization DG space discretization Discrete counterpart of weak form for an element K T h K u h φ h dx h φ h F(u h, h u h ) dx t K + φ h F(u h K, h u h K ) n dσ + φ h s(u h, u h ) dx = 0 K Flux function of boundary integral not uniquely defined A consistent and stable discretization of viscous fluxes must account for the effect of interface jumps on h u h Thus we introduce suitable numerical fluxes and modified gradients K
16 DG space discretization DG space discretization Sum above equation over the elements Using the trace operators the DG dicretization reads: find u h1,..., u hm Φ h such that u h φ h dx h φ h F (u h, h u h + r([[u h ]])) dx K T K t h K T K h + [[φ h ]] F ( u ± h, ( hu h + η f r f ([[u h ]])) ±) dσ f F f + φ h s (u h, h u h + r([[u h ]])) dx = 0 K T h K for all φ h Φ h
17 DG space discretization DG space discretization The lifting operator r f, applied to the jumps of u h componentwise, is defined as the solution of the problem φ h r f (v) dx = {φ h } v dσ, φ h [Φ h ] N, v [ L 2 (f ) ] N Ω h f (1) The function r is related to r f by r(v) def = r f (v) (2) f F Inviscid flux: Godunov or van Leer/Hänel flux vector splitting Viscous flux ( F u ± h, ( hu h + η f r f ([[u h ]])) ±) def = {F (u h, h u h + η f r f ([[u h ]]))} (3)
18 Time integration Time Integration The DG space discretized equations can be written as M du dt + R (U) = 0 where U is the global vector of DOFs and M is the global block diagonal mass matrix
19 Time integration Time Integration The DG space discretized equations can be written as M du dt + R (U) = 0 where U is the global vector of DOFs and M is the global block diagonal mass matrix Implicit time discretization
20 Time integration Time Integration The DG space discretized equations can be written as M du dt + R (U) = 0 where U is the global vector of DOFs and M is the global block diagonal mass matrix Implicit time discretization Backward Euler method for steady problems Second-order Runge-Kutta scheme for unsteady problems
21 Time integration Backward Euler scheme By linearizing at time level n we obtain [ M t + R ] (Un ) ( U n+1 U n) = R (U n ) U }{{} A where R(Un ) U is the Jacobian matrix of the DG space discretization and A denotes the linear system matrix
22 Time integration Backward Euler scheme By linearizing at time level n we obtain [ M t + R ] (Un ) ( U n+1 U n) = R (U n ) U }{{} A where R(Un ) U is the Jacobian matrix of the DG space discretization and A denotes the linear system matrix The compact DG space discretization here employed minimizes the number of non-zero blocks for each (block) row of A
23 Time integration Backward Euler scheme By linearizing at time level n we obtain [ M t + R ] (Un ) ( U n+1 U n) = R (U n ) U }{{} A where R(Un ) U is the Jacobian matrix of the DG space discretization and A denotes the linear system matrix The compact DG space discretization here employed minimizes the number of non-zero blocks for each (block) row of A To solve the linear system we can use either direct or iterative solvers
24 Time integration Backward Euler scheme By linearizing at time level n we obtain [ M t + R ] (Un ) ( U n+1 U n) = R (U n ) U }{{} A where R(Un ) U is the Jacobian matrix of the DG space discretization and A denotes the linear system matrix The compact DG space discretization here employed minimizes the number of non-zero blocks for each (block) row of A To solve the linear system we can use either direct or iterative solvers Linear algebra and parallelization are handled through PETSc
25 Time integration Implicit Runge-Kutta scheme Second-order scheme (Iannelli and Baker, 1988) [ ] M t + α R (Un ) K 1 = R (U n ) U [ ] M t + α R (Un ) K 2 = R (U n + βk 1 ) U U n+1 U n = Y 1 K 1 + Y 2 K 2, where the constants α, β, Y 1, Y 2 are given by α = 2 ( ) 2 1, β = 8α 2 2 α, Y 1 = 1 1 8α, Y 2 = 1 Y 1
26 Time integration Implicit Runge-Kutta scheme Second-order scheme (Iannelli and Baker, 1988) [ ] M t + α R (Un ) K 1 = R (U n ) U [ ] M t + α R (Un ) K 2 = R (U n + βk 1 ) U U n+1 U n = Y 1 K 1 + Y 2 K 2, where the constants α, β, Y 1, Y 2 are given by α = 2 ( ) 2 1, β = 8α 2 2 α, Y 1 = 1 1 8α, Y 2 = 1 Y 1 For Y 1 = 1, Y 2 = 0, α = 1 we obtain the backward Euler scheme For Y 1 = 1, Y 2 = 0, α = 1/2 we obtain the Crank-Nicolson scheme
27 Stanitz elbow Stanitz elbow Rectangular cross section duct for secondary flows studies Mis = 0.26 Rew,is = 2,738,800 based on channel width w Tu = 0.01, µt/µ = 0.1 Computational grid hexahedral elements hmin /w = Maximum aspect ratio is 2,332 Computational details Restarted GMRES(30) linear solver with 60 iterations DG-P0: ILU prec., 200 time steps, CFL number up to , residuals to machine zero DG-P1: ILU prec., 200 time steps, CFL number up to , residuals to DG-P2: block diag. prec., 200 time steps, CFL number limited to , residuals to
28 Stanitz elbow Stanitz elbow (a) Symmetry plane Mach contours, DG-P 2 solution. (b) Total pressure loss contours, DG-P 2 solution. Figure: Geometry and surface grid.
29 Stanitz elbow Stanitz elbow Figure: Pressure coefficient on suction and pressure walls for two spanwise sections.
30 Stanitz elbow Stanitz elbow Figure: Total pressure loss contours at exit, exp. and DG-P 1 solution.
31 Stanitz elbow Stanitz elbow Figure: Total pressure loss contours at exit, exp. and DG-P 2 solution.
32 Stanitz elbow Stanitz elbow Figure: Detail of total pressure loss contours at exit, exp. and DG-P 2 solution.
33 ONERA M6 wing ONERA M6 wing Validation case for transonic external flow around a wing M = α = 3.06 Rec, = based on the chord c = m Tu = 0.001, µt/µ = 0.1 Computational grid NPARC alliance grid containing 294,912 hexahedral elements hmin /c = Computational details Restarted GMRES(30) linear solver with 60 iterations for P0 and P 1 solutions, Richardson method with 15 iterations for P 2 solution DG-P0: ILU prec., 300 time steps, CFL number up to , residuals to DG-P1: ILU prec., 400 time steps, CFL number up to , residuals to DG-P2: block diag. prec., 800 time steps, CFL number limited to 15, residuals to
34 ONERA M6 wing ONERA M6 wing (a) Surface grid. (b) Pressure coefficient contours, DG-P 2 solution. Figure: ONERA M6 wing.
35 ONERA M6 wing ONERA M6 wing Figure: Pressure coefficient on the wing surface at different spanwise sections.
36 ONERA M6 wing ONERA M6 wing Figure: Pressure coefficient on the wing surface at different spanwise sections.
37 ONERA M6 wing ONERA M6 wing Figure: Pressure coefficient on the wing surface at different spanwise sections.
38 ONERA M6 wing ONERA M6 wing Figure: Pressure coefficient on the wing surface at different spanwise sections.
39 ONERA M6 wing ONERA M6 wing (a) (b) Spanwise sections. Figure: Turbulence intensity contours, DG-P 2 solution.
40 ONERA M6 wing ONERA M6 wing C L C D C Dp C Df DG-P DG-P DG-P Table: Lift and drag coefficients of ONERA M6 wing.
41 Conclusions The paper presented a parallel, high-order implicit DG solver for the RANS and k-ω equations The code proved reliable, robust and accurate for the complex 3D high-reynolds number flows here considered On-going work Improvement of computational efficiency of the implicit DG method Reliable and accurate shock capturing techniques for high-reynolds number flows Assessment of parallel efficiency of the code... further work promised within the ADIGMA project
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationCode MIGALE state- of- the- art
Code MIGALE state- of- the- art A. Colombo HiOCFD4 4th International Workshop on High- Order CFD Method Foundation for Research and Technology Hellas (FORTH), Heraklion (Crete) 4th June 2016 1 with the
More informationAn 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 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 informationCurved grid generation and DG computation for the DLR-F11 high lift configuration
ECCOMAS Congress 2016, Crete Island, Greece, 5-10 June 2016 Curved grid generation and DG computation for the DLR-F11 high lift configuration Ralf Hartmann 1, Harlan McMorris 2, Tobias Leicht 1 1 DLR (German
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
16th International Conference on Finite Elements in Flow Problems DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology
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 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 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 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 informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationInterior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations
Interior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations Will Pazner Brown University, 8 George St., Providence, RI, 9, U.S.A. Per-Olof Persson University of
More informationIs 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 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 informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
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 informationHYBRID RANS-LES SIMULATIONS WITH THE DISCONTINUOUS GALERKIN METHOD
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 5 June 8, Glasgow, UK HYBRID RANS-LES SIMULATIONS WITH THE DISCONTINUOUS GALERKIN
More informationA P-ADAPTIVE IMPLICIT DISCONTINUOUS GALERKIN METHOD FOR THE UNDER-RESOLVED SIMULATION OF COMPRESSIBLE TURBULENT FLOWS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11-15 June 2018, Glasgow, UK A P-ADAPTIVE IMPLICIT DISCONTINUOUS GALERKIN METHOD
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 informationA Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization. Laslo Tibor Diosady
A Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization by Laslo Tibor Diosady B.A.Sc., University of Toronto (2005) Submitted to
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 informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
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 informationCFD 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 informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
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 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 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 informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
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 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 informationAdaptive backward difference formula discontinuous Galerkin finite element method for the solution of conservation laws
Adaptive backward difference formula discontinuous Galerkin finite element method for the solution of conservation laws V. Dolejší and P. Kůs Charles University Prague, Faculty of Mathematics and Physics,
More informationLARGE EDDY SIMULATION OF FLOW OVER NOZZLE GUIDE VANE OF A TRANSONIC HIGH PRESSURE TURBINE
20 th Annual CFD Symposium, August 09-10, 2018, Bangalore LARGE EDDY SIMULATION OF FLOW OVER NOZZLE GUIDE VANE OF A TRANSONIC HIGH PRESSURE TURBINE Bharathan R D, Manigandan P, Vishal Tandon, Sharad Kapil,
More informationAalto 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 informationThe hybridized DG methods for WS1, WS2, and CS2 test cases
The hybridized DG methods for WS1, WS2, and CS2 test cases P. Fernandez, N.C. Nguyen and J. Peraire Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics, MIT 5th High-Order
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationFinite 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 informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationA Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM
206 9th International Conference on Developments in esystems Engineering A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM Hayder Al-Jelawy, Stefan Kaczmarczyk
More informationMODELLING OF INFLUENCE OF TURBULENT TRANSITION ON HEAT TRANSFER CONDITIONS KRZYSZTOF BOCHON, WŁODZIMIERZ WRÓBLEWSKI
TASK QUARTERLY 12 No 3, 173 184 MODELLING OF INFLUENCE OF TURBULENT TRANSITION ON HEAT TRANSFER CONDITIONS KRZYSZTOF BOCHON, WŁODZIMIERZ WRÓBLEWSKI AND SŁAWOMIR DYKAS Institute of Power Engineering and
More informationNumerical Simulation of Newtonian and Non-Newtonian Flows in Bypass
Numerical Simulation of Newtonian and Non-Newtonian Flows in Bypass Vladimír Prokop, Karel Kozel Czech Technical University Faculty of Mechanical Engineering Department of Technical Mathematics Vladimír
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationAerodynamics of Centrifugal Turbine Cascades
ASME ORC 2013 2nd International Seminar on ORC Power Systems October 7th-8th, Rotterdam, The Netherlands Aerodynamics of Centrifugal Turbine Cascades G. Persico, M. Pini, V. Dossena, and P. Gaetani Laboratorio
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
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 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 informationV (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 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 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 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 informationCompressible Flow LES Using OpenFOAM
Compressible Flow LES Using OpenFOAM I. B. Popov supervisor: S. J. Hulshoff promoter: H. Bijl PhD student, TU Delft, Delft, The Netherlands researcher, NEQLab Research B.V., The Hague, The Netherlands
More informationCode Verification of Multiphase Flow with MFIX
Code Verification of Multiphase Flow with MFIX A N I R U D D H A C H O U D H A R Y ( A N I R U D D @ V T. E D U ), C H R I S T O P H E R J. ROY ( C J R O Y @ V T. E D U ), ( V I R G I N I A T E C H ) 23
More informationNear-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Near-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods Axel Probst (now at: C 2 A 2 S 2 E, DLR Göttingen) René Cécora,
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationPhysics of turbulent flow
ECL-MOD 3A & MSc. Physics of turbulent flow Christophe Bailly Université de Lyon, Ecole Centrale de Lyon & LMFA - UMR CNRS 5509 http://acoustique.ec-lyon.fr Outline of the course A short introduction to
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationFourth-Order Implicit Runge-Kutta Time Marching Using A Newton-Krylov Algorithm. Sammy Isono
Fourth-Order Implicit Runge-Kutta Time Marching Using A Newton-Krylov Algorithm by Sammy Isono A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationExtension to moving grids
Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical
More informationDerivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint
Derivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint Charles A. Mader University of Toronto Institute for Aerospace Studies Toronto, Ontario, Canada Joaquim
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 informationHigh Order Discontinuous Galerkin Methods for Aerodynamics
High Order Discontinuous Galerkin Methods for Aerodynamics Per-Olof Persson Massachusetts Institute of Technology Collaborators: J. Peraire, J. Bonet, N. C. Nguyen, A. Mohnot Symposium on Recent Developments
More informationLocal Preconditioning for Low Mach-number Aerodynamic Flow Simulation with a Steady Compressible Navier-Stokes Flow Solver
Local Preconditioning for Low Mach-number Aerodynamic Flow Simulation with a Steady Compressible Navier-Stokes Flow Solver Kurt Sermeus and Eric Laurendeau Advanced Aerodynamics Department Bombardier Aerospace
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 informationNEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS
NEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS P.-O. PERSSON AND J. PERAIRE Abstract. We study preconditioners for the iterative solution of the
More informationSolution Algorithms for Viscous Flow
Solution Algorithms for Viscous Flow Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA Indian Institute of Science Bangalore, India September 20, 2004 Flo3xx Computational
More informationProblem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600
Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 6 Overview This problem is aimed at testing the accuracy and the performance of high-order methods on the direct numerical simulation
More informationCourant Institute of Mathematical Sciences, New York University, New York, NY 10012,
EFFICIENT VARIABLE-COEFFICIENT FINITE-VOLUME STOKES SOLVERS MINGCHAO CAI, ANDY NONAKA, JOHN B. BELL, BOYCE E. GRIFFITH, AND ALEKSANDAR DONEV Abstract. We investigate several robust preconditioners for
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER
THERMAL SCIENCE: Year 2018, Vol. 22, No. 2, pp. 963-972 963 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER by Jitesh RANA, Anshuman SILORI, Rajesh MAITHANI *, and
More informationActive Flux for Advection Diffusion
Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationTURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows
TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows Presented by William A. Sirignano Mechanical and Aerospace Engineering University of California
More informationAPPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM
AIAA 2003-3693 APPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM Vladimir V. Golubev* and Axel Rohde Embry-Riddle Aeronautical University Daytona
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationApplication 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 informationarxiv: v1 [math.na] 23 Jun 2017
On the use of kinetic energy preserving DG-schemes for large eddy simulation David Flad 1 arxiv:1706.07601v1 [math.na] 23 Jun 2017 Institute for Aerodynamics and Gas Dynamics, University of Stuttgart,
More informationImplicit Large Eddy Simulation of Transitional Flow over a SD7003 Wing Using High-order Spectral Difference Method
40th Fluid Dynamics Conference and Exhibit 28 June - 1 July 2010, Chicago, Illinois AIAA 2010-4442 Implicit Large Eddy Simulation of Transitional Flow over a SD7003 Wing Using High-order Spectral Difference
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationTAU Solver Improvement [Implicit methods]
TAU Solver Improvement [Implicit methods] Richard Dwight Megadesign 23-24 May 2007 Folie 1 > Vortrag > Autor Outline Motivation (convergence acceleration to steady state, fast unsteady) Implicit methods
More informationConstrained pseudo-transient continuation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2015; 102:1683 1703 Published online 17 March 2015 in Wiley Online Library (wileyonlinelibrary.com)..4858 Constrained
More informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationUncertainty quantification for RANS simulation of flow over a wavy wall
Uncertainty quantification for RANS simulation of flow over a wavy wall Catherine Gorlé 1,2,3, Riccardo Rossi 1,4, and Gianluca Iaccarino 1 1 Center for Turbulence Research, Stanford University, Stanford,
More informationOptimizing Runge-Kutta smoothers for unsteady flow problems
Optimizing Runge-Kutta smoothers for unsteady flow problems Philipp Birken 1 November 24, 2011 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany. email:
More informationLarge Scale Simulations of Turbulent Flows for Industrial Applications. Lakhdar Remaki BCAM- Basque Centre for Applied Mathematics
Large Scale Simulations o Turbulent Flows or Industrial Alications Lakhdar Remaki BCAM- Basque Centre or Alied Mathematics Outline Flow Motion Simulation Physical Model Finite-Volume Numerical Method CFD
More informationPhysic-based Preconditioning and B-Splines finite elements method for Tokamak MHD
Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD E. Franck 1, M. Gaja 2, M. Mazza 2, A. Ratnani 2, S. Serra Capizzano 3, E. Sonnendrücker 2 ECCOMAS Congress 2016, 5-10
More informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
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 informationDifferential 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 informationRANS 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 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 informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationImplicit kinetic relaxation schemes. Application to the plasma physic
Implicit kinetic relaxation schemes. Application to the plasma physic D. Coulette 5, E. Franck 12, P. Helluy 12, C. Courtes 2, L. Navoret 2, L. Mendoza 2, F. Drui 2 ABPDE II, Lille, August 2018 1 Inria
More information