Numerical modeling for the coupled thermo-mechanical processes and spalling phenomena in sp Pillar Stability Experiment (APSE)

In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both continuum and discontinuum based numerical methods. Two-dimensional （2D） and three- dimensional （3D） finite element method （FEM） and 2D distinct element method （DEM） with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone （EDZ） during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.

Design and Coupled Thermo-Mechanical Analysis of Silicon Carbide Primary Mirror Assembly

Based on the principle that the thermal expansion coefficient of the support structure should match that of the mirror, three schemes of primary mirror assembly were designed. Of them, the first is fused silica mirror plus 4J32 flexible support plus ZTC4 support back plate, the second K9 mirror plus 4J45 flexible support plus ZTC4 support back plate, and the third SiC mirror plus SiC rigid support back plate. A coupled thermo-mechanical analysis of the three primary mirror assemblies was made with finite element method. The results show that the SiC assembly is the best of all schemes in terms of their combination properties due to its elimination of the thermal expansion mismatch between the materials. The analytical results on the cryogenic property of the SiC primary mirror assembly show a higher surface finish of the SiC mirror even under the cryogenic condition.

Numerical modeling of coupled thermo-mechanical response of a rock pillar

Understanding the rock mass response to excavation and thermal loading and improving the capability of the numerical models for simulating the progressive failure process of brittle rocks are important for safety assessment and optimization design of nuclear waste repositories.The international cooperative DECOVALEX-2011 project provides a platform for development,validation and comparison of numerical models,in which the sp pillar stability experiment（APSE） was selected as the modeling target for Task B.This paper presents the modeling results of Wuhan University（WHU） team for stages 1 and 2 of Task B by using a coupled thermo-mechanical model within the framework of continuum mechanics.The rock mass response to excavation is modeled with linear elastic,elastoplastic and brittle-plastic models,while the response to heating is modeled with a coupled thermo-elastic model.The capabilities and limitations of the model for representation of the thermo-mechanical responses of the rock pillar are discussed by comparing the modeling results with experimental observations.The results may provide a helpful reference for the stability and safety assessment of the hard granite host rock in China＇s Beishan preselected area for high-level radioactive waste disposal.

Numerical study on coupled thermo-mechanical processes in sp Pillar Stability Experiment

This paper presents a study of the full three-dimensional thermo-mechanical （TM） behavior of rock pillar in,Aspo Pillar Stability Experiment （APSE） using a self-developed numerical code TM-EPCA3D. The transient thermal conduction function was descritized on space and time scales, and was solved by using cellular automaton （CA） method on space scale and finite difference method on time scale, respectively. The advantage of this approach is that no global, but local matrix is used so that it avoids the need to develop and solve large-scale linear equations and the complexity therein. A thermal conductivity versus stress function was proposed to reflect the effect of stress on thermal field. The temperature evolution and induced thermal stress in the pillar part during the heating and cooling processes were well simulated by the developed code. The factors that affect the modeling results were discussed. It is concluded that, the complex TM behavior of Aspo rock pillar is significantly influenced by the complex boundary and initial conditions.

Coupled thermo-mechanical constitutive damage model for sandstone

Underground rock dynamic disasters are becoming more severe due to the increasing depth of human operations underground.Underground temperature and pressure conditions contribute significantly to these disasters.Therefore,it is important to understand the coupled thermo-mechanical(TM)behaviour of rocks for the long-term safety and maintenance of underground tunnelling and mining.Moreover,investigation of the damage,strength and failure characteristics of rocks under triaxial stress conditions is important to avoid underground rock disasters.In this study,based on Weibull distribution and Lemaitre's strain equivalent principle,a statistical coupled TM constitutive model for sandstone was established under high temperature and pressure conditions.The triaxial test results of sandstone under different temperature and pressure conditions were used to validate the model.The proposed model was in good agreement with the experimental results up to 600℃.The total TM damage was decreased with increasing temperature,while it was increased with increasing confining pressure.The model's parameters can be calculated using conventional laboratory test results.

Numerical Simulation of Coupled Thermo-mechanical Behavior of a Cylinder Billet during Hot-forging Process

A finite volume method is applied to simulate a closed die hot forging process of a cylinder billet. Since variation and distribution of temperature play very important role in hot forging, the code involves a methodology of a coupled system of mechanical and thermal equations. The simulated results are compared with the experimental ones. The distribution of temperature in the billet obtained from the simulation is also discussed.

Coupled Thermo-Mechanical Modeling of Long Cone-Shaped Forging Based on Rigid Visco-Plastic FEM

Based on synthetically considering the coupled thermo mechanical relations between temperature and deforming, a numerical simulation of the forging process for the special long cone shaped workpiece of Al 5.44Mg 2.15Li 0.12Zr alloy at high temperature was conducted by using the rigid visco plastic finite element method. The relations between the total load and the displacement during the forging, and the distributions of stress, strain, temperature and strain rate, which can provide useful information for the process design, are obtained.

Numerical analyses of coupled thermo-mechanical behaviors of cement sheath in geothermal wells

The stability of cement sheath under high temperature and high pressure is one of the most critical issues for the durability of geother-mal well systems.In this study,a two-dimensional plane-strain finite element code was developed to investigate the coupled thermo-mechanical behaviors of the casing-cement-formation system.Different from previous linear elastic analyses,a thermoelasto-plastic con-stitutive model based on the thermodynamic theory was adopted for the cement sheath.It is shown that the finite element simulations using the proposed model provide a more accurate and realistic prediction of stress–strain responses of the cement sheath under high temperature.The results demonstrate that the radial stress concentration and the tensile strain concentration occur at both the cement–casing interface and the cement–formation interface,where the cement sheath is most likely to fail.High strength and low stiff-ness in the cement sheath and the formation are preferred for the integrity of the system.Both large thermal cycles and large differences between the internal fluid pressure and the external pressure should be avoided during operation.The new code is an alternative tool for guiding the geothermal well design.The finite element framework described herein is universal for other thermo-mechanical applications,such as energy foundations and energy tunnels.

Shear behavior of intact granite under thermo-mechanical coupling and three-dimensional morphology of shear-formed fractures

The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.

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.

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.

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.

Dynamic thermo-mechanical responses of road-soft ground system under vehicle load and daily temperature variation

A complete road-soft ground model is established in this paper to study the dynamic responses caused by vehicle loads and/or daily temperature variation.A dynamic thermo-elastic model is applied to capturing the behavior of the rigid pavement,the base course,and the subgrade,while the soft ground is characterized using a dynamic thermo-poroelastic model.Solutions to the road-soft ground system are derived in the Laplace-Hankel transform domain.The time domain solutions are obtained using an integration approach.The temperature,thermal stress,pore water pressure,and displacement responses caused by the vehicle load and the daily temperature variation are presented.Results show that obvious temperature change mainly exists within 0.3 m of the road when subjected to the daily temperature variation,whereas the stress responses can still be found in deeper places because of the thermal swelling/shrinkage deformation within the upper road structures.Moreover,it is important to consider the coupling effects of the vehicle load and the daily temperature variation when calculating the dynamic responses inside the road-soft ground system.

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.

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.

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.

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.

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.