Scribd
Upload
42 views

Lab Sheet For Computational Fluid Dynamics: College of Engineering, Pune

The document describes a lab assignment on computational fluid dynamics involving heat advection and convection simulations. It includes instructions for 1D and 2D simulations using different numerical schemes to model advection and convection, and to analyze the results. Students are asked to develop code to model heat transfer in fluids, run simulations, and produce plots and discussions of temperature contours and profiles for various test cases involving uniform flow, channel flow, and different grid and scheme parameters.

Uploaded by

Abul Lais Nalband
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
42 views

Lab Sheet For Computational Fluid Dynamics: College of Engineering, Pune

The document describes a lab assignment on computational fluid dynamics involving heat advection and convection simulations. It includes instructions for 1D and 2D simulations using different numerical schemes to model advection and convection, and to analyze the results. Students are asked to develop code to model heat transfer in fluids, run simulations, and produce plots and discussions of temperature contours and profiles for various test cases involving uniform flow, channel flow, and different grid and scheme parameters.

Uploaded by

Abul Lais Nalband
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 9

Page 1 of 9

Lab sheet for Computational Fluid Dynamics


College of Engineering, Pune
Lab # 5: Computational Heat Advection and Convection
Date:

Time:

1. 2D Computational Heat Advection (CHA):


Consider a 2D Cartesian (x,y) computational domain of size L=1m and H=1 m, for CHA of a fluid (=1000
kg/m3 and cp=4180 W/m.K) moving with a uniform velocity u=v=1 m/s and an initial temperature of 500C. The
left and top boundary of the domain is subjected to 1000C; and the bottom and right boundary to 00C
Develop the code for the above problem, run the code for three different advection schemes: (a)
FOU, (b) SOU and (c) QUICK. Take the maximum number of grid points in x-and y-direction as
imax=jmax=22 and convergence criteria as 0.000001. (Refer the code: Lab4_1_2DAdvection.sci).
Report the results as
a) Plot and discuss the steady state temperature contours for the different advection schemes (3
figures).
b) Plot and discuss the temperature profile at the vertical centerline(x=0.5), T(y), for the different
advection schemes (3 figures).

Page 2 of 9

Answer Sheet
2D Computational Heat Advection (CHA):
a) Plot and discuss the steady state temperature contours for the different advection schemes (3 figures).
b) Plot and discuss the temperature profile at the vertical centerline(x=0.5), T(y), for the different advection
schemes (3 figures).

Discuss Fig. 4.1 here, limited inside this text box only
(a1)

(a2)

(b1)

(b2)

(c1)

(c2)

Fig. 4.1: Steady state temperature contours using the (a1) FOU scheme (b1) SOU scheme (c1) QUICK
scheme Temperature variation along the vertical centerline using the (a2) FOU scheme (b2)
SOU scheme (c2) QUICK scheme

Page 3 of 9

Lab sheet for Computational Fluid Dynamics

Page 4 of 9

College of Engineering, Pune


Lab # 6: Computational Heat Advection and Convection
Date:

Time:

1D Computational Heat Convection (CHC):


Consider one-dimensional unsteady-state convection, in a domain of length L=1. The fluid has a density =1,
diffusion coefficient =0.02 and a velocity of U=1; resulting in a Peclet number Pe= UL/ = 50. The
boundary condition for the advected variable is =0 at x=0 and =1 at x=1. The initial condition is =0.5.
Develop the code for the above problem, run the code for two different advection schemes: (a) FOU and (b) CD.
Take the convergence criteria as 0.000001.
Plot and discuss the steady state temperature profile T(x) for the different advection schemes on two different
grid sizes:
a) imax=12 (2 figures)
b) imax=42 (2 figures)
Plot and discuss the steady state temperature profile T(x) for the FOU and CD advection schemes on two
different grid sizes: (a) imax=12 (2 figures) and (b) imax=42 (2 figures). (Refer the code:
Lab4_2_1DConvection.sci)

Answer Sheet
1D Computational Heat Convection (CHC):

Page 5 of 9
Plot and discuss the steady state temperature profile T(x) for the FOU and CD advection schemes on two
different grid sizes: (a) imax=12 (2 figures) and (b) imax=42 (2 figures).

Discuss Fig. 4.2 here, limited inside this text box only

(a1)

(a2)

(b1)

(b2)

Fig. 4.2: Steady state temperature profile T(x) using the (a1,a2) FOU scheme (b1,b2) CD scheme, on a
grid size of (a1,b1) 12 and (a2,b2) 42 grid points.

Page 6 of 9

Lab sheet for Computational Fluid Dynamics

College of Engineering, Pune

Page 7 of 9
Lab # 7: Computational Heat Advection and Convection
Date:

Time:

2D Computational Heat Convection (CHC):


Consider a 2D Cartesian computational x-y domain of size L=6 unit and H=1 unit, for CHC with a prescribed
velocity field. This corresponds to a slug flow (u=1, v=0) of a fluid in a channel; subjected to a non-dimensional
temperature of 1 at the inlet and 0 at the walls. At the outlet, fully developed Neumann BC is used. The initial
condition for non-dimensional temperature of the fluid is 0.
Develop the code for the above problem, run the code for two different advection schemes: (a) FOU
and (b) QUICK; at Re=10 and Pr=1. Take the maximum number of grid points in x-and y-direction as imax=62
and jmax=22, respectively; and convergence criteria as 0.000001. (Refer the Code: Lab4_3_2DConvection.sci)
Report the results as
a) Plot and discuss the steady state temperature contours for the different advection schemes (2
figures).
b) Plot and discuss the temperature profile, T(y), at different axial locations (x/L=0.2, 0.4, 0.6, 0.8 and
1), for the different advection schemes (2 figures).

(a1)

Page 8 of 9

Discuss Fig. 4.3 here, limited inside this text box only

(a2)

(b1)

(b2)

Answer Sheet
1D Computational Heat Convection (CHC):
a) Plot and discuss the steady state temperature contours for the different advection schemes (2 figures).
b) Plot and discuss the temperature profile, T(y), at different axial locations (x/L=0.2, 0.4, 0.6, 0.8 and
1), for the different advection schemes (2 figures).

Fig. 4.3: Steady state temperature contours using the (a1) FOU scheme and (b1) QUICK scheme.
Temperature variation at different axial locations using the (a2) FOU scheme and (b2) QUICK
scheme.

Page 9 of 9

You might also like