多孔复合介质周期结构热传导与质扩散问题的多尺度数值方法.pdfVIP

中国工程热物理学会 传热传质学学术会议 编号：993112 多孔复合介质周期结构热传导和质扩散问题的多尺度数值方法 曹礼群 (中国科学院计算数学与科学工程计算所，北京，100080) 罗剑兰 ‘ (湘潭大学物理系，湖南省湘潭市，411105) 摘要 本文针对复杂的多孔复合介质周期结构的热传导和质扩散问题，提出了多尺度 有限元计算公式，它是宏观和细观相结合的数值方法． 关键词：多孔复合介质，热传导，质扩散，多尺度渐近分析，有限元法 1引言 多孔复合介质的热传导和质扩散问题，由于细部几何构造的复杂性和物理参数 的各向异性，研究起来十分困难，曲须对实际问题作一些简化处理，即便如此，一 般而言也无法用解析法求解，只能通过数值方法寻求近似解．然而，如果直接对原 始问题实施有限元法或有限差分法，因为网格生成和线性代数方程组求解的工作量 太大而难以实现．目前工程上普遍采用的一些有效导热系数法，见文[1])能较好地 反映结构的宏观平均特性，但对两相不同介质界面及孔洞附近的温度场奇异性难以 用于多孔复合介质的热传导和质扩散问题，见文[2,3，41，由于均匀化方法也是通过寻 找有效介质得到平均宏观热传导和质扩散方程，它也不能刻划局部物理量的涨落． 当然，它的求解方法与工程上普遍采用的有效导热方法有本质的不同，数学上能证 明原问题的解与均匀化问题的解具有某种收敛性．理论分析与数值实验均表明：多 Analysis 尺度渐近分析方法(MuhiscaleAs》mptoticMethod)对线性定常(稳态)问题 十分有效，见文【2,3．4，5，6，11】，而对非定常问题的讨论，目前理论结果尚不多，许多问 e 题数学上仍是未解决的难题，见[3j．本文对一类多孔复合介质周期结构的热传导和 质扩散问题，给出了具体的多尺度渐近展开公式，并在此基础上设计了有限元算法 格式，本文的结果均有严格的数学汪明． 2．问题的数学描述 根据孔隙的大小，对极细小孔隙实行填充处理，利用统计方法寻找有效导热系数， 其中f为一个基本周期单元(单胞)的细观尺寸，三为结构的宏观尺寸．设。是n维欧氏 c 空问R“中具有l一周期结构的羌界区域，如躅1所示，sloR“是整周期无孔洞区域， ；∈L v～4T 其中Z是整数集，口是单位方体，即0=“：0白1，j=1，2，…n} E∽．总的边界8阱=ru是 构外部边界F=0nI、r’，孔洞表面=f20n 热传导方程： 出，》，争掣一去‰㈦争掣)_F㈦巩‰肌胁(0T) u九口 匿 图1．无界区域。 图2整周期无孔洞区域‰ 图3．有孔洞区域ni 图4求解区域n5=f2iunl t 其中p(z，詈)．c(。，詈)，kij(X．吾)分别袁示密度函数，比热函数和导热系数，r(z 是热源，Fj(。6)是温度场，此处 f、 f，·(詈) z∈ni p‘2’；’=1p。 。∈n。 (2) ，。、 f c(；) 。∈瞒