有限差分法中几何守恒律的机理及算法.docVIP

有限差分法中几何守恒律的机理及算法 刘君;韩芳;夏冰 【摘 要】采用有限差分法求解复杂外形物体绕流场时经常进行坐标变换,由此会引入坐标变换系数等几何参数,采用不同的差分格式离散坐标变换系数得到的结果不同,导致在计算过程中可能出现均匀流场不能保持均匀的现象,消除这种误差需要研究几何守恒律.本文对坐标变换过程进行理论分析,发现坐标变换过程中采用的数学恒等式在离散条件下不再成立,这是引起物理量不守恒的本质机理,认为增加坐标变换系数恒等式作为源项的方程形式才是曲线贴体坐标系下的离散等价方程,提出只要源项和对流项的离散格式相同就能满足几何守恒律的构造准则.按照上述理论准则建立了基于离散等价方程的几何守恒律算法,通过AUSM、HLLC、Roe、VanLeer四种分裂格式的算例,表明这种新的几何守恒律算法适用于通量差分裂格式(Flux-Difference Splitting,FDS)和矢通量分裂格式(Flux-Vector Splitting,FVS),且均能消除由坐标变换(包括网格运动)引起的误差,保持流场的均匀特性. 【期刊名称】《空气动力学学报》 【年(卷),期】2018(036)006 【总页数】10页(P917-926) 【关键词】几何守恒律;有限差分法;贴体坐标系;离散近似方程;离散等价方程;源项耦合算法 【作 者】刘君;韩芳;夏冰 【作者单位】大连理工大学航空航天学院,辽宁大连 116024;大连理工大学航空航天学院,辽宁大连 116024;大连理工大学航空航天学院,辽宁大连 116024 【正文语种】中 文 【中图分类】V211.3 0 引 言美国在《2030年CFD远景规划》[1]中认为，自适应网格技术和误差估计是复杂流动和复杂几何外形数值模拟的重大瓶颈技术。网格是计算流体力学(CFD)解决复杂外形工程应用的关键，是该学科的重要内容。在网格理论研究方面，除了传统的生成算法外，近些年几何守恒律(Geometric Conservation Law, GCL)问题也受到关注。它直接影响计算结果精度，是解决“误差估计”的基础。1 几何守恒律的研究现状几何守恒律的概念最早是Thomas和Lombard于1978年在AIAA会议上提出[2]，随后在AIAA期刊正式发表[3]。但对几何守恒律问题的研究，最早可见于1961年Trulio和Trigger的工作报告中对一维坐标下“体积守恒(Conservation of Volume)”问题的讨论[4]。对于笛卡尔坐标系(t,x,y,z)的三维守恒型Euler方程：(1)应用于复杂外形时，需要进行坐标变换，在计算空间中(τ,ξ,η,ζ)形式如下：(2)其中：对于流动参数为常数的均匀流场，以上方程变成：(3)这是Thomas针对有限差分法给出的微分形式的GCL方程，但未对此进行深入讨论，因为他们主要关注有限体积法在网格变形时涉及到的GCL问题。对于有限体积法，基于ALE(Arbitrary Lagrange-Euler)描述下的三维Euler方程：(4)其中，Ω为控制体，?Ω为控制边界，n为控制体边界外法向向量，dσ为体积微元，ds为面积微元，xc为网格面的运动速度。把均匀流场参数代入，得到积分形式的GCL方程：(5)根据式(5)，对于静止网格xc=0，有限体积法不存在GCL问题，即使是运动网格，如果没有变形，式(5)也自动满足。国内采用结构网格与物体刚性固连运动开展动导数研究的文献[5-8]和相关学术专著[9]采纳了以上结论。但是，从本文推导来看，这种结论仅对梯度重构界面通量的有限体积法成立，对于采用有限差分法重构界面通量的高精度格式不一定正确。本文后面的例子表明，所有包含有四边形的网格均存在几何参数计算不相容引入误差的风险。下面考察有限差分法产生GCL问题的数学本质。推导物理空间(t,x,y,z)的守恒控制方程变换到计算空间(τ,ξ,η,ζ)得到如下原始形式：=UIt+FIx+GIy+HIz=0(6)只是在时间和空间导数连续条件下，使用了式(3)及以下恒等式：(7)从理论上讲，带有源项的方程(6)才是原始方程(1)的等价形式，但坐标变换系数大多采用网格位置的差分来近似求解，此时恒等式(3)和式(7)可能不成立，在这种情况下方程(2)与方程(1)不等价，把它作为曲线坐标系下的等价形式进行CFD计算，得到的数值解必然存在误差。研究GCL算法的目的是消除由于几何参数计算引起的误差。在GCL研究领域，把涉及网格变形和时间导数的式(3)称为体积守恒律，把仅有空间导数的式(7)称为面积守恒律[10]。动网格和静网格均须考虑面积守恒律，这点已受到广泛关注，亦是目前GCL研究的主要内容。下面对相关研究文献进行综述。有限