小位移运动物体三维贴体直角网格的流场数值近似计算-期刊.pdfVIP

小位移运动物体三维贴体直角网格的流场 数值近似计算 李盾 航夭空气动力技术研究院，北京100074 摘要：利用直角弼格的空阊不变的特点，提出了运动豹磷辩，采用静态空间网格、杨面边界运动钓三维 非结构贴体直角网格的数值方法进行非定常流场的近似计算。提出的诗算方法在窑向上完全满足运动刚格 的N．S方程的通量守恒条件．和几何守恒条件，在小位移假{殳的条件下对物面网格进行运动边界的近似处 理，逼近{#定常过程。通过单体运动和多体运动的数值计算的计算对托，说明提出的新方法可以反映物体 遥动的瞬间趋势，在小鞋闻步长的情况下，计算的气动力特性更加接近真实情况。为变形体运动等小位移 运动过程的三维粘性流场计算提供了一个可能的解决方案。 关键词：粘性直角网格相对运动数值模拟运动边界 1、引 言 物体的运动可以分为大位移运动过程和小位移运动过程，对于物体的小位移过程，运动 阿格的生成主要有插值方法、网格拉伸的弹簧近似方法、弹性体方法等。 运动物体非定常过程静流场数值模拟计算量很大，一般难以实际成用，实际应用主要采 取近似方法处理，对于相对运动的非定卷流场计算问题，大致可以分为几种近似：动嘲格的 完全非定常的数值计算；本文提出的三维非结构贴体直角固定潮格、运动边界的数值计算， 这种方法针对非结构矗角网格的特点，保证绝大部分流场的物理守恒和网格的几何守恒，在 运动物体的瞬间，考虑了流场每一点螅情况，对边界流场进行了一定盼近似，大大筒仡了非 定常过程，适用于三维复杂物形的快速计算，非定常过程的近似程度较高；铡定嬲格固定边 界的攻角修正的定常计算，这种方法采用流场平均的方法进行攻角修正，没有考虑物体每一 点的流场情况，主要适用于简单流场的快速计算，可以近似非定常过程； 固定网格固定边 界的完全定常计算，这种方法是采用定常{{。算代替非定常过程，非定常过程的近似程度较低。 本文研究了三维非结构贴体直角固定网格、运动边界的数值计算方法，并进行了数值分 析，说明这种近似方法可以在一定程度上可以反映流场的运动情况，为变形体运动等小位移 运动过程的三维物体流场计算提供了一个可能的解决方案。 2、运动物体三维非结构直角网格的流场数值计算 2．1相对运动物体流场的三维非定常流动控制方程 笛卡尔坐标系下，流体相对坐标系的速度为U，物体相对坐标系是非静止的，某一微小 边界运动的控制体口rf，也在运动，控翻体的边界a口ff，有质量以及这些质量携带的动量、 能爨等等物理嶷通量流过，因此考虑的控制体的运动，由质量守恒、动量守恒及能量守恒 可以推得网格及控制体相对运动的N-S方程的积分形式为： 怯黔帆册cU州∥麟‘0 陆鼢UdV+野叭U_u卜P}nd8。翳mds U)】-趟鼯础f』，[(z-U)-q]．n瓠 lI昙黔州+册(U州。“p 其中P、q为流体静压张量和热流量，U。为单元的边界速度。已知边界的 运动速度，则流体通过边界的相对运动速度为u一％。 考虑了物体的耜对运动，则无量纲的N—s方程可以改写为： 昙腑dv+骶t赫=俨D．．．dS (2) 。n a0 孙 其中： p(V—Va)p(w—wn)- p(V—Va)up(w～w。)ui WD= K= p(v—Vd)V+pp(w—wd)Vl (3) p(v—Vd)wp(w—Wd)w1+d p凹∥眇E (琶十p)(V—Vd) 孵”陋陋啪茹=卅 (E十p)(w—wd)U 蚶叶们咖一 粘性项： 0 0 t瓢 tw