Answer of Assignment
Answer of Assignment
2 Algorithm for unsteady state heat conduction equation: Start Divide the domain with number of control volumes Discretized general form of the governing equation Initialize T Write the discretized equation for all grid points Set Told = T Step forward in time, t = t + t Calculate source terms
NO
Problem 1: Grid Details: Number of grid points taking x = 0.02 m is 100 Middle length temperature reaches 100 C at 0.1776 hr a. Explicit Scheme:
Legends: *- Analytical solution and solid line- Numerical solution Figure: Temperature profile using explicit scheme b. Crank Nicolson scheme (Semi implicit scheme)
Although Crank Nicolson is unconditionally stable in mathematical sense, it may not guarantee realistic solution. To obtain realistic solution,
c. Implicit Scheme
Initial condition: At t = 0; T(x,0) = -1 C Boundary condition: a. At t > 0; T(x = 0, t) = -15 C and b. At t > 0; T(x = L, t) = 20 C Where L = total length of the composite slab
Time step: As the problem is solved using implicit scheme, hence the time step will not influence the result.
Figure: Temperature distribution in the slab using an implicit approach till steady-state is reached
T = T(x,y,t) Initial condition: At t = 0; T(x,0) = 100 C Boundary condition: A. B. C. D. At t > 0; T(0,y,t) = 100 C and At t > 0; T(L,y,t) = 100 C At t > 0; T(x,0,t) = 200 C and At t > 0; T(x,H,t) = 150 C Where L and H are the length and height of the plate
Number of time step Figure: Time step-Temperature curve [Number of time step will be replaced by (Number of time step t)] b. 0.5 s and 0.01 m Temperature (C)
Number of time step Figure: Time step-Temperature curve [Number of time step will be replaced by (Number of time step t)]
Number of time step Figure: Time step-Temperature curve [Number of time step will be replaced by (Number of time step t)] d. 1 s and 0.005 m Temperature (C)
Number of time step Figure: Time step-Temperature curve [Number of time step will be replaced by (Number of time step t)]
Comparison of temperature (C) at given grid points grid size 0.01 0.3 0.5 1 grid size 0.005 1 x = 0.15, y = 0.2 60.23 61.27 62.03 x = 0.1, y = 0.15 69.78 70.39 70.84 x = 0.2, y = 0.05 45.1 45.53 45.84
61.92
70.77
45.79