11.10.2013 Views

Conference, Proceedings

Conference, Proceedings

Conference, Proceedings

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.

1a. Three‐dimensional temperature field inside a block<br />

The border conditions mentioned above are<br />

( 1 = 2, 2, 3, ) = ( 1 =− 2, 2, 3, ) = ( 2 = 2 2, 1, 3,<br />

) =<br />

( ) ( ) ( )<br />

T x l x x t T x l x x t T x l x x t<br />

T x2 =− l2 2, x1, x3, t = T x3 = l3 2, x1, x2, t = T x3 =− l3 2, x1, x2, t = Ts,<br />

t > 0 (1.1)<br />

and the initial condition is<br />

( )<br />

T x1, x2, x3, t = 0 = T0 , -li 2 < xi < li<br />

2<br />

,<br />

( 1, 2, 3,<br />

)<br />

Solution<br />

T x x x t of 3D heat conduction equation<br />

∂T 2 2<br />

−a ∇ T =<br />

i = 1, 2, 3<br />

, t > 0 (1.2)<br />

0<br />

∂ t<br />

(1.3)<br />

by the method of variable separation assuming it obeys boundary (1.1), and initial (1.2)<br />

conditions is given as [1]<br />

( )<br />

T x1, x2, x3, t −Ts<br />

=<br />

T −T<br />

0<br />

s<br />

(<br />

n1+ n2+ n3+<br />

)<br />

( 2 −1)( 2 −1)( 2 −1)<br />

( n ) ( n ) ( n )<br />

3<br />

3 2 2 2<br />

∞<br />

⎛ 4 ⎞<br />

−1 ⎡ ⎛<br />

2 2 2 1−1 2 2 −1 2 3−1<br />

⎞ ⎤<br />

= ⎜ ⎟ ∑<br />

exp ⎢− π a ⎜ + + ⎟t⎥×<br />

2 2 2<br />

⎝π⎠ n1, n2, n3=<br />

1 n1 n2 n3 ⎢ ⎜ l1 l2 l ⎟<br />

⎝ 3<br />

⎣ ⎠ ⎥<br />

⎦<br />

x x<br />

cos 2 1 cos 2 1 cos 2 1<br />

1 2<br />

3<br />

( n − ) π ( n − ) π ( n − ) π<br />

1 2 3<br />

l1 l2 l3<br />

Threefold infinite sum contains terms (in the second line behind the sum) which with increase<br />

n<br />

of trinity natural numbers 1, n2, n3as well as with increase of the time t decrease. Space<br />

distribution of temperature inside the block at the instant of time t is determined by the product<br />

x, 1,2,3<br />

of three cosine functions. This product represents even function of space variables i i =<br />

.<br />

We see that each term in infinite sum exponentially decreases with characteristic time equal to<br />

τ<br />

i<br />

i = ⎨ ⎬<br />

aπ( 2ni−1) 2<br />

⎧⎪ l ⎫⎪<br />

1<br />

⎛ li<br />

⎞<br />

τ =<br />

⎪⎩ ⎪⎭<br />

2n1 and amplitude decrease is proportional to i −<br />

⎜ ⎟<br />

. The time ⎝2a⎠ is called the relaxation time. Here 2<br />

a = λ cρ<br />

is the diffusivity, λ ‐thermal conductivity and cρ is<br />

the heat capacity density of a body.<br />

Because, the centre of the sample coincides with coordinate origin the time development of<br />

cos( γ i )<br />

temperature at that point is given by the same formula (1.4) in which last three<br />

functions are equal to one.<br />

l<br />

In a case of a cube 1 = l2 = l3 = l<br />

at the center of a cube triple infinite sum can by reduce to one<br />

sum. The time development temperature at that place then is given by the relation<br />

x<br />

2<br />

(1.4)<br />

27

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

Saved successfully!

Ooh no, something went wrong!