ME 6151: Computational Heat and Fluid Flow Jan – May 2025
Assignment 5
Assigned date: 22 April 2025; Due date: 01 May 2025; Max marks: 100
1. Consider a one-dimensional steady incompressible flow in a porous medium described by the following
equations:
𝑑
(𝜌𝑢) = 0 (1)
𝑑𝑥
𝑑𝑝 1
− − 𝐶 (𝜌𝑢 2 ) = 0 (2)
𝑑𝑥 2
where 𝜌 = 1000 kg/m3 and 𝐶 = 0.1 m−1 everywhere in the domain. The domain has a length 𝐿 = 0.3
m and is divided into 3 equal control volumes as shown in figure 1. A co-located mesh is used to store
the velocities and pressures as indicated in the figure. The velocity at the inlet is given as 𝑢 1 = 1 m/s
and the pressure at the outlet is given as 𝑝 4 = 10 N/m2 .
(a) Consider the term corresponding to porosity, given by, 12 𝐶 (𝜌𝑢 2 ) as a source term, perform
source-term linearization and use this linearized source in all further calculations.
(b) Derive the discrete momentum equations for the velocities 𝑢 𝐴 , 𝑢 𝐵 and 𝑢𝐶 .
(c) Consider half-cell control volumes centered at nodes 1 and 4 and derive the discrete momentum
equations for face velocities 𝑢 1 and 𝑢 4 .
(d) Derive the pressure correction equations for the cells 𝐴, 𝐵 and 𝐶
(e) Taking initial guess everywhere for velocity as 0.5 m/s and for pressure as 10 N/m2 perform
one iteration of the SIMPLE algorithm and report the values of unknown velocity and pressure,
𝑢 𝐴 , 𝑢 𝐵 , 𝑢𝐶 , 𝑢 4 and 𝑝 1 , 𝑝 2 , 𝑝 3 , 𝑝 𝐴 , 𝑝 𝐵 , 𝑝 𝐶 ?
(f) Explain whether you need to perform further iterations to obtain a converged solution?
1 A 2 B 3 C 4
p4
u1
uA , pA uB , p B uC , pC
Figure 1: Computational domain and mesh for problem 1.
PLEASE TURN OVERLEAF
ME 6151: Computational Heat and Fluid Flow Jan – May 2025
2. Consider one-dimensional heat conduction as shown in figure 2. The domain has a length of 4 m, and is
discretized into 4 uniform cells. The thermal conductivity of the material is 𝑘 = 10 W/m.K, and there is
a source given by 𝑆 = (3000 − 20𝑇) W/m3 . The left-boundary of the domain is subjected to a constant
heat flux of 𝑞 ′′𝑎 = 1000 W/m2 and the right-boundary of the boundary is maintained at a constant
temperature of 𝑇𝑏 = 400 K. Two virtual coarse levels are constructed as shown in figure, to enable
the solution using algebraic multigrid method. Take an initial guess of 𝑇𝑖 = 300 K everywhere in the
domain and consider the initial error to be zero on the levels 1 and 2. Perform one V-cycle of algebraic
multigrid method using 𝑛1 = 0 and 𝑛2 = 1 as number of pre- and post-sweeps for Guass-Seidel. Obtain
and indicate the following during the solution:
(a) Report the values for all the coefficients for all cells on the fine mesh (level-0).
(b) Write the relevant equations that connect the level-1 coefficients with the level-0 coefficients.
(c) Report the numerical values for all the coefficients for cells 1 and 2 on level-1.
(d) Write the relevant equations that connect the level-2 coefficients with the level-0 coefficients.
(e) Report the numerical value for cell 1 on level-2.
(f) Perform one V-cycle of algebraic multigrid and indicate relevant values and the value of tem-
perature on the fine mesh at the end of one cycle.
qa′′ 1 2 3 4
Level 0 Tb
1 2
Level 1
1
Level 2
Figure 2: Computational domain and mesh for problem 2.
END OF THE DOCUMENT
