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.

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.

Coupled Numerical Simulation of Electromagnetic and Flow Fields in a Magnetohydrodynamic Induction Pump

Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.

Multi-scenario Simulation and Spatial-temporal Analysis of LUCC in China's Coastal Zone Based on Coupled SD-FLUS Model

Increased human activities in China's coastal zone have resulted in the depletion of ecological land resources.Thus,conducting current and future multi-scenario simulation research on land use and land cover change(LUCC)is crucial for guiding the healthy and sustainable development of coastal zones.System dynamic(SD)-future land use simulation(FLUS)model,a coupled simulation model,was developed to analyze land use dynamics in China's coastal zone.This model encompasses five scenarios,namely,SSP1-RCP2.6(A),SSP2-RCP4.5(B),SSP3-RCP4.5(C),SSP4-RCP4.5(D),and SSP5-RCP8.5(E).The SD model simulates land use demand on an annual basis up to the year 2100.Subsequently,the FLUS model determines the spatial distribution of land use for the near term(2035),medium term(2050),and long term(2100).Results reveal a slowing trend in land use changes in China's coastal zone from 2000–2020.Among these changes,the expansion rate of construction land was the highest and exhibited an annual decrease.By 2100,land use predictions exhibit high accuracy,and notable differences are observed in trends across scenarios.In summary,the expansion of production,living,and ecological spaces toward the sea remains prominent.Scenario A emphasizes reduced land resource dependence,benefiting ecological land protection.Scenario B witnesses an intensified expansion of artificial wetlands.Scenario C sees substantial land needs for living and production,while Scenario D shows coastal forest and grassland shrinkage.Lastly,in Scenario E,the conflict between humans and land intensifies.This study presents pertinent recommendations for the future development,utilization,and management of coastal areas in China.The research contributes valuable scientific support for informed,long-term strategic decision making within coastal regions.

Simulation of Underground Reservoir Stability of Pumped Storage Power Station Based on Fluid-Structure Coupling

Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.

Multi-scale simulation of diffusion behavior of deterrent in propellant

Concentration distribution of the deterrent in single-base propellant during the process of firing plays an important role in the ballistic properties of gun propellant in weapons. However, the diffusion coefficient calculated by molecular dynamics(MD) simulation is 6 orders of magnitude larger than the experimental values. Meanwhile, few simple and comprehensive theoretical models can explain the phenomenon and accurately predict the concentration distribution of the propellant. Herein, an onion model combining with MD simulation and finite element method of diffusion in propellants is introduced to bridge the gap between the experiments and simulations, and correctly predict the concentration distribution of deterrent. Furthermore, a new time scale is found to characterize the diffusion process. Finally, the time-and position-depended concentration distributions of dibutyl phthalate in nitrocellulose are measured by Raman spectroscopy to verify the correctness of the onion model. This work not only provides guidance for the design of the deterrent, but could be also extended to the diffusion of small molecules in polymer with different crystallinity.

Fluid simulation of the effect of a dielectric window with high temperature on plasma parameters in inductively coupled plasma

To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature.This high temperature can modulate the plasma characteristics to a large degree.We thus study the effect of dielectric window temperature on plasma parameters in two different ICP structures based on COMSOL software.The distributions of various plasma species are examined at different dielectric window temperatures.The concentration of neutral gas is found to be largely modulated at high dielectric window temperature,which further affects the electron collision probability with neutrals and the electron temperature.However,the electron density profiles are barely affected by the dielectric window temperature,which is mainly concentrated at the center of the reactor due to the fixed power input and pressure.

Experimental and simulation studies on similitude design method for shock responses of beam-plate coupled structure

The similitude theory helps to understand the physical behaviors of large structures through scaled models. Several papers have studied the similitude of shock issues. However, the dynamic similitude for shock responses of coupled structures is rarely incorporated in open studies. In this paper, scaling laws are derived for the shock responses and spectra of coupled structures. In the presented scaling laws, the geometric distortion and energy loss are considered. The ability of the proposed scaling laws is demonstrated in the simulation and experimental cases. In both cases, the similitude prediction for the prototype's time-domain waveform and spectrum is conducted with the scaled model and scaling laws. The simulation and experimental cases indicate that the predicted shock responses and spectra agree well with those of the prototype, which verifies the proposed scaling laws for predicting shock responses.

Modeling and Simulation of a Transmission Line Response to a 400 kV/400V Capacitor Coupled Substation

The access to electricity in rural areas is extremely limited, but it is crucial for all citizens. The population in rural areas of sub-Saharan African (SSA) countries is generally low, making it economically unfeasible to implement traditional rural electrification (CRE) projects due to the high cost of establishing the necessary distribution infrastructure. To address this cost issue, one alternative technology for rural electrification (URE) that can be explored is the Capacitor Coupled Substation (CCS) technology. CCS is a cost-effective solution for supplying electricity to rural areas. The research is necessitated by the need to offer a cost-effective technology for supplying electricity to sparsely populated communities. This paper examines the impact on the transmission network when a 400 kV/400V CCS is connected to it. The system response when a CCS is connected to the network was modeled using MATLAB/Si-mulink. The results, based on the fixed load of 80 kW, showed negligible interference on the transmission line voltage. However, there was minor impact on the parameters downstream of the tapping point. These findings were further supported by introducing a fault condition to the CCS, which showed that interferences with the CCS could affect the overall stability of the transmission network downstream of the tapping node, similar to the behavior of an unstable load.

Multi-scale coupling simulation in directional solidification of superalloy based on cellular automaton-finite difference method

Casting microstructure evolution is difficult to describe quantitatively by only a separate simulation of dendrite scale or grain scale, and the numerical simulation of these two scales is difficult to render compatible. A three-dimensional cellular automaton model couplling both dendritic scale and grain scale is developed to simulate the microstructure evolution of the nickel-based single crystal superalloy DD406. Besides, a macro–mesoscopic/microscopic coupling solution algorithm is proposed to improve computational efficiency. The simulation results of dendrite growth and grain growth of the alloy are obtained and compared with the results given in previous reports. The results show that the primary dendritic arm spacing and secondary dendritic arm spacing of the dendritic growth are consistent with the theoretical and experimental results. The mesoscopic grain simulation can be used to obtain results similar to those of microscopic dendrites simulation. It is indicated that the developed model is feasible and effective.

Phase field simulation for non-isothermal solidification of multicomponent alloys coupled with thermodynamics database

In order to quantitively model the real solidification process of industrial multicomponent alloys, a non-isothermal phase field model was studied for multicomponent alloy fully coupled with thermodynamic and diffusion mobility database, which can accurately predict the phase equilibrium, solute diffusion coefficients, specific heat capacity and latent heat release in the whole system. The results show that these parameters are not constants and their values depend on local concentration and temperature. Quantitative simulation of solidification in multicomponent alloys is almost impossible without such parameters available. In this model, the interfacial region is assumed to be a mixture of solid and liquid with the same chemical potentials, but with different composition. The anti-trapping current is also considered in the model. And this model was successfully applied to industrial A1-Cu-Mg alloy for the free equiaxed dendrite solidification process.

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.

Effect of coil and chamber structure on plasma radial uniformity in radio frequency inductively coupled plasma

Enhancing plasma uniformity can be achieved by modifying coil and chamber structures in radio frequency inductively coupled plasma(ICP)to meet the demand for large-area and uniformly distributed plasma in industrial manufacturing.This study utilized a two-dimensional self-consistent fluid model to investigate how different coil configurations and chamber aspect ratios affect the radial uniformity of plasma in radio frequency ICP.The findings indicate that optimizing the radial spacing of the coil enhances plasma uniformity but with a reduction in electron density.Furthermore,optimizing the coil within the ICP reactor,using the interior point method in the Interior Point Optimizer significantly enhances plasma uniformity,elevating it from 56%to 96%within the range of the model sizes.Additionally,when the chamber aspect ratio k changes from 2.8 to 4.7,the plasma distribution changes from a center-high to a saddleshaped distribution.Moreover,the plasma uniformity becomes worse.Finally,adjusting process parameters,such as increasing source power and gas pressure,can enhance plasma uniformity.These findings contribute to optimizing the etching process by improving plasma radial uniformity.

The Impact of Deformable Surrounding Rocks on Naturally-Fractured Reservoir Simulation

For the case of a fractured reservoir surrounded by deformable rocks, the appropriateness and applicability of the two common methods of coupling of flow and deformation, explicit (coupled) and implicit (uncoupled) methods are investigated. The explicit formulation is capable of modelling surrounding media;while the implicit coupling is unable to do so as deformation vector does not appear as a primary variable in the formulation. The governing differential equations and the finite element approximation of the governing equations for each of the methods are presented. Spatial discretization is achieved using the Galerkin method, and temporal discretisation using the finite difference technique. In the explicit model, coupling between flow and deformation is captured through volumetric strain compatibility amongst the phases within the system. In the implicit model, this is achieved by defining the pore space storativity as a function of the formation compressibility and the compressibility of the fluid phases within the pore space. The impact of rock deformability on early, intermediate and late time responses of fractured reservoir is investigated through several numerical examples. Salient features of each formulation are discussed and highlighted. It is shown that the implicit model is unable to capture the constraining effects of a non-yielding, surrounding rock, leading to incorrect projections of reservoir production irrespective of the history matching strategy adopted.

Multi-dimensional Simulation of Phase Change by a 0D-2D Model Coupling via Stefan Condition

Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.

COUPLED SIMULATION OF 3D ELECTRO-MAGNETO-FLOW FIELD IN HALL-HEROULT CELLS USING FINITE ELEMENT METHOD

Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.

Wave-tide-surge coupled model simulation for Typhoon Maemi

Reasonably accurate predictions of wave heights, current and elevations during storm events are vital information for marine operations and design of offshore and coastal structures in the surrounding seas of Korea Peninsula. Ocean circulation and wind-wave models have traditionally been run separately, but recent researches have identified potentially important interactions between current and wave motions. The coupled tide-surge and the WAM wave models at the atmospheric boundary layer and bottom boundary layer around the Korea Peninsula are applied for the Typhoon Maemi （0314） event. Communication between the models is aehievod using MPI. Results are compared with coastal tide gauges and moored wave buoys and comparisons are also made between wave computations from the coupled model and the independent third generation wave models. Results suggest that applying the fide-surge-coupled model can be an effective means of obtaining wave and storm surge predictions simultaneously.

A COUPLED MODEL FOR MUILTIPHASE FLUID FLOW AND SEDIMENTATION DEFORMATION IN OIL RESERVOIR AND ITS NUMERICAL SIMULATION

A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.

Numerical simulation of aluminum holding furnace with fluid-solid coupled heat transfer

To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.

A Hybrid Coupled Ocean–Atmosphere Model and Its Simulation of ENSO and Atmospheric Responses

A new hybrid coupled model(HCM) is presented in this study, which consists of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model. The ocean component is the intermediate ocean model(IOM)of the intermediate coupled model(ICM) used at the Institute of Oceanology, Chinese Academy of Sciences(IOCAS). The atmospheric component is ECHAM5, the fifth version of the Max Planck Institute for Meteorology atmospheric general circulation model. The HCM integrates its atmospheric and oceanic components by using an anomaly coupling strategy. A100-year simulation has been made with the HCM and its simulation skills are evaluated, including the interannual variability of SST over the tropical Pacific and the ENSO-related responses of the global atmosphere. The model shows irregular occurrence of ENSO events with a spectral range between two and five years. The amplitude and lifetime of ENSO events and the annual phase-locking of SST anomalies are also reproduced realistically. Despite the slightly stronger variance of SST anomalies over the central Pacific than observed in the HCM, the patterns of atmospheric anomalies related to ENSO,such as sea level pressure, temperature and precipitation, are in broad agreement with observations. Therefore, this model can not only simulate the ENSO variability, but also reproduce the global atmospheric variability associated with ENSO, thereby providing a useful modeling tool for ENSO studies. Further model applications of ENSO modulations by ocean–atmosphere processes, and of ENSO-related climate prediction, are also discussed.