摘要：利用MacPherson等人的理论分析计算三维液态铝合金泡沫体在不同的泡沫结构模型、不同的表面张力下泡沫内气泡体积分布的演化过程。模拟结果表明，球型泡沫体平均体积近似以v∝ t1. 482方式增长，与MacPherson等人的理论结果相符合。结果表明，泡沫模型对泡沫的胀大、缩小有直接的影响。正六面体型泡沫体平均体积随时间减小，而正二十面体型的泡沫体平均体积随时间增长。表面张力是加速泡沫体演化的重要因素。

由于良好的热学、力学、声学性能及可回收再利用等优点，泡沫金属材料，尤其以泡沫铝合金为代表，已成为21世纪材料科学与技术的研究热点。固态金属泡沫体内部孔隙率、胞元尺寸大小、分布及拓扑结构对泡沫金属的各项性能有直接的影响。然而泡沫金属制备中，泡沫体内的金属相一般要经历由液态向固态转变的过程，故液态金属泡沫体的孔隙率、胞元尺寸及分布将影响到凝固后泡沫金属成品的各项指标参数。所以，很有必要研究液态金属泡沫体的演化过程。目前，大量研究集中于泡沫的析液现象。Gergely等研究了液态泡沫铝的析液过程，他们基于等径五边形十二面体模型计算了液态金属铝泡沫体内孔隙率随时间的变化，分析了表面张力、泡径大小等因素对析液过程的影响。Brunke等人通过实验研究了泡沫铝合金的内部泡沫尺寸分布以及析液过程，并考虑了气泡的合并。Magrabi等人基于Lemlich的泡沫胀缩理论研究了压缩空气泡沫体中泡沫尺寸随时间的分布。

研究大多数基于等径泡沫体模型，但实际生产过程中泡沫体气泡尺寸大小并不均匀，在一定的尺度范围内存在一定的分布规律。1952年，著名的数学家Von Neumann揭开了二维气泡演化的规律，证明二维情况下描述气泡的表面积随时间变化与其边数相关：A∝ ta, a = 1。其中，A为气泡表面积，t为时间。由于该问题在数学以及许多科学及工程部门都具有十分重要的意义，半个多世纪以来，科学家都力图将Von Neumann的结论推广到三维空间。但直到2007年，美国数学家MacPherson和材料学家Stolovitz才解决了这一问题，成为受到学术界广泛关注的重大突破。他们指出：泡沫是大量的单个小空间区域相互作用形成的，这些小区域不断变化——张大、收缩或者合并。其研究还表明，在三维空间中，气泡不同表面之间交界边缘的总长度如果超过平均宽度的6倍，那么气泡将会膨胀；反之，气泡则会收缩。

泡沫铝的孔径大小、分布对于泡沫铝的性能有着至关重要的影响。在生产过程中，由于过程的高温及不可视性，给研究带来了相当大的难度。为了突破这一局限，本文基于美国数学家MacPherson和材料学家Stolovitz在2007年取得的重大理论突破——三维空间泡沫演化理论，对三维液态铝合金泡沫体内气泡尺寸分布的演化过程进行数值模拟与深度分析。本文将重点探讨在不同泡沫结构模型和不同表面张力条件下，泡沫内气泡体积分布的演变规律，特别是揭示表面张力在其中的核心驱动作用。

理论基础：从二维到三维的泡沫演化规律

泡沫演化的数学描述是一个跨越了半个多世纪的难题。早在1952年，著名数学家Von Neumann揭示了二维气泡的演化规律：在二维情况下，气泡表面积随时间的变化仅与其边数相关。这一发现为后续研究奠定了基础，但将这一规律推广到三维空间却耗费了科学家们数十年的努力。

直到2007年，MacPherson和Stolovitz才成功解决了这一难题。他们指出，泡沫是由大量单个小空间区域相互作用形成的复杂系统，这些小区域处于不断的变化之中——张大、收缩或合并。他们的核心结论是：在三维空间中，气泡不同表面之间交界边缘的总长度如果超过平均宽度的6倍，气泡将会膨胀；反之，则会收缩。这一理论为三维泡沫的尺寸分布演化提供了坚实的数学框架。

在本文的研究中，作者采用了多面体泡沫模型，并基于以下几点假设进行简化分析：

1. 泡沫体内气体为理想气体，液体为不可压缩流体；

2. 不考虑相邻气泡间的合并现象，气泡总数保持不变；

3. 不考虑气体在气液界面的扩散现象；

4. 系统保温良好，过程中温度不变，铝合金始终处于液态。

为了描述泡沫尺寸分布随时间的演化，研究定义了一个气泡尺寸（体积）分布函数。该函数表示某一时刻泡沫体内部气泡体积在特定范围内的总个数。由于过程中气泡总数恒定，由此建立守恒方程。在此基础上，根据MacPherson的理论，定义了气泡体积随时间的变化率，该变化率与表面张力、单个气泡Plateau边界长度之和以及欧拉特征长度密切相关。

其中，欧拉特征长度与气泡的拓扑结构直接相关。研究列举了正四面体、正六面体、正八面体、正十二面体和正二十面体五种不同形状气泡模型的欧拉特征数。例如，正四面体具有6条Plateau边界，欧拉特征长度为182倍的Plateau边界长度；而正十二面体和正二十面体均有30条边界，但欧拉特征长度分别为529倍和348倍的Plateau边界长度。这些结构参数的差异，直接决定了泡沫演化的不同走向。

表面张力：泡沫演化的核心驱动力

在液态铝合金泡沫的演化过程中，当孔隙率大于95%时，析液量可以忽略不计；同时，对于液态金属泡沫，由气体扩散造成的尺寸变化也被排除在外。在这种情况下，表面张力成为了影响泡沫尺寸分布演化过程的绝对主导因素。

表面张力在气液界面上起着至关重要的作用。改变铝合金的化学成分，会导致气液界面上的表面张力发生变化。数值模拟结果清晰地表明，表面张力的大小直接决定了泡沫演化的速度。

具体而言，在相同的初始泡沫尺寸分布下，表面张力越大，泡沫平均体积的变化速度越快，尤其是在泡沫演化的后期，这种加速效应更为显著。然而，需要明确的是，表面张力虽然极大地加速了泡沫演化的进程，但它并不直接决定泡沫平均体积是增加还是减小。决定泡沫平均体积增减的根本原因在于泡沫的拓扑结构模型。

这一发现对于理解实际生产中的泡沫行为具有重要意义。表面张力作为驱动气泡合并、收缩或膨胀的能量来源，其数值的微小变化都可能通过欧拉特征长度等拓扑参数被放大，从而影响整个泡沫体系的稳定性与最终结构。

泡沫结构模型对演化的决定性影响

泡沫的拓扑结构是决定其演化方向的内在因素。不同的多面体模型，由于Plateau边界条数和欧拉特征长度的不同，会导致截然不同的演化结果。研究通过对比球形（近似为五边形十二面体）、正六面体以及正二十面体等模型，揭示了这一规律。

球形泡沫体（五边形十二面体模型）

以五边形十二面体模型近似代替球形气泡，这是研究中最为详尽的部分。计算结果表明，此类泡沫体的平均体积随时间增长，增长方式近似符合幂律关系。研究指出，泡沫平均体积的增长指数约为1.48至1.50之间，这与MacPherson等人的理论结果高度吻合。

在演化过程中，泡沫尺寸分布沿着尺寸增大的方向发展。随着时间的推移，气泡尺寸分布的区间逐渐变宽，分布函数的峰值逐渐变小。这意味着小气泡在相对缩小，而大气泡在持续胀大，整体分布向大尺寸方向偏移。这种演化趋势在不同的初始泡沫尺寸分布下均保持一致，证明了该模型下泡沫胀大的内在稳定性。

正六面体泡沫体模型

与球形模型截然相反，正六面体泡沫体的平均体积随时间减小。在正六面体模型下，Plateau边界条数为12，欧拉特征长度为3倍的Plateau边界长度。数值模拟显示，此类泡沫体的尺寸分布演化过程比十二面体模型更快，但方向却是收缩的。

这种收缩特性意味着，如果泡沫中的气泡主要呈现正六面体拓扑结构，那么整个泡沫体系将趋于塌陷或致密化。这解释了为什么在某些制备条件下，泡沫无法稳定存在而迅速衰变。

正二十面体泡沫体模型

正二十面体泡沫体同样具有30条Plateau边界，但其欧拉特征长度为348倍的Plateau边界长度。与正六面体不同，正二十面体泡沫体的平均体积随时间增长。这表明，面数较多的拓扑结构更倾向于膨胀和粗化。

大小气泡的拓扑差异与双峰分布现象

在实际的泡沫铝合金中，气泡尺寸并不均一。研究提出了一个极具洞察力的观点：小气泡和大气泡可能具有不同的泡沫结构。小气泡由于面数较少，可能接近正六面体型；而大气泡面数较多，可能接近正二十面体。Kraynik等人的研究也支持这一观点，指出泡沫体中小气泡的面数可以是四面、五面或六面体，而大气泡的面数甚至可达170个。

基于这一假设，研究构建了混合模型：对于等效直径在1毫米以下的小气泡采用正六面体结构模型，对于等效直径大于1毫米的大气泡采用正二十面体结构模型。模拟结果出现了有趣的现象：在一段演化时间后，泡沫孔径大小分布出现了双峰。这一结果不仅解释了“小气泡缩小、大气泡胀大”的实验观察，也与Brunke等人近期通过实验证明的泡沫铝合金制备过程中会出现双峰或多峰分布的结论相吻合。

这种双峰分布的产生，正是不同拓扑结构对演化方向产生相反影响的直接证据。小气泡因结构倾向于收缩而逐渐减小，大气泡因结构倾向于膨胀而不断增大，两者共同作用导致了分布曲线的分裂。

数值模拟方法与初始条件

为了准确模拟泡沫尺寸分布的演化，研究采用了韦伯分布函数来近似泡沫初始尺寸分布。韦伯分布能够灵活描述不同形状的单峰分布，通过调整参数可以适应不同工艺条件下的初始状态。

在求解过程中，研究将复杂的偏微分方程转化为有限差分格式进行数值求解。对于不同的泡沫体模型，定义了不同的函数表达式来反映体积变化率与体积的关系。例如，在十二面体模型中，体积变化率与体积的1/3次方成正比；而在正六面体模型中，相关参数则呈现负值，导致收缩行为。

模拟中使用的液态铝合金物性参数均取自文献，确保了计算的可靠性。通过改变表面张力数值和初始分布参数，研究系统地考察了这些因素对演化过程的影响，为后续分析提供了丰富的数据支持。

演化过程中的关键现象与机制分析

平均体积增长的幂律特征

研究特别关注了泡沫平均体积随时间的变化规律。Lemlich的研究表明，二维泡沫体平均尺寸的增长方式符合幂律关系，指数因泡沫体物性而异。对于三维泡沫体，研究同样认为平均体积的增长方式为幂律形式。

通过数值求解，发现球形泡沫体的平均体积增长指数在1.48至1.50之间。这一数值显著高于二维情况下的典型值，反映了三维空间中拓扑约束的复杂性。这种幂律增长表明，泡沫演化是一个自相似过程，表面张力驱动的粗化机制在长时间尺度上维持着特定的统计规律。

表面张力的加速效应

如前所述，表面张力是加速泡沫演化的重要因素。在不同表面张力条件下，泡沫平均体积随时间的变化曲线显示，表面张力越大，曲线斜率越大，演化越快。这种加速效应在演化后期尤为明显，可能是因为随着气泡体积增大，表面能的变化更为显著。

然而，表面张力并不直接决定演化方向。无论是增长还是收缩，表面张力都起着“催化剂”的作用，而真正的“方向盘”则是泡沫的拓扑结构。这一结论纠正了以往可能过于简化地将表面张力视为唯一决定因素的观点。

双峰分布的理论解释

实际生产中观察到的“小气泡缩小、大气泡胀大”现象，长期以来缺乏合理的理论解释。以往的观点多将其归因于气体在气泡间的扩散，但对于液态金属泡沫而言，这种扩散机制是否为主要原因尚无定论。

本研究通过引入拓扑结构差异，提供了一个全新的解释框架：小气泡和大气泡由于尺寸不同，其稳定存在的多面体形状也不同。小气泡倾向于形成面数较少的几何形状（如正六面体），这种结构在表面张力作用下趋于收缩；大气泡倾向于形成面数较多的形状（如正二十面体），这种结构则趋于膨胀。因此，即使不考虑气体扩散，仅由拓扑结构差异就能导致双峰分布的出现。

这一理论不仅与Brunke等人的实验观察一致，也为控制泡沫孔径分布提供了新的思路：通过调整合金成分或工艺参数，改变气泡的拓扑结构，可能实现对泡沫分布的主动调控。

结论

液态金属泡沫演化过程极其复杂，气泡间的合并机理以及泡沫的稳定性与气泡拓扑结构的关系至今都没有一个通用、合理的规律可循。金属泡沫体内部泡沫结构、尺寸分布等由于实验条件(高温、不可视)的限制也很难准确测得，所以，理论及模拟研究十分必要。

借助MacPherson等人的理论，研究分析了金属泡沫体内泡沫尺寸(体积)分布的演化过程，对不同表面张力、不同气泡拓扑结构条件下泡沫体的演化进行了计算模拟，结果表明，在忽略泡间气体扩散和析液作用的前提下，表面张力是泡沫尺寸演化的主要因素。由于泡沫体内部存在大小不一的气泡，而大小气泡的拓扑结构模型有可能不同，从而导致泡沫的演化过程存在差别，大气泡的平均体积增长而小气泡则减小。这些都与MacPherson等人的理论结果相吻合，对金属泡沫材料制备工艺的改进和优化具有参考价值。鉴于未考虑由于液膜破裂导致气泡间的合并而对泡沫尺寸分布的影响，这一假设在泡沫体发展初期是可以接受的，因而，此模型仅适用于泡沫发展初期。