Scribd
Upload
71 views

CFD - Assignment 1: Report

The document provides instructions to solve the 2D steady-state heat equation on a rectangular plate using finite difference methods and compares the results to an analytical solution. It requires: 1) Solving the heat equation using Gauss-Seidel and PSOR methods on a grid with specified dimensions and convergence criteria. 2) Computing the analytical solution by summing the Fourier series up to n=110 terms. 3) Presenting the temperature values at grid points from each method in a table for comparison. 4) Plotting temperature contours from the Gauss-Seidel results.

Uploaded by

Atharva Shinde
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
71 views

CFD - Assignment 1: Report

The document provides instructions to solve the 2D steady-state heat equation on a rectangular plate using finite difference methods and compares the results to an analytical solution. It requires: 1) Solving the heat equation using Gauss-Seidel and PSOR methods on a grid with specified dimensions and convergence criteria. 2) Computing the analytical solution by summing the Fourier series up to n=110 terms. 3) Presenting the temperature values at grid points from each method in a table for comparison. 4) Plotting temperature contours from the Gauss-Seidel results.

Uploaded by

Atharva Shinde
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 2

CFD - Assignment 1

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.
It is required to obtain the steady-state temperature distribution on a 2D rectangular plate
as shown below.

Dirichlet B.C.: y = 0, T = T1; x = 0, T = T2; y = H, T = T3; x = L, T = T4;


(T1 = 100oC, T2 = T3 = T4 = 0)
Governing equation for steady 2D heat conduction (which if there is no heat generation is
a Laplace equation):
 2T  2T
 2T  2  2  0
x y
It is intended to solve this elliptic PDE subject to imposed boundary conditions. Select
constant spatial step sizes ∆x = ∆y = 0.05 m, in which case imax = 21 and jmax = 41. The
temperature distribution is to be computed for a total of (imax -2)x(jmax -2) = 741 grid
points. Condition of convergence to be used is if ERROR < ERROR max, the solution has
converged where
j  jmax 1
i imax 1

ERROR = 
i 2
Ti ,nj1  Ti ,nj
j 2

ERRORmax = 0.01
(a) 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:
(b)
  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

i j T (Gauss-Seidel) T (PSOR) T (Analytical)

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.

You might also like