Solution of 2D Euler Equations and Application to Airfoil Design
|
|
- Mark Doyle
- 5 years ago
- Views:
Transcription
1 WDS'6 Proceedings of Contributed Papers, Part I, 47 52, 26. ISBN MATFYZPRESS Solution of 2D Euler Equations and Application to Airfoil Design J. Šimák Charles University, Faculty of Mathematics and Physics, Prague, Czech Republic. Abstract. This paper deals with a numerical method for an airfoil design. It is shown how to create an airfoil from a given velocity distribution along a mean camber line. The method is based on searching a fixed point of a contractive operator. We need to have a fast solver of the Euler equations. The Newton method for solving implicit finite volume scheme is described. The resulting system of linear algebraic equations is solved by GMRES, the Jacobian-free version is described. Numerical results are presented. Introduction We present here a method how to get an airfoil from a given velocity distribution. This method can be used in the case of a subsonic inviscid compressible flow described by the Euler equations. The method was presented by [Pelant, 1998]. The main idea is to use a couple of a direct and an inverse operator. The direct operator represents an assignment of a velocity distribution to a shape of an airfoil and the second operator is its inversion. Since we are not able to construct an exact inversion we have to use some iteration method. We put the direct and the approximate inverse operator together to create a contractive operator and find its fixed point. The method can be used to modify existing airfoils. We get a velocity distribution of some airfoil, then we change it and we get a new airfoil with required quality. In the next sections these two operators are described. The theory of a numerical solution of flow dynamics can be found in [Feistauer et al., 23]. Solution of Euler equations using Newton-Krylov method Our goal is to find a steady-state solution of the Euler equations describing inviscid compressible flow. The nonstationary equations can be written in the form where w t + f(w) x + g(w) =, (1) y w = (ρ,ρu,ρv,e) T, f(w) = (ρu,ρu 2 + p,ρuv,(e + p)u) T, g(w) = (ρv,ρuv,ρv 2 + p,(e + p)v) T. The symbol ρ denotes the density, u,v denote the velocity components, E denotes the energy and p denotes the pressure. The Pressure p can be expressed by ( p = (γ 1) E 1 2 ρ( u 2 + v 2)), (2) where γ is the Poisson adiabatic constant. The problem is solved by the implicit finite volume method. We assume a standard finite volume mesh D h (in experiments a quadrilateral structured grid was used). By w h we shall denote a vector of 4-dimensional blocks w i of the values of approximate solution on finite 47
2 volumes D i D h. For w h R n the vector Φ(w h ) consists of 4-dimensional blocks Φ i (w h ) given by Φ i (w h ) = 1 H (w i,w j,n ij ) Γ ij, (3) D i j S(i) where D i denotes the cell area, Γ ij denotes the length of the edge between D i and D j, n ij denotes the outer normal to D i and H(u,v,n) denotes the numerical flux. In this case the Osher- Solomon numerical flux is used. In order to have a higher order method we apply the Van Leer κ-scheme or the Van Albada limiter. By S(i) we denote a set of indices of neighbouring elements and by symbol τ we shall denote the time step. So we get an implicit FV scheme w k+1 i = w k i τ k Φ i (w k+1 h ). (4) Since we are interested in a steady-state solution, we can deal with the nonlinear system Φ(w h ) = and apply the Newton method to this system: DΦ(w k h ) Dw ( ) w k+1 h wh k = Φ(wh k ), k =,1,..., (5) where wh is some suitable initial approximation. Since the employed numerical flux is differentiable, it is possible to compute the Jacobian matrix DΦ(w)/Dw. We will refer to these Newton iterations as outer iterations. We get a large system of linear algebraic equations Ax = b. The matrix A is sparse and consists of four-by-four blocks. Although the structure is symmetric, the matrix itself is nonsymmetric. This leads us to the use of iterative Krylov subspace methods, namely GMRES. As the GMRES solver the program SPARSKIT2 was used. The algorithm can be found in [Saad, 23]. Since we use an iterative solver, it is necessary to stop the process after a sufficient precision is reached and avoid oversolving of the system. The stopping criterion can be residuum < 1 2 Φ(w k h ). To these iterations we will refer as inner iterations. Since the Jacobian matrix is very ill-conditioned, it is necessary to use some preconditioning. ILU(p) preconditioning is used. The choice of p influences the speed of GMRES and also the amount of memory necessary to store the preconditioning matrix. The optimum value appears to be p = 2. In the GMRES algorithm the Jacobian matrix A(w h ) is required only to form a matrixvector product. This product can be approximated with the use of the Taylor series expansion, Av Φ(w h + ǫv) Φ(w h ). (6) ǫ The parameter ǫ is computed in every inner iteration as a function of a machine precision and the norm of v: ǫmach ǫ = v. (7) Thus, we can avoid the generation of the Jacobian matrix and save some memory. In the matrix version the Jacobian matrix is generated at the beginning of iterations and in each iteration there is one matrix-vector multiplication. On the other hand, in the matrix-free version one numerical flux is computed in each iteration. The computing time depends on the number of inner iterations. Since the average number of iterations is small (about 8) the matrix-free version is slightly faster. However, we need some Jacobian matrix in order to generate preconditioning matrix. But this matrix is not generated in each outer iteration, it is update only when it is needed. This Jacobian matrix can be simpler than the true Jacobian matrix. Some results of the Jacobian-free method are presented by [Blanco, Zingg, 1997] and [Knoll, Keyes, 24]. 48
3 Since the Newton method is not globally convergent, we cannot start the method from an arbitrary initial condition. From this reason we solve the time dependent problem w h t + Φ(w h ) = (8) until residuum is smaller than a given tolerance. For subsonic cases the tolerance could be 1 1, for transonic 1 3. After that we switch to the Newton iterations. This problem leads to a system of linear algebraic equations in the form (I + τa)x = τb, where A and b are the same as in the time independent problem. The matrix (I +τa) is better conditioned than the matrix A and for small time steps only few inner iterations are needed. Inverse problem The solution of the inverse problem is based on the application of a direct operator and an approximate inverse operator. The direct operator P is the solution of the Euler equations described above. The approximate inverse operator L represents an assignment of an airfoil shape to a velocity distribution on the surface of the airfoil. If we put these operators together we get an equation PL(u) = f, (9) where f denotes the required velocity distribution and u is the unknown fictitious velocity (it has no physical meaning). We transform this equation into a contractive operator equation and search for its fixed point. So we get a sequence {u k } k=,u k+1 = u k + α(f PLu k ), (1) where α (,1 is such that the sequence tends to a limit u. The first term u of the sequence is suitably chosen (for example u = f). The desired airfoil ψ is obtained by the application of L to the limit u. The velocity distribution is given on the upper and lower side of the airfoil along its mean camber line. This curve can be different from the standard mean camber line, because it connects the stagnation points on the leading and trailing edge. The position of the stagnation point on the leading edge depends on the angle of attack. If the chord line lies on the x-axis and the stagnation point on the leading edge is at the origin and the point on the trailing edge has the x-coordinate b, then the velocity distribution f = {f upper,f lower } must satisfy f upper () = f lower () = and f upper (b) = f lower (b). (11) We consider the approximate inverse operator L : u (ψ 1,ψ 2 ), where u is a fictitious velocity distribution and the (ψ 1,ψ 2 ) coordinates of the airfoil are given by the relations s (x) ψ 1 (x) = x t(x), (12) 1 + s 2 (x) ψ 2 (x) = s(x) ± t(x) 1, x,b. (13) 1 + s 2 (x) The function s(x) describes the mean camber line and the function t(x) describes the thickness of the airfoil. The derivation of these functions can be found in [Pelant, 1998]. 49
4 Since there is some restriction on the velocity distribution, it is necessary to find a suitable angle of attack to ensure that the conditions (11) are satisfied. Thus, we look for a root α of a function describing the position of the stagnation point. Each evaluation requires a steadystate solution of the Euler equations. If we remember the angle from previous inverse problem iteration, we are close enough and only a few evaluations are needed. Moreover, the change of the angle is so small that the flow remains nearly the same. Numerical experiments The described method for the solution of the inverse problem was applied to the velocity distribution along a known airfoil in order to show that our method yields this airfoil. In the first case we start from the airfoil NACA321 and get a velocity distribution satisfying the conditions for the stagnation points (the angle of attack 1.94 ). The Mach number at infinity is M =.52. In the second case we start from the airfoil Eppler E337. The Mach number at infinity is M =.5 and the angle of attack α = In both cases the parameter α in (1) is set to 1. By the symbol v we shall denote the computed velocity distribution on Mach.5 f-v Iterations Figure 1: Given velocity distribution, Figure 2: Convergence history of the L 2 -norm NACA321, M =.52, α = 1.94 f v 2 in a log-scale Y f-v Figure 3: Resulting airfoil shape (NACA321) Figure 4: Distribution of the error f v 5
5 the airfoil and by f the given velocity distribution, both along the mean camber line. Case 1 The given velocity distribution corresponding to NACA321 is shown in Figure 1. The convergence history of the error f v 2 is shown in Figure 2. After 2 iterations of the inverse problem the L 2 -norm of error f v 2 = Resulting shape shows Figure 3 and the error distribution shows Figure 4. Case 2 The given velocity distribution corresponding to E337 is shown in Figure 5. The convergence history of the error is shown in Figure 6. After 2 iterations of the inverse problem the L 2 -norm of error f v 2 = Resulting shape is shown in Figure 7 and the error distribution is shown in Figure 8..7 Mach.5 f-v Iterations Figure 5: Given velocity distribution, E337, Figure 6: Convergence history of the L 2 -norm M =.5, α = 2.48 f v 2 in a log-scale Y f-v Figure 7: Resulting airfoil shape (E337) Figure 8: Distribution of the error f v 51
6 Conclusion We have described a method, which allows us to get an airfoil shape for a given velocity distribution. It turns out that the ability of fast solving the Euler equations is essential. We need to compute at least one steady-state solution in each iteration of the inverse problem. From this reason explicit methods are unsuitable and we use an implicit method which needs the solution of a large ill-conditioned system of linear algebraic equations. GMRES with ILU preconditioning is used. Apparently, it follows from Figures 4 and 8 that the error distribution is not uniform. This is due to the behaviour of the velocity function and also the inaccuracy of the inverse operator. But the influence of the error near the leading edge is small. This method could be modified to the use of the Navier-Stokes equations or to compute an airfoil from the pressure distribution. Acknowledgments. Results have been achieved with the support of the Ministry of Education, Youth and Sports of the Czech Republic, the project MSM of Aeronautical Research and Test Institute and the project MSM of the Charles University. References Blanco M., Zingg D.W., A Fast Solver for the Euler Equations on Unstructured Grids Using a Newton- GMRES Method, AIAA Paper , Feistauer M., Felcman J., Straškraba I., Mathematical and Computational Methods for Compressible Flow, Clarendon Press, Oxford, 23. Knoll D.A., Keyes D.E., Jacobian-free Newton-Krylov Methods: A Survey of Approaches and Applications, Journal of Computational Physics 193, , 24. Pelant J., Inverse Problem for Two-dimensional Flow around a Profile, Report VZLU Z-69, Prague, Saad Y., Iterative Methods for Sparse Linear Systems, 2nd ed., SIAM,
On the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More informationA High-Order Accurate Unstructured GMRES Solver for Poisson s Equation
A High-Order Accurate Unstructured GMRES Solver for Poisson s Equation Amir Nejat * and Carl Ollivier-Gooch Department of Mechanical Engineering, The University of British Columbia, BC V6T 1Z4, Canada
More informationComparison of parallel preconditioners for a Newton-Krylov flow solver
Comparison of parallel preconditioners for a Newton-Krylov flow solver Jason E. Hicken, Michal Osusky, and David W. Zingg 1Introduction Analysis of the results from the AIAA Drag Prediction workshops (Mavriplis
More informationCase C1.3: Flow Over the NACA 0012 Airfoil: Subsonic Inviscid, Transonic Inviscid, and Subsonic Laminar Flows
Case C1.3: Flow Over the NACA 0012 Airfoil: Subsonic Inviscid, Transonic Inviscid, and Subsonic Laminar Flows Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts
More informationCase C3.1: Turbulent Flow over a Multi-Element MDA Airfoil
Case C3.1: Turbulent Flow over a Multi-Element MDA Airfoil Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts Institute of Technology I. Code Description ProjectX
More informationLimiters for Unstructured Higher-Order Accurate Solutions of the Euler Equations
Limiters for Unstructured Higher-Order Accurate Solutions of the Euler Equations Krzysztof Michalak and Carl Ollivier-Gooch Advanced Numerical Simulation Laboratory University of British Columbia Higher-order
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationNumerical Simulations of Fluid-Structure Interaction Problems using MpCCI
Numerical Simulations of Fluid-Structure Interaction Problems using MpCCI François Thirifay and Philippe Geuzaine CENAERO, Avenue Jean Mermoz 30, B-6041 Gosselies, Belgium Abstract. This paper reports
More informationEFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES
EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES B. Treidler, J.A. Ekaterineris and R.E. Childs Nielsen Engineering & Research, Inc. Mountain View, CA, 94043 Abstract Preliminary
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationCase C2.2: Turbulent, Transonic Flow over an RAE 2822 Airfoil
Case C2.2: Turbulent, Transonic Flow over an RAE 2822 Airfoil Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts Institute of Technology I. Code Description ProjectX
More informationEXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS
EXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS Edisson Sávio de Góes Maciel, edissonsavio@yahoo.com.br Mechanical
More informationPROGRESS IN NEWTON-KRYLOV METHODS FOR AERODYNAMIC CALCULATIONS. University of Toronto Institute for Aerospace Studies
PROGRESS IN NEWTON-KRYLOV METHODS FOR AERODYNAMIC CALCULATIONS Alberto Pueyo David W. Zingg y University of Toronto Institute for Aerospace Studies 4925 Duerin Street, Downsview, Ontario M3H 5T6 Canada
More informationOptimization with Gradient and Hessian Information Calculated Using Hyper-Dual Numbers
Optimization with Gradient and Hessian Information Calculated Using Hyper-Dual Numbers Jeffrey A. Fike and Juan J. Alonso Department of Aeronautics and Astronautics, Stanford University, Stanford, CA 94305,
More informationImprovements to a Newton-Krylov Adjoint Algorithm for Aerodynamic Optimization
Improvements to a Newton-Krylov Adjoint Algorithm for Aerodynamic Optimization David W. Zingg, Timothy M. Leung, Laslo Diosady, Anh H. Truong, and Samy Elias Institute for Aerospace Studies, University
More informationAn Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm
An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm P.A. Sherar, C.P. Thompson, B. Xu, B. Zhong Abstract A new method is presented to deal with shape optimization problems.
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationHigh-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method
High-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method Per-Olof Persson University of California, Berkeley, Berkeley, CA 9472-384, U.S.A. We study some of the properties
More informationEfficient Implicit Time-Marching Methods Using a Newton-Krylov Algorithm
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-164 Efficient Implicit Time-Marching Methods Using a Newton-Krylov
More informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
More informationAn efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment
The Aeronautical Journal November 2015 Volume 119 No 1221 1451 An efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment J. A. Camberos
More informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
More informationSubsonic Airfoils. W.H. Mason Configuration Aerodynamics Class
Subsonic Airfoils W.H. Mason Configuration Aerodynamics Class Typical Subsonic Methods: Panel Methods For subsonic inviscid flow, the flowfield can be found by solving an integral equation for the potential
More informationA NEW MIXED PRECONDITIONING METHOD BASED ON THE CLUSTERED ELEMENT -BY -ELEMENT PRECONDITIONERS
Contemporary Mathematics Volume 157, 1994 A NEW MIXED PRECONDITIONING METHOD BASED ON THE CLUSTERED ELEMENT -BY -ELEMENT PRECONDITIONERS T.E. Tezduyar, M. Behr, S.K. Aliabadi, S. Mittal and S.E. Ray ABSTRACT.
More informationComputation of Fictitious Gas Flow with Euler Equations
1 Computation of Fictitious Gas Flow with Euler Equations Pei Li and Helmut Sobieczky DLR Göttingen, Germany Abstract The Fictitious Gas Concept supports some computational design methods to construct
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationSHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY
SHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY A. N. Ryabinin Department of Hydroaeromechanics, Faculty of Mathematics and Mechanics, Saint-Petersburg State University, St. Petersburg, Russia E-Mail: a.ryabinin@spbu.ru
More informationCONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION USED AS A TEST CASE FOR THE INVISCID PART OF RANS SOLVERS
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More information39th AIAA Aerospace Sciences Meeting and Exhibit January 8 11, 2001/Reno, NV
AIAA 1 717 Static Aero-elastic Computation with a Coupled CFD and CSD Method J. Cai, F. Liu Department of Mechanical and Aerospace Engineering University of California, Irvine, CA 92697-3975 H.M. Tsai,
More informationAirfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition
Airfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition.5.5.5.5.5.5.5..5.95.9.85.8.75.7 Patrick A. LeGresley and Juan J. Alonso Dept. of Aeronautics & Astronautics
More informationDiscontinuous Galerkin Spectral Element Approximations for CFD
Discontinuous Galerkin Spectral Element Approimations for CFD D.A. Kopriva Florida State Universit Tallahassee, FL 3236 G.B. Jacobs San Diego State Universit San Diego, CA 92182 September 3, 211 1 Code
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationSubsonic Airfoils. W.H. Mason Configuration Aerodynamics Class
Subsonic Airfoils W.H. Mason Configuration Aerodynamics Class Most people don t realize that mankind can be divided into two great classes: those who take airfoil selection seriously, and those who don
More informationA Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids
A Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids Ruo Li, Xin Wang and Weibo Zhao January 26, 27 Abstract We proposed an efficient and robust algorithm to solve the steady Euler
More informationComparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils
Comparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils Moritz Kompenhans 1, Esteban Ferrer 2, Gonzalo Rubio, Eusebio Valero E.T.S.I.A. (School of Aeronautics)
More informationImplicit/Multigrid Algorithms for Incompressible Turbulent Flows on Unstructured Grids
AIAA 95-1740 Implicit/Multigrid Algorithms for Incompressible Turbulent Flows on Unstructured Grids W. Kyle Anderson,* Russ D. Rausch, Daryl L. Bonhaus NASA Langley Research Center, Hampton, Virginia 23681
More informationA higher-order finite volume method with collocated grid arrangement for incompressible flows
Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.
More informationGrid Dependence Study of Transonic/Supersonic Flow Past NACA Air-foil using CFD Hemanth Kotaru, B.Tech (Civil Engineering)
Grid Dependence Study of Transonic/Supersonic Flow Past NACA Air-foil using CFD Hemanth Kotaru, B.Tech (Civil Engineering) Abstract Computational fluid dynamics is a relatively young field in engineering.
More informationAil implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations
Ail implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations Vimala Dutta Computational and Theoretical Fluid Dynamics Division National Aerospace
More informationAIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS
AIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS Emre Alpman Graduate Research Assistant Aerospace Engineering Department Pennstate University University Park, PA, 6802 Abstract A new methodology
More informationHigh-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering
High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering Daniel J. Garmann and Miguel R. Visbal Air Force Research Laboratory, Wright-Patterson
More informationA Knowledge Based Approach to Mesh Optimization in CFD Domain: ID Euler Code Example
A Knowledge Based Approach to Mesh Optimization in CFD Domain: ID Euler Code Example Tharini Santhanam, J.C. Browne, J. Kallinderis and D. Miranker Department of Computer Science The University of Texas
More informationA robust WENO type finite volume solver for steady Euler equations on unstructured grids
A robust WENO type finite volume solver for steady Euler equations on unstructured grids Guanghui Hu,, Ruo Li 2, Tao Tang Department of Mathematics, Hong Kong Baptist University, Kowloon Tong, Hong Kong.
More informationMATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES
MATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES Goals of research Dr. Uldis Raitums, Dr. Kārlis Birģelis To develop and investigate mathematical properties of algorithms
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationAerodynamic Shape Optimization in Three-Dimensional Turbulent Flows Using a Newton-Krylov Approach
Aerodynamic Shape Optimization in Three-Dimensional Turbulent Flows Using a Newton-Krylov Approach Lana M. Osusky and David W. Zingg University of Toronto Institute for Aerospace Studies, 4925 Dufferin
More information41st AIAA Aerospace Sciences Meeting and Exhibit Jan 6 9, 2003/Reno, Nevada
AIAA 23 68 Aerodynamic Design of Cascades by Using an Adjoint Equation Method Shuchi Yang, Hsiao-Yuan Wu, and Feng Liu Department of Mechanical and Aerospace Engineering University of California, Irvine,
More informationIntroduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workshop 2 Transonic Flow Over a NACA 0012 Airfoil. Introduction to CFX WS2-1 Goals The purpose of this tutorial is to introduce the user to modelling flow in high speed external aerodynamic applications.
More informationAxisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
More informationApplication of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution
Application of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution P. BURDA a,, J. NOVOTNÝ b,, J. ŠÍSTE a, a Department of Mathematics Czech University of Technology
More informationNumerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1
Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1 Overview These notes are meant to be an overview of the various memory access patterns
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationSensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method
2th AIAA Computational Fluid Dynamics Conference 27-3 June 2, Honolulu, Hawaii AIAA 2-348 Sensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method Li Wang
More informationApplication of Wray-Agarwal Turbulence Model for Accurate Numerical Simulation of Flow Past a Three-Dimensional Wing-body
Washington University in St. Louis Washington University Open Scholarship Mechanical Engineering and Materials Science Independent Study Mechanical Engineering & Materials Science 4-28-2016 Application
More informationAn Implicit Hermite WENO Reconstruction-Based Discontinuous Galerkin Method on Tetrahedral Grids
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-4205 An Implicit Hermite WENO Reconstruction-Based Discontinuous Galerkin Method on
More informationAbstract. Introduction
EULER SOLUTIONS AS LIMIT OF INFINITE REYNOLDS NUMBER FOR SEPARATION FLOWS AND FLOWS WITH VORTICES Wolfgang Schmidt and Antony Jameson Dornier GmbH, D-7990 Friedrichshafen, FRG and Princeton University,
More informationAdjoint Solver Workshop
Adjoint Solver Workshop Why is an Adjoint Solver useful? Design and manufacture for better performance: e.g. airfoil, combustor, rotor blade, ducts, body shape, etc. by optimising a certain characteristic
More informationA linear solver based on algebraic multigrid and defect correction for the solution of adjoint RANS equations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2014; 74:846 855 Published online 24 January 2014 in Wiley Online Library (wileyonlinelibrary.com)..3878 A linear solver
More informationAn Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid
An Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid J. V. Lassaline Ryerson University 35 Victoria St, Toronto, ON, M5B 2K3, Canada D. W. Zingg University
More informationCOMPARISON OF SHOCK WAVE INTERACTION FOR THE THREE-DIMENSIONAL SUPERSONIC BIPLANE WITH DIFFERENT PLANAR SHAPES
26 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMPARISON OF SHOCK WAVE INTERACTION FOR THE THREE-DIMENSIONAL SUPERSONIC BIPLANE WITH DIFFERENT PLANAR SHAPES M. Yonezawa*, H. Yamashita* *Institute
More informationPost Stall Behavior of a Lifting Line Algorithm
Post Stall Behavior of a Lifting Line Algorithm Douglas Hunsaker Brigham Young University Abstract A modified lifting line algorithm is considered as a low-cost approach for calculating lift characteristics
More informationIntroduction to Aerodynamic Shape Optimization
Introduction to Aerodynamic Shape Optimization 1. Aircraft Process 2. Aircraft Methods a. Inverse Surface Methods b. Inverse Field Methods c. Numerical Optimization Methods Aircraft Process Conceptual
More informationFOR P3: A monolithic multigrid FEM solver for fluid structure interaction
FOR 493 - P3: A monolithic multigrid FEM solver for fluid structure interaction Stefan Turek 1 Jaroslav Hron 1,2 Hilmar Wobker 1 Mudassar Razzaq 1 1 Institute of Applied Mathematics, TU Dortmund, Germany
More informationproblems on. Of course, additional interface conditions must be
Thirteenth International Conference on Domain Decomposition Methods Editors: N. Debit, M.Garbey, R. Hoppe, J. Périaux, D. Keyes, Y. Kuznetsov c 2001 DDM.org 53 Schur Complement Based Preconditioners for
More informationA Hybrid Cartesian Grid and Gridless Method for Compressible Flows
rd AIAA Aerospace Sciences Meeting and Exhibit,, January 5, Reno, Nevada A Hybrid Cartesian Grid and Gridless Method for Compressible Flows Hong Luo and Joseph D. Baum Science Applications International
More informationAERODYNAMIC SHAPES DESIGN ON THE BASE OF DIRECT NEWTON TYPE OPTIMIZATION METHOD
AERODYNAMIC SHAPES DESIGN ON THE BASE OF DIRECT NEWTON TYPE OPTIMIZATION METHOD A.V. Grachev*, A.N. Kraiko**, S.A. Takovitskii* *Central Aerohydrodynamic Institute (TsAGI), **Central Institute of Aviation
More informationMultigrid Solvers in CFD. David Emerson. Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK
Multigrid Solvers in CFD David Emerson Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK david.emerson@stfc.ac.uk 1 Outline Multigrid: general comments Incompressible
More informationSolver Settings. Introductory FLUENT Training ANSYS, Inc. All rights reserved. ANSYS, Inc. Proprietary
Solver Settings Introductory FLUENT Training 2006 ANSYS, Inc. All rights reserved. 2006 ANSYS, Inc. All rights reserved. 5-2 Outline Using the Solver Setting Solver Parameters Convergence Definition Monitoring
More informationSeed Point. Agglomerated Points
AN ASSESSMENT OF LINEAR VERSUS NON-LINEAR MULTIGRID METHODS FOR UNSTRUCTURED MESH SOLVERS DIMITRI J. MAVRIPLIS Abstract. The relative performance of a non-linear FAS multigrid algorithm and an equivalent
More informationAIAA A Computational Study of Subsonic Flows over A Medium Range Cargo Aircraft
S AIAA 23-366 A Computational Study of Subsonic Flows over A Medium Range Cargo Aircraft 2st Applied Aerodynamic Conference 23-26 June 23 Orlando, Florida For permission to copy or to republish, contact
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationParallel Semi Meshless Stencil Selection for Moving Geometry Simulations
Parallel Semi Meshless Stencil Selection for Moving Geometry Simulations J. Angulo, D. J. Kennett, S. Timme and K. J. Badcock University of Liverpool, Liverpool, England L69 3BX, United Kingdom Computational
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationPrecise FEM solution of corner singularity using adjusted mesh applied to 2D flow
Precise FEM solution of corner singularity using adjusted mesh applied to 2D flow Jakub Šístek, Pavel Burda, Jaroslav Novotný Department of echnical Mathematics, Czech echnical University in Prague, Faculty
More informationSupersonic Wing Design Method Using an Inverse Problem for Practical Application
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-1465 Supersonic Wing Design Method Using an Inverse Problem for Practical
More informationAN INVERSE DESIGN METHOD FOR ENGINE NACELLES AND WINGS
24th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES AN INVERSE DESIGN METHOD FOR ENGINE NACELLES AND WINGS Roland Wilhelm German Aerospace Center DLR, Lilienthalplatz 7, D-388 Braunschweig, Germany
More informationPerformance Evaluation of Multiple and Mixed Precision Iterative Refinement Method and its Application to High-Order Implicit Runge-Kutta Method
Performance Evaluation of Multiple and Mixed Precision Iterative Refinement Method and its Application to High-Order Implicit Runge-Kutta Method Tomonori Kouya Shizuoa Institute of Science and Technology,
More informationTHE EFFECT OF REPLACING THE JOUKOWSKI MAP WITH THE GENERALIZED KARMAN-TREFFTZ MAP IN THE METHOD OF ZEDAN
GSJ: VOLUME 6, ISSUE 2, FEBRUARY 2018 1 GSJ: Volume 6, Issue 2, February 2018, Online: ISSN 2320-9186 THE EFFECT OF REPLACING THE JOUKOWSKI MAP WITH THE GENERALIZED KARMAN-TREFFTZ MAP IN THE METHOD OF
More informationA Higher-Order Accurate Unstructured Finite Volume Newton-Krylov Algorithm for Inviscid Compressible Flows
A Higher-Order Accurate Unstructured Finite Volume Newton-Krylov Algorithm for Inviscid Compressible Flows by AMIR NEJAT B.Sc. (Aerospace Engineering), AmirKabir University of Technology, 1996 M.Sc. (Aerospace
More informationTVD Flux Vector Splitting Algorithms Applied to the Solution of the Euler and Navier-Stokes Equations in Three-Dimensions Part II
TVD Flux Vector Splitting Algorithms Applied to the Solution of the Euler and Navier-Stokes Equations in Three-Dimensions Part II EDISSON SÁVIO DE GÓES MACIEL IEA- Aeronautical Engineering Division ITA
More informationA STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS
HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS Muthukumaran.C.K.
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationDevelopment of a Consistent Discrete Adjoint Solver for the SU 2 Framework
Development of a Consistent Discrete Adjoint Solver for the SU 2 Framework Tim Albring, Max Sagebaum, Nicolas Gauger Chair for Scientific Computing TU Kaiserslautern 16th Euro-AD Workshop, Jena December
More informationApplication of Jetstream to a Suite of Aerodynamic Shape Optimization Problems. Karla Telidetzki
Application of Jetstream to a Suite of Aerodynamic Shape Optimization Problems by Karla Telidetzki A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate
More informationAerofoil Optimisation Using CST Parameterisation in SU2
Aerofoil Optimisation Using CST Parameterisation in SU2 Marques, S., & Hewitt, P. (2014). Aerofoil Optimisation Using CST Parameterisation in SU2. Paper presented at Royal Aeronautical Society Applied
More informationA Parallel Newton Krylov Method for Navier Stokes Rotorcraft Codes
International Journal of Computational Fluid Dynamics, 2003 Vol. 17 (3), pp. 225 230 Short Communication A Parallel Newton Krylov Method for Navier Stokes Rotorcraft Codes KIVANC EKICI a and ANASTASIOS
More informationINTERACTIVE AERODYNAMICS ANALYSIS AND DESIGN PROGRAMS FOR USE IN THE UNDERGRADUATE ENGINEERING CURRICULUM
INTERACTIVE AERODYNAMICS ANALYSIS AND DESIGN PROGRAMS FOR USE IN THE UNDERGRADUATE ENGINEERING CURRICULUM Ralph Latham, Kurt Gramoll and L. N. Sankar School of Aerospace Engineering Georgia Institute of
More informationExpress Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil
Express Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil Dimitrios Sofialidis Technical Manager, SimTec Ltd. Mechanical Engineer, PhD PRACE Autumn School 2013 -
More informationNUMERICAL SOLUTION OF INVISCID AND VISCOUS FLOWS USING MODERN SCHEMES AND QUADRILATERAL OR TRIANGULAR MESH
Proceedings of the Czech Japanese Seminar in Applied Mathematics 24 Czech Technical University in Prague, August 4-7, 24 pp. 8 NUMERICAL SOLUTION OF INVISCID AND VISCOUS FLOWS USING MODERN SCHEMES AND
More informationKeywords: CFD, aerofoil, URANS modeling, flapping, reciprocating movement
L.I. Garipova *, A.N. Kusyumov *, G. Barakos ** * Kazan National Research Technical University n.a. A.N.Tupolev, ** School of Engineering - The University of Liverpool Keywords: CFD, aerofoil, URANS modeling,
More informationAERODYNAMIC DESIGN OF FLYING WING WITH EMPHASIS ON HIGH WING LOADING
AERODYNAMIC DESIGN OF FLYING WING WITH EMPHASIS ON HIGH WING LOADING M. Figat Warsaw University of Technology Keywords: Aerodynamic design, CFD Abstract This paper presents an aerodynamic design process
More informationA COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS
A COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS L. Mangani Maschinentechnik CC Fluidmechanik und Hydromaschinen Hochschule Luzern Technik& Architektur
More informationComparison of h- and p- Adaptations for Spectral Difference Methods
40th Fluid Dynamics Conference and Exhibit 28 June - 1 July 2010, Chicago, Illinois AIAA 2010-4435 Comparison of h- and p- Adaptations for Spectral Difference Methods Yi Li, Sachin Premasuthan and Antony
More informationA parallel direct/iterative solver based on a Schur complement approach
A parallel direct/iterative solver based on a Schur complement approach Gene around the world at CERFACS Jérémie Gaidamour LaBRI and INRIA Bordeaux - Sud-Ouest (ScAlApplix project) February 29th, 2008
More informationNUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING
Review of the Air Force Academy No.3 (35)/2017 NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING Cvetelina VELKOVA Department of Technical Mechanics, Naval Academy Nikola Vaptsarov,Varna, Bulgaria (cvetelina.velkova1985@gmail.com)
More informationA parallel Newton-Krylov-Schur flow solver for the Navier-Stokes equations using the SBP-SAT approach
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January, Orlando, Florida AIAA -6 A parallel Newton-Krylov-Schur flow solver for the Navier-Stokes equations
More informationA COMPARISON OF VARIOUS NODAL DISCONTINUOUS GALERKIN METHODS FOR THE 3D EULER EQUATIONS
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 information