多孔材料渗透系数性能设计.docVIP

多孔材料渗透系数性能设计 徐胜利1 程耿东1 贾海昌2 1大连理工大学工程力学系，工业装备结构分析国家重点实验室，116023 2Department of Mechanical and Aerospace Engineering, Rutgers University Email: shengli.xu.xu@ 摘要：本文将流体拓扑优化技术应用到具有周期性微结构多孔材料等效渗透系数性能设计，建立了相应优化问题的数学模型和灵敏度分析方法。宏观多孔材料中的流体流动服从Darcy定律，微观单胞内流体流动服从Stokes方程。Darcy定律中的渗透系数联系了多孔材料宏观流动和材料单胞微观流动，多孔材料的渗透系数采用均匀化方法对单胞的分析得到。本文采用逆均匀化方法设计固体材料在材料单胞中的分布，得到了一种具有各向同性最大渗透系数性能的多孔材料微结构形式。 关键词：渗透系数，多孔介质，拓扑优化 1 引言 多孔介质材料是自然界和工业界普遍存在的载体形式，多孔介质中质量动量能量传递问题涉及到许多领域，是构成众多自然现象的基本过程，包括核反应堆冷却，多孔介质强化传热传质，石油开采，植物与土壤间的水分养分传递等。H.Darcy[1]于1856年提出了著名的Darcy定律，Darcy定律描述了一定条件下多孔介质中流体流动的规律。Darcy方程中的渗透系数反映了多孔介质对其内部流体的输运能力，是多孔介质一个重要的物理参数。渗透系数的大小是由多孔介质本身决定的，它与多孔介质内部空腔的大小，几何形状以及分布形式有关。多孔介质的许多应用领域都对其渗透性能提出了要求，因此设计具有特定渗透性能的多孔介质材料的微结构形式是一项有意义的工作。 设计具有特定性能的材料结构一直是科研人员研究的热点，一直以来，材料的力学性能，热性能得到了广泛的研究。近来材料的声学性能，电磁学性能以及材料多功能性设计成为研究的热点[2]。材料性能设计可以归结为在一定设计域上的材料布局问题，拓扑优化方法是进行材料性能设计的有力工具[3]。连续体拓扑优化的目的是设计在一定区域上如何布置材料，通常可以描述为实体材料和空心材料（无材料）两相材料的布局设计。著名学者Bendsoe和Kikuchi[4]首先研究了最大化刚度设计下材料的布局问题。丹麦学者Sigmund首先将拓扑优化方法和均匀化方法用于设计桁架类单胞的力学性能[5]。单胞是这类具有周期性微结构材料的最小单位，利用均匀化方法对单胞进行分析可以得到这类具有周期性微结构材料的等效力学性能。逆均匀化方法是通过设计材料在单胞中的分布使得由单胞组成的周期性材料具有指定性能的方法。逆均匀化方法成功的应用到了指定热弹性性能，极端弹性和热膨胀性能[6]，压电复合材料[7]的单胞拓扑设计。近来，逆均匀化方法在热电耦合[8]，刚度和流体耦合[9]等多功能材料设计上得到了应用。 流体流道拓扑设计是拓扑优化新的应用领域。最近，拓扑优化被用来设计Stokes流的流体流道分布以减小流体流动的功耗[10]。流体拓扑优化技术也被用来设计中低雷诺数的流体流道[11]。Guest[9]利用逆均匀化技术研究了考虑渗流性能的多孔材料设计，得到了在一定材料体分比下具有最大各向同性渗透系数周期性多孔介质的单胞微结构形式。Guest采用了速度场和压力场同阶插值有限元求解Stokes方程，提出了用Darcy流近似模拟设计单胞上固体－流体边界，给出了三种固体材料初始密度分布形式：1，固体材料密度从单胞中心向四周逐渐增大；2，固体材料密度从单胞中心向四周逐渐减小；3，固体材料密度均匀分布。对于第三种密度初值分布方式，为了得到材料的拓扑结构，需要固定特定位置的速度值，不同初始密度分布得到的单胞拓扑形式是不同的。文章[9]指出对于随机分布的密度初值，经过几次迭代之后各单元密度趋于均匀分布，无法得到单胞拓扑。 本文提出了一种材料渗透系数性能设计的优化模型，并且给出了相应的伴随灵敏度分析方法。采用混合Galerkin有限元求解相应的控制方程。计算结果表明上面提到的三种密度分布初值比较容易得到单胞的拓扑形式，对于随机分布密度初值我们也进行了计算，并且得到了相同的优化结果，这说明本文提出的优化模型不依赖于设计变量初值。本文的结构如下：第二部分给出了计算等效渗透系数的均匀化方法；第三部分介绍了优化的数学模型；第四部分是相应的灵敏度分析方法；第五部分是算例；最后是结论。 2 Stokes流的均匀化方法 由于多孔介质具有多尺度的特点，直接模拟流体在多孔介质流道里的流动是非常困难的。一般采用Darcy定律等宏观等效方法描述多孔介质内部流体流动。通常多孔介质的渗透系数是不知道的，必须通过实验或者计算的方法得到。采用均匀化的办法，我们可以通过分析微观尺度上的一个单胞得到宏观多孔介质等效的渗透系数。S