MM210A6
MM210A6
2. Solve the temperature distribution in the plate for the problem given in L17b, slide 15. Use
Gauss-Seidel method to solve iteratively for the temperatures until the percent relative error
falls below 𝜀𝑠 = 2%.
3. Read Section 32.1 in the textbook (Chapra and Canale, 8th edition) carefully.
The advection-diffusion equation is used to compute the distribution of concentration along the
length of a rectangular chemical reactor:
𝜕𝑐 𝜕2𝑐 𝜕𝑐
=𝐷 2−𝑈 − 𝑘𝑐
𝜕𝑡 𝜕𝑥 𝜕𝑥
1
where, 𝑐 is the concentration (𝑚𝑔/𝑚3 ), 𝑡 is the time (𝑚𝑖𝑛), 𝐷 is a diffusion coefficient
(𝑚2 /𝑚𝑖𝑛), 𝑥 is the distance along the tank’s longitudinal axis (𝑚), where 𝑥 = 0 at the tank’s
inlet, 𝑈 is the velocity along the 𝑥 direction (𝑚/𝑚𝑖𝑛), and 𝑘 is a reaction rate (𝑚𝑖𝑛−1 ) whereby
the chemical decays to another form. Develop an explicit scheme to solve this equation
numerically. Test it for 𝑘 = 0.15, 𝐷 = 100, and 𝑈 = 1 for a tank of length 10 𝑚. Use ∆𝑥 =
1𝑚, and a step size ∆𝑡 = 0.005. Assume that the in flow concentration is 100 and that the
initial concentration in the tank is zero. Perform the simulation from 𝑡 = 0 to 100 and plot the
final resulting concentrations versus 𝑥. You may use any software of your choice to solve this
problem. In your submission, clearly write the equations and submit the code along with it.
4. Solve the problem given in L19, slide 12 using the Crank-Nicolson method and ∆𝑡 = 2𝑠.
Clearly show all the equations and steps in your calculation. You may write a code to solve for
the final values. Please submit your code along with the submission.
Extra questions for students taking this as replacement for MA214. Regular students may treat
these as additional practice questions:
5. Determine the temperature distribution and fluxes for the plate shown below. The plate is
60x60x1 cm and is made out of aluminum (𝑘 ′ = 0.49 𝑐𝑎𝑙/(𝑠. 𝑐𝑚. ℃)), with an input of
10 𝑐𝑎𝑙/𝑠 into the middle node (2,2).