Numerical Simulation of Particle Separation in An Oil-Sand Separator

The gathering systems of crude oil are greatly endangered by the fine sand and soil in oil. Up to now , how to separate sand from the viscid oil is still a technical problem for oil production home or abroad. Recently , Institute of Mechanics in Chinese Academy of Sciences has developed a new type of oil-sand separator , which has been applied successfully in oil field in situ. In this paper, the numerical method of vortex-stream function is used to predict the liquid-solid separating course and the efficiency for this oil-sand separator. Results show that the viscosity and particle diameter have much influence on the particle motion. The calculating separating efficiency is compared with that of experiment and indicates that this method can be used to model the complex two-phase flow in the separator.

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes, equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tolimien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

Experimental Investigation and Numerical Simulation on the Effect of Fissure Water Pressure in Vertical Sliding Surface

This paper studies the effect of fissure water pressure in different fractures on the critical angle of landslide by laboratory investigation and numerical simulation in order to understand the mechanisms of fissure water pressure on landslide stability. Laboratory observations show that the effect of fissure water pressure on the critical angle of landslide is little when the distance between water-holding fracture and slope toe is three times greater than the depth of fissure water. These experimental results are also simulated by a three-dimensional face-to-face contact discrete element method. This method has included the fissure water pressure and can accurately calculate the critical angle of jointed slope when fissure water pressure in vertical sliding surface exists. Numerical results are in good agreement with experimental observations. It is revealed that the location of water-holding structural surface is important to landslide stability. The ratio of the distance between water-holding fissure and slope toe to the depth of fissure water is a key parameter to justify the effect of fissure water pressure on the critical angle of landslide.

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.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Numerical Investigation on Laser Transformation Hardening with Different Temporal Pulse Shapes

A 3-D numerical model for pulsed laser transformation hardening (LTH) is developed using the finite element method. In this model, laser spatial and temporal intensity distribution, temperature-dependent thermophysical properties of material, and multi-phase transformations are considered. The influence of laser temporal pulse shape on connectivity of hardened zone, maximum surface temperature of material and hardening depth is numerically investigated at different pulse energy levels. Results indicate that these hardening parameters are strongly dependent on the temporal pulse shape. For the rectangular temporal pulse shape, the temperature field obtained from this model is in excellent agreement with analytical solution, and the predicted hardening depth is favorably compared with experimental one. It should be pointed out that appropriate temporal pulse shape should be selected according to pulse energy level in order to achieve desirable hardening quality under certain laser spatial intensity distribution.

Numerical Analysis of LEC Growth of GaAs with an Axial Magnetic Field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Numerical Simulation of 3D-Thermal Field with Phase Transformation

以激光熔凝表面强韧化处理为背景,应用空间的弹塑性有限单元和高精度的数值算法、同时考虑材料组织性能的变化来模拟材料的温度场。主要研究激光熔凝加工中瞬时温度场数值模拟,同时考虑相变潜热的影响,为第二步热应力场及残余应力的数值模拟做准备。最后用算例验证了模型的正确性,并给出了不同时刻温度场的分布。

Teracluster LSSC-II - Its Designing Principles and Applications in Large Scale Numerical Simulations

The teracluster LSSC-II installed at the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences is one of the most powerful PC clusters in China. It has a peek performance of 2Tflops. With a Linpack performance of 1.04Tflops, it is ranked at the 43rd place in the 20th TOP500 List (November 2002), 51st place in the 21st TOP500 List (June 2003), and the 82nd place in the 22nd TOP500 List (November 2003) with a new Linpack performance of 1.3Tflops. In this paper, we present some design principles of this cluster, as well as its applications in some largescale numerical simulations.

Numerical Investigation of Toroidal Shock Wave Focusing in a Cylindrical Chamber

In this paper, focusing of a toroidal shock wave propagating from an annular shock tube into a cylindrical chamber was investigated numerically with the dispersion controlled dissipation (DCD) scheme. The first case for an incident Mach number of 1.5 was conducted and compared with experiments for validation. Then, several cases were calculated for higher incident Mach numbers varying from 2.0 to 5.0, and complicated flow structures were observed. The numerical study was mainly focused on two aspects: focusing process and flow structures. The process, including diffraction, focusing, and reflection, is displayed to reveal the focusing mechanism, and the flow structures at different incident. Mach numbers are used to demonstrate shock reflection styles and focusing characteristics.

Numerical Simulation of Transient Thermal Field in Laser Melting Process

Numerical simulation of thermal field was studied in laser processing. The 3-D finite element model of transient thermal calculation is given by thermal conductive equation. The effects of phase transformation latent are considered. Numerical example is given to verify the model. Finally the real example of transient thermal field is given.

Numerical Simulation of Pulsed Liquid Fluidized Bed and Its Experimental Validation

The pulsed liquid fluidized bed was studied using numerical simulation and experimental methods, The area-averaged two-fluid model (TFM) was used to simulate the pulsed fluidization. The bed expansion and collapse processes were simulated first and the phenomena obtained from the calculation were consistent with our previous experiments and observations. In the pulsed fluidization, the variation of bed height, the variations of particle velocity and concentration distribution were obtained and analyzed. Experiments were carried out to validate the simulation results. The pressure variation with time at different locations was measured using pressure transducers and compared with the simulated results. The variations of bed height and particle concentration distribution were recorded using a digital video camera recorder. The results were consistent with the simulation results as a whole.

Numerical Simulation on the Micro-Scale Bending Induced by Pulsed Laser Beam

利用自行研制的含热传导、冲击动力学大、变形有限元程序,模拟了小尺寸梁在脉冲激光加热条件下的变形过程。在此基础上,利用商用程序模拟了冷却及残余应力的产生,研究了激光参数（强度及分布）等对于微弯曲的影响。数值模拟结果与文献中的实验观察相吻合。

Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.