CFD - Assignment
CFD - Assignment
CFD - Assignment
Submit by date..
Note:
(i) You are required to enclose a printout of the code(s) you have developed.
(ii) Please write down the date of submission conspicuously on the cover page of your
report.
j jmax 1
i imax 1
ERROR =
i2
Ti ,nj1 Ti ,nj
j 2
ERRORmax = 0.01
Use an initial guess T = 0.0 at all interior points.
(a) Use point Gauss-Seidel method to compute the temperature distribution and write
down the number of iterations required. Present the results in the format given
by TABLE 1.
(b) Do the same using PSOR with a relaxation parameter given by the optimal
relation:
2
cos
2
cos
opt
2 2 1
, where = imax 1 jmax 1 , x
1 2 y
(c) Perform numerical experiments with taking values from 0.8 to 2.0 at intervals
of 0.1 and find the number of iterations required for convergence with each value
of . Give an (vertical axis) vs Number of iterations for convergence
(horizontal axis) plot.
(d) Find the temperature distribution from the analytical solution given below at the
grid nodes and present the result in the format given in TABLE 1:
(e)
n ( H y )
1 (1) n sinh
T T1 2 L sin n x
n 1 n n H L
sinh
L
Carry out the summation for n = 1, , 110 and use `long double
TABLE 1
11 1 - - -
11 2 - - -
- - - - - - - - - - - - - - - - - - - - - - - --
- - - - - - - - - - - - - - - - - - - - - - - --
11 41 - - -
Draw contours of constant temperature for the rectangular plate from results obtained by
Gauss-Seidel method. You may use TECPLOT (preferred) or GNU plot.