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