叶栅激波流场的数值模拟与高分辨率差分格式

一、前言 近年来,激波流场数值模拟技术的一个重要进展,是各种高分辨率总变差减小(TVD)格式的建立和发展。这类格式与常用的中心差分方法相比,其突出的优点在于,不需要人为地在差分方程中附加耗散项就能给出激波上、下游无波动的数值解,保证了数值流场具有很高的空间分辨率。目前这种方法已引起人们的高度重视,并逐渐发展成为应用计算中的流行格式。本文的目的,是探讨TVD格式在叶栅激波流场数值分析     

具有声激波的跨音流管道中声传播的数值方法

本文采用四阶MacCormack格式和附加四阶粘性项方法,求解具有声激波的跨音流变差分格壁管中的声传播问题,比前人结果有明显改善。本文详细介绍了这种差分方法,特别是关于截面硬式,人工粘性项和计算可靠性判据。这种方法省内存,省机时,可以在微机上实现计算。     

双曲守恒律系统的分段边界层函数方法(PBLM)

本文提出了一个新的高阶Godunov格式。此格式放弃了[9]、[10]中关于参数在格子中满足多项式分布近似及在格子边界上存在间断的假设,直接引入了一个分段边界层型函数分布假设。由于引入的分布函数是单调可微的,因此PBLM格式无需进行如同MUSCL、PPM等格式中的单调性校核。该格式由于不进行单调性修正,在PPM格式中需进行修正而精度降阶的点上仍保持原有精度。对一维激波管的计算表明PBLM格式对激波的展开比PPM格式还要小,计算时间相当。同PPM格式一样,PBLM格式在激波后存在有2%~4%的皮后伪振荡,应加上适当的人工粘性。     

高分辨率熵相容格式及其在激波流动中的应用

为解决熵守恒格式在激波附近出现数值振荡的问题,本文将熵相容格式与MUSCL格式相结合,提出一种既能适合于激波问题、又不依赖于传统人工黏性经验模型的高分辨率熵相容格式,通过对多个激波问题的数值计算,并对比二阶中心格式、熵守恒格式、熵相容格式和高分辨率熵相容格式的计算结果,发现:熵相容格式具有较好的激波捕捉能力,有效解决了熵守恒格式在激波附近的数值振荡问题;MUSCL重构格式进一步提高了熵相容格式的数值模拟能力,既能精确捕捉激波附近的流动细节,又在光滑区保持二阶精度;在对比的四种格式中,本文提出的高分辨率熵相容格式对激波问题的预测性能最佳。该项工作对发展激波湍流相互作用模型、提高跨/超音速叶轮机械流动预测精度具有理论价值和应用潜力。     

再入体头部粘性流场高分辨率格式数值模拟

TVD格式是目前数值研究以激波为主要特征之一的超声速、高超声速流场的最先进的算法之一。本文用二阶迎风TVD格式,对三种烧蚀外形的轴对称粘性流场和10°钝锥有攻角三维粘性流场进行了数值模拟,得到了高质量的头部脱体激波和与实验结果及直线推进法计算一致的物面压力分布,表明了TVD格式在再入体粘性绕流计算中的独特优势。     

迎风TVD格式在粘性流计算中的应用研究与改进

首先以高超声速钝楔绕汉热流计算和Ｍ∞＝２．０、激波角θ＝３２．５８°的平板激波附面层干扰分离流动计算,研究了迎风ＴＶＤ格式的粘性特性、发现即使在很小的物面网格下,原有格式计算驻点热流仍较实验和文献结果低,激波附面层干扰的分离区小于实验值,原因主要在于熵修正公式。针对这些不足提出了新的熵修正公式,用改进后的格式重新计算,得到了与实验符合较好的结果。最后用改进的格式对分离流动的壁面温度控制效应进行了研     

粘性跨音速叶栅绕流的数值模拟及前驻点处理

一、前言 对于粘性跨音速叶栅绕流,直接用N-S方程数值求解仍有很多困难,因此利用边界层理论考虑流体的粘性效应,进行有粘—无粘迭代计算是粘性流体叶栅绕流的很好模拟。采用多重网格有限体积法进行跨音速叶栅的无粘流计算,具有易于适应复杂几何形状,保证差分格式的守恒性、能较准确地计算激波。加速收敛等良好特性。为了提高精度和简化计算,在求解边界层微分方程时,采用了Illingworth变换方法,用驻点方程的相似性解给定边值速度剖面,使用全湍流Prandtl混合长度粘性模式并且其混合长度选用Plctcher     

用TVD及MmB方法模拟理想非线性色谱过程

用高精度、高分辨率的激波捕捉法TVD(全变差衰减)及MmB(局部保最大、最小、有界)差分格式,对理想非线性色谱过程进行了数值模拟,以较高的分辨率确定了非线性色谱过程和激波效应,并计算了流出过程的保留时间,所得结果与理论分析及实验结果作了比较.     

拉氏数值模拟中的人工粘性方法

较好的人工粘性需要满足较小的计算开销、不能去除真实具有的涡运动等条件.提出一种应用于拉氏数值模拟中基于Lew人工粘性,同时增加了限制器的人工粘性方法.可以有效减少数值模拟结果对网格的依赖;采用特征值限制器控制施加的人工粘性大小,通过限制器能够区分激波压缩和等熵压缩;方便应用在二维、三维,结构网格或者非结构网格上.     

一种改进的HLLEM格式及其激波稳定性分析

使用低耗散激波捕捉格式对高超声速流动问题进行数值模拟时经常会遭受不同形式的激波不稳定性.本文基于二维无黏可压缩Euler方程,对低耗散HLLEM格式进行激波稳定性分析.结果表明:激波面横向通量中切向速度的扰动增长诱发了格式的不稳定性.通过增加耗散来治愈HLLEM格式的激波不稳定性.为了避免引入过多的耗散进而影响剪切层的分辨率,定义激波探测函数和亚声速区探测函数,使得只有在计算激波层亚声速区的横向数值通量时才增加耗散,其余地方的数值通量依然采用低耗散的HLLEM格式来计算.稳定性分析和数值模拟的结果表明,改进的HLLEM格式不仅保留了原格式高分辨率的优点,还大大提高了格式的鲁棒性,在计算强激波问题时能够有效地抑制不稳定现象的发生.     

High-order well-balanced schemes and applications to non-equilibrium flow

The appearance of the source terms in modeling non-equilibrium flow problems containing finite-rate chemistry or combustion poses additional numerical difficulties beyond that for solving non-reacting flows. A well-balanced scheme, which can preserve certain non-trivial steady state solutions exactly, may help minimize some of these difficulties. In this paper, a simple one-dimensional non-equilibrium model with one temperature is considered. We first describe a general strategy to design high-order well-balanced finite-difference schemes and then study the well-balanced properties of the high-order finite-difference weighted essentially non-oscillatory (WENO) scheme, modified balanced WENO schemes and various total variation diminishing (TVD) schemes. The advantages of using a well-balanced scheme in preserving steady states and in resolving small perturbations of such states will be shown. Numerical examples containing both smooth and discontinuous solutions are included to verify the improved accuracy, in addition to the well-balanced behavior.     

高精度TVD有限差分格式的自由梯度条件

1自由流条件守恒性气动方程组的高精度TVD有限差分格式一般通过使用通量限制器加入非线性的人工粘性来构成。在曲线坐标系中，通量限制器比较的两个因子与各自所在单元的界面雅可比转换行列式相关。为了消除畸变网格单元界面雅可比转换行列式的变化对通量限制器比较结果的影响，本文提出了高精度TVD格式的自由梯度性质，并分别就Osher和HartenTVD格式讨论了同时满足自由流和自由梯度的度量张量计算方法，使高精度TVD格式在畸变网格上得以实现。任意曲线坐标系中守恒型二维欧拉方程表示为其中，E、F为笛卡尔坐标系中的通量；J为坐标变…     

An efficient finite volume method for electric field–space charge coupled problems

The goal of this paper is to introduce some recently developed finite volume schemes to enable numerical simulation of electric field–space charge coupled problems. The key features of this methodology are the possibility of handling problems with complex geometries and accurately capturing the charge density distribution. The total variation diminishing (TVD) scheme and the improved deferred correction (IDC) scheme are used to compute the convective and diffusive fluxes respectively. Our technique is firstly verified with the computation of hydrostatic solutions in a two coaxial cylinders configuration. The homogeneous and autonomous injection from the inner or outer electrode is considered. Comparison has been made with the analytical solution. The numerical technique is also applied to the problem of corona discharge in a blade-plane configuration. The good agreement between our numerical solution and the one obtained with a combination approach of Finite Element Method (FEM) and Method of Characteristics (MoC) is shown.     

A comparison of spatial discretization schemes for differential solution methods of the radiative transfer equation

A comparison of discretization schemes required to evaluate the radiation intensity at the cell faces of a control volume in differential solution methods of the radiative transfer equation is presented. Several schemes developed using the normalized variable diagram and the total variation diminishing formalisms are compared along with essentially non-oscillatory schemes and genuinely multidimensional schemes. The calculations were carried out using the discrete ordinates method, but the analysis is equally valid for the finite-volume method. It is shown that the S schemes of the genuinely multidimensional family perform quite well, particularly in problems with discontinuous radiation intensity fields. However, they are time consuming, and so they do not always become more attractive regarding the trade-off between accuracy and computational requirements, in comparison with other high-order schemes that, although being less accurate, are also more economical.     

Application of the Generalized Differential Quadrature Method in Solving Burgers'Equations

The aim of this paper is to obtain numerical solutions of the one-dimensional, two-dimensional and coupled Burgers' equations through the generalized differential quadrature method (GDQM). The polynomial-based differential quadrature (PDQ) method is employed and the obtained system of ordinary differential equations is solved via the total variation diminishing Runge-Kutta (TVD-RK) method. The numerical solutions are satisfactorily coincident with
the exact solutions. The method can compete against the methods applied in the literature.     

抗爆容器内爆炸流场数值模拟

 采用计算流体动力学中的二阶精度TVD差分格式和特殊算子分裂法,按轴对称问题,对半球顶圆柱筒密闭式抗爆容器内部爆炸流场进行了数值模拟。计算得到的容器壁面载荷分布与实验结果基本一致。不同时刻爆炸流场压力分布图像清晰地描述了容器壁面的冲击波加载过程。     

Explicit Multi-Symplectic Splitting Methods for the Nonlinear Dirac Equation

In this paper, we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac (NLD) equation. Based on its multi-symplectic formulation, the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system. Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize the linear subproblem, respectively. And the nonlinear subsystem is solved by a symplectic scheme. Finally, a composition method is applied to obtain the final schemes for the NLD equation. We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly. Numerical experiments are presented to show the effectiveness of the proposed methods.     

A highly accurate method to solve Fisher’s equation

In this study, we present a new and very accurate numerical method to approximate the Fisher’s-type equations. Firstly, the spatial derivative in the proposed equation is approximated by a sixth-order compact finite difference (CFD6) scheme. Secondly, we solve the obtained system of differential equations using a third-order total variation diminishing Runge–Kutta (TVD-RK3) scheme. Numerical examples are given to illustrate the efficiency of the proposed method.     

Several teleportation schemes of an arbitrary unknown multi-particle state via different quantum channels

We first provide four new schemes for two-party quantum teleportation of an arbitrary unknown multi-particle state by using three-, four-, and five-particle states as the quantum channel, respectively. The successful probability and fidelity of the four schemes reach 1. In the first two schemes, the receiver can only apply one of the unitary transformations to reconstruct the original state, making it easier for these two schemes to be directly realized. In the third and fourth schemes, the sender can preform Bell-state measurements instead of multipartite entanglement measurements of the existing similar schemes, which makes real experiments more suitable. It is found that the last three schemes may become tripartite controlled teleportation schemes of teleporting an arbitrary multi-particle state after a simple modification. Finally, we present a new scheme for three-party sharing an arbitrary unknown multi-particle state. In this scheme, the sender first shares three three-particle GHZ states with two agents. After setting up the secure quantum channel, an arbitrary unknown multi-particle state can be perfectly teleported if the sender performs three Bell-state measurements, and either of two receivers operates an appropriate unitary transformation to obtain the original state with the help of other receiver's three single-particle measurements. The successful probability and fidelity of this scheme also reach 1. It is demonstrated that this scheme can be generalized easily to the case of sharing an arbitrary unknown multi-particle state among several agents.     

Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations

In this paper, we propose a finite difference AMR-WENO method for hyperbolic conservation laws. The proposed method combines the adaptive mesh refinement (AMR) framework  and  with the high order finite difference weighted essentially non-oscillatory (WENO) method in space and the total variation diminishing (TVD) Runge–Kutta (RK) method in time (WENO-RK)  and  by a high order coupling. Our goal is to realize mesh adaptivity in the AMR framework, while maintaining very high (higher than second) order accuracy of the WENO-RK method in the finite difference setting. The high order coupling of AMR and WENO-RK is accomplished by high order prolongation in both space (WENO interpolation) and time (Hermite interpolation) from coarse to fine grid solutions, and at ghost points. The resulting AMR-WENO method is accurate, robust and efficient, due to the mesh adaptivity and very high order spatial and temporal accuracy. We have experimented with both the third and the fifth order AMR-WENO schemes. We demonstrate the accuracy of the proposed scheme using smooth test problems, and their quality and efficiency using several 1D and 2D nonlinear hyperbolic problems with very challenging initial conditions. The AMR solutions are observed to perform as well as, and in some cases even better than, the corresponding uniform fine grid solutions. We conclude that there is significant improvement of the fifth order AMR-WENO over the third order one, not only in accuracy for smooth problems, but also in its ability in resolving complicated solution structures, due to the very low numerical diffusion of high order schemes. In our work, we found that it is difficult to design a robust AMR-WENO scheme that is both conservative and high order (higher than second order), due to the mass inconsistency of coarse and fine grid solutions at the initial stage in a finite difference scheme. Resolving these issues as well as conducting comprehensive evaluation of computational efficiency constitute our future work.