Effect of dry-wet cycles on dynamic properties and microstructures of sandstone:Experiments and modelling

Underground pumped storage power plant(UPSP)is an innovative concept for space recycling of abandoned mines.Its realization requires better understanding of the dynamic performance and durability of reservoir rock.This paper conducted ultrasonic detection,split Hopkinson pressure bar(SHPB)impact,mercury intrusion porosimetry(MIP),and backscatter electron observation(BSE)tests to investigate the dynamical behaviour and microstructure of sandstone with cyclical dry-wet damage.A coupling FEM-DEM model was constructed for reappearing mesoscopic structure damage.The results show that dry-wet cycles decrease the dynamic compressive strength(DCS)with a maximum reduction of 39.40%,the elastic limit strength is reduced from 41.75 to 25.62 MPa.The sieved fragments obtain the highest crack growth rate during the 23rd dry-wet cycle with a predictable life of 25 cycles for each rock particle.The pore fractal features of the macropores and micro-meso pores show great differences between the early and late cycles,which verifies the computational statistics analysis of particle deterioration.The numerical results show that the failure patterns are governed by the strain in pre-peak stage and the shear cracks are dominant.The dry-wet cycles reduce the energy transfer efficiency and lead to the discretization of force chain and crack fields.

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.

Determination of energy coupling to material in laser welding by a novel “sandwich” method

A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the keyhole and the inverse Bremsstrahlung absorption coefficient in the keyhole plasma are obtained from the experiments;2） instead of using a parallel incident beam,a focused laser beam with real Gaussian intensity distribution is implemented;3） both Fresnel absorption and inverse Bremsstrahlung absorption during multiple reflections are considered.The calculation results show that the distribution of absorbed laser intensity by the keyhole wall is not uniform.The maximum laser energy is absorbed by the bottom of the keyhole,although no rays irradiate directly onto the bottom.According to analysis of beam focusing characteristics,the location of the focal plane plays a more important role in the laser energy absorption by the front wall than by the rear wall.

COUPLING MODEL OF EXTENDED MANUFACTURING ORGANIZATION AND ITS APPLICATION

For the feature of complex weapon manufacturing on internet,a coupling model is proposed.By using the model,the correlation between manufacturing cells in an extended manufacturing organization can be evaluated quantitatively,so an appropriate control plan is determined.A strategy to improve and reduce the coupling relationship of the organization is studied.A correlation matrix of extended tasks is built to analyze the relationship between sub-tasks and manufacturing resources.An optimization method for manufacturing resource configuration is presented based on the coupling model.Finally,a software system for analyzing coupling model about manufacturing organization on internet is developed,and the result shows that the coupling model is effective.

Application of 1D/3D finite elements coupling for structural nonlinear analysis

An element coupling model (ECM) method was proposed to simulate the global behavior and local damage of a structure.In order to reflect the local damage and improve the computational efficiency,three-dimensional (3D) solid elements and one-dimensional (1D) beam element were coupled by the multi-point constraint equations.A reduced scale 1?8 model test was simulated by the ECM and a full three dimensional model (3DM) contrastively.The results show that the global behavior and local damages of ECM agree well with the test and 3DM.It is indicated that the proposed method can be used in the structural nonlinear analysis accurately and efficiently.

The dynamic coupling model and its application of urbanization and eco-environment in Hexi Corridor

This paper, taking Hexi Corridor as an example, analyzes the altemating intimidation and the dynamic evolving relation between urbanization and eco-environment in arid area of West China. We argue that the harmonious development system of the urbanization and eco-environment would go through four phases： rudimentary symbiotic phase, harmonious developmental phase, utmost increasing phase and spiral type rising phase. Throughout the four phases, the elements of the system would influence each other, coerce each other, and complete the spiral type rising process from low-grade symbiosis to high-grade harmony together. The study on Hexi Corridor shows that the urbanization level in Hexi Corridor has increased gradually from 1985 to 2003 accompanied with the fluctuations of eco-environment state. The response of eco-environment to urbanization has been evident, but lagged behind the urbanization course. At present, the harmonious development system in Hexi Corridor was in its harmonious developmental phase. However, the coupling degree has increased quickly and approached 90 yet, which is signaling that the system is about to enter the utmost increasing phase, and the ecological crisis will enter the latent period. We have found that the coupling degree can well reflect the interactive coercing and dynamic evolving situation between urbanization and eco-environment in Hexi Corridor. From the temporal change of the coupling degree, it can be concluded that urbanization sometimes needs to pay a certain cost for the damage of the eco-environment in its initial stages, but as the urbanization continues, the state of the eco-environment would be meliorated.

COUPLED PHYSICAL-ECOLOGICAL MODELLING IN THE CENTRAL PART OF JIAOZHOU BAY Ⅱ. COUPLED WITH AN ECOLOGICAL MODEL

Sharples’ 1 D physical model employing tide wind driven turbulence closure and surface heating cooling physics, was coupled with an ecological model with 9 biochemical components: phytoplankton, zooplankton, shellfish, autotrophic and heterotrophic bacterioplankton, dissolved organic carbon (DOC), suspended detritus and sinking particles to simulate the annual evolution of ecosystem in the central part of Jiaozhou Bay. The coupled modeling results showed that the phytoplankton shading effect could reduce seawater temperature by 2℃, so that photosynthesis efficiency should be less than 8%; that the loss of phytoplankton by zooplankton grazing in winter tended to be compensated by phytoplankton advection and diffusion from the outside of the Bay; that the incident irradiance intensity could be the most important factor for phytoplankton growth rate; and that it was the bacterial secondary production that maintained the maximum zooplankton biomass in winter usually observed in the 1990s, indicating that the microbial food loop was extremely important for ecosystem study of Jiaozhou Bay.

Study on inhomogeneous cooling behavior of extruded profile with unequal and large thicknesses during quenching using thermo-mechanical coupling model

The interfacial heat transfer coefficient between hot profile surface and cooling water was determined by using inverse heat conduction model combined with end quenching experiment. Then, a Deform-3 D thermo-mechanical coupling model for simulating the on-line water quenching of extruded profile with unequal and large thicknesses was developed. The temperature field, residual stress field and distortion of profile during quenching were investigated systematically. The results show that heat transfer coefficient increases as water flow rate increases. The peak heat transfer coefficient with higher water flow rates appears at lower interface temperatures. The temperature distribution across the cross-section of profile during quenching is severe nonuniform and the maximum temperature difference is 300 ℃ at quenching time of 3.49 s. The temperature difference through the thickness of different parts of profile first increases sharply to a maximum value, and then gradually decreases. The temperature gradient increases obviously with the increase of thickness of parts. After quenching, there exist large residual stresses on the inner side of joints of profile and the two ends of part with thickness of 10 mm. The profile presents a twisting-type distortion across the cross-section under non-uniform cooling and the maximum twisting angle during quenching is 2.78°.

A Numerical Study on Effects of Land-Surface Heterogeneity from' Combined Approach' on Atmospheric ProcessPart II: Coupling-Model Simulations

Two land surface schemes, one the standard Biosphere / Atmosphere Transfer Scheme Version ie (BOZ) and the other B1Z based on B0Z and heterogeneously-treated by' combined approach' , were co 'pled to the meso-scale model MM4, respectively. Through the calculations of equations from the companion paper, parameters representing land surface heterogeneity and suitable for the coupling models were found out. Three cases were simulated for heavy rainfalls during 36 hours, and the sensitivity of short-term weather modeling to the land surface heterogeneity was tested. Through the analysis of the simulations of the three heavy rainfalls, it was demonstrated that BIZ, compared with BOZ, could more realistically reflect the features of the land surface heterogeneity, therefore could more realistically reproduce the circulation and precipitation amount in the heavy rainfall processes of the three cases. This shows that even short-term weather is sensitive to the land surface heterogeneity, which is more obvious with time passing, and whose influence is more pronounced in the lower layer and gradually extends to the middle and upper layer. Through the analysis of these simulations with BlZ, it is suggested that the bulk effect of smaller-scale fluxes (i.e., the momentum, water vapor and sensible heat fluxes) near the s ig nificantly-heterogeneous land surface is to change the larger-scale (i.e., meso-scale) circulation, and then to influence the development of the low-level jets and precipitation. And also, the complexity of the land-atmosphere interaction was shown in these simulations.

Coupling effects of influence factors on probability of corrosion initiation time of reinforced concrete

The coupling model of major influence factors such state affecting the chloride diffusion process in concrete is as environmental relative humidity, load-induced crack and stress discussed. The probability distributions of the critical chloride concentration Cc, the chloride diffusion coefficient D, and the surface chloride concentration Cs were determined based on the collected natural exposure data. And the estimation of probability of corrosion initiation considering the coupling effects of influence factors is presented. It is found that the relative humidity and curing time are the most effective factors affecting the probability of corrosion initiation before and after 10 years of exposure time. At the same exposure time, the influence of load-induced crack and stress state on the probability of corrosion initiation is obvious, in which the effect of crack is the most one

A model for coupling reservoir inflow and wellbore flow in fishbone wells

A coupling model is proposed in this paper by using the Green Function and Newman＇s product principle, and the solution method is provided here as well. This model can be used to describe the reservoir inflow and wellbore flow for fishbone wells in an unsteady flow or pseudo-steady flow state. A case study indicates that the bottom hole pressure declines quickly in the unsteady flow period which is very short. The pressure drop per unit time remains unchanged under the pseudo-steady flow conditions. The distribution of flow rate along the main wellbore shows a wave shape under the unsteady flow condition, and the flow rate distribution in each branch is similar. The flow rate distribution along the main wellbore is irregular ＂U＂ shaped under the pseudo-steady flow condition, and the space-symmetrical branches have the same flow distribution pattern. In the initial production period, the flow rate increases significantly as the length of branches and the angle between branches and the main wellbore increase. As the production continues, the length and angle of branches have only a slight effect on the flow in fishbone wells.

Roll System and Stock's Multi-parameter Coupling Dynamic Modeling Based on the Shape Control of Steel Strip

The existence of rolling deformation area in the rolling mill system is the main characteristic which dis- tinguishes the other machinery. In order to analyze the dynamic property of roll system＇s flexural deformation, it is necessary to consider the transverse periodic movement of stock in the rolling deformation area which is caused by the flexural deformation movement of roll system simul- taneously. Therefore, the displacement field of roll system and flow of metal in the deformation area is described by kinematic analysis in the dynamic system. Through intro- ducing the lateral displacement function of metal in the deformation area, the dynamic variation of per unit width rolling force can be determined at the same time. Then the coupling law caused by the co-effect of rigid movement and flexural deformation of the system structural elements is determined. Furthermore, a multi-parameter coupling dynamic model of the roll system and stock is established by the principle of virtual work. More explicitly, the cou- pled motion modal analysis was made for the roll system. Meanwhile, the analytical solutions for the flexural defor- mation movement＇s mode shape functions of rolls are discussed. In addition, the dynamic characteristic of the lateral flow of metal in the rolling deformation area has been analyzed at the same time. The establishment ofdynamic lateral displacement function of metal in the deformation area makes the foundation for analyzing the coupling law between roll system and rolling deformation area, and provides a theoretical basis for the realization of the dynamic shape control of steel strip.

Coupling non-isothermal trickle-bed reactor with catalyst pellet models to understand the reaction and diffusion in gas oil hydrodesulfurization

In this work,a trickle-bed reactor coupled with catalyst pellet model is employed to understand the effects of the temperature and catalyst pellet structures on the reaction-diffusion behaviors in gas oil hydrodesulfurization(HDS).The non-isothermal reactor model is determined to be reasonable due to non-negligible temperature variation caused by the reaction heat.The reaction rate along the reactor is mainly dominated by the temperature,and the sulfur concentration gradient in the catalyst pellet decreases gradually along the reactor,leading to the increased internal effectiveness factor.For the fixed catalyst bed volume,there exists a compromise between the catalyst reaction rate and effectiveness factor.Under commonly studied catalyst pellet size of 0.8-3 mm and porosity of 0.4-0.8,an optimization of the temperature and catalyst pellet structures is carried out,and the optimized outlet sulfur content decreases to 7.6 wppm better than the commercial level at 0.96 mm of the catalyst pellet size and 0.40 of the catalyst porosity.

An integrated method to calculate the spatial distribution of overburden strata failure in longwall mines by coupling GIS and FLAC^(3D)

The spatial distribution of overburden strata failure is of significant importance to affect the safety of underground mining. Because the traditional methods cannot be applied in all coal mines due to geological conditions or mining structures, a method of coupling FLAC3D with GIS was presented to calculate the spatial distribution of overburden strata failure in longwall coal mines. After building the spatio-temporal database from the calculation results of FLAC3D, the height of the mining-induced fractured zone in the overburden strata can be calculated by using the given height function. The results of case study show that the height of the fractured zone reached the maximum value at the face advance equal to about the panel width. The outcome of the work presented will be helpful in practice to improve safety in the production.

Characteristic analysis of mechanical thermal coupling model for bearing rotor system of high-speed train

Based on Newton’s second law and the thermal network method,a mechanical thermal coupling model of the bearing rotor system of high-speed trains is established to study the interaction between the bearing vibration and temperature.The influence of lubrication on the vibration and temperature characteristics of the system is considered in the model,and the real-time relationship between them is built up by using the transient temperature field model.After considering the lubrication,the bearing outer ring vibration acceleration and node temperature considering grease are lower,which shows the necessity of adding the lubrication model.The corresponding experiments for characteristics of vibration and temperature of the model are respectively conducted.In the envelope spectrum obtained from the simulation signal and the experimental signal,the frequency values corresponding to the peaks are close to the theoretical calculation results,and the error is very small.In the three stages of the temperature characteristic experiment,the node temperature change of the simulation model is consistent with the experiment.The good agreement between simulation and experiments proves the effectiveness of the model.By studying the influence of the bearing angular and fault size on the system node temperature,as well as the change law of bearing lubrication characteristics and temperature,it is found that the worse the working condition is,the higher the temperature is.When the ambient temperature is low,the viscosity of grease increases,and the oil film becomes thicker,which increases the sliding probability of the rolling element,thus affecting the normal operation of the bearing,which explains the phenomenon of frequent bearing faults of high-speed trains in the low-temperature area of Northeast China.Further analysis shows that faults often occur in the early stage of train operation in the low-temperature environment.

Analytical solution of vacuum preloading technology combined with electroosmosis coupling considering impacts of distribution of soil’s electrical potential

Combining vacuum preloading technology and electroosmosis can improve the treatment effect of soft soil foundation by utilizing the advantages of both methods.Many studies indicate that the soil electrical potential is non-linearly distributed in the treatment process by the combined method.However,in the previous theoretical study,the non-linear-distribution impacts of soil’s electrical potential on soft soil foundation treatment have not been considered.It is always assumed to be linear distribution,which is different from the experimental results.In this paper,the coupling consolidation model of this technology under the two-dimensional plane strain condition is initially established;and the well resistance effect,the vacuum load decreasing along the soil depth and the non-linear variation of electrical potential in the soil are considered.Then,the analytical solutions of the average excess pore water pressure and soil’s consolidation degree in the anode affected area are acquired based on the soil’s electrical potential distribution.Finally,the rationality of the analytical solution is testified by conducting an experimental model test,which proves the scientificity of the analytical solution.The analytical solution is adopted to better predict the dissipation of excess pore water pressure and soil consolidation degree when using the combined technology.This study can provide a reference with more accuracy for the engineering practices of this combined technology in the future.

Coupling numerical simulation with remotely sensed information for the study of frozen soil dynamics

The acquisition of spatial-temporal information of frozen soil is fundamental for the study of frozen soil dynamics and its feedback to climate change in cold regions.With advancement of remote sensing and better understanding of frozen soil dynamics,discrimination of freeze and thaw status of surface soil based on passive microwave remote sensing and numerical simulation of frozen soil processes under water and heat transfer principles provides valuable means for regional and global frozen soil dynamic monitoring and systematic spatial-temporal responses to global change.However,as an important data source of frozen soil processes,remotely sensed information has not yet been fully utilized in the numerical simulation of frozen soil processes.Although great progress has been made in remote sensing and frozen soil physics,yet few frozen soil research has been done on the application of remotely sensed information in association with the numerical model for frozen soil process studies.In the present study,a distributed numerical model for frozen soil dynamic studies based on coupled water-heat transferring theory in association with remotely sensed frozen soil datasets was developed.In order to reduce the uncertainty of the simulation,the remotely sensed frozen soil information was used to monitor and modify relevant parameters in the process of model simulation.The remotely sensed information and numerically simulated spatial-temporal frozen soil processes were validated by in-situ field observations in cold regions near the town of Naqu on the East-Central Tibetan Plateau.The results suggest that the overall accuracy of the algorithm for discriminating freeze and thaw status of surface soil based on passive microwave remote sensing was more than 95%.These results provided an accurate initial freeze and thaw status of surface soil for coupling and calibrating the numerical model of this study.The numerically simulated frozen soil processes demonstrated good performance of the distributed numerical model based on the coupled water-heat transferring theory.The relatively larger uncertainties of the numerical model were found in alternating periods between freezing and thawing of surface soil.The average accuracy increased by about 5%after integrating remotely sensed information on the surface soil.The simulation accuracy was significantly improved,especially in transition periods between freezing and thawing of the surface soil.

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.

Simplification and Experimental Verification for Temperature and Humidity Field Coupling Model of Conservatory Soil

Based on the surface energy balance model which is widely used abroad, a temperature and humidity field coupling model of conservatory soil without crop vegetation in full illumination was established. Considering the relatively closed environment in conservatory, weak solar radiation and little surface evaporation of soil, the daily variation of water content in different soil layers may be neglected, then the temperature and humidity field coupling model was simplified to a one-dimensional thermal diffusion model. The simplified model and the temperature and humidity field coupling model adopted the same computational method of soil physical parameters and discrete format of heat diffusion differential equations, and were applied to the continuous simulation of temperature field in conservatory soil without crop vegetation in full illumination. Through the comparison between simulation results and experimental data, the precision of the simplified model was verified. The typical rule of soil heat flux variation in a 24 h cycle was also obtained.

Modelling the response of larch growth to age,density,and elevation and the implications for multifunctional management in northwest China

Plantations of Rupprecht's larch(Larix principis-rupprechtii)have been widely established in the drylands of northwest and north China under traditional fastgrowing plantation management strategies.These strategies and the long-term logging ban have led to over-populated stands with lower structural and functional stability,less economic benefit and higher water consumption.To guide the sustainable management of larch plantations,field surveys and historical data compilation were undertaken in the Liupan Mountains of northwest China.The main influencing factors(stand structure and site condition)and their effects on mean tree height,mean DBH and timber volumes were determined based on up-boundary line analysis.Tree growth models coupling the effects of tree age,stand density,and elevation were established.Both height and DBH markedly increased initially and then slowly with tree age,decreased with stand density,and showed unimodal change with elevation.The coupled growth models accounted for72-78%of the variations in tree height,DBH and timber growth.Recommendations for future plantation management are:(1)prolong the rotation to at least 60 years to produce large-diameter,high-quality timber and maintain greater carbon stocks;(2)zone the target functions of stands by elevation;and,(3)reduce stand density for balanced supply of multiple ecosystem services.The growth models developed can predict growth response of larch plantations to density alteration under given ages and elevations,and assist the transformation from traditional management for maximum timber production to site-specific and multifunctional management with longer rotations and moderate tree density.