Averaged Dynamics of Fluids near the Oscillating Interface in a Hele-Shaw Cell

The steady flow in a Hele-Shaw cell filled with fluids with a high viscosity contrast in the presence of fluid oscillations is experimentally studied.The control of oscillatory dynamics of multiphase systems with interfaces is a challenging technological problem.We consider miscible(water and glycerol)and immiscible(water and high-viscosity silicone oil PMS-1000)fluids under subsonic oscillations perpendicular to the interface.Observations show that the interface shape depends on the amplitude and frequency of oscillations.The interface is undisturbed only in the absence of oscillations.Under small amplitudes,the interface between water and glycerol widens due to mixing.When the critical amplitude is reached,the interface becomes unstable to the fingering instability:Aqueous fingers penetrate the high-viscosity glycerol and induce intensive mixing of miscible fluids and associated decay of the instability.After the disappearance of the fingers,the interface takes a U-shape in the central part of the cell.A similar effect is observed for immiscible fluids:The oscillating interface tends to bend to the side of a high-viscosity fluid.Again,when the critical amplitude is reached,the fingering instability arises at the convex interface.This paper focuses on the causes of bending of the initially undisturbed interface between miscible or immiscible fluids.For this purpose,we measure the steady flow velocity near the interface and in the bulk of a high-viscosity fluid using Particle Image Velocimetry(PIV).

Mechanisms of salt rejection at the ice-liquid interface during the freezing of pore fluids in the seasonal frozen soil area

Seasonal frozen soil accounts for about 53.50%of the land area in China.Frozen soil is a complex multiphase system where ice,water,soil,and air coexist.The distribution and migration of salts in frozen soil during soil freezing are notably different from those in unfrozen soil areas.However,little knowledge is available about the process and mechanisms of salt migration in frozen soil.This study explores the mechanisms of salt migration at the ice-liquid interface during the freezing of pore fluids through batch experiments.The results are as follows.The solute concentrations of liquid and solid phases at the ice-liquid interface(C*_(L),C*_(S))gradually increased at the initial stage of freezing and remained approximately constant at the middle stage.As the ice-liquid interface advanced toward the system boundary,the diffusion of the liquid phase was blocked but the ice phase continued rejecting salts.As a result,C*_(L)and C*_(S)rapidly increased at the final stage of freezing.The distribution characteristics of solutes in ice and the liquid phases before C*_(L)and C*_(S)became steady were mainly affected by the freezing temperature,initial concentrations,and particle-size distribution of media(quartz sand and kaolin).In detail,the lower the freezing temperature and the better the particle-size distribution of media,the higher the solute proportion in the ice phase at the initial stage of freezing.Meanwhile,the increase in concentration first promoted but then inhibited the increase of solutes in the ice phase.These results have insights and scientific significance for the tackling of climate change,the environmental protection of groundwater and soil,and infrastructure protection such as roads,among other things.

Dynamics of Low-Viscosity Liquids Interface in an Unevenly Rotating Vertical Layer

The behavior of two immiscible low-viscosity liquids differing in density and viscosity in a vertical flat layer undergoing modulated rotation is experimentally studied.The layer has a circular axisymmetric boundary.In the absence of modulation of the rotation speed,the interphase boundary has the shape of a short axisymmetric cylinder.A new effect has been discovered,under the influence of rotation speed modulation,the interface takes on a new dynamic equilibrium state.A more viscous liquid covers the end boundaries of the layer in the form of thin films,which have the shape of round spots of almost constant radius;with increasing amplitude of the velocity modulation,the wetting boundary expands.It is found that upon reaching the critical amplitude of oscillations,the film of a viscous liquid loses stability,and the outer edge of the wetting spot collapses and takes on a feathery structure.It is shown that this threshold is caused by the development of the Kelvin-Helmholtz oscillatory instability of the film.The spreading radius of a spot of light viscous liquid and its stability are studied depending on the rotation rate,amplitude,and frequency of rotation speed modulation.The discovered averaged effects are determined by different oscillatory interaction of fluids with the end-walls of the cell,due to different viscosities.The effect of films forming can find application in technological processes to intensify mass transfer at interphase boundaries.

NUMERICAL SIMULATION ON MULTI-FLUID INTERFACES

Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eulerian equations. The boundary conditions at interface of two fluids are determined by Ghost fluid method (GFM). The distributions of fluid parameters, such as pressure and density, were got at different time steps. The results show that the method presented in this paper can track the density discontinuity perfectly. Superior to previous results, the density discontinuity remains sharper. Also, the mixing of fluids can be greatly enhanced by setting disturbances along the initial fluid interfaces.

Monolithic approach to thermal fluid-structure interaction withnonconforming interfaces

This paper presents a monolithic approach to the thermal fluidstructure interaction （FSI） with nonconforming interfaces. The thermal viscous flow is governed by the Boussinesq approximation and the incompressible NavierStokes equations. The motion of the fluid domain is accounted for by an arbitrary LagrangianEulerian （ALE） strategy. A pseudosolid formulation is used to manage the deformation of the fluid do main. The structure is described by the geometrically nonlinear thermoelastic dynamics. An efficient data transfer strategy based on the Gauss points is proposed to guarantee the equilibrium of the stresses and heat along the interface. The resulting strongly coupled set of nonlinear equations for the fluid, solution procedure. A numerical example efficiency of the methodology. structure, and heat is solved by a monolithic is presented to demonstrate the robustness and

Sealing mechanism of magnetic fluids

Sealing is one of the most successful apphcatious of magnetic fluids. However, the sealing pressure difference is not satisfactory. This paper theoretically analyzes the mechanism of magnetic fluids sealing. Main factors that have significant effects on the sealing ability include viscous stress on the interracial surface, magnetic surface tension, and the shape of the interracial surface. The sealing pressure with magnetic fluids decreases with increase of rotational speed. Experiments were carried out to study the stability of the interface between magnetic fluids and water. It has been shown that stability of the interface will be damaged by washing of water when the relative flow between water and magnetic fluid becomes turbulent.

Progress in molecular-simulation-based research on the effects of interface-induced fluid microstructures on flow resistance

In modern chemical engineering processes, solid interface involvement is the most important component of process intensification techniques, such as nanoporous membrane separation and heterogeneous catalysis. The fundamental mechanism underlying interfacial transport remains incompletely understood given the complexity of heterogeneous interfacial molecular interactions and the high nonideality of the fluid involved. Thus, understanding the effects of interface-induced fluid microstructures on flow resistance is the first step in further understanding interfacial transport. Molecular simulation has become an indispensable method for the investigation of fluid microstructure and flow resistance. Here, we reviewed the recent research progress of our group and the latest relevant works to elucidate the contribution of interface-induced fluid microstructures to flow resistance.We specifically focused on water, ionic aqueous solutions, and alcohol–water mixtures given the ubiquity of these fluid systems in modern chemical engineering processes. We discussed the effects of the interfaceinduced hydrogen bond networks of water molecules, the ionic hydration of ionic aqueous solutions, and the spatial distributions of alcohol and alcohol–water mixtures on flow resistance on the basis of the distinctive characteristics of different fluid systems.

Physicochemical properties and structure of fluid at nano-/micro-interface:Progress in simulation and experimental study

In modern chemical engineering processes, the involvement of solid/fluid interface is the most important component of process intensification techniques, such as confined membrane separation and catalysis. In the review, we summarized the research progress of the latest theoretical and experimental works to elucidate the contribution of interface to the fluid properties and structures at nano-and micro-scale. We mainly focused on water, alcohol aqueous solution, and ionic liquids, because they are classical systems in interfacial science and/or widely involved in the industrialization process. Surface-induced fluids were observed in all reviewed systems and played a critical role in physicochemical properties and structures of outside fluid. It can even be regarded as a new interface, when the adsorption layer has a strong interaction with the solid surface. Finally, we proposed a perspective on scientific challenges in the modern chemical engineering processes and outlined future prospects.

Convective transport of nanoparticles in multi-layer fluid flow

Technologically, multi-layer fluid models are important in understanding fluid-fluid or fluid-nanoparticle interactions and their effects on flow and heat transfer characteristics. However, to the best of the authors＇ knowledge, little attention has been paid to the study of three-layer fluid models with nanofluids. Therefore, a three-layer fluid flow model with nanofluids is formulated in this paper. The governing coupled nonlinear differential equations of the problem are non-dimensionalized by using appropriate fundamental quantities. The resulting multi-point boundary value problem is solved numerically by quasi-linearization and Richardson＇s extrapolation with modified boundary conditions. The effects of the model parameters on the flow and heat transfer are obtained and analyzed. The results show that an increase in the nanoparticle concentration in the base fluid can modify the fluid-velocity at the interface of the two fluids and reduce the shear not only at the surface of the clear fluid but also at the interface between them. That is, nanofluids play a vital role in modifying the flow phenomena. Therefore, one can use nanofluids to obtain the desired qualities for the multi-fluid flow and heat transfer characteristics.

Polar interface and surface optical vibration spectra in multi-layer wurtzite quantum wires： transfer matrix method

The polar interface optical （IO） and surface optical （SO） phonon modes and the corresponding Froehlich electron phonon-interaction Hamiltonian in a freestanding multi-layer wurtzite cylindrical quantum wire （QWR） are derived and studied by employing the transfer matrix method in the dielectric continuum approximation and Loudon＇s uniaxial crystal model. A numerical calculation of a freestanding wurtzite GaN/AlN QWR is performed. The results reveal that for a relatively large azimuthal quantum number m or wave-number kz in the free z-direction, there exist two branches of IO phonon modes localized at the interface, and only one branch of SO mode localized at the surface in the system. The degenerating behaviours of the IO and SO phonon modes in the wurtzite QWR have also been clearly observed for a small kz or m. The limiting frequency properties of the IO and SO modes for large kz and m have been explained reasonably from the mathematical and physical viewpoints. The calculations of electron-phonon coupling functions show that the high-frequency IO phonon branch and SO mode play a more important role in the electron phonon interaction.

TRACKING AND RECONSTRUCTION METHODS FOR MOVING-INTERFACES

A simple,efficient and accurate high resolution method to tracking moving-interfaces-the characteristic integral-averaging finite volume method on unstructured meshes is proposed. And some numerical tests and evaluation of six main efficient methods for interface reconstruction are made. Through strict numerical simulation,their characters,advantages and shortcomings are compared,analyzed and commended in particular.

Influence of Drilling Fluid Components on Shear Strength at Cement-aquifuge Interface in Coalbed Methane Wells

The hydraulic fracturing is still an effective technology for the exploitation of coalbed methane （CBM）. However, after the hydraulic fracturing operation, the high water cut or sudden water flooding of CBM well usually occurs due to upward migration of bottom water, which is called water channeling （water inrush）. This problem has been severely limiting the hydraulic fracturing effect of CBM wells. Some studies show that the aquifuge and cement paste themselves will not crush under hydraulic fracturing pressure. Water channeling often occurs at cement- aquifuge interface （CAI）.

NUMERICAL STUDY OF TWO INTERFACE METHODS FOR COMPRESSIBLE MULTI-FLUID FLOWS

Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations written in a quasi conservative form, which is solved by a standard high resolution piecewise parabolic method (PPM) with multi fluid Riemann solver. The level set interface capturing method uses a narrow band ghost fluid method (GFM) with no numerical smearing. Several examples are presented and compared for one and two dimensions, which show the feasibility of the two methods applied to various multi fluid problems.

Transport diffuse interface model for simulation of solid-fluid interaction

For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. This is because the conservative variables in phasedensity equations include volume fractions. The degeneracy can be avoided by adding an artificial quantity of another material into the pure phase. However, nonphysical waves,such as shear waves in fluids, are introduced by the artificial treatment. In this paper,a transport diffuse interface model, which is able to treat zero/unity volume fractions, is presented for solid-fluid interaction. In the proposed model, a new formulation for phase densities is derived, which is unrelated to volume fractions. Consequently, the new model is able to handle zero/unity volume fractions, and nonphysical waves caused by artificial volume fractions are prevented. One-dimensional and two-dimensional numerical tests demonstrate that more accurate results can be obtained by the proposed model.

Acoustic Modeling of a 3-Layered Panel Incorporating Electro-Rheological/ Magneto-Rheological (EMR) Fluids

Applications of Electro-Rheological (ER) or Magneto-Rheological (MR) fluids as typical smart materials have been widely investigated over the past decades (since their introduction in 40’s). The special applications of these materials as a means of noise suppression are not yet investigated. Constrained Layer Damping (CLD) sheets can be realized by incorporating EMR (ER/MR) materials. In this way, a multilayered damping sheet is obtained with adaptive (tunable) stiffness and damping characteristics. These properties are easily changed in proportion to the electric (magnetic) field applied upon the EMR layer. This notion has been introduced for semi-active vibration control problems. Herein, such panels incorporating EMR material are proposed for adaptive acoustic treatments. Modeling (simulation) of a 3-layered panel with the middle layer being EMR with adjustable properties is carried out in this paper. The tunability of transmission/absorption characteristics of these composite sheets enables us making smart panels for adaptive noise and acoustic treatments. An adaptive performance can be achieved via changing the properties of such panels, on line, according to some sensor outputs. The main objective is to develop proper models to predict the Transmission Loss (TL) of such panels. Also, the TL of this panel is compared with the middle layer of a Newtonian fluid.

Application of Fluid Phase Recovery Method in the Early Development of BZ26-6 Volatile Oil Field in Bohai Sea

BZ26-6 Oilfield is a kind of deep metamorphic rock buried-hill volatile oilfield in Bohai Sea, China. Its early development plan is restricted due to the simultaneous production of oil and gas in large sections of reservoirs, unclear understanding of formation fluid properties and uncertainty of gas-oil interface. Through theoretical research on phase recovery and experimental analysis of crude oil phase characteristics in the original formation, characteristic parameters of the equilibrium condensate gas fluid are restored and calculated. Through the superimposed phase diagram of volatile oil and condensate gas, BZ26-6 Oilfield is determined to be a volatile oil reservoir with a condensate gas cap, with formation pressure and saturation pressure of 36.1 MPa, respectively. Based on the research results of oil-gas phase behavior characteristics, the thermodynamic equations and equation of state are jointly used to solve the problem, and the content change curves of each component at different depths are drawn. Combined with the sensitivity analysis of numerical simulation, the gas-oil interface is determined to be -3726 m above sea level. The fluid phase analysis software, Fluidmodeler, is used to simulate volatile oil degassing and condensate gas separation experiments. In combination with oil and gas production data obtained through the production test, the specific oil recovery index and the specific gas recovery index are determined to be 0.408 m3/(MPa·d·m) and 1195 m3/(MPa·d·m), respectively. And the reasonable production capacity prediction is conducted on the early development of BZ26-6 Oilfield. The research results can provide a theoretical basis for the efficient development of similar complex oil and gas reservoirs.

亲油水泥浆界面封隔性能评价研究

页岩气井油基钻井液条件下固井时,套管与井壁表面附着油膜或油基钻井液易导致水泥环界面封隔能力下降,形成窜流通道,影响后期增产改造作业。虽然前置液能有效提高界面润湿反转,但受用量、冲洗效率等因素限制,界面仍会出现油膜附着情况。为此,在水泥浆中加入亲油表面活性剂,制备形成亲油水泥浆,赋予水泥环亲油能力。采用接触角、剪切胶结强度、界面水力封隔测试等评价了亲油水泥石的亲油能力及界面封隔效果。研究发现,非极性溶剂在亲油水泥石表面接触角远低于常规水泥石表面,具备良好的亲油性能;亲油水泥石与含白油、油基钻井液的套管、页岩岩芯胶结后,一、二界面抗流体窜流压力分别提升500%和400%,胶结强度分别提高205%和122%;亲油表面活性剂的加入不会对水泥水化程度、水泥浆工程性能及水泥石力学性能产生负面影响。结果表明,亲油表面活性剂掺入后可有效提高水泥环与含油界面的封隔能力,具备提高油基钻井液条件下水泥环套管地层界面胶结与封隔性能的潜力。

深水气田长距离混输模式下水相运动规律研究

陵水25-1深水气田混输管道距离长、管道爬坡段坡度大,初始投产阶段混输流体水相到达平台的时间长且预测困难,给管道安全运行带来风险。基于多相流计算模型,运用LedaFlow软件,对陵水25-1深水气田初始投产阶段的水相运动规律进行研究,并深入分析含水率、油量、气量、油水相间摩擦系数等因素对水相运行的影响。结果表明:初始投产阶段水相在爬坡段分离和积聚,形成长时间的停滞,管道坡度越大,水相波动越剧烈;含水率越高的物流,水相在爬坡段聚集时间越短;气量增大则会降低水相滞液程度,提高水相运行速度,有助于水相爬坡;油量变化对水相运行速度的影响较小,油水相间摩擦系数对爬坡段水相运动有直接影响,油水相间摩擦系数修正拟合,对实际管道运行有重要指导意义。研究结果可为深水气田开发研究提供借鉴。

锦州25-1南X井区潜山油气藏精细评价技术与实践

锦州25-1南是渤海首个规模性潜山油气藏开发油田,目前已累计投产36口井,是渤海最大的生产潜山裂缝性油藏。随着勘探区域的不断扩大及勘探井网的多次加密,在精确判定潜山内幕带优势储层发育位置及流体界面等问题上仍存在诸多困扰。依托该区域已钻井及围区新增勘探井为数据体,通过对已钻井的随钻气测、工程参数及地化录井等数据进行深入分析,引入储层评价系数,建立参数模型及评价图版,实现了储层快速评价及潜山流体类型的识别。该方法成功实现了锦州25-1南太古界潜山评价目标,X井区油底加深136 m,同时该方法为后续渤海油田锦州25-3区块及渤中26-6区块太古界潜山综合高效开发评价提供关键技术支持。

Mathematical Modelling of Particle Movement Ahead of the Solid-liquid Interface in Continuous Casting

Whether the particle will be trapped by the solid-liquid interface or not is dependent on its moving behavior ahead of the interface, so a mathematical model has been developed to investigate the movement of the particle ahead of the solid-liquid interface. Based on the theory for the boundary layer, the fluid velocity field near the solid-liquid interface was obtained, and the trajectories of particles were calculated by the equations of motion for particles. In this model, the drag force, the added mass force, the buoyance force, the gravitational force, the Saffman force and the Basset history force are considered. The results show that the behavior of the particle ahead of the solid-liquid interface is affected by the physical property of the particle and fluid flow. And in the continuous casting process, if it moves in the stream directed upward or downward near vertical solid-liquid interface or in the horizontal flow under the solid-liquid interface, the particle with the diameter from 5 um to 60um can reach the solid-liquid interface. But if it moves in horizontal flow above the solid-liquid interface, only the particle with the diameter from 5 um to 10 um can reach the solid-liquid interface.