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.

Asymptotic analysis of thin plates under normal load and horizontal edge thrust

We consider the radially symmetric nonlinear von Kármán plate equations for circular or annular plates in the limit of small thickness. The loads on the plate consist of a radially symmetric pressure load and a uniform edge load. The dependence of the steady states on the edge load and thickness is studied using asymptotics as well as numerical calculations. The von Kármán plate equations are a singular perturbation of the Fӧppl membrane equation in the asymptotic limit of small thickness. We study the role of compressive membrane solutions in the small thickness asymptotic behavior of the plate solutions.

We give evidence for the existence of a singular compressive solution for the circular membrane and show by a singular perturbation expansion that the nonsingular compressive solution approach this singular solution as the radial stress at the center of the plate vanishes. In this limit, an infinite number of folds occur with respect to the edge load. Similar behavior is observed for the annular membrane with zero edge load at the inner radius in the limit as the circumferential stress vanishes.

We develop multiscale expansions, which are asymptotic to members of this family for plates with edges that are elastically supported against rotation. At some thicknesses this approximation breaks down and a boundary layer appears at the center of the plate. In the limit of small normal load, the points of breakdown approach the bifurcation points corresponding to buckling of the nondeflected state. A uniform asymptotic expansion for small thickness combining the boundary layer with a multiscale approximation of the outer solution is developed for this case. These approximations complement the well known boundary layer expansions based on tensile membrane solutions in describing the bending and stretching of thin plates. The approximation becomes inconsistent as the clamped state is approached by increasing the resistance against rotation at the edge. We prove that such an expansion for the clamped circular plate cannot exist unless the pressure load is self-equilibrating.

Studies in vortex motion

This thesis covers four different problems in the understanding of vortex sheets, and these are presented in four chapters.

In Chapter 1, free streamline theory is used to determine the steady solutions of an array of identical, hollow or stagnant core vortices in an inviscid, incompressible fluid. Assuming the array is symmetric to rotation through π radians about an axis through any vortex centre, there are two solutions or no solutions depending on whether A^(1/2)/L is less than or greater than 0.38 where A is the area of the vortex and L is the separation distance. Stability analysis shows that the more deformed shape is unstable to infinitesimal symmetric disturbances which leave the centres of the vortices undisplaced.

Chapter 2 is concerned with the roll-up of vortex sheets in homogeneous fluid. The flow over conventional and ring wings is used to test the method of Fink and Soh (1974). Despite modifications which improve the accuracy of the method, unphysical results occur. A possible explanation for this is that small scales are important and an alternate method based on "Cloud-in-Cell" techniques is introduced. The results show small scale growth and amalgamation into larger structures.

The motion of a buoyant pair of line vortices of opposite circulation is considered in Chapter 3. The density difference between the fluid carried by the vortices and the fluid outside is considered small, so that the Boussinesq approximation may be used. A macroscopic model is developed which shows the formation of a detrainment filament and this is included as a modification to the model. The results agree well with the numerical solution as developed by Hill (1975b) and show that after an initial slowdown, the vortices begin to accelerate downwards.

Chapter 4 reproduces completely a paper that has already been published (Baker, Barker, Bofah and Saffman (1974)) on the effect of "vortex wandering" on the measurement of velocity profiles of the trailing vortices behind a wing.

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.

Topics in linear and nonlinear dispersive waves

The various singularities and instabilities which arise in the modulation theory of dispersive wavetrains are studied. Primary interest is in the theory of nonlinear waves, but a study of associated questions in linear theory provides background information and is of independent interest.

The full modulation theory is developed in general terms. In the first approximation for slow modulations, the modulation equations are solved. In both the linear and nonlinear theories, singularities and regions of multivalued modulations are predicted. Higher order effects are considered to evaluate this first order theory. An improved approximation is presented which gives the true behavior in the singular regions. For the linear case, the end result can be interpreted as the overlap of elementary wavetrains. In the nonlinear case, it is found that a sufficiently strong nonlinearity prevents this overlap. Transition zones with a predictable structure replace the singular regions.

For linear problems, exact solutions are found by Fourier integrals and other superposition techniques. These show the true behavior when breaking modulations are predicted.

A numerical study is made for the anharmonic lattice to assess the nonlinear theory. This confirms the theoretical predictions of nonlinear group velocities, group splitting, and wavetrain instability, as well as higher order effects in the singular regions.

Numerical methods for ill-posed, linear problems

A means of assessing the effectiveness of methods used in the numerical solution of various linear ill-posed problems is outlined. Two methods: Tikhonov' s method of regularization and the quasireversibility method of Lattès and Lions are appraised from this point of view.

In the former method, Tikhonov provides a useful means for incorporating a constraint into numerical algorithms. The analysis suggests that the approach can be generalized to embody constraints other than those employed by Tikhonov. This is effected and the general "T-method" is the result.

A T-method is used on an extended version of the backwards heat equation with spatially variable coefficients. Numerical computations based upon it are performed.

The statistical method developed by Franklin is shown to have an interpretation as a T-method. This interpretation, although somewhat loose, does explain some empirical convergence properties which are difficult to pin down via a purely statistical argument.

Numerical solution of parabolic equations by the box scheme

The box scheme proposed by H. B. Keller is a numerical method for solving parabolic partial differential equations. We give a convergence proof of this scheme for the heat equation, for a linear parabolic system, and for a class of nonlinear parabolic equations. Von Neumann stability is shown to hold for the box scheme combined with the method of fractional steps to solve the two-dimensional heat equation. Computations were performed on Burgers' equation with three different initial conditions, and Richardson extrapolation is shown to be effective.