Effect of the Jacobian evalution on direct solutions of the Euler equations

Download
2003
Onur, Ömer
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. However, numerical method may have accuracy problem and may need longer execution time. In order to improve the accuracy of numerical method detailed error analyses were performed. It was demonstrated that the finite-difference perturbation magnitude and computer precision are the most important parameters that affect the accuracy of numerical Jacobians. A relation was developed for optimum perturbation magnitude that can minimize the error in numerical Jacobians. Results show that very accurate numerical Jacobians can be calculated with optimum perturbation magnitude. The effects of the accuracy of numerical Jacobians on the convergence of flow solver are also investigated. In order to reduce the execution time for numerical Jacobian evaluation, flux vectors with perturbed flow variables are calculated for only related cells. A sparse matrix solver based on LU factorization is used for the solution, and to improve the Jacobian matrix solution some strategies are considered. Effects of different flux splitting methods, higher-order discretizations and several parameters on the performance of the solver are analyzed.

Suggestions

Effects of the Jacobian evaluation on Newton's solution of the Euler equations
Onur, O.; Eyi, Sinan (null; 2005-12-01)
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 accur...
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 cyclic thermal loading on the performance of steel H-piles in integral bridges with stub-abutments
Dicleli, Murat (Elsevier BV, 2004-02-01)
In this paper, analytical equations are developed to estimate the lateral displacement capacity of steel-H piles in integral bridges with stub abutments subjected to cyclic thermal variations. First, steel H piles that are capable of sustaining large plastic deformations are identified based on their local buckling strength. The normalized moment-curvature relationships of these piles are then obtained for various axial load levels. Next, a low-cycle fatigue damage model is employed to determine the maximum...
Differential Quadrature Solution of Hyperbolic Telegraph Equation
Pekmen, B.; Tezer, Münevver (2012-01-01)
Differential quadrature method (DQM) is proposed for the numerical solution of one- and two-space dimensional hyperbolic telegraph equation subject to appropriate initial and boundary conditions. Both polynomial-based differential quadrature (PDQ) and Fourier-based differential quadrature (FDQ) are used in space directions while PDQ is made use of in time direction. Numerical solution is obtained by using Gauss-Chebyshev-Lobatto grid points in space intervals and equally spaced and/or GCL grid points for th...
Two-dimensional unsteady Navier-Stokes solution method with moving overset grids
Tuncer, İsmail Hakkı (American Institute of Aeronautics and Astronautics (AIAA), 1997-03-01)
A simple numerical algorithm to localize intergrid boundary points and to interpolate unsteady solution variables across two-dimensional, structured overset grids is presented. Overset grids are allowed to move in time relative to each other. Intergrid boundary points are localized in a triangular stencil on the donor grid by a directional search algorithm. The final parameters of the search algorithm give the interpolation weights at the intergrid boundary point. Numerical results are presented for steady ...
Citation Formats
Ö. Onur, “Effect of the Jacobian evalution on direct solutions of the Euler equations,” M.S. - Master of Science, Middle East Technical University, 2003.