Boundary value problem

In mathematics, a boundary value problem consists of a differential equation to be satisfied at all points in the interior of an interval or a region and a set of boundary conditions specifying the values of the solution or some of its derivatives everywhere on the boundary of the interval or region. Boundary value problems may be posed for ordinary differential equations as well as partial differential equations.

Boundary value problems arise in several branches of physics. Problems involving the wave equation, such as the determination of normal modes, are often stated as boundary value problems.

A large class of important boundary value problems are the Sturm-Liouville problems. The analysis of these problems involves the eigenfunctions of a differential operator.

Contents

1 Example (ordinary) 2 Example (partial) 3 See also 4 External links 5 Bibliography

Example (ordinary)

We wish to find a function y(x) which solves the following Sturm-Liouville problem:

<math>y'' = \lambda y\,<math>

and satisfy the boundary conditions

<math>y(0)=0, \quad y(\pi) = 0<math>.

We will use k to denote the square root of the absolute value of <math>\lambda<math>.

If <math>\lambda = 0<math> then

<math>y(x) = Ax + B\,<math>

solves the ODE.

Substituted boundary conditions give that both A and B are equal to zero.

For positive <math>\lambda<math> we obtain that

<math>y(x) = A e^{kx} + B e^{-kx}\,<math>

solves the ODE.

Substitution of boundary conditions again yields A = B = 0.

For negative <math>\lambda<math> it is easy to show that

<math>y(x) = A \sin(kx) + B \cos(kx)\,<math>

solves the ODE.

From the first boundary condition,

<math>0 = y(0) = A \sin(0) + B \cos(0) = B\,<math>.

Now, after the cosine is gone, we will substitute the second boundary condition:

<math>y(\pi) = A \sin(k \pi) = 0\,<math>.

So either A = 0 or k is an integer.

Thus we get that the eigenfunctions which solve the "boundary value problem" are

<math>y_n(x) = A \sin(nx) \quad n = 0,1,2,3,...<math>.

One may easily check that they satisfy the boundary conditions.

Example (partial)

Consider the elliptic eigenvalue problem (boundary value problem)

<math>\nabla^2 v + \lambda v = {\partial^2 v \over \partial x^2} + {\partial^2 v \over \partial y^2} + \lambda v = 0,\,0

with boundary conditions

<math>v(x,0) = v(x,1) = 0,\,0

<math>{\partial\over\partial x}v(0,y) = {\partial\over\partial x}v(1,y) = 0,\,0

We suppose the solution is of the form

<math>v = X(x)Y(y)\,<math>

substituting,

<math>{\partial^2\over\partial x^2}(X(x)Y(y))+{\partial^2\over\partial y}(X(x)Y(y))+\lambda X(x)Y(y)\,<math>

<math>= X''(x)Y(y)+X(x)Y''(y)+\lambda X(x)Y(y)= 0\,<math>.

Divide throughout by X(x):

<math>= {X''(x)Y(y) \over X(x)}+{X(x)Y''(y)\over X(x)}+{\lambda X(x)Y(y)\over X(x)} = {X''(x)Y(y) \over X(x)}+Y''(y)+\lambda Y(y) = 0<math>

and then by Y(y):

<math>= {X''(x)\over X(x)}+{Y''(y)+\lambda Y(y)\over Y(y)} = 0<math>.

Now X′′(x)/X(x) is a function of x only, as is (Y′′(y)+λY(y))/Y(y), so there are separation constants so

<math>{X''(x)\over X(x)} = k = {Y''(y)+\lambda Y(y)\over Y(y)}.<math>

From our boundary conditions we have

<math>v(x,0) = X(x)Y(0) = 0,\ v(x,1) = X(x)Y(1) = 0<math>,

<math>{\partial\over\partial x}v(0,x) = X'(0)Y(y) = 0,\ {\partial\over\partial x}v(1,y)=X'(1)Y(y)=0<math>

we want

<math>Y(0) = 0,\ Y(1) = 0,\ X'(0) = 0,\ X'(1) = 0<math>

which splits up into ordinary differential equations

<math>{X''(x)\over X(x)} = k<math>,

<math>X''(x) - k X(x) = 0\,<math>

and

<math>{Y''(y)+\lambda Y(y)\over Y(y)} = k<math>,

<math>Y''(y)+(\lambda-k) Y(y) = 0\,<math>

which we can evaluate the boundary conditions and solutions accordingly.

See also

initial value problem

External links

• Linear Partial Differential Equations: Exact Solutions and Boundary Value Problems (http://eqworld.ipmnet.ru/en/solutions/lpde.htm) at EqWorld: The World of Mathematical Equations.

Bibliography

• A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, 2002.de:randwertproblem