A coupled model for multiphase fluid flow and sedimentation deformation in oil reservoir and its numerical simulation

A mathematical model for coupled multiphase fluid flow and sedimentation deformation is developed based on fluid-solid interaction mechanism. A finite difference-finite element numerical approach is presented. The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances, and the coupled model has practical significance for oilfield development.

Experimental study of thermocapillary convection on a liquid bridge consisting of fluid with low Prandtl number in a floating half-zone

A device of mercury liquid bridge of floating half-zone is designed to experimentally explore thermocapillary convection and its instability of a low Prandtl number liquid. Noncontacted diagnostic techniques were developed to monitor surface flow and surface deformation. The surface flow and the influence of a growing surface film (or skin) on the flow were observed. It is shown that the film is a key factor in changing the behavior associated with the thermocapillary convection. The experiment indicates that the critical Marangoni number should be much higher than that expected by the numerical simulation. The condition and process of surface film growth are discussed. The surface oscillation of the mercury bridge wrapped with ''dirt-film'' was observed, and the characteristics and the frequency associated with this oscillation are given.

Two-dimensional fluid model simulation of bell jar top inductively coupled plasma

In the present paper, argon (Ar) plasmas in a bell jar inductively coupled plasma (ICP) source are systematically studied over pressures from 5 to 20 mtorr and power inputs from 0.2 to 0.5 kW. In this study, both a two-dimensional (2-D) fluid model simulation and global model calculation are compared, The 2-D fluid model simulation with a self-consistent power deposition is developed to describe the Ar plasma behavior as well as predict the plasma parameter distributions, Finally, a quantitative comparison between the global model and the fluid model is made to test their validity.

Compact finite difference - Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

Potential energy model for turbulent entrainment across density interface in a linearly stratified fluid

A potential energy model is developed for turbulent entrainment in the absence of mean shear in a linearly stratified fluid. The relation between the entrainment distance D and the time t and the relation between dimensionless entrainment rate E and the local Richardson number are obtained. An experiment is made for examination. The experimental results are in good agreement with the model, in which the dimensionless entrainment distance D is given by DBAR = A(i)(SBAR)-1/4(fBAR)1/2(tBAR)1/8, where A(i) is the proportional coefficient, S is the dimensionless stroke, fBAR is the dimensionless frequency of the grid oscillation, tBAR the dimensionless time.

Growth of zero-mean-shear turbulent-mixed layer with suspended particles in a 2-layer fluid system

The growth behaviour of zero-mean-shear turbulent-mixed layer containing suspended solid particles has been studied experimentally and analysed theoretically in a two-layer fluid system. The potential model for estimating the turbulent entrainment rate of the mixed layer has also been suggested, including the results of the turbulent entrainment for pure two-layer fluid. The experimental results show that the entrainment behaviour of a mixed layer with the suspended particles is well described by the model. The relationship between the entrainment distance and the time, and the variation of the dimensionless entrainment rate E with the local Richardson number Ri1 for the suspended particles differ from that for the pure two-layer fluid by the factors-eta-1/5 and eta-1, respectively, where eta = 1 + sigma-0-DELTA-rho/DELTA-rho-0.

Hierarchal structure of the simplified Navier-Stokes equations (SNSE) and its mechanical connotation and application

The hierarchial structure and mathematical property of the simplified Navier-Stokesequations (SNSE) are studied for viscous flow over a sphere and a jet of compressible flu-id. All kinds of the hierarchial SNSE can be divided into three types according to theirmathematical property and also into five groups according to their physical content. Amultilayers structure model for viscous shear flow with a main stream direction is pre-sented. For the example of viscous incompressible flow over a flat plate there existthree layers for both the separated flow and the attached flow; the character of thetransition from the three layers of attached flow to those of separated flow is elucidated.A concept of transition layer being situated between the viscous layer and inviscidlayer is introduced. The transition layer features the interaction between viscous flow andinviscid flow. The inner-outer-layers-matched SNSE proposed by the present author inthe past is developed into the layers matched (LsM)-SNSE.

Fission of solitary wave in stratified fluid

In this paper, we study the fission of a solitary wave in the stratified fluid with a free surface. It has been discovered that there is no difference between the fissions of the internal solitary waves in odd or even modes, and the effect of the stratification on the fission of a surface solitary wave can almost be neglected

Visualization of turbulence and internal waves in a stratified fluid

A new method for measuring the density, temperature and velocity of N2 gas flow by laser induced biacetyl phosphorescence is proposed. The characteristics of the laser induced phosphorescence of biacetyl mixed with N2 are investigated both in static gas and in one-dimensional flow along a pipe with constant cross section. The theoretical and experimental investigations show that the temperature and density of N2 gas flow could be measured by observing the phosphorescence lifetime and initial intensity of biacetyl triplet (3Au) respectively. The velocity could be measured by observing the time-of-flight of the phosphorescent gas after pulsed laser excitation. The prospect of this method is also discussed.