Double exponential model of corpse cooling under conditions of linearly varying ambient temperature

Cover Page

Cite item

Full Text

Abstract

BACKGROUND: The main condition for the correctness of determining the postmortem interval by the method of thermometry of the deep tissues of the corpse is the constancy of the ambient temperature. This condition significantly limits the range of application of the method. The priority of thermometry of the core of the body in the diagnosis of prescription of death is explained by the slower cooling of deep tissues, which allows to increase the duration of the postmortem interval available for diagnosis, and less exposure to the influence of various random factors on the cooling process. The finite element models proposed recently can take into account almost all essential cooling conditions, including changes in ambient temperature, however, due to their high complexity, they require serious physical and mathematical training and technical skills, expensive software and postmortem computed tomography. For these reasons, they have not yet found wide application in expert practice.

In this article, a mathematical model of cooling the core of a corpse at a linearly varying ambient temperature is proposed.

AIMS: Construction a mathematical model of cooling the core of a corpse based on the Marshall-Hoare phenomenological law under conditions of linearly varying external temperature, to find a numerical algorithm for solving the model and to develop a computer program that implements it.

MATERIAL AND METHODS: A direct analytical modeling of the corpse cooling under conditions of linearly varying ambient temperature was carried out, performed on the basis of the Marshall-Hoare phenomenological cooling law and focused on solving the problem of determination of the postmortem interval by rectal or cranioencephalic temperature.

RESULTS: A mathematical model of cooling the core of a corpse under conditions of linearly varying ambient temperature has been developed. The chord method is proposed as a numerical algorithm for solving this model. The developed mathematical model and an iterative algorithm for its solution, as well as procedures for calculating interval estimates of the postmortem interval, are implemented in the C# language in the format of the Warm Bodies MHNH computer program.

CONCLUSIONS: It is advisable to use the proposed model and the program implementing it in forensic medical expert practice when determining the postmortem interval by the rectal or cranioencephalic temperature of a corpse in conditions of linearly varying ambient temperature.

About the authors

German V. Nedugov

Samara State Medical University

Author for correspondence.
Email: nedugovh@mail.ru
ORCID iD: 0000-0002-7380-3766
SPIN-code: 3828-8091
Scopus Author ID: 25947646500
ResearcherId: ABH-5590-2020

MD, Dr. Sci. (Med.), Associate Professor

Russian Federation, 89 Chapaevskaya street, Samara, 443099

References

  1. Marshall TK, Hoare FE. Estimating the time of death. The rectal cooling after death and its mathematical expression. J Forensic Sci. 1962;7(1):56–81.
  2. Henssge C. Death time estimation in case work. I. The rectal temperature time of death nomogram. Forensic Sci Int. 1988;38(3-4):209–236. doi: 10.1016/0379-0738(88)90168-5
  3. Henssge C. Rectal temperature time of death nomogram: dependence of corrective factors on the body weight under stronger thermic insulation conditions. Forensic Sci Int. 1992;54(1):51–66. doi: 10.1016/0379-0738(92)90080-g
  4. Althaus L, Henssge C. Rectal temperature time of death nomogram: sudden change of ambient temperature. Forensic Sci Int. 1999;99(3):171–178. doi: 10.1016/s0379-0738(98)00188-1
  5. Bisegna P, Henssge C, Althaus L, Giusti G. Estimation of the time since death: sudden increase of ambient temperature. Forensic Sci Int. 2008;176(2-3):196–199. doi: 10.1016/j.forsciint.2007.09.007
  6. Nedugov GV. Mathematical modeling of the corpse cooling under conditions of varying ambient temperature. Russian Journal of Forensic Medicine. 2021;7(1):29–35. (In Russ.) doi: 10.17816/fm360
  7. Nedugov GV. Numerical method for solving double exponential models of corpse cooling in the determination of the time of death. Forensic Medical Expertise. 2021;64(6):25–28. (In Russ.) doi: 10.17116/sudmed20216406125
  8. Madea B. Methods for determining time of death. Forensic Sci Med Pathol. 2016;12(4):451–485. doi: 10.1007/s12024-016-9776-y
  9. Henssge C, Madea B. Estimation of the time since death in the early post-mortem period. Forensic Sci Int. 2004;144(2-3):167–175. doi: 10.1016/j.forsciint.2004.04.051
  10. Nedugov GV. New computer technologies for determining the postmortem interval by the Henssge method. Russian Journal of Forensic Medicine. 2021;7(3):152–158. (In Russ.) doi: 10.17816/fm406
  11. Mall G, Eisenmenger W. Estimation of time since death by heat-flow Finite-Element model. Part I: method, model, calibration and validation. Leg Med (Tokyo). 2005;7(1):1–14. doi: 10.1016/j.legalmed.2004.06.006
  12. Mall G, Eisenmenger W. Estimation of time since death by heat-flow Finite-Element model. Part II: application to non-standard cooling conditions and preliminary results in practical casework. Leg Med (Tokyo). 2005;7(2):69–80. doi: 10.1016/j.legalmed.2004.06.007
  13. Schenkl S, Muggenthaler H, Hubig M, et al. Automatic CT-based finite element model generation for temperature-based death time estimation: feasibility study and sensitivity analysis. Int J Legal Med. 2017;131(3):699–712. doi: 10.1007/s00414-016-1523-0

Supplementary files

Supplementary Files
Action
1. JATS XML
2. Fig. 1. Temperature curves from examples 1 and 2.

Download (242KB)
3. Fig. 2. The geometry of the solution of equation (9) by the chord method for the data of Example 1. The first 3 iterations are shown at t0=10 and t1=0,1 h.

Download (210KB)
4. Fig. 3. The geometry of the solution of equation (10) by the chord method for the data of Example 2. The first 5 iterations are shown at t0=0,1 and t1=20 h.

Download (191KB)
5. Fig. 4. The Excel workspace with the data of Example 1 for solving the function (9) in relation to rectal temperature by the chord method. The cells intended for the introduction of initial physical quantities are highlighted in color. The remaining numeric data is generated automatically.

Download (403KB)

Copyright (c) 2021 Nedugov G.V.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies