Numerical Approximations of the Dynamical System Generated by Burgers' Equation with Neumann-Dirichlet Boundary Conditions
ESAIM: Mathematical Modelling and Numerical Analysis , Tome 47 (2013) no. 5, pp. 1465-1492

Voir la notice de l'article provenant de la source Numdam

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.

DOI : 10.1051/m2an/2013084
Classification : 37L05, 65P30, 35B32, 35B41, 65M99, 65P40, 34B15
Keywords: nonlinear dynamical system, finite precision arithmetic, bifurcation, asymptotic behavior, numerical approximation, stability, nonlinear partial differential equation, boundary value problem
@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 },
     pages = {1465--1492},
     publisher = {EDP-Sciences},
     volume = {47},
     number = {5},
     year = {2013},
     doi = {10.1051/m2an/2013084},
     mrnumber = {3100771},
     zbl = {1283.37072},
     language = {en},
     url = {http://geodesic.mathdoc.fr/articles/10.1051/m2an/2013084/}
}
TY  - JOUR
AU  - Allen, Edward J.
AU  - Burns, John A.
AU  - Gilliam, David S.
TI  - Numerical Approximations of the Dynamical System Generated by Burgers' Equation with Neumann-Dirichlet Boundary Conditions
JO  - ESAIM: Mathematical Modelling and Numerical Analysis 
PY  - 2013
SP  - 1465
EP  - 1492
VL  - 47
IS  - 5
PB  - EDP-Sciences
UR  - http://geodesic.mathdoc.fr/articles/10.1051/m2an/2013084/
DO  - 10.1051/m2an/2013084
LA  - en
ID  - M2AN_2013__47_5_1465_0
ER  - 
%0 Journal Article
%A Allen, Edward J.
%A Burns, John A.
%A Gilliam, David S.
%T Numerical Approximations of the Dynamical System Generated by Burgers' Equation with Neumann-Dirichlet Boundary Conditions
%J ESAIM: Mathematical Modelling and Numerical Analysis 
%D 2013
%P 1465-1492
%V 47
%N 5
%I EDP-Sciences
%U http://geodesic.mathdoc.fr/articles/10.1051/m2an/2013084/
%R 10.1051/m2an/2013084
%G en
%F 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 , Tome 47 (2013) no. 5, pp. 1465-1492. doi: 10.1051/m2an/2013084

Cité par Sources :