项目名称不可压缩流动问题高效数值方法研究及应用-宝鸡文理学院.doc

项目名称：不可压缩流动问题高效数值方法研究及应用 完成单位：西安交通大学、 宝鸡文理学院 完成人：何银年、李剑 项目简介： 本项目属于计算数学及计算流体动力学的交叉学科领域。 Navier-Stokes（N-S）方程组是描述流体运动的典型非线性方程组，在工程应用和非线性科学研究方面非常重要，并对人们认识和控制湍流及复杂流体动力学行为具有重大意义。然而人们对非线性现象本质的认识有限，这导致数值模拟成为一种十分重要的研究手段。但直接数值模拟N-S方程组有四大难题，即如何处理不可压缩约束条件 div u=0、强非线性、长时间区间积分和小粘性参数的问题。本项目针对在求解不可压缩N-S方程组中出现的上述四大难题，用有限元和有限体积元方法设计出高效稳定的数值格式，并进行了全面和系统的分析研究，建立了一套有关不可压缩N-S方程组数值模拟求解的完整理论分析和数值计算方法。主要发现点如下： 提出了求解N-S方程组的低等阶元局部高斯积分稳定化方法，证明了数值解的适定性和最优收敛阶，解决了处理不可压缩约束条件div u=0的难题。并首次给出对应的有限体积方法系统的最优阶理论分析。欧洲科学院院士Alfio Quarteroni教授发表文章指出：“Li和Chen提出的特殊的Petrov-Galerkin方法，推进了有限体积元方法的发展”。 提出了求解非定常N-S方程组的全离散时空多水平算法，将非定常N-S方程组在时空细网格上求解归结为时空粗网格上的半隐格式和时空细网格上的隐式/显式格式的求解，证明了数值解的几乎无条件稳定性和收敛性，降低了计算工作量，部分解决了如何处理强非线性和长时间区间积分的难题。法国科学院院士Temam教授评价指出：“He等人提出的多水平方法由于降低计算量的优点，被成功应用于求解非定常N-S方程组”。 发现并证明了在初值光滑条件下求解非定常N-S方程组的线性化隐式/显式全离散有限元算法是几乎无条件收敛的，突破了传统意义上时间步长依赖网格尺度的苛刻约束条件限制，解决了如何处理强非线性和长时间区间积分的难题。美国三院院士Hughes教授、美国工程院院士Oden教授、法国科学院和欧洲科学院Glowinski教授及巴西科学院院士袁锦韵教授指出：“He等人提出的线性化隐式/显式大时间步长格式属于当前非线性科学研究的热点领域，克服了显式格式对时间步长苛刻的约束条件，并首先给出了线性化隐式/显式全离散有限元算法几乎无条件稳定和收敛的证明”。 提出了求解定常N-S方程组的粘性相关迭代法和时空隐式/显式迭代算法，证明了数值解的稳定性和最优收敛阶，解决了求解定常N-S方程组时如何处理小粘性参数的难题。加拿大阿尔伯塔大学的资深教授Yaushu Wong评价指出：“该方法解决了如何处理小粘性参数的难题”。 本项目在国际知名期刊发表SCI文章63篇， 其中大多数发表在计算数学和力学类国际顶尖期刊SIAM J. Numer. Anal.， Numerische Mathematik，Math. Comp.及Comput. Methods Appl. Mech. Engrg.和J. Comput. Phys.上，其中20篇主要论著被SCI他引601次，8篇代表作被SCI引用421次，其中3篇主要论著[1,2,15]为ESI高被引论文。项目第一完成人在ESI全球数学领域的高被引作者排名居前18%（281/1519），H指数为21，并于2014-2015两年连续入选爱思唯尔中国高被引学者榜单。 主要论文专著目录 A stabilized finite element method based on two local Gauss integrations for the Stokes equations，作者: Li, Jian; He, Yinnian JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (IF：1.266) 卷: 214 期: 1 页: 58-65 出版年: APR 15 2008， A stabilized finite element method based on local polynomial pressure projection for the stationary Navier-Stokes equations ，作者: He, Yinnian; Li, Jian APPLIED NUMERICAL MATHEMATICS (IF:1.221) 卷: 58 期: 10 页: 1503-1514 出版年: OCT 2008 3. On large time-stepping methods for the Cahn-Hilli