Fully coupled fluid-solid productivity numerical simulation of multistage fractured horizontal well in tight oil reservoirs

A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.

A semi-analytical model for coupled flow in stress-sensitive multi-scale shale reservoirs with fractal characteristics

A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.

Coupled multiphysical model for investigation of influence factors in the application of microbially induced calcite precipitation

The study presents a comprehensive coupled thermo-bio-chemo-hydraulic(T-BCH)modeling framework for stabilizing soils using microbially induced calcite precipitation(MICP).The numerical model considers relevant multiphysics involved in MICP,such as bacterial ureolytic activities,biochemical reactions,multiphase and multicomponent transport,and alteration of the porosity and permeability.The model incorporates multiphysical coupling effects through well-established constitutive relations that connect parameters and variables from different physical fields.It was implemented in the open-source finite element code OpenGeoSys(OGS),and a semi-staggered solution strategy was designed to solve the couplings,allowing for flexible model settings.Therefore,the developed model can be easily adapted to simulate MICP applications in different scenarios.The numerical model was employed to analyze the effect of various factors,including temperature,injection strategies,and application scales.Besides,a TBCH modeling study was conducted on the laboratory-scale domain to analyze the effects of temperature on urease activity and precipitated calcium carbonate.To understand the scale dependency of MICP treatment,a large-scale heterogeneous domain was subjected to variable biochemical injection strategies.The simulations conducted at the field-scale guided the selection of an injection strategy to achieve the desired type and amount of precipitation.Additionally,the study emphasized the potential of numerical models as reliable tools for optimizing future developments in field-scale MICP treatment.The present study demonstrates the potential of this numerical framework for designing and optimizing the MICP applications in laboratory-,prototype-,and field-scale scenarios.

Nonlinearly induced entanglement in dissipatively coupled optomechanical system

Nonlinearly induced steady-state photon–phonon entanglement of a dissipative coupled system is studied in the bistable regime. Quantum dynamical characteristics are analysed by solving the mean-field and fluctuation equations of the system. It is shown that dissipative coupling can induce bistable behaviour for the effective dissipation of the system.Under suitable parameters, one of the steady states significantly reduces the dissipative effect of the system. Consequently,a larger steady-state entanglement can be achieved compared to linear dynamics. Furthermore, the experimental feasibility of the parameters is analysed. Our results provide a new perspective for the implementation of steady-state optomechanical entanglement.

Seismic performance evaluation of hybrid coupled shear wall system with shear and flexural fuse-type steel coupling beams

Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically evaluate the seismic behavior of flexural and shear fuse beams situated in short-,medium-and high-rise RC buildings that have HCSWs.Three building groups hypothetically located in a high seismic hazard zone were studied.A series of 2D nonlinear time history analyses was accomplished in OpenSees,using the ground motion records scaled at the design basis earthquake level.It was found that the effectiveness of fuses in HCSWs depends on various factors such as size and scale of the building,allowable rotation value,inter-story drift ratio,residual drift quantity,energy dissipation value of the fuses,etc.The results show that shear fuses better meet the requirements of rotations and drifts.In contrast,flexural fuses dissipate more energy,but their sectional stiffness should increase to meet other requirements.It was concluded that adoption of proper fuses depends on the overall scale of the building and on how associated factors are considered.

Effects of asymmetric coupling and boundary on the dynamic behaviors of a random nearest neighbor coupled system

This study investigates the dynamical behaviors of nearest neighbor asymmetric coupled systems in a confined space.First, the study derivative analytical stability and synchronization conditions for the asymmetrically coupled system in an unconfined space, which are then validated through numerical simulations. Simulation results show that asymmetric coupling has a significant impact on synchronization conditions. Moreover, it is observed that irrespective of whether the system is confined, an increase in coupling asymmetry leads to a hastened synchronization pace. Additionally, the study examines the effects of boundaries on the system's collective behaviors via numerical experiments. The presence of boundaries ensures the system's stability and synchronization, and reducing these boundaries can expedite the synchronization process and amplify its effects. Finally, the study reveals that the system's output amplitude exhibits stochastic resonance as the confined boundary size increases.

Customized modulation on plasma uniformity by non-uniform magnetic field in capacitively coupled plasma

A two-dimensional fluid model based on COMSOL Multiphysics is developed to investigate the modulation of static magnetic field on plasma homogeneity in a capacitively coupled plasma(CCP)chamber. To generate a static magnetic field, direct current is applied to a circular coil located at the top of the chamber. By adjusting the magnetic field's configuration, which is done by altering the coil current and position, both the plasma uniformity and density can be significantly modulated. In the absence of the magnetic field, the plasma density exhibits an inhomogeneous distribution characterized by higher values at the plasma edge and lower values at the center. The introduction of a magnetic field generated by coils results in a significant increase in electron density near the coils. Furthermore, an increase in the sets of coils improves the uniformity of the plasma. By flexibly adjusting the positions of the coils and the applied current,a substantial enhancement in overall uniformity can be achieved. These findings demonstrate the feasibility of using this method for achieving uniform plasma densities in industrial applications.

A Coupled Thermomechanical Crack Propagation Behavior of Brittle Materials by Peridynamic Differential Operator

This study proposes a comprehensive,coupled thermomechanical model that replaces local spatial derivatives in classical differential thermomechanical equations with nonlocal integral forms derived from the peridynamic differential operator(PDDO),eliminating the need for calibration procedures.The model employs a multi-rate explicit time integration scheme to handle varying time scales in multi-physics systems.Through simulations conducted on granite and ceramic materials,this model demonstrates its effectiveness.It successfully simulates thermal damage behavior in granite arising from incompatible mineral expansion and accurately calculates thermal crack propagation in ceramic slabs during quenching.To account for material heterogeneity,the model utilizes the Shuffle algorithm andWeibull distribution,yielding results that align with numerical simulations and experimental observations.This coupled thermomechanical model shows great promise for analyzing intricate thermomechanical phenomena in brittle materials.

Coexisting fast–slow dendritic traveling waves in a 3D-array electric field coupled neuronal network

Coexistence of fast and slow traveling waves without synaptic transmission has been found in hhhippocampal tissues,which is closely related to both normal brain activity and abnormal neural activity such as epileptic discharge. However, the propagation mechanism behind this coexistence phenomenon remains unclear. In this paper, a three-dimensional electric field coupled hippocampal neural network is established to investigate generation of coexisting spontaneous fast and slow traveling waves. This model captures two types of dendritic traveling waves propagating in both transverse and longitude directions: the N-methyl-D-aspartate(NMDA)-dependent wave with a speed of about 0.1 m/s and the Ca-dependent wave with a speed of about 0.009 m/s. These traveling waves are synaptic-independent and could be conducted only by the electric fields generated by neighboring neurons, which are basically consistent with the in vitro data measured experiments. It is also found that the slow Ca wave could trigger generation of fast NMDA waves in the propagation path of slow waves whereas fast NMDA waves cannot affect the propagation of slow Ca waves. These results suggest that dendritic Ca waves could acted as the source of the coexistence fast and slow waves. Furthermore, we also confirm the impact of cellular spacing heterogeneity on the onset of coexisting fast and slow waves. The local region with decreasing distances among neighbor neurons is more liable to promote the onset of spontaneous slow waves which, as sources, excite propagation of fast waves. These modeling studies provide possible biophysical mechanisms underlying the neural dynamics of spontaneous traveling waves in brain tissues.

Multi-soliton solutions of coupled Lakshmanan-Porsezian-Daniel equations with variable coefficients under nonzero boundary conditions

This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.

A mechanically coupled MEMS filter with high-Q width extensional mode resonators

This work presents a novel radio frequency(RF)narrowband Si micro-electro-mechanical systems(MEMS)filter based on capacitively transduced slotted width extensional mode(WEM)resonators.The flexibility of the plate leads to multiple modes near the target frequency.The high Q-factor resonators of around 100000 enable narrow bandwidth filters with small size and simplified design.The 1-wavelength and 2-wavelength WEMs were first developed as a pair of coupled modes to form a passband.To reduce bandwidth,two plates are coupled with aλ-length coupling beam.The 79.69 MHz coupled plate filter(CPF)achieved a narrow bandwidth of 8.8 kHz,corresponding to a tiny 0.011%.The CPF exhibits an impressive 34.84 dB stopband rejection and 7.82 dB insertion loss with near-zero passband ripple.In summary,the RF MEMS filter presented in this work shows promising potential for application in RF transceiver front-ends.

Coupled effects of climate change and human activities on vegetation dynamics in the Southwestern Alpine Canyon Region of China

The driving effects of climate change and human activities on vegetation change have always been a focal point of research.However,the coupling mechanisms of these driving factors across different temporal and spatial scales remain controversial.The Southwestern Alpine Canyon Region of China(SACR),as an ecologically fragile area,is highly sensitive to the impacts of climate change and human activities.This study constructed a vegetation cover dataset for the SACR based on the Enhanced Vegetation Index(EVI)from 2000 to 2020.Spatial autocorrelation,Theil-Sen trend,and Mann-Kendall tests were used to analyze the spatiotemporal characteristics of vegetation cover changes.The main drivers of spatial heterogeneity in vegetation cover were identified using the optimal parameter geographic detector,and an improved residual analysis model was employed to quantify the relative contributions of climate change and human activities to interannual vegetation cover changes.The main findings are as follows:Spatially,vegetation cover exceeds 60%in most areas,especially in the southern part of the study area.However,the border area between Linzhi and Changdu exhibits lower vegetation cover.Climate factors are the primary drivers of spatial heterogeneity in vegetation cover,with temperature having the most significant influence,as indicated by its q-value,which far exceeds that of other factors.Additionally,the interaction q-value between the two factors significantly increases,showing a relationship of bivariate enhancement and nonlinear enhancement.In terms of temporal changes,vegetation cover shows an overall improving trend from 2000 to 2020,with significant increases observed in 68.93%of the study area.Among these,human activities are the main factors driving vegetation cover change,with a relative contribution rate of 41.31%,while climate change and residual factors contribute 35.66%and 23.53%,respectively.By thoroughly exploring the coupled mechanisms of vegetation change,this study provides important references for the sustainable management and conservation of the vegetation ecosystem in the SACR.

Abnormal transition of the electron energy distribution with excitation of the second harmonic in low-pressure radio-frequency capacitively coupled plasmas

The self-excited second harmonic in radio-frequency capacitively coupled plasma was significantly enhanced by adjusting the external variable capacitor.At a lower pressure of 3 Pa,the excitation of the second harmonic caused an abnormal transition of the electron energy probability function,resulting in abrupt changes in the electron density and temperature.Such changes in the electron energy probability function as well as the electron density and temperature were not observed at the higher pressure of 16 Pa under similar harmonic changes.The phenomena are related to the influence of the second harmonic on stochastic heating,which is determined by both amplitude and the relative phase of the harmonics.The results suggest that the self-excited high-order harmonics must be considered in practical applications of lowpressure radio-frequency capacitively coupled plasmas.

Numerical simulation study on the mold strength of magnetic mold casting based on a coupled electromagnetic-structural method

The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.

Contrasts of bimodal tropical instability waves(TIWs)-induced wind stress perturbations in the Pacific Ocean among observations,ocean models,and coupled climate models

The coupling between wind stress perturbations and sea surface temperature(SST)perturbations induced by tropical instability waves(TIWs)in the Pacific Ocean has been revealed previously and proven crucial to both the atmosphere and ocean.However,an overlooked fact by previous studies is that the loosely defined“TIWs”actually consist of two modes,including the Yanai wave-based TIW on the equator(hereafter eTIW)and the Rossby wave-based TIW off the equator(hereafter vTIW).Hence,the individual feedbacks of the wind stress to the bimodal TIWs remain unexplored.In this study,individual coupling relationships are established for both eTIW and v TIW,including the relationship between the TIW-induced SST perturbations and two components of wind stress perturbations,and the relationship between the TIW-induced wind stress perturbation divergence(curl)and the downwind(crosswind)TIW-induced SST gradients.Results show that,due to different distributions of eTIW and vTIW,the coupling strength induced by the eTIW is stronger on the equator,and that by the vTIW is stronger off the equator.The results of any of eTIW and vTIW are higher than those of the loosely defined TIWs.We further investigated how well the coupling relationships remained in several widely recognized oceanic general circulation models and fully coupled climate models.However,the coupling relationships cannot be well represented in most numerical models.Finally,we confirmed that higher resolution usually corresponds to more accurate simulation.Therefore,the coupling models established in this study are complementary to previous research and can be used to refine the oceanic and coupled climate models.

Multi-Physics Coupled Acoustic-Mechanics Analysis and Synergetic Optimization for a Twin-Fluid Atomization Nozzle

Fine particulate matter produced during the rapid industrialization over the past decades can cause significant harm to human health.Twin-fluid atomization technology is an effective means of controlling fine particulate matter pollution.In this paper,the influences of the main parameters on the droplet size,effective atomization range and sound pressure level(SPL)of a twin-fluid nozzle(TFN)are investigated,and in order to improve the atomization performance,a multi-objective synergetic optimization algorithm is presented.A multi-physics coupled acousticmechanics model based on the discrete phase model(DPM),large eddy simulation(LES)model,and Ffowcs Williams-Hawkings(FW-H)model is established,and the numerical simulation results of the multi-physics coupled acoustic-mechanics method are verified via experimental comparison.Based on the analysis of the multi-physics coupled acoustic-mechanics numerical simulation results,the effects of the water flow on the characteristics of the atomization flow distribution were obtained.A multi-physics coupled acoustic-mechanics numerical simulation result was employed to establish an orthogonal test database,and a multi-objective synergetic optimization algorithm was adopted to optimize the key parameters of the TFN.The optimal parameters are as follows:A gas flow of 0.94 m^(3)/h,water flow of 0.0237 m^(3)/h,orifice diameter of the self-excited vibrating cavity(SVC)of 1.19 mm,SVC orifice depth of 0.53 mm,distance between SVC and the outlet of nozzle of 5.11 mm,and a nozzle outlet diameter of 3.15 mm.The droplet particle size in the atomization flow field was significantly reduced,the spray distance improved by 71.56%,and the SPL data at each corresponding measurement point decreased by an average of 38.96%.The conclusions of this study offer a references for future TFN research.

CAS-ESM2.0 Successfully Reproduces Historical Atmospheric CO_(2) in a Coupled Carbon−Climate Simulation

The atmospheric carbon dioxide(CO_(2))concentration has been increasing rapidly since the Industrial Revolution,which has led to unequivocal global warming and crucial environmental change.It is extremely important to investigate the interactions among atmospheric CO_(2),the physical climate system,and the carbon cycle of the underlying surface for a better understanding of the Earth system.Earth system models are widely used to investigate these interactions via coupled carbon-climate simulations.The Chinese Academy of Sciences Earth System Model version 2(CAS-ESM2.0)has successfully fixed a two-way coupling of atmospheric CO_(2)with the climate and carbon cycle on land and in the ocean.Using CAS-ESM2.0,we conducted a coupled carbon-climate simulation by following the CMIP6 proposal of a historical emissions-driven experiment.This paper examines the modeled CO_(2)by comparison with observed CO_(2)at the sites of Mauna Loa and Barrow,and the Greenhouse Gases Observing Satellite(GOSAT)CO_(2)product.The results showed that CAS-ESM2.0 agrees very well with observations in reproducing the increasing trend of annual CO_(2)during the period 1850-2014,and in capturing the seasonal cycle of CO_(2)at the two baseline sites,as well as over northern high latitudes.These agreements illustrate a good ability of CAS-ESM2.0 in simulating carbon-climate interactions,even though uncertainties remain in the processes involved.This paper reports an important stage of the development of CAS-ESM with the coupling of carbon and climate,which will provide significant scientific support for climate research and China’s goal of carbon neutrality.

Stress-corrosion coupled damage localization induced by secondary phases in bio-degradable Mg alloys:phase-field modeling

In this study,a phase-field scheme that rigorously obeys conservation laws and irreversible thermodynamics is developed for modeling stress-corrosion coupled damage(SCCD).The coupling constitutive relationships of the deformation,phase-field damage,mass transfer,and electrostatic field are derived from the entropy inequality.The SCCD localization induced by secondary phases in Mg is numerically simulated using the implicit iterative algorithm of the self-defined finite elements.The quantitative evaluation of the SCCD of a C-ring is in good agreement with the experimental results.To capture the damage localization,a micro-galvanic corrosion domain is defined,and the buffering effect on charge migration is explored.Three cases are investigated to reveal the effect of localization on corrosion acceleration and provide guidance for the design for resistance to SCCD at the crystal scale.

Uncertainty quantification of mechanism motion based on coupled mechanism—motor dynamic model for ammunition delivery system

In this paper,a dynamic modeling method of motor driven electromechanical system is presented,and the uncertainty quantification of mechanism motion is investigated based on this method.The main contribution is to propose a novel mechanism-motor coupling dynamic modeling method,in which the relationship between mechanism motion and motor rotation is established according to the geometric coordination of the system.The advantages of this include establishing intuitive coupling between the mechanism and motor,facilitating the discussion for the influence of both mechanical and electrical parameters on the mechanism,and enabling dynamic simulation with controller to take the randomness of the electric load into account.Dynamic simulation considering feedback control of ammunition delivery system is carried out,and the feasibility of the model is verified experimentally.Based on probability density evolution theory,we comprehensively discuss the effects of system parameters on mechanism motion from the perspective of uncertainty quantization.Our work can not only provide guidance for engineering design of ammunition delivery mechanism,but also provide theoretical support for modeling and uncertainty quantification research of mechatronics system.