fd1d_heat_steady_test

fd1d_heat_steady_test, a Fortran90 code which calls fd1d_heat_steady(), which applies the finite difference method to estimate the solution of the steady state (time independent) heat equation over a one dimensional region, which can be thought of as a thin metal rod.

Licensing:

The information on this web page is distributed under the MIT license.

Related Data and codes:

fd1d_heat_steady, a Fortran90 code which uses the finite difference method (FDM) to solve the steady (time independent) heat equation in 1D.

Source Code:

• fd1d_heat_steady_test.sh, runs all the tests.

• problem1.f90, uses K(X) = 1, F(X) = 0, so the solution should be the straight line that connects the boundary values.
• problem1_nodes.txt, the coordinates of the nodes.
• problem1_values.txt, the computed temperatures at the nodes.
• problem1.png, a PNG image of the solution.

• problem2.f90, uses K(X) which is set to different constants over three subregions, and F(X) = 0.0, so the solution will be a piecewise linear function that connects the boundary values.
• problem2_nodes.txt, the coordinates of the nodes.
• problem2_values.txt, the computed temperatures at the nodes.
• problem2.png, a PNG image of the solution.

• problem3.f90, uses K(X) = 1, F(X) defines a heat source, so the solution can rise above the boundary values.
• problem3_nodes.txt, the coordinates of the nodes.
• problem3_values.txt, the computed temperatures at the nodes.
• problem3.png, a PNG image of the solution.

• problem4.f90, uses K(X) = 1, F(X) defines a heat source and a heat sink, so the solution can go above and below the boundary values.
• problem4_nodes.txt, the coordinates of the nodes.
• problem4_values.txt, the computed temperatures at the nodes.
• problem4.png, a PNG image of the solution.