Efficient and generic algorithm for rigorous integration forward in time of dPDEs. I. (English) Zbl 1296.65138

Summary: We propose an efficient and generic algorithm for rigorous integration forward in time of partial differential equations written in the Fourier basis. By rigorous integration we mean a procedure which operates on sets and return sets which are guaranteed to contain the exact solution. The presented algorithm generates, in an efficient way, normalized derivatives which are used by the Lohner algorithm to produce a rigorous bound. The algorithm has been successfully tested on several partial differential equations (PDEs) including the Burgers equation, the Kuramoto-Sivashinsky equation, and the Swift-Hohenberg equation. The problem of rigorous integration in time of partial differential equations is a problem of large computational complexity and efficient algorithms are required to deal with PDEs on higher dimensional domains, like the Navier-Stokes equation. Technicalities regarding the various optimization techniques implemented in the software used in this paper will be reported elsewhere.


65M99 Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems
35B40 Asymptotic behavior of solutions to PDEs
35B41 Attractors
65G20 Algorithms with automatic result verification
65M70 Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs
65Y20 Complexity and performance of numerical algorithms
65T50 Numerical methods for discrete and fast Fourier transforms


FADBAD++; CAPD; Taylor
Full Text: DOI


[1] Arioli, G; Koch, H, Integration of dissipative partial differential equations: a case study, SIAM J. Appl. Dyn. Syst., 9, 1119-1133, (2010) · Zbl 1298.37071
[2] Bendtsen, C., Stauning, O.: FADBAD, a Flexible C++ Package for Automatic Differentiation. Technical University of Denmark, Department of Mathematical Modelling, Copenhagen (1996)
[3] Breuer, B; McKenna, PJ; Plum, M, Multiple solutions for a semilinear boundary value problem: a computational multiplicity proof, J. Differ. Equ., 195, 243-269, (2003) · Zbl 1156.35359
[4] CAPD—computer assisted proofs in dynamics, a package for rigorous numeric, http://capd.ii.uj.edu.pl · Zbl 0513.65092
[5] Canuto, C., Hussaini, M.Y., Quarteroni, A., Zang, T.A.: Spectral Methods. Fundamentals in Single Domains. Springer, Berlin (2006) · Zbl 1093.76002
[6] Cooley, JW; Tukey, JW, An algorithm for the machine calculation of complex Fourier series, Math. Comp., 19, 297-301, (1965) · Zbl 0127.09002
[7] Cyranka, J.: Existence of Globally Attracting Fixed Points of Viscous Burgers Equation with Constant Forcing. A Computer Assisted Proof, preprint, 2011, available on-line http://www.ii.uj.edu.pl/ cyranka · Zbl 1365.65220
[8] Cyranka, J.: Efficient Algorithms for Rigorous Integration Forward in Time of dPDEs. Existence of Globally Attracting Fixed Points of Viscous Burgers Equation with Constant Forcing, a Computer Assisted Proof, PhD dissertation, Jagiellonian University, Cracow, 2013, available on-line http://www.ii.uj.edu.pl/ cyranka · Zbl 1106.65315
[9] Day, S; Hiraoka, Y; Mischaikow, K; Ogawa, T, Rigorous numerics for global dynamics: a study of Swift-Hohenberg equation, SIAM J. Appl. Dyn. Syst., 4, 1-31, (2005) · Zbl 1058.35050
[10] Day, S; Lessard, JP; Mischaikow, K, Validated continuation for equilibria of pdes, SIAM J. Numer. Anal., 45, 1398-1424, (2007) · Zbl 1151.65074
[11] Gameiro, M; Lessard, JP; Mischaikow, K, Validated continuation over large parameter ranges for equilibria of pdes, Math. Comput. Simul., 79, 1368-1382, (2008) · Zbl 1166.65379
[12] Gameiro, M; Lessard, JP, Analytic estimates and rigorous continuation for equilibria of higher-dimensional pdes, J. Differ. Equ., 249, 2237-2268, (2010) · Zbl 1256.35196
[13] Gameiro, M; Lessard, JP, Rigorous computation of smooth branches of equilibria for the three dimensional Cahn-Hilliard equation, Numer. Math., 117, 753-778, (2011) · Zbl 1216.65145
[14] Griewank, A.: Evaluating Derivatives. Principles and Techniques of Algorithmic Differentiation, Frontiers in Applied Mathematics, vol. 19. SIAM, Philadelphia (2000) · Zbl 0958.65028
[15] Heywood, JG; Nagata, W; Xie, W, A numerically based existence theorem for the Navier-Stokes equations, J. Math. Fluid. Mech., 1, 5-23, (1999) · Zbl 0934.35115
[16] Hiraoka, Y; Ogawa, T, Rigorous numerics for localized patterns to the quintic Swift-Hohenberg equation, Japan, J. Ind. Appl. Math., 22, 57-75, (2005) · Zbl 1067.65146
[17] Hiraoka, Y; Ogawa, T, An efficient estimate based on FFT in topological verification method, J. Comput. Appl. Math., 199, 238-244, (2007) · Zbl 1109.65110
[18] Jorba, Á; Zou, M, A software package for the numerical integration of ODEs by means of high-order Taylor methods, Exp. Math., 14, 99-117, (2005) · Zbl 1108.65072
[19] Kapela, T; Zgliczyński, P, A lohner-type algorithm for control systems and ordinary differential inclusions, Discret. Cont. Dyn. Sys. B, 11, 365-385, (2009) · Zbl 1185.65079
[20] Kim, M; Nakao, MT; Watanabe, Y; Nishida, T, A numerical verification method of bifurcating solutions for 3-dimensional Rayleigh-Bénard problems, Numer. Math., 111, 389-406, (2009) · Zbl 1155.76024
[21] Lohner, RJ; Cash, JR (ed.); Gladwell, I (ed.), Computation of guaranteed enclosures for the solutions of ordinary initial and boundary value problems, (1992), Oxford · Zbl 0767.65069
[22] Maier-Paape, S; Miller, U; Mischaikow, K; Wanner, T, Rigorous numerics for the Cahn-Hilliard equation on the unit square, Rev. Mat. Complut., 21, 351-426, (2008) · Zbl 1160.37418
[23] Moore, R.E.: Interval Analysis. Prentice-Hall, New Jersey (1966) · Zbl 0176.13301
[24] Moore, R.E.: Methods and Applications of Interval Analysis. SIAM, Philadelphia (1979) · Zbl 0417.65022
[25] Nakao, MT, Numerical verification methods for solutions of ordinary and partial differential equations, Numer. Funct. Anal. Optim., 22, 321-356, (2001) · Zbl 1106.65315
[26] Orszag, SA; Patterson, GS, Spectral calculations of isotropic turbulence: efficient removal of aliasing interactions, Phys. Fluids, 14, 2538-2541, (1971) · Zbl 0225.76033
[27] Simó, C.: Taylor Method for the Integration of ODE. Lectures given at the LTI07 Advanced School on Long Time Integrations, available online http://www.maia.ub.es/dsg/2007/ · Zbl 0513.65094
[28] Software package and the data from tests, http://www.ii.uj.edu.pl/ cyranka/FFT
[29] Temperton, C, Self-sorting mixed-radix fast Fourier transforms, J. Comput. Phys., 52, 1-23, (1983) · Zbl 0513.65092
[30] Temperton, C, Fast mixed-radix real Fourier transforms, J. Comput. Phys., 52, 340-350, (1983) · Zbl 0513.65094
[31] Warmus, M.: Calculus of Approximations. Bull. de l’Académie Polonaise des Sciences Cl. III 4, 253-259 (1956)
[32] Zgliczyński, P, Attracting fixed points for the Kuramoto-Sivashinsky equation—a computer assisted proof, SIAM J. Appl. Dyn. Syst., 1, 215-288, (2002) · Zbl 1004.35017
[33] Zgliczyński, P, \(C^1\)-lohner algorithm, Found. Comput. Math., 2, 429-465, (2002) · Zbl 1049.65038
[34] Zgliczyński, P, Rigorous numerics for dissipative partial differential equations II. periodic orbit for the Kuramoto-Sivashinsky PDE—a computer assisted proof, Found. Comput. Math., 4, 157-185, (2004) · Zbl 1066.65105
[35] Zgliczyński, P, Rigorous numerics for dissipative PDEs III. an effective algorithm for rigorous integration of dissipative pdes, Topol. Methods Nonlinear Anal., 36, 197-262, (2010) · Zbl 1230.65113
[36] Zgliczyński, P; Mischaikow, K, Rigorous numerics for partial differential equations: the Kuramoto-Sivashinsky equation, Found. Comput. Math., 1, 255-288, (2001) · Zbl 0984.65101
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.