Duhamel's formula: Difference between revisions

From DispersiveWiki
Jump to navigationJump to search
No edit summary
 
m (Minor correction)
 
(5 intermediate revisions by the same user not shown)
Line 5: Line 5:
for some spatial operator ''L'', is given by
for some spatial operator ''L'', is given by
<center><math> u(t) = e^{tL} u_0 + \int_0^t e^{(t-t')L} F(t')\ dt',</math></center>
<center><math> u(t) = e^{tL} u_0 + \int_0^t e^{(t-t')L} F(t')\ dt',</math></center>
provided that ''L'' has enough of a functional calculus, and ''u'', ''u_0'', ''F'' have enough regularity, to justify all computations.  (If ''L'' is constant coefficient, then the Fourier transform can usually be used to justify everything so long as one works in the category of tempered distributions.)  Note that the case ''L=0'' is simply the fundamental theorem of calculus, indeed one can view Duhamel's formula as the fundamental theorem of calculus twisted (conjugated) by the free propagator <math>e^{tL}</math>.
provided that ''L'' has enough of a functional calculus, and <math> u </math>, <math> u_0 </math>, <math> F </math> have enough regularity, to justify all computations.  (If ''L'' is constant coefficient, then the Fourier transform can usually be used to justify everything so long as one works in the category of tempered distributions.)  Note that the case ''L=0'' is simply the fundamental theorem of calculus, indeed one can view Duhamel's formula as the fundamental theorem of calculus twisted (conjugated) by the free propagator <math>e^{tL}</math>.


For equations which are second order in time, the formula is slightly more complicated.  For instance, the solution to the inhomogeneous initial value problem
For equations which are second order in time, the formula is slightly more complicated.  For instance, the solution to the inhomogeneous initial value problem
Line 11: Line 11:
is given (formally, at least) by
is given (formally, at least) by
<center><math> u(t) = \cos(t\sqrt{L}) u_0 + \frac{\sin(t\sqrt{L})}{\sqrt{L}} u_1 + \int_0^t \frac{\sin((t-t')\sqrt{L})}{\sqrt{L}} F(t')\ dt'.</math></center>
<center><math> u(t) = \cos(t\sqrt{L}) u_0 + \frac{\sin(t\sqrt{L})}{\sqrt{L}} u_1 + \int_0^t \frac{\sin((t-t')\sqrt{L})}{\sqrt{L}} F(t')\ dt'.</math></center>
Anyhow, we note that in this case the solution can be cast in the standard form. So, let us introduce the vectors
<center><math>{\underline y}=\left[\begin{matrix} u \\ w \end{matrix}\right]</math></center>
and
<center><math>{\underline \Phi}=\left[\begin{matrix} 0 \\ F(t) \end{matrix}\right]</math></center>
with the matrix
<center><math>{\hat M}=\begin{bmatrix} 0 & 1 \\
                                L & 0 \end{bmatrix}</math></center>.
We can write the second order equation as
<center><math>{\underline y}_{t}-{\hat M}{\underline y}={\underline \Phi(t)}</math></center>
and write down the solution as expected in the original Duhamel's formula, that is
<center><math> {\underline y}(t) = e^{t\hat M} {\underline y}_0 + \int_0^t e^{(t-t')\hat M} {\underline \Phi}(t')\ dt',</math></center>.
Useful applications of this approach can be found for systems having a Hamiltonian flow.


[[Category:methods]]
[[Category:methods]]

Latest revision as of 14:26, 17 September 2007


Duhamel's formula expresses the solution to a general inhomogeneous linear equation as a superposition of free solutions arising from both the initial data and the forcing term. For instance, the solution to the inhomogeneous initial value problem

for some spatial operator L, is given by

provided that L has enough of a functional calculus, and , , have enough regularity, to justify all computations. (If L is constant coefficient, then the Fourier transform can usually be used to justify everything so long as one works in the category of tempered distributions.) Note that the case L=0 is simply the fundamental theorem of calculus, indeed one can view Duhamel's formula as the fundamental theorem of calculus twisted (conjugated) by the free propagator .

For equations which are second order in time, the formula is slightly more complicated. For instance, the solution to the inhomogeneous initial value problem

is given (formally, at least) by

Anyhow, we note that in this case the solution can be cast in the standard form. So, let us introduce the vectors

and

with the matrix

.

We can write the second order equation as

and write down the solution as expected in the original Duhamel's formula, that is

.

Useful applications of this approach can be found for systems having a Hamiltonian flow.