有限差分适应性参数法解二阶常微分方程式之探讨.docVIP

第十五屆全國計算流體力學學術研討會　　　　　　　　　　　　　　　　　　高雄，中華民國九十七年八月 The 15th National Computational Fluid Dynamics Conference　　 　　　　 　Kaohsiung,　August,2008 PAGE 有限差分適應性參數法解二階常微分方程式之探討 Adaptive Parameter Scheme of Finite Difference Method in Solving Second Order Ordinary Differential Equations 李志良 郭鴻森 德霖技術學院機械工程系助理教授 大同大學機械工程系教授 摘要 本研究利用有限差分適應性參數法探討含非齊次項之二階常微分方程式，研究首先將有限差分通式代入微分方程式中，所得差分方程式中的非齊次項與各格點係數則含有一適應性參數。利用對原微分方程式加以微分的增階動作，如此可以將截斷誤差轉換至高階，且同時可得到差分式的修正項。有限差分適應性參數法之適應性參數的最佳值可以由原微分方程式的齊次解析解推導得來，但是在本研究中由於左右差分格點係數的表示式於尋求適應性參數最佳值時不能一致而無法求得。本研究最後以三個算例探討，研究結果顯示，本研究之數值解比中央差分法準確但是有穩定性的問題。研究中發現若將有限差分法與有限解析法結合，針對二階常微分方程式可發展出具有兩個適應性參數的有限差分法。 關鍵字: 有限差分適應性參數法，有限差分通式，有限解析法 一、前言 由於現代工程的問題愈行複雜，所以其中大部份問題已經無法單靠物理直覺以及過去經驗來解決。雖然以實驗處理問題的方式，在過去一直很成功地解決許多問題，但是一但遭遇到高速、高負荷、高溫或其他不正常狀況時則告失敗。而且實驗工作不但複雜又費時昂貴，此時數值方法的產生及應用，便可以用電腦計算模擬工程問題，以降低尋找答案的工作量與費用。而微分方程式在此便有其基本重要性，因為許多物理定律與關係在數學上常以微分方程式的形式出現。而微分方程式的數值方法至今一直在發展當中，有限差分法便是其中被人常用的一種數值方法。 目前在有限差分法中將微分式變為差分式的差分方法已經很多，以一階導數的差分方法為例，有以差分格點求曲線斜率的方式，如中央差分法Central Difference Scheme (CDS)、上風差分法Upwind Difference Scheme (UDS)以及Leonard Difference Scheme [1-4]等等。又Leonard [5]亦利用差分格點先求出前後各半格的格點之表示式，再以中央差分法求出一階導數差分式，即Quadratic Upstream Interpolation for Convection Kinematics (QUICK) 的方法。此外甚至有用有限元素法導出的一階導數差分式[6]，如Quadratic Upwiind Petro-Galerkin scheme (QUPG) 的方法，還有Leonard 又將LDS 及LUDS 兩種方法利用Weighting factor method 組成LSUDS 的一階導數差分式[1]等等，而且新的差分方法一直在被研究中。然而這麼多的差分方法中，對於其截斷誤差的探討研究卻並不是十分清楚，一般僅知道其截斷誤差的階數，這也是為什麼會發展出那麼多差分方法的原因，以及希望能盡量減少誤差達到更準確的地步。上述各相關文獻之研究，大多著重於各階導數之差分展開式其截斷誤差的盡量減低，而非考量整個統御方程差分式之截斷誤差的有效降低，尤其是在非齊次之微分方程式時。Adaptive Parameter Scheme適應性參數有限差分法[7]，便是一種可以補回整個差分式截斷誤差的方法，本研究即應用此法探討二階非齊次之微分方程式。 二、公式推導 一般二階非齊次微分方程式之形式如下： (1) 在此先以常係數問題探討，對(1)式中之一階與二階導數以三個差分格點描述之有限差分通式[8]做差分式處理，且在此特別強調零階導數依循傳統有限差分法貫用的作法，不做任何差分式的處理，差分式的推導如下。在此將以三點差分格點描述之一階與二階導數有限差分通式代入微分方程式(1) 式中，且其中一階導數係數與零階導數係數分別為常數A與B，其初步得到的差分方程式如下： (2) 整理(2)式後可變為： (3) 且其差分方程式的完整截斷誤差為： (4) 接著依適應性參數法處理截斷誤差的步驟，將微分方程式(1)式加以微分，如此能得到增階之關係式，再將其代入(4)式的截斷誤差中，則可消去二階截斷誤差。(1)式經微分後結果為： (5)