CFD study of a sudden-expanding coal combustor using Euler-Euler and Euler-Lagrange models

In this study, the Euler-Euler (E-E) and Euler-Lagrange (E-L) models designed for the same chemical mechanism of heterogeneous reactions were used to predict the performance of a typical sudden-expanding coal combustor. The results showed that the current E-E model underestimated the coal burnout rate because the particle temperature fluctuation on char combustion is not adequately considered. A comparison of the E-E and E-L simulations showed the underestimation of heterogeneous chemical reaction rates by the E-E model. (C) 2010 Elsevier Ltd. All rights reserved.

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.

激光诱发活塞的热负荷

活塞的热疲劳性能对柴油发动机的全寿命至关重要。由于能量有限和可控性差等缺点，现有实验系统均不能满意地进行活塞热负荷模拟实验。为此，提出并建立了一套激光诱发活塞热负荷的实验系统。该系统通过对激光束的空间整形，使之投射到活塞表面后诱发的温度场分布满足特定要求。基于PROFIBUS-DP现场总线技术实现了系统集成和实验过程的全反馈控制。针对活塞的典型热负荷条件，即高周热疲劳和热冲击分别进行实验，以模拟正常工作循环和“启动一停车”等热负荷或转速突变工况。通过设置加热一冷却周期或上限下限温度，可以获得相应的热负荷模拟实验效果。研究结果表明，采用经光束整形的激光进行活塞热负荷模拟实验具有周期短、可控性好等优点。 

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Step Structure in Cold-Rolled Deformed Nanocrystalline Nickel

The microstructure characteristic of the cold-rolled deformed nanocrystalline nickel metal is studied by transmission electron microscopy. The results show that there are step structures nearby the grain boundary (GB), and the contrast of stress field in front of the step corresponds to the step in the shape. It is indicated that the interaction between twins and dislocations is not a necessary condition to realizing the deformation. In the later stage of the deformation when the grain size becomes about 100nm, the deformation can depend upon the moving of the boundary of the stack faults (SFs) which result from the partial dislocations emitted from GBs. However, when the size of SFs grows up, the local internal stress which is in front of the step gradually becomes higher. When this stress reaches a critical value which stops the gliding of the partial dislocations, the SFs will stop to grow up and leave a step structure behind.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.

Computer Simulation Of Random Sphere Packing In An Arbitrarily Shaped Container

Most simulations of random sphere packing concern a cubic or cylindric container with periodic boundary, containers of other shapes are rarely studied. In this paper, a new relaxation algorithm with pre-expanding procedure for random sphere packing in an arbitrarily shaped container is presented. Boundaries of the container are simulated by overlapping spheres which covers the boundary surface of the container. We find 0.4 similar to 0.6 of the overlap rate is a proper value for boundary spheres. The algorithm begins with a random distribution of small internal spheres. Then the expansion and relaxation procedures are performed alternately to increase the packing density. The pre-expanding procedure stops when the packing density of internal spheres reaches a preset value. Following the pre-expanding procedure, the relaxation and shrinking iterations are carried out alternately to reduce the overlaps of internal spheres. The pre-expanding procedure avoids the overflow problem and gives a uniform distribution of initial spheres. Efficiency of the algorithm is increased with the cubic cell background system and double link data structure. Examples show the packing results agree well with both computational and experimental results. Packing density about 0.63 is obtained by the algorithm for random sphere packing in containers of various shapes.

突扩圆管内液-固两相流固体颗粒运动特性的DPM数值模拟

采用考虑颗粒碰撞的欧拉-拉格朗日数值模拟方法(DPM),对水平突扩圆管中液固两相流固体颗粒的碰撞过程进行了数值计算。在模型中,对液相采用欧拉法建立控制方程,对离散颗粒采用拉格朗日方法模拟。采用硬球模型描述颗粒间的碰撞作用。计算结果表明,该模型可以真实地模拟液固两相流中固体颗粒运动的动态变化过程以及颗粒的非均匀分布特征,从单颗粒层次上提供颗粒的运动信息,这有助于深入研究液固两相流中固体颗粒的运动规律

突然扩张流槽中血管内皮细胞表面粘附蛋白表达的研究

动脉粥样硬化的非随机分布与当地的血流动力环境有关，为了研究复杂的流体动力学条件对血和皮细胞生理功能的影响，构建了平行板式平直流槽和突然扩张流槽，通过数值模拟分析了流型的特征，并探讨流型改变对人脐静脉血管内皮细胞表面粘附蛋白表达的影响，发现突然扩张流槽中流动的空间变化使得总体细胞表面粘附蛋白ICAM-1的表达显著高于平直流槽中的均匀定常剪切作用，表明局部流动空间变化的性质可以影响血管内皮细胞的功能。

Transition Between Detonation and Deflagration in a Tube with a Cavity

In this paper, the transition of a detonation from deflagration was investigated numerically while a detonation wave propagates in a tube with a sudden change in cross section, referred to as the expansion cavity. The dispersion-controlled scheme was adopted to solve Euler equations of axis-symmetric flows implemented with detailed chemical reaction kinetics of hydrogen-oxygen (or hydrogen-air) mixture. The fractional step method was applied to treat the stiff problems of chemical reaction flow. It is observed that phenomena of detonation quenching and reigniting appear when the planar detonation front diffracts at the vertex of the expansion cavity entrance. Numerical results show that detonation front in mixture of higher sensitivity keeps its substantial coupled structure when it propagates into the expansion cavity. However, the leading shock wave decouples with the combustion zone if mixture of lower sensitivity was set as the initial gas.

On the origin of shear banding instability in metallic glasses

To uncover the physical origin of shear-banding instability in metallic glass (MG), a theoretical description of thermo-mechanical deformation of MG undergoing one-dimensional simple shearing is presented. The coupled thermo-mechanical model takes into account the momentum balance, the energy balance and the dynamics of free volume. The interplay between free-volume production and temperature increase being two potential causes for shear-banding instability is examined on the basis of the homogeneous solution. It is found that the free-volume production facilitates the sudden increase in the temperature before instability and vice versa. A rigorous linear perturbation analysis is used to examine the inhomogeneous deformation, during which the onset criteria and the internal length and time scales for three types of instabilities, namely free-volume softening, thermal softening and coupling softening, are clearly revealed. The shear-banding instability originating from sole free-volume softening takes place easier and faster than that due to sole thermal softening, and dominates in the coupling softening. Furthermore, the coupled thermo-mechanical shear-band analysis does show that an initial slight distribution of local free volume can incur significant strain localization, producing a shear band. During such a localization process, the local free-volume creation occurs indeed prior to the increase in local temperature, indicating that the former is the cause of shear localization, whereas the latter is its consequence. Finally, extension of the above model to include the shear-induced dilatation shows that such dilatation facilitates the shear instability in metallic glasses.

压缩真空环境中无相互作用两量子比特系统的纠缠动力学

本文研究了压缩真空中无相互作用双原子的纠缠动力学特性。通过分析不同初始纠缠态的演化，发现压缩真空中纠缠原子失去纠缠的速度比在普通真空中更快，并且压缩越大纠缠衰减越快。可以用concurrence和可分性“距离”Lambda的时间演化来解释这种不同寻常的纠缠突然消失现象。

Using drift nets to capture early life stages and monitor spawning of the Yangtze River Chinese sturgeon (Acipenser sinensis)

P>A sampling system for capturing sturgeon eggs using a D-shaped bottom anchored drift net was used to capture early life stages (ELS) of Chinese sturgeon, Acipenser sinensis, and monitor annual spawning success at Yichang on the Yangtze River, 1996-2004, before and just after the Three Gorges Dam began operation. Captured were 96 875 ELS (early life stages: eggs, yolk-sac larvae = eleuthero embryos, and larvae); most were eggs and only 2477 were yolk-sac larvae. Most ELS were captured in the main river channel and inside the bend at the Yichang spawning reach. Yolk-sac larvae were captured for a maximum of 3 days after hatching began, indicating quick dispersal downstream. The back-calculated day of egg fertilization over the eight years indicated a maximum spawning window of 23 days (20 October-10 November). Spawning in all years was restricted temporally, occurred mostly at night and during one or two spawning periods, each lasting several days. The brief temporal spawning window may reduce egg predation by opportunistic predators by flooding the river bottom with millions of eggs. During 1996-2002, the percentage of fertilized eggs in an annual 20-egg sample was between 63.5 to 94.1%; however, in 2003 the percentage fertilized was only 23.8%. This sudden decline may be related to the altered environmental conditions at Yichang caused by operation of the Three Gorges Dam. Further studies are needed to monitor spawning and changes in egg fertilization in this threatened population.