Scribd
Upload
8 views

Assignment#4

Uploaded by

Pooja Agarwal
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
8 views

Assignment#4

Uploaded by

Pooja Agarwal
Copyright
© © All Rights Reserved
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