求解双曲守恒律方程的三阶修正模板WENO格式北大核心CSCD

为了降低经典的三阶加权本质无振荡(WENO)格式的数值耗散,提出了一种新的三阶WENO格式的修正模板近似方法.改进了经典WENO-JS格式中各候选模板上数值通量的一阶多项式逼近,通过加入二次项使模板逼近达到三阶精度.计算了相应的候选通量,并且通过引入可调函数φ(x),使得新的格式具有ENO性质.最后给出了一系列数值算例,证明了该方法的有效性.     

基于新光滑因子的WENO5格式

针对WENO格式的构造,本文给出了一个WENO为5阶的充分条件,降低了Henrick等人提出的充分条件对于权因子的精度要求.另外,对于Jiang和Shu提出的WENO5中的光滑因子中两项的系数做出了调整,并结合Borges等人的方法得到了新的WENO权因子计算方法.从数值试验的结果可以看出,新的WENO格式对于连续波形...     

激波捕捉差分方法研究

在迎风型格式和矢通量分裂技术的基础之上,对捕捉激波方法进行一种新的尝试．该方法首先对原始格式在特征方向上进行投影,然后用限制器对这些特征分量的变化幅值进行限制以抑止非物理波动,最后再把它转换成守恒形式,得到了基本上无振荡的激波捕捉格式．用该方法对两种迎风显示格式(二阶和三阶)和3种迎风紧致格式(三阶、五阶和七阶)进行处理,并在一维和二维的情况下进行了应用测试．通过与高阶WENO、MP、Compact-WENO等格式的比较,表明该方法在光滑捕捉激波的前提下仍有较高精度和分辨率．     

基于高阶空间精度格式和嵌套网格对旋翼尾迹捕捉的改进

发展了一种基于高阶迎风格式和嵌套网格捕捉直升机悬停旋翼涡尾迹的方法．无粘通量采用Roe Reimann求解器,使用改进的5阶加权基本无振荡(WENO)格式对交界面左右状态进行高阶插值,并与MUSCL插值进行比较．为便于捕捉尾迹和实施周期性边界条件,计算采用结构嵌套网格,其中高质量的旋翼网格完全嵌套于背景网格中．当解达到近似收敛后在桨尖涡分布区域对背景网格进行加密,如此经过3次得到优化的背景网格．考虑到WENO格式插值的特点,提出了搜索3层洞边界和人工外边界的方法以便插值的直接进行．用该方法对一跨音速和一亚音速悬停旋翼粘性流场进行了数值计算．数值结果表明:所发展方法对涡尾迹具有很高的捕捉能力;与MUSCL格式相比,WENO格式具有较低的数值耗散．     

求解一维对流方程的高精度紧致差分格式

本文提出数值求解一维对流方程的一种两层隐式紧致差分格式,采用泰勒级数展开法以及对截断误差余项中的三阶导数进行修正的方法对时间和空间导数进行离散.格式的截断误差为O(τ~4+τ~2h~2+h~4),即该格式在时间和空间上均可达到四阶精度.利用von Neumann方法分析得到该格式是无条件稳定的.通过数值实验验证了本文格式的精确性和稳定性.     

一个改进的三阶WENO-Z型格式

在WENO-Z型格式框架下,基于高阶全局光滑因子,在非线性权建立过程中引入参数,通过收敛性分析确定参数取值范围,兼顾精确性与不振荡性,得到参数最佳取值.最终得到一个低耗散、高分辨率的三阶WENO差分格式,该格式在函数一阶极值点处仍保持预期三阶精度.最后通过精确解算例验证了格式在各种类型极值点处精度恢复情况,并通过一、二维Euler方程组经典算例测试了格式的低耗散、高分辨特性.结果表明,该文格式是一个性能优良的激波捕捉格式.     

非线性发展方程的小模板简化Pade格式

在有理逼近的紧致格式的理论基础上，采用特别的统一的Pade逼近形式，构造了针对高阶非线性发展方程的、简单小模板的差商格式．不仅保持了格式的四阶精度，而且还可以采用追赶法求解得到的3对角矩阵，或者采用三阶Runge-Kuna法直接求解积分．计算效果通过多种算例表明是十分令人满意的．相对于其他差分格式，此方法具有模板较小而精度保持四阶的优点．     

基于映射函数的三阶WENO改进格式及其应用

低耗散、高分辨率激波捕捉格式对含激波流场的数值模拟具有重要意义.在传统三阶WENO格式（WENO-JS3）和三阶WENO Z格式（WENO-Z3）基础上,基于映射函数,给出WENO-M3、WENO MZ3格式.选用Sod激波管、激波与熵波相互作用、双爆轰波碰撞及双Mach(马赫)反射等经典算例,考察上述格式的计算性能.数值结果表明,WENO-MZ3格式相较其他格式具有耗散低、对流场结构分辨率高的特性.为了进一步扩展WENO-MZ3格式的应用范围,采用该格式数值研究封闭方形舱室内柱形高压、高密度气体爆炸波传播过程,波系演化规律以及壁面典型测点压力载荷.数值计算结果表明WENO-MZ3格式能够较好地模拟包含高压比、高密度比的爆炸波且给出数值耗散较小的壁面压力载荷.     

求解时间分布阶扩散方程的两个高阶有限差分格式

基于复化Simpson公式和复化两点Gauss-Legendre公式，构造了两个求解时间分布阶扩散方程的高阶有限差分格式.不同于以往文献中提出的时间一阶或二阶格式，这两种格式在时间方向都具有三阶精度，而在分布阶和空间方向可达到四阶精度.数值结果表明，两种算法都是稳定且收敛的，从而是有效的.两种格式的收敛速率也通过数值实验进行了验证，并且通过和文献中的算法对比可以得出其更为高效,     

高精度隐式参差光滑格式的构造

参照Lax-Wendroff格式的构造方法,就双曲型方程、抛物型方程和双曲-抛物型方程,构造了一种新的IRS(implicit residual smoothing)格式。该IRS格式有二阶或三阶时间精度且大大地拓宽了解的稳定区域和CFL数。这种新的IRS格式有中心加权型和迎风偏向型两种,并用von-Neumann分析方法分析了格式的稳定范围。讨论了在透平机械中广泛应用的Dawes方法的局限性,发现该方法对稳态问题得出的解与时间步长的选取有关,对粘性问题求解时,时间步长受严格限制。最后,结合TVD(total variation diminishing)格式和四阶Runge-Kutta技术,用IRS格式和Dawes方法对二维反射激波场进行了数值模拟,数值结果支持本文的分析结论。     

一种新的求解双曲型守恒律的三阶中心差分格式

本文提出了一种求解双曲型守恒律新的三阶中心差分格式，主要是引入了一种推广的三阶重构，并证明了这种重构在网格边界无振荡．所提的格式保持了中心差分格式简单的优点，不需用Riemann解算器，避免了进行特征解耦．数值试验结果表明本文格式是高精度、高分辨率的。     

An adaptive mesh method for 1D hyperbolic conservation laws

An adaptive method is developed for solving one-dimensional systems of hyperbolic conservation laws, which combines the rezoning approach with the finite volume weighted essentially non-oscillatory (WENO) scheme. An a posteriori error estimate, used to equidistribute the mesh, is obtained from the differences between respective numerical solutions of 5th-order WENO (WENO5) and 3rd-order ENO (ENO3) schemes. The number of grids can be adaptively readjusted based on the solution structure. For higher efficiency, mesh readjustment is performed every few time steps rather than every time step. In addition, a high order conservative interpolation is used to compute the physical solutions on the new mesh from old mesh based on the finite volume ENO reconstruction. Extensive examples suggest that this adaptive method exhibits more accurate resolution of discontinuities for a similar level of computational time comparing with that on a uniform mesh.     

基于新的参考光滑性指示子的改进的三阶WENO格式北大核心CSCD

针对计算流体力学对高精度高分辨率的需求,基于降低经典的三阶加权本质无振荡(WENO)格式的数值耗散特性,该文提出了一种新的参考光滑性指示子.其构造方法与经典的WENO-Z格式不同,它是通过候选子模板上重构多项式的导数的线性组合与整个全局模板上重构多项式的导数的L^(2)范数逼近获得的.采用该计算方法可以得到比WENO-Z格式更高阶的参考光滑性指示子,另外改变自由参数φ的取值,可以获得不同的参考光滑性指示子.该文通过一系列数值算例证明了该参考光滑性指示子的有效性.     

On the use of Mühlbach expansions in the recovery step of ENO methods

Summary. The recovery step is the most expensive algorithmic ingredient in modern essentially non-oscillatory (ENO) shock capturing methods on triangular meshes for the numerical simulation of compressible fluid flow. While recovery polynomials in Newton form are used in one-dimensional ENO schemes it is a priori not clear whether such useful as well as numerically stable form of polynomials exists in multiple dimensions. As was observed in [1] a very general answer to this question was provided by Mühlbach in two subsequent papers [15] and [16]. We generalise his interpolation theory further to the general recovery problem and outline the use of Mühlbach's expansion in ENO schemes. Numerical examples show the usefulness of this approach in the problem of recovery from cell average data. Received August 24, 1995 / Revised version received December 14, 1995     

WENO schemes with Lax–Wendroff type time discretizations for Hamilton–Jacobi equations

In this paper, a class of weighted essentially non-oscillatory (WENO) schemes with a Lax–Wendroff time discretization procedure, termed WENO-LW schemes, for solving Hamilton–Jacobi equations is presented. This is an alternative method for time discretization to the popular total variation diminishing (TVD) Runge–Kutta time discretizations. We explore the possibility in avoiding the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong discontinuous derivative. As a result, comparing with the original WENO with Runge–Kutta time discretizations schemes (WENO-RK) of Jiang and Peng [G. Jiang, D. Peng, Weighted ENO schemes for Hamilton–Jacobi equations, SIAM J. Sci. Comput. 21 (2000) 2126–2143] for Hamilton–Jacobi equations, the major advantages of WENO-LW schemes are more cost effective for certain problems and their compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method.     

Third order nonoscillatory central scheme for hyperbolic conservation laws

Summary. A third-order accurate Godunov-type scheme for the approximate solution of hyperbolic systems of conservation laws is presented. Its two main ingredients include: 1. A non-oscillatory piecewise-quadratic reconstruction of pointvalues from their given cell averages; and 2. A central differencing based on staggered evolution of the reconstructed cell averages. This results in a third-order central scheme, an extension along the lines of the second-order central scheme of Nessyahu and Tadmor \cite{NT}. The scalar scheme is non-oscillatory (and hence – convergent), in the sense that it does not increase the number of initial extrema (– as does the exact entropy solution operator). Extension to systems is carried out by componentwise application of the scalar framework. In particular, we have the advantage that, unlike upwind schemes, no (approximate) Riemann solvers, field-by-field characteristic decompositions, etc., are required. Numerical experiments confirm the high-resolution content of the proposed scheme. Thus, a considerable amount of simplicity and robustness is gained while retaining the expected third-order resolution. Received April 10, 1996 / Revised version received January 20, 1997     

Mapped WENO schemes based on a new smoothness indicator for Hamilton–Jacobi equations

In this paper, we introduce an improved version of mapped weighted essentially non-oscillatory (WENO) schemes for solving Hamilton–Jacobi equations. To this end, we first discuss new smoothness indicators for WENO construction. Then the new smoothness indicators are combined with the mapping function developed by Henrick et al. (2005) [31]. The proposed scheme yields fifth-order accuracy in smooth regions and sharply resolve discontinuities in the derivatives. Numerical experiments are provided to demonstrate the performance of the proposed schemes on a variety of one-dimensional and two-dimensional problems.     

Generalizing the ENO-DB2<Emphasis Type="Italic">p</Emphasis> transform using the inverse wavelet transform

The essentially non-oscillatory (ENO)-wavelet transform developed by Chan and Zhou (SIAM J. Numer. Anal. 40(4), 1369–1404, 2002) is based on a combination of the Daubechies-2p wavelet transform and the ENO technique. It uses extrapolation methods to compute the scaling coefficients without differencing function values across jumps and obtains a multiresolution framework (essentially) free of edge artifacts. In this work, we present a different way to compute the ENO-DB2p wavelet transform of Chan and Zhou which allows us to simplify the process and easily generalize it to other families of orthonormal wavelets.     

无波动,无自由参数,耗散的隐式差分格式

本文建立了求解NS方程和Euler方程无波动、无自由参数、耗散的隐式差分格式.该格式是TVD的和无条件稳定的.其隐式部分在1,2,3维情况下仅分别依赖于3,5,9个点,且系数矩阵是主对角占优的.计算例题表明,该方法可获得和显式方法相同的精度,能很好地捕捉激波和剪切层,且计算时间比显式有较多的节省.