PNNL

Abstract

In this work, we develop an hp-adaptive minimum action method (MAM). MAM plays an important role to minimize the Freidlin-Wentzell action functional, which is the central object of the Freidlin-Wentzell theory of large deviations for transitions induced by small random perturbations. Because of the demanding computation cost, especially in spatially extended systems, numerical efficiency is a critical issue for MAM. Difficulties come from discretizations of both time and space. One hurdle for the application of MAM to large scale systems is the global reparametrization in time direction, which is required by most versions of MAM. We recently introduced a new version of MAM in [25], called tMAM, where we used some simple heuristic criteria to demonstrate that tMAM can be effectively coupled with h-adaptivity. The target of this paper is to integrate hp-adaptivity into tMAM using a posterior error estimation techniques. More specifically, we use the arc length constraint given by geometric minimum action method (gMAM) to define an indicator of the effect of linear time scaling in tMAM, and the derivative recovery technique to construct an error indicator and a regularity indicator for hp-refinement. Numerical results are presented.