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.

Numerical Simulation and Dynamical Analysis for Low Salinity Water Lens in the Expansion Area of the Changjiang Diluted Water

The low salinity water lenses（LSWLes） in the expansion area of the Changjiang diluted water（CDW） exist in a certain period of time in some years. The impact of realistic river runoff, ocean currents and weather conditions need to be taken into account in the dynamical analysis of LSWL, which is in need of research. In this paper, the POM-σ-z model is used to set up the numerical model for the expansion of the CDW. Then LSWL in summer 1977 is simulated, and its dynamic mechanism driven by wind, tide, river runoff and the Taiwan Warm Current is also analyzed. The simulated results indicate that the isolated LSWL detaches itself from the CDW near the river mouth, and then moves towards the northeast region outside the Changjiang Estuary. Its maintaining period is from July 26 to August 11. Its formation and development is mainly driven by two factors. One is the strong southeasterly wind lasting for ten days. The other is the vertical tidal mixing during the transition from neap tide to spring tide.

Numerical simulation of flow and heat transfer in a random packed bed

Random packed beds have more complex interior structure than structured beds and are widely used in industry and engineering. CFD simulation was carried out to investigate and analyze the local flow and heat transfer in a 120-sphere random packed bed. 3D Navier-Stokes equation was solved with a finite volume formulation based on the Chimera meshing technique. Investigation was focused on low Reynolds number flow （Re=4.6-56.2）, which typically occurs in packed bed reactors in bio-chemical fields. Detailed temperature field information was obtained. Inhomogeneity of flow and heat transfer due to the non-uniform distribution of void fraction was discussed and analyzed.

The Erosion Effect of Kapton Film in a Ground-based Atomic Oxygen Irradiation Simulator

In order to investigate the effect of space environmental factors on spacecraft materials, a ground-based simulation facility for space atomic oxygen（AO） irradiation was developed in our laboratory. Some Kapton film samples were subjected to AO beam generated by this facility. The Kapton films before and after AO exposure were analyzed comparatively using optical microscopy, scanning electronic microscopy, atomic force microscopy, high-precision microbalance, and X-ray photoelectron spectroscopy. The experimental results indicate that the transmittance of Kapton film will be reduced by AO irradiation notably, and its color deepens from pale yellow to brown. Surface roughness of the AO-treated sample is already increased obviously after AO irradiation for 5 hours, and exhibits a flannel-like appearance after 15 hours’ exposure in AO beam. The imide rings and benzene rings in kapton molecule are partially decomposed, and some new bonds form during AO irradiation. The mass loss of kapton film increases linearly with the increase of AO fluence, which is resulted from the formation of volatile products, such as CO, CO2 and NOx. The breakage in structure and degradation in properties of AO-treated Kapton film can be attributed to the integrated effect ofimpaction and oxidization of AO beam. The test results agree well with the space flight experimental data.

Numerical study of separation on the trailing edge of a symmetrical airfoil at a low Reynolds number

This study focuses on the trailing-edge separation of a symmetrical airfoil at a low Rey-nolds number. Finite volume method is adopted to solve the unsteady Reynolds-averaged Navier-Stokes （RANS） equation. Flow of the symmetrical airfoil SD8020 at a low Reynolds number has been simulated. Laminar separation bubble in the flow field of the airfoil is observed and process of unsteady bubble burst and vortex shedding from airfoil surfaces is investigated. The time-dependent lift coefficient is characteristic of periodic fluctuations and the lift curve varies nonlinearly with the attack of angle. Laminar separation occurs on both surfaces of airfoil at small angles of attack. With the increase of angle of attack, laminar separation occurs and then reattaches near the trailing edge on the upper surface of airfoil, which forms laminar separation bubble. When the attack of angle reaches certain value, the laminar separation bubble is unstable and produces two kinds of large scale vortex, i.e. primary vortex and secondary vortex. The periodic processes that include secondary vortex production, motion of secondary vortex and vortex shedding cause fluctuation of the lift coefficient. The periodic time varies with attack of angle. The secondary vortex is relatively stronger than the primary vortex, which means its influence is relatively stronger than the primary vortex.

Particle-in-Cell with Monte CarloCollisions Gun Code Simulations of a Surface-Conversion H^(−)Ion Source

We present an extended update on the status of a particle-in-cellwithMonte Carlo collisions(PIC-MCC)gun code developed at LosAlamos for the study of surfaceconverter H−ion sources.The programis fully kinetic.Some of the program’s features include:solution of arbitrary electrostatic and magnetostatic fields in an axisymmetric(r,z)geometry to describe the self-consistent time evolution of a plasma;simulation of a multi-species(e^(−),H^(+),H^(+)_(2),H^(+)_(3),H^(−))plasma discharge from a neutral hydrogen gas and filament-originated seed electrons;full 2-dimensional(r,z)3-velocity(v_(r),v_(z),vφ)dynamics for all species;detailed collision physics between charged particles and neutrals and the ability to represent multiple smooth(not stair-stepped)electrodes of arbitrary shape and voltage whose surfaces may be secondary-particle emitters(H^(−)and e^(−)).The status of this development is discussed in terms of its physics content and current implementation details.