Introduction to Mathematical Physics/Some mathematical problems and their solution/Boundary, spectral and evolution problems

From Wikibooks, open books for an open world
Jump to: navigation, search

In order to help the reading of the next chapters, a quick classification of various mathematical problems encountered in the modelization of physical phenomena is proposed in the present chapter. More precisely, the problems considered in this chapter are those that can be reduced to the finding of the solution of a partial differential equation (PDE). Indeed, for many physicists, to provide a model of a phenomenon means to provide a PDE describing this phenomenon. They can be boundary problems, spectral problems, evolution problems. General ideas about the methods of exact and approximate solving of those PDE is also proposed[1]. This chapter contains numerous references to the "physical" part of this book which justify the interest given to those mathematical problems.

In classical books about PDE, equations are usually classified into three categories: hyperbolic, parabolic and elliptic equation. This classification is connected to the proof of existence and unicity of the solutions rather than to the actual way of obtaining the solution. We present here another classification connected to the way one obtains the solutions: we distinguish mainly boundary problems and evolution problems.

Let us introduce boundary problems:

proeq

Problem: (Boundary problem) Find u such that:

  1. Lu=f,\ u\in E, x\in\Omega

  2. u satisfies the boundary conditions.

This class of problems can be solved by integral methods (see section chapmethint) and by variational methods (see section chapmetvar). Let us introduce a second class of problems: evolution problems. Initial conditions are usually implied by time variables. One knows at time t_0 the function u(x,y,t_0) and try to get the values of u(x,y,t) for t greater than t_0. This our second class of problem:


probevollin

Problem: (Evolution problem) Find u\in V such that:

  1. eqevollinspec

    \frac{\partial u}{\partial t}=Lu+f,\ u\in E, x\in\Omega

  2. u satisfies the boundary conditions.
  3. u satisfies initial conditions.

Of course in some problems time can act "like" a space variable. This is the case when shooting problems where u should satisfy a condition at t=t_0 and another condition at t=t_1.

The third class of problem (spectral problem) often arises from the solving (see section chapmethspec) of linear evolution problems (evolution problems where the operator L is linear):

prospec1

Problem: (Spectral problem) Find u and \lambda such that:


  1. 
Lu-\lambda u=f,\ u\in E, x\in\Omega

  2. u satisfies the boundary conditions.


The difference between boundary problems and evolution problem rely on the different role played by the different independent variables on which depend the unknown function. Space variables usually implies boundary conditions. For instance the elevation u(x,y) of a membrane which is defined for each position x,y in a domain \Omega delimited by a boundary \partial \Omega should be zero on the boundary. If the equation satisfied by u at equilibrium and at position (x,y) in \Omega is:

Lu=f,

then L should be a differential operator acting on space variables which is at least of the second order.

Let us now develop some ideas about boundary conditions. In the case of ordinary differential equations (ODE) the unicity of solution is connected to the initial conditions via the Cauchy theorem. For instance to determine fully a solution of an equation:

\frac{d^2u}{dt^2}=f(t)

one needs to know both u(0) and \frac{du}{dt}(0). Boundary conditions are more subtle, since the space can have a dimension greater than one. Let us consider a 1-D boundary problem. The equation is then an ODE that can be written:

eqgenerwrit

Lu=f

Boundary conditions are imposed for two (at least) different values of the space variables x. Thus the operator should be at least of second order. Let us take L=\frac{d^2}{dx^2}. Equation eqgenerwrit is then called Laplace equation\index{Laplace equation}. The elevation of a string obey to such an equation, where f is the distribution of weight on the string. A clamped string (see Fig. figcordef) corresponds to the case where the boundary conditions are:


u(0)=0,\ u(L)=0;

Those conditions are called Dirichlet conditions. A sliding string (see Fig. figcordeg) corresponds to the case where the boundary conditions are:


\frac{du}{dx}(0)=0,\ \frac{du}{dx}(L)=0;

Those conditions are called Neumann conditions. Let us recall here the definition of adjoint operator\index{adjoint operator}:

Definition: The adjoint operator in a Hilbert space \mathcal H of an operator L is the operator L^* such that for all u and h in \mathcal H, one has  \mathrel{<} Lu,h\mathrel{>} = \mathrel{<} u,L^*h\mathrel{>} .

figcordef

The elevation of a clamped string verifies a Laplace equation with Dirichlet boundary conditions.

One can show that \mathcal H, the space of functions zero in 0 and L is a Hilbert space for the scalar product \mathrel{<} u,v \mathrel{>}=\int_0^L u(x)v(x) dx. (One speaks of Sobolev space, see section secvafor). Moreover, one can show that the adjoint of L is L^*=\frac{\partial^{2}}{\partial x^{2}}:

eqadjoimq

<Lu,v>=\int_0^l \frac{d^2u}{dx^2}vdx=
\int_0^lu\frac{d^2v}{dx^2}dx+[v\frac{du}{dx}-u\frac{dv}{dx}]^l_0

As L^*=L, L is called self-adjoint .

figcordeg

The elevation of a sliding string verifies a Laplace equation with Neumann boundary conditions.

One can also show that the space \mathcal H_2 of function with derivative zero at x=0 and at x=L is a Hilbert space for the scalar product: \mathrel{<} u,v
\mathrel{>}=\int_0^Lu(x)v(x) dx. Using equation eqadjoimq, one shows that L is also self-adjoint.

The form \index{form (definite)}


<Lu,u>=\int_o^l\frac{d^2u}{dx^2}udx=
[u\frac{du}{dx}]_0^l-\int_0^l(\frac{df}{dx})^2dx=
-\int_0^l(\frac{df}{dx})^2dx

is negative definite in the case of the clamped string and negative (no definite) in the case of the sliding string. Indeed, in this last case, function u=\mbox{ constant } makes <Lu,u> zero.

Let us conclude those remarks about boundary conditions by the compatibility between the solution u of a boundary problem


Lu=f

and the right hand member f ([#References|references]).

Theorem:

If there exists a solution v_0 for the homogeneous adjoint problem :


L^*v_0=0

then a necessary condition that the problem Lu=f has a solution is that f is orthogonal to v_0.

Proof:

Let us assume that u_1 is solution of:


Lu=f

Then


\forall v,\ <v,Lu_1>=<v,f>

thus


\forall v,\ <L^*v,u_1>=<v,f>

If there exists a function v_0 such that L^*v_0=0, then previous equation implies:


<v_0,f>=0

Example:

In the case of the sliding string, function v_0=1 over [0,L] is solution of the homogeneous adjoint problem. The previous compatibility equation is then: f should be of average zero:


<v_0,f>=\int_0^lf dx.

For the clamped string the homogeneous adjoint problem has no non zero solution.


  1. The reader is supposed to have a good knowledge of the solving of ordinary differential equations. A good reference on this subject is ([ma:equad:Arnold83]).