Burgers' equation or Bateman–Burgers equation is a fundamental partial differential equation and convection–diffusion equation^[1] occurring in various areas of applied mathematics, such as fluid mechanics,^[2] nonlinear acoustics,^[3] gas dynamics, and traffic flow.^[4] The equation was first introduced by Harry Bateman in 1915^[5][6] and later studied by Johannes Martinus Burgers in 1948.^[7]

For a given field ${\displaystyle u(x,t)}$ and diffusion coefficient (or kinematic viscosity, as in the original fluid mechanical context) ${\displaystyle \nu }$, the general form of Burgers' equation (also known as viscous Burgers' equation) in one space dimension is the dissipative system:

$${\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=\nu {\frac {\partial ^{2}u}{\partial x^{2}}}.}$$

When the diffusion term is absent (i.e. ${\displaystyle \nu =0}$), Burgers' equation becomes the inviscid Burgers' equation:

$${\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=0,}$$

which is a prototype for conservation equations that can develop discontinuities (shock waves). The previous equation is the advective form of the Burgers' equation. The conservative form is found to be more useful in numerical integration

$${\displaystyle {\frac {\partial u}{\partial t}}+{\frac {1}{2}}{\frac {\partial (u^{2})}{\partial x}}=0.}$$

Terms

There are 4 parameters in Burgers' equation: ${\displaystyle u,x,t}$ and ${\displaystyle \nu }$. In a system consisting of a moving viscous fluid with one spatial (${\displaystyle x}$) and one temporal (${\displaystyle t}$) dimension, e.g. a thin ideal pipe with fluid running through it, Burgers' equation describes the speed of the fluid at each location along the pipe as time progresses. The terms of the equation represent the following quantities:^[8]

• ${\displaystyle x}$: spatial coordinate
• ${\displaystyle t}$: temporal coordinate
• ${\displaystyle u(x,t)}$: speed of fluid at the indicated spatial and temporal coordinates
• ${\displaystyle \nu }$: viscosity of fluid

The viscosity is a constant physical property of the fluid, and the other parameters represent the dynamics contingent on that viscosity.

Inviscid Burgers' equation

This is a numerical simulation of the inviscid Burgers Equation in two space variables up until the time of shock formation.

The inviscid Burgers' equation is a conservation equation, more generally a first order quasilinear hyperbolic equation. The solution to the equation and along with the initial condition

$${\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=0,\quad u(x,0)=f(x)}$$

can be constructed by the method of characteristics. The characteristic equations are

$${\displaystyle {\frac {dx}{dt}}=u,\quad {\frac {du}{dt}}=0.}$$

Integration of the second equation tells us that ${\displaystyle u}$ is constant along the characteristic and integration of the first equation shows that the characteristics are straight lines, i.e.,

$${\displaystyle u=c,\quad x=ut+\xi }$$

where ${\displaystyle \xi }$ is the point (or parameter) on the x-axis (t = 0) of the x-t plane from which the characteristic curve is drawn. Since ${\displaystyle u}$ at ${\displaystyle x}$-axis is known from the initial condition and the fact that ${\displaystyle u}$ is unchanged as we move along the characteristic emanating from each point ${\displaystyle x=\xi }$, we write ${\displaystyle u=c=f(\xi )}$ on each characteristic. Therefore, the family of trajectories of characteristics parametrized by ${\displaystyle \xi }$ is

$${\displaystyle x=f(\xi )t+\xi .}$$

Thus, the solution is given by

$${\displaystyle u(x,t)=f(\xi )=f(x-ut),\quad \xi =x-f(\xi )t.}$$

This is an implicit relation that determines the solution of the inviscid Burgers' equation provided characteristics don't intersect. If the characteristics do intersect, then a classical solution to the PDE does not exist and leads to the formation of a shock wave. Whether characteristics can intersect or not depends on the initial condition. In fact, the breaking time before a shock wave can be formed is given by^[9]

$${\displaystyle t_{b}=\inf _{x}\left({\frac {-1}{f'(x)}}\right)}$$

Inviscid Burgers' equation for linear initial condition

Subrahmanyan Chandrasekhar provided the explicit solution in 1943 when the initial condition is linear, i.e., ${\displaystyle f(x)=ax+b}$, where a and b are constants.^[10] The explicit solution is

$${\displaystyle u(x,t)={\frac {ax+b}{at+1}}.}$$

This solution is also the complete integral of the inviscid Burgers' equation because it contains as many arbitrary constants as the number of independent variables appearing in the equation.^[11] Using this complete integral, Chandrasekhar obtained the general solution described for arbitrary initial conditions from the envelope of the complete integral.

Viscous Burgers' equation

This is a numerical solution of the viscous two dimensional Burgers equation using an initial Gaussian profile. We see shock formation, and dissipation of the shock due to viscosity as it travels.

The viscous Burgers' equation can be converted to a linear equation by the Cole–Hopf transformation,^[12][13][14]

$${\displaystyle u=-2\nu {\frac {1}{\phi }}{\frac {\partial \phi }{\partial x}},}$$

which turns it into the equation

$${\displaystyle {\frac {\partial }{\partial x}}\left({\frac {1}{\phi }}{\frac {\partial \phi }{\partial t}}\right)=\nu {\frac {\partial }{\partial x}}\left({\frac {1}{\phi }}{\frac {\partial ^{2}\phi }{\partial x^{2}}}\right)}$$

which can be integrated with respect to ${\displaystyle x}$ to obtain

$${\displaystyle {\frac {\partial \phi }{\partial t}}=\nu {\frac {\partial ^{2}\phi }{\partial x^{2}}}+g(t)\phi }$$

where ${\displaystyle g(t)}$ is a function that depends on boundary conditions. If ${\displaystyle g(t)=0}$ identically (e.g. if the problem is to be solved on a periodic domain), then we get the diffusion equation

$${\displaystyle {\frac {\partial \phi }{\partial t}}=\nu {\frac {\partial ^{2}\phi }{\partial x^{2}}}.}$$

The diffusion equation can be solved, and the Cole–Hopf transformation inverted, to obtain the solution to the Burgers' equation:

$${\displaystyle u(x,t)=-2\nu {\frac {\partial }{\partial x}}\ln \left\{(4\pi \nu t)^{-1/2}\int _{-\infty }^{\infty }\exp \left[-{\frac {(x-x')^{2}}{4\nu t}}-{\frac {1}{2\nu }}\int _{0}^{x'}f(x'')dx''\right]dx'\right\}.}$$

Other forms

Generalized Burgers' equation

The generalized Burgers' equation extends the quasilinear convective to more generalized form, i.e.,

$${\displaystyle {\frac {\partial u}{\partial t}}+c(u){\frac {\partial u}{\partial x}}=\nu {\frac {\partial ^{2}u}{\partial x^{2}}}.}$$

where ${\displaystyle c(u)}$ is any arbitrary function of u. The inviscid ${\displaystyle \nu =0}$ equation is still a quasilinear hyperbolic equation for ${\displaystyle c(u)>0}$ and its solution can be constructed using method of characteristics as before.^[15]

Stochastic Burgers' equation

Added space-time noise ${\displaystyle \eta (x,t)={\dot {W}}(x,t)}$, where ${\displaystyle W}$ is an ${\displaystyle L^{2}(\mathbb {R} )}$ Wiener process, forms a stochastic Burgers' equation^[16]

$${\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=\nu {\frac {\partial ^{2}u}{\partial x^{2}}}-\lambda {\frac {\partial \eta }{\partial x}}.}$$

This stochastic PDE is the one-dimensional version of Kardar–Parisi–Zhang equation in a field ${\displaystyle h(x,t)}$ upon substituting ${\displaystyle u(x,t)=-\lambda \partial h/\partial x}$.

See also

• Euler–Tricomi equation
• Chaplygin's equation
• Conservation equation
• Fokker–Planck equation

References

External links

• Burgers' Equation at EqWorld: The World of Mathematical Equations.
• Burgers' Equation at NEQwiki, the nonlinear equations encyclopedia.