|
|
| Line 1: |
Line 1: |
| In [[mathematics]], and more specifically in [[partial differential equations]], '''Duhamel's principle''' is a general method for obtaining solutions to [[homogeneous differential equation|inhomogeneous]] linear evolution equations like the [[heat equation]], [[wave equation]], and [[vibrating plate]] equation. It is named after [[Jean-Marie Duhamel]] who first applied the principle to the inhomogeneous heat equation that models, for instance, the distribution of heat in a thin plate which is heated from beneath. For linear evolution equations without spatial dependency, such as a [[harmonic oscillator]], Duhamel's principle reduces to the method of [[variation of parameters]] technique for solving linear inhomogeneous [[ordinary differential equations]].<ref> Fritz John, "Partial Differential Equations', New York, Springer-Verlag, 1982, 4th ed., 0387906096 </ref>
| | Hello. Allow me introduce the author. Her title is Refugia Shryock. To perform baseball is the hobby he will never stop doing. For a while I've been in South Dakota and my mothers and fathers reside close by. In her expert life she is a payroll clerk but she's usually needed her own business.<br><br>my web blog; [http://premium.asslikethat.com/blog/9567 home std test kit] |
| | |
| The philosophy underlying Duhamel's principle is that it is possible to go from solutions of the [[Cauchy problem]] (or initial value problem) to solutions of the inhomogeneous problem. Consider, for instance, the example of the heat equation modeling the distribution of heat energy ''u'' in '''R'''<sup>''n''</sup>. The initial value problem is
| |
| :<math>
| |
| \begin{cases}
| |
| u_t(x,t) - \Delta u(x,t) = 0 &(x,t)\in \mathbf{R}^n\times (0,\infty)\\
| |
| u(x,0) = g(x) & x\in \mathbf{R}^n
| |
| \end{cases}
| |
| </math>
| |
| where ''g'' is the initial heat distribution. By contrast, the inhomogeneous problem for the heat equation is
| |
| :<math>
| |
| \begin{cases}
| |
| u_t(x,t) -\Delta u(x,t) = f(x,t) &(x,t)\in \mathbf{R}^n\times (0,\infty)\\
| |
| u(x,0) = 0 & x\in \mathbf{R}^n
| |
| \end{cases}
| |
| </math>
| |
| corresponds to adding an external heat energy ''ƒ''(''x'',''t'')''dt'' at each point. Intuitively, one can think of the inhomogeneous problem as a set of homogeneous problems each starting afresh at a different time slice ''t'' = ''t''<sub>0</sub>. By linearity, one can add up (integrate) the resulting solutions through time ''t''<sub>0</sub> and obtain the solution for the inhomogeneous problem. This is the essence of Duhamel's principle.
| |
| | |
| ==General considerations==
| |
| Formally, consider a linear inhomogeneous evolution equation for a function
| |
| :<math>u:D\times(0,\infty)\to \mathbf{R}</math>
| |
| with spatial domain ''D'' in '''R'''<sup>''n''</sup>, of the form
| |
| :<math>
| |
| \begin{cases}
| |
| u_t(x,t) -Lu(x,t) = f(x,t) &(x,t)\in D\times (0,\infty)\\
| |
| u|_{\partial D} = 0 &\\
| |
| u(x,0) = 0 & x\in D,
| |
| \end{cases}</math>
| |
| where ''L'' is a linear differential operator that involves no time derivatives.
| |
| | |
| Duhamel's principle is, formally, that the solution to this problem is
| |
| :<math>u(x,t) = \int_0^t (P^sf)(x,t)\,ds</math>
| |
| where ''P''<sup>''s''</sup>''ƒ'' is the solution of the problem
| |
| :<math>
| |
| \begin{cases}
| |
| u_t - Lu = 0 & (x,t)\in D\times (s,\infty)\\
| |
| u|_{\partial D} = 0 &\\
| |
| u(x,s) = f(x,s) & x\in D.
| |
| \end{cases}
| |
| </math>
| |
| | |
| Duhamel's principle also holds for linear systems (with vector-valued functions ''u''), and this in turn furnishes a generalization to higher ''t'' derivatives, such as those appearing in the wave equation (see below). Validity of the principle depends on being able to solve the homogeneous problem in an appropriate function space and that the solution should exhibit reasonable dependence on parameters so that the integral is well-defined. Precise analytic conditions on ''u'' and ''f'' depend on the particular application.
| |
| | |
| ==Examples==
| |
| ===Wave equation===
| |
| Given the inhomogeneous wave equation:
| |
| | |
| :<math>u_{tt}-c^2u_{xx}=f(x,t)\,</math>
| |
| | |
| with initial conditions
| |
| | |
| :<math>u(x,0)=u_t(x,0)=0.\,</math>
| |
| | |
| A solution is
| |
| | |
| :<math>u(x,t) = \frac{1}{2c}\int_0^t\int_{x-c(t-s)}^{x+c(t-s)} f(\xi,s)\,d\xi\,ds.\,</math>
| |
| | |
| ===Constant-coefficient linear ODE===
| |
| Duhamel's principle is the result that the solution to an inhomogeneous, linear, partial differential equation can be solved by first finding the solution for a step input, and then superposing using [[Duhamel's integral]].
| |
| Suppose we have a constant coefficient, m<sup>th</sup> order inhomogeneous [[ordinary differential equation]].
| |
| | |
| :<math> P(\partial_t)u(t) = F(t) \,</math>
| |
| | |
| :<math> \partial_t^j u(0) = 0, \; 0 \leq j \leq m-1 </math>
| |
| | |
| where
| |
| | |
| :<math> P(\partial_t) := a_m \partial_t^m + \cdots + a_1 \partial_t + a_0,\; a_m \neq 0. </math>
| |
| | |
| We can reduce this to the solution of a homogeneous ODE using the following method. All steps are done formally, ignoring necessary requirements for the solution to be well defined.
| |
| | |
| First let ''G'' solve
| |
| | |
| :<math> P(\partial_t)G = 0, \; \partial^j_t G(0) = 0, \quad 0\leq j \leq m-2, \; \partial_t^{m-1} G(0) = 1/a_m. </math>
| |
| | |
| Define <math> H = G \chi_{[0,\infty)} </math>, with <math>\chi_{[0,\infty)}</math> being the [[Indicator function|characteristic function]] of the interval <math>[0,\infty)</math>. Then we have
| |
| | |
| :<math> P(\partial_t) H = \delta </math>
| |
| | |
| in the sense of [[Distribution (mathematics)|distributions]]. Therefore
| |
| | |
| :<math> u(t) = (H \ast F)(t) </math> | |
| | |
| :<math> = \int_0^\infty G(\tau)F(t-\tau)\,d\tau </math>
| |
| | |
| :<math> = \int_{-\infty}^t G(t-\tau)F(\tau)\, d\tau </math>
| |
| | |
| solves the ODE.
| |
| | |
| ===Constant-coefficient linear PDE===
| |
| More generally, suppose we have a constant coefficient inhomogeneous [[partial differential equation]]
| |
| | |
| :<math> P(\partial_t,D_x)u(t,x) = F(t,x) \,</math>
| |
| | |
| where
| |
| | |
| :<math> D_x = \frac{1}{i} \frac{\partial}{\partial x}. \,</math>
| |
| | |
| We can reduce this to the solution of a homogeneous ODE using the following method. All steps are done formally, ignoring necessary requirements for the solution to be well defined.
| |
| | |
| First, taking the [[Fourier transform]] in ''x'' we have
| |
| | |
| :<math> P(\partial_t,\xi)\hat u(t,\xi) = \hat F(t,\xi).</math>
| |
| | |
| Assume that <math> P(\partial_t,\xi) </math> is an m<sup>th</sup> order ODE in ''t''. Let <math>a_m </math> be the coefficient of the highest order term of <math> P(\partial_t,\xi) </math>.
| |
| Now for every <math>\xi </math> let <math>G(t,\xi) </math> solve
| |
| | |
| :<math> P(\partial_t,\xi)G(t,\xi) = 0, \; \partial^j_t G(0,\xi) = 0 \; \mbox{ for } 0\leq j \leq m-2, \; \partial_t^{m-1} G(0,\xi) = 1/a_m. </math>
| |
| | |
| Define <math>H(t,\xi) = G(t,\xi) \chi_{[0,\infty)}(t) </math>. We then have
| |
| | |
| :<math> P(\partial_t,\xi) H(t,\xi) = \delta(t) </math>
| |
| | |
| in the sense of [[Distribution (mathematics)|distributions]]. Therefore
| |
| | |
| :<math> \hat u(t,\xi) = (H(\cdot,\xi) \ast \hat F(\cdot,\xi))(t) </math>
| |
| | |
| :<math> = \int_0^\infty G(\tau,\xi)F(t-\tau,\xi)\,d\tau </math>
| |
| | |
| :<math> = \int_{-\infty}^t G(t-\tau,\xi)F(\tau,\xi)\, d\tau </math>
| |
| | |
| solves the PDE (after transforming back to ''x'').
| |
| | |
| ==See also==
| |
| *[[Retarded potential]]
| |
| | |
| == References==
| |
| {{reflist}}
| |
| | |
| [[Category:Wave mechanics]]
| |
| [[Category:Partial differential equations]]
| |
| | |
| [[Category:Mathematical principles]]
| |