An application of the finite volume method to the bio-heat-transfer-equation in premature infants

M. Ludwig and J. Koch and B. Fischer

ETNA
28
136--148
(2008)

In this report the development of a ﬁnite volume method for the time-accurate simulation of the temperature distribution in a premature infant inside an incubator or in an open radiant warmer is described. The real geometry of a premature infant is obtained from MRT-images. The infants thermoregulation is modeled by the so-called bio-heat-transfer-equation incorporating source terms and Neumann boundary conditions. The source terms describe the metabolic heat production, the blood ﬂow and the respiratorical water loss whereas the Neumann boundary conditions model the heat transfer by transepidermal water loss, radiation, convection and conduction. The numerical solution is carried out by the developed ﬁnite volume method whose spatial discretization is done by a 3D-mesh-generator from CFD. For the time integration a semi-implicit multistep method is used. The arising large, sparse linear systems are eﬃciently solved with a Krylov subspace method. Some successful test runs using real life data are presented.