Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

This is with regards to the first section of's_principle. I want to see if I am understanding this.

Basically the inhomogeneous equation says that heat is being added at a rate of $f(x,t)$, so at each point $x$, $f(x,t)\,dt$ of heat is being added to the system.

And you can think of this as a set of systems $u_t(x,t)-\Delta u(x,t)=0$, $u(x,t_0)=f(x,t_0)$ for all $t_0\in(0,\infty)$, which says the initial heat distribution is $f(x,t)$, and no heat is being added, and if we integrate these solutions over $t_0$ we get the inhomogeneous solution...

It's the third paragraph I am a bit uncomfortable with. Can anyone offer some help?

share|cite|improve this question
up vote 14 down vote accepted

Don't hurt your brain too much in trying to understand the physical interpretation. Start from the mathematics! Let me do the ODE version first:

Let $X: \mathbb{R}\to \mathbb{R}^N$ be a vector valued function of one (time) variable. Let $A:\mathbb{R}^N \to\mathbb{R}^N$ be a linear transformation that is independent of the variable $t$. Let $R_t$ denote the solution operator: that is $ X(t) = R_tX_0$ where $X_0$ is a fixed vector solves the equation $$ \dot{X} = AX $$ with initial data $X(0) = X_0$. Now let $Y(t)$ be another vector valued function. Consider the expression

$$ X(t) = \int_0^t R_{t-s}Y(s) \mathrm{d}s $$

Computing it explicitly we have that

$$ \left(\frac{d}{dt} - A\right)X(t) = \int_0^t\left(\frac{d}{dt} - A\right)\left[R_{t-s}Y(s)\right]~ \mathrm{d}s + R_{t-s}Y(s)\Big|_{s = t} $$

The second term comes from the fundamental theorem of calculus when the $\frac{d}{dt}$ derivative hits the integral sign.

By definition the term under the integral sign evaluates to 0, since $Z(t) = R_{t-s}Y(s)$ solves the homogeneous equation with initial data $Z(s) = Y(s)$. So we are left with

$$ (\frac{d}{dt} - A)X(t) = R_0Y(s) = Y(s) $$

i.e. that $X(t)$ solves the inhomogeneous problem.

For the PDE version, you just replace $\mathbb{R}^N$ with a Banach or Hilbert space, and the computation formally carries through in exactly the same way.

Coming back to the physical interpretation: you see that Wikipedia is being a bit imprecise: it is not the adding up of solutions that matters really; what drives Duhamel's principle is the fact that you are adding only a limited number of solutions.

Let me clarify: imagine you have a family of initial data $Y(s)$. And you write down the expression

$$ X(t) = \int_{a}^{b} R_{t-s}Y(s) \mathrm{d}s $$

which would be what we do if we were to just add (integrate) the contributions from all of the linear waves coming from the "inhomogeneity", by linearity it is clear that $X(t)$ will still solve the homogeneous equation, since it is a fixed (as in the limits of the integral) sum of may solutions.

The magic of Duhamel's principle is in that the upper-limit of the integral is time! That is, we defined

$$ X(t) = \int_0^t R_{t-s}Y(s) \mathrm{d}s $$

As you see from the derivation above, it is this upper-limit which, when acted on using the fundamental theorem of calculus, give you the inhomogeneous term. So what is the physical interpretation then? That the upper-limits also changes represents the fact that the solution at $X(t + \Delta t)$ consists of the forward time evolution of the solution at $X(t)$ plus a new contribution from the data in $(t,t+\Delta t)$ which was not included in the computation in $X(t)$. This "adding a new contribution" is precisely what we imagine the inhomogeneous term as, that is, a source term!

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.