An Efficient Scheme for Approximating Long-time Dynamics of a Class of Non-linear Models

Abstract

We propose a novel, highly efficient, second-order accurate, longtime unconditionally stable numerical scheme for a class of finitedimensional nonlinear models important in geophysical fluid dynamics. The scheme is highly efficient in that only a fixed symmetric positive definite linear problem (with varying right-hand sides) is solved at each time step. The solutions remain uniformly bounded for all time. We show that the scheme accurately captures the long-time dynamics of the underlying geophysical model, with the global attractors and invariant measures of the scheme converging to those of the original model as the step size approaches zero. In our numerical experiments, we adopt an indirect approach, using statistics from long-time simulations to approximate the invariant measures. Our results suggest that the convergence rate of the long-term statistics, as a function of terminal time, is approximately first-order under the Jensen-Shannon metric and half-order under the total variation metric. This implies that extremely long simulations are required to achieve high-precision approximations of the invariant measure (or climate). Nevertheless, the second-order scheme significantly outperforms its first-order counterpart, requiring far less time to reach a small neighborhood of statistical equilibrium for a given step size.

Department(s)

Mathematics and Statistics

Keywords and Phrases

global attractors; invariant measures; Jensen-Shannon entropy; Long-time behavior; Lorenz 96

International Standard Serial Number (ISSN)

2163-4548; 1526-7555

Document Type

Article - Journal

Document Version

Citation

File Type

text

Language(s)

English

Rights

© 2025 International Press, All rights reserved.

Publication Date

01 Jan 2025

Link to Full Text

Share

 
COinS