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.