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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through poro...In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.展开更多
The level of deformation development of surrounding rocks is a vital predictor to evaluate impending coal mine disasters and it is important to establish accurate measurements of the deformed status to ensure coal min...The level of deformation development of surrounding rocks is a vital predictor to evaluate impending coal mine disasters and it is important to establish accurate measurements of the deformed status to ensure coal mine safety. Traditional deformation monitoring methods are mostly based on single parameter, in this paper, multiple approaches are integrated: firstly, both electric and elastic models are established,from which electric field distribution and seismic wave recording are calculated and finally, the resistivity profiles and source position information are determined using inversion methods, from which then the deformation and failure of mine floor are evaluated. According to the inversion results of both electric and seismic field signals, multiple-parameter dynamic monitoring of surrounding rock deformation in deep mine can be performed. The methodology is validated using numerical simulation results which shows that the multi-parameter dynamic monitoring methods have better results for surrounding rock deformation in deep mine monitoring than single parameter methods.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Rock burst in a circular tunnel under high in-situ stress conditions was investigated with a numerical method coupled the rock failure process theory (RFPA) and discontinuous deformation theory (DDA). Some numerical t...Rock burst in a circular tunnel under high in-situ stress conditions was investigated with a numerical method coupled the rock failure process theory (RFPA) and discontinuous deformation theory (DDA). Some numerical tests were carraied out to investigate the failuer patterns of circular tunnel under unloading conditions. Compared the results under loading conditions,the shapes of failure zones are more regular under the unloading conditions. The failure pat-terns in the same type of rock mass are clearly different because of non-homogeneity of the rock material. The extension of cracks shows some predictability with an increasing of in-situ stress. When the homogeneity index of rocks (m) is ei-ther relatively high or low and lateral pressure coefficients (λ) is high,the number of regular shear slide cracks decreases and the probability of a rock burst also becomes lower. Our numerical simulation results show that the stability of sur-face rock and the natural bedding stratification of rock material greatly affect rock bursts. Installing bolts with due dili-gence and suitably can effectively prevent rock bursts. However,it is not effective to control rock bursts by releasing the strain energy with normal pre-boreholes.展开更多
The present study examines the results of the researches related to the gob bulking factor carried out at home and abroad.A mathematical function of a three-dimensional gob bulking factor is described based on a three...The present study examines the results of the researches related to the gob bulking factor carried out at home and abroad.A mathematical function of a three-dimensional gob bulking factor is described based on a three-dimensional gob model.The method of taking value for interstice and permeability ratios is also proposed.The law of air leakage of fully mechanized top coal is researched in this study.The results show that the speed of air flow near the upper and lower crossheadings is higher than that in the central section of the gob at the same distance from the working face.When the amount of air at the working face exceeds a critical amount,the width of the spontaneous combustion zone in the upper and lower crossheadings is also larger than that in the central section.In this situation,the key is preventing the coal left in the upper and lower crossheadings from self-igniting.Reducing the amount of air at the working face can decrease the width of the spontaneous combustion zone,especially the width near the upper and lower crossheadings.This also moves the spontaneous combustion zone in the direction of the working face.It can prevent the coal in the gob from self-igniting by making the coal left in the crossheadings to be inert and by effectively controlling the amount of air at the working face.展开更多
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.展开更多
基金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.
基金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.
基金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 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.
基金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.
基金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.
文摘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.
文摘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.
基金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.
基金Projects 50534090 and 50674090 supported by the National Natural Science Foundation of China2005CB221503 by the National Key Basic ResearchDevelopment Program (973 Program)
文摘In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.
基金financial support from the Fundamental Research Funds for the Central Universities of China (No. 2015QNB19)the financial support from the Open Fund of Key Laboratory of Safety and High-efficiency Coal Mining, Ministry of Education of China (No. JYBSYS2015107)+2 种基金the National Natural Science Foundation of China (Nos. 51404254, 41430317 and U1261202)the China Postdoctoral Science Foundation of China (No. 2014M560465)the Jiangsu Planned Projects for Postdoctoral Research Funds of China (No. 1302050B)
文摘The level of deformation development of surrounding rocks is a vital predictor to evaluate impending coal mine disasters and it is important to establish accurate measurements of the deformed status to ensure coal mine safety. Traditional deformation monitoring methods are mostly based on single parameter, in this paper, multiple approaches are integrated: firstly, both electric and elastic models are established,from which electric field distribution and seismic wave recording are calculated and finally, the resistivity profiles and source position information are determined using inversion methods, from which then the deformation and failure of mine floor are evaluated. According to the inversion results of both electric and seismic field signals, multiple-parameter dynamic monitoring of surrounding rock deformation in deep mine can be performed. The methodology is validated using numerical simulation results which shows that the multi-parameter dynamic monitoring methods have better results for surrounding rock deformation in deep mine monitoring than single parameter methods.
基金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.
基金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.
文摘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.
文摘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.
文摘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.
基金Projects 50639100 supported by the National Natural Science Foundation of China50539100 by the New Century Talents Plan of Education Depart- ment
文摘Rock burst in a circular tunnel under high in-situ stress conditions was investigated with a numerical method coupled the rock failure process theory (RFPA) and discontinuous deformation theory (DDA). Some numerical tests were carraied out to investigate the failuer patterns of circular tunnel under unloading conditions. Compared the results under loading conditions,the shapes of failure zones are more regular under the unloading conditions. The failure pat-terns in the same type of rock mass are clearly different because of non-homogeneity of the rock material. The extension of cracks shows some predictability with an increasing of in-situ stress. When the homogeneity index of rocks (m) is ei-ther relatively high or low and lateral pressure coefficients (λ) is high,the number of regular shear slide cracks decreases and the probability of a rock burst also becomes lower. Our numerical simulation results show that the stability of sur-face rock and the natural bedding stratification of rock material greatly affect rock bursts. Installing bolts with due dili-gence and suitably can effectively prevent rock bursts. However,it is not effective to control rock bursts by releasing the strain energy with normal pre-boreholes.
基金Financial support for this work was obtained from the National Natural Science Foundation of China(No.51074059)
文摘The present study examines the results of the researches related to the gob bulking factor carried out at home and abroad.A mathematical function of a three-dimensional gob bulking factor is described based on a three-dimensional gob model.The method of taking value for interstice and permeability ratios is also proposed.The law of air leakage of fully mechanized top coal is researched in this study.The results show that the speed of air flow near the upper and lower crossheadings is higher than that in the central section of the gob at the same distance from the working face.When the amount of air at the working face exceeds a critical amount,the width of the spontaneous combustion zone in the upper and lower crossheadings is also larger than that in the central section.In this situation,the key is preventing the coal left in the upper and lower crossheadings from self-igniting.Reducing the amount of air at the working face can decrease the width of the spontaneous combustion zone,especially the width near the upper and lower crossheadings.This also moves the spontaneous combustion zone in the direction of the working face.It can prevent the coal in the gob from self-igniting by making the coal left in the crossheadings to be inert and by effectively controlling the amount of air at the working face.
基金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.
