Burgers' equation explained

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

u(x,t) and diffusion coefficient (or kinematic viscosity, as in the original fluid mechanical context) \nu , the general form of Burgers' equation (also known as viscous Burgers' equation) in one space dimension is the dissipative system:(\partial u)/(\partial t) + u (\partial u)/(\partial x) = \nu (\partial ^{2} u)/(\partial x ^{2}) .

The term

u\partial u/\partial x can also rewritten as \partial(u ^{2/2)/\partial} x . When the diffusion term is absent (i.e. \nu =0 ), Burgers' equation becomes the inviscid Burgers' equation:(\partial u)/(\partial t) + u (\partial u)/(\partial x) = 0 ,

which is a prototype for conservation equations that can develop discontinuities (shock waves).

The reason for the formation of sharp gradients for small values of

\nu becomes intuitively clear when one examines the left-hand side of the equation. The term \partial/\partial t + u \partial/\partial x is evidently a wave operator describing a wave propagating in the positive x -direction with a speed u . Since the wave speed is u , regions exhibiting large values of u will be propagated rightwards quicker than regions exhibiting smaller values of u ; in other words, if u is decreasing in the x -direction, initially, then larger u 's that lie in the backside will catch up with smaller u 's on the front side. The role of the right-side diffusive term is essentially to stop the gradient becoming infinite.

Inviscid Burgers' equation

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

(\partial u)/(\partial t) + u (\partial u)/(\partial x) = 0 , u(x, 0 ) = f(x)

can be constructed by the method of characteristics. Let

t be the parameter characterising any given characteristics in the x -t plane, then the characteristic equations are given by(dx)/(dt) = u, (du)/(dt) =0 .

Integration of the second equation tells us that

u is constant along the characteristic and integration of the first equation shows that the characteristics are straight lines, i.e.,u = c, x = ut + \xi

where

\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 u at x -axis is known from the initial condition and the fact that u is unchanged as we move along the characteristic emanating from each point x = \xi , we write u = c = f(\xi) on each characteristic. Therefore, the family of trajectories of characteristics parametrized by \xi isx = f(\xi) t + \xi.

Thus, the solution is given by

u(x,t) = f(\xi) = f(x - ut), \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^{[8] [9]}

t _{b} = (-1)/(inf _{x} \left(f ^{\prime(x)\right)}) .

Complete integral of the inviscid Burgers' equation

The implicit solution described above containing an arbitrary function

f is called the general integral. However, the inviscid Burgers' equation, being a first-order partial differential equation, also has a complete integral which contains two arbitrary constants (for the two independent variables).^[10] Subrahmanyan Chandrasekhar provided the complete integral in 1943,^[11] which is given byu(x,t) = (ax + b)/(at +1) .

where

a and b are arbitrary constants. The complete integral satisfies a linear initial condition, i.e., f(x) = ax + b . One can also construct the geneal integral using the above complete integral.

Viscous Burgers' equation

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

u(x,t) = -2 \nu (\partial)/(\partial x) ln \varphi(x,t),

which turns it into the equation

2 \nu (\partial)/(\left[ \partial x) (1)/(\left( \varphi) (\partial\varphi)/(\partial t) - \nu (\partial ^{2\varphi})/(\partial x ^{2}) \right)\right] =0 ,

which can be integrated with respect to

x to obtain(\partial\varphi)/(\partial t) - \nu (\partial ^{2\varphi})/(\partial x ^{2}) = \varphi (df(t))/(dt) ,

where

df/dt is an arbitrary function of time. Introducing the transformation \varphi\to \varphi e ^{f} (which does not affect the function u(x,t) ), the required equation reduces to that of the heat equation^[15] (\partial\varphi)/(\partial t) = \nu (\partial ^{2\varphi})/(\partial x ^{2}) .

The diffusion equation can be solved. That is, if

\varphi(x, 0 ) = \varphi _{0(x)} , then\varphi(x,t) = (1)/(\sqrt{ 4 \pi \nu t) }\int_^\infty \varphi_0(x') \exp \left[-\frac{(x-x')^2}{4\nu t}\right]dx'.

The initial function

\varphi _{0(x)} is related to the initial function u(x, 0 ) = f(x) byln \varphi _{0(x)} = - (1)/(2 \nu) \int ^{x} _{0} f(x') dx',

where the lower limit is chosen arbitrarily. Inverting the Cole–Hopf transformation, we have

u(x,t) =-2 \nu (\partial)/(ln \left\{ \partial x) (1)/(\sqrt{ 4 \pi \nu t) }\int_^\infty \exp\left[-\frac{(x-x')^2}{4\nu t} - \frac{1}{2\nu}\int_0^{x'}f(x'')dx''\right]dx'\right\}

which simplifies, by getting rid of the time-dependent prefactor in the argument of the logarthim, to

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

This solution is derived from the solution of the heat equation for

\varphi that decays to zero as x\to\pm inf ty ; other solutions for u can be obtained starting from solutions of \varphi that satisfies different boundary conditions.

Some explicit solutions of the viscous Burgers' equation

Explicit expressions for the viscous Burgers' equation are available. Some of the physically relevant solutions are given below:^[16]

Steadily propagating traveling wave

If

u(x, 0 ) = f(x) is such that f( -inf ty) = f ^{+} and f( +inf ty) = f ^{-} and f'(x)< 0 , then we have a traveling-wave solution (with a constant speed c = (f ^{++f} ^{-)/2} ) given byu(x,t) = c - (f ^{+-f} ^{-})/(\tanh\left[ 2) (f ^{+-f} ^{-})/(4 \nu) (x - ct)\right].

This solution, that was originally derived by Harry Bateman in 1915, is used to describe the variation of pressure across a weak shock wave. When

f ^{+=2} and f ^{-=0} tou(x,t) = (2)/(1+ e ^{(x - t)/(\nu)})

with

c =1 .

Delta function as an initial condition

If

u(x, 0 ) = 2 \nu Re \delta(x) , where Re (say, the Reynolds number) is a constant, then we have^[17] u(x,t) = \sqrt{ (\nu)/(\pi t) } \left[\frac{(e^{Re}-1)e^{-x^2/4\nu t}}{1 + (e^{Re}-1) \mathrm{erfc}(x/\sqrt{4\nu t})/\sqrt{2}}\right].

In the limit

Re\to 0 , the limiting behaviour is a diffusional spreading of a source and therefore is given byu(x,t) = (2 \nu Re)/(\sqrt{ 4 \pi \nu t) } \exp\left(-\frac\right).

On the other hand, In the limit

Re\to inf ty , the solution approaches that of the aforementioned Chandrasekhar's shock-wave solution of the inviscid Burgers' equation and is given byu(x,t) = \begin{cases} (x)/(t) , 0 <x<\sqrt{ 2 \nu Ret},\\ 0 , otherwise .\end{cases}

The shock wave location and its speed are given by

x = \sqrt{ 2 \nu Re t} and \sqrt{\nu Re/t}.

N-wave solution

The N-wave solution comprises a compression wave followed by a rarafaction wave. A solution of this type is given by

u(x,t) = (x)/(t) \left[ 1 + (1)/(\sqrt{ e ^{Re _{0-1}}) (t)/(t _{0}) }\exp\left(-\frac\right)\right]^

where

R _{0} may be regarded as an initial Reynolds number at time t = t _{0} and Re(t) = ( 1 / 2 \nu) \int ^{inf ty} _{0} udx =ln ( 1+ \sqrt{\tau/t}) with \tau = t _{0} \sqrt{e ^{Re _{0}} -1 } , may be regarded as the time-varying Reynold number.

Other forms

Multi-dimensional Burgers' equation

In two or more dimensions, the Burgers' equation becomes

(\partial u)/(\partial t) + u ⋅ \nabla u = \nu \nabla ^{2} u.

One can also extend the equation for the vector field

u , as in(\partial u)/(\partial t) + u ⋅ \nabla u = \nu \nabla ^{2} u.

Generalized Burgers' equation

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

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

where

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

Stochastic Burgers' equation

Added space-time noise

η (x,t) = • W(x,t) , where W is an L ^{2(R)} Wiener process, forms a stochastic Burgers' equation^[19] (\partial u)/(\partial t) + u (\partial u)/(\partial x) = \nu (\partial ^{2} u)/(-λ \partial x ^{2}) (\partial η)/(\partial x) .

This stochastic PDE is the one-dimensional version of Kardar–Parisi–Zhang equation in a field

h(x,t) upon substituting u(x,t) =-λ \partial h/\partial x .

See also

• Chaplygin's equation
• Conservation equation
• Euler–Tricomi equation
• Fokker–Planck equation
• KdV-Burgers equation

External links

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

Notes and References

1. Misra . Souren . Raghurama Rao . S. V. . Bobba . Manoj Kumar . 2010-09-01 . Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation . International Journal of Computational Fluid Dynamics . 24 . 8 . 303–315 . 10.1080/10618562.2010.523518 . 2010IJCFD..24..303M . 123001189 . 1061-8562.
2. It relates to the Navier–Stokes momentum equation with the pressure term removed : here the variable is the flow speed y=u
3. It arises from Westervelt equation with an assumption of strictly forward propagating waves and the use of a coordinate transformation to a retarded time frame: here the variable is the pressure
4. Musha . Toshimitsu . Higuchi . Hideyo . 1978-05-01 . Traffic Current Fluctuation and the Burgers Equation . Japanese Journal of Applied Physics . en . 17 . 5 . 811 . 10.1143/JJAP.17.811 . 1978JaJAP..17..811M . 121252757 . 1347-4065.
5. Bateman, H.. 1915. Some recent researches on the motion of fluids. Monthly Weather Review. 43. 4. 163–170. 10.1175/1520-0493(1915)43<163:SRROTM>2.0.CO;2. 1915MWRv...43..163B . free.
6. Whitham, G. B. (2011). Linear and nonlinear waves (Vol. 42). John Wiley & Sons.
7. Burgers, J. M.. A Mathematical Model Illustrating the Theory of Turbulence. Advances in Applied Mechanics. 1. 171–199. 10.1016/S0065-2156(08)70100-5. 1948. 9780123745798.
8. Book: Olver, Peter J.. Introduction to Partial Differential Equations. Undergraduate Texts in Mathematics . Springer. 2013. 978-3-319-02098-3. Online. 37. 10.1007/978-3-319-02099-0 . 220617008 .
9. Web site: Cameron . Maria . February 29, 2024 . Notes on Burger's Equation . February 29, 2024 . University of Maryland Mathematics Department, Maria Cameron's personal website.
10. Book: Forsyth, A. R. . Andrew Forsyth . 1903 . A Treatise on Differential Equations . London . Macmillan .
11. Chandrasekhar . S. . 1943 . . Report No. 423 . Ballistic Research Laboratories .
12. Julian . Cole . Julian Cole . 1951 . On a quasi-linear parabolic equation occurring in aerodynamics . Quarterly of Applied Mathematics . 9 . 3 . 225–236 . 10.1090/qam/42889 . 43633894 . free .
13. Eberhard Hopf. The partial differential equation u_t + uu_x = μu_xx . Communications on Pure and Applied Mathematics. September 1950. 3. 3. 201–230 . 10.1002/cpa.3160030302 . Eberhard Hopf.
14. Book: Kevorkian, J. . Partial Differential Equations: Analytical Solution Techniques . Belmont . Wadsworth . 1990 . 0-534-12216-7 . 31–35 .
15. Landau, L. D., & Lifshitz, E. M. (2013). Fluid mechanics: Landau And Lifshitz: course of theoretical physics, Volume 6 (Vol. 6). Elsevier. Page 352-354.
16. Salih, A. "Burgers’ Equation." Indian Institute of Space Science and Technology, Thiruvananthapuram (2016).
17. Whitham, Gerald Beresford. Linear and nonlinear waves. John Wiley & Sons, 2011.
18. Courant, R., & Hilbert, D. Methods of Mathematical Physics. Vol. II.
19. W. . Wang . A. J. . Roberts . Diffusion Approximation for Self-similarity of Stochastic Advection in Burgers' Equation . Communications in Mathematical Physics . 333 . 1287–1316 . 2015 . 3 . 10.1007/s00220-014-2117-7 . 1203.0463 . 2015CMaPh.333.1287W . 119650369 .