ASSIGNMENT-3

Advance Heat and Mass Transfer

IMRAN SAJID SHAHID 12SP-MS-ME-002

Problem#1 Solve by analytical method for temperature distribution in given rectangular plate subject to given boundary conditions shown in figure-1?

200C

(1.0, 1.5)

100C

100C

(0, 0)

X 200C Figure-1

Solution: Assumption: i) 2D problem; ii) Steady State Problem; iii) No internal heat generation Mathematical Model: To find the steady state temperature u(x,y) in a rectangular plate whose edges are subject to above mentioned boundary conditions, as shown in figure. The boundary value problem is: ( ) ( ) ( ) ( ( ) ) ( ( ) )

In order to solve 2D Laplaces Heat Equation subject to inhomogenous boundary conditions, we first develop general solution for given problem. Consider the figure-2 to solve 2D Laplaces Equation for rectangular plate subject to the Dirichlet boundary conditions. [1]
Y g(x) (a,b)

F(y)

G(y)

(0, 0)

f(x)

Figure-2: General 2D Laplaces Heat Problem The method of separation of variables is not applicable to a Dirichlet problem when the boundary conditions on all four sides of the rectangular are inhomogeneous. To get around this difficulty, we break the problem ( ) ( ) ( ) ( ( ) ) ( ) ( ) ( ( ) ) ( ) ( )

into two problems, each of which has homogeneous boundary conditions on parallel boundaries.

Problem-1

Problem-2

( (

) ) ( )

( (

) ) ( )

( (

) )

( )

( (

) )

( )

Suppose u1 and u2 are the solutions of Problems 1 and 2, respectively. If we define ( ) ( ) ( )

u(x,y) satisfies all boundary conditions of the problems. By solving Problem-1 and Problem-2 and adding their solutions, we have solved the original problem. This additive property of solutions is known as the superposition principle. See Figure-3.

Figure-3: Solution u = Solution u1 of Problem-1 + Solution u2 of Problem-2 Consider one by one each problem and determine the solution of each Boundary Value Problem (BVP). Problem-1 ( ) ( ) ( ) Suppose; ( ) ( ) ( ) ( ( ) ) ( ) ( ( ) ) ( )

Since the boundary conditions at y=0 and y=b are functions of x we choose to separate Laplaces equation as

So that we have system of ODEs:

And the solution of ODEs is; ( ) ( ) Now ( ) ( ) Implies gives A = 0 and or

Implies that, ( ) ( ) Thus the proposed solution ( And ( ) ( ) ( ) ( ) ) becomes, ( ) ( ) ( )

Determine the At ( )

Constants by imposing remaining two boundary conditions. ( ) we have, ( )

And by using Fourier Expansion we have, ( At ) ( ) ( ) we have, ( )

( )

Indicates that the entire expression in the parenthesis is given by: We can now solve for : Implies that; ( ) ( ( ) ) given in (12) and ( ) ( )

A solution to the given boundary-value problem consist of series (11) with coefficients (13), respectively. Problem-2 ( ( ( ) ) ) Suppose; ( ) ( ) ( ) ( ( ) ) ( ) ( ( ) ) ( )

Since the boundary conditions at y=0 and y=b are functions of x we choose to separate Laplaces equation as

So that we have system of ODEs:

And the solution of ODEs is; ( ) ( ) Now ( ) gives C = 0 and

( )

Implies

or

Implies that, ( ) ( ) Thus the proposed solution u2(x,y) becomes, ( And ( ) ( ) ( ) ) ( ) ( ) ( )

Determine the At ( )

Constants by imposing remaining two boundary conditions. ( ) we have, ( )

And consequently by using Fourier Expansion we have, ( At ) ( ) ( ) we have, ( ) ( ) ( )

Indicates that the entire expression in the parenthesis is given by: We can now solve for : Implies that; ( ) ( )

( )

) given in (19) and

A solution to the given boundary-value problem consist of series (18) with coefficients (20), respectively. Summarize the results; Problem-1 Problem-2

( (

) ) ( )

( (

) ) ( )

( (

) )

( )

( (

) )

( )

Solution u1 of Problem-1
( ) ( ) ( ) (

Solution u2 of Problem-2
)

Where, ( ) ( ( ) )

Where, ( ) ( ( ) )

Now, one can have solution (

) after substituting the solutions of both problems, ( ) ( ) ( )

In the given problem statement we identify ( ) Problem-1 takes the form: ( ) ( ) ( ) ( )

Subject to; ( ( ) ) ( ( ) ) ( ( ) )

Previously we have discovered the solution u1 of Problem-1, then; ( ) ( ) ( )

And Constants takes the form, Implies that, ( ) [ ( ) ] ( ) [ ]

And , ( [ ( ) ) ] [ ) ] ) ( ( )

) ]

Substitute the value of constants in eq(24); ( ) ( [ ( ) ] [ ( ) ] ( ) )

) ]

) ]

) ]

)]

Problem-2 takes the form: ( )

Subject to; ( ) ( ) ( )

Previously we have discovered the solution u2 of Problem-2, then; ( ) ( ) ( )

And Constants takes the form, ( )

Implies that; ( ) [ ( ) ]

And , ( ( ) [ ( ) ] ) ( )

Substitute the value of constants in eq(31); ( ) ( [ ( ) ] [ ( ) ] ( ) )

) ]

) ]

)]

So the Solution u(x,y) in rectangular plate is obtained by substituting the values of

from eq(27) & eq(34):

) ]

)]

Problem#2 Solve by analytical method for temperature distribution in given rectangular plate subject to given boundary conditions shown in figure-4?

T1 + T2 Sin (x/W)

(1.0, 1.5)

T1

T1

(0, 0) T1

Figure-4 Where, T1 = 50OC T2=100OC WIDTH,W =1m HEIGHT,H=1.5m

Solution: Assumption: i) 2D problem; ii) Steady State Problem; iii) No internal heat generation Mathematical Model: To find the steady state temperature u(x,y) in a rectangular plate whose edges are subject to above mentioned boundary conditions, as shown in figure. The boundary value problem is: ( ( ( ) ) ) ( ( ) ) ( ( ) ) ( )

After applying the superposition principle, we can break the problem into two problems, each of which has homogeneous boundary conditions on parallel boundaries, as shown in Figure-3 Suppose u1 and u2 are the solutions of Problems 1 and 2, respectively. If we define ( ( ) ( ) ( )

) satisfies all boundary conditions of the problem.

In the given problem statement we identify ( ) Problem-1 takes the form: ( ) ( ) ( ) ( ) ( )

Subject to; ( ( ) ) ( ( ) ) ( ( ) ) ( ( ))

Previously we have discovered the general solution u1 of Problem-1, then; ( ) ( ) ( )

And Constants takes the form, ( ) [ ]

Implies that, ( ) [ ( ) ]

And , ( ( ) ) ( ( ) ).

Consider the second integral term. From trigonometry we have learnt that; In our above equation ( ( ( [ ( ) ] [ (
( ( ) ) ( ( ) )

and taking limits we have; ) ) ( ( ) ) ] ) ) ( ( ) ) ( ( ) )

) ]

Substitute the value of constants in eq(42); ( ) ( [ ( ) ] [ ( ) ] ( ) )

) ]

) ]

))

) ]

)]

Problem-2 takes the form: ( )

Subject to; ( ( ) ) ( ( ) ) ( ( ) )

Previously we have discovered the solution u2 of Problem-2, then; ( ) ( ) ( )

And Constants takes the form, ( )

Implies that; ( ) [ ( ) ]

And , ( ( ) [ ( ) ] ) ( )

Substitute the value of constants in eq(31); ( ) ( [ ( ) ] [ ( ) ] ( ) )

) ]

) ]

)]

So the Solution u(x,y) in rectangular plate is obtained by substituting the values of

from eq(45) & eq(52): )+

(
( ) ]

) ]

(
( )]

MATLAB PDETool was utilized to plot temperature plots and heat flux for both problems 1& 2. Following plots illustrate the model and results obtained from PDETool, Heat transfer module [2]. Problem-1

Figure-5: Temperature isothermal lines along with Heat flux arrows plot

Figure-6: Triangular Mesh Elements

Figure-7: PDE Equation is defines with k of steel [3]

Figure-8: 3D Temperature plot Problem-2

Figure-9: Temperature isothermal lines along with Heat flux arrows plot

Figure-10: Boundary condition on top face

Figure-11: 3D Temperature plot