Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.

NUMERICAL SIMULATION OF 3-D DENSE SOLID-LIQUID TWO-PHASE TURBULENT FLOW IN A NON-CLOGGING MUD PUMP

A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to TECPLOT as the post-processor, the simulation results isvisualized. The results show the main flow characteristics: There exists backflow and aberrantvelocities at inlet area and a relative velocity slip between two phases; A jet-wake flow pattern isdiscerned around the shroud-suction side area; The relative velocity vector of solid phase iscloser to the pressure surface than that of liquid phase and the trend is more obvious with theincrease of diameter; The kinetic energy of turbulence k and the dissipation rate e reach theirpeaks at the corner of pressure and suction surface. The simulation results show a good agreementwith the experimental flow features in the impeller channel, which prove the turbulent model used isvalid and provide a theoretical design basis to non-clogging pumps.

NUMERICAL SIMULATION FOR INCOMPRESSIBLE TURBULENT FLOW IN THE IMPELLER OF SEWAGE PUMP

The 3-D turbulent flow in the impeller of sewage pump is simulated. Thetime-averaged N-S equations and the kappa-epsilon turbulent model is modified. The calculation iscarried out in body-fitted coordinated grid by applying SIMPLE-C algorithm. The calculated velocity,pressure distributions of the turbulent flow in the sewage pump are obtained for the first time,which will be helpful for the optimal design and performance prediction of sewage pumps on the basisof flow field simulation.

Numerical simulation on turbulent flow field in convergentdiveroent nozzle

Because of the complication of turbulence＇s mechanism and law as well as the jet pressure in nozzle is difficult to test by experiment, five turbulent models were applied to numerically simulate the turbulent flow field in convergent-divergent nozzle. Theory analysis and experiment results of mass flow rates conclude that the RNG k-ε model is the most suitable model. The pressure distribution in the convergent-divergent nozzle was revealed by computational fluid dynamic （CFD） simulating on the turbulent flow field under different pressure conditions. The growing conditions of cavitation bubbles were shown; meanwhile, the phenomena in the experiment could be explained. The differential pres- sure between the upstream and downstream in nozzle throat section can improve the cavitating effect of cavitation water jet.

Numerical Simulation on Gas-Solid Two-Phase Turbulent Flow in FCC Riser Reactors(Ⅱ) Numerical Simulation on the Gas-Solid Two-Phase Turbulent Flow

Numerical simulation on the flow,heat transfer and cracking reactions in commercial fluid catalyticcracking(FCC)riser reactors were carried out employing the developed turbulent gas-solid two-phase flow-reac-tion model for FCC riser reactors given in Part Ⅰ of the present paper.Detailed information about the turbulentflow fields in the riser reactor obtained revealed the basic characteristics of the gas-solid two-phase turbulentflows when heat transfer and catalytic cracking reactions were co-existing in the riser.Results showed that thedistributions of the flow,the turbulence kinetic energy and the catalyst particle concentration are not uniform inthe axial,radial and tangential directions.The most complicated part of the riser reactor is the feed injectingzone.The complicated configuration of the turbulent gas-solid two-phase flows would exert a great influence onthe results of interphase heat transfer and cracking reactions.

Numerical Simulation of the Effect of Air Distribution on Turbulent Flow and Combustion in a Tubular Heating Furnace

A three-dimension full-size numerical simulation of the effect of air distribution on turbulent flow and combustion in a tubular heating furnace was carried out. A standard k –ε turbulent model, a simplified PDF combustion model and a discrete ordinate transfer radiation model were used. The hybrid grid combining a structured and a non-structured grid was generated without any simplification of the complicated geometric configuration around the burner. It was found that the multistage combustion could reduce and control the peak value of temperature. At the same time, it was concluded that the amount of primary air had little effect on the global distribution of velocity and temperature in the furnace, but a great effect on that around the burner. It is recommended that 45% - 65% of the total amount of air be taken in in primary air inlets in the furnace. All the results are important to optimize the combustion progress.

Numerical Simulation of Two-Dimensional Flow over Three Cylinders by Lattice Boltzmann Method

The numerical simulation using the multiple relaxation time lattice Boltzmann method （MRT-LBM） is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the moving cylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/ D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D = 0.5 and fv= 0.4 （where A is the forced amplitude, D is the cylinder diameter, and fv stands for the vibrating frequency, respectively）. The calculated results not only indicate that the spacing ratios T/D （T is the center-to-center spacing between the two upstream cylinders） have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D （where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder）.

Direct numerical simulation on relevance of fluctuating velocities and drag reduction in turbulent channel flow with spanwise space-dependent electromagnetic force

Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.

THE NUMERICAL SIMULATION OF TURBULENT FLOWS OF NEWTONIAN AND A SORT OF NON-NEWTONIAN FLUID IN 180-DEGREE SQUARE SECTIONED BEND

Using k- model of turbulence and measured wall functions, turbulent flows of Newtonian (pure water) andasort of non-Newtonian fluid (dilute, drag-reduction solution of polymer) in a 180-degree curved bend were simulated numerically. The calculated results agreed well with the measured velocity profiles. On the basis of calculation and measurement, appropriateness of turbulence model to complicated flow in which the large-scale vortex exists was analyzed and discussed.

NUMERICAL SIMULATION OF 3-D TURBULENT FLOWS OVER DREDGED TRENCHES

A 3- D free surface flow in open channels based on the Reynolds equations with the k-ε turbulence closure model is presented in this paper. Insted of the 'rigid lid' approximation, the solution of the free surface equation is implemented in the velocity-pressure iterative procedure on the basis of the conventional SIMPLE method. This model was used to compute the flow in rectangular channels with trenches dredged across the bottom. The velocity, eddy viscosity coefficient, turbulent shear stress, turbulent kinetic energy and elevation of the free surface can be obtained. The computed results are in good agreement with previous experimental data.

NUMERICAL SIMULATION OF THREE DIMENSIONAL TURBULENT FLOW IN SUDDENLY EXPANDED RECTANGULAR DUCT

In this paper,to simulate the three dimensional turbulent flow in suddenly expanded rectangular duct numerically,the SIMPLEC algorithm is employed to solve the incompressible Navier-Stckes equation with k-εturbulence model.The numerical resulis show well the three dimensional turbulent flow field in the rectangular duct behind the sudden expansion cross-section. and agree.fairly well with the experimental result of the length of the main circumfluence.The numerical method of this paper can be applied to numerical analysis of this kind of turbulent flow.

NUMERICAL SIMULATION OF CHARGED GAS-LIQUID TWO PHASE JET FLOW IN ELECTROSTATIC SPRAYING

Multi-fluid k-e-kp, two phase turbulence model is used to simulate charged gas-liquid two phase coaxial jet, which is the transorting flow field in electrostatic spraying. Compared with the results of experiment, charged gas-liquid two-phase turbulence can be well predicted by this model.

Numerical Research on the Fiber Suspensions in a Turbulent T-shaped Branching Channel Flow

The concentration and orientation of fiber in a turbulent T-shaped branching channel flow are investi-gated numerically. The Reynolds averaged Navier-Stokes equations together with the Reynolds stress turbulent model are solved for the mean flow field and the turbulent kinetic energy. The fluctuating velocities of the fluid are assumed as a random variable with Gaussian distribution whose variance is related to the turbulent kinetic energy. The slender-body theory is used to simulate the fiber motion based on the known mean and fluctuating velocities of the fluid. The results show that at low Reynolds number, fiber concentration is high in the flow separation regions, and fiber orientation throughout the channel is widely distributed with a slight preference of aligning along the horizontal axis. With increasing of Re, the high concentration region disappears, and fiber orientation becomes ho-mogeneous without any preferred direction. At high Reynolds number, fiber concentration increases gradually along the flow direction. The differences in the distribution of concentration and orientation between different fiber aspect ratio are evident only at low Re. Both Re and fiber aspect ratio have small effect on the variance of orientation angle.

Validation of the numerical simulation of fluid flow in the continuous casting tundish

The 3-dimensional fluid flow in a water model of the continuous casting tundish is simulated with the k-ε two-equation turbulence model. The methods to decide the size of wall-adjacent grids and the effects of residuals and mesh size on the simulation accuracy are discussed. The current investigation concludes that the following condition should be satisfied to get accurate enough simulation results: (1) If the dimension of the domain is in the order of cubic meters, the average size of cells in the mesh system should be at least smaller than 30 mm; (2) The normalized non-scaled residual should be reached at least smaller than 10-4.

Numerical Simulation of Interaction Between Laminar Flow and Elastic Sheet

A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.

Computational modeling of free-surface slurry flow problems using particle simulation method

The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed.

Numerical simulation of flow around square cylinder using different low-Reynolds number turbulence models

ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.

Computer Simulation of Turbulent Flow through a Hydraulic Turbine Draft Tube

Based on the Naviev-Stokes equations and the standard κ-ε turbulence model, this paper presents the derivation of the governing equations for the turbulent flow field in a draft tube. The mathematical model for the turbulent flow through a draft tube is set up when the boundary conditions, including the inlet boundary conditions, the outlet boundary conditions and the wall boundary conditions, have been implemented. The governing equations are formulated in a discrete form on a staggered grid system by the finite volume method. The second-order central difference approximation and hybrid scheme are used for discretization. The computation and analysis on internal flow through a draft tube have been carried out by using the simplee algorithm and cfx-tasc flow software so as to obtain the simulated flow fields. The calculation results at the design operating condition for the draft tube are presented in this paper. Thereby, an effective method for simulating the internal flow field in a draft tube has been explored.

Numerical simulation of three-dimensional combustion flows

A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.