在現(xiàn)代科技的迅速發(fā)展中，磁性液體作為一種特殊的智能材料，因其能在外加磁場下迅速改變自身力學性質(zhì)而受到廣泛關(guān)注，特別是在精密制動系統(tǒng)、震動控制及目標導向藥物輸送系統(tǒng)等領(lǐng)域。然而，傳統(tǒng)的磁性液體數(shù)值模型多依賴于均勻分布的球形顆粒進行模擬，這雖然在初期為我們理解磁性液體的基本性質(zhì)提供了便利，但因其忽略了磁性顆粒本身形狀在磁場作用下可能形成的鏈狀或更復(fù)雜的結(jié)構(gòu)，限制了對磁性液體微觀復(fù)雜行為的準確預(yù)測和深入理解。因此，亟需開發(fā)一種新的模型，能夠準確捕捉非球型磁性顆粒在磁場下的復(fù)雜相互作用和材料本身磁流變性能的變化。
中國科學技術(shù)大學龔興龍教授團隊報告了一種基于修正的分子動力學和磁偶極子相互作用構(gòu)建的針對納米球鏈磁性流體的數(shù)值模型。通過多組模擬，該研究不僅揭示了納米球鏈的長細比和幾何尺寸對宏觀磁流變效應(yīng)的影響，還建立了納米球鏈的幾何參數(shù)與其介觀結(jié)構(gòu)形成之間的聯(lián)系。最終，該研究確定的納米球鏈最佳幾何參數(shù)是長細比為7、粒徑為20 nm。另外，與傳統(tǒng)的球顆粒模型相比，該模型成功地將預(yù)測誤差減少了3.883倍。相關(guān)研究以“Nanosphere Chain Model of Magnetic Fluid and Its Dynamic Performance”為題發(fā)表于最新一期《Langmuir》上，論文第一作者為中國科學技術(shù)大學工程科學學院博士生王康，通訊作者為龔興龍教授和鄧華夏教授。
【數(shù)值模擬方法的構(gòu)建】
與傳統(tǒng)的基于球形顆粒的模型不同，本模型采用納米球鏈作為基礎(chǔ)計算單元。在同一納米球鏈內(nèi)，納米球通過磁偶極相互作用力結(jié)合在一起，忽略了排斥力和范德華力，整個球鏈被視為一個不可破壞的剛體。而在不同納米球鏈之間，則考慮了它們的旋轉(zhuǎn)和橫向碰撞交互。在每次模擬中，通過調(diào)整球鏈中納米球的數(shù)量，來改變球鏈顆粒的長細比。模擬初期，所有納米球鏈均沿z軸排列，隨后在z軸方向施加穩(wěn)定的剪切流，促使結(jié)構(gòu)開始變化，直到整個系統(tǒng)達到動態(tài)平衡。最終，該模型綜合考慮了磁偶極力、范德華力、彈性擠出力和流體粘滯阻力，在x-y平面的上下邊界處應(yīng)用了周期性邊界條件以模擬連續(xù)體。
圖1 模擬的基本設(shè)置參數(shù)和納米球鏈的示意圖
【數(shù)值模型的誤差分析】
為了驗證模型的準確性，分別模擬了含有細長的DTM/MS2磁性顆粒和棒狀顆粒的磁性流體。DTM/MS2磁性顆粒的平均長度被設(shè)定為631 nm，長細比為5；棒狀顆粒的平均長度被設(shè)定為5 μm，長細比為8。模擬在300K溫度和7%體積分數(shù)下進行，磁場強度從0.1到10 kOe進行變化。對于棒狀顆粒，在1.1、2.2和3.3 kOe的磁場下進行了模擬。此外，還額外構(gòu)建了一個基于球形顆粒的傳統(tǒng)模型進行比較。為了保持模擬體積和質(zhì)量分數(shù)，模型中的顆粒數(shù)量與鏈模型相同，其他參數(shù)也保持與鏈模型相同，唯一改變的是模型的計算單元，一個將顆粒視為球形，另一個將顆粒視為鏈條。兩個模型的模擬結(jié)果與文獻報告的實驗結(jié)果進行了比較，球鏈模型成功地將平均預(yù)測誤差減少了3.883倍。
圖2 模型的誤差分析
【模擬參數(shù)對材料磁流變性能的影響】
通過修改球鏈的長細比，該工作首先研究了長細比對材料磁流變性能的影響。如圖3a所展示，長細比在2至10的范圍內(nèi)只有一個峰值。當長細比達到7時，基于納米球鏈的磁性流體顯示出了最大的剪切應(yīng)力。根據(jù)長細比不同，球鏈可分為三種類型：短鏈（長細比2至4）、中長鏈（長細比5至7）和長鏈（長細比大于7）。在短鏈區(qū)間，隨著長細比的增加，剪切應(yīng)力穩(wěn)定增長，其中直徑20nm的球鏈剪切應(yīng)力增長最顯著，平均剪切應(yīng)力從235Pa提升至562Pa，增幅達139%。同時，介觀結(jié)構(gòu)的傾斜角度減少，其徑向分布函數(shù)緩慢上升（圖3c）。對于中長鏈，剪切應(yīng)力的增加速度減緩，甚至略有下降，而徑向和傾斜角度分布函數(shù)的增長緩慢。而在長鏈情況下，剪切應(yīng)力開始降低，平均傾斜角度也隨之減少。當長細比為7時，球鏈的傾斜角度最接近理想傾斜角度（圖3b），此時納米球鏈達到了最大剪切應(yīng)力。此外，不同長細比對平均分子勢能的影響也進行了分析，結(jié)果顯示不同長細比下平均分子勢能沒有顯著變化（圖3d）。
圖3 顆粒長細比對宏觀性能和介觀結(jié)構(gòu)參數(shù)的影響
同時分析顆粒長細比對自身屬性的影響時發(fā)現(xiàn)，隨著長細比的增加，磁偏轉(zhuǎn)力矩呈現(xiàn)近線性增長，而旋轉(zhuǎn)速度以二次函數(shù)曲線減。▓D4a）。在低長細比下，旋轉(zhuǎn)速度大于剪切速率。隨著模擬的進行，具有更強磁偏轉(zhuǎn)力矩的納米球鏈會構(gòu)建更垂直的微觀結(jié)構(gòu)。在高長細比下，剪切速率強于旋轉(zhuǎn)速度。此時，具有較大旋轉(zhuǎn)速度的納米球鏈有更小傾斜角度的更高概率。同時，流體分層現(xiàn)象也直接影響材料的宏觀性質(zhì)。為了將分層現(xiàn)象進行定量分析，模擬區(qū)域被平均劃分為100個子區(qū)域。通過計算每個子區(qū)域的體積分數(shù)，得到磁性流體的分層程度參數(shù)。如圖4b所示，隨著長細比的增加，分層程度線性增加。高度分層程度會導致局部體積分數(shù)過小從而引起材料磁流變效應(yīng)的降低，這可以解釋當長細比超過7時剪切應(yīng)力的減小。
圖4 顆粒長細比對顆粒自身力學參數(shù)和介觀結(jié)構(gòu)分層程度的影響
之后，通過修改顆粒的粒徑大小，對不同幾何尺寸的磁性液體顆粒進行了研究。如圖5a所示，隨著納米球的尺寸增加，剪切應(yīng)力增加并在長細比1和2時達到穩(wěn)定值。在較大長細比的情況下，剪切應(yīng)力在20nm處達到峰值，然后下降至穩(wěn)定值。通過觀察平均傾斜角度的變化圖（圖5b），可以看出平均傾斜角度隨著納米球直徑的增加而減小。這個結(jié)果對于20nm和30nm的納米球鏈來說，是最接近理論最佳角度的。圖5c和5d顯示，隨著直徑的增加，聚集度逐漸降低？梢钥闯，三個峰值在20nm時達到最大，證明了此時的微觀結(jié)構(gòu)具有良好的緊湊性。在20nm時，由于傾斜角度和微觀結(jié)構(gòu)的緊湊性都處于良好狀態(tài)，剪切應(yīng)力達到了最大值，與靜止狀態(tài)相比有103%的提升（從322Pa提高到654Pa）。
圖5 顆粒粒徑對宏觀性能和介觀結(jié)構(gòu)參數(shù)的影響
最終，文章還探討了剪切速率對材料磁流變性能的影響。納米球鏈的剪切應(yīng)力隨剪切率的增加而增加，顯示出典型賓漢流體的性質(zhì)，響應(yīng)時間則隨剪切率的增加而減少。介觀結(jié)構(gòu)方面，納米球鏈在不同剪切率下均為柱狀結(jié)構(gòu)，顆粒的平均傾斜角度隨速率的增大而增大，該變化趨勢與剪切應(yīng)力保持一致，而粒子的聚集程度則受到剪切速率的影響較小。
圖6 剪切速率對宏觀性能和介觀結(jié)構(gòu)參數(shù)的影響
本研究成功開發(fā)了一種新型的基于納米球鏈的顆粒動力學模型，深入討論了納米磁性顆粒自身幾何參數(shù)對宏觀力學性能和介觀結(jié)構(gòu)形成的影響。通過深入探索納米球鏈的長細比和幾何尺寸對磁性液體性能的影響，本研究揭示了優(yōu)化這些參數(shù)以提升磁流變效應(yīng)的可能性。最后，本研究確定了最優(yōu)的長細比和幾何尺寸參數(shù)，為設(shè)計高性能磁流變材料提供了重要的參考。
原文鏈接
Kang Wang, Bing Liu, Xinyu Lian, Shouhu Xuan, Huaxia Deng and Xinglong Gong. Nanosphere Chain Model of Magnetic Fluid and Its Dynamic Performance.Langmuir, 2024, 40, 6187-6197.
DOI: 10.1021/acs.langmuir.3c03538
https://doi.org/10.1021/acs.langmuir.3c03538
