二阶无振荡特征格式在对流扩散方程中的应用.pdfVIP

二阶无振荡特征格式 在对流扩散方程中的应用 于恒 忻孝康 沈一帆 (复旦大学应用力学系， 摘要本文就已发展的一致二阶精度的无振荡特征格式…，进一步捐展丁其无振荡区 域．并把此原来适用于取曲型方程的方法推广到对流扩散方程的汁算。利用适当的边界条件延 拓和时间相关法，成功地解决了定常及非定常对流扩散方程的初、边值问题的计算，对一系列 较典型而困难的对流扩散问题，给出了较为满意的数值结果。 关键词 无振荡格式，对流～扩散方程，数值计算 1引言 众所周知，在对双曲型守恒方程进行差分计算时，二阶(或更高阶)数值格式在间 断面附近总是伴有非物理的数值振荡，它可以通过添加人工耗散来加以克服，但这样一 来便大大降低了数值解的精度和对间断的分辨率。Harten的高分辨率TVD格式“’，开 创了寻求无振荡高阶数值格式的新篇章。但具有间断高分辨的压缩性TVD格式，常常将 梯度较大的光滑解压缩成为间断解，并且几乎所有的TVD格式都在局部极值点退化为一 阶精度。 为了寻求无振荡的一致高精度格式，对于双曲型方程的初值问题 “，+，(“)，=0，u(x，O)=”o(x) (1) 机制作了详尽的分析之后得出：只要二阶特征插值的抛物线极值点不落入插值区域，则 该格式此时不会产生过靶现象。于是得到二阶特征格式的无振荡区域(口z谚／护“≥0 时)： Lax-wendr。ff格式：陲l，2—1；≤2 (2．1) 二阶迎风格式：{t。，：一1j≤2 (22) ll 其嘞=鬻， 一■ 瓣。由此启发得到了～种二阶无振荡特征格 式，它不带有复杂的TVD限制器，但却是无振荡的，并且具有～致二阶精度。该格式用 于双曲型方程的数值计算，得到了非常满意的结果。 ·30· 本文将该格式的无振荡区域进一步拓展，并分析了其TVD性质，然后将此原来适用 于双曲型方程的格式，推广到定常及非定常对流扩散方程的数值计算。特别对定常问 题，由于该格式往往涉及四点或五点格式，如何处理边界条件是一大问题，本文通过采 用合理的边界条件延拓，利用时间相关法，成功地解决了定常和非定常对流扩散方程的 初、边值问题。文中给出了一系列较典型而困难问题的数值结果，结果是令人满意的。 2拓展的无振荡区域 找出了二阶特征格式数值振荡的几何机制之后，受此启发，我们可进～步得到如下 结论：对于双曲型方程(1)的任何二阶特征数值格式，只要数值解满足 G一=(彬“一《)(彬“一“?～，)≤0 当口≥0时 (3．1) G+=(“；“一¨?)(彬“一≯鼻，)≤o当口0时 (3．2) 则此时格式不会产生过靶现象。 公式(3．I)和(3．2)的几何意义是显而易见的，它要求插值解落在左右插值节点的被 插函数值之间，不大予昂大值，也不小于最小值。将Lax—wendroff格式和二阶迎风格 式分别代入条件(3．I)和(3．2)，容易得到相应的最大无振荡区域(见表．I)。 表1二阶特征格式的最大无振荡区间 风向 dO d0 格式 LUP L—Ⅵ L-W RIP 判断参数 ～ — 0—1，2 r二l，2 0+】，2 t1，2 7 ’ 熊凝刚礴付嘶徽h一南9一s萄 显见它对任意CFL数盯(o≤J仃l≤1)都满足的最小无振荡区间(口≥0时)正是 (2．1)和(2．