An Hp-Adaptive Minimum Action Method Based on a Posteriori Error Estimate

In this work, we develop an hp-adaptivity strategy for the minimum action method (MAM) using a posteriori error estimate. MAM plays an important role in minimizing the Freidlin-Wentzell action functional, which is the central object of the Freidlin-Wentzell theory of large deviations for noise-induced transitions in stochastic dynamical systems. Because of the demanding computation cost, especially in spatially extended systems, numerical efficiency is a critical issue for MAM. Difficulties come from both temporal and spatial discretizations. One severe hurdle for the application of MAM to large scale systems is the global reparametrization in time direction, which is needed in most versions of MAM to achieve accuracy. We recently introduced a new version of MAM in [22], called tMAM, where we used some simple heuristic criteria to demonstrate that tMAM can be effectively coupled with h-adaptivity, i.e., the global reparametrization can be removed. The target of this paper is to integrate hpadaptivity into tMAM using a posteriori error estimation techniques, which provides a general adaptive MAM more suitable for parallel computing. More specifically, we use the zero-Hamiltonian constraint to define an indicator to measure the error induced by linear time scaling, and the derivative recovery technique to construct an error indicator and a regularity indicator for the transition paths approximated by finite elements. Strategies for hp-adaptivity have been developed. Numerical results are presented. AMS subject classifications: 60H35, 65C20, 65N20, 65N30

[1]  M. Freidlin,et al.  Random Perturbations of Dynamical Systems , 1984 .

[2]  J. Zabczyk,et al.  Stochastic Equations in Infinite Dimensions , 2008 .

[3]  Ira B. Schwartz,et al.  An iterative action minimizing method for computing optimal paths in stochastic dynamical systems , 2012, 1210.5153.

[4]  Jinchao Xu,et al.  Asymptotically Exact A Posteriori Error Estimators, Part II: General Unstructured Grids , 2003, SIAM J. Numer. Anal..

[5]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[6]  Guang Lin,et al.  Hybrid parallel computing of minimum action method , 2013, Parallel Comput..

[7]  Weiqing Ren,et al.  Adaptive minimum action method for the study of rare events. , 2008, The Journal of chemical physics.

[8]  Stephen J. Wright,et al.  Numerical Optimization (Springer Series in Operations Research and Financial Engineering) , 2000 .

[9]  Randolph E. Bank,et al.  hp Adaptive finite elements based on derivative recovery and superconvergence , 2011, Comput. Vis. Sci..

[10]  E. Weinan,et al.  STUDY OF NOISE-INDUCED TRANSITIONS IN THE LORENZ SYSTEM USING THE MINIMUM ACTION METHOD , 2010 .

[11]  Xiaoliang Wan,et al.  A minimum action method for small random perturbations of two-dimensional parallel shear flows , 2013, J. Comput. Phys..

[12]  E Weinan,et al.  Model the nonlinear instability of wall-bounded shear flows as a rare event: a study on two-dimensional Poiseuille flow , 2015 .

[13]  Haijun Yu,et al.  A dynamic-solver-consistent minimum action method: With an application to 2D Navier-Stokes equations , 2017, J. Comput. Phys..

[14]  M. Freidlin.  Random perturbations of reaction-difiusion equations: the quasi de-terministic approximation , 1988 .

[15]  E. Vanden-Eijnden,et al.  The geometric minimum action method: A least action principle on the space of curves , 2008 .

[16]  Xiaoliang Wan,et al.  An adaptive high-order minimum action method , 2011, J. Comput. Phys..

[17]  Jinchao Xu,et al.  Superconvergent Derivative Recovery for Lagrange Triangular Elements of Degree p on Unstructured Grids , 2007, SIAM J. Numer. Anal..

[18]  Xiaoliang Wan,et al.  A Minimum Action Method with Optimal Linear Time Scaling , 2015 .

[19]  Giovanna Jona-Lasinio,et al.  Large fluctution for a non linear heat equation with noise , 1982 .

[20]  E. Vanden-Eijnden,et al.  Rare Event Simulation of Small Noise Diffusions , 2012 .

[21]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.