龙格库塔间断有限元方法在计算爆轰问题中的应用

构造求解带源项守恒律方程组的龙格库塔间断有限元(RKDG)方法,并分别结合源项的Strang分裂法和无分裂法数值求解模型守恒律方程和反应欧拉方程.为了和有限体积型WENO方法进行比较,设计计算源项的WENO重构格式.对一维带源项守恒律的计算表明,对于非刚性问题,RKDG方法比有限体积型WENO方法的误差更小;对于刚性问题,RKDG方法对于间断面位置的捕捉更为精确.对于一二维爆轰波问题的计算结果表明,RKDG方法对爆轰波结构的分辨和爆轰波位置的捕捉能力更强.     

一种基于HLLC算法的Mie-Grüneisen多介质混合模型

建立一种修正的HLLC （Harten-Lax-Van Leer-Contact）格式下稳定的Mie-Grüneisen多介质混合计算模型.Mie-Grüneisen混合模型中的通量包括守恒和非守恒两个部分,原始的HLLC格式对守恒部分适用,但是原始的HLLC格式直接用于非守恒部分,很难控制住数值振荡的产生.在原始格式中,间断面的移动速度为间断网格的左侧或右侧速度,修正后替换为网格中的平均速度,经过修正后,对HLLC格式重新进行推导,并随之扩展到二维问题.数值实验表明,利用修改后的HLLC格式,Mie-Grüneisen混合模型可以取得较好的稳定性和准确性.     

带副翼的翼身组合体绕流的Euler和N-S方程解

将对接分区网格与分区求解算法结合,有效地求解了带副翼偏转的翼身组合体绕流的N S方程.数值方法中选用VanLeer分裂格式离散无粘通量项,采用中心差分格式来离散粘性通量项.分区交界面采用了一种满足通量守恒的内边界耦合条件.数值算例表明该方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

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

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

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

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

二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法

构造矩形网格下求解Lagrangian坐标系下气动方程组的单元中心型格式. 空间离散采用控制体积间断Petrov-Galerkin方法,时间离散采用二阶TVD Runge-Kutta方法. 利用限制器来抑制非物理震荡并保证RKCV算法的稳定性. 构造的算法可以保证物理量的局部守恒. 与Runge-Kutta间断Galerkin（RKDG）方法相比较,RKCV方法的计算公式少一项积分项使得计算较简单. 给出一些数值算例验证了算法的可靠性及效率.     

双曲守恒律的Taylor-Galerkin有限元方法

利用双曲守恒律的Hamilton-Jacobi方程形式,应用Taylor公式与Galerkin有限元给出了求解双曲守恒律的计算方法。采用TVD差分格式的构造思想,对数值通量作修正,在等距网格情形下有限元方法得到的计算格式满足TVD性质,并给出了数值例子。     

多介质流的高分辨率Euler方法

在多介质流动问题中,不同介质有不同的状态方程。这使通量成为间断函数,从而没有通量的Jacobi矩阵。而用Euler坐标系描述的方程组的很多高分辨率格式都要用到Jacobi矩阵及其特征值和特征向量,即要求通量连续可微。因此必须重新处理整个守恒律方程组。对于γ气体问题将γ看作一个新未知量并增加一个守恒方程,从而使整个方程组的通量成为光滑函数,为高分辨率格式的构造铺平了道路。由于真实流动只遵守三个守恒律,多加的一个守恒律虽然对偏微分方程组没有影响,但对差分方程数值解有影响。这一点在数值实验中已有表现。提出了一个方案将这一影响尽量消除。所用格式可完全照搬单介质流动的任何现有格式。对一维多介质流动Euler方程组的激波管问题的数值实验表明这样处理所构造的格式具有同单介质流动问题同样的效果。     

二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法（英文）

构造矩形网格下求解Lagrangian坐标系下气动方程组的单元中心型格式.空间离散采用控制体积间断Petrov-Galerkin方法,时间离散采用二阶TVD Runge-Kutta方法.利用限制器来抑制非物理震荡并保证RKCV算法的稳定性.构造的算法可以保证物理量的局部守恒.与Runge-Kutta间断Galerkin(RKDG)方法相比较,RKCV方法的计算公式少一项积分项使得计算较简单.给出一些数值算例验证了算法的可靠性及效率.     

可压缩多介质流体数值模拟中的Level-Set间断跟踪方法

针对可压缩多介质流体的数值模拟,发展了一种Level-Set间断追踪技术,用LS(Level-Set)函数追踪激波和捕捉界面,用Riemann问题解构造带状区域内的虚拟流体状态,对物理量的外推方法、间断附近虚拟流体的构造、间断推进速度的计算等问题进行了研究.最后对可压缩多介质流体一维和二维守恒律方程组进行数值模拟,数值计算采用通量重构的高精度WENO格式,计算结果令人满意.     

A gas kinetic scheme for the Baer–Nunziato two-phase flow model

Numerical methods for the Baer–Nunziato (BN) two-phase flow model have attracted much attention in recent years. In this paper, we present a new gas kinetic scheme for the BN two-phase flow model containing non-conservative terms in the framework of finite volume method. In the view of microscopic aspect, a generalized Bhatnagar–Gross–Krook (BGK) model which matches with the BN model is constructed. Based on the integral solution of the generalized BGK model, we construct the distribution functions at the cell interface. Then numerical fluxes can be obtained by taking moments of the distribution functions, and non-conservative terms are explicitly introduced into the construction of numerical fluxes. In this method, not only the complex iterative process of exact solutions is avoided, but also the non-conservative terms included in the equation can be handled well.     

A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows

We propose a new model and a solution method for two-phase compressible flows. The model involves six equations obtained from conservation principles applied to each phase, completed by a seventh equation for the evolution of the volume fraction. This equation is necessary to close the overall system. The model is valid for fluid mixtures, as well as for pure fluids. The system of partial differential equations is hyperbolic. Hyperbolicity is obtained because each phase is considered to be compressible. Two difficulties arise for the solution: one of the equations is written in non-conservative form; non-conservative terms exist in the momentum and energy equations. We propose robust and accurate discretisation of these terms. The method solves the same system at each mesh point with the same algorithm. It allows the simulation of interface problems between pure fluids as well as multiphase mixtures. Several test cases where fluids have compressible behavior are shown as well as some other test problems where one of the phases is incompressible. The method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state.     

A fully discrete,kinetic energy consistent finite-volume scheme for compressible flows

A robust, implicit, low-dissipation method suitable for LES/DNS of compressible turbulent flows is discussed. The scheme is designed such that the discrete flux of kinetic energy and its rate of change are consistent with those predicted by the momentum and continuity equations. The resulting spatial fluxes are similar to those derived using the so-called skew-symmetric formulation of the convective terms. Enforcing consistency for the time derivative results in a novel density weighted Crank–Nicolson type scheme. The method is stable without the addition of any explicit dissipation terms at very high Reynolds numbers for flows without shocks. Shock capturing is achieved by switching on a dissipative flux term which tends to zero in smooth regions of the flow. Numerical examples include a one-dimensional shock tube problem, the Taylor–Green problem, simulations of isotropic turbulence, hypersonic flow over a double-cone geometry, and compressible turbulent channel flow.     

An all-speed relaxation scheme for interface flows with surface tension

We consider interface flows where compressibility and capillary forces (surface tension) are significant. These flows are described by a non-conservative, unconditionally hyperbolic multiphase model. The numerical approximation is based on finite-volume method for unstructured grids. At the discrete level, the surface tension is approximated by a volume force (CSF formulation). The interface physical properties are recovered by designing an appropriate linearized Riemann solver (Relaxation scheme) that prevents spurious oscillations near material interfaces. For low-speed flows, a preconditioning linearization is proposed and the low Mach asymptotic is formally recovered. Numerical computations, for a bubble equilibrium, converge to the required Laplace law and the dynamic of a drop, falling under gravity, is in agreement with experimental observations.     

A roe-average algorithm for a granular-gas model with non-conservative terms

A Roe-average algorithm has been derived for a granular-gas model, proposed by Goldshtein and Shapiro [Goldshtein, Shapiro, Mechanics of collisional motion of granular materials: Part 1. General hydrodynamic equations, J. Fluid Mech. 282 (1995) 75–114], which contains non-conservative terms in the Euler-like hyperbolic governing equations apart from sink terms, which arise from inelastic collision of granules and are present only in the energy equation. The non-conservative terms introduce non-isentropic effects in acoustic-wave propagation within granular media and they also contribute to the Rankine–Hugoniot relations across a discontinuity. A Roe-average algorithm, based on the same granular-gas model, was derived in the literature [V. Kamenetsky, A. Goldshtein, M. Shapiro, D. Degani, Evolution of a shock wave in a granular gas, Phys. Fluids, 12 (2000) 3036–3049] which then required the implementation of a shock-fitting technique at a discontinuity. In the present work, Roe-averaged variables have been obtained from the Rankine–Hugoniot jump relations and the non-conservative terms have been incorporated in the numerical flux formula consistent with upwind principles associated with the granular speed of sound. Results for unsteady one-dimensional granular flows, colliding with a wall, demonstrate the capability of the proposed algorithm to capture strong shocks in addition to flow features not found in molecular gases, such as a fluidized region downstream of the shock and a compacted solid-block region adjacent to the wall.     

连续能量蒙特卡罗方法组件均匀化中的临界效应

研究蒙特卡罗均匀化中的临界效应,应用B_N理论修正临界效应.简化压水堆模型和临界试验装置的验证表明:在连续能量蒙特卡罗方法均匀化中,B_N方法很好的考虑了群常数的临界效应,曲率基模成功度量了堆芯的有限性.     

Vibration and stability of a free circular plate subjected to non-conservative loading

This paper is concerned with the stability and vibration of a completely free circular plate subjected to a non-conservative edge loading. The eigencurves and mode shapes of the plate are obtained for various values of the non-conservativeness parameter. Numerical results are presented for the asymmetrical mode shapes of the plate. Interesting conclusions are drawn from these results some of which are verified analytically.     

Power flow and energy distribution in conservatively or non-conservatively coupled oscillators under correlative excitations

As a foundation of SEA for consevatively or non-conservativelycoupled systems under correlative excitations,power flow and energy distribu-tion in two coupled osillators under proportionally correlative excitations arestudied mainly in this paper.Theoretical relationship among power flow andtime averaged vibational energy of the osillators,as well as power balance equa-tions and energy ratio or them,are derived.Further along,power flow,energydistribution and power balance equations of the systems in conservative ornon-conservative coupling and under proportional correlative excitations arediscussed in detail,with comparing with characteristics of those innon-correlative excitations.Numerical resu1ts supporting the basic theory aregiven,too.Results show that the power flow and energy distribution of theosillators under correlative excitations,either conservative or non-conservativecoupling,are much different to those in the condition of non-correlativeexcitations.     

LWDG Method for a Multi-Class Traffic Flow Model on an Inhomogeneous Highway

In this paper, we apply the discontinuous Galerkin method with Lax-Wendroff type time discretizations (LWDG) using the weighted essentially non-oscillatory (WENO) limiter to solve a multi-class traffic flow model for an inhomogeneous highway, which is a kind of hyperbolic conservation law with spatially varying fluxes. The numerical scheme is based on a modified equivalent system which is written as a "standard" hyperbolic conservation form. Numerical experiments for both the Riemann problem and the traffic signal control problem are presented to show the effectiveness of these methods.     

原子分子离化截面分波法的数值计算

讨论了原子分子离化截面分波方法的数值计算,特别是高分波项的处理问题。利用外推法,在包括任意大的高分波项的同时,节省了大量计算机时,而且得到较满意的计算结果.