摘要：利用MacPherson等人的理論分析計(jì)算三維液態(tài)鋁合金泡沫體在不同的泡沫結(jié)構(gòu)模型、不同的表面張力下泡沫內(nèi)氣泡體積分布的演化過(guò)程。模擬結(jié)果表明，球型泡沫體平均體積近似以v∝ t1. 482方式增長(zhǎng)，與MacPherson等人的理論結(jié)果相符合。結(jié)果表明，泡沫模型對(duì)泡沫的脹大、縮小有直接的影響。正六面體型泡沫體平均體積隨時(shí)間減小，而正二十面體型的泡沫體平均體積隨時(shí)間增長(zhǎng)。表面張力是加速泡沫體演化的重要因素。

由于良好的熱學(xué)、力學(xué)、聲學(xué)性能及可回收再利用等優(yōu)點(diǎn)，泡沫金屬材料，尤其以泡沫鋁合金為代表，已成為21世紀(jì)材料科學(xué)與技術(shù)的研究熱點(diǎn)。固態(tài)金屬泡沫體內(nèi)部孔隙率、胞元尺寸大小、分布及拓?fù)浣Y(jié)構(gòu)對(duì)泡沫金屬的各項(xiàng)性能有直接的影響。然而泡沫金屬制備中，泡沫體內(nèi)的金屬相一般要經(jīng)歷由液態(tài)向固態(tài)轉(zhuǎn)變的過(guò)程，故液態(tài)金屬泡沫體的孔隙率、胞元尺寸及分布將影響到凝固后泡沫金屬成品的各項(xiàng)指標(biāo)參數(shù)。所以，很有必要研究液態(tài)金屬泡沫體的演化過(guò)程。目前，大量研究集中于泡沫的析液現(xiàn)象。Gergely等研究了液態(tài)泡沫鋁的析液過(guò)程，他們基于等徑五邊形十二面體模型計(jì)算了液態(tài)金屬鋁泡沫體內(nèi)孔隙率隨時(shí)間的變化，分析了表面張力、泡徑大小等因素對(duì)析液過(guò)程的影響。Brunke等人通過(guò)實(shí)驗(yàn)研究了泡沫鋁合金的內(nèi)部泡沫尺寸分布以及析液過(guò)程，并考慮了氣泡的合并。Magrabi等人基于Lemlich的泡沫脹縮理論研究了壓縮空氣泡沫體中泡沫尺寸隨時(shí)間的分布。

研究大多數(shù)基于等徑泡沫體模型，但實(shí)際生產(chǎn)過(guò)程中泡沫體氣泡尺寸大小并不均勻，在一定的尺度范圍內(nèi)存在一定的分布規(guī)律。1952年，著名的數(shù)學(xué)家Von Neumann揭開(kāi)了二維氣泡演化的規(guī)律，證明二維情況下描述氣泡的表面積隨時(shí)間變化與其邊數(shù)相關(guān)：A∝ ta, a = 1。其中，A為氣泡表面積，t為時(shí)間。由于該問(wèn)題在數(shù)學(xué)以及許多科學(xué)及工程部門(mén)都具有十分重要的意義，半個(gè)多世紀(jì)以來(lái)，科學(xué)家都力圖將Von Neumann的結(jié)論推廣到三維空間。但直到2007年，美國(guó)數(shù)學(xué)家MacPherson和材料學(xué)家Stolovitz才解決了這一問(wèn)題，成為受到學(xué)術(shù)界廣泛關(guān)注的重大突破。他們指出：泡沫是大量的單個(gè)小空間區(qū)域相互作用形成的，這些小區(qū)域不斷變化——張大、收縮或者合并。其研究還表明，在三維空間中，氣泡不同表面之間交界邊緣的總長(zhǎng)度如果超過(guò)平均寬度的6倍，那么氣泡將會(huì)膨脹；反之，氣泡則會(huì)收縮。

泡沫鋁的孔徑大小、分布對(duì)于泡沫鋁的性能有著至關(guān)重要的影響。在生產(chǎn)過(guò)程中，由于過(guò)程的高溫及不可視性，給研究帶來(lái)了相當(dāng)大的難度。為了突破這一局限，本文基于美國(guó)數(shù)學(xué)家MacPherson和材料學(xué)家Stolovitz在2007年取得的重大理論突破——三維空間泡沫演化理論，對(duì)三維液態(tài)鋁合金泡沫體內(nèi)氣泡尺寸分布的演化過(guò)程進(jìn)行數(shù)值模擬與深度分析。本文將重點(diǎn)探討在不同泡沫結(jié)構(gòu)模型和不同表面張力條件下，泡沫內(nèi)氣泡體積分布的演變規(guī)律，特別是揭示表面張力在其中的核心驅(qū)動(dòng)作用。

理論基礎(chǔ)：從二維到三維的泡沫演化規(guī)律

泡沫演化的數(shù)學(xué)描述是一個(gè)跨越了半個(gè)多世紀(jì)的難題。早在1952年，著名數(shù)學(xué)家Von Neumann揭示了二維氣泡的演化規(guī)律：在二維情況下，氣泡表面積隨時(shí)間的變化僅與其邊數(shù)相關(guān)。這一發(fā)現(xiàn)為后續(xù)研究奠定了基礎(chǔ)，但將這一規(guī)律推廣到三維空間卻耗費(fèi)了科學(xué)家們數(shù)十年的努力。

直到2007年，MacPherson和Stolovitz才成功解決了這一難題。他們指出，泡沫是由大量單個(gè)小空間區(qū)域相互作用形成的復(fù)雜系統(tǒng)，這些小區(qū)域處于不斷的變化之中——張大、收縮或合并。他們的核心結(jié)論是：在三維空間中，氣泡不同表面之間交界邊緣的總長(zhǎng)度如果超過(guò)平均寬度的6倍，氣泡將會(huì)膨脹；反之，則會(huì)收縮。這一理論為三維泡沫的尺寸分布演化提供了堅(jiān)實(shí)的數(shù)學(xué)框架。

在本文的研究中，作者采用了多面體泡沫模型，并基于以下幾點(diǎn)假設(shè)進(jìn)行簡(jiǎn)化分析：

1. 泡沫體內(nèi)氣體為理想氣體，液體為不可壓縮流體；

2. 不考慮相鄰氣泡間的合并現(xiàn)象，氣泡總數(shù)保持不變；

3. 不考慮氣體在氣液界面的擴(kuò)散現(xiàn)象；

4. 系統(tǒng)保溫良好，過(guò)程中溫度不變，鋁合金始終處于液態(tài)。

為了描述泡沫尺寸分布隨時(shí)間的演化，研究定義了一個(gè)氣泡尺寸（體積）分布函數(shù)。該函數(shù)表示某一時(shí)刻泡沫體內(nèi)部氣泡體積在特定范圍內(nèi)的總個(gè)數(shù)。由于過(guò)程中氣泡總數(shù)恒定，由此建立守恒方程。在此基礎(chǔ)上，根據(jù)MacPherson的理論，定義了氣泡體積隨時(shí)間的變化率，該變化率與表面張力、單個(gè)氣泡Plateau邊界長(zhǎng)度之和以及歐拉特征長(zhǎng)度密切相關(guān)。

其中，歐拉特征長(zhǎng)度與氣泡的拓?fù)浣Y(jié)構(gòu)直接相關(guān)。研究列舉了正四面體、正六面體、正八面體、正十二面體和正二十面體五種不同形狀氣泡模型的歐拉特征數(shù)。例如，正四面體具有6條Plateau邊界，歐拉特征長(zhǎng)度為182倍的Plateau邊界長(zhǎng)度；而正十二面體和正二十面體均有30條邊界，但歐拉特征長(zhǎng)度分別為529倍和348倍的Plateau邊界長(zhǎng)度。這些結(jié)構(gòu)參數(shù)的差異，直接決定了泡沫演化的不同走向。

表面張力：泡沫演化的核心驅(qū)動(dòng)力

在液態(tài)鋁合金泡沫的演化過(guò)程中，當(dāng)孔隙率大于95%時(shí)，析液量可以忽略不計(jì)；同時(shí)，對(duì)于液態(tài)金屬泡沫，由氣體擴(kuò)散造成的尺寸變化也被排除在外。在這種情況下，表面張力成為了影響泡沫尺寸分布演化過(guò)程的絕對(duì)主導(dǎo)因素。

表面張力在氣液界面上起著至關(guān)重要的作用。改變鋁合金的化學(xué)成分，會(huì)導(dǎo)致氣液界面上的表面張力發(fā)生變化。數(shù)值模擬結(jié)果清晰地表明，表面張力的大小直接決定了泡沫演化的速度。

具體而言，在相同的初始泡沫尺寸分布下，表面張力越大，泡沫平均體積的變化速度越快，尤其是在泡沫演化的后期，這種加速效應(yīng)更為顯著。然而，需要明確的是，表面張力雖然極大地加速了泡沫演化的進(jìn)程，但它并不直接決定泡沫平均體積是增加還是減小。決定泡沫平均體積增減的根本原因在于泡沫的拓?fù)浣Y(jié)構(gòu)模型。

這一發(fā)現(xiàn)對(duì)于理解實(shí)際生產(chǎn)中的泡沫行為具有重要意義。表面張力作為驅(qū)動(dòng)氣泡合并、收縮或膨脹的能量來(lái)源，其數(shù)值的微小變化都可能通過(guò)歐拉特征長(zhǎng)度等拓?fù)鋮?shù)被放大，從而影響整個(gè)泡沫體系的穩(wěn)定性與最終結(jié)構(gòu)。

泡沫結(jié)構(gòu)模型對(duì)演化的決定性影響

泡沫的拓?fù)浣Y(jié)構(gòu)是決定其演化方向的內(nèi)在因素。不同的多面體模型，由于Plateau邊界條數(shù)和歐拉特征長(zhǎng)度的不同，會(huì)導(dǎo)致截然不同的演化結(jié)果。研究通過(guò)對(duì)比球形（近似為五邊形十二面體）、正六面體以及正二十面體等模型，揭示了這一規(guī)律。

球形泡沫體（五邊形十二面體模型）

以五邊形十二面體模型近似代替球形氣泡，這是研究中最為詳盡的部分。計(jì)算結(jié)果表明，此類(lèi)泡沫體的平均體積隨時(shí)間增長(zhǎng)，增長(zhǎng)方式近似符合冪律關(guān)系。研究指出，泡沫平均體積的增長(zhǎng)指數(shù)約為1.48至1.50之間，這與MacPherson等人的理論結(jié)果高度吻合。

在演化過(guò)程中，泡沫尺寸分布沿著尺寸增大的方向發(fā)展。隨著時(shí)間的推移，氣泡尺寸分布的區(qū)間逐漸變寬，分布函數(shù)的峰值逐漸變小。這意味著小氣泡在相對(duì)縮小，而大氣泡在持續(xù)脹大，整體分布向大尺寸方向偏移。這種演化趨勢(shì)在不同的初始泡沫尺寸分布下均保持一致，證明了該模型下泡沫脹大的內(nèi)在穩(wěn)定性。

正六面體泡沫體模型

與球形模型截然相反，正六面體泡沫體的平均體積隨時(shí)間減小。在正六面體模型下，Plateau邊界條數(shù)為12，歐拉特征長(zhǎng)度為3倍的Plateau邊界長(zhǎng)度。數(shù)值模擬顯示，此類(lèi)泡沫體的尺寸分布演化過(guò)程比十二面體模型更快，但方向卻是收縮的。

這種收縮特性意味著，如果泡沫中的氣泡主要呈現(xiàn)正六面體拓?fù)浣Y(jié)構(gòu)，那么整個(gè)泡沫體系將趨于塌陷或致密化。這解釋了為什么在某些制備條件下，泡沫無(wú)法穩(wěn)定存在而迅速衰變。

正二十面體泡沫體模型

正二十面體泡沫體同樣具有30條Plateau邊界，但其歐拉特征長(zhǎng)度為348倍的Plateau邊界長(zhǎng)度。與正六面體不同，正二十面體泡沫體的平均體積隨時(shí)間增長(zhǎng)。這表明，面數(shù)較多的拓?fù)浣Y(jié)構(gòu)更傾向于膨脹和粗化。

大小氣泡的拓?fù)洳町惻c雙峰分布現(xiàn)象

在實(shí)際的泡沫鋁合金中，氣泡尺寸并不均一。研究提出了一個(gè)極具洞察力的觀(guān)點(diǎn)：小氣泡和大氣泡可能具有不同的泡沫結(jié)構(gòu)。小氣泡由于面數(shù)較少，可能接近正六面體型；而大氣泡面數(shù)較多，可能接近正二十面體。Kraynik等人的研究也支持這一觀(guān)點(diǎn)，指出泡沫體中小氣泡的面數(shù)可以是四面、五面或六面體，而大氣泡的面數(shù)甚至可達(dá)170個(gè)。

基于這一假設(shè)，研究構(gòu)建了混合模型：對(duì)于等效直徑在1毫米以下的小氣泡采用正六面體結(jié)構(gòu)模型，對(duì)于等效直徑大于1毫米的大氣泡采用正二十面體結(jié)構(gòu)模型。模擬結(jié)果出現(xiàn)了有趣的現(xiàn)象：在一段演化時(shí)間后，泡沫孔徑大小分布出現(xiàn)了雙峰。這一結(jié)果不僅解釋了“小氣泡縮小、大氣泡脹大”的實(shí)驗(yàn)觀(guān)察，也與Brunke等人近期通過(guò)實(shí)驗(yàn)證明的泡沫鋁合金制備過(guò)程中會(huì)出現(xiàn)雙峰或多峰分布的結(jié)論相吻合。

這種雙峰分布的產(chǎn)生，正是不同拓?fù)浣Y(jié)構(gòu)對(duì)演化方向產(chǎn)生相反影響的直接證據(jù)。小氣泡因結(jié)構(gòu)傾向于收縮而逐漸減小，大氣泡因結(jié)構(gòu)傾向于膨脹而不斷增大，兩者共同作用導(dǎo)致了分布曲線(xiàn)的分裂。

數(shù)值模擬方法與初始條件

為了準(zhǔn)確模擬泡沫尺寸分布的演化，研究采用了韋伯分布函數(shù)來(lái)近似泡沫初始尺寸分布。韋伯分布能夠靈活描述不同形狀的單峰分布，通過(guò)調(diào)整參數(shù)可以適應(yīng)不同工藝條件下的初始狀態(tài)。

在求解過(guò)程中，研究將復(fù)雜的偏微分方程轉(zhuǎn)化為有限差分格式進(jìn)行數(shù)值求解。對(duì)于不同的泡沫體模型，定義了不同的函數(shù)表達(dá)式來(lái)反映體積變化率與體積的關(guān)系。例如，在十二面體模型中，體積變化率與體積的1/3次方成正比；而在正六面體模型中，相關(guān)參數(shù)則呈現(xiàn)負(fù)值，導(dǎo)致收縮行為。

模擬中使用的液態(tài)鋁合金物性參數(shù)均取自文獻(xiàn)，確保了計(jì)算的可靠性。通過(guò)改變表面張力數(shù)值和初始分布參數(shù)，研究系統(tǒng)地考察了這些因素對(duì)演化過(guò)程的影響，為后續(xù)分析提供了豐富的數(shù)據(jù)支持。

演化過(guò)程中的關(guān)鍵現(xiàn)象與機(jī)制分析

平均體積增長(zhǎng)的冪律特征

研究特別關(guān)注了泡沫平均體積隨時(shí)間的變化規(guī)律。Lemlich的研究表明，二維泡沫體平均尺寸的增長(zhǎng)方式符合冪律關(guān)系，指數(shù)因泡沫體物性而異。對(duì)于三維泡沫體，研究同樣認(rèn)為平均體積的增長(zhǎng)方式為冪律形式。

通過(guò)數(shù)值求解，發(fā)現(xiàn)球形泡沫體的平均體積增長(zhǎng)指數(shù)在1.48至1.50之間。這一數(shù)值顯著高于二維情況下的典型值，反映了三維空間中拓?fù)浼s束的復(fù)雜性。這種冪律增長(zhǎng)表明，泡沫演化是一個(gè)自相似過(guò)程，表面張力驅(qū)動(dòng)的粗化機(jī)制在長(zhǎng)時(shí)間尺度上維持著特定的統(tǒng)計(jì)規(guī)律。

表面張力的加速效應(yīng)

如前所述，表面張力是加速泡沫演化的重要因素。在不同表面張力條件下，泡沫平均體積隨時(shí)間的變化曲線(xiàn)顯示，表面張力越大，曲線(xiàn)斜率越大，演化越快。這種加速效應(yīng)在演化后期尤為明顯，可能是因?yàn)殡S著氣泡體積增大，表面能的變化更為顯著。

然而，表面張力并不直接決定演化方向。無(wú)論是增長(zhǎng)還是收縮，表面張力都起著“催化劑”的作用，而真正的“方向盤(pán)”則是泡沫的拓?fù)浣Y(jié)構(gòu)。這一結(jié)論糾正了以往可能過(guò)于簡(jiǎn)化地將表面張力視為唯一決定因素的觀(guān)點(diǎn)。

雙峰分布的理論解釋

實(shí)際生產(chǎn)中觀(guān)察到的“小氣泡縮小、大氣泡脹大”現(xiàn)象，長(zhǎng)期以來(lái)缺乏合理的理論解釋。以往的觀(guān)點(diǎn)多將其歸因于氣體在氣泡間的擴(kuò)散，但對(duì)于液態(tài)金屬泡沫而言，這種擴(kuò)散機(jī)制是否為主要原因尚無(wú)定論。

本研究通過(guò)引入拓?fù)浣Y(jié)構(gòu)差異，提供了一個(gè)全新的解釋框架：小氣泡和大氣泡由于尺寸不同，其穩(wěn)定存在的多面體形狀也不同。小氣泡傾向于形成面數(shù)較少的幾何形狀（如正六面體），這種結(jié)構(gòu)在表面張力作用下趨于收縮；大氣泡傾向于形成面數(shù)較多的形狀（如正二十面體），這種結(jié)構(gòu)則趨于膨脹。因此，即使不考慮氣體擴(kuò)散，僅由拓?fù)浣Y(jié)構(gòu)差異就能導(dǎo)致雙峰分布的出現(xiàn)。

這一理論不僅與Brunke等人的實(shí)驗(yàn)觀(guān)察一致，也為控制泡沫孔徑分布提供了新的思路：通過(guò)調(diào)整合金成分或工藝參數(shù)，改變氣泡的拓?fù)浣Y(jié)構(gòu)，可能實(shí)現(xiàn)對(duì)泡沫分布的主動(dòng)調(diào)控。

結(jié)論

液態(tài)金屬泡沫演化過(guò)程極其復(fù)雜，氣泡間的合并機(jī)理以及泡沫的穩(wěn)定性與氣泡拓?fù)浣Y(jié)構(gòu)的關(guān)系至今都沒(méi)有一個(gè)通用、合理的規(guī)律可循。金屬泡沫體內(nèi)部泡沫結(jié)構(gòu)、尺寸分布等由于實(shí)驗(yàn)條件(高溫、不可視)的限制也很難準(zhǔn)確測(cè)得，所以，理論及模擬研究十分必要。

借助MacPherson等人的理論，研究分析了金屬泡沫體內(nèi)泡沫尺寸(體積)分布的演化過(guò)程，對(duì)不同表面張力、不同氣泡拓?fù)浣Y(jié)構(gòu)條件下泡沫體的演化進(jìn)行了計(jì)算模擬，結(jié)果表明，在忽略泡間氣體擴(kuò)散和析液作用的前提下，表面張力是泡沫尺寸演化的主要因素。由于泡沫體內(nèi)部存在大小不一的氣泡，而大小氣泡的拓?fù)浣Y(jié)構(gòu)模型有可能不同，從而導(dǎo)致泡沫的演化過(guò)程存在差別，大氣泡的平均體積增長(zhǎng)而小氣泡則減小。這些都與MacPherson等人的理論結(jié)果相吻合，對(duì)金屬泡沫材料制備工藝的改進(jìn)和優(yōu)化具有參考價(jià)值。鑒于未考慮由于液膜破裂導(dǎo)致氣泡間的合并而對(duì)泡沫尺寸分布的影響，這一假設(shè)在泡沫體發(fā)展初期是可以接受的，因而，此模型僅適用于泡沫發(fā)展初期。