（计算数学专业论文）海水入侵问题的间断有限元方法.pdf

j1 原创性声明 本人郑重声明：所呈交的学位论文，是本人在导师的指导下， 独立进行研究所取得的成果。除文中已经注明引用的内容外，本 论文不包含任何其他个人或集体己经发表或撰写过的科研成果。 -对本文的研究作出重要贡献的个人和集体，均已在文中以明确方 _ 式标明。本声明的法律责任由本人承担。 论文作者签名：i 靡厘缝 关于学位论文使用授权的声明 本人完全了解山东大学有关保留、使用学位论文的规定，同意学校 保留或向国家有关部门或机构送交论文的复印件和电子版，允许论文 被查阅和借阅；本人授权山东大学可以将本学位论文的全部或部分内 容编入有关数据库进行检索，可以采用影印、缩印或其他复制手段保 存论文和汇编本学位论文。 ( 保密论文在解密后应遵守此规定) 论文作者签名：似导师签名：孵日 期：翟丛业 山东大学博士学位论文 中文摘要 目录 英文摘要 第一章预备知识 1 1 海水入侵问题的数学模型 1 2 定义和记号 1 2 1 空间和范数 1 2 2 空间剖分 1 2 3 间断s o b o l e v 空间 1 3 定理和不等式。 第二章间断有限元方法和混合有限元方法相结合求解海水人侵问题 2 1 引言 2 2 混合元和问断有限元相结合的组合方法 2 2 1 水头方程 2 2 2 浓度方程 2 2 3 组合半离散格式 2 3 先验误差估计 2 3 1 咒1 范数下的误差估计 2 3 2c 2 范数下的误差估计 2 4 时间离散格式。 2 5 结论 第三章迎风间断有限元方法和混合元方法相结合求解海水人侵问题 3 1 引言 3 2 混合元与迎风间断有限元相结合的组合方法 3 2 1 水头方程 3 2 2 浓度方程 3 2 3 组合全离散格式 3 3 辅助结果 v x l l 3 3 5 6 7 0 0 o 0 l 2 3 4 2 4 9 0 0 1 l l 3 4 - 一 1 1 1 1 1 1 1 1 2 2 2 翻3 3 3 3 3 3 i 山东大学博士学位论文 3 4 误差分析 3 7 3 5 一个半隐格式4 6 3 6 数值算例4 9 3 7 结论 5 3 第四章抛物方程的间断有限元区域分解方法 5 4 4 1 引言 5 4 4 2 区域分解算法 5 5 4 3 误差分析 5 7 4 4 数值实验6 4 4 5 结论 6 9 参考文献 7 0 致谢 7 9 攻读博士学位期间完成论文情况 8 0 作者简介 8 1 山东大学博士学位论文 插图目录 4 1 稳定性实验的初值。6 5 4 2 当7 = 3 5 时，i l i u l i i 随时间的变化情况 6 5 4 3 当7 = 0 4 9 5 时，川u 川随时间的变化情况 6 6 4 4 当2 = 0 2 5 和l 时，i l i u i f i 随时间的变化情况 6 6 4 5 当a t h 2 = 1 5 时，i i i u i l l 随时间的变化情况 6 7 4 6 不匹配网格6 8 山东大学博士学位论文 表格目录 1 1 数学模型中的记号 2 3 1 最低阶混合元和问断分片线性n i p g 5 0 3 21 阶混合元和间断分片二次n i p g5 1 3 3 最低阶混合元和间断分片线性s i p g 5 1 3 41 阶混合元和间断分片二次s i p g 。5 2 3 5 基于n i p g 的半离散格式 5 2 3 6 基于s i p g 的半离散格式 5 3 4 1 收敛性结果6 8 山东大学博士学位论文 海水人侵问题的间断有限元方法 廉西猛 ( 山东大学数学学院，济南，2 5 0 1 0 0 ) 指导老师：芮洪兴教授 摘要 海水入侵是在许多沿海国家和地区都普遍存在的问题，它是指海水侵入地下 淡水蓄水层的复杂运动，一般足由自然水环境的变化和社会经济发展所造成的 由于淡水的密度小于海水的密度，所以在地下蓄水层中，淡水总足存在于海水的 上方，并且在咸淡水界面上达到一种动态的压力甲衡沿海国家和地区为了满足 人类生活和经济发展的需要，往往会过度地开采地下水( 淡水) ，这就会引发海水 入侵问题当淡水的开采量超过了它的补给量，海岸带附近的淡水水头( 水位) 就会随之下降，进而导致淡水在咸淡水界面处的产生压力减小，压力平衡遭到破 坏在压力作用下，海水会侵染淡水层，以达到新的甲衡当海水侵染到淡水开 采井时，开采出的水就不再适合饮用和生产，这将给沿海地区的人类生活和社会 经济发展带来严重的危害目前，全世界有几十个国家和地区的几百个地方都发 现了海水入侵问题，因此，这一问题受到了国际社会的广泛关注，许多国家都在 积极开展海水入侵问题的研究和治理工作 近些年来，对海水入侵问题的数学模型及数值解法的研究已经取得了一些成 果基于咸淡水界面为两种流体混合的过渡带的假设，h e n r y 5 1 1 得到了稳定流 在匀质介质、简化边界条件情况下的解析解1 9 7 0 年，p i n d e r 6 7 】最早提出了海 水入侵过渡带模型，并提出了h e n r y 模型的有限元数值解法s e g o l 等人f 7 7 ，7 8 】 发展了地下水位和浓度相互依赖的剖面二维有限元模型h u y a k o r n 等【5 4 】提出 了依赖浓度的地下水流方程和溶质运移方程，研究了三维情形的有限元模型但 是，上述海水入侵的研究工作都是在一些具体的假设条件得到的，并不能真实的 反映海水入侵过程9 0 年代初，薛禹群等【8 5 ，8 6 】研制我国第一个三维有限元海 水入侵模型他们考虑了速度随盐分变化的情形，但是没有给出对防止海水入侵 的各类工程的数值预测；参见【4 7 】及其引文自上世纪末，袁益让，梁栋和芮洪 兴教授深入研究了海水入侵数学模型及其强对流占优和非线性的特点，提出了处 理海水入侵模型的特征有限元【8 8 】、迎风分数步有限差分 8 9 ，9 0 ，9 1 ，9 2 】等方法， 山东大学博士学位论文 对海水入侵防治工程实施后咸淡水变化运移的真实过程进行了定量模拟类似于 文献【3 9 ，3 8 ，4 4 1 中处理油藏问题的方法，这些方法克服了普通差分或有限元方 法容易产生数值振荡的缺点，但也存在缺点：迎风方法引入的额外的扩散会影响 结果的精度，而特征方法不便于处理非线性问题【4 5 】另外，对于出现间断解的 情形，上述方法不能很好的解决本文在上述研究成果的基础上，充分考虑了海 水入侵模型中相互耦合的两个方程的特性，选用混合元和间断有限元相结合的方 法来处理这一问题，并相应地提出半离散和全离散格式 1 9 7 3 年，r e e d 和h i l l 6 9 在处理中子输运方程一一不依赖时间的线性双曲 型方程时，提出了第一个间断有限元方法自那时起，求解双曲型方程的间断伽 略金方法得到了广泛的关注和积极的发展闻断的思想与其他数值技巧相结合， 发展出了许多不同的数值方法，如j o h n s o n 5 5 1 等人将迎风的思想融入到问断伽 略金方法中，用来处理双曲型方程；c o c k b u r n 等人将高精度有限差分和有限体 积方法的技术( 如迎风、数值流通量、斜率限制器等，参见【5 7 】) 与伽略金方法 相结合，提出了r u n g e - k u t t a 间断伽略金方法【2 4 ，2 3 ，2 2 ，2 0 ，2 5 】，这也是双曲 型方程间断伽略金方法的主要进展【2 1 】为了能够处理椭圆型和对流扩散型方 程，c o c k b u r n 等人又提出了局部间断伽略金方法【2 6 ，1 5 】同样在7 0 年代，求 解椭圆型和抛物型方程的间断伽略金有限元方法( 也被称为内罚方法) 也被独立 地提出来，并得到了深入的研究和发展 3 7 ，7 ，8 4 ，2 】近几年来，许多新型的内 罚方法相继被提出和研究，包括对称内罚伽略金方法【2 】，b a u m a n n - o d c n 方法 【9 ，6 4 1 ，非对称内罚伽略金方法【7 4 ，7 2 ，7 5 】，以及不完全内罚伽略金方法【3 4 】等； 更多细节可参见综述性文献 3 ，1 4 ，4 ，7 0 由于与连续的伽略金有限元方法相比，间断伽略金方法具有局部守恒，易于 实现却自适应，可以更好的逼近闻断解等优点【4 0 】，所以它们被广泛的应用于许 多实际问题。尤其以在计算流体力学领域的应用最为广泛，涉及不可压缩8 0 ，8 1 】 和可压缩【1 6 ，3 0 ，3 1 】混溶驱动问题，n a v i e r - s t o k c s 问题【8 ，4 8 ，4 9 ，5 0 ，6 5 1 ，无粘 性可压流体【3 5 】，两相流【4 3 ，4 1 ，4 2 】，气体动力学【17 】，半导体等器件的数值模拟 【1 7 ，6 0 】，多孔介质流体【8 2 ，7 3 】和多孔弹性介质流体 6 6 】等另外，它们还被应用 于弹性力学【6 1 ，7 l 】等其他物理问题中注意到上述实际问题的模型大多足对流 扩散方程，所以处理此类问题的数值格式往往足通过将两类问断伽略金方法相结 合得到的，即采用内罚方法处理扩散项，而用双曲型间断伽略金方法处理对流项 不仅如此，在处理对流扩散问题，特别是对流占优的情形时，往往在间断伽略金 方法中引入一些有限差分和有限元方法的技巧，如迎风【5 5 ，7 9 ，5 2 ，8 1 ，8 2 ，3 1 ，7 0 1 ， 山东大学博士学位论文 斜率限制器【7 0 】和数值流通量【3 5 ，3 6 】等，以取得更好的逼近效果 混合元方法自出现以来，就得到了广泛和深入地研究【1 1 ，6 8 ，6 2 ，6 3 】，主要 用于求解椭圆型和抛物型偏微分方程【5 6 ，1 3 ，1 2 】混合元方法首先通过引入一个 新的变量( 一般足未知函数导数的某种线性组合) ，将最初的方程化成混合形式， 然后对混合变分形式使用有限元方法进行求解这样可以将未知函数和新引入的 变量同时求出，并且对这两个函数的逼近都可以达到最优阶基于这一优点，早 在上世纪八十年代，混合元方法就被用于与其他方法相结合，共同求解耦合的偏 微分方程系统，如多孔介质混溶驱动模型【3 8 ，4 4 】近几年来，随着间断伽略金方 法的发展，这种新型的方法也经常与混合元联系在一起使用，来处理一些复杂的 耦合系统，例如不可压缩【8 0 ，8 1 】和可压缩【1 6 ，3 0 ，3 1 】混溶驱动问题，多孔弹性 流体问题【6 6 】等由于海水入侵的数学模型是一个由一个抛物方程和一个对流 扩散方程组成的耦合系统，本论文也采用类似于上述的想法来处理这一问题，即 利用混合元方法求解抛物方程，用间断伽略金方法求解对流扩散方程 在芮洪兴教授的精心指导下，作者完成博士论文的写作在本论文中，作者 充分研究了袁益让等人在文献【8 8 ，8 9 ，9 0 ，9 1 ，9 2 】中考虑的海水入侵数学模型， 提出了几种处理该模型问题的数值格式本文所采用的数值方法是混合元方法 和间断伽略金有限元方法，又由于此海水入侵的数学模型是耦合的偏微分方程系 统，所以本文中还用到了迭代方法应用这些方法，作者提出了连续时间格式和 隐式离散时问格式对这些格式，我们在理论上分析了格式的收敛性，得到了先 验的误差估计为了提高计算的速度和效率，作者改进了全隐离散时问格式，提 出了一个半隐的离散格式，并在最后给出了并行求解此格式的区域分解方法本 文共分为四个章节，下面简要介绍一下各章的主要内容 第一章，我们做一些准备工作首先，给出了在文献郾，8 9 ，9 0 ，9 1 ，9 2 】中提 出并研究的海水入侵的数学模型，该模型足由一个抛物方程和一个对流扩散方程 组成的耦合系统对于此系统，我们给出了两种不同的边值条件，并在后面几章 分别给出相应的处理格式接下来，给出了一些定义和记号，包括s o b o l e v 空间 和分片s o b o l e v 空间及它们的范数、空间剖分及其属性、标量和矢量跳跃和平均 的定义等最后，以引理的形式给出了在误差估计中必需的定理和不等式，这些 引理在许多的书籍和文献中都曾被提到或证明 第二章，对于前一章给出的数学模型和第一个边值条件，我们提出一个连续 时间格式混合元方法被用来处理数学模型中的抛物方程，而间断伽略金有限元 方法( 对称内罚方法) 被用来处理系统的对流扩散方程，其中对流项只是单纯做 山东大学博士学位论文 了内积处理在给出了上述方程的变分形式的基础上，我们提出了迭代的连续时 间格式对于该格式，得到了c 2 范数和咒1 范数下的最优误差估计这一部分 内容主要出自文章【5 8 】，该文章已被j o u r n a lo fs y s t e m ss c i e n c ea n dc o m p l e x i t y ( s c i ) 接收。 接着，针对上述连续时间格式，提出了一些时问离散方法，涉及显式和隐式 方法以及一些高阶方法对于这些方法，我们简要地研究了它们的稳定性和收敛 性，并讨论了其优缺点 第三章，我们主要考虑附加了第二个边值条件的海水入侵模型模型中的抛 物方程仍然采用混合有限元方法来求解。而对于对流扩散方程，为了更好的逼近 方程中的对流项，我们使用迎风间断伽略金有限元方法来处理另外，通过引入 一个参数a 和口，我们在同一个表达式中考虑了两种不同的内罚方法：对称内 罚和不对称内罚方法以及两种时间离散方法：向后e u l e r 方法和c r a n k - n i c o l s o n 方法在此基础上，我们得到了全隐式的迭代格式理论上的坳先验误差估计 和设计的数值实验结果都证明此格式的有效性和可靠性。这些结果主要出自文章 【5 9 】 众所周知，非线性问题的全隐格式的求解效率很低，而显格式是条件稳定的， 时间步长会受到限制为了改变这种两难处境，我们提出了一种半隐格式，即对 于线性的扩散项采用隐式处理，而非线性的对流项采用显式处理理论分析和数 值实验都表明此格式的计算精度与全隐格式相差无几 第四章，我们讨论了一种求解抛物方程的基于间断有限元方法的区域分解算 法由于上一章中的半隐格式可以看成足对抛物方程的全隐格式，所以这一章的 并行方法可以有效的提高半隐格式的计算速度我们以一个标准的抛物方程为 例，提出了区域分解算法，并且考虑它的稳定性和收敛性，最后给出了数值算例 以验证我们的结论这一章的结果出自文章【7 6 】 关键词：海水人侵问题；间断有限元方法；混合有限元方法；先验误差估计； 区域分解 山东大学博士学位论文 d i s c o n t i n u o u sf i n i t ee l e m e n t 1 。1 h o d sforinuousi ne m ei v l e t l 2 0si o r s e a w a t e ri n t r u s i o np r o b l e m x i m e n gl i a n ( s c h o o lo fm a t h e m a t i c s ，s h a n d o n gu n i v e r s i t y , j i n a n2 5 0 1 0 0 ，c h i n a ) s u p e r v i s o r ：p r o f h o n g x i n gr u i a b s t r a c t s e a w a t e ri n t r u s i o no c c u r sw i d e l yi nm a n yc o a s t a lc o u n t r i e sa n dr e g i o n s ，i t i st h ec o m p l i c a t e dm o v e m e n to fs a l i n ew a t e ri n t oaf r e s h w a t e ra q u i f e r ，w h i c hi s g e n e r a l l yc a u s e db yc h a n g e so fn a t u r a lw a t e re n v i r o n m e n ta n d t h es o c i a la n d e c o n o m i c a ld e v e l o p m e n t i ti sw e l lk n o w nt h a tt h ed e n s i t yo ft h ef r e s h w a t e ri s l o w e rt h a nt h a to ft h es e a w a t e r ，s ot h ef r e s hw a t e rc o l u m nw o u l df l o wo v e rt h e c o l u m no ft h es e a w a t e ri nt h ea q u i f e r ，a n dt h e yw i l la r r i v ed y n a m i ce q u i l i b r i u m o fp r e s s u r e so nt h ei n t e r f a c e i no r d e rt os a t i s f yt h ed e m a n do ft h ew a t e rf o r p e o p l e sl i v e sa n de c o n o m i c a lu s e s ，t h ec o a s t a lc o u n t i e sa n dr e g i o n so f t e nm i n e t h eg r o u n d w a t e re x c e s s i v e l y , w h i c hw i l li n d u c et h cp r o b l c mo fs e a w a t e ri n t r u s i o n m e nf r e s hw a t e ri sw i t h d r a w na taf a s t e rr a t et h a ni tc a nb er e p l e n i s h e d ，t h e w a t e rt a b l ei sd r a w nd o w na sar c s u l t t h ed e c l i n eo ft h ew a t e rt a b l er e d u c e st h e p r e s s u r eo nt h ei n t e r f a c e ，w h i c hl e a d st ot h ed e s t r u c t i o no fb a l a n c e t or e a c ht h e n e wb a l a n c eo fp r e s s u r e s 。s e a w a t e ri sp u l l e di n t ot h ef r e s h w a t e ra q u i f e r 册l e n t h ep u m p i n gw e l li sc o n t a m i n a t e do w i n gt ot h es e a w a t e ri n v a s i o n ，t h ew a t e rc a n n o tb ep o t a b l ea n du s e df o rp r o d u c t i o n ，a n dt h i sb r i n g sg r e a td a m a g e st op e o p l e s l i v e sa n ds o c i a la n de c o n o m i c a ld e v e l o p m e n to ft h ec o a s t a lc o u n t r i e sa n dr e g i 0 1 1 8 s of a rs e a w a t e ri n t r u s i o np r o b l e mh a sb e e nf o u n di nh u n d r e d so fa r e a si ns c o r e s o fc o u n t r i e sa n dr e g i o i l s ，t h u si ti sw i d e l yc o n c e r n e db yi n t e r n a t i o n a lc o m m u n i t y , a n di su n d e ra c t i v er e s e a r c ha n dg o v e r n m e n t i nr e c e n td e c a d e s ，t h e r ei ss o m el i t e r a t u r eo nt h em a t h e m a t i c a lm o d e la n d i t sn u m e r i c a la p p r o x i m a t i o no fs e a w a t e ri n t r u s i o np r o b l e m b a s e do nt h e 瓣 s u m p t i o nt h a tt h ei n t e r f a c eb e t w e e nt h et w of l u i di sa t r a n s i t i o nz o n e ，h e n r y 【5 1 】 o b t a i n e dt h ea n a l y t i cs o l u t i o nf o rt h es t e a d yf l o wi nt h ec a s eo fh o m o g e n e o u s 山东大学博士学位论文 m e d i u ma n ds i m p l i f i e db o u n d a r yc o n d i t i o n i n1 9 7 0 ，p i n d e r 【6 7 f i r s td e v e l o p e d t h et r a n s i t i o nz o n em o d e lo fs e a w a t e ri n t r u s i o n ，a n dt h e ya l s op r o p o s e daf i n i t e e l e m e n tm e t h o d sf o rh e n r y sm o d e l s e g o le ta 1 【7 7 ，7 8 】d e v e l o p e da c u t p l a n e t w od i m e n s i o n a lf i n i t e - e l e m e n tm o d e lw i t ht h ew a t e rt a b l ea n dc o n c e n t r a t i o nc o i l - p l i n gw i t he a c ho t h e r h u y a k o me ta 1 【5 4 】p r o p o s e dt h eg r o u n d w a t e re q u a t i o n w h i c hd e p e n do nc o n c e n t r a t i o na n dt h es o l u t et r a n s p o r te q u a t i o n ，a n dc o n s i d - e r e df i n i t e - e l e m e n tm o d e lo ft h r e e - d i m e n s i o n a lc a s e h o w e v e r ，t h er e s e a r c hr e s u l t s a b o v ew h i c ha r eo b t a i n e du n d e rs o m ep r e c i s ea s s u m p t i o n s ，c a n n o tr e a l l yr e f e c t t h es e a w a t c ri n t r u s i o np r o c e s s e a r l yi n1 9 9 0 s ，x u ee ta 1 【8 5 ，8 6 】d e s i g n e dt h e f i r s tt h r e e - d i m e n s i o n a lf i n i t e - e l e m e n ts e a w a t e ri n t r u s i o nm o d e l t h e yc o n s i d e r e d t h ec a s et h a tt h ef l u i dv e l o c i t yv a r i e sw i t ht h ec o n c e n t r a t i o n ，b u tt h e yd i dn o t g i v et h en u m e r i c a lp r e d i c t i o nf o rt h ep r o t e c t i o np r o j c c t sa v o i d i n gs e a w a t e ri n t r u - s i o n ；f o rd e t a i ls e e 4 7 1a n dc i t a t i o n st h e r e i n s i n c et h ee n do ft h el a s tc e n t u r y , y u a ne ta 1 d e e p l ys t u d i e dt h em a t h e m a t i c a lm o d e lo ft h es e a w a t e ri n t r u s i o n ， w h i c hh a st h ec h a r a c t e r i s t i co fc o n v e c t i o n - d o m i n a t e dd i f f u s i o na n dn o n l i n e a r i t y , a n dp r o p o s e dc h a r a c t e r i s t i c sf i n i t ee l e m e n tm e t h o d s 8 8 】a n dm o d i f i e dm e t h o d o fu p w i n dw i t hf i n i t ed i f f e r e n c ef r a c t i o n a ls t e p sp r o c e d u r e 8 9 ，9 0 ，9 1 ，9 2 】f o r s e a w a t e ri n t r u s i o nt os i m u l a t et h et r a n s p o r to ft h em i m n gf l u i da f t e rt h ep r o t e c - t i o np r o j e c t si sc a r r i e do u t s i m i l a rt ot h em e t h o d sd e v e l o p e df o ro i lr c s e r v i o r p r o b l e m si n 【8 9 ，9 0 ，9 1 ，9 2 】，t h e s em e t h o d sa v o i dn u m e r i c a lo s c i l l a t i o nw h i c hi s t h ed i s a d v a n t a g eo ff i n i t ed i f f e r e n c eo rf i n i t ee l e m e n tm e t h o d s h o w e v e r ，t h e y a l s oh a v ed r a w b a c k s ：e x c e s s i v ed i f f u s i o ni n t r o d u c e db yu p w i n dm e t h o d sa f f e c t s t h ea c c u r a c y ；t h ec h a r a c t e r i s t i cm e t h o d sa r ei n c o n v e n i e n tf o rn o n l i n e a rp r o b l e m s 【4 5 】f u r t h e r m o r e ，t h e yc a n n o td ow e l lw i t ht h ep r o b l e m sw i t hd i s c o n t i n u o u s s o l u t i o n s i no u rp a p e r ，a f t e ro b s e r v i n gt h ec h a r a c t e r i s t i co ft h ec o u p l e de q u a - t i o n s ，w ed e v e l o p e ds c h e m e sc o m b i n i n gd i s c o n t i n u o u sg a l e r k i nw i t hm i x e df i n i t e e l e m e n tm e t h o d sf o rs e a w a t e ri n t r u s i o nt oo v e r c o m et h ed r a w b a c k sa b o v e t h ef i r s td i s c o n t i n u o u sg a l e r k i nm e t h o dw a si n t r o d u c e di n1 9 7 3b yr e e d a n dh i l l 【6 9 】f o rt h en e u t r o nt r a n s p o r te q u a t i o nw h i c hi sat i m e - i n d e p e n d e n t l i n e a rh y p e r b o l i ce q u a t i o n ，a n ds i n c et h e nd i s c o n t i n u o u sg a l e r k i nm e t h o d sf o r h y p e r b o l i ce q u a t i o nh a v eb e e nw i d e l yc o n c e r n e da n da c t i v e l yd e v e l o p e d m a n y d i f f e r e n tn u m e r i c a lm e t h o d sw e r ep u tf o r w a r db yi n c o r p o r a t i n gd i s c o n t i n u o u s x 山东大学博士学位论文 g a l e r k i nm e t h o d sw i t ho t h e rm e t h o d s f o re x a m p l e ，j o h n s o na n dp i t k k r a n t a 5 5 】d e v e l o p e dm e t h o d sf o rh y p e r b o l i ce q u a t i o n ，w h i c hi so b t a i n e db yc o m b i n i n g d i s c o n t i n u o u sg a l e r k i nm e t h o d sw i t hu p w i n dm e t h o d ；c o c k b u r ne la 1 i n c o r p o - r a t e ds oc a l l e d h i g hr e s o l u t i o n t e c h n i q u e s ( s u c ha su p w i n d ，n u m e r i c a lf l u xa n d s l o p el i m i t c r s ；s e e 【5 7 1t ok n o wt h ei d e a sf o rd e t a i l ) f o rf i n i t ed i f f e r e n c eo rf i _ n i t ev o l u m em e t h o d si n t ot h et h ef r a m e w o r ko fd i s c o n t i n u o u sg a l e r k i nm e t h o d s a n dp r o p o s e dr u n g e - k u t t ad i s c o n t i n u o u sg a l e r k i n ( r k d g ) m e t h o d si nt h e i r s e r i a lp a p e r sf 2 4 ，2 3 ，2 2 ，2 0 ，2 5 】，w h i c ha r et h em a j o rd e v e l o p m e n to fd i s c o n - t i n u o u sg a l e r k i nm e t h o d sf o rh y p e r b o l i ce q u a t i o n s 【2 1 】t oh a n d l ee l l i p t i ca n d c o n v e c t i o n - d i f f u s i o ne q u a t i o n s ，b a s e do nt h ei d e a so fr u n g e - k u t t ad i s c o n t i n u o u s g a l e r k i nm e t h o d s ，c o c k b u r na n ds h u 【2 6 】a n dc a s t i l l oa n dc o c k b u mf 1 5 1 p r e - s e n t e dl o c a ld i s c o n t i n u o u sg a l c r k i n ( l d g ) m e t h o d s a l s oi n1 9 7 0 s ，b u ti n d c p e n - d e n t l y , d i s c o n t i n u o u sg a l e r k i nf i n i t ee l e m e n tm e t h o d sf o re l l i p t i ca n dp a r a b o l i c e q u a t i o n s ( a l s oc a l l e di n t e r i o rp e n a l t y ( i p ) m e t h o d s ) w e r ep r o p o s e da n ds t u d i e d 【3 7 ，7 ，8 4 ，2 】i nr e c e n ty e a r s ，m a n yn e wt y p e so fi n t e r i o rp e n a l t ym e t h o d sa r ed e - v e l o p e d ，i n c l u d i n gs y m m e t r i ci n t e r i o rp e n a l t yg a l e r k i n ( s i p g ) m e t h o d s 【2 】，t h e m e t h o do fb a u m a n na n do d e n 【9 ，6 4 】，n o n s y m m e t r i ci n t e r i o rp e n a l t yg a l e r k i n ( n i p g ) m e t h o d s 【7 4 ，7 2 ，7 5 1 ，i n c o m p l e t ei n t e r i o rp e n a l t yg a l e r k i n ( i i p g ) m e t h - o d s ( 3 4 】，a n ds oo n ；s e es u m m a r yl i t e r a t u r e 3 ，1 4 ，4 ，7 0 lf o rm o r ei n f o r m a t i o n b e c a u s ed i s c o n t i n u o u sg a l e r k i nm e t h o d sc a nc o n s e r v em a s s l o c a l l y , b cs u i t a b l ef o rh pa d a p t i v i t ya n da p p r o x i m a t es o l u t i o nw i t hd i s c o n t i n u i t i e sb e t t e rc o m - p a r i n gw i t hc o n t i n u o u sg a l e r k i nm e t h o d s ，t h e ya r eo f t e ne m p l o y e dt od e a lw i t h p r a c t i c a lp r o b l e m s c o m p u t a t i o n a lf l u i dd y n a m i ci st h em o s tw i d e l yu s e df i e l d ， f o re x a m p l e s ，i n c o m p r e s s i b l e 8 0 ，8 1 】a n dc o m p r e s s i b l e 【1 6 ，3 0 ，3 1 】m i s c i b l ed i s - p l a c e m e n tp r o b l e m s ，g a sd y n a m i cp r o b l e m s 17 ，s e m i c o n d u c t o rd e v i c es i m u l a t i o n 【17 ，6 0 ，n a v i e r - s t o k e sp r o b l e m s 【8 ，4 8 ，4 9 ，5 0 ，6 5 ，t w o - p h a s ef l o w 【4 3 ，4 1 ，4 2 ， i n v i s c i dc o m p r e s s i b l ef l u i d 【3 5 】，f l u i di np o r o u sm e d i a 7 3 ，8 2 】a n dp o r o e l a s t i c i t y m e d i a 【6 6 ，a n ds oo n a d d i t i o n a l l y , t h e ya r ea l s oa p p l i e dt oo t h e rp h y s i c a lp r o b - l e m s ，f o re x a m p l ee l a s t i c i t y 【6 1 ，7 1 o b s e r v i n gt h a tt h em o d e l sm e n t i o n e da b