Fundamental solution

Fundamental solution

In mathematics, a fundamental solution for a linear partial differential operator "L" is a formulation in the language of distribution theory of the older idea of a Green's function. In terms of the Dirac delta function δ("x"), a fundamental solution "f" is the solution of the inhomogeneous equation

:"Lf" = δ("x").

Here "f" is "a priori" only assumed to be a Schwartz distribution.

This concept was long known for the Laplacian in two and three dimensions. It was investigated for all dimensions for the Laplacian by Marcel Riesz. The existence of a fundamental solution for any operator with constant coefficients — the most important case, directly linked to the possibility of using convolution to solve an arbitrary right hand side — was shown by Malgrange and Ehrenpreis.

Example

Consider the following differential equation "Lf" = sin(x) with L as,

: L=frac{partial^2}{partial x^2}

The fundamental solutions can be obtained by solving "Lf" = δ("x"), explicitly,

: frac{partial^2}{partial x^2} F(x) = delta(x)

Since for the Heaviside function "H" we have

:"H"′("x") = δ("x")

there is a solution

:"F"′("x") = "H"(x) + "C".

Here "C" is an arbitrary constant. For convenience, set

:"C" = − 1/2.

After integrating and taking the integration constant as zero, we get

: F(x)=frac{1}{2} |x|

Fundamental solutions for some partial differential equations

Laplace equation

: [- abla^2] Phi(mathbf{x},mathbf{x}') = delta(mathbf{x}-mathbf{x}')

The fundamental solutions in two and three dimensions are

: Phi_{2D}(mathbf{x},mathbf{x}')=-frac{1}{2pi}ln|mathbf{x}-mathbf{x}'|,quad Phi_{3D}(mathbf{x},mathbf{x}')=frac{1}{4pi|mathbf{x}-mathbf{x}'
Helmholtz equation where the parameter "k" is real and the fundamental solution a modified Bessel function.

: [- abla^2+k^2] Phi(mathbf{x},mathbf{x}') = delta(mathbf{x}-mathbf{x}')

The two and three dimensional Helmholtz equations have the fundamental solutions : Phi_{2D}(mathbf{x},mathbf{x}')=frac{1}{2pi}K_0(k|mathbf{x}-mathbf{x}'|),quadPhi_{3D}(mathbf{x},mathbf{x}')=frac{1}{4pi|mathbf{x}-mathbf{x}'exp(-k|mathbf{x}-mathbf{x}'|)

Biharmonic equation

: [- abla^4] Phi(mathbf{x},mathbf{x}') = delta(mathbf{x}-mathbf{x}')

The biharmonic equation has the fundamental solutions

:Phi_{2D}(mathbf{x},mathbf{x}')=-frac{8pi}

Motivation

The motivation to find the fundamental solution is because once one finds the fundamental solution, it is easy to find the desired solution of the original equation. In fact, this process is achieved by convolution.

Fundamental solutions also play an important role in the numerical solution of partial differential equations by the boundary element method.

Application to the example

Consider the operator L, mentioned in the example.

: frac{partial^2}{partial x^2} f(x) = sin(x)

Since we have found the fundamental solution, we can easily find the solution of the original equation by convolution,

: int_{-infty}^{infty} frac{1}{2}|x - y|sin(y)dy

Proof that the convolution is the desired solution

Denote the convolution operation as

:"f"*"g".

Say we are trying to find the solution of

:"Lf" = "g"("x").

When applying the differential operator, "L", to the convolution it is known that

:"L"("f"*"g")=("Lf")*"g",

provided "L" has constant coefficients.

If "f" is the fundamental solution, the RHS reduces to

:δ*"g".

It is straightforward to verify that this is in fact "g"("x") (in other words the delta function acts as identity element for convolution). Summing up,

: L(F*g)=(LF)*g=delta(x)*g(x)=int_{-infty}^{infty} delta (x-y) g(y) dy=g(x)

Therefore, if "F" is the fundamental solution, the convolution "F"*"g" is the solution of "Lf" = "g"("x").

See also

* parametrix


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • Fundamental theorem of algebra — In mathematics, the fundamental theorem of algebra states that every non constant single variable polynomial with complex coefficients has at least one complex root. Equivalently, the field of complex numbers is algebraically closed.Sometimes,… …   Wikipedia

  • Fundamental unit (number theory) — In algebraic number theory, a fundamental unit is a generator for the torsion free unit group of the ring of integers of a number field, when that group is infinite cyclic. See also Dirichlet s unit theorem.For rings of the form mathbb{Z} [sqrt… …   Wikipedia

  • Fundamental theorem of linear algebra — In mathematics, the fundamental theorem of linear algebra makes several statements regarding vector spaces. These may be stated concretely in terms of the rank r of an m times; n matrix A and its LDU factorization::PA=LDUwherein P is a… …   Wikipedia

  • Fundamental theorem of combinatorial enumeration — The fundamental theorem of combinatorial enumeration is a theorem in combinatorics that solves the enumeration problem of labelled and unlabelled combinatorial classes. The unlabelled case is based on the Pólya enumeration theorem.This theorem is …   Wikipedia

  • Solución fundamental — En matematicas, una solución fundamental para un operador diferencial con derivadas parciales L es una formulación en el lenguaje de la teoría de distribuciones proveniente de la antigua idea de la función de Green. En terminos de la función… …   Wikipedia Español

  • Ideal solution — In chemistry, an ideal solution or ideal mixture is a solution in which the enthalpy of solution is zero; [ A to Z of Thermodynamics Pierre Perrot ISBN 0198565569] the closer to zero the enthalpy of solution, the more ideal the behavior of the… …   Wikipedia

  • bandaid solution — /bændeɪd səˈluʃən/ (say bandayd suh loohshuhn) noun an answer to a problem which is superficial by nature and does not address fundamental issues. Also, bandaid remedy …  

  • Heat equation — The heat equation is an important partial differential equation which describes the distribution of heat (or variation in temperature) in a given region over time. For a function of three spatial variables ( x , y , z ) and one time variable t ,… …   Wikipedia

  • Pell's equation — is any Diophantine equation of the form:x^2 ny^2=1,where n is a nonsquare integer and x and y are integers. Trivially, x = 1 and y = 0 always solve this equation. Lagrange proved that for any natural number n that is not a perfect square there… …   Wikipedia

  • Dirac delta function — Schematic representation of the Dirac delta function by a line surmounted by an arrow. The height of the arrow is usually used to specify the value of any multiplicative constant, which will give the area under the function. The other convention… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”