計算機在濃度場中的數(shù)值模擬研究

1、計算機在濃度場中數(shù)值模擬研究蔣新宇 材料科學與工程學院 班級：b130606學號：b13060637摘要：從物理場、濃度場以及數(shù)值模擬的定義出發(fā)，了解計算機模擬在濃度場分 析中的應用，以及就數(shù)值模擬而言，討論其自身的特點。關鍵字：物理場濃度場數(shù)值模擬引言：計算機在多個方面已冇其廣泛的應用，而在材料科學中，對于濃度場的數(shù)值分析, 相關軟件、方程的應川能夠更加精準的計算、分析在一定條件下的物理規(guī)律，并冃最終投入 實際的運用中去。濃度場屬丁物理場的一種。對丁物理場，場指物體在空間中的分 布情況。場是用空間位置函數(shù)來表征的。當溶液間存在濃度差時，在 介而允許溶質(zhì)自由通過的條件下，高濃度側與低濃度側的溶

2、質(zhì)在空間 上的分布是均勻遞減的，此種濃度差在空間上的遞減稱為濃度梯度。 該濃度梯度的作用稱為濃度場。數(shù)值模擬數(shù)值模擬也叫計算機模擬。依靠電子計算機，結合有限元或有限 容積的概念，通過數(shù)值計算和圖像顯示的方法，達到對工程問題和物 理問題乃至自然界各類問題研究的目的。在計算機上實現(xiàn)一個特定的 計算，非常類似于履行一個物理實驗。這時分析人員已跳出了數(shù)學方 程的圈子來對待物理現(xiàn)象的發(fā)生，就像做一次物理實驗。數(shù)值模擬實際上應該理解為用計算機來做實驗。比如某一特定機 翼的繞流，通過計算并將其計算結果在熒光屏上顯示，可以看到流場 的各種細節(jié)：如激波是否存在，它的位置、強度、流動的分離、表面 的壓力分布、受力

3、大小及其隨吋間的變化等。通過上述方法，人們可 以清楚地看到激波的運動、渦的生成與傳播。總之數(shù)值模擬可以形象 地再現(xiàn)流動情景，與做實驗沒有什么區(qū)別。從上面的例子可以看到，數(shù)值模擬包含以下幾個步驟：首先要建 立反映問題（工程問題、物理問題等）木質(zhì)的數(shù)學模型。具體說就是 要建立反映問題各量之間的微分方程及相應的定解條件。這是數(shù)值模 擬的出發(fā)點。沒有正確完善的數(shù)學模型，數(shù)值模擬就無從談起。牛頓 型流體流動的數(shù)學模型就是著名的納維一斯托克斯方程（簡稱方程） 及其相應的定解條件。數(shù)學模型建立之后，需要解決的問題是尋求高效率、高準確度的 計算方法。由于人們的努力，目前已發(fā)展了許多數(shù)值計算方法。計算 方法不僅

4、包括微分方程的離散化方法及求解方法，還包括貼體坐標的 建立，邊界條件的處理等。這些過去被人們忽略或回避的問題，現(xiàn)在 受到越來越多的重視和研究。在確定了計算方法和坐標系后，就可以開始編制程序和進行計 算。實踐表明這一部分工作是整個工作的主體，占絕大部分時間。由 于求解的問題比較復雜，比如方程就是一個非線性的十分復雜的方 程，它的數(shù)值求解方法在理論上不夠完善，所以需要通過實驗來加以 驗證。正是在這個意義上講，數(shù)值模擬又叫數(shù)值試驗。應該指出這部 分工作決不是輕而易舉的。在計算工作完成后，大量數(shù)據(jù)只能通過圖像形象地顯示出來。因 此數(shù)值的圖像顯示也是一項十分重要的工作。目前人們已能把圖作得 像相片一樣逼

5、真。利用錄像機或電影放映機可以顯示動態(tài)過程，模擬 的水平越來越高，越來越逼真。計算機模擬技術在我國工程技術方面的應用已經(jīng)達到一個成熟 的境界。如長江三峽截流，計算機模擬技術為長江三峽截流成功起到 了舉足輕重的作用。在截流前期，科學家們采用計算機模擬的方式將 三峽截流的相關數(shù)據(jù)和資料輸入到計算機中，采用模擬試驗進行截流 模擬，為長江三峽截流的實際成功奠定了良好的基礎。不僅如此，在 后期，三峽工程的專家們?nèi)匀徊捎糜嬎銠C模擬的方式對三峽工程將帶 來的社會效益、經(jīng)濟效益進行了模擬，通過計算機模擬在這方面取得 了可靠地資料，為三峽工程決策提供了重要依據(jù)。黃河小浪底工程、萬家寨引黃工程等，這些工程也都采用

6、了計算 機模擬技術，為其實際工程的運營提供了良好的數(shù)據(jù)依據(jù)。除此之外, 在建筑業(yè)中，采用計算機模擬技術及吋發(fā)現(xiàn)工程中的問題，對其進行 及時修改，以免造成嚴重后果，大大降低了修改成本，同時也推動了 建筑業(yè)的信息化發(fā)展。計算機模擬技術的優(yōu)點在于其具有反復模擬功能，其很容易改變 系統(tǒng)的結構和系數(shù)。在解決一些實際問題時，用分析方法建模的形式 來處理一些帶隨機因素的復雜系統(tǒng)，此吋需要作很多簡化假設，這可 能與實際問題有相當大的差距，最終導致無法應用所得到的答案。這 時，計算機模擬將成為解決這一問題的有效方法；計算機模擬具有較 強的靈活性和任意性，其不受時間、空間、條件的限制，這一優(yōu)勢使 其成為一個解決實

7、際問題的一個強大的工具；計算機模擬只需從實際 數(shù)據(jù)出發(fā)來模仿描述系統(tǒng)，逐步求精，準確地將系統(tǒng)反映出來，此吋 根本不需要對變動機理進行深究；計算機模擬可以建立更高層次的模 型，進行智能化模型處理，為決策提供數(shù)據(jù)支持。當然數(shù)值模擬也存在著一些不足：計算機模擬缺乏統(tǒng)一可信的測 量尺度一一主要是由于模型都是由建模者建立的，其與這些建模者的 水平有很大關系。另外，建模者在觀測的過程中很容易被自身預先設 計大的結果所影響，導致結果主觀性較強，因此導致模擬的結果不能 被大眾所認可。建模者本身的能力是有限的，要想超越系統(tǒng)、超越自 身是非常困難的。這也是計算機模擬的弊端。 =j在固體屮，擴散是物質(zhì)傳輸?shù)奈ㄒ环绞?/p>

8、，通常是指從材料內(nèi)部原 子遷移的微觀過程。以及由于大量原子遷移而造成的物質(zhì)的宏觀流 動。擴散與材料在生產(chǎn)和使用過程中的許多重要的物理化學過程密切 相關，因此對于擴散的濃度場計算具有重要意義。菲克定律，是描述氣體擴散現(xiàn)象的宏觀規(guī)律，這是生理學家菲克 (fick)于1855年發(fā)現(xiàn)的。早在1855年，菲克就提出了：在單位 時間內(nèi)通過垂直于擴散方向的單位截面積的擴散物質(zhì)流量(稱為擴散 通量diffusion flux,用j表示)與該截面處的濃度梯度(concentration gradient)成正比，也就是說，濃度梯度越大，擴散通量越大。這就是 菲克第一定律。菲克第二定律是在第一定律的基礎上推導出來

9、的。菲克第二定 律指出，在非穩(wěn)態(tài)擴散過程中，在距離x處，濃度隨時間的變化率等 于該處的擴散通量隨距離變化率的負值。其中菲克第一定律只適應于j和c不隨時間變化穩(wěn)態(tài)擴散(steady-state diffusion)的場合。對于穩(wěn)態(tài)擴散也可以描述為：在擴 散過程中，各處的擴散組元的濃度c只隨距離x變化，而不隨時間t 變化，每一吋刻從前邊擴散來多少原子，就向后邊擴散走多少原子， 沒有盈虧，所以濃度不隨時間變化。實際上，大多數(shù)擴散過程都是在 非穩(wěn)態(tài)條件下進行的。非穩(wěn)態(tài)擴散(nonsteady-state diffusion)的特點 是：在擴散過程中，j隨時間和距離變化。通過各處的擴散通量j隨 著距離x

10、在變化，而穩(wěn)態(tài)擴散的擴散通量則處處相等，不隨時間而發(fā) 生變化。對于非穩(wěn)態(tài)擴散，就要應用菲克第二定律了。擴散微分方程在形態(tài)上與熱傳導微分方程非常接近，因此也可以 采用與溫度場求解相同的數(shù)值方法來求解。帶氣幕的局部潔凈室數(shù)值模擬因為潔凈室的換氣次數(shù)較大，且氣流組織基本可視為強制對流流 型，所以本次的計算模型可采用標準的高雷諾數(shù)二方程模型。為了簡 化計算，對潔凈室的實際條件進行了如下假設：室內(nèi)氣流流動為穩(wěn)態(tài)流動，室內(nèi)氣流不可壓縮流體，物性為常數(shù), 忽略質(zhì)量力；室內(nèi)無內(nèi)熱源，圍護結構絕熱，對于潔凈室來說，可假 設為無溫差送風，而且將室內(nèi)溫度場視為均勻溫度場；忽略污染粒了 的質(zhì)量，并假定它是被動量對氣流

11、無作用，室內(nèi)污染源的發(fā)塵速率恒 定?？拷诿嫣幉捎脡毫Ρ诿婧瘮?shù)。離散方法采用有限差分法。在劃分網(wǎng)格時，使用交錯網(wǎng)格，且在 氣幕風口號上方設置不均勻風格。在方程組求解時，對耦合方程組使 用simple算法，單個方程組使用adi逐行迭代法。綜合各算例，主送風速度與氣幕風速對不均勻度的作用有所不 同。當主送風速度相同時，氣幕風速（或?qū)挾龋鶆蚨葏^(qū)域而積影 響不大。而當主送風速度增人時，可滿足不均勻度在±20%內(nèi)的區(qū)域 寬度增大。模擬中將點污染源按所處位置分為兩類：按橫向分為主潔凈區(qū)內(nèi) 澎湃和環(huán)境塵源；按高度分為工作區(qū)高度塵源的非工作區(qū)高度塵源。 釆用散發(fā)率2.5x104粒/s的點污染源，這與一般潔凈室設計規(guī)定在 人穿著潔凈服劇烈活動時，全身散塵率為5.6x105粒/min相比還是 比較安全。從模擬結果可以看出：污染場是沿流線分布的，近似看成是二維分布場。丄潔凈區(qū)內(nèi)塵源相同位置時所產(chǎn)生的污染場主要隨主送風風速 發(fā)牛變化。主送風風速增大，污染場濃度略有降低。位置不同時，例 如當污染源靠近潔凈室中線，高濃度都集中于污染源下風側，工作面 以上擴散很小，不影響工作區(qū)潔凈度。當塵源在主潔凈區(qū)內(nèi)時，氣幕寬度和風速對該塵源的污染場基本 無作用。而當環(huán)境塵源位于工作區(qū)高度時，對風速相同的寬口氣幕比 窗口氣幕隔斷效果顯著。污染源的位置對局部