arrow
Return

A monolithic homotopy continuation algorithm with application to computational fluid dynamics

delete2016-09-01
delete21
PRE
AI
D
David A. Brown *
D
David W. Zingg
DOI:10.1016/j.jcp.2016.05.031delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
A new class of homotopy continuation methods is developed suitable for globalizing quasi-Newton methods for large sparse nonlinear systems of equations. The new continuation methods, described as monolithic homotopy continuation, differ from the classical predictor-corrector algorithm in that the predictor and corrector phases are replaced with a single phase which includes both a predictor and corrector component. Conditional convergence and stability are proved analytically. Using a Laplacian-like operator to construct the homotopy, the new algorithm is shown to be more efficient than the predictor-corrector homotopy continuation algorithm as well as an implementation of the widely-used pseudo-transient continuation algorithm for some inviscid and turbulent, subsonic and transonic external aerodynamic flows over the ONERA M6 wing and the NACA 0012 airfoil using a parallel implicit Newton-Krylov finite-difference flow solver. (C) 2016 Elsevier Inc. All rights reserved.
Keywords:
Homotopy
Continuation
Globalization
Newton-Krylov
Computational fluid dynamics

Journal

Journal of Computational Physics cover
Journal of Computational Physics
IF:
3.8
Papers:
1.5W
Citations:
7.4W

Organization

U
university of toronto
Scholars:
14.7W
Papers: 12.0W
Citations: 165