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.

A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

Effect of microalloying on wettability and interface characteristics of Zr-based bulk metallic glasses with W substrate

The infiltration casting method is widely employed for the preparation of ex-situ composite materials.However,the production of composite materials using this method must necessitates a comprehensive understanding of the wettability and interface characteristics between the reinforcing phase and the bulk metallic glasses(BMGs).This work optimized the composition of Zr-based BMGs through microalloying methods,resulting in a new set of Zr-based BMGs with excellent glass-forming ability.Wetting experiments between the Zr-based BMGs melts and W substrates were conducted using the traditional sessile drop method,and the interfaces were characterized utilizing a scanning electron microscope(SEM)equipped with energy dispersive X-ray spectroscopy(EDS).The work demonstrates that the microalloying method substantially enhances the wettability of the Zr-based BMGs melt.Additionally,the incorporation of Nb element impedes the formation of W-Zr phases,but the introduction of Nb element does not alter the extent of interdiffusion between the constituent elements of the amorphous matrix and W element,indicating that the influence of Nb element on the diffusion of individual elements is minute.

Comparison of the interface reaction behaviors of CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) solid-state systems based on the diffusion couple method

The formation mechanism of calcium vanadate and manganese vanadate and the difference between calcium and manganese in the reaction with vanadium are basic issues in the calcification roasting and manganese roasting process with vanadium slag.In this work,CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples were prepared and roasted for different time periods to illustrate and compare the diffusion reaction mechanisms.Then,the changes in the diffusion product and diffusion coefficient were investigated and calculated based on scanning electron microscopy (SEM) with energy dispersive X-ray spectroscopy (EDS) analysis.Results show that with the extension of the roasting time,the diffusion reaction gradually proceeds among the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples.The regional boundaries of calcium and vanadium are easily identifiable for the CaO–V_(2)O_(5) diffusion couple.Meanwhile,for the MnO_(2)–V_(2)O_(5) diffusion couple,MnO_(2) gradually decomposes to form Mn_(2)O_(3),and vanadium diffuses into the interior of Mn_(2)O_(3).Only a part of vanadium combines with manganese to form the diffusion production layer.CaV_(2)O_(6) and MnV_(2)O_(6) are the interfacial reaction products of the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples,respectively,whose thicknesses are 39.85 and 32.13μm when roasted for 16 h.After 16 h,both diffusion couples reach the reaction equilibrium due to the limitation of diffusion.The diffusion coefficient of the CaO–V_(2)O_(5) diffusion couple is higher than that of the MnO_(2)–V_(2)O_(5) diffusion couple for the same roasting time,and the diffusion reaction between vanadium and calcium is easier than that between vanadium and manganese.

The Behavior of Amphiphile at Oil-Water Interface by Monte Carlo Simulation

A novel simple two-dimensional square-lattice model of amphiphile at oil-water interface is developed,in which oil and water act as solvent and occupy empty sites and amphiphile occupies chains of sites. In this model, the oil-water interface is fixed, And amphiphile molecules will be enriched at the oil-water interface. The interfacial concentration of amphiphile calculated by Monte Carlo method shows that it is easier for the hydrophilic-hydrophobic balanced amphiphile to stay at the interface. And the adsorption of amphiphile increases with the increase of amphiphile concentration and the decrease with temperature.

Synergistic Effect of Amphiphiles at Oil-Water Interface: By Monte Carlo Simulation

Amphiphile-oil-water system is complicated. The real behavior of amphiphile in the interface is still undnown despite that this behavior is very important in determining the stability of emulsion system. In this paper, the interface properties of amphiphile at oil-water interface were investigated by a square-lattice model Monte Carlo simulation method. The synergistic effect was found for hydrophobic and hydrophilic amphiphile mixture systems; and the synergistic effect disappears or was weakened as the amphiphile at the interface region became dilute with the increasing of temperature.

Formation of Uniform Oil-Soluble Fe_3O_4 Nanoparticles via Oil-Water Interface System

Nowadays,the novel oil water interface method has attracted a considerable attention owing to the advantages of mild reaction conditions,simple operation,low cost,and high efficiency.In this paper,uniform oil-soluble Fe_3O_4 nanoparticles(NPs) were synthesized by oil-water interface method from mixing iron tristearate of 0.067mol/L in cyclohexane with ferrous sulfate in water.The as-prepared products were characterized by X-ray diffraction(XRD),transmission electron microscopy(TEM),vibrating sample magnetometer(VSM),Fourier transform infrared spectroscopy(FT-IR) and thermogravimetric analyzer(TGA).TEM images and XRD profiles showed that the size of the oil-soluble products ranged in 1.7-6.9 nm.VSM indicated that the Fe_3O_4 NPs were superparamagetic.FT-IR and TGA proved that oleic acid was combined to the surface of Fe_3O_4 NPs closely.TEM images and XRD profiles revealed that the most suitable reaction concentration of NH_3·H_2O,oleic acid/water in volume,reaction temperature and reaction time were 4.5 mol/L,50:1 000,80℃ and 6 h,respectively.The formation mechanism of the nearly monodispersed Fe_3O_4 NPs was that the preformed Fe_3O_4 nuclei were capped by oleic acid as early as the nucleation occurred in oil-water interface and subsequently entered into oil phase to stop growing.

Growth interface of CdMnTe crystal by traveling heater method

The growth interfaces of CdMnTe（CMT） crystals grown by traveling heater method（THM） were studied. Two types of polycrystalline CMT feed ingots synthesized in a traditional rocking furnace and vertical Bridgman（VB） furnace were adopted in THM growth, and the effects of the polycrystalline feed on the growth interface were revealed. The morphology of the growth interface of CMT crystal（CMT2） grown from the feed by vertical Bridgman was smoother with lower curvature compared with that of CMT crystal（CMT1） from the feed by rocking furnace. The radial Mn composition and Te inclusion distribution of the CMT wafers were analyzed and correlated to the growth interface. The Mn segregation along the radial direction and Te inclusion density of CMT2 were lower than those of CMT1. The VB method synthesized polycrystalline feed could improve the growth interface morphology, which is beneficial for decreasing the Te inclusions and Mn segregation in CMT wafers.

Two-grid methods for semi-linear elliptic interface problems by immersed finite element methods

In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.

Studying the mechanical properties of the soil-root interface using the pullout test method

It is important to quantify the effect of the root diameter, the embedment length of the root and load speed on the soil-root interface mechanical properties for studying the root anchorage. The soilroot interface mechanical properties can be obtained through the pullout force and root slippage curve(F-S curve). About 120 Pinus tabulaeformis single roots whose diameters ranged from 1 mm to 10 mm divided into 6 groups based on different root embedment length(50 mm, 100 mm and 150 mm) and different load velocity(10 mm·min^(-1), 50 mm·min^(-1), 100 mm·min^(-1) and 300 mm·min^(-1)) were investigated using the pullout method. This study aims to explore the mechanical properties of the soil-root interface in the real conditions using the pullout test method. The results showed two kinds of pullout test failure modes during the experimental process: breakage failure and pullout failure. The results showed that the roots were easier to be broken when the root diameter was smaller or the loading speed was larger. The relationship between the maximum anchorage force and root diameter was linear and the linearly dependent coefficient(R^2) was larger than 0.85. The anchorage force increased with the root embedment length. An increase of 10%^(-1)5% for the maximumanchorage force was found when load speed increased from 10 to 300 mm.min^(-1). The mean peak slippage of the root was from 13.81 to 35.79 mm when the load velocity varied from 10 to 300 mm.min^(-1). The study will be helpful for the design of slopes reinforced by vegetation and in predicting risk of uprooting of trees, and will have practical benefits for understanding the mechanism of landslide.

Experimental investigation of interface curing stresses between PMMA and composite using digital speckle correlation method

This paper studies the interface curing stresses between polymethyl methacrylate (PMMA) and composite by means of digital speckle correlation method (DSCM).A new method by combining DSCM with the marker points is developed to measure the interface curing stresses,and the measurement principle is introduced.The interface curing stresses between PMMA and composite with different curing bonding conditions are measured and analyzed,this indicates that the residual stress for furnace heating and furnace cooling is the smallest.Finally,the measurement error is discussed by means of finite element method,the influences of glass microsphere between adhesive and PMMA can be ignored.

LOCAL AND PARALLEL FINITE ELEMENT METHOD FOR THE MIXED NAVIER-STOKES/DARCY MODEL WITH BEAVERS-JOSEPH INTERFACE CONDITIONS

In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.

Immersed Interface Finite Element Methods for Elasticity Interface Problems with Non-Homogeneous Jump Conditions

In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.

P_1-nonconforming triangular finite element method for elliptic and parabolic interface problems

The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.

A NEW TREATMENT FOR INTERFACE OF METAL AND DIELECTRIC IN FDTD METHOD

A new treatment for interface of metal and dielectric in finite-difference time domain (FDTD) method was presented, further how and why this treatment should be as it is was explained. The FDTD method was developed to compute the low-frequency electromagnetic fields in the metal instead of the traditional high-frequency ones as it was originally proposed. The computational results agreed well with the experimental ones.

Modeling Methods and Test Verification of Root Insert Contact Interface for Wind Turbine Blade

Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the metal contact interface,LMOM is proposed to analyze the load path and stress distribution characteristics,while GMM is used to calculate and analyze the stress distribution characteristics of the resin layer established between the bushing and composite layers of root insert.To validate the GMM,a tension test is carried out.The result successfully shows that the shear strain expresses a similar strain distribution tendency with the GMM′s results.

AN EMBEDDED BOUNDARY METHOD FOR ELLIPTIC AND PARABOLIC PROBLEMS WITH INTERFACES AND APPLICATION TO MULTI-MATERIAL SYSTEMS WITH PHASE TRANSITIONS

The embedded boundary method for solving elliptic and parabolic problems in geometrically complex domains using Cartesian meshes by Johansen and Colella （1998, J. Comput. Phys. 147, 60） has been extended for elliptic and parabolic problems with interior boundaries or interfaces of discontinuities of material properties or solutions. Second order accuracy is achieved in space and time for both stationary and moving interface problems. The method is conservative for elliptic and parabolic problems with fixed interfaces. Based on this method, a front tracking algorithm for the Stefan problem has been developed. The accuracy of the method is measured through comparison with exact solution to a two-dimensional Stefan problem. The algorithm has been used for the study of melting and solidification problems.

Effect of SDBS on interfacial electron transfer at the liquid/liquid interface by thin layer method

The effect of an adsorbed anionic surfactant sodium dodecyl benzene sulfonate （SDBS） on electron transfer （ET） reaction between ferricyanide aqueous solution and decamethylferrocene （DMFc） located on the adjacent organic phase was investigated for the first time by thin layer method. The adsorption of SDBS at the interface resulted in a decay in the cathodic plateau current of bimolecular reaction with increasing concentrations of SDBS in aqueous phase. However, the rate constant of electron transfer （ket） increased monotonically as the SDBS concentrations increased from 0 to 200 p, moFL. The experimental results showed that SDBS formed patches on the interface and influenced the structure of electrical double layer. 2009 Xiao Quan Lu. Published by Elsevier B.V. on behalf of Chinese Chemical Society. All rights reserved.

Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients

The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.

Immersed Interface Method for Fokker-Planck Equation with Discontinuous Drift

The Immersed Interface Method (IIM) is derived to solve the corresponding Fokker-Planck equation of Brownian motion with pure dry friction, which is one of the simplest models of piecewise-smooth stochastic systems. The IIM is capable of treating a discontinuity in the drift of Fokker-Planck equation and it is readily extended to the dry and viscous friction model. Analytic results of the considered model are used to confirm the effectiveness and design accuracy of the scheme.