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.

Simulation of Gas-Water Two-Phase Flow in Tight Gas Reservoirs Considering the Gas Slip Effect

A mathematical model for the gas-water two-phase flow in tight gas reservoirs is elaborated.The model can account for the gas slip effect,stress sensitivity,and high-speed non-Darcy factors.The related equations are solved in the framework of a finite element method.The results are validated against those obtained by using the commercial software CMG(Computer Modeling Group software for advanced recovery process simulation).It is shown that the proposed method is reliable.It can capture the fracture rejection characteristics of tight gas reservoirs better than the CMG.A sensitivity analysis of various control factors(initial water saturation,reservoir parameters,and fracturing parameters)affecting the production in tight gas wells is conducted accordingly.Finally,a series of theoretical arguments are provided for a rational and effective development/exploitation of tight sandstone gas reservoirs.

Simulation of liquid cone formation on the tip apex of indium field emission electric propulsion thrusters

Field emission electric propulsion(FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities,making them ideal propulsion devices for micro-nano satellites. However, the detection of certain aspects, such as the evolution process of the liquid cone and the physical quantities at the cone apex, proves challenging due to the minute size of the needle tip and the vacuum environment in which they operate. Consequently, this paper introduces a computational fluid dynamics(CFD) model to gain insight into the formation process of the liquid cone on the tip apex of indium FEEP. The CFD model is based on electrohydrodynamic(EHD) equations and the volume of fluid(VOF) method. The entire cone formation process can be divided into three stages, and the time-dependent characteristics of the physical quantities at the cone apex are investigated. The influences of film thickness, apex radius size and applied voltage are compared.The results indicate a gradual increase in the values of electrostatic stress and surface tension stress at the cone apex over an initial period, followed by a rapid escalation within a short duration.Apex configurations featuring a small radius, thick film and high voltage exhibit a propensity for liquid cone formation, and the cone growth time decreases as the film thickness increases.Moreover, some unstable behavior is observed during the cone formation process.

Numerical Simulation and Analysis of Solid-liquid Two-phase Flow in Centrifugal Pump

The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the two-phase flow pumps, i.e., low overall efficiency and severe abrasion. In this study, the three-dimensional, steady, incompressible, and turbulent solid-liquid two-phase flows in a low-specific-speed centrifugal pump are numerically simulated and analyzed by using a computational fluid dynamics （CFD） code based on the mixture model of the two-phase flow and the RNG k-~ two-equation turbulence model, in which the influences of rotation and curvature are fully taken into account. The coupling between impeller and volute is implemented by means of the frozen rotor method. The simulation results predicted indicate that the solid phase properties in two-phase flow, especially the concentration, the particle diameter and the density, have strong effects on the hydraulic performance of the pump. Both the pump head and the efficiency are reduced with increasing particle diameter or concentration. However, the effect of particle density on the performance is relatively minor. An obvious jet-wake flow structure is presented near the volute tongue and becomes more remarkable with increasing solid phase concentration. The suction side of the blade is subject to much more severe abrasion than the pressure side. The obtained results preliminarily reveal the characteristics of solid-liquid two-phase flow in the centrifugal pump, and are helpful for improvement and empirical correction in the hydraulic design of centrifugal pumps.

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.

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.

Numerical simulation of two-phase flow in fractured porous media using streamline simulation and IMPES methods and comparing results with a commercial software

Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.

Numerical Simulation and Experimental Analysis of the Influence of Asymmetric Pressure Conditions on the Splitting of a Gas-Liquid Two-Phase Flow at a T-Junction

Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid velocities.The flow rates in the two branches were measured accurately to determine how the two considered phases distribute in the two outlets.The experimental results have shown that when the two outlet pressures are asymmetric,the two-phase flow always tends to flow into the outlet which has a lower pressure.As the inlet liquid velocity increases,however,the two-phase flow gradually tends to split evenly.Compared with the experiment results,the pressure difference between the two outlets can be determined more accurately by means of numerical simulation.The trends of experimental results and simulations are in very good agreement.

Numerical simulation of predicting and reducing solid particle erosion of solid-liquid two-phase flow in a choke

Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.

Two-dimensional Dynamics Simulation of Two-phase Debris Flow

To investigate the movement mechanism of debris flow, a two-dimensional, two-phase, depthintegrated model is introduced. The model uses Mohr-Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second-order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.

Numerical simulation of the gas-solid two-phase flow inside the multi-channel nozzle for the surface nanocrystallization induced by the ultrasonic particulate peening

Using a gas-solid two-phase model(a discrete phase model),the authors investigated the flow field inside the multi-channel nozzle for surface nanocrystallization(SNC)induced by the ultrasonic particulate peening(USPP).By computation,the velocity fields of both the gas and the solid phases were simulated and the track of the solid phase was analyzed in detail.It can be found that the velocities of the two phases are able to reach an ultrasonic level;meanwhile,the dispersion width of the solid phase at the nozzle exit is less than that of the gas phase.When particle diameters are less than 5 μm,there is a decreasing trend in the dispersion width of the solid phase with an increase in particle diameters.The trend becomes stable as the particle diameters are greater than 5 μm;in the meantime,the distribution of solid particles is near the axis of the jet flow.The optimal standoff distance between the nozzle and the substrate in the process of USPP is about 120 mm.Simulation results can help improve the design of mass-production-oriented multi-channel nozzles for SNC induced by USPP.

Numerical Simulation of the Gas-solid Two-phase Flows in a Precalciner

The gas-solid two-phase flous of the precalciner were simulated by different multiphase models,such as mixture model,the Enderium model,including mixture and dispersed,and discrete phase model（DPM）,The results of the different multiphase models were analyzed and compared.showing the rationality of the diffusion and mixture of the cenment raic meals and coal poroder some extent Moreover,the results also shose the rationality of the given inlets parameters of actual process of the precalciner.

Finite element equations and numerical simulation of elastic wave propagation in two-phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.

Thermal Simulation for Two-Phase Liquid Cooling 3D-ICs

This work presents an algorithm for simulating more accurate temperature distribution in two-phase liquid cooling for three-dimensional integrated circuits than the state of-the-art methods by utilizing local multi-linear interpolation techniques on heat transfer coefficients between the microchannel and silicon substrate, and considering the interdependence between the thermal conductivity of silicon and temperature values. The experimental results show that the maximum and average errors are only 9.7% and 6.7% compared with the measurements, respectively.

SIMULATION OF GAS AND LIQUID TWO-PHASE FLOW THROUGH THE BLAST FURNACE DROPPING ZONE

A three-dimensional mathematical model,based on differential balances of mass and momentum,hasbeen developed to describe the two-phase flow of gas and liquid through the dropping zone of the blast fur-nace.Agreement between observed and calculated values verifies the validity of this model.On the basis of this model,various parameters for the surrounding of the dry zone of Blast FurnaceNo.I-BF of the Beijing Iron and Steel Company have been computed,from which a diagram for demar-cation of fluidization of coke and flooding of slag has been proposed.

Liquid-Gas Two-Phase Flow Simulation in Mold Filling Process with Level Set Method

Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect. The model is based on the combination of SOLA and Level Set Method. The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment. Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.

Affine particle-in-cell method for two-phase liquid simulation

Background The interaction of gas and liquid can produce many interesting phenomena,such as bubbles rising from the bottom of the liquid.The simulation of two-phase fluids is a challenging topic in computer graphics.To animate the interaction of a gas and liquid,MultiFLIP samples the two types of particles,and a Euler grid is used to track the interface of the liquid and gas.However,MultiFLIP uses the fluid implicit particle(FLIP)method to interpolate the velocities of particles into the Euler grid,which suffer from additional noise and instability.Methods To solve the problem caused by fluid implicit particles(FLIP),we present a novel velocity transport technique for two individual particles based on the affine particle-in-cell(APIC)method.First,we design a weighed coupling method for interpolating the velocities of liquid and gas particles to the Euler grid such that we can apply the APIC method to the simulation of a two-phase fluid.Second,we introduce a narrowband method to our system because MultiFLIP is a time-consuming approach owing to the large number of particles.Results Experiments show that our method is well integrated with the APIC method and provides a visually credible two-phase fluid animation.Conclusions The proposed method can successfully handle the simulation of a two phase fluid.

基于CFD的LNG储罐瞬时大风环境数值模拟研究

为分析不同大风风速下储罐迎风侧和背风侧压力分布情况和风场特点,采用k-ε湍流模型,建立液化天然气(LNG)储罐区瞬时大风环境模型。利用计算流体动力学(CFD)软件Fluent对LNG储罐风环境进行模拟,利用软件分析模块生成储罐压力云图和风速矢量图。研究结果表明:LNG储罐在迎风侧承受压力高于背风侧,压力呈现对称分布;储罐前端与两侧易产生气体漩涡,但总体速度曲线较为规律,无混乱复杂的速度曲线。研究结果可为LNG场站抗灾韧性增强、极端大风天气应急疏散路径规划等多方面提供理论参考,对LNG储罐区安全建设具有重要指导意义。

变容积密集烤房的CFD分析与试验研究

为保障密集烤房装烟密度,降低烘烤的能源消耗,研发一套适用于密集烤房的变容积系统。在完成变容积装置的设计后,基于CFD方法模拟分析装置与烟叶的不同距离对烤房内部气体分布均匀性的影响。通过烘烤试验获取实际烘烤数据,对模拟值加以验证。试验结果表明:当隔板与烟叶距离为0 cm时,流速不均匀系数Kv为0.40,温度不均匀系数Kt为0.41,距离为10 cm时,Kv=0.41,Kt=0.43;距离为20 cm时,Kv=0.42,Kt=0.49。装烟区9个测量点的温度模拟值与实测值基本吻合,误差在6%以内。变容积烤房在装烟量为一半时,相比常规烤房的燃料消耗可节约13.4%。研究结果表明:当隔板与烟叶距离为0 cm时烤房内部的气体分布最均匀;CFD模型与数值模拟结果具有可靠性;变容积装置具有较好的保温效果,可保证装烟密度,降低烤烟能耗。

基于CFD与风洞试验的边主梁涡振气动措施

为能够快捷且经济地完成开口类钝体桥梁断面涡振制振气动措施的选型,以一座边主梁叠合梁斜拉桥为背景,采用“CFD(computation fluid dynamics)数值模拟选型+风洞试验验证”的思路对其涡振制振气动措施选型进行研究.该桥原设计主梁断面在常遇风速下存在显著涡激振动,为完成气动措施的初步选型,采用CFD数值计算对原设计断面的流场进行模拟,通过研究原设计断面的旋涡脱落状态确定主要旋涡抑制对象,进而有针对性地模拟了3种气动措施(下中央稳定板、导流板与风嘴)对主要脱落旋涡的抑制作用,通过将各断面旋涡脱落状态与三分力系数进行对比分析,得到各断面涡振性能的相对优劣关系,并最终选取风嘴与下中央稳定板结合而成的组合气动措施进行风洞验证试验.试验结果表明:该组合气动措施能够有效抑制梁体在各风攻角下的涡激振动,且在+5°风攻角下,通过风洞试验得到的导流板、下中央稳定板、风嘴组合气动3种措施对原设计断面涡振振幅的减小作用依此递增,分别为2.7%、27.7%与87.4%,制振能力高低关系与数值模拟结果相一致;本次数值模拟结果符合预期要求,未来可针对不同类型桥梁断面进一步扩展数值模拟与风洞试验结果对比的数据集,以期更为准确、快捷地完成气动措施的选型.