气相爆轰波基元反应模型数值模拟

应用频散可控格式和基元反应模型计算了氢气-氧气和氢气-空气混合物的爆轰波一维结构，并对爆轰波在二维轴对称突扩截面管内的传播过程进行了数值模拟。采用算子分裂的方法来处理化学反应的刚性问题。计算结果表明，频散可控格式可以精确捕捉爆轰波及其反射过程。爆轰波在突扩截面处绕射，由于稀疏波的作用而引起局部熄爆现象。数值结果还反应出，在一定条件下，熄爆区域可以二次起爆，而且二次起爆过程与可燃混合气体的敏感度有关，对于敏感度高的可燃气体，由于爆轰波具有较强的自持能力，其波阵面在绕射过程中所产生的局部熄爆区后会自动再次起爆

超声速气流中横向煤油射流的数值模拟

为了对超声速气流中液体横向射流进行数值模拟，采用Eulerian-Lagmngian方法描述气液两相流动现象，用随机轨道模型追踪液滴的运动历程，并考虑气流可压缩性、流场不均匀性及液滴变形对液滴运动的影响，开发了相应的计算程序。计算结果发现，液体射流与气流之间存在着强烈的相互作用，液滴在进入气流中不久就破碎成很小的子液滴，受煤油液雾的影响气流速度、温度急剧下降，同时在喷孔上游出现～道弓形激波。与液雾结构的纹影图像和液雾穿透的实验结果对比显示，数值计算结果基本合理。

Effects of the Spanwise Characteristic Length on the Transition Feature in the Wake of a Cylinder

The transition features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations of 3-D incompressible Navier-Stokes equations. The spanwise characteris-tic length determines the transition features and global properties of the wake.

Simulations of oil spread on a water surface

Numerical simulation of an oil slick spreading on still and wavy surfaces is described in this paper. The so-called sigma transformation is used to transform the time-varying physical domain into a fixed calculation domain for the water wave motions and, at the same time, the continuity equation is changed into an advection equation of wave elevation. This evolution equation is discretized by the forward time and central space scheme, and the momentum equations by the projection method. A damping zone is set up in front of the outlet boundary coupled with a Sommerfeld-Orlanski condition at that boundary to minimize the wave reflection. The equations for the oil slick are depth-averaged and coupled with the water motions when solving numerically. As examples, sinusoidal and solitary water waves, the oil spread on a smooth plane and on still and wavy water surfaces are calculated to examine the accuracy of simulating water waves by Navier-Stokes equations, the effect of damping zone on wave reflection and the precise structures of oil spread on waves.

Formation Mechanism of Cracks in Saturated Sand

The formation mechanism of “water film” (or crack) in saturated sand is analyzed theoretically and numerically. The theoretical analysis shows that there will be no stable “water film” in the saturated sand if the strength of the skeleton is zero and no positions are choked. It is shown by numerical simulation that stable water films initiate and grow if the choking state keeps unchanged once the fluid velocities decrease to zero in the liquefied sand column. The developments of “water film” based on the model presented in this paper are compared with experimental results.

Shock-Induced Near-Wall Two-Phase Flow Structure Over a Micron-Sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micronsized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed- phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail.The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.

Two bifurcation transitions of the floating half zone convection in a fat liquid bridge of larger Pr

The transient process of the thermocapillary convection was obtained for the large Pu floating half zone by using the method of three-dimensional and unsteady numerical simulation. The convection transits directly from steady and axisymmetric state to oscillatory flow for slender liquid bridge, and transits first from steady and axisymmetric convection to the steady and non-axisymmetric convection, then, secondly to the oscillatory convection for the fatter liquid bridge. This result implies that the volume of liquid bridge is not only a sensitive critical parameter for the onset of oscillation, but also relates to the new mechanism for the onset of instability in the floating half zone convection even in case of large Prandtl number fluid.

A new model of thermocapillary convection in floating half zone

A new model of thermocapillary convection in floating half zone was suggested in the present paper. The liquid bridge floats between two co-axis rods, the lower rod consists of metal with constant temperature and the upper rod consists of thermal insulating materials, where the normal gradient of temperature is nearly zero. In this case, the new model is relatively closer to simulate a half part of floating full zone in microgravity environment in comparison with the usual model of floating half zone. Basic features of the new model were studied by both numerical simulation and experiments, and the comparisons with the usual model were also discussed.

Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust

Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here the authors address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. $En Echelon$ crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality------from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that the authors do not understand earthquake mechanisms very well but that they still know little about our earth's crust.

加卸载响应比的数值模拟

用数值模型来模拟和验证加卸载响应比理论。编制非均匀各向异性材料宏细观损伤破坏程度，模拟固体材料在损伤破坏时的加卸载过程，证实了加卸载响应比理论的正确性。结果表明，ＬＵＲＲ这一反映地球介质临近破坏的指标，也可以用来预测材料的破坏。

Optimization of dynamic properties and link performance for directly-modulated gain-coupled DFB lasers

This paper demonstrates the respective roles that combined gain- and index-coupling play in the dynamic properties and overall link performance of DFB lasers. It is shown that for datacommunication applications, modest gain-coupling enables optimum transmission at 10Gbit/s.

Oscillatory Instabilities of Two-Layer Rayleigh-Marangoni-Benard Convection

The Rayleigh–Marangoni–Bénard convective instability (R–M–B instability) in the two-layer systems such as Silicone oil (10cSt)/Fluorinert (FC70) and Silicone oil (2cSt)/water liquids are studied. Both linear instability analysis and nonlinear instability analysis (2D numerical simulation) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. The results show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. The secondary instability phenomenon found in the real two-layer system of Silicone oil over water could explain the difference in the comparison of the Degen’s experimental observation with the previous linear stability analysis results of Renardy et al.

油水分离技术研究之一--直管和螺旋管的数值模拟

本文采用数值模拟的方法研究油水混合物在直管和螺旋管中的流动状况.计算采用Euler-Euler法和Euler模型:支配油水分离两相流动的基本方程包括连续方程和动量方程,湍流模型采用多相流中混合型k-ε模式,基本方程的离散和求解采用SIMPLE算法.利用Fluent软件,以直管和螺旋管为例进行了计算,获得了初步计算结果.计算表明,本文所用方法可以较好地模拟直管中油水在重力作用下的分离,以及螺旋分离器中油水在重力和离心力共同作用下的分离现象.并可为实验研究提供参考。

The inner Pressure of the Atoms in a Solid in Refined TFD Theory

The relation between the inner pressure of an atom in a solid and the density of energy of electrons under Refined TFD theory is given.

Two Bifurcation Transition Processes in Floating Half Zone Convection of Larger Prandtl Number Fluid

Processes of the onset oscillation in the thermocapillaxy convection under the Earth's gravity are investigated by the numerical simulation and experiments in a floating half zone of large Prandtl number with different volume ratio. Both computational and experimental results show that the steady and axisymmetric convection turns to the oscillatory convection of m=1 for the slender liquid bridge, and to the oscillatory convection before a steady and 3D asymmetric state for the case of a fat liquid bridge. It implies that, there are two critical Marangoni numbers related, respectively, to these two bifurcation transitions for the fat liquid bridge. The computational results agree with the results of ground-based experiments.