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.

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.

Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Numerical simulation of hypersonic viscous flow around space shuttle with method of diffusion analogy

Hypersonic viscous flow around a space shuttle with M(infinity) = 7, Re = 148000 and angle of attack alpha = 5-degrees is simulated numerically with the special Jacobian matrix splitting technique and simplified diffusion analogy method. With the simplified diffusion analogy method the efficiency of computation and resolution of the shock can be improved.

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.