求解浅水波方程的熵相容格式

提出了一种求解浅水波方程组的熵相容格式．在熵稳定通量中添加特征速度差分绝对值的项来抵消解在跨过激波时所产生的熵增,从而实现熵相容．新的数值差分格式具有形式简单、计算效率高、无需添加任何的人工数值粘性的特点．数值算例充分说明了其显著的优点．利用新格式成功地模拟了不同类型溃坝问题的激波、稀疏波传播及溃坝两侧旋涡的形成,是求解浅水波方程组较为理想的方法.     

不可压磁流体力学方程组的高精度数值解法

为数值求解低雷诺数下不可压流体在电磁场作用下的流动,提出一种四阶紧致差分方法.由二维原始变量的MHD方程组出发,推导出具有较少未知量的电流密度-涡量-流函数形式MHD方程组.建立了求解二维非定常不可压MHD方程组的电流密度-涡量-流函数形式的四阶精度紧致差分格式.为验证本文提出的高精度紧致差分方法的精确性和可靠性,对有...     

迎风紧致格式与驱动方腔流动问题的直接数值模拟

本文给出了一种求解不可压缩流动问题的高精度差分格式,即迎风紧致格式.出发方程采用二维非定常原始变量Naiver-Stokes方程组.在差分方程中,对流项采用三阶精度的迎风紧致差分,其余空间导数项采用四阶紧致差分.本文利用该差分格式在等距网格上数值模拟了驱动方腔流动中的分离涡运动.在257×257的细网格上,Re数最高计算到10000.Re≤5000时的计算结果与前人结果符合得很好.当Re≥7500时发现流动不存在定常层流解而为非定常周期性解,并首次给出了非定常解的结果。     

太阳风中的Alfvén孤波和中间激波 *

提出中间激波由动力Alfv啨ｎ波经波形连续变陡演化而成 .由动力Alfv啨ｎ波的色散作用形成Alfv啨ｎ孤波再进而演化成为激波 .由导出的非线性自治方程组计算出各种稳态结构 .它们分别对应中间激波(I_b 类) ,Switch-on激波和Alfvén孤波 .另有一类由稀疏孤波演化成的稀疏型间断面 .由于方程中没有强耗散机制 ,所得到的激波峰面上残留着孤波系列 .     

一种治愈强激波数值不稳定性的混合方法

HLLC（Harten-Lax-Leer-contact）格式是一种高分辨率格式,能够准确捕捉激波、接触间断和稀疏波.但是使用HLLC格式计算多维问题时,在强激波附近会出现激波不稳定现象.FORCE（first-order centred）格式在强激波附近表现出很好的稳定性,并且其数值耗散比HLL（Harten-Lax-Leer）格式小.分析了HLLC格式和FORCE格式在特定流动条件下的稳定性,构造了HLLC-FORCE混合格式并且进一步结合开关函数来消除HLLC格式的激波不稳定现象.数值试验表明新构造的混合格式不仅能够消除HLLC格式的激波不稳定现象,还最大程度地保留HLLC格式高分辨率的优点.     

一维非线性弹性力学方程组的Riemann问题

本文研究了一维非线性弹性力学方程组的Riemann问题.根据左右状态所处的相对位置,分情况构造了问题的唯一整体解.由于激波条件退化,系统的基本波除了稀疏波和激波还包含退化激波.     

多等级交通流LWR模型中的非线性波描述与WENO数值逼近

证明了交通流多等级LWR(Lighthill-Whitham-Richards)模型的双曲性质,并根据交通流的特征给出关于其非线性波的描述,主要包括车流通过激波和稀疏波时密度和速度的单调性变化．由于方程组没有显式的特征分解,所以引入具有高分辨和高精度的WENO(weighted essentially non~oscillatory)格式作数值模拟,得到与理论描述完全一致的数值结果．     

Riemann问题的间断有限元

用时空全离散间断零次有限元对Riemarm问题进行了数值求解，没有出现振荡，很好的模拟了稀疏波的逐渐稀疏化和激波的剧烈变化。     

带有摩擦项的广义Chaplygin气体非对称Keyfitz-Kranzer方程组的Riemann问题

研究了带有摩擦项的广义Chaplygin气体非对称Keyfitz-Kranzer方程组的Riemann问题,并得到其Riemann解的整体结构.Riemann解中包含激波,稀疏波,接触间断和δ-激波.与齐次非对称Keyfitz-Kranzer方程组不同的是非齐次非对称Keyfitz-Kranzer方程组的Riemann解是非自相似的.     

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

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

超音速探测器-刚性盘-缝-带型降落伞系统的大涡模拟研究

研究了Mach数为2时,流场不同块结构自适应网格加密精度对探测器-刚性盘-缝-带型降落伞系统的气动减速性能以及流场结构特性的影响.对于非定常可压缩流体流动,采用了兼顾激波与湍流的WENO(weighted essentially non-oscillatory)和TCD(tuned center difference)混合计算格式以及拉伸涡亚格子模型的大涡模拟方法.结果表明:在较低的流场块结构自适应网格分辨率下,是难以准确模拟计算降落伞系统重要的气动阻力系数和捕捉流场流动特征细节的.随后验证了流场自适应网格的收敛性.     

Shock boundary layer interactions in dense gases

The concept of triple deck theory is applied to study laminar interacting boundary layers of dense gases in external purely supersonic flow. An impinging shock is generated at distances which are large compared to the upper deck thickness. As predicted by weakly nonlinear theory such a discontinuity may disintegrate into a sonic shock and an associated wave fan depending on its amplitude and the magnitude of the so called fundamental derivative. Incoming and outgoing waves are computed analytically by means of the method of multiple scales taking into account that mutual interaction effects between them are restricted to the linear upper deck region. The lower deck problem is solved numerically. The results show that it is possible to reduce the size of the separation bubble or even to avoid the occurrence of flow separation by choosing an optimal thermodynamic state. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Riemann problem for the pressure-gradient system in three pieces

The Riemann problem for a two-dimensional pressure-gradient system is considered. The initial data are three constants in three fan domains forming different angles. Under the assumption that only a rarefaction wave, shock wave or contact discontinuity connects two neighboring constant initial states, it is proved that the cases involving three rarefaction waves are impossible. For the cases involving one shock (rarefaction) wave and two rarefaction (shock) waves, only the combinations when the three elementary waves have the same sign are possible (impossible).     

On the Riemann problem for 2D compressible Euler equations in three pieces

The Riemann problem for two-dimensional isentropic Euler equations is considered. The initial data are three constants in three fan domains forming different angles. Under the assumption that only a rarefaction wave, shock wave or contact discontinuity connects two neighboring constant initial states, it is proved that the cases involving three shock or rarefaction waves are impossible. For the cases involving one rarefaction (shock) wave and two shock (rarefaction) waves, only the combinations when the three elementary waves have the same sign are possible (impossible).     

Fine structures for the solutions of the two-dimensional Riemann problems by high-order WENO schemes

The two-dimensional Riemann problem with polytropic gas is considered. By a restriction on the constant states of each quadrant of the computational domain such that there is only one planar centered wave connecting two adjacent quadrants, there are nineteen genuinely different initial configurations of the problem. The configurations are numerically simulated on a fine grid and compared by the 5th-order WENO-Z5, 6th-order WENO-??6, and 7th-order WENO-Z7 schemes. The solutions are very well approximated with high resolution of waves interactions phenomena and different types of Mach shock reflections. Kelvin-Helmholtz instability-like secondary-scaled vortices along contact continuities are well resolved and visualized. Numerical solutions show that WENO-??6 outperforms the comparing WENO-Z5 and WENO-Z7 in terms of shock capturing and small-scaled vortices resolution. A catalog of the numerical solutions of all nineteen configurations obtained from the WENO-??6 scheme is listed. Thanks to their excellent resolution and sharp shock capturing, the numerical solutions presented in this work can be served as reference solutions for both future numerical and theoretical analyses of the 2D Riemann problem.     

一维不定常流方程的Harten解

为了利用一个二阶五点差分格式,对一维不定常流方程,计算强平面爆炸波遇不动固壁正反射的解,本文采用一定的技巧构造了边界和次边界的差分格式,并推算了初值的自模拟解析解.同时对初值的奇性提出了处理方法,从而满意地得到了这个较为困难问题的计算结果.本文提出和所利用的差分格式在处理激波间断问题方面具有实用意义.     

A composite semi-conservative scheme for hyperbolic conservation laws

In this work a first order accurate semi-conservative composite scheme is presented for hyperbolic conservation laws. The idea is to consider the non-conservative form of conservation law and utilize the explicit wave propagation direction to construct semi-conservative upwind scheme. This method captures the shock waves exactly with less numerical dissipation but generates unphysical rarefaction shocks in case of expansion waves with sonic points. It shows less dissipative nature of constructed scheme. In order to overcome it, we use the strategy of composite schemes. A very simple criteria based on wave speed direction is given to decide the iterations. The proposed method is applied to a variety of test problems and numerical results show accurate shock capturing and higher resolution for rarefaction fan.     

可压缩平面流中尖角处不定常涡产生的初始阶段的数值模拟

本文采用有限体积法结合高分辨率的ＴＶＤ（总变差不增）格式，对平面激波遇矩形、三角形障碍物问题进行数值模拟，采用Ｓｃｈｗａｒｚ变换生成计算网格，可使网格线在尖角处较密集且网格设置符合流场物理特性。初始阶段的流场复杂且粘性影响的尺度很小，模型中可以忽略流体的粘性。计算结果描述了平面无粘可压流中，由于弓形激波后熵与温度的不均匀性引起的尖角下游不定常集中涡的产生过程。数值计算与应我们要求所做激波管实验的光测结果相符合。     

A difference scheme with anti-diffusion terms to improve the computation of contact discontinuities and the study of weak discontinuities in compressible flow

The investigation of Mach reflection formed after the impingement of a weak plane shock wave on a wedge with shock Mach number M_s near 1, is still an open problem[12]. It's difficult for shock tube experiments with interferometer to detect contact discontinuities if it is too weak; also difficult to catch with due accuracy the transition condition between Mach reflection and regular reflection. The interest to this phenomenon is continuing, especially for weak shocks, because there was systematic discrepancy between simplified three shock theory of von Neumann [8] and shock tube results [15] which was named by G. Birkhoff as “von Neumann Paradox on three shock theory” [18].In 1972, K.O.Friedrichs called for more computational efforts on this problem. Recently it is known that for weak impinging shocks it's still difficult to get contact discontinuities and curved Mach stem with satisfactory accuracy. Recent numerical computation sometimes even fails to show reflected shock wave[6]. These explain why von Neumann paradox of the three shock theory in case of weak discontinuities is still a problem of interesting [9,12,14]. In this paper, on one hand, we investigate the numerical methods for Euler's equation for compressible inviscid flow, aiming at improving the computation of contact discontinuities, on the other hand, a methodology is suggested to correctly plot flow data from the massive information in storage. On this basis, all the reflected shock wave , contact discontinuities and the curved Mach stem are determined. We get Mach reflection under the condition when over-simplified shock theory predicts no such configuration[5].     

On the structure of internal dissipation of composite compact schemes for gasdynamic simulation

The structure of the internal dissipation terms in composite compact schemes intended for gasdynamic simulation is considered. The main cause of the insufficient stability of high-order accurate schemes is indicated. A method for controlling the dissipative properties of schemes is proposed that makes it possible to compute compressible gas flows with strong shock waves. The supersonic turbulent unsteady flow past a two-dimensional cavity directed toward the stream is computed.