Parabolic partial differential equation

Last updated September 10, 2024

A parabolic partial differential equation is a type of partial differential equation (PDE). Parabolic PDEs are used to describe a wide variety of time-dependent phenomena in, i.a., engineering science, quantum mechanics and financial mathematics. Examples include the heat equation, time-dependent Schrödinger equation and the Black–Scholes equation.

Definition

To define the simplest kind of parabolic PDE, consider a real-valued function $u(x,y)$ of two independent real variables, $x$ and $y$ . A second-order, linear, constant-coefficient PDE for $u$ takes the form

Au_{xx}+2Bu_{xy}+Cu_{yy}+Du_{x}+Eu_{y}+F=0,

where the subscripts denote the first- and second-order partial derivatives with respect to $x$ and $y$ . The PDE is classified as parabolic if the coefficients of the principal part (i.e. the terms containing the second derivatives of $u$ ) satisfy the condition^[1]

B^{2}-AC=0.

Usually $x$ represents one-dimensional position and $y$ represents time, and the PDE is solved subject to prescribed initial and boundary conditions. Equations with $B^{2}-AC<0$ are termed elliptic while those with $B^{2}-AC>0$ are hyperbolic. The name "parabolic" is used because the assumption on the coefficients is the same as the condition for the analytic geometry equation $Ax^{2}+2Bxy+Cy^{2}+Dx+Ey+F=0$ to define a planar parabola.

The basic example of a parabolic PDE is the one-dimensional heat equation

u_{t}=\alpha \,u_{xx},

where $u(x,t)$ is the temperature at position $x$ along a thin rod at time $t$ and $\alpha$ is a positive constant called the thermal diffusivity.

The heat equation says, roughly, that temperature at a given time and point rises or falls at a rate proportional to the difference between the temperature at that point and the average temperature near that point. The quantity $u_{xx}$ measures how far off the temperature is from satisfying the mean value property of harmonic functions.

The concept of a parabolic PDE can be generalized in several ways. For instance, the flow of heat through a material body is governed by the three-dimensional heat equation

u_{t}=\alpha \,\Delta u,

where

\Delta u:={\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}+{\frac {\partial ^{2}u}{\partial z^{2}}},

denotes the Laplace operator acting on $u$ . This equation is the prototype of a multi-dimensional parabolic PDE.^[2]

Noting that $-\Delta$ is an elliptic operator suggests a broader definition of a parabolic PDE:

u_{t}=-Lu,

where $L$ is a second-order elliptic operator (implying that $L$ must be positive; a case where $u_{t}=+Lu$ is considered below).

A system of partial differential equations for a vector $u$ can also be parabolic. For example, such a system is hidden in an equation of the form

\nabla \cdot (a(x)\nabla u(x))+b(x)^{\text{T}}\nabla u(x)+cu(x)=f(x)

if the matrix-valued function $a(x)$ has a kernel of dimension 1.

Solution

Under broad assumptions, an initial/boundary-value problem for a linear parabolic PDE has a solution for all time. The solution $u(x,t)$ , as a function of $x$ for a fixed time $t>0$ , is generally smoother than the initial data $u(x,0)=u_{0}(x)$ .

For a nonlinear parabolic PDE, a solution of an initial/boundary-value problem might explode in a singularity within a finite amount of time. It can be difficult to determine whether a solution exists for all time, or to understand the singularities that do arise. Such interesting questions arise in the solution of the Poincaré conjecture via Ricci flow.^{[ citation needed ]}

Backward parabolic equation

One occasionally encounters a so-called backward parabolic PDE, which takes the form $u_{t}=Lu$ (note the absence of a minus sign).

An initial-value problem for the backward heat equation,

{\begin{cases}u_{t}=-\Delta u&{\textrm {on}}\ \ \Omega \times (0,T),\\u=0&{\textrm {on}}\ \ \partial \Omega \times (0,T),\\u=f&{\textrm {on}}\ \ \Omega \times \left\{0\right\}.\end{cases}}

is equivalent to a final-value problem for the ordinary heat equation,

{\begin{cases}u_{t}=\Delta u&{\textrm {on}}\ \ \Omega \times (0,T),\\u=0&{\textrm {on}}\ \ \partial \Omega \times (0,T),\\u=f&{\textrm {on}}\ \ \Omega \times \left\{T\right\}.\end{cases}}

Similarly to a final-value problem for a parabolic PDE, an initial-value problem for a backward parabolic PDE is usually not well-posed (solutions often grow unbounded in finite time, or even fail to exist). Nonetheless, these problems are important for the study of the reflection of singularities of solutions to various other PDEs.^[3]

Notes

Related Research Articles

<span class="mw-page-title-main">Wave equation</span> Differential equation important in physics

The wave equation is a second-order linear partial differential equation for the description of waves or standing wave fields such as mechanical waves or electromagnetic waves. It arises in fields like acoustics, electromagnetism, and fluid dynamics.

<span class="mw-page-title-main">Partial differential equation</span> Type of differential equation

In mathematics, a partial differential equation (PDE) is an equation which computes a function between various partial derivatives of a multivariable function.

In mathematics, mathematical physics and the theory of stochastic processes, a harmonic function is a twice continuously differentiable function $where U is an open subset of ⁠ ⁠ that satisfies Laplace's equation, that is, everywhere on U . This is usually written as or$

<span class="mw-page-title-main">Heat equation</span> Partial differential equation describing the evolution of temperature in a region

In mathematics and physics, the heat equation is a certain partial differential equation. Solutions of the heat equation are sometimes known as caloric functions. The theory of the heat equation was first developed by Joseph Fourier in 1822 for the purpose of modeling how a quantity such as heat diffuses through a given region.

In mathematics, Frobenius' theorem gives necessary and sufficient conditions for finding a maximal set of independent solutions of an overdetermined system of first-order homogeneous linear partial differential equations. In modern geometric terms, given a family of vector fields, the theorem gives necessary and sufficient integrability conditions for the existence of a foliation by maximal integral manifolds whose tangent bundles are spanned by the given vector fields. The theorem generalizes the existence theorem for ordinary differential equations, which guarantees that a single vector field always gives rise to integral curves; Frobenius gives compatibility conditions under which the integral curves of r vector fields mesh into coordinate grids on r-dimensional integral manifolds. The theorem is foundational in differential topology and calculus on manifolds.

Second-order linear partial differential equations (PDEs) are classified as either elliptic, hyperbolic, or parabolic. Any second-order linear PDE in two variables can be written in the form

In mathematics, the method of characteristics is a technique for solving partial differential equations. Typically, it applies to first-order equations, although more generally the method of characteristics is valid for any hyperbolic and parabolic partial differential equation. The method is to reduce a partial differential equation to a family of ordinary differential equations along which the solution can be integrated from some initial data given on a suitable hypersurface.

In the mathematical fields of differential equations and geometric analysis, the maximum principle is one of the most useful and best known tools of study. Solutions of a differential inequality in a domain D satisfy the maximum principle if they achieve their maxima at the boundary of D.

In mathematics, a hyperbolic partial differential equation of order $is a partial differential equation (PDE) that, roughly speaking, has a well-posed initial value problem for the first derivatives. More precisely, the Cauchy problem can be locally solved for arbitrary initial data along any non-characteristic hypersurface. Many of the equations of mechanics are hyperbolic, and so the study of hyperbolic equations is of substantial contemporary interest. The model hyperbolic equation is the wave equation. In one spatial dimension, this is The equation has the property that, if u and its first time derivative are arbitrarily specified initial data on the line t = 0, then there exists a solution for all time t .$

In mathematics, Harnack's inequality is an inequality relating the values of a positive harmonic function at two points, introduced by A. Harnack. Harnack's inequality is used to prove Harnack's theorem about the convergence of sequences of harmonic functions. J. Serrin, and J. Moser generalized Harnack's inequality to solutions of elliptic or parabolic partial differential equations. Such results can be used to show the interior regularity of weak solutions.

In the mathematical study of heat conduction and diffusion, a heat kernel is the fundamental solution to the heat equation on a specified domain with appropriate boundary conditions. It is also one of the main tools in the study of the spectrum of the Laplace operator, and is thus of some auxiliary importance throughout mathematical physics. The heat kernel represents the evolution of temperature in a region whose boundary is held fixed at a particular temperature, such that an initial unit of heat energy is placed at a point at time $t = 0$ .

In numerical analysis, finite-difference methods (FDM) are a class of numerical techniques for solving differential equations by approximating derivatives with finite differences. Both the spatial domain and time domain are discretized, or broken into a finite number of intervals, and the values of the solution at the end points of the intervals are approximated by solving algebraic equations containing finite differences and values from nearby points.

In mathematics, an elliptic boundary value problem is a special kind of boundary value problem which can be thought of as the stable state of an evolution problem. For example, the Dirichlet problem for the Laplacian gives the eventual distribution of heat in a room several hours after the heating is turned on.

In mathematics, the Babuška–Lax–Milgram theorem is a generalization of the famous Lax–Milgram theorem, which gives conditions under which a bilinear form can be "inverted" to show the existence and uniqueness of a weak solution to a given boundary value problem. The result is named after the mathematicians Ivo Babuška, Peter Lax and Arthur Milgram.

In mathematics, the p-Laplacian, or the p-Laplace operator, is a quasilinear elliptic partial differential operator of 2nd order. It is a nonlinear generalization of the Laplace operator, where $is allowed to range over . It is written as$

In mathematics, the Hopf lemma, named after Eberhard Hopf, states that if a continuous real-valued function in a domain in Euclidean space with sufficiently smooth boundary is harmonic in the interior and the value of the function at a point on the boundary is greater than the values at nearby points inside the domain, then the derivative of the function in the direction of the outward pointing normal is strictly positive. The lemma is an important tool in the proof of the maximum principle and in the theory of partial differential equations. The Hopf lemma has been generalized to describe the behavior of the solution to an elliptic problem as it approaches a point on the boundary where its maximum is attained.

The Kansa method is a computer method used to solve partial differential equations. Its main advantage is it is very easy to understand and program on a computer. It is much less complicated than the finite element method. Another advantage is it works well on multi variable problems. The finite element method is complicated when working with more than 3 space variables and time.

In mathematics, the walk-on-spheres method (WoS) is a numerical probabilistic algorithm, or Monte-Carlo method, used mainly in order to approximate the solutions of some specific boundary value problem for partial differential equations (PDEs). The WoS method was first introduced by Mervin E. Muller in 1956 to solve Laplace's equation, and was since then generalized to other problems.

<span class="mw-page-title-main">Stokes problem</span>

In fluid dynamics, Stokes problem also known as Stokes second problem or sometimes referred to as Stokes boundary layer or Oscillating boundary layer is a problem of determining the flow created by an oscillating solid surface, named after Sir George Stokes. This is considered one of the simplest unsteady problems that has an exact solution for the Navier–Stokes equations. In turbulent flow, this is still named a Stokes boundary layer, but now one has to rely on experiments, numerical simulations or approximate methods in order to obtain useful information on the flow.

The Cole–Hopf transformation is a change of variables that allows to transform a special kind of parabolic partial differential equations (PDEs) with a quadratic nonlinearity into a linear heat equation. In particular, it provides an explicit formula for fairly general solutions of the PDE in terms of the initial datum and the heat kernel.

References

Taylor, Michael E. (1975). "Reflection of singularities of solutions to systems of differential equations". Communications on Pure and Applied Mathematics. 28 (4): 457–478. CiteSeerX 10.1.1.697.9255 . doi:10.1002/cpa.3160280403. ISSN 0010-3640.
Zauderer, Erich (2006). Partial Differential Equations of Applied Mathematics. Hoboken, N.J: Wiley-Interscience. ISBN 978-0-471-69073-3. OCLC 70158521.