Finite difference solver for two-phase heat transfer in porous media with Brinkman flow

This dataset contains the source code of a standalone Python application for simulating transient heat transfer in a porous medium composed of a moving fluid (liquid) and a stationary solid matrix. The solver accounts for local thermal non-equilibrium between phases (separate energy equations for solid and fluid), anisotropic effective thermal conductivity, internal heat sources, and convective heat exchange between phases.

Four velocity profiles are implemented: constant, parabolic, cosine, and Brinkman (Darcy–Brinkman flow with prescribed flow rate). The code solves the Brinkman equation via a tridiagonal matrix algorithm to obtain the velocity profile used in the fluid energy equation.

The graphical user interface (Tkinter) provides interactive contour plots of temperature fields (solid/liquid), velocity profiles, and time–temperature diagrams. The solver includes stability checks (NaN/inf detection, bounds). Results can be saved as compressed NumPy .npz files.