Scribd
Upload
19 views2 pages

MM210A6

Uploaded by

navneetacharya07
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
19 views2 pages

MM210A6

Uploaded by

navneetacharya07
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

MM210: Assignment 6

Due: 11:59 pm of April 15, 2024


Please show all steps in your calculations. No marks will be given if only the final answers are
given, with incomplete steps.
While you are encouraged to solve these problems using computer codes, the MM210
evaluation will be only based on written mathematical expressions and their calculated values.
You can use a calculator/computer to calculate the numerical values for the mathematical
expressions.
Please upload your scanned/typed solutions on the moodle submission portal.

1. The following system of equations is designed to determine concentrations (the 𝑐’s in


𝑔/𝑐𝑚3 ) in a series of coupled reactors as a function of the amount of mass input to each reactor
(the right hand sides in 𝑔/𝑑𝑎𝑦):
15𝑐1 − 3𝑐2 − 𝑐3 = 3800
−3𝑐1 + 18𝑐2 − 6𝑐3 = 1200
−4𝑐1 − 𝑐2 + 12𝑐3 = 2350
Solve for the concentrations in the three reactors, using:
(a) Gauss elimination,
(b) Gauss-Seidel method until the percent relative error falls below 𝜀𝑠 = 5%.

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).

You might also like