CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki
|
|
- Giles Weaver
- 5 years ago
- Views:
Transcription
1 CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing Lakhdar Remaki
2 Outline What we doing in CFD? CFD in Industry Shock-filter model for mesh adaptation Motivation The model Brief insight in the mathematical analysis Numerical tests Conclusions
3 What we doing in CFD? Navier-Stokes Equations t ρ + ( ρv ) = 0 Continuity equation t ( ρv ) + ( ρ v v ) + p = ρf + ( λ v ) + ( 2 µ D ) Momentum equation t ( ρe ) + [ ( ρe + p )v ] + p = ρf v + (κ θ ) + (λ ( v )v ) + (2 µdv ) Energy equation
4 Bloodhound SSC Project Constructing a vehicle to take the World Land Speed Record to 1000 mph
5 Bloodhound SSC Project Thrust SSC: 763mph Bloodhound SSC: 1000mph
6 Governing Equations d α α Qd V = F ( Q ) n d S G ( Q ) n ds α α dt V S S ρ ρu 1 Q = ρu 2 ρ u 3 E 2 u τ βα = µ k δ βα 3 xk ρu α ρu1uα + pδ α 1 Fα = ρu 2 uα + pδ α 2 ρu 3uα + pδ α 3 ( E + p ) u α u β uα + µ + x xβ α Stress tensor α = 1,2,3 0 τ 1α Gα = τ 2α τ 3α u βτ 3α qα qα = k T xβ Average heat flux
7 Equations Discretization Finite Volume Method (FLITE system) Cell vertex finite volume solution (Dual mesh) Time discretization: explicit multi-stage Runge Kutta Convergence acceleration to steady state by local time stepping and an agglomerated multigrid process.
8 Inviscid flux Approximation α F ( Q ) n ds α S ~ F IJ J Λ I ρqij 1 ρu1qij + pnij ~ FIJ = ρu2 qij + pnij2 ρu q + pn31 IJ 3 IJ ( E + p )qij Solve a Riemann problem ~ Q + FIJ (Q) = 0 t s QI if s < 0 Q (0, s ) = QJ if s > 0 qij = uα nαij
9 Numerical Results Transonic viscous flow over an ONERA M6 wing M= , AoA=3.06, Re=11.720e6, tetrahedral, 10 prismatic layers Transonic flow over ONRA M6 wing. Used mesh and Cp profile
10 Numerical Results Comparison to Experimental Results ONERA M6: Cp chordwise profile at 44% of wing semispan ONERA M6: Cp chordwise profile at 65% of wing semispan
11 Numerical Results Comparison to Experimental Results ONERA M6: Cp chordwise profile at 80% of wing semispan ONERA M6: Cp chordwise profile at 90% of wing semispan
12 Numerical Results Supersonic inviscid flow over an F15 fighter M=2.0, AoA=3.0, tetrahedral
13 Numerical Results Residual curve
14 Numerical Results M=3.0, AoA=3.0 Cp Profile Residual curve
15 Spay Drag Model for Bloodhound SSC Supersonic Vehicle: Motivation Pressure Shocks around the vehicle
16 Governing Equations Fluid phase (Navier stokes-equations) tφ f ρ f + (φ f ρ f U f ) = 0 on Ω [ 0.T ] t (φ f ρ f U f ) + (φ f ρ f U f U f ) = φ f Pf +.τ β (U p U f ) + φ f ρ f g Solid phase tφ p + (φ pu p ) = 0 on Ω [ 0.T ] t (φ pu p ) + (φ pu p U p ) = φp +φ f = 1 φp ρf Pf + β (U p U f ) + φ p (1 ) g ρp ρp φpρ f U p U f φ p2 µ f φ f D p2 Dp β = 3 φ p ρ f U p U f 2.65 φp 4 Cd D p R ( Re p ) if Re p < 1000 Cd = e p Drag coefficient 0.4 else Re p = if φ f < 0.8 if φ f 0.8 D pφ U p U f f νf Particle Reynolds number g Gravity
17 Validation Experimental Apparatus: (a) General Flow System, (b) Geometry of the curved square duct Y. Kliafas, M. Holt, LDV measurements of a turbulent airsolid twophase fow in a 90 bend. Experiments in Fluids 5, 7385.
18 Validation Hybrid used mesh and Volume Fraction profile for the curved duct
19 Validation Residual convergence
20 Validation Station θ=0 Station θ=15 Station θ=45 Station θ=30 90 Bend case: Mean Stream fluid and particles velocity comparison to experimental results for different stations
21 Application to Bloodhound SSC Bloodhound SSC Supersonic Car: Hybrid mesh The delimited area using normal velocity gradient criterion
22 Application to Bloodhound SSC Volume fraction profile
23 Application to Bloodhound SSC Residual convergence (a) Drag convergence before and after injecting sand particles: (a) Volume Fraction= 1.5e-3, Drag increases by 5%. (b) Volume Fraction = 5e-3, Drag increases by 10%. (b)
24 Bloodhound SSC Project Aerodynamic Optimization Initial Geometry Final Geometry
25 Tricycle recreational vehicle (CTA-BRPUDES)
26 CAD Mesh
27 CAD Mesh
28 CFD Model CAD/CATIA Meshing/CCM + CAD Fix Flow simulation/ccm+ Suitable mesh parameters Optimizatio n Validation/Experimental 28
29 Velocity profile 29
30 Objective: Minimizing Drag by Fender Shape modification Drag drop 26.15% Drag drop13.18%
31 Shock-filter model for mesh adaptation The CFD solutions are mesh dependant Poor mesh quality leads to a non-accurate CFD solution Generating good mesh quality is a difficult task Solution: Automatic mesh adaptation
32 Example: supersonic flow around a NACA0012
33 Process of an Anisotropic SolutionMesh Adaptation OptiGrid Initial grid Solver OptiGrid Solver OptiGrid Solver Final grid Solver
34 Mesh adaptation overview, 1 Adaptation principle: Objective: Estimation Minimize the computational error u u h u uh X < C1 u rh u X over the domain X Interpolation error For linear interpolation we have: u rh u X < C 2 h 2 D 2u h X Procedure: Build a new mesh from the previous one by Equidistributing the interpolation error bound h 2 D 2u h ( x ) of a selected scalar variable
35 Mesh adaptation overview, 3 To obtain an anisotropic mesh we try to equi-distribute the error along the edges The v second derivative along an edge 1 D u ( x) = 2 < Hv, v > v2 2 v h Derive The M = H = P 1 D P a metric error along and edge v 1 T ε (v ) = v M ( s )v ds 0 H Hessian
36 Mesh adaptation overview, 4 Optimization operations include: nodes movement, refinement by cutting edges with large error and coarsening by concatenating edges with small errors
37 Multiple shocks problem Problem: the aim of mesh adaptation is to start with an arbitrary and coarse initial mesh, taking less time and effort to generate. However, when the initial grid is too coarse discontinuities are badly defined, especially in the case of multiple shocks, whatever the solver used. The proposed solution: Enhance the discontinuities of the selected scalar variable before adaptation. How? Using a shock filter.
38 Shock-filter model, 1 Rudin filter : Nonlinear hyperbolic equation ut + F (u xx ) u x = 0 ( x, t ) (Ω IR) IR + u (0, t ) = u0 ( x) Initial condition sf ( s ) > 0 s Ω Steady state
39 Shock-filter model, 2 The proposed shock-filter ut + a1 ( x, y ) F1 ( u, u x ) x f1 (u ) + a2 ( x, y ) F2 ( u, u y ) y f 2 (u ) = 0 ( S ) (t, x, y ) IR Ω u (0, x, y ) = u0, srf1 ( s, r ) > 0 srf2 ( s, r ) > 0 ( s, r ) Ω F1, F2, f1, f 2 C a1, a2 L (Ω)
40 Mathematical Analysis Existence and uniqueness of a solution of equation (S) in a suitable space Classical frame of Distribution space Obstacle: The product of two distribution is not defined in general (Schwartz theorem) Generalized functions framework (Colombeau algebra): More suitable
41 Insight into the generalized functions algebra, 1 Intuitive idea Microscopic aspect ε u U ε Macroscopic aspect { } U = Calss u ε, 0 < ε < 1 Generalized function
42 Insight into generalized functions algebra, 2 Construction of the algebra The reservoir of representatives : c Es (Ω) = R C (] 0,1] Ω, IR) / D, n IN, c > 0, x Ω, ( DR)(ε, x) n ε The ideal I s (Ω) = { R C (] 0,1] Ω, IR) / D, q IN, cq > 0 / x Ω, ( DR)(ε, x) cqε q } Definition of the Algebra Gs (Ω ) = E s (Ω) I s (Ω )
43 Insight into generalized functions algebra, 3 Some definitions and properties g1, g 2 Gs (Ω), R1, R 2 Es (Ω) / g1 = R1 et g 2 = R 2 The association operation g1, g 2 are said to be associated ( g1 g 2 ) if and only if for two representatives R1, R2 [ ( R (ε, x) R (ε, x))ϕ ( x)] dx 0 quand ε 0 Ω 1 2 D(Ω) The space of ϕ D(Ω) C function with a compact support
44 Insight into generalized functions algebra, 4 Some definitions and properties Product of two generalized functions g1 g 2 = R1 R 2 Derivatives of a generalized function g = R (ε,.)
45 Insight into generalized functions algebra, 5 Some definitions and properties Regularized derivatives g = ( R(ε,.)) ρ h (ε ) ρ D(Ω) or step functions / ρdx = 1, ρ h (ε ) = 1 ρ ( x / h(ε )) n h (ε ) h : ] 0,1] ] 0,1] / h(ε ) 0 qd ε 0
46 Insight into generalized functions algebra, 6 Some definitions and properties Injective mapping L (Ω) Gs (Ω) f L (Ω) f ρε ( x) Es (Ω) How distributions are related to generalized functions T D' (Ω), g Gs (Ω) / g T J.F. Colombeau, Multiplication of Distribution, Lecture Notes in Mathematics, Vol. 1532, Springer-Verlag, 1992.
47 An existence and uniqueness result of the proposed PDE within the generalized functions algebra, 1 A rough Idea of the proof Consider the 1D case ut + a ( x) F (u xx, u x ) x f (u ) = 0 ( S ) (t, x) IR Ω u (0, x) = u0, Interpretation of (S) in u0ε = u0 ρε a ε = a ρε Gs (Ω) U 0 = u0ε, A = a ε U t + A( x) F ( xxu, xu ) x f (U ) 0 ( SG ) (t, x) IR Ω U (0, x) = U 0,
48 An existence and uniqueness result of the proposed PDE within the generalized functions algebra, 2 Proving the existence of a generalized function by proving the existence of a representative u ε t + a ε ( x) F ( xx u ε, x u ε ) x f (u ε ) = 0 ( S E ) ε u (0, x) = u0ε, (t, x) IR Ω Show that exists and ε Show the uniqueness u U1 and U 2 and then u ε Es are two solutions of ODE equation U = uε ( S G ), then U1 U 2
49 An existence and uniqueness result of the proposed PDE within the generalized functions algebra, 3 Theorem [1,2]: Under the assumptions: u L ( Ω ) L (Ω ) 0 1 a L (Ω) f, G C (Ω) L (Ω) Equation (Sg) admits a unique solution in Gs ( IR + Ω) Proposition [3]: The shock-filter equation (S) admits a steady-state u ε L (Ω)
50 Numerical Analysis of the PDE, 1 Results Finite difference scheme u i +1ε, n 2u i ε, n + u i 1ε, n ε, n +1 n h2 ui =u i rmax 0, a i F ε, n ε,n u u i 1 i h u i +1ε, n 2u i ε, n + u i 1ε, n 2 h rmin 0, a i F ε,n ε,n u u i i + 1 h ui ε,0 = u ε,0 (ih), ε,n ε,n ε,n f ' ( u ) u u i 1 i i (, ε,n ε,n ε,n f ' ( u ) u u i i i +1 ( ) )
51 Numerical Analysis of the PDE, 1 Results Proposition : Under the assumptions a BV (Ω) L (Ω) 1 2 The previous scheme is stable for the norm L, for the total variation in space and for the total variation in time in Tonnelli-Cesari sense (i ) u iε,n u 0 L ( ℜ ), i Z, n IN ε,n ε,n (ii ) u i +1 u i TV (u 0 ), n IN i Z (iii ) u ε,n +1 u ε,n TV (u ), n IN i i 0 i Z CFL condition r aff ' <
52 References [1] L. Remaki and M. Cheriet, Enhanced and Restored Signal as a Generalized Solution for Shock Filter Models: Part-I: Existence and Uniqueness result of the Cauchy Problem, Journal of Mathematical Analysis And Applications. Vol. 279, No 1, pp , [2] L. Remaki and M. Cheriet, Enhanced and Restored Signal as a Generalized Solution for Shock Filter Models: Part-II: Numerical Study, Journal of Mathematical Analysis And Applications. Vol. 279, No 2, pp , [3] L. Remaki and W.G. Habashi 3D-mesh adaptation on multiple weak shocks and boundary layers, SIAM J. Sci. Comupt. Vol. 28, No. 4, pp
53 Behaviour of the Shock-filter : 1D case
54 Application to Mesh Adaptation, 1 The scalar solution of adaptation The steady-state of the shock filter equation is a piecewise constant function. Initial solution Steady state Solution of adaptation = Average solution ((original + steady state)/2)
55 Application to Mesh Adaptation, 2 Shock-filter effect on a transonic profile. Transonic Flow around NACA0012: Mach profile The shock filter effect
56 Numerical tests: Inviscid supersonic flow on a ramp: AoA=0, Mach =2 Initial grid Mach profil
57 Numerical tests: Inviscid supersonic flow on a ramp: AoA=0, Mach =2 Grid after 2 cycles of adaptation without pre-processing Grid after 2 cycles of adaptation with preprocessing Mach profile Mach profile
58 Numerical tests: Inviscid supersonic flow on a ramp: AoA=0, Mach =2 Corresponding residuals curves
59 Impact of the Shock filter pre-processing Inviscid test case: NACA 0012, AoA=1, Mach = 0.85 Original grid Cp profile
60 Adaptation without shock-filter preprocessing. The different cycles and corresponding grids The different adapted grids
61 Adaptation without shock-filter preprocessing. The different cycles and corresponding Cp profiles, 1 Different Cp profiles
62 Adaptation without shock-filter preprocessing. The different cycles and corresponding Cp profiles, 2 The chordwise pressure distributions
63 Adaptation with shock-filter pre-processing. The different cycles and corresponding grids and Cp profiles, 1 Different grids and corresponding Cp profiles
64 Adaptation with shock-filter pre-processing. The different cycles and corresponding grids and Cp profiles, 2 The chordwise pressure distributions
65 Inviscid test case: NACA 0012, AoA=1.25, Mach = 0.8 Initial grid Cp profile
66 Adaptation without shock-filter preprocessing. The different cycles and corresponding grids Figure 10-a: Different adapted grids Different adapted grids
67 Adaptation without shock-filter preprocessing. The different cycles and corresponding Cp profiles, 1 Different Cp profiles
68 Adaptation without shock-filter preprocessing. The different cycles and corresponding Cp profiles, 2 The chordwise pressure distributions
69 Adaptation with shock-filter pre-processing. The different cycles and corresponding grids and Cp profiles
Large 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 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 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 informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
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
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
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 informationFinite 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 informationNumerical Solution of Partial Differential Equations governing compressible flows
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
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 informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
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 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 informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
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 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 informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationImplicit 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 informationGetting 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 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 informationBasic 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 informationChapter 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 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 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 informationA 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 informationChapter 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 informationA Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization
A Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization Juan J. Alonso Department of Aeronautics & Astronautics Stanford University jjalonso@stanford.edu Lecture 19 AA200b Applied Aerodynamics
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
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 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 informationProgress in Mesh-Adaptive Discontinuous Galerkin Methods for CFD
Progress in Mesh-Adaptive Discontinuous Galerkin Methods for CFD German Aerospace Center Seminar Krzysztof Fidkowski Department of Aerospace Engineering The University of Michigan May 4, 2009 K. Fidkowski
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
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 informationA 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 informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
More informationChapter 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 informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
More informationApplied 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 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 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 informationChapter 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 informationAA210A 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 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 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 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 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 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 informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
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 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 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 informationNumerical 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 informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationLecture: 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 informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
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 informationDiscrete and continuous adjoint method for compressible CFD J. Peter ONERA
Discrete and continuous adjoint method for compressible CFD J. Peter ONERA J. Peter 1 1 ONERA DMFN October 2, 2014 J. Peter (ONERA DMFN) Adjoint method for compressible CFD October 2, 2014 1 / 60 Outline
More informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
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 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 informationAnisotropic mesh refinement for discontinuous Galerkin methods in two-dimensional aerodynamic flow simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2008; 56:2111 2138 Published online 4 September 2007 in Wiley InterScience (www.interscience.wiley.com)..1608 Anisotropic
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 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 informationShape Optimisation of Axisymmetric Scramjets
Shape Optimisation of Axisymmetric Scramjets Hideaki Ogawa 11 th International Workshop on Shock Tube Technology Australian Hypersonics Workshop 2011 University of Queensland, Brisbane 23 rd March 2011
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 information3.5 Vorticity Equation
.0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous
More informationFundamentals of compressible and viscous flow analysis - Part II
Fundamentals of compressible and viscous flow analysis - Part II Lectures 3, 4, 5 Instantaneous and averaged temperature contours in a shock-boundary layer interaction. Taken from (Pasquariello et al.,
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationCompressible 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 informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationComputational 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 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 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 informationScalar conservation laws with moving density constraints arising in traffic flow modeling
Scalar conservation laws with moving density constraints arising in traffic flow modeling Maria Laura Delle Monache Email: maria-laura.delle monache@inria.fr. Joint work with Paola Goatin 14th International
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
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 information5. 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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
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 informationContinuity 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 informationThe Riemann problem. The Riemann problem Rarefaction waves and shock waves
The Riemann problem Rarefaction waves and shock waves 1. An illuminating example A Heaviside function as initial datum Solving the Riemann problem for the Hopf equation consists in describing the solutions
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationDrag of a thin wing and optimal shape to minimize it
Drag of a thin wing and optimal shape to minimize it Alejandro Pozo December 21 st, 211 Outline 1 Statement of the problem 2 Inviscid compressible flows 3 Drag for supersonic case 4 Example of optimal
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 informationHPC enabling of OpenFOAM R for CFD applications
HPC enabling of OpenFOAM R for CFD applications Numerical acoustic analysis of a turbulent flow around a bluff body 06-08 April 2016, Casalecchio di Reno, BOLOGNA. Marta Cianferra, University of Trieste
More informationEntropy-Based Drag Error Estimation and Mesh. Adaptation in Two Dimensions
Entropy-Based Drag Error Estimation and Mesh Adaptation in Two Dimensions Krzysztof J. Fidkowski, Marco A. Ceze, and Philip L. Roe University of Michigan, Ann Arbor, MI 48109 This paper presents a method
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 information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
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 informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
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 informationPIPE FLOWS: LECTURE /04/2017. Yesterday, for the example problem Δp = f(v, ρ, μ, L, D) We came up with the non dimensional relation
/04/07 ECTURE 4 PIPE FOWS: Yesterday, for the example problem Δp = f(v, ρ, μ,, ) We came up with the non dimensional relation f (, ) 3 V or, p f(, ) You can plot π versus π with π 3 as a parameter. Or,
More informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 014 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
More information