热冲击作用下层合皮肤组织热传导过程数值模拟  

Numerical analysis for non-Fourier thermal behavior of biological multilayer skin tissue subjected to impulsive heat source

作　　者：郭攀 武文华[2] 赵军 李倩 

机构地区：[1]郑州大学力学与工程科学学院工程力学系,微纳成型技术国家级国际联合研究中心,郑州450001 [2]大连理工大学运载工程与力学学部工程力学系,工业装备结构分析国家重点实验室,大连116024 [3]郑州大学力学与工程科学学院,郑州450001

出　　处：《计算力学学报》2017年第6期683-689,共7页Chinese Journal of Computational Mechanics

基　　金：国家自然科学基金(11572072);国家重点基础研究发展计划(973计划)(2014CB046803);国家科技重大专项(2016ZX05028-02-05);国家国际科技合作专项(2015DFA30550)资助项目

摘　　要：在急剧温度变化等强间断温度冲击作用下的生物层合组织非傅里叶热传导分析中,经典时域连续有限元方法(如Newmark等方法)会在波阵面以后的和层合组织界面附近的区域表现出强烈的数值振荡。这类数值振荡会影响问题求解精度,并带来较大不确定性。针对这类现象,本文发展了改进时域间断Galerkin有限元方法,进一步开展了相关问题的数值模拟。其控制方程的基本未知数(温度)及其时间导数在指定时间间隔内假设存在间断且独立插值。在有限元离散列式中引入比例刚度阵人工阻尼,以成功消除波前位置的虚假数值振荡行为。通过算例对比分析,相比Newmark方法和传统间断Galerkin方法,所提出的改进时域间断Galerkin有限元方法较好消除了波前、波后以及组织界面处的数值振荡,有效捕捉了波阵面的间断行为,提高了计算的精度。At present,the traditional time-stepping numerical method,such as Newmark method,produces numerical oscillations at the wave-after,wave-front stage in solving the non-Fourier bio-heat problems of multilayer tissue under high frequency heat load. The spurious non-physical oscillations obviously decrease the accuracy and stability of the simulated results. This paper presents a modified time discontinuous Galerkin finite element method( MDGFEM). The basic unknown variable( temperature) and its time derivative are assumed to be discontinuous and interpolated individually at each time level in time domain respectively. The artificial damping term is induced to eliminate spurious numerical oscillations the wave-front stage. Numerical simulations show the advantages of the present MDGFEM in eliminating the spurious numerical oscillations at the wave-after,wave-front stage and the interface between layers and in catching the discontinuous behavior in wave front surface than that of Newmark and former DGFEM.

关 键 词：热冲击 时域间断 有限元方法 层合皮肤组织 非傅里叶热传导 

分 类 号：TK121[动力工程及工程热物理—工程热物理] O242.21[动力工程及工程热物理—热能工程]