有限差分法。它是計算機數值模擬最早采用的方法,至今仍被廣泛運用。
考慮時間因素的影響,差分格式可以分為顯格式、隱格式、顯隱交替格式等。其基本思想是用離散的、只含有限個未知數的差分方程去代替連續變量的微分方程和定解條件,求出差分方程的解作為求偏微分方程的近似解。該方法將求解域劃分為差分網格,用有限個網格節點代替連續的求解域,以泰勒級數展開等方法,把控制方程中的導數用網格節點上的函數值的差商代替進行離散,建立以網格節點上的值為未知數的代數方程組,進而獲得所需的數值解。
到目前為止,國內外專家學者已采用有限差分法對多種典型形狀食品的凍結時間做了大量數值模擬,取得了許多成果。Allada和Quan(1966),Bailey
(1974),Mannapperuma和Singh(1988)先后在以食品熱物理性質為溫度的函數的假設前提下,推導出了一些預測食品凍結時間的顯式差分公式;Cleland和Earle
(1984)確定了計算凍結時間和針對各種相變問題的有限差分解決方案,認為三階隱式有限差分的Lees(1966)式模擬預測最為精確;Sheen和Hayakawa(I990)
應用交替隱式差分法和有限容積法,建立了一種適用于食品不同導熱過程的新型有限差分模型來求解凍結時間;關志強(1999)等通過建立平板狀食品凍結過程的溫度模型和焓模型,應用完全隱式差分法獲得凍結時間的數值解,并以牛肉為例實驗驗證:劉寶林(2000)用隱式差分法求解了無限大平板食品的一維速凍過程,討論了介質溫度、表面傳熱系數和食品厚度對凍結時間的影響;謝晶等
(2001)利用有限差分的焓方程數值求解鱈魚在肋板鼓風凍結裝置中的凍結時間,計算結果與實驗的誤差小于10%;汪政高(Zhengfu Wang)等(2007)利用VB編程,有限差分離散模型方程,其模擬與實驗結果的相對誤差分別僅為4.69%(平板狀卡爾斯魯厄測試物)、8.03%(平板狀牛肉)、0.55%(球狀牛肉)、5.61%
(圓柱狀卡爾斯魯厄測試物)。
有限差分法是一種直接將微分問題變為代數問題的近似數值解法,數學概念直觀,表達簡單,是發展較早且比較成熟的數值方法。它特別適合處理形狀規則食品的邊界條件及線性間題,但對復雜區域的適應性較差及數值解的守恒性難以保證。