Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods...Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.展开更多
Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experim...Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.展开更多
The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distanc...The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distance flight airframe power.In this paper,a computational model of the scramjet magnetohydrody⁃namic channel is developed and verified by using the commercial software Fluent.It is found that when the mag⁃netic induction intensity is 1,2,3,4 T,the power generation efficiency is 22.5%,22.3%,22.0%,21.5%,and decreases with the increase of the magnetic induction intensity,and the enthalpy extraction rate is 0.026%,0.1%,0.21%,0.34%,and increases with the increase of the magnetic induction intensity.The deceleration ef⁃fect of electromagnetic action on the airflow in the power channel increases with the increase of magnetic induc⁃tion intensity.The stronger the magnetic field intensity,the more obvious the decreasing effect of fluid Mach num⁃ber in the channel.The power generation efficiency decreases as the magnetic induction intensity increases and the enthalpy extraction rate is reversed.As the local currents gathering at inlet and outlet of the power generation area,total temperature and enthalpy along the flow direction do not vary linearly,and there are maximum and minimum values at inlet and outlet.Increasing the number of electrodes can effectively regulate the percentage of Joule heat dissipation,which can improve the power generation efficiency.展开更多
In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommend...In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.展开更多
To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyur...To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing...This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing,logging calculation,and seismic inversion technology,we obtained the current insitu stress characteristics of a single well and rock mechanical parameters.Simultaneously,significant controlling factors of rock mechanical properties were analyzed.Subsequently,by coupling hydraulic fracturing physical experiments with finite element numerical simulation,three different fracturing models were configured:single-cluster,double-cluster,and triple-cluster perforations.Combined with acoustic emission technology,the fracture initiation mode and evolution characteristics during the loading process were determined.The results indicate the following findings:(1)The extension direction and length of the fracture are significantly controlled by the direction of the maximum horizontal principal stress.(2)Areas with poor cementation and compactness exhibit complex fracture morphology,prone to generating network fractures.(3)The interlayer development of fracturing fractures is controlled by the strata occurrence.(4)Increasing the displacement of fracturing fluid enlarges the fracturing fracture length and height.This research provides theoretical support and effective guidance for hydraulic fracturing design in tight oil and gas reservoirs.展开更多
A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radia...A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radiation, and convection. A formula for the calculation of the heat conductivity of a sintering system containing solid phase, liquid phase, and gas phase is given. Due to the continuous movement of the laser beam, a local coordinate system centered on the laser beam is used to simplify the analytical calculation. Assuming that it is approximately a Gaussian laser beam, the heat conduction equation is resolved based on the assumption of the thermal insulating boundary conditions and the fixed thermal physical parameters. The FORTRAN language is employed to compile the program to simulate the temperature field in the direct copper powder sintering process. It shows a good agreement with the preliminary experimental results.[KH3/4D]展开更多
The 3 D non isothermal flow of non Newtonian viscous polymer melt in a co rotating twin screw extruder is modeled. The distributions of the velocity, temperature, pressure and the viscous dissipation in the fl...The 3 D non isothermal flow of non Newtonian viscous polymer melt in a co rotating twin screw extruder is modeled. The distributions of the velocity, temperature, pressure and the viscous dissipation in the flow domain are presented by using a fluid dynamics analysis package (Polyflow). The numerical results show that the temperatures are high in the intermeshing region and on the screw surface, the maximum pressure and the minimum pressure occur in the intermeshing region, and the flow rate is almost proportional to the screw speed.展开更多
Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the...Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.展开更多
In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the sof...In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.展开更多
Aiming at the impaction among granules of non obstructive particle damping(NOPD), the vibration absorption model for vertical impact of granules is established by adopting Hertz contact theory. The numerical simulati...Aiming at the impaction among granules of non obstructive particle damping(NOPD), the vibration absorption model for vertical impact of granules is established by adopting Hertz contact theory. The numerical simulation of the granules movement process is proceeded, and the vibration response of a free free uniform beam is obtained for the case when all granules act on it. Through this method, the effect on vibration absorption of impaction is investigated. The simulational data show that multi gra nule vertical impaction is not sensitive to the movement clearance. The vibration absorption is also very well when the clearance changes within a large range. Therefore, the phenomenon that the vibration magnitude may increase if the clearance in a single impact body is improperly selected will not happen. The effect of vibration suppression in the range of middle and high frequencies(2 500~6 000 Hz) is better than that in the range of low frequency(<2 500 Hz). It indicates that the effect on vibration absorption of multi granule can well restrain the vibration of middle and high frequencies.展开更多
The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investig...The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.展开更多
Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-materia...Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.展开更多
LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation ha...LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.展开更多
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.展开更多
A Double Shear Model(DSM) was used in a numerical simulation on bolted rock joint shearing performance.An entire bolt deformed as the letter'U'under a shear load between two joints.Near the bolt-joint intersec...A Double Shear Model(DSM) was used in a numerical simulation on bolted rock joint shearing performance.An entire bolt deformed as the letter'U'under a shear load between two joints.Near the bolt-joint intersection,the bolt partly deformed as the letter'Z'.There were two critical points along the bolt:one was at the bolt-joint intersection with zero bending moment and the other at the maximum bending moment(plastic hinge) with zero shear stress.The blocks on two sides slid along the bolt as it deformed. A separation area was found between the two joint contact surfaces of the middle rock block and sided block.This area of separation was related to bolt diameter and external forces.We assume that this area is related to the work of external forces.Further research is needed.展开更多
In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, re...In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group （RNG） model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.展开更多
Gas extraction practice has been proven for the clear majority of coal mines in China to be unfavorable using drill holes in the coal seam. Rather, mining-induced fractures in the goaf should be utilized for gas extra...Gas extraction practice has been proven for the clear majority of coal mines in China to be unfavorable using drill holes in the coal seam. Rather, mining-induced fractures in the goaf should be utilized for gas extraction. To study gas migration in mining-induced fractures, one mining face of 10 th Mine in Pingdingshan Coalmine Group in Henan, China, has been selected as the case study for this work. By establishing the mathematical model of gas migration under the influence of coal seam mining, discrete element software UDEC and Multiphysics software COMSOL are employed to model gas migration in mining-induced fractures above the goaf. The results show that as the working face advances, the goaf overburden gradually forms a mining-induced fracture network in the shape of a trapezoid, the size of which increases with the distance of coal face advance. Compared with gas migration in the overburden matrix, the gas flow in the fracture network due to mining is far greater. The largest mining-induced fracture is located at the upper end of the trapezoidal zone, which results in the largest gas flux in the network. When drilling for gas extraction in a mining-induced fracture field, the gas concentration is reduced in the whole region during the process of gas drainage, and the rate of gas concentration drops faster in the fractured zone. It is shown that with gas drainage, the gas flow velocity in the mininginduced fracture network is faster.展开更多
A numerical program is built to simulate the performance of a spark ignited two-stroke free-piston engine coupled with a linear generator. The computational model combines a series of dynamic and thermodynamic equatio...A numerical program is built to simulate the performance of a spark ignited two-stroke free-piston engine coupled with a linear generator. The computational model combines a series of dynamic and thermodynamic equations that are solved simultaneously to predict the performances of the engines. The dynamic analysis performed consists of an evaluation of the frictional force and load force introduced by the generator. The thermodynamic analysis used a single zone model to describe the engine＇ s working cycle which includes intake, scavenging, compression, combustion and expansion, and to evaluate the effect of heat transfer based on the first law of thermodynamics and the ideal gas state equation. Because there is no crankshaft, a time based Wiebe equation was used to express the fraction of fuel burned in the combustion. The calculated results were validated by using the experimental data from another research group. The results indicate that the free-piston generator has some advantages over conventional engines.展开更多
基金National Natural Science Foundation of China(No.11575080)Hunan Provincial Natural Science Foundation of China(No.2022JJ30482)Hunan Provincial Innovation Foundation for Postgraduate(No.QL20220206).
文摘Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.
基金This work was supported by the National Natural Science Foundation of China(No.11575080)the National Natural Science Foundation of Hunan Province,China(No.2022JJ30482)the Hunan Provincial Innovation Foundation for Postgraduates(No.QL20220206).
文摘Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.
文摘The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distance flight airframe power.In this paper,a computational model of the scramjet magnetohydrody⁃namic channel is developed and verified by using the commercial software Fluent.It is found that when the mag⁃netic induction intensity is 1,2,3,4 T,the power generation efficiency is 22.5%,22.3%,22.0%,21.5%,and decreases with the increase of the magnetic induction intensity,and the enthalpy extraction rate is 0.026%,0.1%,0.21%,0.34%,and increases with the increase of the magnetic induction intensity.The deceleration ef⁃fect of electromagnetic action on the airflow in the power channel increases with the increase of magnetic induc⁃tion intensity.The stronger the magnetic field intensity,the more obvious the decreasing effect of fluid Mach num⁃ber in the channel.The power generation efficiency decreases as the magnetic induction intensity increases and the enthalpy extraction rate is reversed.As the local currents gathering at inlet and outlet of the power generation area,total temperature and enthalpy along the flow direction do not vary linearly,and there are maximum and minimum values at inlet and outlet.Increasing the number of electrodes can effectively regulate the percentage of Joule heat dissipation,which can improve the power generation efficiency.
基金the funding supports of the National Key Research and Development Plan,China(Grant No.2022YFC3801800)National Natural Science Foundation of China(Grant Nos.52038010 and 52078368)。
文摘In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.
基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200494)China Postdoctoral Science Foundation(Grant No.2021M701725)+3 种基金Jiangsu Postdoctoral Research Funding Program(Grant No.2021K522C)Fundamental Research Funds for the Central Universities(Grant No.30919011246)National Natural Science Foundation of China(Grant No.52278188)Natural Science Foundation of Jiangsu Province(Grant No.BK20211196)。
文摘To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
基金supported by the Major Scientific and Technological Projects of CNPC under grant ZD2019-183-006the National Science and Technology Major Project of China(2016ZX05014002-006)the National Natural Science Foundation of China(42072234)。
文摘This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing,logging calculation,and seismic inversion technology,we obtained the current insitu stress characteristics of a single well and rock mechanical parameters.Simultaneously,significant controlling factors of rock mechanical properties were analyzed.Subsequently,by coupling hydraulic fracturing physical experiments with finite element numerical simulation,three different fracturing models were configured:single-cluster,double-cluster,and triple-cluster perforations.Combined with acoustic emission technology,the fracture initiation mode and evolution characteristics during the loading process were determined.The results indicate the following findings:(1)The extension direction and length of the fracture are significantly controlled by the direction of the maximum horizontal principal stress.(2)Areas with poor cementation and compactness exhibit complex fracture morphology,prone to generating network fractures.(3)The interlayer development of fracturing fractures is controlled by the strata occurrence.(4)Increasing the displacement of fracturing fluid enlarges the fracturing fracture length and height.This research provides theoretical support and effective guidance for hydraulic fracturing design in tight oil and gas reservoirs.
文摘A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radiation, and convection. A formula for the calculation of the heat conductivity of a sintering system containing solid phase, liquid phase, and gas phase is given. Due to the continuous movement of the laser beam, a local coordinate system centered on the laser beam is used to simplify the analytical calculation. Assuming that it is approximately a Gaussian laser beam, the heat conduction equation is resolved based on the assumption of the thermal insulating boundary conditions and the fixed thermal physical parameters. The FORTRAN language is employed to compile the program to simulate the temperature field in the direct copper powder sintering process. It shows a good agreement with the preliminary experimental results.[KH3/4D]
文摘The 3 D non isothermal flow of non Newtonian viscous polymer melt in a co rotating twin screw extruder is modeled. The distributions of the velocity, temperature, pressure and the viscous dissipation in the flow domain are presented by using a fluid dynamics analysis package (Polyflow). The numerical results show that the temperatures are high in the intermeshing region and on the screw surface, the maximum pressure and the minimum pressure occur in the intermeshing region, and the flow rate is almost proportional to the screw speed.
基金Supported by the National Natural Defense Basic Scientific Research Program of China(A262006-1288)the Key Disciplines Program of Shanghai Municipal Commission of Education(J50501)~~
文摘Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.
文摘In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.
文摘Aiming at the impaction among granules of non obstructive particle damping(NOPD), the vibration absorption model for vertical impact of granules is established by adopting Hertz contact theory. The numerical simulation of the granules movement process is proceeded, and the vibration response of a free free uniform beam is obtained for the case when all granules act on it. Through this method, the effect on vibration absorption of impaction is investigated. The simulational data show that multi gra nule vertical impaction is not sensitive to the movement clearance. The vibration absorption is also very well when the clearance changes within a large range. Therefore, the phenomenon that the vibration magnitude may increase if the clearance in a single impact body is improperly selected will not happen. The effect of vibration suppression in the range of middle and high frequencies(2 500~6 000 Hz) is better than that in the range of low frequency(<2 500 Hz). It indicates that the effect on vibration absorption of multi granule can well restrain the vibration of middle and high frequencies.
基金Supported by the National Natural Science Foundation of China(10672052)the Advanced TalentStart-Up Foundation of Jiangsu University(08JDG018)~~
文摘The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.
基金Supported by the National Natural Science Foundation of China(10476011)~~
文摘Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.
文摘LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.
基金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.
基金Financial support from the National Natural Science Foundation of China(No.50978251)the National Government Building High-Level University Graduate Programs of the China Scholarship Council
文摘A Double Shear Model(DSM) was used in a numerical simulation on bolted rock joint shearing performance.An entire bolt deformed as the letter'U'under a shear load between two joints.Near the bolt-joint intersection,the bolt partly deformed as the letter'Z'.There were two critical points along the bolt:one was at the bolt-joint intersection with zero bending moment and the other at the maximum bending moment(plastic hinge) with zero shear stress.The blocks on two sides slid along the bolt as it deformed. A separation area was found between the two joint contact surfaces of the middle rock block and sided block.This area of separation was related to bolt diameter and external forces.We assume that this area is related to the work of external forces.Further research is needed.
基金Project supported by the National Natural Science Foundation of China （Grant Nos 40233030, 40405004 and 40405014） and the Special Program of the Scientific and Social Practices for Graduate Students in Chinese Academy of Sciences, China.
文摘In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group （RNG） model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.
基金financially supported by the National Key Research and Development Program (No. 2016YFC0801402)the National Natural Science Foundation of China (No. 51374236)Chongqing Research Program of Basic Research and Frontier Technology of China (No. cstc2015jcyj BX0076)
文摘Gas extraction practice has been proven for the clear majority of coal mines in China to be unfavorable using drill holes in the coal seam. Rather, mining-induced fractures in the goaf should be utilized for gas extraction. To study gas migration in mining-induced fractures, one mining face of 10 th Mine in Pingdingshan Coalmine Group in Henan, China, has been selected as the case study for this work. By establishing the mathematical model of gas migration under the influence of coal seam mining, discrete element software UDEC and Multiphysics software COMSOL are employed to model gas migration in mining-induced fractures above the goaf. The results show that as the working face advances, the goaf overburden gradually forms a mining-induced fracture network in the shape of a trapezoid, the size of which increases with the distance of coal face advance. Compared with gas migration in the overburden matrix, the gas flow in the fracture network due to mining is far greater. The largest mining-induced fracture is located at the upper end of the trapezoidal zone, which results in the largest gas flux in the network. When drilling for gas extraction in a mining-induced fracture field, the gas concentration is reduced in the whole region during the process of gas drainage, and the rate of gas concentration drops faster in the fractured zone. It is shown that with gas drainage, the gas flow velocity in the mininginduced fracture network is faster.
文摘A numerical program is built to simulate the performance of a spark ignited two-stroke free-piston engine coupled with a linear generator. The computational model combines a series of dynamic and thermodynamic equations that are solved simultaneously to predict the performances of the engines. The dynamic analysis performed consists of an evaluation of the frictional force and load force introduced by the generator. The thermodynamic analysis used a single zone model to describe the engine＇ s working cycle which includes intake, scavenging, compression, combustion and expansion, and to evaluate the effect of heat transfer based on the first law of thermodynamics and the ideal gas state equation. Because there is no crankshaft, a time based Wiebe equation was used to express the fraction of fuel burned in the combustion. The calculated results were validated by using the experimental data from another research group. The results indicate that the free-piston generator has some advantages over conventional engines.
