An inverse analysis of the brain cooling process in neonates using the particle filter method

This study deals with the computational simulation and inverse analysis of the cooling treatment of the hypoxic-ischemic encephalopathy in neonates. A reduced-order model is implemented for real-time monitoring of the internal body temperatures. The purpose of this study is to sequentially estimate the transient temperatures of the brain and other body regions with reduced uncertainties.

Pennes’ model was applied in each body element, and Fiala’s blood pool concept was used for the solution of the forward bioheat transfer problem. A state estimation problem was solved with the Sampling Importance Resampling (SIR) algorithm of the particle filter method.

The particle filter method was stable and accurate for the estimation of the internal body temperatures, even in situations involving large modeling and measurement uncertainties.

The proposed reduced-order model was verified with the results of a high-fidelity model available in the literature. Validation of the proposed model and of the solution of the state estimation problem shall be pursued in the future.

The solution of the state estimation problem with the reduced-order model presented in this paper has great potential to perform as an observer of the brain temperature of neonates, for the analysis and control of the systemic cooling treatment of neonatal hypoxic-ischemic encephalopathy.

The main treatment for hypoxic-ischemic encephalopathy in neonates is the cooling of affected regions. Accurate and fast models might allow the development of individualized protocols, as well as control strategies for the cooling treatment.

This paper presents the application of the SIR algorithm for the solution of a state problem during the systemic cooling of a neonate for the treatment of the hypoxic-ischemic encephalopathy.