Scribd
Upload
8 views

Assignment#4

Uploaded by

Pooja Agarwal
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
8 views

Assignment#4

Uploaded by

Pooja Agarwal
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

Department of Metallurgical & Materials Engineering

Indian Institute of Technology Kharagpur

Assignment#4

Computer Programming and numerical methods in Metallurgical Engineering (MT60003)


Autumn, 2024-2025

1. Solve the following ODE from t=0 to 2 by (a) Euler with h =0.5 and 0.25, ii) With h=0.5
Heun’s method with h=0.5, iii) 4th order RK method with h=0.5 and 0.25. Euler: 0.1135
Heun: 1.884
RK4=2.513
dy
 yt 3  1.5 y
dt
y(0)=1;

Compare with exact solution.

2. Solve the following second order differential equation by 2nd order RK (Ralston): Convert the
2 second order
d y
t  y  0 ODE to two first
dt 2 order ODES.
Where, y0  4; y0  0 at x=4
y=45.615
Take h=0.1 and solve for x=0 to 4. Plot y versus x. Plot y  versus x. dy/dt=
-6.2211

3 A non-isothermal batch reactor can be described by the following equation: Maximum


temperature
=278oC
𝑑𝐶 5
(− ) Concentration
= −𝑒 𝑇 𝐶
𝑑𝑡 after 5 secs =
0.00359284
𝑑𝑇 5
(− )
g-mole/L
= 2000𝑒 𝑇 𝐶 − 2(𝑇 − 30)
𝑑𝑡
where, C is the concentration of the reactant and T is the temperature in K of the
reactor. Initially the reactor is at 30oC and has a concentration of reactant C of 0.5
g- mole/L. Calculate the evolution concentration and temperature of the reactor as
a function of time up to 5 secs (time step = 0.1) using 4th order RK method.

Use finer time step of 0.02s upto 0.5 secs and then coarser time step of 0.05s upto
5 seconds. Plot concentration and temperature Vs time, using separate scale for
temperature and concentration.

4 Isothermal reduction of iron ore-coal composite pellets can be represented Since we


start with 1
by the following coupled equations. mole of
hematite, we
will end up
with 2 moles
of Fe and at
each time
step the iron
balance
should be
maintained
at 2.

Where, H, M, W, F represent the concentrations for Hematite, Magnetite, Wustite


and Iron, respectively. Assume, the first order rate constant at 1300oC for H to
M, M to W, W to F are 1, 0.5, and 0.1/s respectively. Take H(0)=1.0, C(0)=3,
and M(0) = W(0) = F(0) = 0.

Calculate evolution of phases with time 2nd order RK (mid-point method) and
plot the evolution of H,M,W,F.C with time. Take time step =0.5 s and calculate
for 50 secs.

You might also like