Computational Assignment AE603
Computational Assignment AE603
1. Finite Difference (FD) Code for Solving the Heat Equation using direct solution
method (Gaussian elimination followed by back substitution):
Write a computer code that implements the Finite Difference (FD) method for solving the
two-dimensional steady state heat conduction equation with heat source on a square
plate. Essentially this results in a system of linear equations, Ax=b form. For simplicity
use the following assumptions:
a) The mesh is uniform.
b) A uniformly distributed source.
c) The material composition is homogeneous, i.e. constant thermal conductivity: k.
d) Temperature of the walls prescribed, on all four faces of the square plate.
Write your code in fortran/C/matlab only. Include a README text file that states:
i. How to compile the code: the line command necessary to compile.
ii. How to execute the code: the line command necessary to execute.
iii. Status of your code: operational/compiles-doesnt-run/doesnt compile.
iv. Briefly describe the problem solved by your code, expected input, resulting
output, any limitations or restrictions.
A sample input file (using the values listed below) and corresponding output file
produced by your code. Format of the input file may be in the following
sequence (all values in SI units):
Line
Number
Value
Meaning
8/12/16/20
300.0
Tw = Wall temperature
100.0
50.0
k = Thermal conductivity
Note: your code should function for general (but permissible) input values not just the above;
the sample input serves only as a check of your successful completion of the assignment.
Line
Number
Value
Meaning
8/12/16/20
300.0
Tw = Wall temperature
100.0
50.0
k = Thermal conductivity
6.
1.e-6
7.
0/1/2
2
, where I is as
1 sin / I
defined above.
c. Compare the direct solution method with the GS method with respect to the storage
requirement (as a function of I).
d. Compare the three iterative methods with respect to the number of iterations needed for
the solution to converge to the pre-specified tolerance, i.e., Max( ) c , where is the
relative convergence error.