非均匀网格三点四阶紧致差分格式.pdfVIP

第十届全国计算流体力学会议论文集 非均匀网格三点四阶紧致差分格式 于欣 中国科学院力学研究所，北京100080 摘要本文对一阶、二阶导数给出非均匀网格上的三点四阶精度帑致差分离散。求解 二阶微分方程二维或三维问题．精度高、耗散色敬低、计算稳定．适用于湍流直接数值模 拟等(见川，f21)。 0引言 随着电子计算机的发展．越来越多的实际问题数值模拟成为现实，但同时，实际问题 对计算的要求也越来越高．有些问题可能要计算多日甚至更长。在保证计算精度的前提下， 减少计算时间的一个方法是提高离散精度。 例如对直接数值模拟，计算机容量和速度需要极大，计算方法受到挑战，由于湍流 中包含不同尺度的结构，而大小尺度结构相差十分悬殊，这就要求数值计算方法同时具有 捕捉大小尺度流动结构的能力。具有高精度+低耗教低色散，计算稳定，具有抑制混淆误 差的能力。 近十年来，紧致差分格式得到更广泛的应用，她比传统差分格式具有使用的网格点 少，精度高，低耗敬低色敖，稳定性好等优点．更适合于湍流直接数值模拟。 我们在【3】中给出了交错网格紧致差分格式，他比非交错网格紧致格式计算更稳定。 过去的紧致差分格式通常采用均匀网格。对非均匀网格，通过Jacobi变换使计算区 域网格均匀。这在物理域网格变化剧烈的情况下会产生较大误差。另外，有时由于计算需 要，网格变化很大。LC．mmet等人于1999年提出了基于非均匀网格的对称紧致格式．该 格式甚至在随机网格上仍能保持很高精度，见【l】。 J．奉项研究最新进展见： ed蛐un“ b却：／／mvwoer∞．psu cfm．brown 个人主页：lattp：t／ww,v edu／faculty／yu／ 2．作者的Email：”@cerse．psu．edu 3．作者感谢中国科学院回国人员择优基金的资助 李新亮等A构造了基于非均匀网格的迎风紧致格式【21。文章用数值实验说明，该方法 精度高，对网格要求宽松，能抑制高额非物理振荡，可推广到三维情况，是进行槽道湍流 直接数值模拟的有效方法。 我们这里给出的基于非均匀阿格的紧致差分格式- l 对于一阶导数：与LGAMET文(见【1])中的四阶格式相比．精度相同(四阶)t 系数简单，差分点少【三个差分点，而LGAMET格式五个差分点；李新亮文[2】 中的迎风紧致格式为五点五阶)； 2．对于二阶导数：与[1】中的三阶格式相比，我们这里给出的基于非均匀网格的紧致 格式，精度高一阶(四阶)，差分点少(三个差分点，而L．GAMET格式五个差 分点)。 本文给出的紧致差分格式是对一阶和二阶导数的离散．因此可以应用到二阶微分方程 和偏微分方程二维和三维问题。 1非均匀网格三点四阶紧致差分离散 对—阶导蠡，我们给出如下基于非均匀阿格的紧致差分格式 db““鼎铲黼莆+黼等 其中 bm=k一1／^=(ze一“．)／^ 岵=肌，h=(zP—zc)lh u。，u0．u麓在q—l=￡。点取值 u，“7，u”在Zi=‰点取值 邯，嵋，“；在 2i+i=。口点取值 对角占优- 1 蚱 6矗 2b。b 一砾了耵一瓦可了2砾葡 。： 与L．GAMET丈【ll中的四阶格式相比，精度相同(四阶)，系数俺单，差分点少(三 十差分点。而L．GAMET格式五个差分点)． 上述一阶导敷的差分式中u。，u，邯也可记为Ⅵ。Ⅵ，啦+。，用光滑爵数口代入该式得