自适应网格局部加密方法在可压缩流场模拟中的应用

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/     

局部网格加密方法计算不可压N—S方程

本文发展了二维非定常不可压N-S方程的局部网格加密方法(LMR):需加密区域预先指定;仅对空间方向加密,分粗细二层网格,细网格覆盖在粗网格上;粗和细网格上分别用显式和隐式差分格式。结合压力修正法类的余量型差分格式,恰当地处理了粗细网格之间的信息传递,使得粗细网格交接面上质量守恒,即满足连续性条件。粗网格通过插值给细网格边界值影响细网格;细网格通过粗细网格压力Poisson方程耦合求解影响粗网格,并且由于压力Poisson方程从动量方程形成,在细网格覆盖下的粗网格上,压力Poisson方程是细网格方程的组合,从而耦合求解时粗细网格压力Poisson方程不需迭代。 本文中计算了二维方腔内的自然对流模型问题,Rayleigh数是10~6。粗网格用显式ULWC格式,细网格用隐式余量型近似因式分解格式。Poisson方程用快速直接算法PO1STG(在FISHPACK中)。粗网格距△=1/16,细网格距△=1/64,加密1/4区域,计算到定常局部网格加密方法所需CPU时间比同等均匀网格(△=1/64)节约一倍以上。