In mathematics, a Cauchy (pronounced "coe-she") boundary condition imposed on an ordinary differential equation or a partial differential equation specifies both the values a solution of a differential equation is to take on the boundary of the domain and the normal derivative at the boundary. It 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 to have a particular solution one has to specify the value of the function and the value of the derivative at a given initial or boundary point, i.e.,
- <math>y(a)=\alpha \ ,</math>
and
- <math>y'(a)=\beta \ . </math>
where <math> a \ </math> 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.
- <math>u_x = {\part u \over \part x} \ </math>
- <math>u_{xy} = {\part^2 u \over \part y\, \part x} \ </math>
and let us now define a simple, second order, partial differential equation:
- <math>\psi_{xx} + \psi_{yy}= \psi(x,y) \ </math>
We have a two dimensional domain whose boundary is a boundary line, which in turn can be described by the following parametric equations
- <math>x=\xi (s) \ </math>
- <math>y=\eta (s) \ </math>
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, and its normal derivative in order to solve the partial differential equation, that is to say, both
- <math>\psi (s) \ </math>
and
- <math>\frac{d\psi}{dn}(s)=\mathbf{n}\cdot\nabla\psi \ </math>
are specified at each point on the boundary of the domain of the given partial differential equation (PDE), where <math>\nabla\psi(s) \, </math> 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 <math>s \ </math> 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 Robin or impedance boundary condition, a mixture of Dirichlet and Neumann boundary conditions are given by
- <math>\alpha (s)\psi (s)+ \beta (s) \frac{d\psi }{dn}(s)=f(s) \ </math>
here <math>\alpha (s) \ </math>, <math>\beta (s) \ </math>, and <math> f(s) \ </math> 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
- <math>u_t = k \nabla^2 u \ </math>
where <math>k \ </math> is a material-specific constant called thermal conductivity. and suppose that such equation is applied over the region <math>G \ </math>, which is the upper semidisk centered at the origin of radius <math>a \ </math>. 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
- <math>u=0 \ \forall (x,y) \in r=a, 0\leq \theta \leq \pi \ </math>
and
- <math>u_y = 0, y = 0 \ </math>
We can use separation of variables by considering the function as composed by the product of the spatial and the temporal part
- <math>u(x,y,t)= \phi (x,y) \psi (t)\ </math>
applying such product to the original equation we obtain
- <math>\phi (x,y) \psi ' (t)= k \phi (x,y) \psi (t) \ </math>
whence
- <math> \frac{\psi '(t)}{k \psi (t)} = \frac{\phi (x,y)}{\phi (x,y)} </math>
Since the left hand side (l.h.s.) depends only on <math>t \ </math>, and the right hand side (r.h.s) depends only on <math>(x,y) \ </math>, we conclude that both should be equal to the same constant
- <math>\frac {\psi '(t)}{k \psi (t)}= - \lambda = \frac {\phi (x,y)}{\phi (x,y)} </math>
Thus we are led to two equations: the first in the spatial variables
- <math>\phi_{xx}+\phi_{yy}+\lambda \phi (x,y)=0 \ </math>
and a second equation in the <math>t \ </math> variable,
- <math>\psi '(t) +\lambda k \psi (t)=0 \ </math>
Once we impose the boundary conditions, the solution of the temporal ODE is
- <math>\psi (t) =A e^{-\lambda k t}\ </math>
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 <math>\phi (x,y) = X(x)Y(y) \ </math> into the PDE and dividing by <math> X(x) Y(y) \ </math> from which we obtain (after reorganizing terms)
- <math>\frac {Y}{Y}+\lambda =-\frac {X}{X} </math>
since the l.h.s depends only on y and r.h.s only depends on <math>x \ </math>, both sides must equal a constant, say <math>\mu \ </math>,
- <math>\frac {Y}{Y}+ \lambda =- \frac {X}{X} = \mu </math>
so we obtain a pair of ODE's upon which we can impose the boundary conditions that we defined above
External links
- http://mathworld.wolfram.com/CauchyBoundaryConditions.html
- http://www.math.ohio-state.edu/~gerlach/math/BVtypset/node133.html
- http://ocw.mit.edu/NR/rdonlyres/Mathematics/18-303Fall-2004/88A7832B-A69E-48B9-A566-2F3B810164F0/0/pde3d.pdf
References
Cooper, Jeffery M. "Introduction to Partial Differential Equations with MATLAB". ISBN 0-8176-3967-5
