fd1d_heat_steady_test

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

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Related Data and Programs:

fd1d_heat_steady, an Octave code which uses the finite difference method to solve the steady heat equation in 1D.

Source Code:

• fd1d_heat_steady_test.m, calls all the tests;
• fd1d_heat_steady_test.sh, runs all the tests;
• fd1d_heat_steady_test.txt, the output file;

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

• test2.m, 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.
• test2_nodes.txt, the coordinates of the nodes.
• test2_values.txt, the computed temperatures at the nodes.
• test2.png, a PNG image of the solution.

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

• test4.m, 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.
• test4_nodes.txt, the coordinates of the nodes.
• test4_values.txt, the computed temperatures at the nodes.
• test4.png, a PNG image of the solution.