Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately c...Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.展开更多
To investigate the reason for Mn segregation in TC2 titanium alloy bars,a multiphysics-coupled mathematical model was established using the BMPS-VAR numerical simulation software,incorporating electro-magnetic,thermal...To investigate the reason for Mn segregation in TC2 titanium alloy bars,a multiphysics-coupled mathematical model was established using the BMPS-VAR numerical simulation software,incorporating electro-magnetic,thermal,and flow fields.Numerical simulation was performed to analyze the dynamic evolution of varying-mass electrode fragments during the vacuum arc remelting(VAR)of a∅508-mm TC2 titanium alloy in-got.The results indicate that Mn segregation caused by 15-kg electrode fragmentation during the VAR process of a TC2 titanium alloy ingot corresponds to the segregation observed in the TC2 titanium alloy bar.The numerical simulation of the VAR process provides effective result prediction and technical support for solving practical problems in smelting.展开更多
Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic th...Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.展开更多
Since the General Channel designed for the South-to-North Water Transfer Project in China has to cross many rivers and streams flowing from west to east, there are potentially serious effects additional flooding on th...Since the General Channel designed for the South-to-North Water Transfer Project in China has to cross many rivers and streams flowing from west to east, there are potentially serious effects additional flooding on the western side of the project alignment. Therefore, a 2-D numerical model for forecasting basin flood disasters was established and verified using historical flood data. The model was applied to researching the interaction between the proposed Project and flooding events for 5 streams in the Anyang River reach as a representative case study. Simulated results indicate that the model could correctly forecast the flood, submerged area and depths, and water surface elevations along the left side of the channel. The discharge capacity and location of hydraulic structures in the transfer canal alignment were analyzed. Then adjustments to the dimensions and positioning of proposed hydraulic structures were recommended at intersections, especially the addition of a channel to transfer flood water from one stream to another, which can effectively limit the sluice and protect the Anyang City from flooding.展开更多
In this paper, the erosion-resisting coefficient was introduced to computebed deformation in a crush bedrock river. In the case of crush bedrock, there has been no propercontrol equation to describe bed stability, whi...In this paper, the erosion-resisting coefficient was introduced to computebed deformation in a crush bedrock river. In the case of crush bedrock, there has been no propercontrol equation to describe bed stability, which leads to difficulty in calculation of the beddeformation with conventional methods. The data from field survey were used to give thee-rosion-resisting capability with an appropriate coefficient. After the determination oflongitudinal distribution expressed by polynomial regression and transversal distribution expressedby normal distribution function, the plane distribution of erosion-resisting coefficient in a crushbedrock river was obtained. With the computational results from a 2-D horizontal flow mathematicalmodel, the erosion-resisting coefficient and controlling condition of local stability were employedto compute the values of bed deformation when riverbed is stable. The above method was applied in acase study, and the computational results of flow and bed deformations are in good a-greement withphysical model test data.展开更多
The assessment of the radiological impact of the liquid discharges from nuclear power plants is a major issue for the envi- ronmental protection. In this study, a numerical model for the radionuclide transport in the ...The assessment of the radiological impact of the liquid discharges from nuclear power plants is a major issue for the envi- ronmental protection. In this study, a numerical model for the radionuclide transport in the aquatic environment is built, based on the hydrodynamic equations, including the complete set of Saint-Venant equations, the sediment transport equations, with consideration of several different particle sizes and the deposition and erosion of the suspended sediments, and the radionuclide transport equations The exchanges of radionuclides between water, suspended matter and bed sediments are described in terms of kinetic transfer coe- fficients. The model is used to simulate the transport of the radionuclides discharged from a planned nuclear power plant project to be sited along the lower Yangtze River. From the model results, one may see the detailed temporal-spatial evolution of the radio- nuclide contamination in the solution, in the suspended matter and in the bed sediments. The model can be used as a basic tool for studying the environmental impacts of the liquid discharges from nuclear facilities on a river system.展开更多
Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vor...Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson＇s work, the maximum amplitude of the cylinder with low mass ratio m＊=2.6 can reach as large as 1.5D to be called as the ＂super-upper branch＂, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.展开更多
Based on the fully three-dimensional(3-D)and two-dimensional(2-D)comprehensive CFD(Computational Fluid Dynamics)combustion models for a circulating fluidized bed boiler,a simplified 3-D computational domain considerin...Based on the fully three-dimensional(3-D)and two-dimensional(2-D)comprehensive CFD(Computational Fluid Dynamics)combustion models for a circulating fluidized bed boiler,a simplified 3-D computational domain considering the corrections of furnace side wall openings is proposed.It aims to compensate for the deficiencies of the large amount of computation in the fully 3-D model and improve the air and gas flow treatments at the openings in the simplified 2-D model.Three different computational domains,named as the fully 3-D model,simplified 3-D model and 2-D model,were implemented to perform a comparative CFD analysis in an ultra-supercritical circulating fluidized bed boiler including the hydrodynamics,penetration depth of secondary air,temperature and species distribution.The simulation results computed by the simplified3-D model yield better agreement with the fully 3-D simulation results than those of the 2-D model.The simplified 3-D model is recommended as an alternative computational domain for the conventional 2-D model in the numerical simulation of large-scale circulating fluidized bed boiler.展开更多
The construction and operation of sulfur-containing gas storage are often more difficult than a non-sulfur storage facility due to the need to prevent environmental contamination from H_(2)S leaks,as well as the corro...The construction and operation of sulfur-containing gas storage are often more difficult than a non-sulfur storage facility due to the need to prevent environmental contamination from H_(2)S leaks,as well as the corrosive effects of H_(2)S on production facilities.Rapid elutriation of H_(2)S from the reservoir during the construction of the gas storage is an effective way to avoid these problems.However,the existing H_(2)S elutriation method has low efficiency and high economic cost,which limits the development of reconstructed gas storage of sulfur-containing gas reservoirs.To improve the efficiency of H_(2)S elutriation in sulfur-containing gas reservoirs and enhance the economic benefits,a numerical simulation model of multiphase flow components was established to study the migration law of H_(2)S in the multi-cycle operation of gas storage.Based on the H_(2)S migrate law,the displacement H_(2)S elutriation method was developed,and the elutriation mechanism and elutriation efficiency of the two methods were compared and analyzed.In addition,the main controlling factors affecting the H_(2)S elutriation efficiency were investigated,and the H_(2)S elutriation scheme of H gas storage was optimized.The results indicate that H_(2)S migrates between near-well and far-well regions under pressure differentials.The traditional H_(2)S elutriation method relies on concentration gradient diffusion,whereas the displacement elutriation approach leverages pressure differentials with higher H_(2)S elutriation efficiency.For the displacement elutriation method,higher reservoir permeability enhances the peak-shaving capacity of the gas storage but has a minor impact on H_(2)S elutriation when the formation permeability is between 30 and 100 mD.The elutriation efficiency is significantly higher when wells are drilled in the high structural parts of the reservoir compared to the low structural parts.Longer displacement elutriation time within a cycle improves H_(2)S elutriation efficiency but reduces the working gas volume of the storage.Therefore,the optimal displacement time for H gas storage is 60 days.An optimized H_(2)S elutriation scheme enabled the working gas to meet the national first-class natural gas standard within 10 cycles.This study elucidates H_(2)S migration patterns,H_(2)S elutriation mechanisms,and key influence factors on H_(2)S elutriation efficiency,offering valuable technical insights for sour gas storage operations.展开更多
Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately ...Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.展开更多
The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-E...The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-ECBM technology should be perused. First, the coupling equations of the di usion–adsorption–seepage–heat transfer fields of gas are established. The displacement processes under di erent pressures and temperatures are simulated by COMSOL. Finally, the displacement effects, a comparison of the CO_2 storage capacity with the CH_4 output and the e ective influencing radius of CO_2 injection are analyzed and discussed. The results show that(1) the displacement pressure and temperature are two key factors influencing the CH_4 output and the CO_2 storage capacity, and the increase in the CO_2 storage capacity is more sensitive to temperature and pressure than the CH_4 output.(2) The gas flow direction is from the injection hole to the discharge hole during the displacement process, and the regions with high velocity are concentrated at the injection hole and the discharge hole.(3) A reduction in the CH_4 concentration and an increase in the CO_2 concentration are obvious during the displacement process.(4) The e ective influencing radius of injecting CO_2 with heat increases with the increase in time and pressure. The relationship between the e ective influencing radius and the injection time of CO_2 has a power exponential function, and there is a linear relationship between the functional coe cient and the injection pressure of CO_2. This numerical simulation study on enhancing CBM recovery by injecting CO_2 with heat can further promote the implementation of CO_2-ECBM project in deep coal seams.展开更多
Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surroundi...Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.展开更多
2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock ...2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock mass and controlling the stability of surrounding rock. To achieve an accurate simulation of bolted joint shear tests, we have studied a numerical simulation method that takes into account the 2G-NPR bolt's tensile–shear fracture criterion. Firstly, the indoor experimental study on the tensile–shear mechanical properties of 2G-NPR bolt is carried out to explore its mechanical properties under diferent tensile–shear angles, and the fracture criterion of 2G-NPR bolt considering the tensile–shear angle is established. Then, a three-dimensional numerical simulation method considering the tensile–shear mechanical constitutive and fracture criterion of 2G-NPR bolt, the elastoplastic mechanical behavior of surrounding rock and the damage and deterioration of grouting body is proposed. The feasibility and accuracy of the method are verifed by comparing with the indoor shear test results of 2G-NPR bolt anchorage joints. Finally, based on the numerical simulation results, the deformation and stress of the bolt, the distribution of the plastic zone of the rock mass, the stress distribution and the damage of the grouting body are analyzed in detail. The research results can provide a good reference value for the practical engineering application and shear mechanical performance analysis of 2G-NPR bolt.展开更多
Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to ex...Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to explore the following issues of CGS:fluid and heat flow,solute transport;stresses,displacements and rock failures related to geo-mechanical effects;equilibrium and kinetic chemical reactions;chemical damage to mechanical properties of the rock.Then,the coupling program was applied to the Ordos CGS Project to study the formation response under the multifield interaction caused by CO_(2)injection.The simulation results show that the mechanical process dominates the short CO_(2)injection period.Specifically,the formation’s permeability near the injection well increases by 43%,due to the reduction of effective stress,which significantly promotes the lateral migration of CO_(2).When the injection rate exceeds 0.15 million tons per year,the cohesion of the reservoir rock is not enough to resist the shear force inside the rock and rock failure may occur.During the subsequent long-term sequestration period(200 years),the influence of mineral reactions gradually increases.Due to calcite dissolution,the shear modulus of caprock is predicted to decrease by 7.6%,which will to some extent increase the risk of rock failure.展开更多
The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limit...The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limited understanding of the resources quantity and storage gained from scientific researches.In this study,using a typical cross section across the basin and taking into account its geothermal and geological conditions,a new waterheat coupled model was built and associated modelling was done by the software TOUGH2.During modelling process,the accuracy and applicability of the model was confirmed through the calibration of relevant parameters for modelling the heat and water transport and the formation of geothermal reservoir across the basin,with particular focus on the Neogene geothermal field.Results show that the groundwater that flows from the basin margins to the center is heated by the Neogene and Paleogene sedimentary rocks with high geothermal gradients.Since the east-west extending fault F1 is conductive,it acts as preferential flow paths which on one hand provide additional and rapid flows to the thermal reservoir;and on the other hand,cool down the thermal water to a certain extent due to the infiltration of shallower water sources in the vicinity of the fault.Furthermore,the estimated geothermal resources quantity is close to that of previous studies.In comparison with the Paleogene rock formations,the Neogene geothermal reservoir shows a better nature in terms of water content,aquifer permeability and resources exploitability,although the resource quantity of the Paleogene reservoir is considerable.展开更多
Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster managem...Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.展开更多
Gas breakthrough pressure is a key parameter to evaluate the sealing capacity of caprock,and it also plays important roles in safety and capacity of CO_(2)geological storage.Based on the published experimental results...Gas breakthrough pressure is a key parameter to evaluate the sealing capacity of caprock,and it also plays important roles in safety and capacity of CO_(2)geological storage.Based on the published experimental results,we present numerical simulations on CO_(2)breakthrough pressure in unsaturated low-permeability rock under 9 multiple P-T conditions(which can keep CO_(2)in gaseous,liquid and supercritical states)and thus,a numerical method which can be used to accurately predict CO_(2)breakthrough pressure on rock-core scale is proposed.The simulation results show that CO_(2)breakthrough pressure and breakthrough time are exponential correlated with P-T conditions.Meanwhile,pressure has stronger effects on experimental results than that of temperature.Moreover,we performed sensitivity studies on the pore distribution indexλ(0.6,0.7,0.8,and 0.9)in van Genuchten-Muale model.Results show that with the increase ofλ,CO_(2)breakthrough pressure and breakthrough time both show decreasing trends.In other words,the larger the value ofλis,the better the permeability of the caprock is,and the worse the CO_(2)sealing capacity is.The numerical method established in this study can provide an important reference for the prediction of gas breakthrough pressure on rock-core scale and for related numerical studies.展开更多
According to the influence of the combination of short-distance coal seam group on mining roadway, using numerical simulation software FLAG2D to draw the abutment pressure distribution ahead the working face and the a...According to the influence of the combination of short-distance coal seam group on mining roadway, using numerical simulation software FLAG2D to draw the abutment pressure distribution ahead the working face and the area of influence in fully-mechanized mining conditions, the variation rules of surrounding rock supporting pressure of floor roadway and the deformation rules were summarized. GYS-300 anchor dynamometer was used to measure the roadway surface displacement, and the conclusions of numerical simulation were verified.展开更多
Hot compression tests of 3Cr2NiMnMo steel were performed at temperatures in the range of 850 to 1 100 ℃ and with strain rates of 10 ^-2 s^- 1 to 1 s ^-1. Both the constitutive equations and the hot deformation activa...Hot compression tests of 3Cr2NiMnMo steel were performed at temperatures in the range of 850 to 1 100 ℃ and with strain rates of 10 ^-2 s^- 1 to 1 s ^-1. Both the constitutive equations and the hot deformation activation energy were derived from the correlativity of flow stress, strain rate and temperature. The mathematical models of the dynamic recrystaIiization of 3Cr2NiMnMo steel, which inelude the dynamic recrystallization kinetics model and the crystalliza- tion grain size model, are based on Avrami＇s law and the results of thermosimulation experiments. By integrating de- rived dynamic recrystallization models with the thermal mechanical coupled finite element method, the microstruc ture evolution in hot compressive deformation was simulated. The distribution of dynamic recrystallization grains and grain sizes were determined through a comparison of the simulation results with the experimental results. The distri- bution of strain and dynamic recrystallization grain is also discussed. The similarity between the experimental results and the simulated results indicates that the derived dynamic recrystallization models can be applied effectively to pre diet and analyze the microstructure evolution in hot deformed 3Cr2NiMnMo steel.展开更多
In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubi...In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.52404155 and 52304111)State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining&Technology,Beijing(Grant No.XD2024006).
文摘Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.
文摘To investigate the reason for Mn segregation in TC2 titanium alloy bars,a multiphysics-coupled mathematical model was established using the BMPS-VAR numerical simulation software,incorporating electro-magnetic,thermal,and flow fields.Numerical simulation was performed to analyze the dynamic evolution of varying-mass electrode fragments during the vacuum arc remelting(VAR)of a∅508-mm TC2 titanium alloy in-got.The results indicate that Mn segregation caused by 15-kg electrode fragmentation during the VAR process of a TC2 titanium alloy ingot corresponds to the segregation observed in the TC2 titanium alloy bar.The numerical simulation of the VAR process provides effective result prediction and technical support for solving practical problems in smelting.
基金The National Key R&D Program of China under contract No.2023YFE0126300the National Natural Science Foundation of China under contract Nos 42066002 and U20A2099.
文摘Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.
基金the National Natural Science Foundation of China (Grant No. 50779019).
文摘Since the General Channel designed for the South-to-North Water Transfer Project in China has to cross many rivers and streams flowing from west to east, there are potentially serious effects additional flooding on the western side of the project alignment. Therefore, a 2-D numerical model for forecasting basin flood disasters was established and verified using historical flood data. The model was applied to researching the interaction between the proposed Project and flooding events for 5 streams in the Anyang River reach as a representative case study. Simulated results indicate that the model could correctly forecast the flood, submerged area and depths, and water surface elevations along the left side of the channel. The discharge capacity and location of hydraulic structures in the transfer canal alignment were analyzed. Then adjustments to the dimensions and positioning of proposed hydraulic structures were recommended at intersections, especially the addition of a channel to transfer flood water from one stream to another, which can effectively limit the sluice and protect the Anyang City from flooding.
文摘In this paper, the erosion-resisting coefficient was introduced to computebed deformation in a crush bedrock river. In the case of crush bedrock, there has been no propercontrol equation to describe bed stability, which leads to difficulty in calculation of the beddeformation with conventional methods. The data from field survey were used to give thee-rosion-resisting capability with an appropriate coefficient. After the determination oflongitudinal distribution expressed by polynomial regression and transversal distribution expressedby normal distribution function, the plane distribution of erosion-resisting coefficient in a crushbedrock river was obtained. With the computational results from a 2-D horizontal flow mathematicalmodel, the erosion-resisting coefficient and controlling condition of local stability were employedto compute the values of bed deformation when riverbed is stable. The above method was applied in acase study, and the computational results of flow and bed deformations are in good a-greement withphysical model test data.
文摘The assessment of the radiological impact of the liquid discharges from nuclear power plants is a major issue for the envi- ronmental protection. In this study, a numerical model for the radionuclide transport in the aquatic environment is built, based on the hydrodynamic equations, including the complete set of Saint-Venant equations, the sediment transport equations, with consideration of several different particle sizes and the deposition and erosion of the suspended sediments, and the radionuclide transport equations The exchanges of radionuclides between water, suspended matter and bed sediments are described in terms of kinetic transfer coe- fficients. The model is used to simulate the transport of the radionuclides discharged from a planned nuclear power plant project to be sited along the lower Yangtze River. From the model results, one may see the detailed temporal-spatial evolution of the radio- nuclide contamination in the solution, in the suspended matter and in the bed sediments. The model can be used as a basic tool for studying the environmental impacts of the liquid discharges from nuclear facilities on a river system.
基金financially supported by the National Natural Science Foundation of China(Grant No.51509045)
文摘Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson＇s work, the maximum amplitude of the cylinder with low mass ratio m＊=2.6 can reach as large as 1.5D to be called as the ＂super-upper branch＂, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.
基金by the Key Project of the National Research Program of China(Grant No.2020YFB0606201)。
文摘Based on the fully three-dimensional(3-D)and two-dimensional(2-D)comprehensive CFD(Computational Fluid Dynamics)combustion models for a circulating fluidized bed boiler,a simplified 3-D computational domain considering the corrections of furnace side wall openings is proposed.It aims to compensate for the deficiencies of the large amount of computation in the fully 3-D model and improve the air and gas flow treatments at the openings in the simplified 2-D model.Three different computational domains,named as the fully 3-D model,simplified 3-D model and 2-D model,were implemented to perform a comparative CFD analysis in an ultra-supercritical circulating fluidized bed boiler including the hydrodynamics,penetration depth of secondary air,temperature and species distribution.The simulation results computed by the simplified3-D model yield better agreement with the fully 3-D simulation results than those of the 2-D model.The simplified 3-D model is recommended as an alternative computational domain for the conventional 2-D model in the numerical simulation of large-scale circulating fluidized bed boiler.
基金supported by the Science and Technology Research Program of Chongqing Municipal Education Commission(KJQN202401501,KJZD-M202401501).
文摘The construction and operation of sulfur-containing gas storage are often more difficult than a non-sulfur storage facility due to the need to prevent environmental contamination from H_(2)S leaks,as well as the corrosive effects of H_(2)S on production facilities.Rapid elutriation of H_(2)S from the reservoir during the construction of the gas storage is an effective way to avoid these problems.However,the existing H_(2)S elutriation method has low efficiency and high economic cost,which limits the development of reconstructed gas storage of sulfur-containing gas reservoirs.To improve the efficiency of H_(2)S elutriation in sulfur-containing gas reservoirs and enhance the economic benefits,a numerical simulation model of multiphase flow components was established to study the migration law of H_(2)S in the multi-cycle operation of gas storage.Based on the H_(2)S migrate law,the displacement H_(2)S elutriation method was developed,and the elutriation mechanism and elutriation efficiency of the two methods were compared and analyzed.In addition,the main controlling factors affecting the H_(2)S elutriation efficiency were investigated,and the H_(2)S elutriation scheme of H gas storage was optimized.The results indicate that H_(2)S migrates between near-well and far-well regions under pressure differentials.The traditional H_(2)S elutriation method relies on concentration gradient diffusion,whereas the displacement elutriation approach leverages pressure differentials with higher H_(2)S elutriation efficiency.For the displacement elutriation method,higher reservoir permeability enhances the peak-shaving capacity of the gas storage but has a minor impact on H_(2)S elutriation when the formation permeability is between 30 and 100 mD.The elutriation efficiency is significantly higher when wells are drilled in the high structural parts of the reservoir compared to the low structural parts.Longer displacement elutriation time within a cycle improves H_(2)S elutriation efficiency but reduces the working gas volume of the storage.Therefore,the optimal displacement time for H gas storage is 60 days.An optimized H_(2)S elutriation scheme enabled the working gas to meet the national first-class natural gas standard within 10 cycles.This study elucidates H_(2)S migration patterns,H_(2)S elutriation mechanisms,and key influence factors on H_(2)S elutriation efficiency,offering valuable technical insights for sour gas storage operations.
基金supported by National Natural Science Foundation of China(Nos.41204077,41372290,41572244,51034003,51174210,and 51304126)natural science foundation of Shandong Province(Nos.ZR2011EEZ002 and ZR2013EEQ019)State Key Research Development Program of China(No.2016YFC0600708-3)
文摘Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.
基金financially supported by the National Natural Science Foundation of China(No.41330638)
文摘The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-ECBM technology should be perused. First, the coupling equations of the di usion–adsorption–seepage–heat transfer fields of gas are established. The displacement processes under di erent pressures and temperatures are simulated by COMSOL. Finally, the displacement effects, a comparison of the CO_2 storage capacity with the CH_4 output and the e ective influencing radius of CO_2 injection are analyzed and discussed. The results show that(1) the displacement pressure and temperature are two key factors influencing the CH_4 output and the CO_2 storage capacity, and the increase in the CO_2 storage capacity is more sensitive to temperature and pressure than the CH_4 output.(2) The gas flow direction is from the injection hole to the discharge hole during the displacement process, and the regions with high velocity are concentrated at the injection hole and the discharge hole.(3) A reduction in the CH_4 concentration and an increase in the CO_2 concentration are obvious during the displacement process.(4) The e ective influencing radius of injecting CO_2 with heat increases with the increase in time and pressure. The relationship between the e ective influencing radius and the injection time of CO_2 has a power exponential function, and there is a linear relationship between the functional coe cient and the injection pressure of CO_2. This numerical simulation study on enhancing CBM recovery by injecting CO_2 with heat can further promote the implementation of CO_2-ECBM project in deep coal seams.
基金supported by the National Natural Science Foundation of China(NSFC)(41941018)the Program of China Scholarship Council(202106430031)。
文摘Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.
基金supported by the National Natural Science Foundation of China(NSFC)(41941018).
文摘2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock mass and controlling the stability of surrounding rock. To achieve an accurate simulation of bolted joint shear tests, we have studied a numerical simulation method that takes into account the 2G-NPR bolt's tensile–shear fracture criterion. Firstly, the indoor experimental study on the tensile–shear mechanical properties of 2G-NPR bolt is carried out to explore its mechanical properties under diferent tensile–shear angles, and the fracture criterion of 2G-NPR bolt considering the tensile–shear angle is established. Then, a three-dimensional numerical simulation method considering the tensile–shear mechanical constitutive and fracture criterion of 2G-NPR bolt, the elastoplastic mechanical behavior of surrounding rock and the damage and deterioration of grouting body is proposed. The feasibility and accuracy of the method are verifed by comparing with the indoor shear test results of 2G-NPR bolt anchorage joints. Finally, based on the numerical simulation results, the deformation and stress of the bolt, the distribution of the plastic zone of the rock mass, the stress distribution and the damage of the grouting body are analyzed in detail. The research results can provide a good reference value for the practical engineering application and shear mechanical performance analysis of 2G-NPR bolt.
基金jointly supported by the National Natural Science Foundation of China(Grant No.42141013)the CNPC Innovation Found(2021 D002-1102)the China Postdoctoral Science Foundation(No.2021M701379)。
文摘Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to explore the following issues of CGS:fluid and heat flow,solute transport;stresses,displacements and rock failures related to geo-mechanical effects;equilibrium and kinetic chemical reactions;chemical damage to mechanical properties of the rock.Then,the coupling program was applied to the Ordos CGS Project to study the formation response under the multifield interaction caused by CO_(2)injection.The simulation results show that the mechanical process dominates the short CO_(2)injection period.Specifically,the formation’s permeability near the injection well increases by 43%,due to the reduction of effective stress,which significantly promotes the lateral migration of CO_(2).When the injection rate exceeds 0.15 million tons per year,the cohesion of the reservoir rock is not enough to resist the shear force inside the rock and rock failure may occur.During the subsequent long-term sequestration period(200 years),the influence of mineral reactions gradually increases.Due to calcite dissolution,the shear modulus of caprock is predicted to decrease by 7.6%,which will to some extent increase the risk of rock failure.
基金This work was supported by grants from National Natural Science Foundation of China(No.41402231)geothermal survey project of China Geological Survey(No.DD20190128)+1 种基金Natural Science Foundation of Hebei Province China(No.D2019330003)S&T Program of Hebei China(No.20374201D).
文摘The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limited understanding of the resources quantity and storage gained from scientific researches.In this study,using a typical cross section across the basin and taking into account its geothermal and geological conditions,a new waterheat coupled model was built and associated modelling was done by the software TOUGH2.During modelling process,the accuracy and applicability of the model was confirmed through the calibration of relevant parameters for modelling the heat and water transport and the formation of geothermal reservoir across the basin,with particular focus on the Neogene geothermal field.Results show that the groundwater that flows from the basin margins to the center is heated by the Neogene and Paleogene sedimentary rocks with high geothermal gradients.Since the east-west extending fault F1 is conductive,it acts as preferential flow paths which on one hand provide additional and rapid flows to the thermal reservoir;and on the other hand,cool down the thermal water to a certain extent due to the infiltration of shallower water sources in the vicinity of the fault.Furthermore,the estimated geothermal resources quantity is close to that of previous studies.In comparison with the Paleogene rock formations,the Neogene geothermal reservoir shows a better nature in terms of water content,aquifer permeability and resources exploitability,although the resource quantity of the Paleogene reservoir is considerable.
基金financially supported by Department of Space,India(Grant No.ISRO/RES/4/663/18-19)。
文摘Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.
基金supported by Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.SKLGME021010)funded by the National Natural Science Foundation of China(Grant No.41702251 and 42141010)the MOE Key Laboratory of Groundwater Circulation and Environmental Evolution。
文摘Gas breakthrough pressure is a key parameter to evaluate the sealing capacity of caprock,and it also plays important roles in safety and capacity of CO_(2)geological storage.Based on the published experimental results,we present numerical simulations on CO_(2)breakthrough pressure in unsaturated low-permeability rock under 9 multiple P-T conditions(which can keep CO_(2)in gaseous,liquid and supercritical states)and thus,a numerical method which can be used to accurately predict CO_(2)breakthrough pressure on rock-core scale is proposed.The simulation results show that CO_(2)breakthrough pressure and breakthrough time are exponential correlated with P-T conditions.Meanwhile,pressure has stronger effects on experimental results than that of temperature.Moreover,we performed sensitivity studies on the pore distribution indexλ(0.6,0.7,0.8,and 0.9)in van Genuchten-Muale model.Results show that with the increase ofλ,CO_(2)breakthrough pressure and breakthrough time both show decreasing trends.In other words,the larger the value ofλis,the better the permeability of the caprock is,and the worse the CO_(2)sealing capacity is.The numerical method established in this study can provide an important reference for the prediction of gas breakthrough pressure on rock-core scale and for related numerical studies.
文摘According to the influence of the combination of short-distance coal seam group on mining roadway, using numerical simulation software FLAG2D to draw the abutment pressure distribution ahead the working face and the area of influence in fully-mechanized mining conditions, the variation rules of surrounding rock supporting pressure of floor roadway and the deformation rules were summarized. GYS-300 anchor dynamometer was used to measure the roadway surface displacement, and the conclusions of numerical simulation were verified.
基金Sponsored by Innovation Program of Shanghai Municipal Education Commission of China(12ZZ183)Shanghai IndustryUniversity Coorperation Projects of China(12QY20)
文摘Hot compression tests of 3Cr2NiMnMo steel were performed at temperatures in the range of 850 to 1 100 ℃ and with strain rates of 10 ^-2 s^- 1 to 1 s ^-1. Both the constitutive equations and the hot deformation activation energy were derived from the correlativity of flow stress, strain rate and temperature. The mathematical models of the dynamic recrystaIiization of 3Cr2NiMnMo steel, which inelude the dynamic recrystallization kinetics model and the crystalliza- tion grain size model, are based on Avrami＇s law and the results of thermosimulation experiments. By integrating de- rived dynamic recrystallization models with the thermal mechanical coupled finite element method, the microstruc ture evolution in hot compressive deformation was simulated. The distribution of dynamic recrystallization grains and grain sizes were determined through a comparison of the simulation results with the experimental results. The distri- bution of strain and dynamic recrystallization grain is also discussed. The similarity between the experimental results and the simulated results indicates that the derived dynamic recrystallization models can be applied effectively to pre diet and analyze the microstructure evolution in hot deformed 3Cr2NiMnMo steel.
文摘In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.
