Parabolic partial differential equation

Last updated

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, including heat conduction, particle diffusion, and pricing of derivative investment instruments.

Contents

Definition

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

and this PDE is classified as being parabolic if the coefficients satisfy the condition

Usually represents one-dimensional position and represents time, and the PDE is solved subject to prescribed initial and boundary conditions.

The name "parabolic" is used because the assumption on the coefficients is the same as the condition for the analytic geometry equation to define a planar parabola.

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

where is the temperature at time and at position along a thin rod, and is a positive constant (the thermal diffusivity). The symbol signifies the partial derivative of with respect to the time variable , and similarly is the second partial derivative with respect to . For this example, plays the role of in the general second-order linear PDE: , , and the other coefficients are zero.

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 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,

where

denotes the Laplace operator acting on . This equation is the prototype of a multi-dimensional parabolic PDE.

Noting that is an elliptic operator suggests a broader definition of a parabolic PDE:

where is a second-order elliptic operator (implying that must be positive; a case where is considered below).

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

if the matrix-valued function has a kernel of dimension 1.

Parabolic PDEs can also be nonlinear. For example, Fisher's equation is a nonlinear PDE that includes the same diffusion term as the heat equation but incorporates a linear growth term and a nonlinear decay term.

Solution

Under broad assumptions, an initial/boundary-value problem for a linear parabolic PDE has a solution for all time. The solution , as a function of for a fixed time , is generally smoother than the initial data .

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 (note the absence of a minus sign).

An initial-value problem for the backward heat equation,

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

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. [1]

Examples

See also

Related Research Articles

<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.

<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.

<span class="mw-page-title-main">Elliptic operator</span> Type of differential operator

In the theory of partial differential equations, elliptic operators are differential operators that generalize the Laplace operator. They are defined by the condition that the coefficients of the highest-order derivatives be positive, which implies the key property that the principal symbol is invertible, or equivalently that there are no real characteristic directions.

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 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.

<span class="mw-page-title-main">Differential equation</span> Type of functional equation (mathematics)

In mathematics, a differential equation is an equation that relates one or more unknown functions and their derivatives. In applications, the functions generally represent physical quantities, the derivatives represent their rates of change, and the differential equation defines a relationship between the two. Such relations are common; therefore, differential equations play a prominent role in many disciplines including engineering, physics, economics, and biology.

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

In mathematics, the viscosity solution concept was introduced in the early 1980s by Pierre-Louis Lions and Michael G. Crandall as a generalization of the classical concept of what is meant by a 'solution' to a partial differential equation (PDE). It has been found that the viscosity solution is the natural solution concept to use in many applications of PDE's, including for example first order equations arising in dynamic programming, differential games or front evolution problems, as well as second-order equations such as the ones arising in stochastic optimal control or stochastic differential games.

In mathematics, constraint counting is counting the number of constraints in order to compare it with the number of variables, parameters, etc. that are free to be determined, the idea being that in most cases the number of independent choices that can be made is the excess of the latter over the former.

In mathematics, Harnack's inequality is an inequality relating the values of a positive harmonic function at two points, introduced by A. Harnack (1887). Harnack's inequality is used to prove Harnack's theorem about the convergence of sequences of harmonic functions. J. Serrin (1955), 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.

<span class="mw-page-title-main">Navier–Stokes existence and smoothness</span> Millennium Prize Problem

The Navier–Stokes existence and smoothness problem concerns the mathematical properties of solutions to the Navier–Stokes equations, a system of partial differential equations that describe the motion of a fluid in space. Solutions to the Navier–Stokes equations are used in many practical applications. However, theoretical understanding of the solutions to these equations is incomplete. In particular, solutions of the Navier–Stokes equations often include turbulence, which remains one of the greatest unsolved problems in physics, despite its immense importance in science and engineering.

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 interval are discretized, or broken into a finite number of steps, and the value of the solution at these discrete points is approximated by solving algebraic equations containing finite differences and values from nearby points.

<span class="mw-page-title-main">Elliptic boundary value problem</span>

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 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 numerical analysis, the FTCS method is a finite difference method used for numerically solving the heat equation and similar parabolic partial differential equations. It is a first-order method in time, explicit in time, and is conditionally stable when applied to the heat equation. When used as a method for advection equations, or more generally hyperbolic partial differential equations, it is unstable unless artificial viscosity is included. The abbreviation FTCS was first used by Patrick Roache.

In the finite element method for the numerical solution of elliptic partial differential equations, the stiffness matrix is a matrix that represents the system of linear equations that must be solved in order to ascertain an approximate solution to the differential equation.

In scientific computation and simulation, the method of fundamental solutions (MFS) is a technique for solving partial differential equations based on using the fundamental solution as a basis function. The MFS was developed to overcome the major drawbacks in the boundary element method (BEM) which also uses the fundamental solution to satisfy the governing equation. Consequently, both the MFS and the BEM are of a boundary discretization numerical technique and reduce the computational complexity by one dimensionality and have particular edge over the domain-type numerical techniques such as the finite element and finite volume methods on the solution of infinite domain, thin-walled structures, and inverse problems.

In numerical mathematics, the boundary knot method (BKM) is proposed as an alternative boundary-type meshfree distance function collocation scheme.

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.

<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 method of solving parabolic partial differential equations (PDEs) with a quadratic nonlinearity of the form:

References

  1. Taylor, M. E. (1975), "Reflection of singularities of solutions to systems of differential equations", Comm. Pure Appl. Math., 28 (4): 457–478, CiteSeerX   10.1.1.697.9255 , doi:10.1002/cpa.3160280403

Further reading