一维Stefan问题之有限差分解.PDFVIP

東海科學 19 卷:17-47 17 民國一⃝六年七月出版 一維Stefan 問題之有限差分解 邱靖宏* 黃皇男* 摘 要 本文以有限差分法計算一維二相Stefan 問題的數值解，文中分別採用顯式、Crank-Nicolson 與隱式差分 法來進行，除推導數值計算的穩定條件外，並討論這三種方法對於此類移動邊界問題數值解的精確度以及計 算效率等。由於隱式差分算法為無條件穩定，本文進一步採用此法，以模擬邊界是固定熱源之 Stefan 問題， 而此熱源分別為單邊或雙邊週期時，週期與振幅對於冰水交界面移動的影響。 一、前 言 隨著科學技術和工業快速發展，許多領域中不斷出現具有移動邊界(moving boundary) 的問題，如 物質結晶的成長控制、冰川的溶解- 凝固，化學反應的氣體擴散... 等等。由於這類問題中都包含 有待定的移動介面，而且介面形狀與位置的決定就是所的求解之一部分，使此問題的求解比解一 般的偏微分方程更難。這類移動邊界問題通常稱為 Stefan 問題，最早由J. Stefan 在 1890 年提出 [1]，主要是探討冰在水中融解或是水固化結冰的問題。在一些文獻當中，Stefan 指出在一個半無窮 (semi-infinite) 的液體固化的過程中，其固化的邊界的速率是一固定速率，其溫度 可以表示成下 列形式 ( ) 其中、 與 均為常數，固化的速率 等於。Stefan 在 1889 年提出了此問題的通解， 1959 年Carslaw 和Jaeger 根據Stefan 的結果提出了此問題的Neumann 的解[1, 2, 3]。 現在我們考慮Stefan 問題中比較常見的例子- 冰融解成水的例子。它是一個邊界會隨時間移動 的熱傳導方程(heat equation)，其邊界就是冰與水的交界面。我們的目的就是要觀察當冰融解成水 時，水之溫度 隨著位置 和時間 的變化，並令 。在時間 ，溫度 , 也就是整個區域都是攝氏零度的冰；而當 時，在 處有一個 *東海大學應用數學系 18 穩定熱源開始加熱，即 ，則冰吸熱之後會開始融解成水，其中溶解成水的區域之溫度 必須滿足熱傳導方程式。假設 是表示冰與水的交界面之方程式，此交界面（邊界）的移動 （融化）必會滿足 ，且交界面溫度條件為 ，此條件稱為Stefan 條 件，其中係數 稱為Stefan 數。 如果冰的初始溫度在融點以下，則熱傳導現象除了在水的區域會發生外，冰的區域也會發生，此 類問題稱為兩相(two phase) 問題。其中水的溫度為 ，冰的溫度為 ，並設 在時間 時，還未加熱，溫度為冰的初始溫度 ，當在 處以一穩定熱源加熱時，即 。冰及水均會開始發生熱傳導現象，假設 是代表冰與水的交介面方程式，它會滿足 ，同時介面溫度條件為