（计算数学专业论文）抛物型方程的几种可并行的有限差分方法.pdf

山东大学博士学位论文 抛物型方程的几种可并行的有限差分方法 张守慧 ( 山东大学数学与系统科学学院，济南，2 5 0 1 0 0 ) 摘要 抛物型方程是偏微分方程中基本方程之一在自然科学的众多领域中，许多现象 是用抛物型方程或者方程组描述的 1 】例如热传导以及其它扩散现象、化学反应、粒 子的运输等等另外在一些问题的数值处理中也经常出现抛物型偏微分方程在现代数 值方法中，最早为人们所注意且理论分析完善的是有限差分法，因此抛物型偏微分方 程的有限差分方法一直是人们关心的焦点随着并行机的问世和发展，传统的有限差 分方法在不同方面暴露出各自的弱点例如，古典显式虽然适合于并行计算，但它是条 件稳定的，特别是多维问题中计算步长受到严格的限制；古典隐式和c r a n k - n i c o l s o n 格 式是绝对稳定的，但需要求解联立方程组，不便于直接在并行机上应用因此需要构造 具有良好稳定性、并行性和计算精度的新的差分方法 七十年代初，m i r a n k e r 2 l 指出用有限差分逼近偏微分方程时，主要是组织传统差分 方法的并行实现，至于设计新算法，推动力是很小的；之后的十几年高阶差分格式方面 的研究得到了发展p 7 】八十年代初，上述情况由于e v a n s 和a b d u l l a h 的工作 8 1 2 】 而发生了变化，他们设计的分组显式方法保证了数值计算的稳定性，同时由于显式求解 而使该方法具有很好的并行性质它是不同类型s a u l y e v 非对称格式 1 3 】的恰当组合 由于不同的s a u l y e v 格式的截断误差中某些项绝对值相等，符号相反，在同一时间层 和不同时间层上连续交替使用不同的非对称格式，可带来截断误差的部分抵消，从而提 高方法的计算精度这些非对称格式都是隐格式，但由于它们之间的巧妙结合，可以显 式求解，这就是e v a n s - a b d u l l a h 分组显式( g e ) 这项工作说明了建立满足上述要求的 新的差分格式是可能的但是在将分组显式思想应用于变系数问题时，稳定性的证明 遇到了困难 在此基础上，张宝琳等在 1 4 1 6 】中提出利用s a u l y e v 非对称格式构造分段隐式 的思想，并恰当的使用交替技术建立了多种显一隐式和纯隐式交替并行方法，取得了稳 定性和并行兼顾的研究成果之后又将方法推广到变系数问题，并用能量法证明了方 法的绝对稳定性在数值试验中发现，分段或分块并行计算的结果一般都比原来相应的 未加分裂时的结果精确所以通过分而治之的策略来建立新算法，不但可以用于并行， 还可以提高精度之后涌现出大量的并行差分算法的研究成果，韩臻在【1 7 ，1 8 】中详 山东大学博士学位论文 细研究了一类纯显一隐分段和分块交替方法；【1 9 2 1 】中冯慧等通过不同点的隐式差分 格式之间的相互约化来建立新型迭代方法，此方法和j a c o b i 方法同样具有并行性，却 比j a c o b i 收敛快【2 2 】中张志跃等给出变系数抛物型问题的分组显式方法，并用能量 方法给出稳定性证明( 2 3 - 2 6 】中王文洽等针对不同的方程建立了分段的显隐格式，证 明方法的稳定性并给出数值算例上述的方法在并行性和稳定性方面都有其优良的表 现，但是他们都存在一个共同的问题，那就是它们都是基于二阶差分格式建立的，这 直接影响数值计算中空间的误差精度近年来，研究人员开始致力于研究高阶紧致差 分格式 2 7 】中s a n j i v ak l e l e 提出了高阶的紧致差分格式，文中对格式的误差做了 f o u r i e r 分析，并将它与经典差分格式做了比较【2 8 3 2 】中m a r kh c a r p e n t e r 等针对 不同问题提出高阶紧致格式并给出理论分析将这些高阶差分格式与交替分组思想结 合起来，是否可以得到稳定性好，可并行且精度高的数值算法呢? 近年来，涌现出大量 的高阶交替分组格式的研究工作【4 垂5 2 】 本文作者在王文洽教授的精心指导下，就抛物型问题的几类数学模型利用有限差 分方法的技巧，构造了具有良好数值性质和计算效果的迭代方法、分组显式方法和交 替分段方法，对方法做了理论分析并给出算例说明方法的适用性本人拓广了前人的工 作，不具有重复性本文共分为五章 第一章中主要利用 1 9 】中冯慧提出的数值s t e n c i l 的概念，将其应用于二维对流 扩散方程，建立了比j a c o b i 迭代收敛快的新型迭代算法本章首先给出针对对流扩 散方程的数值s t e n c i l 的概念，经过三次消元过程得到最终的数值s t e n c i l ，在此基础 上建立了新型迭代算法；通过分析迭代误差证明了方法的收敛性，并与j a c o b i 迭代比 较收敛阶；最后数值试验说明方法的适用性，证实了理论分析的结论本章内容已被 ( ( i n t e r n a t i o n a lj o u r n a lo fc o m p u t e rm a t h e m a t i c s ) ) 接受 第一章的创新之处在于将 1 9 】中的方法应用到含有时间项的高维抛物型问题中， 建立了收敛速度快、具有并行性质的新型迭代格式，通过分析迭代误差证明了方法的 收敛性以及与古典迭代法之间收敛阶的比较；最后用实际例子说明了算法的有效性 第二章主要运用 2 2 ，2 4 】中的构造思想，将中心差分格式与分组显式思想相结合， 针对含有变系数的对流扩散方程建立了分组显式方法，并用能量方法证明了该格式的 稳定性本章首先给出基于c r a n k - n i c o l s o n 差分格式的四种非对称的逼近方程，通过 它们的巧妙组合建立交替分组显式方法；由于扩散项为变系数，所以采用能量法证明 稳定性；最后数值试验说明方法的适用性本章内容已投到( ( i n t e r n a t i o n a lj o u r n a lo f c o m p u t e rm a t h e m a t i c s ) ) 第二章的创新之处在于对于变系数的抛物问题给出和c r a n k - n i c o l s o n 格式相匹配 的交替分组显格式，并用能量方法证明了稳定性 接下来的三章内容中，主要借鉴了 1 4 ，2 3 2 6 ，4 6 4 8 】中交替分组( 段) 格式的思想， 山东大学博士学位论文 将它与高阶差分格式 2 7 - 3 2 】相结合，建立了高阶的交替格式；经证明方法都是绝对稳 定的，且具有并行性质；在时间步长足够小时，空间的局部截断误差可达到o ( h 4 ) 第三章中首先给出高阶的显、隐差分格式，在隐格式的基础上构造了四种非对称 格式，通过它们之间的巧妙组合建立了交替分段显隐格式；由 3 3 3 4 】中的k e l l o g g 引 理证明了方法的无条件稳定性；得到了方法的局部截断误差；数值算例证实了方法的实 用性，并且可以达到o ( h 4 ) 的误差精度本章内容已被计算物理接受 第三章的创新之处在于将高阶差分格式与交替分段显隐格式的思想相结合，构造 出高阶的交替分段显隐格式方法具有良好的数值稳定性，空间误差阶可以达到o ( h 4 ) 阶 第四、五章是在第三章的基础上，引入高阶c r a n k - n i c o l s o n 差分格式并适当变形， 构造了八个非对称逼近方程；通过交替使用这些差分格式建立了两种不同的数值计算 方法 第四章是单独应用八个非对称差分格式构造了交替分组显格式，通过k e l l o g g 引 理证明方法的稳定性，通过误差分析得到两层抵消部分误差后的误差可以达到o ( t h ) 阶；数值试验说明了方法的实用性，并且时间步长充分小的前提下误差对于空间来说可 以达到四阶本章内容发表在山东大学学报( 理学版) 第五章利用交替分段格式思想，引入了在非对称格式之间插入对称格式的思想，从 而建立交替分段c r a n kn i c o l s o n 格式本格式同样具有数值稳定性、可并行性质，误差 分析时由于插入中心对称格式，所以在这些点处两层之间部分抵消后误差较小，这在数 值试验中得到了证实交替分段c r a n k - n i c o l s o n 格式的结果比相应的交替分组显格式 结果要好本章内容已投稿到应用数学与力学 后两章的创新之处在于将高阶差分格式与交替分组、交替分段c r a n kn i c o l s o n 思 想充分的结合，建立与第三章不同的非对称s a u l y e v 格式，在此基础上构造相应的交 替方法这些方法都具有绝对稳定性、可并行性质；并且在空间上截断误差可以达到 o ( h 4 ) 阶；最后给出数值算例说明方法的适用性 关键词：差分格式；数值s t e n c i l ；迭代法；交替分组格式；并行；稳定性；收敛阶 山东大学博士学位论文 so m ep a r a l l e l i s mn u m e r i c a lm e t h o d s f o rp a r a bo l i ce q u a t i o n s h o u h u iz h a n g s c h o o lo fm a t h e m a t i c sa n ds y s t e ms c i e n c e ，s h a n d o n gu n i v e r s i t y ， j i n a n ，s h a n d o n g ，2 5 0 10 0 , p r c h i n a a b s t r a c t p a r a b o l i ce q u a t i o ni so n eo fb a s i cp a r t i a le q u a t i o n s i nm a n ys c i e n c ef i e l d s ，m a n y p h e n o m e n aa r ed e s c r i b e db yp a r a b o l i ce q u a t i o n ( s ) 1 ，s u c ha st h ep r o c e s so fh e a tc o n - d u c t i o na n dd i f f u s i o n ，t h ec h e m i c a lr e a c t i o ne t c a m o n gm o d e r nn u m e r i c a lm e t h o d s ， t h ef i n i t ed i f f e r e n c em e t h o di st h ee a r l i e s ta n dm o s tp e r f e c tm e t h o d s ot h ef i n i t e d i f f e r e n c em e t h o df o rs o l v i n gp a r a b o l i ce q u a t i o ni sa l w a y saf o c a lw h i c hp e o p l e sc a r e a b o u t a st h ep a r a l l e lc o m p u t e rc o m e si n t ob e i n ga n dd e v e l o p s ，s o m ed i s a d v a n t a g e d i s a p p e a r si nd i f f e r e n tm e a n s f o re x a m p l e ，t h ec l a s s i c a le x p l i c i ts c h e m ei ss u i tf o r p a r a l l e lc o m p u t i n gb u ti t sc o n d i t i o n a ls t a b i l i t y e s p e c i a l l yf o rh i g h d i m e n s i o np r o b - l e m ，t h et i m es t e pi sl i m i t e dv e r ys e v e r e l y t h ec l a s s i c a li m p l i c i ta n dc r a n k - n i c o l s o n s c h e m ei sa b s o l u t e l ys t a b l e ，b u tt h e yc a nb es o l v e do n l yb ys o l v i n gl i n e a re q u a t i o n s o b v i o u s l yt h e ya r en o ts u i tf o rp a r a l l e lc o m p u t i n g s oi t sw o r t h yt oc o n s t r u c t i n go t h e r n e wd i f f e r e n c em e t h o d sw h i c hh a sb e t t e rs t a b i l i t y , p a r a l l e l i s ma n dh i g h - p r e c i s i o n i nt h ee a r l ys e v e n t i e s ，m i r a n k e r 2 】p o i n t e do u tt h a to r g a n i z i n gt h et r a d i t i o n a l d i f f e r e n c em e t h o di no r d e rt op a r a l l e lc o m p u t i n gi st h em a i nm e t h o dw h e nw ea p p r o x - i m a t e dt h ep a r t i a le q u a t i o nb yf i n i t ed i f f e r e n c em e t h o d b e t w e e nt h es e v e n t i e s ，t h e r e s e a r c hw a sm a i n l ya b o u th i g h o r d e rd i f f e r e n c es c h e m ef o rd i f f e r e n te q u a t i o n s 3 7 1 b u ts i n c et h ee i g h t i e s ，t h es i t u a t i o nc h a n g e db e c a u s eo fe v a n sa n da b d u l l a h sw o r k 8 1 2 i nt h ee a r l ye i g h t i e s ，e v a n sa n da b d u l l a hp r o p o s e dt h ei d e aw h i c hc o n s t r u c t e d g r o u pe x p l i c i tm e t h o db ya p p r o p r i a t ec o m b i n a t i o no fd i f f e r e n ts a u l y e va s y m m e t r i c s c h e m e 1 3 t h eg r o u pe x p l i c i tm e t h o dk e e p st h es t a b i l i t yo fn u m e r i c a lc o m p u t i n g ， a n dh a sb e t t e rp a r a l l e l i s mb e c a u s ei tc a nb es o l v e de x p l i c i t l y b e c a u s es o m et e r m si n t h et r u n c a t i o ne r r o ro fd i f f e r e n ts a u l y e vs c h e m ei se q u a lf o rt h e i ra b s o l u t ev a l u ea n d t h es i g hi sc o n t r a s t ，m a k i n gu s eo ft h e ma l t e r n a t i n gi nat i m el a y e ro rd i f f e r e n tl a y e r x 山东大学博士学位论文 m a yc a n c e ls o m et r u n c a t i o ne r r o ra n dt h ec a l c u l a t i o na c c u r a c yc a nb ei m p r o v e d a n d t h e s es a u l y e vs c h e m ew e r ei m p l i c i t ，b u tt h eg r o u ps c h e m ec a nb es o l v e de x p l i c i t l yb e - c a u s eo fa p p r o p r i a t ec o m b i n a t i o n t h i si se v a n s - a b d u u a h sg r o u pe x p l i c i t ( g e ) t h i s w o r ki n d i c a t e dt h a ti t sp o s s i b l et oc o n s t r u c tn e wd i f f e r e n c em e t h o d sw h i c hs a t i s f yt h e a b o v ec o n d i t i o n s b u tw h e nt h e ye x t e n d e dt h em e t h o dt ov a r i a b l ec o e f f i c i e n tp r o b l e m ， t h ep r o v i n go fs t a b i l i t yi sd i f f i c u l t b a s e do nt h i s ，z h a n gb a o l i ne ta 1 p r o p o s e dt h ei d e aw h i c hc o n s t r u c t e dt h es e g - m e n ti m p l i c i ts c h e m eb yu s i n gt h es a u l y e va s y m m e t r i cs c h e m e ，a n ds e tu pav a r i e t y o fe x p l i c i t i m p l i c i ta n dp u r ei m p l i c i ta l t e r n a t i n gp a r a l l e lm e t h o d sb ym a k i n gu s eo f a l t e r n a t et e c h n o l o g y 1 4 - 1 6 t h e s em e t h o d sc a nk e e pt h es t a b i l i t ya n dp a r a l l e l i s m a f - t e rt h a t ，t h e ye x t e n d e di tt ov a r i a b l ec o e f f i c i e n tp r o b l e ma n dp r o v e di t ss t a b i l i t yb y e n e r g ym e t h o d i nt h ec o u r s eo fn u m e r i c a le x p e r i m e n t s ，t h e yf o u n dt h a tt h er e s u l to f s e g m e n to rb l o c kp a r a l l e la l g o r i t h mi sb e t t e rt h a nt h er e s u l to ft h em e t h o dn os p l i t t i n g s oc o n s t r u c t i n gn e wm e t h o db yd i v i d ea n dc o n q u e rs t r a t e g yc a nn o to n l yb eu s e df o r p a r a l l e lc o m p u t i n gb u ta l s oi m p r o v ec a l c u l a t i o na c c u r a c y - a tt h es a m et i m e ，t h e r e a r em a n yr e s e a r c hc o m i n gi n t ob e i n g f o re x a m p l e ，h a nz h e ns t u d i e dak i n do fp u r e e x p l i c i t - i m p l i c i ts e g m e n ta n db l o c ka l t e r n a t i n gm e t h o di nd e t a i l 1 7 ，1 8 f e n g h u ie ta 1 c o n s t r u c t e dt h en e wi t e r a t i o nm e t h o df o re l l i p t i c a le q u a t i o nb ye l i m i n a t i o nb e t w e e n d i f f e r e n c es c h e m eo fd i f f e r e n tn o d e s ，a n dt h em e t h o dh a ds a m ep a r a l l e l i s ma sj a c o b i m e t h o da n dh i g h e rc o n v e r g e n c er a t e 1 9 2 1 z h a n gz h i y u eg a v et h eg r o u pe x p l i c i t s c h e m ef o rp a r a b o l i cp r o b l e mw i t hv a r i a b l ec o e f f i c i e n ta n dp r o v e dt h es t a b i l i t yb ye n - e r g ym e t h o d 2 2 w a n gw e n q i a e ta 1 c o n s t r u c t e da l t e r n a t i n gs e g m e n te x p l i c i t i m p l i c i t s c h e m ef o rd i f f e r e n tp r o b l e m ，p r o v e dt h e i rs t a b i l i t ya n dg a v en u m e r i c a le x p e r i m e n t s 1 3 2 6 r e c e n t l y , s a n j i v ak l e l ep r o p o s e dh i g h - r a t ec o m p a c td i f f e r e n c es c h e m ea n dd o f o u r i e ra n a l y z i n ga b o u te r r o ra n dc o m p a r ei tw i t ht r a d i t i o n a ls c h e m e 2 7 2 8 - 3 2 】m a r k h c a r p e n t e re ta 1 p r o p o s e ds o m eh i g h o r d e rc o m p a c td i f f e r e n c es c h e m ef o rd i f f e r e n t p r o b l e ma n dd i dn u m e r i c a la n a l y s i s t h e nw h a tw i l lb eo b t a i n e db yc o m b i n i n gt h e h i g h o r d e rd i f f e r e n c es c h e m ew i t ha l t e r n a t i n gg r o u p ? t h e r ew e r em a n yr e s e a r c hc o m i n gi n t ob e i n g 4 4 一s 2 i t st h en e w f o c a lt h a tc o m b i n i n gt h eh i g h - r a t es c h e m ea n dt h e s t r a t e g yo fd i v i d i n ga n dc o n q u e r i n g u n d e rp r o f w a n gw e n q i a sc a r e f u l l yi n s t r u c t i n g ，t h ea u t h o rc o n s t r u c t ss o m ep a r a l l e ld i f f e r e n c em e t h o d sf o rs o m ep a r a b o l i cp r o b l e mw h i c hc o n t a i ni t e r a t i o nm e t h o da n d h i g h r a t ea l t e r n a t i n gg r o u ps c h e m e a n dt h es t a b i l i t yo ft h e s em e t h o d sa r ep r o v e da n d s o m en u m e r i c a le x p e r i m e n t si n d i c a t et h e i ra p p l i c a b i l i t y t h ep a p e re x t e n d st h ew o r k 山东大学博士学位论文 o ft h ep r e d e c e s s o r s ，a n dh a sn o n - r e p e a t a b i l i t y t h ep a p e ri sd i v i d e di n t of i v ec h a p t e r s i nc h a p t e r1 ，an e wi t e r a t i o nm e t h o df o r2 dc o n v e c t i o nd i f f u s i o np r o b l e mi sc o n - s t r u c t e rb ym a k i n gu s eo fn u m e r i c a ls t e n c i l 1 9 f i r s t ，i tg i v e st h ed e f i n i t i o no fs t e n c i l f o rp a r a b o l i ce q u a t i o n ，a n do b t a i n st h ef i n a ln u m e r i c a ls t e n c i la f t e rt h r e es t e n c i le l i m - i n a t i o n b a s e do nt h i s ，t h en e wi t e r a t i o ns c h e m ei sc o n s t r u c t e d t h e nt h ec o n v e r g e n c e o fi t e r a t i o ni sp r o v e db ya n a l y z i n gt h ei t e r a t i o ne r r o ra n dt h ec o n v e r g e n c er a t ei sc o m - p a r e dw i t hj a c o b im e t h o d s f i n a l l yt h ep a p e rg i v e sn u m e r i c a le x p e r i m e n tt os h o wi t s a p p l i c a b i l i t y t h ew o r ka b o u t 1i sa c c e p t e db y ( ( i n t e r n a t i n a lj o u r n a lo fc o m p u t e r m a t h e m a t i c s ) ) t h e n e wi d e ao fc h a p t e r1i st h a tn u m e r i c a ls t e n c i li sf i r s ta p p l i e dt o2 dp a r a b o l i c e q u a t i o n ，a n dt h eh i g h - r a t ec o n v e r g e n c ea n dp a r a l l e l i s mi t e r a t i o ni sc o n s t r u c t e d t h e p r o v i n go fs t a b i l i t yi so b t a i n e da n dn u m e r i c a le x p e r i m e n t ss h o w si t sa p p l i c a b i l i t y c h a p t e r2m a i n l yu s e st h ei d e ai n 【2 2 ，2 4 】a n dg i v e st h ea l t e r n a t i n gg r o u pe x p l i c i t f o rc o n v e c t i o nd i f f u s i o ne q u a t i o nw i t hv a r i a b l ec o e f f i c i e n t t h ea l t e r n a t i n gg r o u pe x - p l i c i ts c h e m ei sc o n s t r u c t e db yc o m b i n i n gc r a n k - n i c o l s o nd i f f e r e n ts c h e m ea n dt h ei d e a o fa l t e r n a t i n gg r o u p f i r s ti tg i v e sf o u ra s y m m e t r i cd i f f e r e n c es c h e m eb a s e do nc r a n k - n i c o l s o nd i f f e r e n c es c h e m e ，a n dc o n s t r u c t st h ea l t e r n a t i n gg r o u pe x p l i c i ts c h e m eb y c o m b i n i n gt h e s ea s y m m e t r i cd i f f e r e n c es c h e m e t h e ni t ss t a b i l i t yi sp r o v e db ye n e r g y m e t h o da n dn u m e r i c a le x p e r i m e n ti n d i c a t e si t sv a l i d i t y t h ew o r ka b o u t 2h a sb e e n s u b m i t t e dt o ( ( i n t e r n a t i o n a lj o u r n a lo fc o m p u t e rm a t h e m a t i c s ) ) t h en e wi d e ao fc h a p t e r2i st h a tt h ea l t e r n a t i n gg r o u pe x p l i c i ts c h e m ei sc o n s t r u c t e df o rc o n v e c t i o nd i f f u s i o ne q u a t i o nw i t hv a r i a b l ec o e f f i c i e n tb yc o m b i n i n gt h e i d e ao fa l t e r n a t i n gg r o u pw i t hc r a n k n i c o l s o ns c h e m e t h ef o l l o w i n gc h a p t e r sm a i n l yu s e st h ei d e ai n 【1 4 ，2 3 2 6 ，4 6 4 8 】a n di n t r o d u c e st h e h i g h r a t ed i f f e r e n c es c h e m e 2 7 3 2 b a s e do nt h i s ，s o m eg r o u pm e t h o d sa r ec o n s t r u c t e d f o rc o n v e c t i o nd i f f u s i o ne q u a t i o n t h e s es c h e m ea r ea l la b s o l u t e l ys t a b l ea n dh a v ep a r a l l e l i s m t h er a t eo fl o c a lt r u n c a t i o ne r r o rc a nr e a c ho ( h 4 ) i nc h a p t e r3 ，t h eh i g h r a t ee x p l i c i ta n di m p l i c i td i f f e r e n c es c h e m ea r eg i v e nf i r s t t h e nf o u ra s y m m e t r i cs c h e m e sa r ec o n s t r u c t e db a s e do ni m p l i c i ts c h e m ea n dt h ea l t e r n a t i n gs e g m e n te x p l i c i t i m p l i c i ts c h e m ei sc o n s t r u c t e db ya p p r o p r i a t ec o m b i n a t i o n o fa b o v ed i f f e r e n c es c h e m e s t h ea b s o l u t e l ys t a b i l i t yi sp r o v e db yk e l l o g gl e m m ai n 3 3 3 4 】a n dt h el o c a lt r u n c a t i o ne r r o ri so b t a i n e db yd e r i v a t i o n f i n a l l yn u m e r i c a le x - p e r i m e n t si n d i c a t e st h ea p p l i c a b i l i t ya n dt h er a t eo fl o c a lt r u n c a t i o ne r r o rc a nr e a c h o ( h 4 ) t i l ew o r ka b o u t 3i sa c c e p t e db y ( ( c h i n e s ej o u r n a lo fc o m p u t a t i o n a lp h y 7 s i c s ) ) 山东大学博士学位论文 c h a p t e r4g i v e st h ea l t e r n a t i n gg r o u pe x p l i c i ts c h e m ef o rt h es a m ee q u a t i o n f i r s t ， i tg i v e sh i g h - r a t ec r a n k - n i c o l s o nd i f f e r e n c es c h e m e s b a s e do nt h i s ，e i g h ta s y m m e t r i c d i f f e r e n c es c h e m e sa r ec o n s t r u c t e di no r d e rt oc o n s t r u c tt h ea l t e r n a t i n gg r o u pe x p l i c i t s c h e m e t h e nt h ea b s o l u t es t a b i l i t yi s p r o v e db yt h es a m ew a y , a n dt h et r u n c a t i o n e r r o rc a nr e a c ho ( t h ) b e c a u s es o m ep a r to ft r u n c a t i o ne r r o rc a nb ec a n c e l l e db yu s i n g d i f f e r e n ts c h e m ea l t e r n a t e l yb e t w e e nt w ot i m el a y e r s f i n a l l yn u m e r i c a le x p e r i m e n t s h o w st h em e t h o di sa p p l i c a b l e t h ew o r ka b o u t 4i sp u b l i s h e di n ( ( j o u r n a lo fs h a n - d o n gu n i v e r s i t y ) ) ( n a t u r a ls c i e n c e ) b a s e do nc h a p t e r4 ，c h a p t e r5i n t r o d u c e st h ei d e ao fa l t e r n a t i n gs e g m e n tc r a n k - n i c o l s o ns c h e m e i tc o m b i n e st h ee i g h ta s y m m e t r i cs c h e m ew i t hh i g h - r a t ec r a n k - n i c o l s o ns c h e m ea n dc o n s t r u c t st h eh i g h - r a t ea l t e r n a t i n gs e g m e n tc r a n k - n i c o l s o n s c h e m e i th a sa b s o l u t es t a b i l i t ya n dp a r a l l e l i s m i na d d i t i o n ，t h et r u n c a t i o ne r r o ro f t h ep o i n t si sm u c hl e s sb e c a u s em o r ep a r t sa r ec a n c e l e db yo n l yu s i n gc r a n k - n