Quasi-convergence of an implementation of optimal balance by backward-forward nudging

Optimal balance is a non-asymptotic numerical method to compute a point on the slow manifold for certain two-scale dynamical systems. It works by solving a modified version of the system as a boundary value problem in time, where the nonlinear terms are adiabatically ramped up from zero to the fully nonlinear dynamics. A dedicated boundary value solver, however, is often not directly available. The most natural alternative is a nudging solver, where the problem is repeatedly solved forward and backward in time and the respective boundary conditions are restored whenever one of the temporal end points is visited. In this paper, we show quasi-convergence of this scheme in the sense that the termination residual of the nudging iteration is as small as the asymptotic error of the method itself, i.e., under appropriate assumptions exponentially small. This confirms that optimal balance in its nudging formulation is an effective algorithm. Further, it shows that the boundary value problem formulation of optimal balance is well posed up at most a residual error as small as the asymptotic error of the method itself. The key step in our proof is a careful two-component Gronwall inequality.

[1]  M. Oliver,et al.  Optimal balance for rotating shallow water in primitive variables , 2020 .

[2]  P. De Maesschalck,et al.  Gevrey asymptotic properties of slow manifolds , 2019, Nonlinearity.

[3]  M. Oliver,et al.  A direct construction of a slow manifold for a semilinear wave equation of Klein–Gordon type , 2019, Journal of Differential Equations.

[4]  M. Oliver,et al.  Hs-class construction of an almost invariant slow subspace for the Klein–Gordon equation in the non-relativistic limit , 2018 .

[5]  Georg A. Gottwald,et al.  Optimal Balance via Adiabatic Invariance of Approximate Slow Manifolds , 2017, Multiscale Model. Simul..

[6]  G. Gottwald,et al.  Slow dynamics via degenerate variational asymptotics , 2014, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[7]  Colin Cotter,et al.  Data assimilation on the exponentially accurate slow manifold , 2013, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[8]  K. U. Kristiansen,et al.  Exponential estimates of symplectic slow manifolds , 2012, 1208.4219.

[9]  Didier Auroux,et al.  The Back and Forth Nudging algorithm for data assimilation problems : theoretical results on transport equations , 2008, 0809.4838.

[10]  Georg A. Gottwald,et al.  Long-Time Accuracy for Approximate Slow Manifolds in a Finite-Dimensional Model of Balance , 2007, J. Nonlinear Sci..

[11]  Colin J. Cotter,et al.  Semigeostrophic Particle Motion and Exponentially Accurate Normal forms , 2005, Multiscale Model. Simul..

[12]  A. Viudez,et al.  Optimal potential vorticity balance of geophysical flows , 2004, Journal of Fluid Mechanics.

[13]  David G. Dritschel,et al.  The Contour-Advective Semi-Lagrangian Algorithm for the Shallow Water Equations , 1999 .

[14]  David G. Dritschel,et al.  A contour‐advective semi‐lagrangian numerical algorithm for simulating fine‐scale conservative dynamical fields , 1997 .

[15]  Neil Fenichel Geometric singular perturbation theory for ordinary differential equations , 1979 .

[16]  P. Davis,et al.  Linear generalizations of Gronwall’s inequality , 1976 .