A parallel block multi-level preconditioner for the 3D incompressible Navier--Stokes equations
Journal of Computational Physics
Jacobian-free Newton-Krylov methods: a survey of approaches and applications
Journal of Computational Physics
Faster PDE-based simulations using robust composite linear solvers
Future Generation Computer Systems - Special issue: Selected numerical algorithms
Adaptive Software for Scientific Computing: Co-Managing Quality-Performance-Power Tradeoffs
IPDPS '05 Proceedings of the 19th IEEE International Parallel and Distributed Processing Symposium (IPDPS'05) - Workshop 10 - Volume 11
Journal of Computational Physics
Making automatic differentiation truly automatic: coupling PETSc with ADIC
Future Generation Computer Systems
A component infrastructure for performance and power modeling of parallel scientific applications
Proceedings of the 2008 compFrame/HPC-GECO workshop on Component based high performance
Proceedings of the 2009 Workshop on Component-Based High Performance Computing
An automated component-based performance experiment environment
Proceedings of the 2009 Workshop on Component-Based High Performance Computing
Making automatic differentiation truly automatic: coupling PETSc with ADIC
Future Generation Computer Systems
Journal of Computational Physics
Characterizing the performance and energy attributes of scientific simulations
ICCS'06 Proceedings of the 6th international conference on Computational Science - Volume Part I
Hi-index | 0.02 |
Pseudotransient continuation is a practical technique for globalizing the computation of steady-state solutions of nonlinear differential equations. The technique employs adaptive time-stepping to integrate an initial value problem derived from an underlying ODE or PDE boundary value problem until sufficient accuracy in the desired steady-state root is achieved to switch over to Newton's method and gain a rapid asymptotic convergence. The existing theory for pseudotransient continuation includes a global convergence result for differential equations written in semidiscretized method-of-lines form. However, many problems are better formulated or can only sensibly be formulated as differential-algebraic equations (DAEs). These include systems in which some of the equations represent algebraic constraints, perhaps arising from the spatial discretization of a PDE constraint. Multirate systems, in particular, are often formulated as differential-algebraic systems to suppress fast time scales (acoustics, gravity waves, Alfven waves, near equilibrium chemical oscillations, etc.) that are irrelevant on the dynamical time scales of interest. In this paper we present a global convergence result for pseudotransient continuation applied to DAEs of index 1, and we illustrate it with numerical experiments on model incompressible flow and reacting flow problems, in which a constraint is employed to step over acoustic waves.