This function solves the three-dimensional Pennes Bioheat Transfer (BHT) equation in a homogeneous medium using Alternating Direction Implicit (ADI) method. The code has been developed for High-Intensity Focused Ultrasound (HIFU) treatments in tissue, but it can be applied to other heating problems as well. If desired, the solution takes into account the perfusion rate, thermal conductivity and specific heat capacity of tissue.
ADI_method.pdf - Written down numerical solution to heat equation using ADI method
solve_heat_equation_implicit_ADI.m - Code for the numerical solution using ADI method
thomas_algorithm.m - Fast algorithm for solving tridiagonal matrices
comparison_to_analytical_solution.m - An example code for comparing the solutions from ADI method to an analytical solution with different heating and cooling durations
Visa Suomi (2023). 3-D Heat Equation Numerical Solution (https://www.mathworks.com/matlabcentral/fileexchange/59336-3-d-heat-equation-numerical-solution), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform CompatibilityWindows macOS Linux
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!
- Checked compatibility with Matlab R2018a and above
- The outputs are now given as 4-D matrices for temperature and thermal dose over time