基于有限元分析的死后头部对流热交换在线工具
- 作者: Nedugiv V.G.1, Nedugov G.V.2
-
隶属关系:
- Samara National Research University
- Samara State Medical University
- 期: 卷 10, 编号 4 (2024)
- 页面: 555-565
- 栏目: 技术报告
- ##submission.dateSubmitted##: 07.08.2024
- ##submission.dateAccepted##: 01.10.2024
- ##submission.datePublished##: 05.12.2024
- URL: https://for-medex.ru/jour/article/view/16175
- DOI: https://doi.org/10.17816/fm16175
- ID: 16175
如何引用文章
详细
背景。在确定死亡时间的热测诊断领域,死后热交换的有限元分析方法是一种前景广阔的技术。相比传统的尸体冷却现象学方程,该方法能够克服其局限性。然而,目前支持这一方法的软件通常价格昂贵,并要求用户具备独立设置求解条件的能力。本文提出了一种开放式在线工具,用于基于头部颅脑热测的尸体冷却有限元分析,具有简洁的界面,可用于确定死亡时间。
研究目的。开发一个用于分析死后头部对流热交换的有限元分析在线工具。
材料与方法。构建了一个可缩放的有限元头部模型,该模型以多层球体为近似,包含1311个节点和9277个有限元。开发了用于计算初始和死后头部温度场的算法,并以Python 3编程语言实现。
研究结果。开发了在线应用程序Simple Finite Element Model of Postmortem Convective Heat Transfer of the Head,其功能包括:考虑初始温度场的特性;头部主要解剖层的尺寸和热物理特性;对流热交换强度;诊断点坐标;尸体冷却过程中外部温度的变化。该工具能够生成诊断点和头部表面在死后24小时内的冷却曲线,并输出相关数值结果。此外,还可显示有限元模型的网格几何特性、温度沿从模型中心到诊断点积分路径的分布情况。
结论。该计算在线工具无需高性能计算机或用户具备专业知识,可在法医学实践中用于通过颅脑热测法确定尸体死亡时间。
全文:
Обоснование
Одним из наиболее перспективных и активно развиваемых современных направлений в области термометрической диагностики давности наступления смерти является метод конечно-элементного анализа посмертного теплообмена [1–3]. Данный подход позволяет преодолеть ограничения, присущие феноменологическим уравнениям охлаждения трупа, поскольку непосредственно основывается на физике теплопередачи и принципиально может быть применён к любым условиям охлаждения [4–7].
Метод конечных элементов, используемый при оценке посмертного теплообмена, является валидированным инструментом решения в сложных геометрических объектах и краевых условиях дифференциальных уравнений с частными производными второго порядка, таких как, например, уравнение теплопроводности, которое в пространственной области Ω размерности n с системой координат имеет вид
, (1)
где c ― удельная теплоёмкость; ρ ― плотность; λ ― теплопроводность; r ― координаты точки; ― искомая функция, задающая температуру в точке пространства Ω в момент времени τ на промежутке [0,∞); ― функция плотности источников тепла [6].
В процессе развития компьютерных технологий метод конечных элементов получил широкое распространение при решении различных термодинамических задач в технике, а также в медицине. В настоящее время на рынке присутствует значительное количество различных коммерческих пакетов программного обеспечения метода конечных элементов. Высокая стоимость этих продуктов, даже отечественного производства, существенно ограничивает возможности их использования в диагностике давности наступления смерти. Бесплатные же версии программных пакетов, выпускаемые в образовательных целях, характеризуются наличием ограничений на количество узлов в генерируемых сетках конечно-элементных моделей (КЭМ) или вовсе допускают постановку только двумерных задач. Кроме того, использование коммерческих решателей предполагает наличие у пользователей компетенций по самостоятельному построению геометрии КЭМ, заданию её теплофизических свойств и краевых условий. Это препятствует широкому внедрению метода конечных элементов для диагностики давности наступления смерти, сводя его к единичным эпизодам применения в научных целях, реализуемым к тому же коллективами исследователей, включающими в свой состав помимо судебных медиков математиков и программистов [2, 4–7].
В этой связи актуальность представляет разработка открытых судебно-медицинских онлайн-инструментов, адаптированных к решению конкретизированных задач конечно-элементного анализа теплообмена трупа и требующих от пользователя минимальных навыков в области метода конечных элементов. Соответствие названным условиям накладывает на подобный сервис облачных вычислений ряд следующих ограничений.
Основное требование к онлайн-решателю заключается в ограничении сложности геометрии КЭМ. Это связано с тем, что повышение сложности КЭМ сопровождается ростом вычислительных затрат и времени, необходимого для получения результата. Поэтому при достаточно сложной КЭМ с точной геометрией и конечно-элементной сеткой реализовать метод конечных элементов в онлайн-режиме не представляется возможным. Отсюда для онлайн-реализации КЭМ должна включать не весь труп, а только какую-то одну его часть, желательно с минимальной вариабельностью своего анатомического строения. Подобный подход помимо снижения вычислительных затрат позволяет также автоматизировать процесс создания геометрии КЭМ, избавляя пользователя от названного наиболее сложного этапа её разработки.
Второе требование заключается в лимитировании граничных условий теплообмена, всё многообразие которых желательно свести к какому-либо одному постоянно встречающемуся и преобладающему типу теплопередачи, каковым, в частности, является конвекция. Это также позволило бы автоматизировать процесс задания граничных условий, оставив пользователю только необходимость ввода величины коэффициента теплоотдачи.
С учётом названных ограничений одной из наиболее перспективных для конечно-элементного моделирования областей тела представляется голова, геометрия мозгового отдела которой очень близка к такому телу, как многослойный шар с равномерным распределением слоёв однородной структуры [8]. Поверхность шара имеет с касательной плоскостью только одну общую точку, что позволяет пренебречь потерями тепла путём теплопроводности при контакте головы с опорной поверхностью и не учитывать в КЭМ теплофизические свойства последней.
Цель работы ― разработка онлайн-инструмента конечно-элементного анализа посмертного конвективного теплообмена головы.
Материалы и методы
В работе представлена разработка КЭМ посмертного конвективного теплообмена головы с последующим построением на её основе онлайн-приложения, предназначенного для определения давности наступления смерти путём краниоэнцефальной термометрии трупа.
Геометрия конечно-элементной модели
В качестве геометрической модели головы взрослого человека и прежде всего её мозгового отдела использовали многослойный шар радиусом 98 мм, состоящий из пяти равномерно распределённых однородных слоёв: кожно-апоневротического лоскута (5 мм), костей свода черепа (5 мм), ликвора субарахноидального пространства (2 мм), а также поверхностного (17 мм) и глубокого (69 мм) слоёв головного мозга (86 мм). Два подслоя в блоке головного мозга выделяли с целью моделирования его начального температурного поля у человека, которое в поверхностном слое является неоднородным с монотонным понижением температуры [9, 10]. Непосредственно расчётной моделью являлось сечение шара диаметральной плоскостью с нулевой аппликатой всех его точек (рис. 1). Сгенерированная сетка КЭМ включала 1311 узлов и 9277 конечных элементов, в том числе 2552 треугольника и 6725 тетраэдров (см. рис. 1).
Рис. 1. Конечно-элементная модель мозгового отдела головы, заданная в онлайн-приложении.
Fig. 1. The finite element model of the brain region of the head, set in an online application.
Теплофизические свойства конечно-элементной модели
В онлайн-приложении предусматривали возможность произвольного задания теплофизических свойств выделенных в КЭМ анатомических слоёв головы. Поскольку на практике точное измерение удельной теплоёмкости, теплопроводности и плотности органов и тканей трупа затруднительно, по умолчанию задавали усреднённые теплофизические свойства актуальных тканей (табл. 1), указанные в специальной литературе [1, 6–10].
Таблица 1. Теплофизические свойства, по умолчанию заданные в конечно-элементной модели
Table 1. Thermophysical properties, set by default in the finite element model
Слой КЭМ | Плотность, кг/м3 | Удельная теплоёмкость, Дж/(К∙кг) | Коэффициент теплопроводности, Вт/(м∙К) | Начальная температура, °С* | Радиальная координата, мм* |
Мозг | 1050 | 3700 | 0,503 | 37,2 | 69 |
Мозг | 1050 | 3700 | 0,503 | 35,5 | 86 |
Ликвор | 1003 | 4174 | 0,618 | 35 | 88 |
Череп | 1500 | 2300 | 1,16 | 34 | 93 |
Скальп | 1000 | 2000 | 0,34 | 31 | 98 |
Примечание. * На внешней границе слоя. КЭМ ― конечно-элементная модель. Note. * On the outer boundary of the layer. КЭМ ― finite element model. |
Начальное температурное поле
Поиск начального температурного поля КЭМ сопровождался решением уравнения Лапласа () c краевыми условиями Дирихле, соответствовавшими постоянной температуре на поверхностях выделенных анатомических слоёв:
,
где Δt ― оператор Лапласа в декартовой системе координат трёхмерного пространства с осями ox, oy и oz, определяемый как
,
а символ ∂Ω обозначает границу области Ω.
По умолчанию в качестве условий Дирихле задавали усреднённые прижизненные температуры на внешних границах слоёв (см. табл. 1). Затем, используя слабую формулировку уравнения Лапласа, выводили матричное уравнение [11, 12], из которого находили начальную температуру в каждой точке КЭМ (рис. 2).
Рис. 2. Начальное температурное поле диаметрального сечения конечно-элементной модели.
Fig. 2. The initial temperature field of the diametrical section of the finite element model.
Посмертное температурное поле
На следующем этапе находили посмертное температурное поле расчётной области. Сумму выделяемой внутренними источниками тепловой энергии и теплопотерь за счёт радиационного теплообмена ввиду их небольшой величины и примерно одинаковых количеств принимали равной нулю [8], также предполагали, что в силу однородности выделенных в КЭМ анатомических слоёв их теплофизические свойства являются константами. С учётом названных допущений уравнение теплопроводности (1) принимает вид
,
где а ― коэффициент температуропроводности, м2/с, равный
,
где c ― удельная теплоёмкость, Дж/(К×кг); ρ ― плотность, кг/м3; λ ― коэффициент теплопроводности, Вт/(м×К).
В качестве начального условия поставленной задачи нестационарного теплопереноса принимали найденное на первом этапе стационарное начальное температурное поле
.
Для моделирования охлаждения в воздушной среде на внешнем ребре модели задавали конвективный теплообмен с воздухом, протекающий по закону Ньютона–Рихмана:
,
где n ― нормаль к поверхности охлаждаемой части тела, м; s ― индекс, обозначающий поверхность физического тела; α ― коэффициент теплоотдачи, Вт/(м2×К); ta ― температура внешней среды, К. Коэффициент теплоотдачи в расчётах задавали постоянным и одинаковым для всей поверхности КЭМ. Актуальный расчётный период принимали равным 24,5 ч с шагом интегрирования 600 с.
Валидация конечно-элементной модели
В ходе валидации онлайн-приложения проверяли адекватность выполняемых им расчётов метода конечных элементов и адекватность геометрии КЭМ и заданных теплофизических свойств.
Адекватность метода конечных элементов тестировали путём оценки соответствия посмертных температурных полей КЭМ при различных условиях конвективной теплоотдачи и внешнего температурного режима с таковыми, полученными с помощью бесплатной версии приложения ELCUT 6.6 (https://elcut.ru/downloads_r.htm) для КЭМ с аналогичной геометрией, теплофизическими свойствами и краевыми условиями.
Адекватность геометрии КЭМ и заданных теплофизических свойств подтверждали путём оценки соответствия посмертных температурных кривых для центра и поверхности КЭМ, полученных с помощью разработанного онлайн-приложения, таковым, вычисленным по уравнениям с общепризнанной валидностью, основанным на законах охлаждения Ньютона–Рихмана и Marshall–Hoare [8, 13].
В качестве валидного описания посмертной динамики краниоэнцефальной температуры использовали уравнение Henssge:
, (2)
где t ― краниоэнцефальная температура, °С; tа ― внешняя температура, °С; τ ― давность наступления смерти, ч [8, 14].
В качестве эталона динамики охлаждения поверхности головы применяли уравнение, основанное на законе охлаждения Ньютона-Рихмана:
, (3)
где t ― поверхностная температура, °С; t0 ― начальная поверхностная температура, °С; k ― константа охлаждения, ч−1 [8, 14].
В уравнении (3) начальную температуру поверхности головы задавали равной 31°С, а константу k определяли по формуле:
,
где t1 и t2 ― температура поверхности головы в 23 ч и 24 ч посмертного периода соответственно, °С; Δτ ― промежуток времени между указанными измерениями, равный 1 ч [8].
В ходе валидации КЭМ находили и оценивали её посмертное температурное поле при внешней температуре от 10°С до 22°С и коэффициенте теплоотдачи от 5 до 25 Вт/(м2×К).
Язык программирования
Код онлайн-приложения, предназначенного для операционных систем Windows и Linux, составляли на языке программирования Python 3. При разработке пользовательского интерфейса также использовали библиотеку React из языка JavaScript. В качестве асинхронной очереди задач на расчёт и брокера сообщений применяли инструменты Celery и Redis.
Результаты
На основе разработанной КЭМ было создано приложение Simple Finite Element Model of Postmortem Convective Heat Transfer of the Head (свидетельство о государственной регистрации программы для ЭВМ № 2024668362) объёмом 1,779 Гб, предназначенное для конечно-элементного онлайн-анализа посмертного конвективного теплообмена головы (https://forenscalc.ru/simple-fem).
Для работы с онлайн-приложением пользователю необходимо указать координаты диагностической точки из области мозгового отдела головы, температуру внешней среды в посмертном периоде, коэффициент теплоотдачи, толщину и теплофизические свойства анатомических слоёв головы (начальную температуру, плотность, теплопроводность и удельную теплоёмкость). Результатом работы приложения является визуализация внешнего температурного тренда и кривых охлаждения в диагностической точке и на поверхности головы в первые сутки посмертного периода с генерацией числовых значений указанных температур в виде файла с расширением .csv, доступного в личном кабинете пользователя. Построение графиков выполняется приложением в режиме реального времени (real-time). Помимо числовых данных, в личном кабинете пользователя дополнительно формируются файлы с расширениями .msh и .png, содержащие графическую и табличную информацию о геометрии и свойствах сетки КЭМ, а также начальную термограмму вдоль контура интегрирования от центра КЭМ до диагностической точки.
Сравнение кривых охлаждения, полученных в одинаковых условиях с помощью разработанного онлайн-приложения и программы ELCUT 6.6, установило их практически полное соответствие. Незначительные несовпадения были обусловлены небольшим количеством узлов и элементов в сетке КЭМ бесплатной версии приложения ELCUT 6.6. и, следовательно, меньшей точностью её вычислений. Дальнейшая отладка КЭМ в различных условиях конвективного теплообмена с отсутствием внешней радиации и контакта с другими физическими телами показала максимальную сходимость её вычислительных результатов с таковыми уравнений (2) и (3) при задании коэффициента теплоотдачи, равного 6 Вт/(м2×К).
Помимо ввода константных значений внешней температуры, приложение допускает также задание любых типов её изменений, описываемых функциями из библиотеки Python 3. Для этого в соответствующее текстовое поле достаточно ввести составленную в соответствии с синтаксисом Python-библиотеки NumPy функцию внешней температуры от давности наступления смерти, измеряемой в секундах.
При отсутствии у судебно-медицинского эксперта информации о биометрических и теплофизических свойствах тканей головы конкретного трупа или их начального температурного поля, а также характеристик конвективного теплообмена следует проводить расчёты, используя заданные по умолчанию значения, в том числе коэффициента теплоотдачи, равного 6 Вт/(м2×К). В этом случае при постоянной внешней температуре результаты метода конечных элементов для точки начала координат в КЭМ будут соответствовать таковым уравнения C. Henssge.
Приведём пример. Внешняя температура в посмертном периоде изменялась по синусоидальной кривой, описываемой уравнением
,
где t ― давность наступления смерти, с (рис. 3).
Рис. 3. Графики охлаждения центра и поверхности конечно-элементной модели при синусоидальных изменениях внешней температуры. Условные обозначения: 1 ― температура в точке с нулевой радиальной координатой; 2 ― температура в точках поверхности конечно-элементной модели; 3 ― внешняя температура. ДНС ― давность наступления смерти.
Fig. 3. Graphs of cooling of the center and surface of the finite element model with sinusoidal changes in external temperature. Symbols: 1 ― temperature at a point with zero radial coordinate; 2 ― temperature at points on the surface of the finite element model; 3 ― ambient temperature. ДНС ― postmortem interval.
Необходимо вычислить динамику краниоэнцефальной температуры трупа при усреднённых биометрических и теплофизических параметрах тканей головы и конвективном теплообмене с коэффициентом теплоотдачи 6 Вт/(м2×К).
Для решения задачи в соответствующих текстовых полях онлайн-приложения Simple Finite Element Model of Postmortem Convective Heat Transfer of the Head необходимо указать декартовы координаты диагностической точки, равные нулю. Графические результаты расчётов представлены на рис. 3.
Технически приложение допускает ввод только абсциссы и ординаты точки, аппликата которой всегда равна нулю. По умолчанию приложение всегда производит расчёты для точки центра КЭМ с нулевыми координатами, а также для точек поверхности головы с радиальной координатой, равной 98 мм.
Разработанное приложение сохраняет в личном кабинете пользователя историю всех расчётов. По этой причине, а также для обеспечения конфиденциальности данных для входа на веб-страницу приложения требуется пройти процедуру регистрации.
Обсуждение
Поскольку возможности феноменологических уравнений охлаждения трупа в плане термометрической диагностики давности наступления смерти практически исчерпаны, основным методом решения данной научно-практической задачи остаётся конечно-элементный анализ посмертного теплообмена мёртвого тела в целом или его отдельных частей. Однако существенным препятствием на пути внедрения метода конечных элементов остаются его сложность, необходимость приобретения коммерческих пакетов программного обеспечения и наличие у пользователя умений по их использованию, наиболее важными из которых следует считать построение геометрии и сетки КЭМ, задание её теплофизических свойств и краевых условий.
В ходе проведённого исследования впервые разработано открытое судебно-медицинское онлайн-приложение конечно-элементного анализа посмертного теплообмена, позволяющее в режиме real-time рассчитать кривую охлаждения для любой точки мозговой части головы трупа и не требующее от пользователя компетенций по заданию геометрии, теплофизических свойств и краевых условий КЭМ.
Геометрия КЭМ созданного онлайн-приложения является масштабируемой и аппроксимирована многослойным шаром, поскольку при условии адекватной фиксации общих объёмов тканей геометрическое разрешение модели практически не влияет на точность определения давности наступления смерти [5]. Правомерность использованной аппроксимации уже была подтверждена ранее на двумерной КЭМ с грубой сеткой всего из 254 конечных элементов однотипной треугольной формы [8]. Названное упрощение геометрии, выполненное с повышением дискретизации расчётной области до 9277 конечных элементов, позволило значительно уменьшить вычислительные затраты и автоматизировать процесс создания КЭМ, а также сделало возможным выполнение расчётов в онлайн-режиме, повысив при этом их точность по сравнению с двумерными КЭМ.
При наличии соответствующих данных пользователь может задавать в расчётах толщину каждого из пяти актуальных анатомических слоёв: кожно-апоневротического лоскута, свода черепа, лептоменинкса с ликвором, поверхностного и глубокого блоков мозга, а также их основные теплофизические свойства и начальную температуру их внешних границ. Это позволяет учитывать реальную геометрию и размеры расчётной области и составляющих её тканевых структур. Кроме того, доступно задание величины конвективного теплообмена, значение коэффициента теплоотдачи которого зависит от типа и скорости перемещения охлаждающей труп среды [15–18].
Важным преимуществом разработанного онлайн-приложения, как и всех аналогичных программных пакетов, основанных на методе конечных элементов, является возможность расчёта кривой охлаждения для любой точки мозгового отдела головы. Интерфейс онлайн-приложения позволяет указать декартовы координаты сечения головы диаметральной плоскостью. Однако ввиду осевой симметрии КЭМ достаточно принять одну из декартовых координат равной нулю, а вместо второй задать радиальную координату точки, соответствующую расстоянию от центра КЭМ, расположенного в толще межталамического сращения головного мозга, до конца раневого канала от игольчатого термощупа. Возможность учёта в расчётах координат диагностической точки позволяет устранить такой присущий феноменологическим уравнениям тип погрешностей определения давности наступления смерти, как диагностические промахи, связанные с отклонениями конца термощупа от точки с наивысшей температурой при проведении краниоэнцефальной термометрии [19].
Помимо учёта особенностей начального температурного поля, размеров охлаждаемой части тела, интенсивности конвективного теплообмена и координат диагностической точки, онлайн-приложение при определении давности наступления смерти позволяет также принять во внимание изменения внешней температуры в процессе остывания мёртвого тела.
К факторам, ограничивающим применение разработанного онлайн-инструмента определения давности наступления смерти, следует отнести комбинированные условия охлаждения трупа, наиболее существенными из которых являются комбинации конвекции с внешним тепловым потоком, например с интенсивной прямой солнечной радиацией или теплопередачей за счёт теплопроводности при значительном по площади контакте охлаждаемой части тела с другими предметами. Ограничением является также асимметрия конвективного теплообмена, сопровождаемая неоднородностью значений коэффициента теплоотдачи на различных участках поверхности головы. Такая ситуация может возникнуть вследствие одновременного расположения головы трупа в различных средах, например в водной и воздушной, или из-за наличия головного убора.
В целом перечисленные ограничения могут быть устранены, но ценой необходимости возложения на пользователя обязанности самостоятельного задания геометрии и краевых условий КЭМ. По этой причине выбранный разработчиками приложения диапазон учитываемых при вычислениях условий охлаждения и геометрических характеристик расчётной области является компромиссом между возможностями моделирования последних и простотой интерфейса приложения.
Таким образом, изложенные данные позволяют использовать разработанное приложение в судебно-медицинской экспертной практике для онлайн-установления давности наступления смерти в раннем посмертном периоде путём краниоэнцефальной термометрии трупа при отсутствии неоднородности конвективного теплообмена головы или его комбинации с другими типами теплопередачи.
Заключение
Одним из наиболее перспективных современных направлений в области термометрической диагностики давности наступления смерти является метод конечно-элементного анализа посмертного теплообмена, позволяющий преодолеть ограничения, присущие феноменологическим уравнениям охлаждения трупа. Однако доступные пакеты программного обеспечения метода конечных элементов характеризуются высокой стоимостью, а также предполагают наличие у пользователя навыков самостоятельного задания условий решаемых задач и анализа результатов выполненных вычислений.
В ходе проведённого исследования разработано открытое онлайн-приложение конечно-элементного анализа посмертного конвективного теплообмена головы человека, предназначенное для определения давности наступления смерти путём краниоэнцефальной термометрии трупа. Для достижения возможности выполнения вычислений в онлайн-режиме геометрия головы при расчётах аппроксимирована многослойным шаром. В отличие от феноменологических уравнений разработанное онлайн-приложение при определении давности наступления смерти позволяет учитывать особенности начального температурного поля, размеры и теплофизические свойства основных анатомических слоёв головы, интенсивность конвективного теплообмена, координаты диагностической точки и изменения внешней температуры в процессе охлаждения трупа. Результатом работы онлайн-приложения являются визуализация кривых охлаждения в диагностической точке и на поверхности головы в первые сутки посмертного периода с генерацией соответствующих числовых значений, а также вывод информации о геометрии и свойствах сетки КЭМ и распределении температуры вдоль контура интегрирования от центра КЭМ до диагностической точки. История и результаты расчётов сохраняются в личном кабинете пользователя.
В отличие от других программных пакетов разработанный решатель не предъявляет высоких системных требований к компьютеру, поскольку реализован в виде веб-приложения с выполнением всех вычислений на сервере, и не требует наличия у пользователя специальной математико-инженерной подготовки. Изложенное позволяет использовать метод конечных элементов в формате разработанного онлайн-инструмента в судебно-медицинской экспертной практике при установлении давности наступления смерти путём краниоэнцефальной термометрии трупа.
Дополнительная информация
Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении работы.
Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи.
Вклад авторов. Все авторы подтверждают соответствие своего авторства международным критериям ICMJE (все авторы внесли существенный вклад в разработку концепции, проведение работы и подготовку статьи, прочли и одобрили финальную версию перед публикацией). Наибольший вклад распределён следующим образом: В.Г. Недугов ― разработка вычислительного алгоритма метода конечных элементов, написание исходного кода онлайн-приложения, научное редактирование рукописи, рассмотрение и одобрение окончательного варианта рукописи; Г.В. Недугов ― концепция и дизайн работы, разработка геометрии конечно-элементной модели, сбор данных по теплофизическим свойствам тканей головы, тестирование и валидация онлайн-приложения, написание текста рукописи.
Additional information
Funding source. This article was not supported by any external sources of funding.
Competing interests. The authors declare that they have no competing interests.
Authors’ contribution. All authors made a substantial contribution to the conception of the work, acquisition, analysis, interpretation of data for the work, drafting and revising the work, final approval of the version to be published and agree to be accountable for all aspects of the work. V.G. Nedugov: development of the computational algorithm, writing the source code of the online application, scientific editing of the manuscript, review and approval the final version of the manuscript; G.V. Nedugov: concept and design of the work, development of the geometry of the finite element model, collection of data on the thermophysical properties of head tissues, testing and validation of the online application, writing the manuscript.
作者简介
Vladimir G. Nedugiv
Samara National Research University
Email: nedugovvg@gmail.com
ORCID iD: 0009-0007-7542-7235
SPIN 代码: 2407-7937
俄罗斯联邦, Samara
German V. Nedugov
Samara State Medical University
编辑信件的主要联系方式.
Email: nedugovh@mail.ru
ORCID iD: 0000-0002-7380-3766
SPIN 代码: 3828-8091
MD, Dr. Sci. (Medicine), Assistant Professor
俄罗斯联邦, Samara参考
- 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
- 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
- Smart JL, Kaliszan M. Use of a finite element model of heat transport in the human eye to predict time of death. J Forensic Sci. 2013;58(Suppl 1):S69–S77. doi: 10.1111/1556-4029.12022
- 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. EDN: UUZEXN doi: 10.1007/s00414-016-1523-0
- Weiser M, Erdmann B, Schenkl S, et al. Uncertainty in temperature-based determination of time of death. Heat and Mass Transfer. 2018;54(9):2815–2826. EDN: ILDZMW doi: 10.1007/s00231-018-2324-4
- Ullrich J, Weiser M, Subramaniam SJ, et al. The impact of anatomy variation on temperature based time of death estimation. Int J Legal Med. 2023;137(5):1615–1627. EDN: KDNXGI doi: 10.1007/s00414-023-03026-w
- Subramaniam JS, Hubig M, Muggenthaler H, et al. Sensitivity of temperature-based time since death estimation on measurement location. Int J Legal Med. 2023;137(6):1815–1837. EDN: GXPTML doi: 10.1007/s00414-023-03040-y
- Nedugov GV. Estimation of the postmortem interval by the method of finite element modeling of postmortem heat transfer in human head. Sci Innovations Med. 2022;7(3):179–185. EDN: CIMZZD doi: 10.35693/2500-1388-2022-7-3-179-185
- Nelson DA, Nunneley SA. Brain temperature and limits on transcranial cooling in humans: Quantitative modeling results. Eur J Appl Physiol Occup Physiol. 1998;78(4):353–359. EDN: AUNLVL doi: 10.1007/s004210050431
- Zhu L, Diao C. Theoretical simulation of temperature distribution in the brain during mild hypothermia treatment for brain injury. Med Biol Eng Comput. 2001;39(6):681–687. EDN: OVPRTJ doi: 10.1007/BF02345442
- Duck FA. Physical properties of tissue: A comprehensive reference book. London: Academic Press; 1990. P. 9–42.
- Logg A, Wells G, Mardal KA. Automated solution of differential equations by the finite element method: The FEniCS book. Berlin: Springer-Verlag; 2012. doi: 10.1007/978-3-642-23099-8
- Muggenthaler H, Hubig M, Schenkl S, et al. Calibration and parameter variation using a finite element model for death time estimation: The influence of the substrate. Leg Med (Tokyo). 2017;25:23–28. doi: 10.1016/j.legalmed.2016.12.007
- 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
- Clark RP, Toy N. Forced convection around the human head. J Physiol. 1975;244(2):295–302. doi: 10.1113/jphysiol.1975.sp010798
- Defraeye T, Blocken B, Koninckx E, et al. Computational fluid dynamics analysis of drag and convective heat transfer of individual body segments for different cyclist positions. J Biomech. 2011;44(9):1695–1701. doi: 10.1016/j.jbiomech.2011.03.035
- Kurazumi Y, Fukagawa K, Sakoi T, et al. Convective heat transfer coefficient relating to evaluation of thermal environment of infant. Heliyon. 2022;8(12):e12076. EDN: SQKCMK doi: 10.1016/j.heliyon.2022.e12076
- Jiang S, Zhang M, Wang S, Li J. Numerical investigation of the convective heat transfer coefficient for a sleeping infant in a ventilation room. Indoor Air. 2022;32(10):e13126. EDN: ZBNRHA doi: 10.1111/ina.13126
- Vavilov AJ. Diagnostic ‘blunder’ as a cause of errors in the calculated determination of the age of death by the thermal method. Problemy ekspertizy v meditsine. 2008;8(3-4):8–11. (In Russ.) EDN: OKCJBD
补充文件
