Cauchy boundary condition

From Wikipedia, the free encyclopedia
Jump to: navigation, search

In mathematics, a Cauchy boundary condition /kˈʃ/ augments an ordinary differential equation or a partial differential equation to completely determine the solution. Thereby, both function value and normal derivative are specified on the boundary of the domain. This corresponds to imposing both a Dirichlet and a Neumann boundary condition. It is named after the prolific 19th century French mathematical analyst Augustin Louis Cauchy.

Cauchy boundary conditions can be understood from the theory of second order ordinary differential equations, where, in order to have a particular solution, one must specify the value of the function and the value of the derivative at a given initial or boundary point, i.e.,

y(a)=\alpha \ ,

and

y'(a)=\beta  \ .

where  a \ is a boundary or initial point.

Cauchy boundary conditions are the generalization of these type of conditions. Let us first recall a simplified form for writing partial derivatives.

\begin{align}
u_x    &= {\part u \over \part x} \\
u_{xy} &= {\part^2 u \over \part y\, \part x}
\end{align}

and let us now define a simple, second order, partial differential equation:

\psi_{xx} + \psi_{yy}= \psi(x,y) \

We have a two dimensional domain whose boundary is a boundary line, which in turn can be described by the following parametric equations

\begin{align}
x &= \xi (s) \\
y &= \eta (s)
\end{align}

hence, in a similar manner as for second-order ordinary differential equations, we now need to know the value of the function at the boundary, as well as its normal derivative, in order to solve the partial differential equation, that is to say, both

\psi (s) \

and

\frac{d\psi}{dn}(s)=\mathbf{n}\cdot\nabla\psi \

are specified at each point on the boundary of the domain of the given partial differential equation (PDE), where \nabla\psi(s) \, is the gradient of the function. It is sometimes said that Cauchy boundary conditions are a weighted average of imposing Dirichlet boundary conditions and Neumann boundary conditions. This should not be confused with statistical objects such as the weighted mean, the weighted geometric mean, or the weighted harmonic mean, since no such formulas are used upon imposing Cauchy boundary conditions. Rather, the term weighted average means that while analyzing a given boundary value problem, one should bear in mind all available information for its well-posedness and subsequent successful solution.

Since the parameter s \ is usually time, Cauchy conditions can also be called initial value conditions or initial value data or simply Cauchy data.

Notice that although Cauchy boundary conditions imply having both Dirichlet and Neumann boundary conditions, this is not the same at all as having a Robin or impedance boundary condition. A mixture of Dirichlet and Neumann boundary conditions is given by

\alpha (s)\psi (s)+ \beta (s) \frac{d\psi }{dn}(s)=f(s) \

where \alpha (s) \ , \beta (s) \ , and  f(s) \ are understood to be given on the boundary (this contrasts to the term mixed boundary conditions, which is generally taken to mean boundary conditions of different types on different subsets of the boundary). In this case the function and its derivative must fulfill a condition within the same equation for the boundary condition.

Example

Let us define the heat equation in two spatial dimensions as follows

u_t = \alpha \nabla^2 u \

where \alpha \ is a material-specific constant called thermal diffusivity

and suppose that such equation is applied over the region G \ , which is the upper semidisk of radius a \ centered at the origin. Suppose that the temperature is held at zero on the curved portion of the boundary, while the straight portion of the boundary is insulated, i.e., we define the Cauchy boundary conditions as

u=0 \ \forall (x,y) \in r=a, 0\leq \theta \leq \pi \

and

u_y = 0, y = 0 \

We can use separation of variables by considering the function as composed by the product of the spatial and the temporal part

u(x,y,t)= \phi (x,y) \psi (t)\

applying such product to the original equation we obtain

\phi (x,y) \psi ' (t)= \alpha \phi '' (x,y) \psi (t) \

whence

 \frac{\psi '(t)}{\alpha \psi (t)} = \frac{\phi '' (x,y)}{\phi (x,y)}

Since the left hand side (l.h.s.) depends only on t \ , and the right hand side (r.h.s.) depends only on (x,y) \ , we conclude that both should be equal to the same constant

\frac {\psi '(t)}{\alpha \psi (t)}= - \lambda = \frac {\phi '' (x,y)}{\phi (x,y)}

Thus we are led to two equations: the first in the spatial variables

\phi_{xx}+\phi_{yy}+\lambda \phi (x,y)=0 \

and a second equation in the t \ variable,

\psi '(t) +\lambda \alpha \psi (t)=0 \

Once we impose the boundary conditions, the solution of the temporal ODE is

\psi (t) =A e^{-\lambda \alpha t}\

where A is a constant which could be defined upon the initial conditions. The spatial part can be solved again by separation of variables, substituting \phi (x,y) = X(x)Y(y) \ into the PDE and dividing by  X(x) Y(y) \ from which we obtain (after reorganizing terms)

\frac {Y''}{Y}+\lambda =-\frac {X''}{X}

since the l.h.s. depends only on y and r.h.s. only depends on x \ , both sides must equal a constant, say \mu \ ,

\frac {Y''}{Y}+ \lambda =- \frac {X''}{X} = \mu

so we obtain a pair of ODE's upon which we can impose the boundary conditions that we defined

See also

References

  • Cooper, Jeffery M. "Introduction to Partial Differential Equations with MATLAB". ISBN 0-8176-3967-5

External links








Creative Commons License