💻 Fortran Source Code Library

Run calculations locally on your own machine. View code structure, read technical explanations, and download compilation packages including sample input files.

Composite Plane Wall Conduction

Core Numerical Engine in Fortran 90 • 0 total downloads

steady_heat_conduction_plane_walls.f90
! Error: Source file steady_heat_conduction_plane_walls.f90 not found on server.

Solver Description

Calculates steady-state 1D heat conduction through composite plane walls. It computes conduction resistances ($L_i / k_i A$) and convection resistance ($1 / h A$) to find the total thermal resistance ($R_{total}$). It solves for heat flux ($q'' = (T_0 - T_\infty)/R_{total}$) and uses Fourier's law to resolve all interface temperatures.

Key Numerical Methods & Architecture

  • Input Redirection: Reads parameters sequentially from standard input (`stdin`) using Fortran sequential read (`read(*,*)`), ensuring modular integration.
  • Modular Design: Formulated using pure mathematical routines, separation of equations from output formatting, and precise numerical solvers (e.g. bisection, Newton-Raphson).
  • Standard Compliant: Written in clean, standards-compliant Fortran 90 to ensure cross-compiler compatibility.

🛠️ Local Compilation

To test this code on your machine, compile the source code file(s) using a standard Fortran compiler (e.g., `gfortran`).

Compilation Command:

gfortran -ffree-line-length-none steady_heat_conduction_plane_walls.f90 -o heat_conduction

Execution Command:

Execute the program by feeding the sample input file into the program using stdin redirection:

heat_conduction < input.txt

📥 Downloads & Local Files

Preview of the required input file (input.txt):

! Number of layers (N)
3
! Base Temperature ($T_0$) [°C]
150.0
! Ambient Fluid Temperature ($T_\infty$) [°C]
20.0
! Layer 1 Thickness ($L_1$) [m]
0.05
! Layer 1 Conductivity ($k_1$) [W/m-K]
15.0
! Layer 2 Thickness ($L_2$) [m]
0.03
! Layer 2 Conductivity ($k_2$) [W/m-K]
60.0
! Layer 3 Thickness ($L_3$) [m]
0.02
! Layer 3 Conductivity ($k_3$) [W/m-K]
1.2
! Convection Coefficient ($h$) [W/m²-K]
10.0