Effects of the Jacobian evaluation on Newton's solution of the Euler equations

2005-12-01
Onur, O.
Eyi, Sinan
A Newton's method is developed for solving the 2-D Euler equations. The Euler equations are discretized using a finite-volume method with upwind flux splitting schemes. Both analytical and numerical methods are used for Jacobian calculations. Although the numerical method has the advantage of keeping the Jacobian consistent with the numerical residual vector and avoiding extremely complex analytical differentiations, it may have accuracy problems and need longer execution time. In order to improve the accuracy of numerical Jacobians, detailed error analyses are performed. Results show that the finite-difference perturbation magnitude and computer precision are the most important parameters that affect the accuracy of numerical Jacobians. A method is developed for calculating an optimal perturbation magnitude that can minimize the error in numerical Jacobians. The accuracy of the numerical Jacobians is improved significantly by using the optimal perturbation magnitude. The effects of the accuracy of numerical Jacobians on the convergence of the flow solver are also investigated. In order to reduce the execution time for numerical Jacobian evaluation, flux vectors with perturbed flow variables are calculated only for neighbouring cells. A sparse matrix solver that is based on LU factorization is used. Effects of different flux splitting methods and higher-order discretizations on the performance of the solver are analyzed. © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.
17th AIAA Computational Fluid Dynamics Conference (2005)

Suggestions

Effects of the Jacobian evaluation on Newton's solution of the Euler equations
Onur, O; Eyi, Sinan (Wiley, 2005-09-20)
Newton's method is developed for solving the 2-D Euler equations. The Euler equations are discretized using a finite-volume method with upwind flux splitting schemes. Both analytical and numerical methods are used for Jacobian calculations. Although the numerical method has the advantage of keeping the Jacobian consistent with the numerical residual vector and avoiding extremely complex analytical differentiations, it may have accuracy problems and need longer execution time. In order to improve the accurac...
Effect of the Jacobian evalution on direct solutions of the Euler equations
Onur, Ömer; Eyi, Sinan; Department of Aerospace Engineering (2003)
A direct method is developed for solving the 2-D planar/axisymmetric Euler equations. The Euler equations are discretized using a finite-volume method with upwind flux splitting schemes, and the resulting nonlinear system of equations are solved using Newton̕s Method. Both analytical and numerical methods are used for Jacobian calculations. Numerical method has the advantage of keeping the Jacobian consistent with the numerical flux vector without extremely complex or impractical analytical differentiations...
A quadtree-based adaptively-refined cartesian-grid algorithm for solution of the euler equations
Bulkök, Murat; Aksel, Mehmet Haluk; Department of Mechanical Engineering (2005)
A Cartesian method for solution of the steady two-dimensional Euler equations is produced. Dynamic data structures are used and both geometric and solution-based adaptations are applied. Solution adaptation is achieved through solution-based gradient information. The finite volume method is used with cell-centered approach. The solution is converged to a steady state by means of an approximate Riemann solver. Local time step is used for convergence acceleration. A multistage time stepping scheme is used to ...
The finite element method over a simple stabilizing grid applied to fluid flow problems
Aydın, Selçuk Han; Tezer-Sezgin, Münevver; Department of Scientific Computing (2008)
We consider the stabilized finite element method for solving the incompressible Navier-Stokes equations and the magnetohydrodynamic (MHD) equations in two dimensions. The well-known instabilities arising from the application of standard Galerkin finite element method are eliminated by using the stabilizing subgrid method (SSM), the streamline upwind Petrov-Galerkin (SUPG) method, and the two-level finite element method (TLFEM). The domain is discretized into a set of regular triangular elements. In SSM, the...
Development of discontinuous galerkin method 2 dimensional flow solver
Güngör, Osman; Özgen, Serkan; Department of Aerospace Engineering (2019)
In this work, 2 dimensional flow solutions of Euler equations are presented from the developed discontinuous Galerkin method finite element method (DGFEM) solver on unstructured grids. Euler equations govern the inviscid and adiabatic flows with a set of hyperbolic equations. The discretization of governing equations for DGFEM is given in detail. The DGFEM discretization provides high order solutions on an element-compact stencil hence only elements having common boundary are coupled. The required elementwise ...
Citation Formats
O. Onur and S. Eyi, “Effects of the Jacobian evaluation on Newton’s solution of the Euler equations,” Toronto, Ontario Canada, 2005, Accessed: 00, 2021. [Online]. Available: https://www.scopus.com/inward/record.uri?partnerID=HzOxMe3b&scp=84884938076&origin=inward.