Mechanical behavior and damage constitutive model of sandstone under hydro-mechanical (H-M) coupling

Underground engineering often passes through water-rich fractured rock masses, which are prone to fracture and instability under the long-term coupling of in-situ stress field and pore water(P-W) pressure, ultimately threatening the stability of underground structures. In order to explore the mechanical properties of rocks under H-M coupling, the corresponding damage constitutive(D-C) model has become the focus of attention. Considering the inadequacy of the current research on rock strength parameters,energy evolution characteristics and D-C model under H-M coupling, the mechanical properties of typical sandstone samples are discussed based on laboratory tests. The results show that the variation of characteristic stresses of sandstone under H-M coupling conforms to the normalized attenuation equation and Mohr-Coulomb(M-C) criterion. The P-W pressure mechanism of sandstone exhibits a dynamic change from softening effect to H-M fracturing effect. The closure stress is mainly provided by cohesive strength, while the initiation stress, damage stress, and peak stress are jointly dominated by cohesive strength and friction strength. In addition, residual stress is attributed to the friction strength formed by the bite of the fracture surface. Subsequently, the energy evolution characteristics of sandstone under H-M coupling were studied, and it was found that P-W pressure weakened the energy storage capacity and energy dissipation capacity of sandstone, and H-M fracturing was an important factor in reducing its energy storage efficiency. Finally, combined with energy dissipation theory and statistical damage theory, two types of D-C models considering P-W pressure are proposed accordingly, and the model parameters can be determined by four methods. The application results indicate that the proposed and modified D-C models have high reliability, and can characterize the mechanical behavior of sandstone under H-M coupling, overcome the inconvenience of existing D-C models due to excessive mechanical parameters,and can be applied to the full-range stress–strain process. The results are conducive to revealing the deformation and damage mechanisms of rocks under H-M coupling, and can provide theoretical guidance for related engineering problems.

Aging Characteristics of Lithium-Ion Battery Under Fast Charging Based on Electrochemical-thermalmechanical Coupling Model

The aging characteristics of lithium-ion battery(LIB)under fast charging is investigated based on an electrochemical-thermal-mechanical(ETM)coupling model.Firstly,the ETM coupling model is established by COMSOL Multiphysics.Subsequently,a long cycle test was conducted to explore the aging characteristics of LIB.Specifically,the effects of charging(C)rate and cycle number on battery aging are analyzed in terms of nonuniform distribution of solid electrolyte interface(SEI),SEI formation,thermal stability and stress characteristics.The results indicate that the increases in C rate and cycling led to an increase in the degree of nonuniform distribution of SEI,and thus a consequent increase in the capacity loss due to the SEI formation.Meanwhile,the increases in C rate and cycle number also led to an increase in the heat generation and a decrease in the heat dissipation rate of the battery,respectively,which result in a decrease in the thermal stability of the electrode materials.In addition,the von Mises stress of the positive electrode material is higher than that of the negative electrode material as the cycling proceeds,with the positive electrode material exhibiting tensile deformation and the negative electrode material exhibiting compressive deformation.The available lithium ion concentration of the positive electrode is lower than that of the negative electrode,proving that the tensile-type fracture occurring in the positive material under long cycling dominated the capacity loss process.The aforementioned studies are helpful for researchers to further explore the aging behavior of LIB under fast charging and take corresponding preventive measures.

Hardening mechanism and thermal-solid coupling model of laminar plasma surface hardening of 65 Mn steel

In the present work,the laminar plasma surface hardening method is employed to enhance the service life of metal components fabricated from 65 Mn steel.The mechanical and wear behaviors of the laminar plasma surface hardened 65 Mn steel were analyzed.The martensite transition transformation of the temperature of the laminar plasma-hardened 65 ferrite Mn steel was determined by a thermal-solid coupling model.Based on the orthogonal experimental results,the optimal hardening parameters were confirmed.The scanning velocity,quenching distance and arc current are 130 mm/min,50 mm and 120 A,respectively.The pearlites and ferrites are transformed into martensites in the hardened zone,while the ratio of martensite in the heataffected zone decreases with the increase in the hardening depth.Compared to the untreated 65Mn steel,the average hardness increases from 220 HV_(0.2)to 920 HV_(0.2)in the hardened zone and the corresponding absorbed power increases from 118.7 J to 175.5 J.At the same time,the average coefficient of friction(COF)decreases from 0.763 to 0.546,and the wear rate decreases from 5.39×10^(-6)mm^(3)/(N·m)to 2.95×10^(-6)mm^(3)/(N·m),indicating that the wear resistance of 65Mn steel could be significantly improved by using laminar surface hardening.With the same hardening parameters,the depth and width of the hardened zone predicted by the thermal-solid coupling model are 1.85 mm and 11.20 mm,respectively,which are in accordance with the experimental results;depth is 1.83 mm and width is 11.15 mm.In addition,the predicted hardness distributions of the simulation model are in accordance with the experimental results.These results indicate that the simulation model could effectively predict the microstructure characteristics of 65 Mn steel.

Degree of Freedom Analysis for Holographic MIMO Based on a Mutual-Coupling-Compliant Channel Model

Degree of freedom(DOF)is a key indicator for spatial multiplexing layers of a wireless channel.Traditionally,the channel of a multiple-input multiple-output(MIMO)half-wavelength dipole array has a DOF that equals the antenna number.However,recent studies suggest that the DOF could be less than the antenna number when strong mutual coupling is considered.We utilize a mutual-coupling-compliant channel model to investigate the DOF of the holographic MIMO(HMIMO)channel and give a upper bound of the DOF with strong mutual coupling.Our numerical simulations demonstrate that a dense array can support more DOF per unit aperture as compared with a half-wavelength MIMO system.

Electro-mechanical coupling properties of band gaps in an with periodically attached “spring-mass” resonators

The model of a locally resonant (LR) epoxy/PZT-4 phononic crystal (PC)nanobeam with “spring-mass” resonators periodically attached to epoxy is proposed. The corresponding band structures are calculated by coupling Euler beam theory, nonlocal piezoelectricity theory and plane wave expansion (PWE) method. Three complete band gaps with the widest total width less than 10GHz can be formed in the proposed nanobeam by comprehensively comparing the band structures of three kinds of LR PC nanobeams with resonators attached or not. Furthermore, influencing rules of the coupling fields between electricity and mechanics,“spring-mass” resonator, nonlocal effect and different geometric parameters on the first three band gaps are discussed and summarized. All the investigations are expected to be applied to realize the active control of vibration in the region of ultrahigh frequency.

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.

Optimization Control of Multi-Mode Coupling All-Wheel Drive System for Hybrid Vehicle

The all-wheel drive(AWD)hybrid system is a research focus on high-performance new energy vehicles that can meet the demands of dynamic performance and passing ability.Simultaneous optimization of the power and economy of hybrid vehicles becomes an issue.A unique multi-mode coupling(MMC)AWD hybrid system is presented to realize the distributed and centralized driving of the front and rear axles to achieve vectored distribution and full utilization of the system power between the axles of vehicles.Based on the parameters of the benchmarking model of a hybrid vehicle,the best model-predictive control-based energy management strategy is proposed.First,the drive system model was built after the analysis of the MMC-AWD’s drive modes.Next,three fundamental strategies were established to address power distribution adjustment and battery SOC maintenance when the SOC changed,which was followed by the design of a road driving force observer.Then,the energy consumption rate in the average time domain was processed before designing the minimum fuel consumption controller based on the equivalent fuel consumption coefficient.Finally,the advantage of the MMC-AWD was confirmed by comparison with the dynamic performance and economy of the BYD Song PLUS DMI-AWD.The findings indicate that,in comparison to the comparative hybrid system at road adhesion coefficients of 0.8 and 0.6,the MMC-AWD’s capacity to accelerate increases by 5.26%and 7.92%,respectively.When the road adhesion coefficient is 0.8,0.6,and 0.4,the maximum climbing ability increases by 14.22%,12.88%,and 4.55%,respectively.As a result,the dynamic performance is greatly enhanced,and the fuel savings rate per 100 km of mileage reaches 12.06%,which is also very economical.The proposed control strategies for the new hybrid AWD vehicle can optimize the power and economy simultaneously.

Time-dependent behavior and permeability evolution of limestone under hydro-mechanical coupling

Excessively high pore water pressure presents unpredictable risks to the safety of rock tunnels in mountainous regions that are predominantly composed of limestone. Investigating the creep characteristics and permeability evolution of limestone under varying hydrated conditions is crucial for a better understanding of the delayed deformation mechanisms of limestone rock tunnels. To this end, this paper initially conducts a series of multi-stage triaxial creep tests on limestone samples under varying pore water pressures. The experiment examines how pore water pressure affects limestone’s creep strain, strain rate, long-term strength, lifespan, and permeability, all within the context of hydraulicmechanical(HM) coupling. To better describe the creep behavior associated with pore water pressure, this paper proposes a new nonlinear fractional creep constitutive model. This constitutive model depicts the initial, steady-state, and accelerated phases of limestone’s creep behavior. Finally, the proposed model is applied to the numerical realization of deformation in limestone tunnel, validating the effectiveness of the proposed constitutive model in predicting tunnel’s creep deformation. This research enhances our understanding of limestone’s creep characteristics and permeability evolution under HM coupling, laying a foundation for assessing the longterm stability of mountain tunnels.

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.

An internal ballistic model of electromagnetic railgun based on PFN coupled with multi-physical field and experimental validation

To accelerate the practicality of electromagnetic railguns,it is necessary to use a combination of threedimensional numerical simulation and experiments to study the mechanism of bore damage.In this paper,a three-dimensional numerical model of the augmented railgun with four parallel unconventional rails is introduced to simulate the internal ballistic process and realize the multi-physics field coupling calculation of the rail gun,and a test experiment of a medium-caliber electromagnetic launcher powered by pulse formation network(PFN)is carried out.Various test methods such as spectrometer,fiber grating and high-speed camera are used to test several parameters such as muzzle initial velocity,transient magnetic field strength and stress-strain of rail.Combining the simulation results and experimental data,the damage condition of the contact surface is analyzed.

Different El Niño Flavors and Associated Atmospheric Teleconnections as Simulated in a Hybrid Coupled Model

A previously developed hybrid coupled model(HCM)is composed of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model(AGCM),denoted as HCMAGCM.In this study,different El Niño flavors,namely the Eastern-Pacific(EP)and Central-Pacific(CP)types,and the associated global atmospheric teleconnections are examined in a 1000-yr control simulation of the HCMAGCM.The HCMAGCM indicates profoundly different characteristics among EP and CP El Niño events in terms of related oceanic and atmospheric variables in the tropical Pacific,including the amplitude and spatial patterns of sea surface temperature(SST),zonal wind stress,and precipitation anomalies.An SST budget analysis indicates that the thermocline feedback and zonal advective feedback dominantly contribute to the growth of EP and CP El Niño events,respectively.Corresponding to the shifts in the tropical rainfall and deep convection during EP and CP El Niño events,the model also reproduces the differences in the extratropical atmospheric responses during the boreal winter.In particular,the EP El Niño tends to be dominant in exciting a poleward wave train pattern to the Northern Hemisphere,while the CP El Niño tends to preferably produce a wave train similar to the Pacific North American(PNA)pattern.As a result,different climatic impacts exist in North American regions,with a warm-north and cold-south pattern during an EP El Niño and a warm-northeast and cold-southwest pattern during a CP El Niño,respectively.This modeling result highlights the importance of internal natural processes within the tropical Pacific as they relate to the genesis of ENSO diversity because the active ocean–atmosphere coupling is allowed only in the tropical Pacific within the framework of the HCMAGCM.

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.

Theoretically predicted innovative palladium stripe dopingcobalt(111) surface with excellent catalytic performance for carbonmonoxide oxidative coupling to dimethyl oxalate

Pd-based catalysts are extensively employed to catalyze CO oxidative coupling to generate DMO,while the expensive price and high usage of Pd hinder its massive application in industrial production.Designing Pd-based catalysts with high efficiency and low Pd usage as well as expounding the catalytic mechanisms are significant for the reaction.In this study,we theoretically predict that Pd stripe doping Co(111)surface exhibits excellent performance than pure Pd(111),Pd monolayer supporting on Co(111)and Pd single atom doping Co(111)surface,and clearly expound the catalytic mechanisms through the density functional theory(DFT)calculation and micro-reaction kinetic model analysis.It is obtained that the favorable reaction pathway is COOCH_(3)-COOCH_(3)coupling pathway over these four catalysts,while the rate-controlling step is COOCH_(3)+CO+OCH_(3)→2COOCH_(3)on Pd stripe doping Co(111)surface,which is different from the case(2COOCH_(3)→DMO)on pure Pd(111),Pd monolayer supporting on Co(111)and Pd single atom doping Co(111)surface.This study can contribute a certain reference value for developing Pd-based catalysts with high efficiency and low Pd usage for CO oxidative coupling to DMO.

Simulation study of supercritical carbon dioxide jet fracturing for carbonate geothermal reservoir based on fluid-thermo-mechanical coupling model

Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon dioxide(SC-CO_(2))jet fracturing is expected to efficiently stimulate the carbonate geothermal reservoirs and achieve the storage of CO_(2) simultaneously.In this paper,we established a transient seepage and fluid-thermo-mechanical coupled model to analyze the impact performance of sc-CO_(2) jet fracturing.The mesh-based parallel code coupling interface was employed to couple the fluid and solid domains by exchanging the data through the mesh interface.The physical properties change of sC-CO_(2) with temperature were considered in the numerical model.Results showed that SC-CO_(2) jet frac-turing is superior to water-jet fracturing with respect to jetting velocity,particle trajectory and pene-trability.Besides,stress distribution on the carbonate rock showed that the tensile and shear failure would more easily occur by SC-CO_(2) jet than that by water jet.Moreover,pressure and temperature control the jet field and seepage field of sC-CO_(2) simultaneously.Increasing the jet temperature can effectively enhance the impingement effect and seepage process by decreasing the viscosity and density of SC-CO_(2).The key findings are expected to provide a theoretical basis and design reference for applying SC-CO_(2) jet fracturing in carbonate geothermal reservoirs.

Hydrologic Response to Future Climate Change in the Dulong-Irra-waddy River Basin Based on Coupled Model Intercomparison Project 6

Within the context of the Belt and Road Initiative(BRI)and the China-Myanmar Economic Corridor(CMEC),the Dulong-Ir-rawaddy(Ayeyarwady)River,an international river among China,India and Myanmar,plays a significant role as both a valuable hydro-power resource and an essential ecological passageway.However,the water resources and security exhibit a high degree of vulnerabil-ity to climate change impacts.This research evaluates climate impacts on the hydrology of the Dulong-Irrawaddy River Basin(DIRB)by using a physical-based hydrologic model.We crafted future climate scenarios using the three latest global climate models(GCMs)from Coupled Model Intercomparison Project 6(CMIP6)under two shared socioeconomic pathways(SSP2-4.5 and SSP5-8.5)for the near(2025-2049),mid(2050-2074),and far future(2075-2099).The regional model using MIKE SHE based on historical hydrologic processes was developed to further project future streamflow,demonstrating reliable performance in streamflow simulations with a val-idation Nash-Sutcliffe Efficiency(NSE)of 0.72.Results showed that climate change projections showed increases in the annual precip-itation and potential evapotranspiration(PET),with precipitation increasing by 11.3%and 26.1%,and PET increasing by 3.2%and 4.9%,respectively,by the end of the century under SSP2-4.5 and SSP5-8.5.These changes are projected to result in increased annual streamflow at all stations,notably at the basin’s outlet(Pyay station)compared to the baseline period(with an increase of 16.1%and 37.0%at the end of the 21st century under SSP2-4.5 and SSP5-8.5,respectively).Seasonal analysis for Pyay station forecasts an in-crease in dry-season streamflow by 31.3%-48.9%and 22.5%-76.3%under SSP2-4.5 and SSP5-8.5,respectively,and an increase in wet-season streamflow by 5.8%-12.6%and 2.8%-33.3%,respectively.Moreover,the magnitude and frequency of flood events are pre-dicted to escalate,potentially impacting hydropower production and food security significantly.This research outlines the hydrological response to future climate change during the 21st century and offers a scientific basis for the water resource management strategies by decision-makers.

Numerical Modelling of Coupled Heat and Mass Transfer in Porous Materials: Application to Cinder Block Bricks

In this work, we present numerical modelling of coupled heat and mass transfer within porous materials. Our study focuses on cinder block bricks generally used in building construction. The material is assumed to be placed in air. Moisture content and temperature have been chosen as the main transfer drivers and the equations governing these transfer drivers are based on the Luikov model. These equations are solved by an implicit finite difference scheme. A Fortran code associated with the Thomas algorithm was used to solve the equations. The results show that heat and mass transfer depend on the temperature of the air in contact with the material. As this air temperature rises, the temperature within the material increases, and more rapidly at the material surface. Also, thermal conductivity plays a very important role in the thermal conduction of building materials and influences heat and mass transfer in these materials. Materials with higher thermal conductivity diffuse more heat.

Mathematical Modeling of Multiple Capacitor Coupled Substations (CCS) Impact on Transmission Lines and Approaches for Ferroresonance Suppression

Rural electrification remains a critical challenge in achieving equitable access to electricity, a cornerstone for poverty alleviation, economic growth, and improved living standards. Capacitor Coupled Substations (CCS) offer a promising solution for delivering cost-effective electricity to these underserved areas. However, the integration of multiple CCS units along a transmission network introduces complex interactions that can significantly impact voltage, current, and power flow. This study presents a detailed mathematical model to analyze the effects of varying distances and configurations of multiple CCS units on a transmission network, with a focus on voltage stability, power quality, and reactive power fluctuations. Furthermore, the research addresses the phenomenon of ferroresonance, a critical issue in networks with multiple CCS units, by developing and validating suppression strategies to ensure stable operation. Through simulation and practical testing, the study provides insights into optimizing CCS deployment, ultimately contributing to more reliable and efficient rural electrification solutions.

Prediction of Sedimentary Microfacies Distribution by Coupling Stochastic Modeling Method in Oil and Gas Energy Resource Exploitation

In view of the problem that a single modeling method cannot predict the distribution of microfacies, a new idea of coupling modeling method to comprehensively predict the distribution of sedimentary microfacies was proposed, breaking the tradition that different sedimentary microfacies used the same modeling method in the past. Because different sedimentary microfacies have different distribution characteristics and geometric shapes, it is more accurate to select different simulation methods for prediction. In this paper, the coupling modeling method was to establish the distribution of sedimentary microfacies with simple geometry through the point indicating process simulation, and then predict the microfacies with complex spatial distribution through the sequential indicator simulation method. Taking the DC block of Bohai basin as an example, a high-precision reservoir sedimentary microfacies model was established by the above coupling modeling method, and the model verification results showed that the sedimentary microfacies model had a high consistency with the underground. The coupling microfacies modeling method had higher accuracy and reliability than the traditional modeling method, which provided a new idea for the prediction of sedimentary microfacies.

Study on damage-stress loss coupling model of rock and prestressed anchor cable in dry-wet environment

The loss of anchoring force is one of the problems to be solved urgently.The anchorage loss is a key factor causing the failure of anchoring engineering,so it is crucial to study the time-dependent variation of anchoring force.Alternating dry and wet(D-W)conditions have a significant effect on deformation of rock.The anchoring system is composed of anchoring components and rock mass,and thus rock deformation has a significant impact on the loss of anchoring force.Quantifying rock deformation under the effects of D-W cycles is a prerequisite to understanding the factors that influence loss of anchoring force in anchor bolts.In this study,we designed an anchoring device that enabled real-time monitoring of the variation in strain during D-W periods and rock testing.Nuclear magnetic resonance(NMR)measurements showed that under D-W conditions,the increment in porosity was smaller for prestressed rock than unstressed rock.The trends of prestress loss and strain variation are consistent,which can be divided into three characteristic intervals:rapid attenuation stage,slow attenuation stage and relatively stable stage.At the same stress level,the rate of stress loss and strain for the soaking specimen was the highest,while that of the dried specimen was the lowest.In the same D-W cycling conditions,the greater the prestress,the smaller the strain loss rate of the rock,especially under soaking conditions.The characteristics of pore structure and physical mechanical parameters indicated that prestress could effectively suppress damage caused by erosion related to D-W cycles.The study reveals the fluctuation behavior of rock strain and prestress loss under D-W conditions,providing a reference for effectively controlling anchoring loss and ideas for inventing new anchoring components.

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.