07.12.2012 Views

Laplace transform isotherm .pdf - University of Hertfordshire ...

Laplace transform isotherm .pdf - University of Hertfordshire ...

Laplace transform isotherm .pdf - University of Hertfordshire ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

simplify the algebra we write this as −µ 2 α. The solutions are then<br />

where T0 is a constant and<br />

T = T0e −µ2 αt<br />

X = Asin µx + B cos µx<br />

with A and B constants. A solution <strong>of</strong> equation(2.3) is then<br />

u = T0(Asin µx + B cos µx)exp � −µ 2 αt �<br />

(2.5)<br />

and since the heat equation is a linear equation, the most general solution<br />

is obtained by summing solutions <strong>of</strong> equation (2.5) type to get<br />

u =<br />

∞�<br />

T0,m(Am sin µmx + Bm cos µmx)exp � −µ 2 mαt� m=1<br />

The constants Am,Bm,T0,m and µm are determined from the initial and<br />

boundary conditions and again, examples are to be found in Crank (1979).<br />

Another analytic method is the use <strong>of</strong> similarity solutions.<br />

We define the dimensionless variable<br />

η = x<br />

√ αt<br />

and then look for solutions <strong>of</strong> the form<br />

u(x,t) = t p g (η) (2.6)<br />

where the number p and the function g(η) are to be found. Substituting<br />

equation (2.6) into equation (2.3) we find<br />

which implies<br />

�<br />

p−1<br />

t pg − η<br />

2 g′ − g ′′�<br />

= 0<br />

g ′′ + η<br />

2 g′ = pg (2.7)<br />

subject to appropriate boundary conditions using equation (2.6). This is<br />

difficult to solve in closed form for arbitrary values <strong>of</strong> p.<br />

10

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!