# American Institute of Mathematical Sciences

September  2018, 23(7): 2859-2877. doi: 10.3934/dcdsb.2018108

## Underlying one-step methods and nonautonomous stability of general linear methods

 1 Sandia National Laboratories, P.O. Box 5800, MS 1320, Albuquerque, NM 87185-1320, USA 2 Department of Mathematics - University of Kansas, 1460 Jayhawk Blvd, Lawrence, KS 66045-7594, USA

* Corresponding author: Andrew J. Steyer

Dedicated to the memory of Timo Eirola.

Received  April 2017 Revised  September 2017 Published  March 2018

Fund Project: This research was supported in part by NSF grant DMS-1419047.

We generalize the theory of underlying one-step methods to strictly stable general linear methods (GLMs) solving nonautonomous ordinary differential equations (ODEs) that satisfy a global Lipschitz condition. We combine this theory with the Lyapunov and Sacker-Sell spectral stability theory for one-step methods developed in [34,35,36] to analyze the stability of a strictly stable GLM solving a nonautonomous linear ODE. These results are applied to develop a stability diagnostic for the solution of nonautonomous linear ODEs by strictly stable GLMs.

Citation: Andrew J. Steyer, Erik S. Van Vleck. Underlying one-step methods and nonautonomous stability of general linear methods. Discrete & Continuous Dynamical Systems - B, 2018, 23 (7) : 2859-2877. doi: 10.3934/dcdsb.2018108
##### References:

show all references

##### References:
Left: Logarithmic plot of the 2-norm of the local truncation error of the numerical solution versus time for various values of $h$. Right: Logarithmic plot of the 2-norm of the numerical solution versus time for various values of $h$. The parameter values used were $a_1=a_2=1.2$, $b_1 = -0.14$, $b_2=-0.15$, $\beta=10.0$, $\omega = 1$ with a final time of $t_f = 40$ and the initial condition $x(0)=(1,0)^T$.
Left: Logarithmic plot of the 2-norm of the local truncation error of the numerical solution versus time for various values of $h$. Right: Logarithmic plot of the 2-norm of the numerical solution versus time for various values of $h$. The parameter values used were using $b_1 = -0.5$, $b_2=-.055$, $\beta=1.0$, $\omega = 1$, and a final time of $t_f = 100$ for various values of $a=a_1=a_2$ using the step-sizes $h=0.05$ and the initial condition $x(0)=(1,0)^T$.
Results of an experiment for the solution of (3) using BDF2, $a_1=a_2=1.2$, $b_1 = -0.14$, $b_2=-0.15$, $\beta=10.0$, $\omega = 1$, and a final time of $t_f = 40$ for various step-sizes $h$ and the initial condition $x(0)=(1,0)^T$. LTEmean is the mean local truncation error, LTEmax is the maximum local truncation error, and ${\mu_{\rm{appr}}(N_f/2,N_f/2)}$ is the value of (26) where $N_f$ is the final step of the approximation.
 $h$ LTEmean LTEmax $\mu_{\rm{appr}}(N_f/2,N_f/2)$ $7.5E-1$ $1.37E10$ $1.51E11$ $7.68E-1$ $7.5E-2$ $3.75E-3$ $9.42E-3$ $9.03E-3$ $7.5E-3$ $3.60E-7$ $6.38E-4$ $-9.70E-2$ $7.5E-4$ $1.95E-9$ $6.24E-5$ $-9.04E-2$
 $h$ LTEmean LTEmax $\mu_{\rm{appr}}(N_f/2,N_f/2)$ $7.5E-1$ $1.37E10$ $1.51E11$ $7.68E-1$ $7.5E-2$ $3.75E-3$ $9.42E-3$ $9.03E-3$ $7.5E-3$ $3.60E-7$ $6.38E-4$ $-9.70E-2$ $7.5E-4$ $1.95E-9$ $6.24E-5$ $-9.04E-2$
Results of an experiment for the solution of (3) using BDF2, using $b_1 = -0.5$, $b_2=-.055$, $\beta=1.0$, $\omega = 1$, and a final time of $t_f = 100$ for various values of $a=a_1=a_2$ using the step-sizes $h=0.05$ and the initial condition $x(0)=(1,0)^T$. LTEmean is the mean local truncation error, LTEmax is the maximum local truncation error, ${\mu_{\rm{appr}}(N_f/2,N_f/2)}$ is the value of (26) where $N_f$ is the final step of the approximation, and ${\tau_{\rm{max}}}$ is the maximum value of $\tau_n$ which denotes the quotient of the local truncation error at time-steps $n+1$ and $n$.
 $a_1=a_2=a$ LTEmean LTEmax $\mu_{\rm{appr}}(N_f/2,N_f/2)$ $\tau_{\rm{max}}$ $1.15$ $5.50E-5$ $4.38E-3$ $-2.33E-2$ $1.068$ $1.45$ $1.18E-4$ $5.02E-3$ $-1.69E-3$ $1.086$ $1.75$ $2.88E-4$ $5.70E-3$ $1.78E-2$ $1.11$ $2.05$ $7.96E-4$ $6.4E-3$ $3.64E-2$ $1.23$
 $a_1=a_2=a$ LTEmean LTEmax $\mu_{\rm{appr}}(N_f/2,N_f/2)$ $\tau_{\rm{max}}$ $1.15$ $5.50E-5$ $4.38E-3$ $-2.33E-2$ $1.068$ $1.45$ $1.18E-4$ $5.02E-3$ $-1.69E-3$ $1.086$ $1.75$ $2.88E-4$ $5.70E-3$ $1.78E-2$ $1.11$ $2.05$ $7.96E-4$ $6.4E-3$ $3.64E-2$ $1.23$
 [1] Lars Grüne. Computing Lyapunov functions using deep neural networks. Journal of Computational Dynamics, 2020  doi: 10.3934/jcd.2021006 [2] Peter Giesl, Sigurdur Hafstein. System specific triangulations for the construction of CPA Lyapunov functions. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020378 [3] Peter Giesl, Zachary Langhorne, Carlos Argáez, Sigurdur Hafstein. Computing complete Lyapunov functions for discrete-time dynamical systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 299-336. doi: 10.3934/dcdsb.2020331 [4] Waixiang Cao, Lueling Jia, Zhimin Zhang. A $C^1$ Petrov-Galerkin method and Gauss collocation method for 1D general elliptic problems and superconvergence. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 81-105. doi: 10.3934/dcdsb.2020327 [5] Gervy Marie Angeles, Gilbert Peralta. Energy method for exponential stability of coupled one-dimensional hyperbolic PDE-ODE systems. Evolution Equations & Control Theory, 2020  doi: 10.3934/eect.2020108 [6] Qiang Fu, Xin Guo, Sun Young Jeon, Eric N. Reither, Emma Zang, Kenneth C. Land. The uses and abuses of an age-period-cohort method: On the linear algebra and statistical properties of intrinsic and related estimators. Mathematical Foundations of Computing, 2020  doi: 10.3934/mfc.2021001 [7] Mehdi Bastani, Davod Khojasteh Salkuyeh. On the GSOR iteration method for image restoration. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 27-43. doi: 10.3934/naco.2020013 [8] Xiaoxiao Li, Yingjing Shi, Rui Li, Shida Cao. Energy management method for an unpowered landing. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020180 [9] Ying Liu, Yanping Chen, Yunqing Huang, Yang Wang. Two-grid method for semiconductor device problem by mixed finite element method and characteristics finite element method. Electronic Research Archive, 2021, 29 (1) : 1859-1880. doi: 10.3934/era.2020095 [10] Hong Niu, Zhijiang Feng, Qijin Xiao, Yajun Zhang. A PID control method based on optimal control strategy. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 117-126. doi: 10.3934/naco.2020019 [11] Qing-Hu Hou, Yarong Wei. Telescoping method, summation formulas, and inversion pairs. Electronic Research Archive, , () : -. doi: 10.3934/era.2021007 [12] Evan Greif, Daniel Kaplan, Robert S. Strichartz, Samuel C. Wiese. Spectrum of the Laplacian on regular polyhedra. Communications on Pure & Applied Analysis, 2021, 20 (1) : 193-214. doi: 10.3934/cpaa.2020263 [13] Qingfang Wang, Hua Yang. Solutions of nonlocal problem with critical exponent. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5591-5608. doi: 10.3934/cpaa.2020253 [14] Li-Bin Liu, Ying Liang, Jian Zhang, Xiaobing Bao. A robust adaptive grid method for singularly perturbed Burger-Huxley equations. Electronic Research Archive, 2020, 28 (4) : 1439-1457. doi: 10.3934/era.2020076 [15] Zexuan Liu, Zhiyuan Sun, Jerry Zhijian Yang. A numerical study of superconvergence of the discontinuous Galerkin method by patch reconstruction. Electronic Research Archive, 2020, 28 (4) : 1487-1501. doi: 10.3934/era.2020078 [16] Yuxia Guo, Shaolong Peng. A direct method of moving planes for fully nonlinear nonlocal operators and applications. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020462 [17] Noah Stevenson, Ian Tice. A truncated real interpolation method and characterizations of screened Sobolev spaces. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5509-5566. doi: 10.3934/cpaa.2020250 [18] Yue Feng, Yujie Liu, Ruishu Wang, Shangyou Zhang. A conforming discontinuous Galerkin finite element method on rectangular partitions. Electronic Research Archive, , () : -. doi: 10.3934/era.2020120 [19] Gang Luo, Qingzhi Yang. The point-wise convergence of shifted symmetric higher order power method. Journal of Industrial & Management Optimization, 2021, 17 (1) : 357-368. doi: 10.3934/jimo.2019115 [20] Xiu Ye, Shangyou Zhang, Peng Zhu. A weak Galerkin finite element method for nonlinear conservation laws. Electronic Research Archive, 2021, 29 (1) : 1897-1923. doi: 10.3934/era.2020097

2019 Impact Factor: 1.27