Energy-conservative implicit integration schemes for the Fokker-Planck transport equation in multidimensional geometries require inverting a dense, non-symmetric matrix (Jacobian), which is very expensive to store and solve using standard solvers. However, these limitations can be overcome with Newt
Using exact Jacobians in an implicit Newton–Krylov method
✍ Scribed by F.D. Bramkamp; H.M. Bücker; A. Rasch
- Publisher
- Elsevier Science
- Year
- 2006
- Tongue
- English
- Weight
- 368 KB
- Volume
- 35
- Category
- Article
- ISSN
- 0045-7930
No coin nor oath required. For personal study only.
✦ Synopsis
In an implicit Newton-Krylov method for inviscid, compressible fluid flow, the derivation of the analytic flux Jacobian can become quite complicated depending on the complexity of the numerical flux calculation. Practically, the derivation of the exact Jacobian by hand is an unrealistic option because of the enormous man-hour investment needed. In this work, automatic differentiation is used to evaluate the exact Jacobian of upwind schemes implemented in the flow solver QUADFLOW. We compare the use of exact Jacobians and Jacobians numerically approximated by first-order forward differences. For a two-dimensional airfoil under three different flight conditions (quasi-incompressible flow, compressible subsonic flow, and transonic flow), we show that the robustness and performance of the present finite volume scheme is significantly improved by using exact Jacobians.
📜 SIMILAR VOLUMES