ME502: Applied Numerical Methods I.I.T.

- Ropar
1st Semester 2022-23
Practice Problem Set #6
General Instructions:
1. Study and revise class notes and the book(s)
2. Write the references for any information (such as figure, table, data, website) that is being cited in your
assignments. Use some standard format to write the references. (For instance ASME format - which is given here
https://www.asme.org/publications-submissions/proceedings/author-guidelines/elements-of-a-paper/references)

Practice Problems/Examples: Students are strongly encouraged to solve these problems on their
own. It is not mandatory for the students to submit the solutions for these ‘practice problems’.

1. (a) Use Euler’s method to solve the following equation.

𝑑𝑦
= 2.7𝑥 2 − 1.8𝑥 + 13.1
𝑑𝑥

The interval over which the solution is required is from x = 0 to x = 4. It is known to us that at
x = 0, y = 2. This problem may be solved using multi step sizes, h = 0.5, h = 0.25, h = 0.1.

(b) Repeat the above problem using Heun’s method.

2.
For a given wall (shown in figure), it is known that 1D steady state
heat transfer is taking place across it (in the x-direction). There is some
heat generation taking place inside the wall, at the rate of egen (W/m3).
The thickness of the wall is L, and its thermal properties are as follows:
thermal conductivity k, density ρ, specific heat cp.

Using the FDM (Finite Different Method) technique for n = 5 nodes,

Condition-A: The left surface of the wall (i.e. at x = 0) is maintained at

Condition-B: The left surface of the wall (i.e. at x = 0) has a constant

[Hint: Do not solve the equations; just write the various equations that
are required to be solved. Start by drawing the nodes boundaries.]

1
3. For a given wall (shown in figure below), it is known that 1D transient state heat transfer is
taking place across it (in the x-direction). There is some heat generation taking place inside the
wall, at the rate of egen (W/m3). The thickness of the wall is L, and its thermal properties are as
follows: thermal conductivity k, density ρ, specific heat cp.
Using the FDM (Finite Different Method) implicit technique for n = 5 nodes, write the
discretized form of equations which will give the solution for the temperature distribution
within the wall at some time t, under the following conditions.

Condition-A: Initially the whole wall is at some uniform temperature, Ti. At time t = 0 sec, the
left surface of the wall (i.e. at x = 0) is maintained at some fixed temperature, TL, and the right
surface of the wall (i.e. at x = L) is maintained at some fixed temperature, TR.

Condition-B: Initially the whole wall is at some uniform temperature, Ti. At time t = 0 sec, the
left surface of the wall (i.e. at x = 0) is heated convectively by the air surrounding it which is at
a temperature of TL, and with a convective heat transfer coefficient of hL. Similarly at time t = 0
sec, the right surface of the wall (i.e. at x = L) is heated convectively by the air surrounding it
which is at a temperature of TR, and with a convective heat transfer coefficient of hR.

[Hint: Do not solve the equations; just write the various equations that are required to be
solved.]

2
4. For a given object (shown in figure below), it is known that 2D steady state heat transfer is
taking place across it (in the x-direction & y-direction). There is some heat generation taking
place inside this object, at the rate of egen (W/m3). The height of the object is L, and its thickness
is 1.5L. Its thermal properties are as follows: thermal conductivity k, density ρ, specific heat cp.

Using the FDM (Finite Different Method) technique, write the discretized form of equations
which will give the solution for the temperature distribution within this object under the
following conditions.

Condition-A: The north side of the object is convectively heated by air (having temperature TN,
and convective heat transfer coefficient hN). The east side of the object is convectively heated by
air (having temperature TE, and convective heat transfer coefficient hE). The south side of the
object is convectively heated by air (having temperature TS, and convective heat transfer
coefficient hS). The west side of the object is convectively heated by air (having temperature TW,
and convective heat transfer coefficient hW).

Condition-B: The north side of the object is maintained at a fixed temperature of TN. The east
side of the object is convectively heated by air (having temperature TE, and convective heat
transfer coefficient hE). The south side of the object is insulated. The west side of the object is
convectively heated by air (having temperature TW, and convective heat transfer coefficient hW).

[Hint: Do not solve the equations; just write the various equations that are required to be
solved.]

3
5. For a given object (shown in figure below), it is known that 2D transient state heat transfer is
taking place across it (in the x-direction & y-direction). There is some heat generation taking
place inside this object, at the rate of egen (W/m3). The height of the object is L, and its thickness
is 1.5L. Its thermal properties are as follows: thermal conductivity k, density ρ, specific heat cp.

Using the FDM (Finite Different Method) implicit technique, write the discretized form of
equations which will give the solution for the temperature distribution within this object at some
time t, under the following conditions.

Condition-A: Initially the whole wall is at some uniform temperature, Ti. At time t = 0 sec, the
object is exposed to the following conditions on its 4 sides. The north side of the object is
convectively heated by air (having temperature TN, and convective heat transfer coefficient hN).
The east side of the object is convectively heated by air (having temperature TE, and convective
heat transfer coefficient hE). The south side of the object is convectively heated by air (having
temperature TS, and convective heat transfer coefficient hS). The west side of the object is
convectively heated by air (having temperature TW, and convective heat transfer coefficient hW).

Condition-B: Initially the whole wall is at some uniform temperature, Ti. At time t = 0 sec, the
object is exposed to the following conditions on its 4 sides. The north side of the object is
maintained at a fixed temperature of TN. The east side of the object is convectively heated by air
(having temperature TE, and convective heat transfer coefficient hE). The south side of the object
is insulated. The west side of the object is convectively heated by air (having temperature TW,
and convective heat transfer coefficient hW).

[Hint: Do not solve the equations; just write the various equations that are required to be
solved.]

4
6. Solve the following differential equation using numerical method.

𝑑𝑦
= 1.2𝑥 3 + 2.3𝑥 2 − 1.3𝑥 + 17.5
𝑑𝑥

The interval over which the solution is required is from x = 0 to x = 5. It is known to us that at
x = 0, y = 3. This problem may be solved using multi step sizes, h = 0.5, h = 0.25, h = 0.1.

Use RK4 method to solve this problem.