New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

Numerical Study on Critical Wedge Angle of Cellular Detonation Reflections

The critical wedge angle (CWA) for the transition from regular reflection (RR) to Mach reflection (MR) of a cellular detonation wave is studied numerically by an improved space-time conservation element and solution element method together with a two-step chemical reaction model. The accuracy of that numerical way is verified by simulating cellular detonation reflections at a 19.3∘ wedge. The planar and cellular detonation reflections over 45∘–55∘ wedges are also simulated. When the cellular detonation wave is over a 50∘ wedge, numerical results show a new phenomenon that RR and MR occur alternately. The transition process between RR and MR is investigated with the local pressure contours. Numerical analysis shows that the cellular structure is the essential reason for the new phenomenon and the CWA of detonation reflection is not a certain angle but an angle range.

An Integrated Evolutionary Algorithm for Expensive Global Optimization

We propose an integrated algorithm named low dimensional simplex evolution extension (LDSEE) for expensive global optimization in which only a very limited number of function evaluations is allowed. The new algorithm accelerates an existing global optimization, low dimensional simplex evolution (LDSE), by using radial basis function (RBF) interpolation and tabu search. Different from other expensive global optimization methods, LDSEE integrates the RBF interpolation and tabu search with the LDSE algorithm rather than just calling existing global optimization algorithms as subroutines. As a result, it can keep a good balance between the model approximation and the global search. Meanwhile it is self-contained. It does not rely on other GO algorithms and is very easy to use. Numerical results show that it is a competitive alternative for expensive global optimization.

A numerical investigation on the flow physics of koi carp (Cyprinus carpio koi) burst-and-coast swimming

Abstract—Burst-and-coast is the most common locomotion type in freely routine swimming of koi carps (Cyprinus carpio koi), which consists of a burst phase and a coast phase in each cycle and mostly leads to a straight-line trajectory. Combining with the tracking experiment, the flow physics of koi carp’s burst-andcoast swimming is investigated using a novel integrated CFD method solving the body-fluid interaction problem. The dynamical equations of a deforming body are formulated. Following that, the loose-coupled equations of the body dynamics and the fluid dynamics are numerically solved with the integrated method. The two burst modes, MT (Multiple Tail-beat) and HT (Half Tail-beat), which have been reported by the experiments, are investigated by numerical simulations in this paper. The body kinematics is predicted and the flow physics is visualized, which are in good agreement with the corresponding experiments. Furthermore, the optimization on the energy cost and several critical control mechanisms in burst-and-coast swimming of koi carps are explored, by varying the parameters in its selfpropelled swimming. In this paper, energetics is measured by the two mechanical quantities, total output power CP and Froude efficiency Fr. Results and discussion show that from the standpoint of mechanical energy, burst-and-coast swimming does not actually save energy comparing with steady swimming at the same average speed, in that frequently changing of speed leads to decrease of efficiency.

Numerical Simulation of Liquefaction Deformation of Sand Layer around Bucket Foundations under Dynamic Loadings

Bucket Foundations under Dynamic Loadings The liquefaction deformation of sand layer around a bucket foundation is simulated under equivalent dynamic ice-induced loadings. A simplified numerical model is presented by taking the bucket-soil interaction into consideration. The development of vertical and horizontal liquefaction deformations are computed under equivalent dynamic ice-induced loadings. Firstly, the numerical model and results are proved to be reliable by comparing them with the centrifuge testing results. Secondly, the factors and the development characteristics of liquefaction deformation are analyzed. Finally, the following numerical simulation results are obtained: the liquefaction deformation of sand layer increases with the increase of loading amplitude and with the decrease of loading frequency and sand skeleton’s strength. The maximum vertical deformation is located on the sand layer surface and 1/4 times of the bucket’s height apart from the bucket’s side wall (loading boundary). The maximum horizontal deformation occurs at the loading boundary. When the dynamic loadings is applied for more than 5 hours, the vertical deformation on the sand layer surface reaches 3 times that at the bottom, and the horizontal deformation at 2.0 times of the bucket height apart from the loading boundary is 3.3% of which on the loading boundary.

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.

NUMERICAL AND EXPERIMENTAL STUDY ON LIQUID-SOLID FLOW IN A HYDROCYCLONE

Hydrocyclones are widely used in industry, of which the geometrical design using CFD techniques is gaining more popularity in recent years. In this study, the Euler-Euler approach and the Reynolds stress model are applied to simulate the liquid-solid flowfield in a hydrocyclone. The methodology is validated by a good agreement between experimental data and numerical results. Within the research range, the simulation indicates that the liquid-solid separation mainly occurs in the conical segment, and increasing conical height or decreasing cylindrical height helps to improve the grade efficiencies of solid particles. Based on these results, two of the same hydrocyclones are designed and installed in series to establish a liquid-solid separation system. Many experiments are then conducted under different conditions, in which the effects of the water cut and the second hydrocyclone on the separation are investigated. The results also confirm that smaller solid particles are more susceptible to the inlet conditions, and the second hydrocyclone plays a more important role as the water cut reduces.

Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments

Smoothed particle hydrodynamics (SPH) is a meshfree particle method based on Lagrangian formulation, and has been widely applied to different areas in engineering and science. This paper presents an overview on the SPH method and its recent developments, including (1) the need for meshfree particle methods, and advantages of SPH, (2) approximation schemes of the conventional SPH method and numerical techniques for deriving SPH formulations for partial differential equations such as the Navier-Stokes (N-S) equations, (3) the role of the smoothing kernel functions and a general approach to construct smoothing kernel functions, (4) kernel and particle consistency for the SPH method, and approaches for restoring particle consistency, (5) several important numerical aspects, and (6) some recent applications of SPH. The paper ends with some concluding remarks.

光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.

Numerical Study On Transient Flow In The Deep Naturally Fractured Reservoir With High Pressure

According to the experimental results and the characteristics of the pressure-sensitive fractured formation, a transient flow model is developed for the deep naturally-fractured reservoirs with different outer boundary conditions. The finite element equations for the model are derived. After generating the unstructured grids in the solution regions, the finite element method is used to calculate the pressure type curves for the pressure-sensitive fractured reservoir with different outer boundaries, such as the infinite boundary, circle boundary and combined linear boundaries, and the characteristics of the type curves are comparatively analyzed. The effects on the pressure curves caused by pressure sensitivity module and the effective radius combined parameter are determined, and the method for calculating the pressure-sensitive reservoir parameters is introduced. By analyzing the real field case in the high temperature and pressure reservoir, the perfect results show that the transient flow model for the pressure-sensitive fractured reservoir in this paper is correct.

A numerical model for predicting the jump of lift force on an oscillating cylinder

The fluid force coefficients on a transversely oscillating cylinder are calculated by applying two- dimensional large eddy simulation method. Considering the ‘‘jump’’ phenomenon of the amplitude of lift coefficient is harmful to the security of the submarine slender structures, the characteristics of this ‘‘jump’’ are dissertated concretely. By comparing with experiment results, we establish a numerical model for predicting the jump of lift force on an oscillating cylinder, providing consultation for revising the hydrodynamic parameters and checking the fatigue life scale design of submarine slender cylindrical structures.

The Experimental and Numerical Studies on Gas Production from Hydrate Reservoir by Depressurization

A set of experimental system to study hydrate dissociation in porous media is built and some experiments on hydrate dissociation by depressurization are carried out. A mathematical model is developed to simulate the hydrate dissociation by depressurization in hydrate-bearing porous media. The model can be used to analyze the effects of the flow of multiphase fluids, the kinetic process and endothermic process of hydrate dissociation, ice-water phase equilibrium, the variation of permeability, convection and conduction on the hydrate dissociation, and gas and water productions. The numerical results agree well with the experimental results, which validate our mathematical model. For a 3-D hydrate reservoir of Class 3, the evolutions of pressure, temperature, and saturations are elucidated and the effects of some main parameters on gas and water rates are analyzed. Numerical results show that gas can be produced effectively from hydrate reservoir in the first stage of depressurization. Then, methods such as thermal stimulation or inhibitor injection should be considered due to the energy deficiency of formation energy. The numerical results for 3-D hydrate reservoir of Class 1 show that the overlying gas hydrate zone can apparently enhance gas rate and prolong life span of gas reservoir.

Numerical and Experimental Investigation On The Prevention of Co Deflagration

The exhaust gases from industrial furnaces contain a huge amount of heat and chemical enthalpy. However, it is hard to recover this energy since exhaust gases invariably contain combustible components such as carbon monoxide (CC). If the CO is unexpectedly ignited during the heat recovery process, deflagration or even detonation could occur, with serious consequences such as complete destruction of the equipment. In order to safely utilize the heat energy contained in exhaust gas, danger of its explosion must be fully avoided. The mechanism of gas deflagration and its prevention must therefore be studied. In this paper, we describe a numerical and experimental investigation of the deflagration process in a semi-opened tube. The results show that, upon ignition, a low-pressure wave initially spreads within the tube and then deflagration begins. For the purpose of preventing deflagration, an appropriate amount of nitrogen was injected into the tube at a fixed position. Both simulation and experimental results have shown that the injection of inert gas can successfully interrupt the deflagration process. The peak value of the deflagration pressure can thereby be reduced by around 50%. (C) 2008 Elsevier Ltd. All rights reserved.

极紫外阿秒脉冲在高次谐波产生过程中引起的非偶极效应

当红外强激光和极紫外（XUV）阿秒脉冲共同作用于原子分子时，电离出去的电子通常会吸收和辐射激光光子而发生能量扩展．讨论了由于XUV阿秒脉冲的短波长与扩展后的电子波包尺度可相比拟时在高次谐波产生过程中引起的非偶极效应．采用彤作为模型分子，并把分子轴置于激光场的传播方向，通过解二维含时薛定谔方程并比较考虑非偶极效应和采用偶极近似两种方法计算得到的结果，两者相比，前者的谐波强度降低，谐波频率向低级次稍有移动，电子能谱的能带内出现了更多的光电子峰．在相同的光电子能量处，两种方法计算得到的信号强度相差2—5倍．并且

Numerical modeling of optical levitation and trapping of the "stuck" particles with a pulsed optical tweezers

We present the theoretical analysis and the numerical modeling of optical levitation and trapping of the stuck particles with a pulsed optical tweezers. In our model, a pulsed laser was used to generate a large gradient force within a short duration that overcame the adhesive interaction between the stuck particles and the surface; and then a low power continuous - wave (cw) laser was used to capture the levitated particle. We describe the gradient force generated by the pulsed optical tweezers and model the binding interaction between the stuck beads and glass surface by the dominative van der Waals force with a randomly distributed binding strength. We numerically calculate the single pulse levitation efficiency for polystyrene beads as the function of the pulse energy, the axial displacement from the surface to the pulsed laser focus and the pulse duration. The result of our numerical modeling is qualitatively consistent with the experimental result. (C) 2005 Optical Society of America.