三维有限体积TVD方法和冲击波的多级扩散研究.pdfVIP

16 三维有限体积TVD方法与冲击波的多级扩散研究 三维有限体积TvD方法与冲击波的多级扩散研究 张海波L2白春华1丁儆1郝保田2 (1北京理工大学-北京100081)P西北核技术研究所，西安 710024) 埔蔓本文给出可压缩流体动力学中的三维有限体积TvD方法，并采用该方法对三维冲击 波的多级坑道扩散问思进行了数值模拟．取得了々人满意的结果． 关t蔼三堆TIfD方怯有限体积法 教值模拟 冲击波 1引言 在计算流体动力学中．TvD方法【1]虽已得到广泛应用．但研究集中在一维、二维情况， 关于三维有限体积m方法的报道，尚未见到。本文在三维直角坐标系中，采用有限体 积法口]，将可压缩流体动力学方程组进行离散化．然后，将ⅣD方法应用到三维坐标， 给出三维有限体积1VD方法。采用有限体积法．差分方程是通过对微分方程在有限体积 的微元(体元)上积分得到的，这使得有限体积法得到的差分方程．适用于不规则网格， 即，可用于计算具有不规则几何形状的流场． 根据以上方法，我们运用面向对象编程技术和C++语言，编制完成TVD3D面向对象程 序，首先进行了平面激波测试，模拟了平面激波的传播、正反射，发现对激波的分辨率 较高，波阵面约l～2个网格。波后平整．激波的反射压强与理论结果的相对误差小于1％。 然后模拟了冲击波的多级坑道扩散这一复杂的三维物理问题．研究了冲击渡的绕射、传 播，马赫反射、正规反射等现象．计算结果表明，三维有限体积TVD方法的间断分辨率 高，光滑区的精度高，且计算稳定。数值结果与实验结果比较，相当一致．这说明三维 有限体积TvD方法是一种根好的计算方法。将TvD方法应用到三维坐标，是成功的。 理想可压缩流体运动的三维直角坐标Euler型守恒方程组为 。(1)’ aU，珧+aP(U)，缸+aS(U)，咖+跏KU)，ez=0 其中U，F(u)，GU)．M(U)均为列向量， 张海激：男．1969年1月生，博士 ——茎三旦笪茎塑墨兰查叁堡 !! Pl l∥ 『 『∥ 舢 绷 P+绷2 pvu pwu U= F=删V G=P+∞2 p~ Mt C例 删 graw pvw P+pw2 DE “(卢匹+，) V(卢目十尸) w(卢匠+P) E为单位质量的总能量，E=F+02+y2+w2)，2 (3) 上式中，工，Y，g为三维直角坐标，P，p，Ⅳ，v，w分别为流场中介质的压 强-密度，x方向速度，Y方向速度，z方向速度，g为单位质量的内能，t为时间。 状态方程P=(y一1)p(E一02+v2+w2)／2) (4) 采用有限体积法，将微分方程(1)在有限体积的微元(体元)V上积分，得 (s) 岳肌鼢啦+m(芸+署+警)dxdydz=0 运用显式时间差分方法，对任意六面体网格．将(5)式离散化，我们得到(1)式对应 的有限体积法差分方程为 (6) u现=u‰一争∑(F4帆+研4Np+研4心) u’JJ^f-I