Modeling and numerical simulation of bioheat transfer and biomechanics in soft tissue

A mathematical model describing the thermomechanical interactions in biological bodies at high temperature is proposed by treating the soft tissue in biological bodies as a thermoporoelastic media. The heat transfer and elastic deformation in soft tissues are examined based on the Pennes bioheat transfer equation and the modified Duhamel-Neuman equations. The three-dimensional governing equations based on the proposed model is discretized using a 19-point finite-difference scheme. The resulting large sparse linear system is solved by a preconditioned Krylov subspace method. Numerical simulations show that the proposed model is valid under our test conditions and the proposed numerical techniques are efficient.