Effects of the Jacobian Evaluation on Direct Solutions of the Euler Equations
|
|
- Steven Lloyd
- 5 years ago
- Views:
Transcription
1 Middle East Technical University Aerospace Engineering Department Effects of the Jacobian Evaluation on Direct Solutions of the Euler Equations by Ömer Onur Supervisor: Assoc. Prof. Dr. Sinan Eyi
2 Outline Introduction Objectives Flow Analysis Accuracy Analysis & Results Performance Analysis & Results Conclusion Recommendations 2
3 Introduction Steady flow computations can be realized either by iterative means,, or using direct methods. Although iterative solvers make an unsteady flow analysis with an advance in time,, with direct methods a steady flow analysis is possible. Low memory requirements made the iterative methods popular untill the development of recent advanced computers. Solving the whole domain at once makes direct solvers stable, and faster convergence is possible due to small number of iterations. In CFD applications like design optimization,, and flutter analysis,, direct methods are preferable. 3
4 Introduction (Cont.) Direct solvers require the calculation of Jacobian matrix. Derivation of analytical Jacobians becomes more difficult as the discretization of governing equations become more complex. The best alternative is to compute the Jacobians numerically as accurate as possible. 4
5 Objectives To develop a direct flow solver code To compare the accuracy of numerical and analytical Jacobians used in direct flow solvers To investigate their effects on the performance (convergence and CPU time) of direct flow solvers To improve the efficiency of Jacobian matrix solution 5
6 Flow Analysis Flow code is a 2-D planar/axisymmetric Euler solver: 1st/2nd order finite-volume discretizations Steger-Warming, Van Leer, Roe upwind flux splitting schemes Newton's method both numerical and analytical Jacobian matrices UMFPACK sparse matrix solver package Different geometries,, grid sizes and BCs are considered. 6
7 Flow Analysis (Cont.) Steady, 2-D planar/axisymmetric Euler equations in generalized coordinates: F(W ˆ ˆ ) G(W ˆ ˆ ) + + σ H(W ˆ ˆ ) = 0 ξ η where Ŵ = J 1 ρ ρu ρv ρet ˆF = J 1 ρu ρuu + ξ x p ρvu + ξ y p ( ρet + p)u Ĝ = J 1 ρv ρuv + ηx p ρvv + ηy p ( ρet + p)v Ĥ = J 1 1 y ρv ρuv 2 ρv ( ρet + p)v 7
8 Flow Analysis (Cont.) Newton s method: R(W ˆ ˆ ) = 0 (Discretized Residual) n ˆR ˆ n ˆ n W R(W ) Ŵ = Wˆ = Wˆ + Wˆ n+ 1 n n [Jacobian Matrix] 8
9 Flow Analysis (Cont.) Analytical Jacobians: Derivation by hand / symbolic manipulators for simple flow models Difficult derivation for complex flow models Numerical Jacobians: Rˆ R(Wˆ ˆ i i + ε e j ) R(W) i = Ŵ ε j Good choice of finite-difference perturbation magnitude ε Usage of higher computer precision 9
10 Accuracy Analysis Error analysis: Condition error [loss of numerical precision] Truncation error [neglected terms in the Taylor series] Total error 2 2 ER f ( ξ ) ε E TOTAL( ε ) = E C( ε ) + E T ( ε ) = + 2 ε x 2 ε<< E C >> & ε >> E T >> ε opt E TOTAL min 10
11 Accuracy Analysis (Cont( Cont.) Optimum perturbation magnitude analysis: E ( ε ) 2 E 1 2 TOTAL R f( ξ ) = + = 0 ε ε 2 2 x 2 ε = 2 OPT 2 ER f( ξ ) 2 x E R can be taken as ε M or simply found by subtraction of double and single precision calculations of flux values. Second derivative of flux can be calculated using a forward/backward finite difference method with double precision or taken as equal to 1. Optimum perturbation magnitude is also found by trialerror 11
12 Accuracy Results Error in numerical flux Jacobians are analyzed for 15º ramp geometry with 33x25 grid M=2.0 free-stream flow inlet-outlet, symmetry and wall BCs Jacobians are calculated using already converged solution 1st order S-W flux splitting Different perturbation magnitude values are investigated for both single and double precision. Optimum perturbation magnitude analysis is performed for single precision. 12
13 Accuracy Results (Cont.) Effect of Control on Total Errors for Residual Jacobians [Single Precision, Forward Differencing] TotalError between ResidualJacobians ( R/ W) w/o Cont. (Max. E) Effect of Control on Total Errors for Residual Jacobians [Single Precision,Backward Differencing] ε =7 wcont.(max.e) w/o Cont. (Avg. E) w Cont. (Avg. E) TotalError between ResidualJacobians ( R/ W) ε = w/o Cont. (Max. E) wcont.(max.e) w/o Cont. (Avg. E) w Cont. (Avg. E) Perturbation Magnitude (ε) Perturbation Magnitude (ε)
14 Accuracy Results (Cont.) Change of Total Maximum Error with Perturbation Magnitude [Single Precision, Forward Differencing] TotalMax. Error between Jacobians F + / W Change of Total Average Error with Perturbation Magnitude [Single Precision, Forward Differencing] ε =7 F - / W G + / W G - / W R pl / W Total Avg. Error between Jacobians F + / W ε =7 F - / W G + / W G - / W R pl / W Perturbation Magnitude (ε) Perturbation Magnitude (ε)
15 Accuracy Results (Cont.) Effect of Axisymmetry on Total Errors for Residual Jacobians [Single Precision, Forward Differencing] Optimum Perturbation Magnitude (ε opt ) Analysis for Single precision TotalError between ResidualJacobians ( R/ W) 10 4 Planar (Avg. E) 10 3 Axis ym. (Avg. E) Planar (Max. E) Axisym. (Max. E) ε =7 Trial-Error Procedure Optimization Method Flux ε opt (max. error) ε opt (avg. error) ε opt (avg.) ε opt ( ε M ) F F G G For single precision; ε opt 7x Perturbation Magnitude (ε)
16 Accuracy Results (Cont.) Effect of Control on Total Errors for Residual Jacobians [Double Precision, Forward Differencing] Effect of Control on Total Errors for Residual Jacobians [Double Precision, Backward Differencing] TotalError between ResidualJacobians ( R/ W) 10 3 w/o Cont. (Max. E) wcont.(max.e) 10 1 w/o Cont. (Avg. E) w Cont. (Avg. E) ε = TotalError between ResidualJacobians ( R/ W) 10 3 w/o Cont. (Max. E) wcont.(max.e) 10 1 w/o Cont. (Avg. E) w Cont. (Avg. E) ε = Perturbation Magnitude (ε) Perturbation Magnitude (ε)
17 Accuracy Results (Cont.) Change of Total Maximum Error with Perturbation Magnitude [Double Precision, Forward Differencing] TotalMax. Error between Jacobians F + / W Change of Total Average Error with Perturbation Magnitude [Double Precision, Forward Differencing] ε = F - / W G + / W G - / W R pl / W TotalAvg. Error between Jacobians F + / W ε = F - / W G + / W G - / W R pl / W Perturbation Magnitude (ε) Perturbation Magnitude (ε)
18 Accuracy Results (Cont.) Effect of Axisymmetry on Total Errors for Residual Jacobians [Double Precision, Forward Differencing] Optimum Perturbation Magnitude (ε opt ) Analysis for Double precision TotalError between ResidualJacobians ( R/ W) 10 0 Planar (Avg. E) 10-1 Axis ym. (Avg. E) Planar (Max. E) Axisym. (Max. E) ε = Trial-Error Procedure Optimization Method Flux ε opt (max. error) ε opt (avg. error) ε opt F F G G For double precision; ε opt 4x Perturbation Magnitude (ε)
19 Performance Analysis Jacobian Matrix structure: a huge square matrix even for a simple geometry composed of either analytical or numerical Jacobians most of the elements are zero Matrix solver: Direct full matrix solvers have very high storage and memory costs. Storage of only non-zero elements improve the efficiency greatly. Sparse Matrix Solvers 19
20 Performance Analysis (Cont.) Matrix Solution Strategies: Frozen Jacobian [using same Jacobian matrix in subsequent iterations] Good initial guess [a time-like term addition to the Jacobian matrix diagonal] 0 n ˆ R(W ˆ ) 1 R n 0 n n 2 [] I + Wˆ = R(W ˆ ) t = t t Wˆ ˆ n R(W ) 2 Δt>> Original Newton s method 20
21 Performance Results Test Case Results: Test Case Results: Convergence and CPU time performance of the direct solver is analyzed for 15º ramp geometry with 33x25 grid M=2.0 free-stream flow inlet-outlet, symmetry and wall BCs Jacobians are calculated using 1st order S-W flux splitting Calculations are realized with forward and backward differencing in single and double precision. Effect of perturbation magnitude ε is investigated Effect of Jacobian freezing and time-like diagonal term addition are analyzed. 21
22 Performance Results (Cont.) 33x25 Grid for Ramp geometry 34x26 cells 4 flow variables total 3536 variables Jacobian matrix has million elements Solver Storage CPU time Full Matrix LU 240 MB 26 h Sparse Matrix 1.2 MB 20 s 22
23 Performance Results (Cont Cont.) Mach Contour for Ramp geometry [2-D planar, Double Precision, 1st order S-W] Mach Contour for Ramp geometry [2-D axisymmetric, Double Precision, 1st order S-W]
24 Performance Results (Cont.) Effect of t onconvergencehistory [2-D planar, Double Precision] Effect of t onconvergencehistory [2-D axisymmetric, Double Precision] 10 0 AnalyticalJac. No frz AnalyticalJac. No frz. Max. Density Residual (R 1 ) Max. Density Res idual (R 1 ) No t t rm =1.5 t rm =5. Full t #ofiterations t CPU Time (s) No t NaN t rm = t rm = Full t No t t rm =1.5 t rm =5. Full t #ofiterations t CPU Time (s) No t 9.79 t rm = t rm = Full t
25 Performance Results (Cont.) 10 1 EffectOfFreezingonConvergenceHistory [2-D planar, Double Precision] 10 1 Effect Of Freezing on Convergence History [2-D axisymmetric, Double Precision] 10-1 Analytical J ac. t rm = Analytical Jac. No t Max. Density Res idual (R 1 ) Max. Density Residual (R 1 ) No frz. R frz = R frz =1 R frz = No frz. R frz = R frz =1 R frz = #ofiterations Tolerance CPU Time (s) No freeze R frz = 1x R frz = 1x R frz = 1x #ofiterations Tolerance CPU Time (s) No freeze 9.79 R frz = 1x R frz = 1x 8.15 R frz = 1x
26 Performance Results (Cont.) 10 0 Convergence History [2-D planar, Double Precision, Forward Differencing] t rm =1.5 R frz = Convergence History [2-D a xis ymmetric, Double P re cis ion, Forward Diffe rencing] t rm =0. R frz = 10-3 Max. Density Residual (R 1 ) Max. Density Residual (R 1 ) An. J a c. ε =4 ε = ε = #ofiterations Jacobian CPU Time (s) Analytic ε = 4x ε = 4x ε = 4x An. J ac. ε =4 ε = ε = #ofiterations Jacobian CPU Time (s) Analytic 8.25 ε = 4x ε = 4x ε = 4x
27 Performance Results (Cont.) Different Flux Splitting Schemes: Different Flux Splitting Schemes: Convergence and CPU time performance of the direct solver is analyzed for the test case using 1st order Van Leer/Roe flux splitting Jacobians are calculated using 1st order Van Leer/Roe flux splitting numerically 1st order S-W flux splitting analytically Optimum perturbation magnitude ε opt is used in the calculations. Calculations are realized with forward and backward differencing in single and double precision. Benefits of using the same flux calculation scheme for both Jacobian and residual calculation are analyzed. 27
28 Performance Results (Cont Cont.) Mach Contour for Ramp geometry [2-D planar, Double Precision, 1st order VL] Mach Contour for Ramp geometry [2-D axisymmetric, Double Precision, 1st order VL]
29 Performance Results (Cont Cont.) Mach Contour for Ramp geometry [2-D planar, Double Precision, 1st order Roe] Mach Contour for Ramp geometry [2-D axisymmetric, Double Precision, 1st order Roe]
30 Performance Results (Cont.) 10 0 Convergence History [2-D planar, Double Precision, S-W An. Jac.] t rm =1.5 R frz = Convergence History [2-D planar, Double Precision, Num. Jac.] t rm =1.5(SW) t rm =2.(VL) t rm =3.(Roe) R frz = ε = Max. Density Residual (R 1 ) Max. Density Residual (R 1 ) Steger-Warming Van Leer Roe #ofiterations FS CPU Time (s) SW VL Roe Steger-Warming Van Leer Roe #ofiterations FS CPU Time (s) SW VL Roe
31 Performance Results (Cont.) 10 0 Convergence History [2-D axisymmetric, Double Precision, S-W An. Jac.] t rm =0. R frz = Convergence History [2-D axisymmetric, Double Precision, Num. Jac.] t rm =0.(SW,VL) t rm =1.5(Roe) R frz = ε = Max. Density Residual (R 1 ) Max. Density Residual (R 1 ) Steger-Warming Van Leer Roe #ofiterations FS CPU Time (s) SW 8.25 VL Roe Steger-Warming Van Leer Roe #ofiterations FS CPU Time (s) SW 9.72 VL Roe
32 Performance Results (Cont.) Higher-Order Discretizations: Higher-Order Discretizations: Convergence and CPU time performance of the direct solver is analyzed for the test case using 2nd order S-W/Van Leer/Roe flux splitting Van Albada limiter Jacobians are calculated using 2nd order S-W/Van Leer/Roe flux splitting numerically Optimum perturbation magnitude ε opt is used in the calculations. Calculations are realized with forward differencing in double precision. 32
33 Ömer Onur - AE500 Presentation Performance Results (Cont Cont.) Mach Contour for Ramp geometry [2-D planar, Double Precision, 2nd order S-W] Mach Contour for Ramp geometry [2-D axisymmetric, Double Precision, 2nd order S-W]
34 Ömer Onur - AE500 Presentation Performance Results (Cont Cont.) Mach Contour for Ramp geometry [2-D planar, Double Precision, 2nd order VL] Mach Contour for Ramp geometry [2-D axisymmetric, Double Precision, 2nd order VL]
35 Ömer Onur - AE500 Presentation Performance Results (Cont Cont.) Mach Contour for Ramp geometry [2-D planar, Double Precision, 2nd order Roe] Mach Contour for Ramp geometry [2-D axisymmetric, Double Precision, 2nd order Roe]
36 Performance Results (Cont.) Convergence History [2-D planar, Double Precision, Num. Jac., Van Albada] Full t No Freezing ε = Convergence History [2-D axisymmetric, Double Precision, Num. Jac., Van Albada] Full t No Freezing ε = Max. Density Residual (R 1 ) Max. Density Residual (R 1 ) Steger-Warming Van Leer Roe #ofiterations FS CPU Time (s) SW VL Roe Steger-Warming Van Leer Roe #ofiterations FS SW VL Roe CPU Time (s) Limit cycle Limit cycle Limit cycle 36
37 Performance Results (Cont.) Different Geometry and Flow Conditions: Different Geometry and Flow Conditions: Convergence and CPU time performance of the direct solver is analyzed for bump geometry with 65x17 grid M=2.0 and M=0.5 free-stream flows inlet-outlet, symmetry and wall BCs 2nd order Roe flux splitting with Van Albada limiter Jacobians are calculated using 2nd order Roe flux splitting numerically Optimum perturbation magnitude ε opt is used in the calculations. Calculations are realized with forward differencing in double precision. 37
38 Performance Results (Cont.) 65x17GridforBumpgeometry 38
39 Performance Results (Cont.) 39
40 Performance Results (Cont.) Convergence History [2-D planar, Double Precision, Num. Jac., Van Albada] Full t No Freezing ε = Convergence History [2-D axisymmetric, Double Precision, Num. Jac., Van Albada] Full t No Freezing ε = Max. Density Res idual (R 1 ) Max. Density Res idual (R 1 ) Roe Roe #ofiterations #ofiterations 40
41 Performance Results (Cont.) 41
42 Performance Results (Cont.) Convergence History [Subsonic 2-D planar, Double Precision, Num. Jac., Van Albada] Full t No Freezing ε = Convergence History [Subsonic 2-D axisymmetric, Double Precision, Num. Jac., Van Albada] Full t No Freezing ε = Max. Density Res idual (R 1 ) Max. Density Res idual (R 1 ) Roe Roe #ofiterations #ofiterations 42
43 Conclusion A direct flow solver code is developed. Accuracy of numerical Jacobians used in the solver is analyzed. A control mechanism is required in 1st order SW Jacobian calculation. Forward/backward differencing does not differ so much. Double precision improves accuracy significantly. Optimum perturbation magnitude is found by an optimization method. 43
44 Conclusion (Cont( Cont.) The effects of the accuracy of numerical Jacobians on the performance (convergence and CPU time) of direct flow solvers is investigated. Double precision improves convergence limits significantly. Optimum perturbation magnitude gives the same convergence with the analytical method. Calculation of fluxes with perturbation only for related cells reduced the execution time greatly. 44
45 Conclusion (Cont( Cont.) The efficiency of Jacobian matrix solution is improved by some strategies. A sparse matrix solver having low storage and memory requiremets is used. A time-like term addition to matrix diagonal stabilizes the solver in case of poor initial conditions. Removing this modification at the right time makes the convergence very fast. Jacobian matrix freezing at the right time decreases the execution time. 45
46 Conclusion (Cont( Cont.) Using the same flux calculation scheme for both Jacobian and residual calculation maintains faster convergence. According to the choice of the flux limiter, higherorder schemes may have convergence problems. Unsteady results like limit-cycle is observed in higherorder schemes. 46
47 Recommendations The control mechanism used in 1st order SW Jacobian calculation can be improved to give the best accuracy for all perturbation magnitude values. More advanced strategies can be considered to improve the Jacobian matrix solution. Higher-order schemes, and especially the choice and application of limiter can be analyzed deeply to obtain fully converged solution. 47
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
More informationParallel Computation of Forced Vibration for A Compressor Cascade
44th AIAA Aerospace Sciences Meeting and Exhibit 9-2 January 26, Reno, Nevada AIAA 26-628 AIAA Paper 26-628 Parallel Computation of Forced Vibration for A Compressor Cascade Zongjun Hu and Gecheng Zha
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 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 informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
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 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 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 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 informationSIMPLE Algorithm for Two-Dimensional Channel Flow. Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Channel Flow Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
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 informationCFD with OpenSource software
CFD with OpenSource software The Harmonic Balance Method in foam-extend Gregor Cvijetić Supervisor: Hrvoje Jasak Faculty of Mechanical Engineering and Naval Architecture University of Zagreb gregor.cvijetic@fsb.hr
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More informationY. Abe, N. Iizuka, T. Nonomura, K. Fujii Corresponding author:
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2012-2801 ICCFD7-2801 Symmetric-conservative metric evaluations for higher-order finite
More 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 informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationEffect of modification to tongue and basic circle diameter on vibration in a double-suction centrifugal pump
5th International Conference on Information Engineering for Mechanics and Materials (ICIMM 2015) Effect of modification to tongue and basic circle diameter on vibration in a double-suction centrifugal
More informationFOR many computational fluid dynamics (CFD) applications, engineers require information about multiple solution
AIAA Aviation 13-17 June 2016, Washington, D.C. 46th AIAA Fluid Dynamics Conference AIAA 2016-3809 Improved Functional-Based Error Estimation and Adaptation without Adjoints William C. Tyson *, Katarzyna
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 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 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 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 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 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 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 informationImplementation and Comparisons of Parallel Implicit Solvers for Hypersonic Flow Computations on Unstructured Meshes
20th AIAA Computational Fluid Dynamics Conference 27-30 June 2011, Honolulu, Hawaii AIAA 2011-3547 Implementation and Comparisons of Parallel Implicit Solvers for Hypersonic Flow Computations on Unstructured
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 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 informationExplicit Incompressible Euler Solver
Explicit Incompressible Euler Solver (Written In Matlab) R D Teja Aerospace Undergraduate IIT Madras Contents: 1. Abstract 2. Introduction 2.1 Euler equations 2.2 Artificial compressibility equations 2.3
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationIntroduction to numerical simulation of fluid flows
Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural
More informationPOSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS. Cheng Wang. Jian-Guo Liu
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 3 Number May003 pp. 0 8 POSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
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 informationComparison of some approximation schemes for convective terms for solving gas flow past a square in a micorchannel
Comparison of some approximation schemes for convective terms for solving gas flow past a square in a micorchannel Kiril S. Shterev and Sofiya Ivanovska Institute of Mechanics, Bulgarian Academy of Sciences,
More informationCFD ANALYSIS OF THE CORK-PHENOLIC HEAT SHIELD OF A REENTRY CUBESAT IN ARC-JET CONDITIONS INCLUDING ABLATION AND PYROLYSIS
CFD ANALYSIS OF THE CORK-PHENOLIC HEAT SHIELD OF A REENTRY CUBESAT IN ARC-JET CONDITIONS INCLUDING ABLATION AND PYROLYSIS 15 th International Planetary Probe Workshop 14/06/2018 Ata Onur Başkaya 1 A. Turchi
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 information8. Introduction to Computational Fluid Dynamics
8. Introduction to Computational Fluid Dynamics We have been using the idea of distributions of singularities on surfaces to study the aerodynamics of airfoils and wings. This approach was very powerful,
More informationDefinition of Ablation Test-case series #3
Tom van Eekelen Alexandre Martin Jean Lachaud Daniele Bianchi Samtech, A Siemens Industry Software Company, Belgium University of Kentucky, USA UARC/Univ. of California Santa Cruz, USA La Sapienza Rome,
More informationHopf Bifurcation Calculations for a Symmetric Airfoil in Transonic Flow
AIAA JOURNAL Vol. 42, No. 5, May 2004 Hopf Bifurcation Calculations for a Symmetric Airfoil in Transonic Flow K. J. Badcock, M. A. Woodgate, and B. E. Richards University of Glasgow, Glasgow, Scotland
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationReduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis
Journal of Mechanical Science and Technology 24 (2) (2010) 601~608 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-010-0106-0 Reduced order model of three-dimensional Euler equations using proper
More information6. Iterative Methods for Linear Systems. The stepwise approach to the solution...
6 Iterative Methods for Linear Systems The stepwise approach to the solution Miriam Mehl: 6 Iterative Methods for Linear Systems The stepwise approach to the solution, January 18, 2013 1 61 Large Sparse
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 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 informationNumerical Analysis: Solutions of System of. Linear Equation. Natasha S. Sharma, PhD
Mathematical Question we are interested in answering numerically How to solve the following linear system for x Ax = b? where A is an n n invertible matrix and b is vector of length n. Notation: x denote
More informationC. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2
On iterative solvers combined with projected Newton methods for reacting flow problems C. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2 1 Delft University of Technology J.M. Burgerscentrum 2 Delft University
More informationThe Euler Equations! Advection! λ 1. λ 2. λ 3. ρ ρu. c 2 = γp/ρ. E = e + u 2 /2 H = h + u 2 /2; h = e + p/ρ. 0 u 1/ρ. u p. t + A f.
http://www.nd.edu/~gtryggva/cfd-course/! Advection! Grétar Tryggvason! Spring! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = x ( / ) ρe ρu E + p ρ E = e + u / H = h + u
More informationSolution Methods In Computational Fluid Dynamics
Solution Methods In Computational Fluid Dynamics T. H. Pulliam NASA Ames Research Center, MS T27B-1 Moffett Field CA, 94035 E-mail: tpulliam@mail.arc.nasa.gov Key Words: Euler, Navier-Stokes, Finite Differences,
More informationOverview of Convection Heat Transfer
Overview of Convection Heat Transfer Lecture Notes for ME 448/548 Gerald Recktenwald March 16, 2006 1 Introduction CFD is a powerful tool for engineering design involving fluid flow and heat transfer.
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationApplication of Chimera Grids in Rotational Flow
CES Seminar Write-up Application of Chimera Grids in Rotational Flow Marc Schwalbach 292414 marc.schwalbach@rwth-aachen.de Supervisors: Dr. Anil Nemili & Emre Özkaya, M.Sc. MATHCCES RWTH Aachen University
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 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 informationAffordable, entropy-consistent, Euler flux functions
Affordable, entropy-consistent, Euler flux functions (with application to the carbuncle phenomenon) Phil Roe Aerospace Engineering University 0f Michigan Ann Arbor Presented at HYP 2006 1 Entropy pairs
More informationC. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2
s: application to chemical vapor deposition C. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2 1 Delft University of Technology J.M. Burgerscentrum 2 Delft University of Technology Department of Multi Scale Physics
More informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More informationA multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction
A multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction J Degroote, S Annerel and J Vierendeels Department of Flow, Heat and Combustion Mechanics, Ghent University,
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 Contours of Pressure.8.6.4...4.6.8 Matthew McMullen Antony Jameson Juan J. Alonso Aero Astro Dept. Stanford
More informationNUMERICAL METHOD FOR THREE DIMENSIONAL STEADY-STATE TWO-PHASE FLOW CALCULATIONS
' ( '- /A NUMERCAL METHOD FOR THREE DMENSONAL STEADY-STATE TWO-PHASE FLOW CALCULATONS P. Raymond,. Toumi (CEA) CE-Saclay DMT/SERMA 91191 Gif sur Yvette, FRANCE Tel: (331) 69 08 26 21 / Fax: 69 08 23 81
More informationNumerical Simulation of Jet Aerodynamics Using the Three-Dimensional Navier-Stokes Code PAB3D
NASA Technical Paper 3596 Numerical Simulation of Jet Aerodynamics Using the Three-Dimensional Navier-Stokes Code PAB3D S. Paul Pao Langley Research Center Hampton, Virginia Khaled S. Abdol-Hamid Analytical
More informationNewton-Krylov-Schwarz Method for a Spherical Shallow Water Model
Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model Chao Yang 1 and Xiao-Chuan Cai 2 1 Institute of Software, Chinese Academy of Sciences, Beijing 100190, P. R. China, yang@mail.rdcps.ac.cn
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf-Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2018/19 Part 4: Iterative Methods PD
More informationScientific Computing
Scientific Computing Direct solution methods Martin van Gijzen Delft University of Technology October 3, 2018 1 Program October 3 Matrix norms LU decomposition Basic algorithm Cost Stability Pivoting Pivoting
More informationComputing Spectra of Linear Operators Using Finite Differences
Computing Spectra of Linear Operators Using Finite Differences J. Nathan Kutz Department of Applied Mathematics University of Washington Seattle, WA 98195-2420 Email: (kutz@amath.washington.edu) Stability
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 informationJACOBI S ITERATION METHOD
ITERATION METHODS These are methods which compute a sequence of progressively accurate iterates to approximate the solution of Ax = b. We need such methods for solving many large linear systems. Sometimes
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 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 informationFinite volume method for CFD
Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor:
More informationMethods for Computing Periodic Steady-State Jacob White
Introduction to Simulation - Lecture 15 Methods for Computing Periodic Steady-State Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Periodic Steady-state problems Application
More informationImproving the Verification and Validation Process
Improving the Verification and Validation Process Mike Fagan Rice University Dave Higdon Los Alamos National Laboratory Notes to Audience I will use the much shorter VnV abbreviation, rather than repeat
More informationIterative Solution methods
p. 1/28 TDB NLA Parallel Algorithms for Scientific Computing Iterative Solution methods p. 2/28 TDB NLA Parallel Algorithms for Scientific Computing Basic Iterative Solution methods The ideas to use iterative
More informationAEROELASTIC SYSTEM IDENTIFICATION USING TRANSONIC CFD DATA FOR A 3D WING
AEROELASTIC SYSTEM IDENTIFICATION USING TRANSONIC CFD DATA FOR A 3D WING G.A. Vio, J.E. Cooper, G. Dimitriadis School of Mechanical, Aerospace and Civil Engineering, The University of Manchester, Manchester
More informationUniversity of Bristol - Explore Bristol Research. Peer reviewed version. Link to published version (if available): /6.
Bagheri, A. K., Jones, D. P., & Gaitonde, A. L. (2016). Linear frequency domain reduced order modelling of aerofoil gust response. In 46th AIAA Fluid Dynamics Conference. American Institute of Aeronautics
More informationCFD Analysis of Mixing in Polymerization Reactor. By Haresh Patel Supervisors: Dr. R. Dhib & Dr. F. Ein-Mozaffari IPR 2007
CFD Analysis of Mixing in Polymerization Reactor By Haresh Patel Supervisors: Dr. R. Dhib & Dr. F. Ein-Mozaffari Introduction Model development Simulation Outline Model Setup for Fluent Results and discussion
More informationChapter 5. The Finite Volume Method for Convection-Diffusion Problems. The steady convection-diffusion equation is
Chapter 5 The Finite Volume Method for Convection-Diffusion Problems Prepared by: Prof. Dr. I. Sezai Eastern Mediterranean University Mechanical Engineering Department Introduction The steady convection-diffusion
More informationA Zooming Approach to Investigate Heat Transfer in Liquid Rocket Engines with ESPSS Propulsion Simulation Tool
A Zooming Approach to Investigate Heat Transfer in Liquid Rocket Engines with ESPSS Propulsion Simulation Tool M. Leonardi, F. Di Matteo, J. Steelant, B. Betti, M. Pizzarelli, F. Nasuti, M. Onofri 8th
More informationArt Checklist. Journal Code: Article No: 6959
Art Checklist Journal Code: JCPH Article No: 6959 Disk Recd Disk Disk Usable Art # Y/N Format Y/N Remarks Fig. 1 Y PS Y Fig. 2 Y PS Y Fig. 3 Y PS Y Fig. 4 Y PS Y Fig. 5 Y PS Y Fig. 6 Y PS Y Fig. 7 Y PS
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationSimulation of Condensing Compressible Flows
Simulation of Condensing Compressible Flows Maximilian Wendenburg Outline Physical Aspects Transonic Flows and Experiments Condensation Fundamentals Practical Effects Modeling and Simulation Equations,
More informationCommissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation. MIT, November 3rd, 2008
ON THE NUMERICAL SIMULATION OF GENERAL MULTI-FIELD TWO-PHASE FLOW MODELS Anela KUMBARO Commissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation MIT, November
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 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 informationMath 471 (Numerical methods) Chapter 3 (second half). System of equations
Math 47 (Numerical methods) Chapter 3 (second half). System of equations Overlap 3.5 3.8 of Bradie 3.5 LU factorization w/o pivoting. Motivation: ( ) A I Gaussian Elimination (U L ) where U is upper triangular
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 informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationDESIGN AND NUMERICAL ANALYSIS OF ASYMMETRIC NOZZLE OF SHCRAMJET ENGINE
ISSN (O): 393-869 DESIGN AND NUMERICAL ANALYSIS OF ASYMMETRIC NOZZLE OF SHCRAMJET ENGINE Lokesh Silwal lokusilwal@gmail.com Sajan Sharma thesajansharma@gm ail.com Nitish Acharya nitishacharya1818@g mail.com
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 informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationVery Simple, Carbuncle-Free, Boundary-Layer-Resolving, Rotated-Hybrid Riemann Solvers
Very Simple, Carbuncle-Free, Boundary-Layer-Resolving, Rotated-Hybrid Riemann Solvers Hiroaki Nishikawa 1 and Keiichi Kitamura 2 1,2 W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department
More informationSIMULATION OF UNSTEADY TURBOMACHINERY FLOWS USING AN IMPLICITLY COUPLED NONLINEAR HARMONIC BALANCE METHOD
Proceedings of ASME Turbo Expo 2011 GT2011 June 6-10, 2011, ancouver, Canada GT2011-46367 SIMULATION OF UNSTEADY TURBOMACHINERY FLOWS USING AN IMPLICITLY COUPLED NONLINEAR HARMONIC BALANCE METHOD Jonathan
More informationA Computational Study on the Thrust Performance of a Supersonic Pintle Nozzle
June 30 - July 3, 2015 Melbourne, Australia 9 P-10 A Computational Study on the Thrust Performance of a Supersonic Pintle Nozzle Ruoyu Deng Department of Mechanical Engineering Andong National University,
More informationEstimation of grid-induced errors in computational fluid dynamics solutions using a discrete error transport equation
Graduate Theses and Dissertations Graduate College 2009 Estimation of grid-induced errors in computational fluid dynamics solutions using a discrete error transport equation Brandon Williams Iowa State
More informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationFinite volume scheme for Euler equations with gravity which is well-balanced for general equations of state and grid systems
Finite volume scheme for Euler equations with gravity which is well-balanced for general equations of state and grid systems Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg,, and Friedrich
More information