摘要：利用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等人的理論結果相吻合，對金屬泡沫材料制備工藝的改進和優化具有參考價值。鑒于未考慮由于液膜破裂導致氣泡間的合并而對泡沫尺寸分布的影響，這一假設在泡沫體發展初期是可以接受的，因而，此模型僅適用于泡沫發展初期。