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