Transient Flows in Low Permeability Porous Media with Composite Start-Up Pressure Gradient

A new mathematical model for the transient flow in the composite low permeability is established. It is solved by FEM with different boundary conditions such as infinite, circular closed and constant pressure boundary conditions. The typical curves for transient wellbore pressure have been presented. It is shown that the pressure and pressure derivative curves with composite start-up pressure gradients have different slopes which are depended on the start-up pressure gradients and the mobility radios in different regions. The boundary effects are the same as the normal reservoirs without start-up pressure gradients. The study provides a new tool to analyze the transient pressure test data in the low permeability reservoir.

Investigation on a pressure-gradient fiber laser hydrophone

In this paper, a pressure-gradient fiber laser hydrophone is demonstrated. Two brass diaphragms are installed at the end of a metal cylinder as sensing elements. A distributed feedback fiber laser, fixed at the center of the two diaphragms, is elongated or shortened due to the acoustic wave. There are two orifices at the middle of the cylinder. So this structure can work as a pressure-gradient microphone in the acoustic field. Furthermore, the hydrostatic pressure is self-compensated and an ultra-thin dimension is achieved. Theoretical analysis is given based on the electro-acoustic theory. Field trials are carried out to test the performance of the hydrophone. A sensitivity of 100 nm MPa-1 has been achieved. Due to the small dimensions, no directivity is found in the test.

Flow Characteristics of Liquid with Pressure-Dependent Viscosities in Microtubes

It is obvious that the pressure gradient alone, the axial direction in a pipe flow keeps constant according to the Haoen-Poiseuille equation. However, recent experiments indicated that the distribution of the pressure seemed no longer linear for liquid flows in microtubes driven by high pressure (1-30MPa). Based on H-P equation with slip boundary condition and Bridgman's relation of viscosity vs. static pressure, the nonlinear distribution of pressure along the axial direction is analyzed in this paper. The revised standard Poiseuille number with the effect of pressure-dependent viscosity taken into account agrees well with the experimental results. Therefore, the dependence of the viscosity on the pressure is one of the dominating, factors under high driven pressure, and is represented by an important property coefficient et of the liquid.

Pressure Drop Studies of Gas-Liquid Bubbly Flow in Vertical Upward Pipeline

This paper describes the experimental and theoretical studies of gas-liquid bubbly flow in vertical upward pipeline carried out at Institute of Mechanics, Chinese Academy of Sciences. Bubbly flow in a vertical pipe with a 3 m long and 5 cm inner diameter plexiglass pipe was experimentally investigated, and studies carried out on the relationship between superficial velocities of the liquid and gas phases and pressure gradient is described. The developed drift-flux model applied to gas-liquid bubbly flow is presented, and the results are compared against the experimental data measured by ours in air/water vertical pipes.

Investigations of phase inversion and frictional pressure gradients in upward and downward oil–water flow in vertical pipes

The present study has attempted to investigate phase inversion and frictional pressure gradients during simultaneous vertical flow of oil and water two-phase through upward and downward pipes. The liquids selected were white oil (44 mPa s viscosity and 860 kg/m3 density) and water. The measurements were made for phase velocities varying from 0 to 1.24 m/s for water and from 0 to 1.87 m/s for oil, respectively. Experiments were carried either by keeping the mixture velocity constant and increasing the dispersed phase fraction or by keeping the continuous phase superficial velocity constant and increasing the dispersed phase superficial velocity. From the experimental results, it is shown that the frictional pressure gradient reaches to its lower value at the phase inversion point in this work. The points of phase inversion are always close to an input oil fraction of 0.8 for upward flow and of 0.75 for downward flow, respectively. A few models published in the literature are used to predict the phase inversion point and to compare the results with available experimental data. Suitable methods are suggested to predict the critical oil holdup at phase inversion based on the different viscosity ratio ranges. Furthermore, the frictional pressure gradient is analyzed with several suitable theoretical models according to the existing flow patterns. The analysis reveals that both the theoretical curves and the experimental data exhibit the same trend and the overall agreement of predicted values with experimental data is good, especially for a high oil fraction.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

DNS of Compressible Turbulent Boundary Layer Over a Blunt Wedge

Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suction is located downstream of the sound-speed point. Statistics are studied and compared with the results in incompressible flat-plate boundary layer. The mean pressure gradient effects on the vortex structure are studied.

Study of drag reduction by gas injection for power-law fluid flow in horizontal stratified and slug flow regimes

In this work, the drag reduction by gas injection for power-law fluid flow in stratified and slug flow regimes has been studied. Experimentswere conducted to measure the pressure gradient within air/CMC solutions in a horizontal Plexiglas pipe that had a diameter of 50mm and a length of 30 m. The drag reduction ratio in stratified flow regime was predicted using the two-fluid model. The results showed that the drag reduction should occur over the large range of the liquid holdup when the flow behaviour index remained at the low value. Furthermore, for turbulent gas-laminar liquid stratified flow, the drag reduction by gas injection for Newtonian fluid was more effective than that for shear-shinning fluid, when the dimensionless liquid height remained in the area of high value. The pressure gradient model for a gas/Newtonian liquid slug flow was extended to liquids possessing the Ostwald–de Waele power law model. The proposed model was validated against 340 experimental data point over a wide range of operating conditions, fluid characteristics and pipe diameters. The dimensionless pressure drop predicted was well inside the 20% deviation region for most of the experimental data. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase slug flows.

Basic equations for suspension flows with interphase mass-transfer

In the case of suspension flows, the rate of interphase momentum transfer M(k) and that of interphase energy transfer E(k), which were expressed as a sum of infinite discontinuities by Ishii, have been reduced to the sum of several terms which have concise physical significance. M(k) is composed of the following terms: (i) the momentum carried by the interphase mass transfer; (ii) the interphase drag force due to the relative motion between phases; (iii) the interphase force produced by the concentration gradient of the dispersed phase in a pressure field. And E(k) is composed of the following four terms, that is, the energy carried by the interphase mass transfer, the work produced by the interphase forces of the second and third parts above, and the heat transfer between phases. It is concluded from the results that (i) the term, (-alpha-k-nabla-p), which is related to the pressure gradient in the momentum equation, can be derived from the basic conservation laws without introducing the "shared-pressure presumption"; (ii) the mean velocity of the action point of the interphase drag is the mean velocity of the interface displacement, upsilonBAR-i. It is approximately equal to the mean velocity of the dispersed phase, upsilonBAR-d. Hence the work terms produced by the drag forces are f(dc) . upsilonBAR-d, and f(cd) . upsilonBAR-d, respectively, with upsilonBAR-i not being replaced by the mean velocity of the continuous phase, upsilonBAR-c; (iii) by analogy, the terms of the momentum transfer due to phase change are upsilonBAR-d-GAMMA-c, and upsilonBAR-d-GAMMA-d, respectively; (iv) since the transformation between explicit heat and latent heat occurs in the process of phase change, the algebraic sum of the heat transfer between phases is not equal to zero. Q(ic) and Q(id) are composed of the explicit heat and latent heat, so that the sum Q(ic) + Q(id)) is equal to zero.

油藏数值模拟中井筒压力梯度的校正方法

有限元法用于油藏数值模拟具有独特的优越性.由于井筒附近流动的特殊性,直接采用数值模拟得到的压力梯度计算油井产量会导致较大的误差.本文分析了数值模拟中井筒压力梯度产生误差的原因,在此基础上提出了计算井筒压力梯度的校正公式.计算结果表明,压力梯度校正公式可显著提高油井产量的计算精度,并能有效地减少井筒附近的网格剖分数量,从而提高了计算效率.

Occurrence of spanning of a submarine pipeline with initial embedment

For better understanding the mechanism of the occurrence of pipeline span for a pipeline with initial embedment, physical and numerical methods are adopted in this study. Experimental observations show that there often exist three characteristic phases in the process of the partially embedded pipeline being suspended: (a) local scour around pipe; (b) onset of soil erosion beneath pipe; and (c) complete suspension of pipe. The effects of local scour on the onset of soil erosion beneath the pipe are much less than those of soil seepage failure induced by the pressure drop. Based on the above observations and analyses, the mechanism of the occurrence of pipeline spanning is analyzed numerically in view of soil seepage failure. In the numerical analyses, the current-induced pressure along the soil surface in the vicinity of the pipe (i.e. the pressure drop) is firstly obtained by solving the N-S equations, thereafter the seepage flow in the soil is calculated with the obtained pressure drop as the boundary conditions along the soil surface. Numerical results indicate that the seepage failure (or piping) may occur at the exit of the seepage path when the pressure gradient gets larger than the critical value. The numerical treatment provides a practical tool for evaluating the potentials for the occurrence of pipe span due to the soil seepage failure.

The influence of Saffman lift force on nanoparticle concentration distribution

The lift force on a spherical nanoparticle near a wall in micro/nanofluidics has not received sufficient attention so far. In this letter the concentration of 200 nm particles is measured at 0.25–2.0 m to a wall in a microchannel with pressure-driven de-ionized water flow pressure gradient 0–2000 kPa/m . The measured data show the influence of the lift force on the nanoparticle concentration distribution. By introducing the Saffman lift force into the Nernst–Planck equation near a wall, we find that the lift force is dominant at the range of 2

Numerical Investigation On The Flowfield Of "Swallowtail" Cavity For Supersonic Mixing Enhancement

A "swallowtail" cavity for the supersonic combustor was proposed to serve as an efficient flame holder for scramjets by enhancing the mass exchange between the cavity and the main flow. A numerical study on the "swallowtail" cavity was conducted by solving the three-dimensional Reynolds-averaged Navier-Stokes equations implemented with a k-epsilon turbulence model in a multi-block mesh. Turbulence model and numerical algorithms were validated first, and then test cases were calculated to investigate into the mechanism of cavity flows. Numerical results demonstrated that the certain mass in the supersonic main flow was sucked into the cavity and moved spirally toward the combustor walls. After that, the flow went out of the cavity at its lateral end, and finally was efficiently mixed with the main flow. The comparison between the "swallowtail" cavity and the conventional one showed that the mass exchanged between the cavity and the main flow was enhanced by the lateral flow that was induced due to the pressure gradient inside the cavity and was driven by the three-dimensional vortex ring generated from the "swallowtail" cavity structure.

Initiation and development of water film by seepage

When water seeps upwards through a saturated soil layer, the soil layer may become instability and water films occur and develop. Water film serves as a natural sliding surface because of its very small friction. Accordingly, debris flow may happen. To investigate this phenomenon, a pseudo-three-phase media is presented first. Then discontinuity method is used to analyze the expansion velocity of water film. Finally, perturbation method is used to analyze the case that a water flow is forced to seep upwards through the soil layer while the movement of the skeleton may be neglected relative to that of water. The theoretical evolutions of pore pressure gradient, effective stress, water velocity, the porosity and the eroded fine grains are obtained. It can be seen clearly that with the erosion and re-deposited of fine grains, permeability at some positions in the soil layer becomes smaller and smaller and, the pore pressure gradient becomes bigger and bigger, while the effective stress becomes smaller and smaller. When the effective stress equals zero, e.f. liquefaction, the water film occurs. It is shown also that once a water film occurs, it will be expanded in a speed of (U)(t)/(1 - E >).

Probability of Sediment Incipient Motion Under Complex Flows

Presented in this paper is a mathematical model to calculate the probability of the sediment incipient motion, in which the effects of the fluctuating pressure and the seepage are considered. The instantaneous bed shear velocity and the pressure gradient on the bed downstream of the backward-facing step flow are obtained according to the PIV measurements. It is found that the instantaneous pressure gradient on the bed obeys normal distribution. The probability of the sediment incipient motion on the bed downstream of the backward-facing step flow is given by the mathematical model. The predicted results agree well with the experiment in the region downstream of the reattachment point while a large discrepancy between the theory and experiment is seen in the region near the reattachment point. The possible reasons for this discrepancy are discussed.