Using Burgers' equation with mixed Neumann-Dirichlet boundary conditions, we highlight a problem that can arise in the numerical approximation of nonlinear dynamical systems on computers with a finite precision floating point number system. We describe the dynamical system generated by Burgers' equation with mixed boundary conditions, summarize some of its properties and analyze the equilibrium states for finite dimensional dynamical systems that are generated by numerical approximations of this system. It is important to note that there are two fundamental differences between Burgers' equation with mixed Neumann-Dirichlet boundary conditions and Burgers' equation with both Dirichlet boundary conditions. First, Burgers' equation with homogenous mixed boundary conditions on a finite interval cannot be linearized by the Cole-Hopf transformation. Thus, on finite intervals Burgers' equation with a homogenous Neumann boundary condition is truly nonlinear. Second, the nonlinear term in Burgers' equation with a homogenous Neumann boundary condition is not conservative. This structure plays a key role in understanding the complex dynamics generated by Burgers' equation with a Neumann boundary condition and how this structure impacts numerical approximations. The key point is that, regardless of the particular numerical scheme, finite precision arithmetic will always lead to numerically generated equilibrium states that do not correspond to equilibrium states of the Burgers' equation. In this paper we establish the existence and stability properties of these numerical stationary solutions and employ a bifurcation analysis to provide a detailed mathematical explanation of why numerical schemes fail to capture the correct asymptotic dynamics. We extend the results in [E. Allen, J.A. Burns, D.S. Gilliam, J. Hill and V.I. Shubov, Math. Comput. Modelling 35 (2002) 1165-1195] and prove that the effect of finite precision arithmetic persists in generating a nonzero numerical false solution to the stationary Burgers' problem. Thus, we show that the results obtained in [E. Allen, J.A. Burns, D.S. Gilliam, J. Hill and V.I. Shubov, Math. Comput. Modelling 35 (2002) 1165-1195] are not dependent on a specific time marching scheme, but are generic to all convergent numerical approximations of Burgers' equation.
@article{M2AN_2013__47_5_1465_0, author = {Allen, Edward J. and Burns, John A. and Gilliam, David S.}, title = {Numerical Approximations of the Dynamical System Generated by Burgers' Equation with Neumann-Dirichlet Boundary Conditions}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique}, volume = {47}, year = {2013}, pages = {1465-1492}, doi = {10.1051/m2an/2013084}, mrnumber = {3100771}, zbl = {1283.37072}, language = {en}, url = {http://dml.mathdoc.fr/item/M2AN_2013__47_5_1465_0} }
Allen, Edward J.; Burns, John A.; Gilliam, David S. Numerical Approximations of the Dynamical System Generated by Burgers' Equation with Neumann-Dirichlet Boundary Conditions. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 47 (2013) pp. 1465-1492. doi : 10.1051/m2an/2013084. http://gdmltest.u-ga.fr/item/M2AN_2013__47_5_1465_0/
[1] Metastability and Transients in Brain Dynamics: Problems and Rigorous Results, in Long-range Interactions, Stochasticity and Fractional Dynamics; Nonlinear Physical Science, edited by Albert C.J. Luo and Valentin Afraimovich. Springer-Verlag (2010) 133-175. | MR 3135776 | Zbl 1244.92004
, and ,[2] The impact of finite precision arithmetic and sensitivity on the numerical solution of partial differential equations. Math. Comput. Model. 35 (2002) 1165-1195. | MR 1910446 | Zbl 1060.65116
, , , and ,[3] On the use of numerical methods for analysis and control of nonlinear convective systems, in Proc. of 47th IEEE Conference on Decision and Control (2008) 197-202.
, and ,[4] Stabilized Finite Element Methods and Feedback Control for Burgers' Equation, in Proc. of the 2000 American Control Conference (2000) 2745-2749.
and ,[5] Exploratory Experimentation and Computation, Notices AMS 58 (2011) 1410-1419. | MR 2884022 | Zbl 1247.00010
and ,[6] Stationary solutions for a boundary controlled Burgers' equation. Math. Comput. Model. 33 (2001) 21-37. | MR 1812539 | Zbl 0967.93051
, and ,[7] Using Global Invariant Manifolds to Understand Metastability in the Burgers Equation With Small Viscosity. SIAM Review 53 (2011) 129-153 [Published originally SIAM J. Appl. Dyn. Syst. 8 (2009) 1043-1065]. | MR 2551255 | Zbl 1201.35044
and ,[8] Control of Turbulent Flows, in Systems Modelling and Optimization, Chapman and Hall CRC, Boca Raton, FL (1999) 3-11. | MR 1678711 | Zbl 0925.93417
, and ,[9] A PDE Sensitivity Equation Method for Optimal Aerodynamic Design. J. Comput. Phys. 136 (1997) 366-384. | MR 1474410 | Zbl 0903.76064
and ,[10] Numerical stationary solutions for a viscous Burgers' equation. J. Math. Syst. Estim. Control 8 (1998) 1-16. | MR 1651458 | Zbl 0892.35134
, , and ,[11] A control problem for Burgers' equation with bounded input/output. Nonlinear Dyn. 2 (1991) 235-262.
and ,[12] A Stabilization problem for Burgers' equation with unbounded control and observation, in Estimation and Control of Distributed Parameter Systems. Int. Ser. Numer. Math. vol. 100, edited by W. Desch, F. Rappel, K. Kunisch. Springer-Verlag (1991) 51-72. | MR 1155636 | Zbl 0767.93046
and ,[13] Optimal fixed-finite-dimensional compensator for Burgers' Equation with unbounded input/output operators. ICASE Report No. 93-19. Institute for Comput. Appl. Sci. Engrg., Hampton, VA. (1993). | MR 1247468 | Zbl 0824.93016
and ,[14] On the Long Time Behavior of Approximating Dynamical Systems, in Distributed Parameter Control, edited by F. Kappel, K. Kunisch and W. Schappacher. Springer-Verlag (2001) 73-86. | Zbl 1043.65079
and ,[15] Boundary control and stabilization for a viscous Burgers' equation. Computation and Control, Progress in Systems Control Theory, vol. 15. Birkhäuser Boston, Boston, MA (1993) 105-120. | MR 1247469 | Zbl 0824.93030
and ,[16] Convergence of trajectories for a controlled viscous Burgers' equation, Control and Estimation of Distributed Parameter Systems: Nonlinear Phenomena. Int. Ser. Numer. Math., vol. 118, edited by W. Desch, F. Rappel, K. Kunisch. Birkhäuser, Basel (1994) 61-77. | MR 1313510 | Zbl 0808.93040
, and ,[17] Steady state response to Burgers' equation with varying viscosity, in Progress in Systems and Control: Computation and Control IV, edited by K. L.Bowers and J. Lund. Birkhäuser, Basel (1995) 75-98. | MR 1349584 | Zbl 0841.35101
, , and ,[18] High gain limits of trajectories and attractors for a boundary controlled viscous Burgers' equation. J. Math. Syst. Estim. Control 6 (1996) 40. | MR 1650140 | Zbl 0856.93048
, and ,[19] Numerical stationary solutions for a viscous Burgers' equation. J. Math. Syst. Estim. Control 8 (1998) 16 (electronic). | MR 1651458 | Zbl 0892.35134
, , and ,[20] On the Global Dynamics of a Controlled Viscous Burgers' Equation. J. Dyn. Control Syst. 4 (1998) 457-519. | MR 1662924 | Zbl 0943.35080
, and ,[21] Boundary Control, Stabilization and Zero-Pole Dynamics for a Nonlinear Distributed Parameter System. Int. J. Robust Nonlinear Control 9 (1999) 737-768. | MR 1711262 | Zbl 0949.93038
, and ,[22] Asymptotic Behavior of Viscous Burgers' Equations with Neumann Boundary Conditions, Third Palestinian Mathematics Conference, Bethlehem University, West Bank. Mathematics and Mathematics Education, edited by S. Elaydi, E. S. Titi, M. Saleh, S. K. Jain and R. Abu Saris. World Scientific (2002) 1-19. | Zbl 1026.35068
and ,[23] Revisiting and extending interface penalties for multi-domain summation-by-parts operators. J. Sci. Comput. 45 (2010) 118-150. | MR 2679793 | Zbl 1203.65176
, and ,[24] Metastable patterns in solutions of ut = ϵ2uxx − f(u). Comm. Pure Appl. Math. 42 (1989) 523-576. | MR 997567 | Zbl 0685.35054
and ,[25] Numerical approximation of a metastable system. IMA J. Numer. Anal. 15 (1995) 505-521. | MR 1355635 | Zbl 0833.65068
, and ,[26] Burgers' equation: A model for all reasons, in Numerical Solutions of J. Partial Differ. Eqns., edited by J. Noye. North-Holland Publ. Co. Amsterdam (1982) 139-225. | MR 649920 | Zbl 0496.76091
,[27] Optimal Neumann Control for the 2D Steady-State Navier-Stokes equations, in New Directions in Math. Fluid Mech. The Alexander. V. Kazhikhov Memorial Volume. Advances in Mathematical Fluid Mechanics, Birkhauser, Berlin (2009) 193-222. | MR 2732011 | Zbl 1200.35214
and ,[28] G. and J. K. Hale, Slow-motion manifolds, dormant instability, and singular perturbations. J. Dyn. Differ. Eqns. 1 (1989) 75-94. | MR 1010961 | Zbl 0684.34055
,[29] Invariant manifolds and the long-time asymptotics of the navier-stokes and vorticity equations on R2. Arch. Rational Mech. Anal. 163 (2002) 209-258. | MR 1912106 | Zbl 1042.37058
and ,[30] Global stability of vortex solutions of the two-dimensional Navier-Stokes equation. Commun. Math. Phys. 255 (2005) 97-129. | MR 2123378 | Zbl 1139.35084
and ,[31] Asymptotic-Numerical Study of Supersensitivity for Generalized Burgers' Equation. SIAM J. Sci. Comput. 22 (2000) 368-385. | MR 1769499 | Zbl 0961.35139
and ,[32] Recovering High-Order Accuracy in WENO Computations of Steady-State Hyperbolic Systems. J. Sci. Comput. 28 (2006) 307-318. | MR 2272634 | Zbl 1158.76365
, and ,[33] Analysis and finite element approximation of optimal control problems for the stationary Navier-Stokes equations with distributed and Neumann controls. Math. Comput. 57 (1991) 123-151. | MR 1079020 | Zbl 0747.76063
, and ,[34] Error estimates of stochastic optimal Neumann boundary control problems. SIAM J. Numer. Anal. 49 (2011) 1532-1552. | MR 2831060 | Zbl 1243.65008
, and ,[35] Spectral Methods for Time Dependent Problems, Cambridge Monographs on Applied and Computational Mathematics, vol. 21. Cambridge University Press (2006). | Zbl 1111.65093
, and ,[36] IEEE Computer Society, IEEE Standard for Binary Floating-Point Arithmetic, IEEE Std 754-1985 (1985).
[37] Application of implicit-explicit high order Runge-Kutta methods to discontinuous-Galerkin schemes. J. Comput. Phys. 225 (2007) 1753-1781. | MR 2349202 | Zbl 1123.65097
, , and ,[38] A high order spectral volume solution to the Burgers' equation using the Hopf-Cole transformation. Int. J. Numer. Meth. Fluids (2011). Available on wileyonlinelibrary.com. DOI: 10.1002/fld.2612. | Zbl 1253.76087
and ,[39] Linear and Quasilinear Equations of Parabolic Type, Translations of the AMS, vol. 23 (1968). | Zbl 0174.15403
, and ,[40] Supersensitive Boundary Value Problems, Asymptotic and Numerical Methods for Partial Differential Equations with Critical Parameters, edited by H.G. Kaper and M. Garbey. Kluwer Publishers (1993) 215-224. | MR 1222424 | Zbl 0788.65088
and ,[41] Distributed and boundary control of the viscous Burgers' equation. Numer. Funct. Anal. Optim. 18 (1997) 143-188. | MR 1442024 | Zbl 0876.93045
, and ,[42] Dynamic Compensators for a Nonlinear Conservation Law, Ph.D. Thesis, Department of Mathematics, Virginia Polytechnic Institute and State University, Blacksburg, VA, 24061 (1993). | MR 2690186
,[43] A Numerical Study of Burgers' Equation With Robin Boundary Conditions, M.S. Thesis. Department of Mathematics, Polytechnic Institute and State University, Blacksburg, VA, 24061 (2001). | Zbl 1156.35450
,[44] Boundary procedures for the time-dependent Burgers' equation under uncertainty. Acta Math. Sci. 30 (2010) 539-550. | MR 2656553 | Zbl 1240.35591
, and ,[45] Slow motion manifolds far from the attractor in multistable reaction-diffusion equations. J. Differ. Eqns. 174 (2001) 101-132. | MR 1844526 | Zbl 0986.35048
,[46] Finite element approximations of Burgers' Equation, M.S. Thesis. Departmant of Mathematics, Polytechnic Institute and State University, Blacksburg, VA, 24061 (1995).
,[47] Dynamics of Evolutionary Equations, vol. 143. Springer-Verlag (2002). | MR 1873467 | Zbl 1254.37002
and ,[48] Exact boundary conditions for the initial value problem of convex conservation laws. J. Comput. Phys. 229 (2010) 3792-3801. | MR 2609753 | Zbl 1190.65127
,[49] Internal layers, small eigenvalues, and the sensitivity of metastable motion. SIAM J. Appl. Math. 55 (1995) 425-445. | MR 1322768 | Zbl 0818.35008
and ,[50] Stabilization of solutions of boundary value problems for a second-order parabolic equation with one space variable, Differentsial'nye Uravneniya 4 (1968) 34D45. | MR 223758 | Zbl 0232.35053
,[51] Qualitative Theory of Parabolic Equations, Part 1, VSP, Utrecht, The Netherlands (1997). | MR 1644092 | Zbl 0918.35001
, and ,