Suppose that we want to find the temperature
in the thin (one dimensional) rod of finite length L extending form
to
extending from
to
. It is governed by the equation
 |
(3.1) |
where
is thermal diffusing (Heat Conduction Coefficient) of the rod. Suppose that end points of the rod are kept at temperature zero. ( e.g. by damping against a large ice block). Thus we get the homogeneous boundary conditions
 |
(3.2) |
Also, Assume that across the length the rod is heated at time
as given by the initial condition
 |
(3.3) |
In order to make (
) consistent with (
), we assume the compatibility condition
 |
(3.4) |
Question: Any restriction on
(If
is piece-wise continuous function, it can be showed that a formal series solution always satisfy the BC's and it is unique solution of the BVP. For part see,chap.6,Churchill and Brawn,Fainter series and solution of BVP.,3rd Ed.)
Solution by method of separation of variables:
Suppose that solution of model problem (
),(
)and (
) can be expressed as
 |
(3.5) |
in which variables
and
are "separated". Substituting (
) in PDE (
), we get
 |
(3.6) |
where ' stands for derivative w.r.t respective variables.
 |
(3.7) |
LHS of (
) is function of
only ,where as RHS. is a function of t only. For a fixed t on RHS and LHS must remain constant as
varies and it should be constant say
 |
(3.8) |
This
is also known as a separation constant. As a result we get two equations
 |
(3.9) |
 |
(3.10) |
in (
) satisfy the BC's (
) iff
 |
(3.11) |
 |
(3.12) |
Since
is not a trivial solution.
Thus the solution
is the solution of BVP.
We can interpret the above system as eigenvalue problem. A non-trivial solution of the system is called eigen-function of the problem with an eigenvalue
. Depending on the sign of
, we have different terms of solution of the above problem
- If
- If
- If
Where
and
are arbitrary constants. Here we have implicitly assumed that
is purely real Why? Answer is hidden in Strum-Liouville theory. In fact it can be shown that all eigenvalues of the problem are real constants.
If
is a strictly positive function
has a unique root at
Thus for
, we have
This eigenvalue problem can not have positive eigenvalues.
If
is not allowed.
trivial solution.
If
For the non-trivial solution
(why not
)
and the corresponding eigen-function is
(They are unique up-to a multiplicative constant.)
Now, Let us look at ODE (
)
 |
(3.14) |
as
: consistent with physical expectation c
choice of sign of
was right. Thus we obtain sequence of separated solutions
By the principle of the superposition of solution,the general solution is given by
But we shall need to find the unknown coefficients
They can be found by using 'initial condition'
 |
(3.17) |
Multiply both sides by
and integrate over 0 to L
 |
(3.18) |
 |
(3.19) |
Using the orthogonality of sin functions
For a function f(x)periodic on interval
instead of
a change of variable gives
These
are the Fourier sine coefficients putting this in eq.( ).We get Fourier expansion of f(x) on [0,L].These Fourier coefficients and the Fourier expansion of f(x) on [0,L] is uniquely determined. Thus the formal series solution of the model problem (1)-(3)
 |
(3.20) |
where
are the Fourier coefficients of
on
- Notice the power all method.For a given initial condition for f, we only have to compute its Fourier sine coefficients in order to obtain the explicit solution.
- The series solution Eq. (
) due to converges rapidly(unless
is very small), due to exponentially decrease solution. Therefore infinite sum can be truncated by finite sum taking sufficiently large
This is used for comparison with numerical solution.
(C.Henny Edwars,Darrd E.Penny ,Differential equations and BVP.s computing and modeling 3rd Ed.)
Solution The BVP for the temperature
of rod is
For this Dirichlet problem, the solution is given by
where
are the Fourier sine coefficients pf the Fourier sine expansion of
on
given by
To compute
This can be written as depend on
as
Thus the rod's temperature is given by
or
Given
and
We want
at
cm,
- With
for iron, we get
- With
for concrete, we get
for
times
and for
times
which implies that concrete is very effective insulator.