Decoupling and recoupling of detonation waves associated with sudden expansion

Detonation propagation behavior associated with sudden expansions has been investigated both experimentally and numerically. Different mechanisms, from sustained propagation to detonation failure and reinitiation including shock and flame front decoupling and recoupling have been observed with the schlieren technique. The shock-induced flame propagation has been modeled with two-step chemistry and structured two-dimensional CFD on arbitrary geometries. The results of the numerical simulations show good correspondence to the variety of phenomena observed in experiments. Thus the numerical simulation can be used to study detonation propagation in complex geometries. It provides a tool for the design of safety devices and aids experimental investigations. Received 4 August 1995 / Accepted 25 April 1996     

Numerical study of wave processes in a pressure-wave refrigerator

The paper provides some results concerning the numerical study of the strongly transient gasdynamic processes in a pressure-wave refrigerator (PWR). A hierarchical set of numerical models from the simplest one-dimensional to a fully three-dimensional formulation is introduced. The computations show that one-dimensional solutions give a reasonable foundation for the understanding of PWR operational principles but cannot satisfactorily predict the refrigeration efficiency. Good agreement with experiment is achieved by considering two- and three-dimensional effects of gas mixing by overlap of rotating nozzles and expansion tubes. Received Received 14 February 1996 / Accepted 1 June 1996     

CH4-O2 混合气中爆燃爆震转捩的数值模拟

运用化学流体力学基本理论和两步燃烧反应模型原理，建立了一维封闭体系可燃气爆燃爆震转变现象的数学模型，利用拉格朗日质量坐标变换下的Ｌａｘ－Ｗｅｎｄｒｏｆ和Ｍｅｃｏｒｍｉｃ气动差分与Ａｄａｍｓ化学差分格式，求解基本方程，成功地完成了过程的数值模拟，清楚地说明了可燃气中ＤＤＴ现象由压缩波到激波达到稳态爆震的发生机制和火焰带引生爆震波的过程行为。     

裂尖超钝化:实验、理论与数值模拟

本文报道了在改性聚丙烯裂尖平面应变段发生超钝化的实验发现。利用环向冷拉模型,可说明裂尖超钝化的形成取决于冷拉传播速率与加载率之比。裂尖超钝化可使材料韧性,尤其是冲击韧性的改善,达到一个新的水平。本文在超弹性-粘塑性各向异性损伤的有限变形本构模型下进行了详尽的数值计算。模拟结果成功地再现了裂尖超钝化的基本特征。     

Lyapunov stable penalty methods for imposing holonomic constraints in multibody system dynamics

This paper presents stability and convergence results on a novel approach for imposing holonomic constraints for a class of multibody system dynamics. As opposed to some recent techniques that employ a penalty functional to approximate the Lagrange multipliers, the method herein defines a penalized dynamical system using penalty-augmented kinetic and potential energies, as well as a penalty dependent constraint violation dissipation function. In as much as the governing equations are not typically cocreive, the usual convergence criteria for linear variational boundary value problems are not directly applicable. Still numerical simulations by various researchers suggest that the method is convergent and stable. Despite the fact that the governing equations are nonlinear, the theoretical convergence of the formulation is guaranteed if the multibody system is natural and conservative. Likewise, stability and asymptotic stability results for the penalty formulation are derived from well-known stability results available from classical mechanics. Unfortunately, the convergence theorem is not directly applicable to dissipative multibody systems, such as those encountered in control applications. However, it is shown that the approximate solutions of a typical dissipative system converge to a nearby collection of trajectories that can be characterized precisely using a Lyapunov/Invariance Principle analysis. In short, the approach has many advantages as an alternative to other computational techniques:

Possibility of acceleration of the threshold processes for multi-component gas in the front of a shock wave

Studies of translational nonequilibrium in the front of a shock wave propagating in a three-component gas were performed by the Monte Carlo simulation method. Simulations were performed for mixtures of components with molecular mass ratios , and shock Mach number . The distribution of relative velocities for pairs of molecules of heavy low-concentration additives 2 and 3 substantially exceeded, in the front, its equilibrium values behind the wave at high values of . The maximum value of this superequilibrium was about for the numerical density ratio: 1000:1:1 and . Calculations showed that high values of the effect of superequilibrium take place up to a ratio of densities 200:1:1. Simulations performed for and a mixture of He, molecular oxygen and Xe with the numerical density ratio 200:1:1 showed also the high value of the superequilibrium effect at corresponding to dissociation threshold of oxygen. Thus, dissociation of oxygen by collisions with Xe in the front of a wave may have a considerably higher rate than total dissociation behind the wave. Received 4 August 1995 / Accepted 25 April 1996     

Real-Space Renormalization Estimates for Two-Phase Flow in Porous Media

We present a spatial renormalization group algorithm to handle immiscibletwo-phase flow in heterogeneous porous media. We call this algorithmFRACTAM-R, where FRACTAM is an acronym for Fast Renormalization Algorithmfor Correlated Transport in Anisotropic Media, and the R stands for relativepermeability. Originally, FRACTAM was an approximate iterative process thatreplaces the L × L lattice of grid blocks, representing the reservoir,by a (L/2) × (L/2) one. In fact, FRACTAM replaces the original L× L lattice by a hierarchical (fractal) lattice, in such a way thatfinding the solution of the two-phase flow equations becomes trivial. Thistriviality translates in practice into computer efficiency. For N=L ×L grid blocks we find that the computer time necessary to calculatefractional flow F(t) and pressure P(t) as a function of time scales as N^1.7 for FRACTAM-R. This should be contrasted with thecomputational time of a conventional grid simulator N^2.3. The solution we find in this way is an accurateapproximation to the direct solution of the original problem.     

纳米压痕过程的三维有限元数值试验研究

采用有限元方法模拟了纳米压痕仪的加、卸载过程，三维有限元模型考虑了纳米压痕仪的标准Berkovich压头．介绍了有限元模型的几何参数、边界条件、材料特性与加载方式，讨论了摩擦、滑动机制、试件模型的大小对计算结果的影响，进行了计算结果与标准试样实验结果的比较，证实了模拟的可靠性．在此基础上，重点研究了压头尖端曲率半径对纳米压痕实验数据的影响．对比分析了尖端曲率半径r＝0与r＝100nm两种压头的材料压痕载荷—位移曲线．结果表明，当压头尖端曲率半径r≠0时，基于经典的均匀连续介质力学本构理论、传统的实验手段与数据处理方法，压痕硬度值会随着压痕深度的减小而升高．     

船用气水分离器惯性级流场分析及阻力特性研究

采用二阶全展开ETG有限元与大涡模拟（LES）相结合的算法，对120－20－35－3型船用气水分离器惯性级在不同雷诺数条下的流动进行了模拟，通过其中一组雷诺数条件下计算所得该实验件阻力系数与物理实验所得阻力系数相比较，确定出该雷诺数条件下采用大涡模拟时所需的亚格子应力模型常数，将该常数带入其它各组雷诺数条件下的计算中，并将计算结果与相同条件下的物理实验结果相比较，证实了该常数的通用性。该常数一经确定，对各雷诺数条件下的流场进行分析，结果反映出采用二阶全展开ETG有限元与大涡模拟（LES）相结合的算法可以捕捉到非常丰富的涡系及涡动的时变过程。在流场分析的基础上本文计算了该实验件内的能耗场，计算结果表明实验件内的能耗主要集中在大涡丰富的区段内。     

泄爆过程中二次爆炸的动力学机理研究

在容积为0.00814m3的柱形泄爆容器中，对泄爆现象进行实验研究. 容器内充满当量比为1的甲烷-空气预混气，采用底端中心点火，泄爆压力为230±15kPa. 基于k-ε湍流模型和EBU燃烧模型，利用同位网格的SIMPLE算法，对该现象进行了数值模拟. 实验和计算获得的外轴线上4个测压点的压力曲线和外流场的阴影和数值照片，形象地描述了高压泄爆时外部流场的变化. 数值结果与实验结果基本一致. 根据实验和数值结果，详细地讨论了泄爆过程中二次爆炸产生的动力学机理. 泄爆的初始阶段，在破膜激波的引导下，泄出的未燃气体因欠膨胀在外流场形成稀疏波低压区和悬激波高压区. 高压区可燃气体密度和温度上升，成为高密度的预热区域. 随后，火焰以射流形式从泄爆口泄出，点燃可燃气云. 受湍流等因素的影响，特别在高密度的预热区域，燃烧速率可能迅速增大，从而导致二次爆炸.