非结构网格的高精度有限体积离散格式

本文针对单元中心非结构网格有限体积算法,通过构造上游虚节点和局部一维坐标,将结构化网格中常用的规正变量(NVSF)和总变差消失通量限制器(TVD FL)两种形式的传统高精度格式移植到非结构网格中。为了提高虚节点的插值精度,使用最小二乘法计算单元中心的变量梯度。经典算例考核并与商业软件FLUENT结果的比较表明,本文对非结构网格的高精度格式移植成功有效。     

不规则对接网格的分区求解计算

采用一种保持通量守恒的不规则对接网格分区求解中交界面耦合条件的计算方法, 结合有限体积法求解了Euler 方程, 无粘通量取用Van Leer 分裂格式, 构造了一种限制器以实现格式的二阶精度和TVD 性质, 并给出数值算例。     

二维欧拉多介质流体力学数值模拟-TVD＋VOF方法

用TVD格式结合VOF界面处理方法编制了二维多介质高分辨欧拉程序，以解决冲击波和多介质界面处理。程序包括单介质网格高精度流体力学计算、多介质网格内界面重构、各种介质输运和压力驰豫平衡过程。其中单介质网格的计算采用Harten二阶TVD格式结合MacCormark方法计算含有源项的非齐次守恒定律方程组，通过4节点限制函数保证格式单调。多介质网格采用Youngs方法构造界面，采用x，y方向分裂格式计算体积份额输运，再根据体积份额输运计算质量、动量和能量的输运，最后利用等熵条件计算各种介质的压力驰豫平衡过程。     

间断Galerkin方法求解跨音压气机转子流场

本文采用间断Galerkin方法求解雷诺平均N-S方程和SA湍流模型方程,采用Roe格式计算无黏通量,混合格式求解黏性通量中的梯度。为防止跨音速流场计算中的数值振荡,采用了TVD限制器和正性修正。计算了跨音速转子NASA Rotor37在设计转速下的变工况流场,得到了与实验吻合良好的结果,表明DG方法在叶轮机械内部流动计算中具有广阔的应用前景。     

求解可压缩流的高精度非结构网格WENO有限体积法

提出-种基于最小二乘重构和WENO限制器的非结构网格高精度有限体积方法.用中心网格的某些邻居网格建立重构多项式,给出-定的原则搜索和存储足够多的邻居网格以建立重构多项式,采用最小二乘法求解重构多项式的系数.用-种通用的方法控制重构邻居个数,以减少存储和计算,采用WENO限制器和旋转Riemann求解器以达到统-的高精度并且抑制守恒律方程求解中的非物理振荡.为检验上述算法,以基于节点的梯度重构,Bath and Jesperson限制器的二阶算法为基准,给出三阶和四阶格式与二阶格式以及高阶格式若干经典算例计算结果的对比和分析.     

反应流体的移动网格动理学格式

在移动网格上构造一种反应流的动理学格式.首先利用BGK模型推导含化学反应的流体力学方程组,并利用其积分形式构造移动网格上离散格式,再利用自适应移动网格方法得到网格速度,最后利用时间精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.一维与二维的数值实验表明这种格式同时具有高精度、高分辨率的特点.     

CMYK防伪印刷微镜阵列

以印刷微镜的反射亮度、反射稳定性和CMYK颜色模式为基础,建立了一种新型的CMYK印刷微镜结构模型。利用微镜表面函数的雅可比行列式和梯度函数,对模型的表面反射亮度和稳定性进行了分析评价。结果表明:该结构模型表面函数的雅可比行列式近似为0,梯度函数随坐标变化而变化。因此,此结构模型表面的反射光具有足够大的亮度和稳定性,适用于创建基于位置和角度变化来改变光变图案反射密度的光学防伪效果。设计了一幅由不同CMYK印刷微镜组成的防伪光变图案,模拟了该防伪图案的印刷微镜阵列结构。     

求解多维欧拉方程的二阶非结构网格混合旋转Riemann求解器

将基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法推广到非结构网格,采用基于网格中心的有限体积法,梯度的计算采用基于节点的方法引入更多的控制体模板,限制器的构造采用与非结构化网格相适应的形式.在求解Riemann问题时,沿具有一定物理意义的两个迎风方向,即控制体界面两侧速度差矢量方向及与之正交的方向.能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或"红斑"现象.为减小计算量,采用HLL和Roe FDS混合旋转格式.     

双曲型守恒律的一种高精度TVD差分格式

构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例.     

一维含化学反应流体力学方程组的ALE间断有限元方法

研究一维含化学反应流体力学方程组的数值模拟方法.结合理想气体状态方程并利用HLLC解法器在各个单元边界处的数值通量,给出ALE间断有限元方法.高阶计算时,使用TVD斜率限制器对数值解可能产生的非物理振荡进行抑制.结果表明：该算法能够保持物理量的守恒性和高精度,并能够清晰地捕捉爆轰波的结构特征.     

Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics

In a pair of earlier papers the author showed the importance of divergence-free reconstruction in adaptive mesh refinement problems for magnetohydrodynamics (MHD) and the importance of the same for designing robust second order schemes for MHD. Second order accurate divergence-free schemes for MHD have shown themselves to be very useful in several areas of science and engineering. However, certain computational MHD problems would be much benefited if the schemes had third and higher orders of accuracy. In this paper we show that the reconstruction of divergence-free vector fields can be carried out with better than second order accuracy. As a result, we design divergence-free weighted essentially non-oscillatory (WENO) schemes for MHD that have order of accuracy better than second. A multi-stage Runge–Kutta time integration is used to ensure that the temporal accuracy matches the spatial accuracy. While this is achieved quite simply up to third order in time, going beyond third order is most simply achieved by using the ADER-WENO schemes that are detailed in a companion paper. (ADER stands for Arbitrary Derivative Riemann Problem.) Accuracy analysis is carried out and it is shown that the schemes meet their design accuracy for smooth problems. Stringent tests are also presented showing that the schemes perform well on those tests.     

对流项离散格式的对比与讨论

本文简单介绍了利用规正变量定义的各种对流项差分格式,给出了利用有限容积法离散粘性对流一扩散问题时的离散方程,其中的对流项采用高阶格式进行离散。以方腔顶盖驱动及圆管突扩区内层流流动考察了各种格式的计算精度与时效。通过对比分析得出:对于常规区域中的流动, QUICK、中心差分(CD)及SMART三种格式的精度与计算时效是比较合理的。     

Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy

In this paper we design a class of numerical schemes that are higher-order extensions of the weighted essentially non-oscillatory (WENO) schemes of G.-S. Jiang and C.-W. Shu (1996) and X.-D. Liu, S. Osher, and T. Chan (1994). Used by themselves, the schemes may not always be monotonicity preserving but coupled with the monotonicity preserving bounds of A. Suresh and H. T. Huynh (1997) they perform very well. The resulting monotonicity preserving weighted essentially non-oscillatory (MPWENO) schemes have high phase accuracy and high order of accuracy. The higher-order members of this family are almost spectrally accurate for smooth problems. Nevertheless, they, have robust shock capturing ability. The schemes are stable under normal CFL numbers. They are also efficient and do not have a computational complexity that is substantially greater than that of the lower-order members of this same family of schemes. The higher accuracy that these schemes offer coupled with their relatively low computational complexity makes them viable competitors to lower-order schemes, such as the older total variation diminishing schemes, for problems containing both discontinuities and rich smooth region structure. We describe the MPWENO schemes here as well as show their ability to reach their designed accuracies for smooth flow. We also examine the role of steepening algorithms such as the artificial compression method in the design of very high order schemes. Several test problems in one and two dimensions are presented. For multidimensional problems where the flow is not aligned with any of the grid directions it is shown that the present schemes have a substantial advantage over lower-order schemes. It is argued that the methods designed here have great utility for direct numerical simulations and large eddy simulations of compressible turbulence. The methodology developed here is applicable to other hyperbolic systems, which is demonstrated by showing that the MPWENO schemes also work very well on magnetohydrodynamical test problems.     

Finite difference methods for second order in space,first order in time hyperbolic systems and the linear shifted wave equation as a model problem in numerical relativity

Motivated by the problem of solving the Einstein equations, we discuss high order finite difference discretizations of first order in time, second order in space hyperbolic systems. Particular attention is paid to the case when first order derivatives that can be identified with advection terms are approximated with non-centered finite difference operators. We first derive general properties of these discrete operators, then we extend a known result on numerical stability for such systems to general order of accuracy. As an application we analyze the shifted wave equation, including the behavior of the numerical phase and group speeds at different orders of approximations. Special attention is paid to when the use of off-centered schemes improves the accuracy over the centered schemes.     

加权ENO格式的构造及数值模拟

将加权ENO格式推广到非结构三角形网格上,构造了一类加权ENO有限体积格式,提出的插值多项式的构造方式,可以减少计算时间.对于出现的病态方程组,给出了解决方法.此外还给出了插值点的选取方式及加权因子的构造方法.结合三阶TVD Runge Kutta时间离散,对二维欧拉方程组进行了数值试验.     

On maximum-principle-satisfying high order schemes for scalar conservation laws

We construct uniformly high order accurate schemes satisfying a strict maximum principle for scalar conservation laws. A general framework (for arbitrary order of accuracy) is established to construct a limiter for finite volume schemes (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO) schemes) or discontinuous Galerkin (DG) method with first order Euler forward time discretization solving one-dimensional scalar conservation laws. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. It is straightforward to extend the method to two and higher dimensions on rectangular meshes. We also show that the same limiter can preserve the maximum principle for DG or finite volume schemes solving two-dimensional incompressible Euler equations in the vorticity stream-function formulation, or any passive convection equation with an incompressible velocity field. Numerical tests for both the WENO finite volume scheme and the DG method are reported.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

Poisson方程的高精度有限差分格式

本文提出数值求解Poisson方程的含选择因子α的预示校正差分格式,它具有四阶精度。第一种格式处理Drichlet边界条件的Poisson方程,它包括Bramble的差分格式和林群、吕涛提出的差分外推格式。第二种格式处理Neumann边界条件的Poisson方程。对于工程计算常用的粗糙网络,作者建议采用α≲0.5的预示校正差分格式。