During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configura...During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.展开更多
The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process ...The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process greatly influences the evolution of confined blast loading and the subsequent structural response,which is crucial in confined blast scenarios.Given the complex nature of the reaction process,accurate analysis of the afterburning effect remains challenging.Previous studies have either overlooked the mechanisms of detonation product combustion or failed to provide experimental validation.This study introduces a three-dimensional model to effectively characterize the combustion of detonation products.The model integrates chemical reaction source terms into the governing equations to consider the combustion processes.Numerical simulations and experimental tests were conducted to analyze the combustion and energy release from the detonation products of fuel-rich explosives in confined spaces.Approximately 50%of the energy was released during the combustion of detonation products in a confined TNT explosion.Although the combustion of these products was much slower than the detonation process,it aligned with the dynamic response of the structure,which enhanced the explosive yield.Excluding afterburning from the analysis reduced the center-point deformation of the structure by 30%.Following the inclusion of afterburning,the simulated quasistatic pressure increased by approximately 45%.Subsequent comparisons highlighted the merits of the proposed approach over conventional methods.This approach eliminates the reliance on empirical parameters,such as the amount and rate of energy release during afterburning,thereby laying the foundation for understanding load evolution in more complex environments,such as ships,buildings,and underground tunnels.展开更多
Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical cha...Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.展开更多
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.展开更多
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.展开更多
Coal mine belt fire develops very rapidly and is difficult to control. If not suppressed quickly, a belt fire could easily lead to airflow disorder and undermine the ventilation system. However, belt fire can be preve...Coal mine belt fire develops very rapidly and is difficult to control. If not suppressed quickly, a belt fire could easily lead to airflow disorder and undermine the ventilation system. However, belt fire can be prevented effectively by establishing fire airflow control system. In this work, the 5th belt roadway of Kongzhuang coal mine was taken as the object of investigation, where geometrical models of this roadway were established firstly. Then, based on mathematical model of fire smoke flow, the CO volume fraction, smoke density distribution, air temperature and pollutant velocity vector in the roadway before and after taking airflow control measures were simulated by using Fluent software. It can be known from the simulation that with the normal ventilation status in 5th belt roadway, the countercurrent of smoke does not happen when a fire occurs; the roadway's section is almost filled with CO at 10 m downstream from the fire source, and with air velocity getting stable gradually, the CO concentration reaches about 15 %. After taking airflow control measures, the effect range of temperature field which are harmful to the miners decreases from 69 m to 30 m; and the distance of the roadway fully filled with CO is 5 m farther than that before taking measures. Finally, according to the numerical simulation results and the actual condition of the belt roadway, the warning and automatic remote airflow control system with short-circuit method for the 5th belt roadway was designed to guarantee the safety production.展开更多
In order to research spread law and distribution law of temperature nearby fire sources on roadway in mine, according to combustion theory and other basic, the theory model of temperature attenuation was determined un...In order to research spread law and distribution law of temperature nearby fire sources on roadway in mine, according to combustion theory and other basic, the theory model of temperature attenuation was determined under unsteady heat-exchange between wind and roadway wall. The full-size roadway fire simulation experiments were carried out in Chongqing Research Institute of China Coal Technology & Engineering Group Corporation. The development processes of mine fire and flow pattern of high temperature gas were analyzed. Experimental roadway is seen as physical model, and through using CFD software, the processes of mine fire have been simulated on computer. The results show that, after fire occurs, if the wind speed is less than the minimum speed which can prevent smoke from rolling back, then the smaller wind speed can cause smoke to roll back easily. Hot plume will lead to secondary disasters in upwind side. Because of roadway wall, hot plume released from roadway fire zone has caused the occurrence of the ceiling jet, and the hot plume has been forced down. Whereas, owing to the higher temperature, buoyancy effect is more obvious. Therefore, smoke rises gradually along the roadway in the flow process, and the hierarchical interface appears wavy.Oxygen-enriched combustion and fuel-enriched combustion are the two kinds of combustion states of fire. The oxygen content of downwind side of fire is maintained at around 15% for oxygen-enriched combustion, and the oxygen content of downwind side of fire is maintained at around 2% for fuel-enriched combustion. Furthermore, fuel-enriched combustion can lead to secondary disasters easily.展开更多
Considering the fact that the temperature distribution in furnace of a tangential fired pulverized coal boiler is difficult to be measured and monitored, two-stage numerical simulation method was put forward. First, m...Considering the fact that the temperature distribution in furnace of a tangential fired pulverized coal boiler is difficult to be measured and monitored, two-stage numerical simulation method was put forward. First, multi-field coupling simulation in typical work conditions was carried out off-line with the software CFX-4.3, and then the expression of temperature profile varying with operating parameter was obtained. According to real-time operating parameters, the temperature at arbitrary point of the furnace can be calculated by using this expression. Thus the temperature profile can be shown on-line and monitoring for combustion state in the furnace is realized. The simul-(ation) model was checked by the parameters measured in an operating boiler, (DG130-9.8/540.) The maximum of relative error is less than 12% and the absolute error is less than 120℃, which shows that the proposed two-stage simulation method is reliable and able to satisfy the requirement of industrial application.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading ef...In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading effect of roof vibration on roadway surrounding rocks as well as the impact on stability. The results show that, on one hand, hard-thick roof will result in high stress concentration on mining surrounding rocks; on the other hand, the breaking of hard-thick roof will lead to mining seismicity, causing dynamic loading effect on coal and rock mass. High stress concentration and dynamic loading combination reaches to the mechanical conditions for the occurrence of rock burst, which will induce rock burst. The mining induced seismic events occurring in the roof breaking act on the mining surrounding rocks in the form of stress wave. The stress wave then has a reflection on the free surface of roadway and the tensile stress will be generated around the free surface. Horizontal vibration of roadway surrounding particles will cause instant changes of horizontal stress of roadway surrounding rocks; the horizontal displacement is directly related to the horizontal stress but is not significantly correlated with the vertical stress; the increase of horizontal stress of roadway near surface surrounding rocks and the release of elastic deformation energy of deep surrounding coal and rock mass are immanent causes that lead to the impact instability of roadway surrounding rocks. The most significant measures for rock burst prevention are controlling of horizontal stress and vibration strength.Key words展开更多
Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was ...Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was usually observed near the excavation boundaries. The formation mechanism of the “11·28” rockburst, which was a typical rockburst and occurred in a drainage tunnel under a deep burial depth, high in-situ stress state and complex geological conditions, has been difficult to explain. Realistic failure process analysis(RFPA3D) software was adopted to numerically simulate the whole failure process of the surrounding rock mass around the tunnel subjected to excavation. The spatial distribution of acoustic emission derived from numerical simulation contributed to explaining the mechanical responses of the process. Analyses of the stress, safety reserve coefficient and damage degree were performed to reveal the effect of faults on the formation of rockbursts in the deep tunnel. The existence of faults results in the formation of stress anomaly areas between the tunnel and the fault. The surrounding rock mass failure propagates toward the fault from the initial failure, to different degrees. The relative positions and angles of faults play significant roles in the extent and development of surrounding rock mass failure, respectively. The increase in the lateral stress coefficient leads to the aggravation of the surrounding rock mass damage, especially in the roof and floor of the tunnel. Moreover, as the rock strength-stress ratio increases, the failure mode of the near-fault tunnel gradually changes from the stress-controlled type to the compound-controlled type. These findings were consistent with the microseismic monitoring results and field observations, which was helpful to understand the mechanical behavior of tunnel excavation affected by faults. The achievements of this study can provide some references for analysis of the failure mechanisms of similar deep tunnels.展开更多
To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mat...To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.展开更多
In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by usi...In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by using FLUENT software,and the numerical simulation analysis of tunnel drilling ventilation and dust removal parameters was carried out.The results show that it is difficult to control the dust pollution of the face by conventional ventilation,and the drilling dust is distributed in the range of 10 m from the face;after the introduction of the long pressure and short suction ventilation scheme,when the ratio of compressed air volume to exhaust air volume is 0.72,the height of the pressure fan is 2.5 m,the distance between the pressure fan and the palm face is 20 m,and the exhaust fan is 12 m away from the palm,the dust concentration control efficiency of the working face is increased by about 60%.Therefore,in the similar long-distance single head tunnel construction,it is appropriate to adopt the dust removal method of long-distance short suction and exhaust fan to ensure the working environment.展开更多
During the compaction of a road subgrade, the mechanical parameters of the soil mass change in real time, but current research assumes that these parameters remain unchanged. In order to address this discrepancy, this...During the compaction of a road subgrade, the mechanical parameters of the soil mass change in real time, but current research assumes that these parameters remain unchanged. In order to address this discrepancy, this paper establishes a relationship between the degree of compaction K and strain ε. The relationship between the compaction degree K and the shear strength of soil(cohesion c and frictional angle φ) was clearly established through indoor experiments. The subroutine UMAT in ABAQUS finite element numerical software was developed to realize an accurate calculation of the subgrade soil compaction quality. This value was compared and analyzed against the assumed compaction value of the model, thereby verifying the accuracy of the intelligent compaction calculation results for subgrade soil. On this basis, orthogonal tests of the influential factors(frequency, amplitude, and quality) for the degree of compaction and sensitivity analysis were carried out. Finally, the ‘acceleration intelligent compaction value’, which is based on the acceleration signal, is proposed for a compaction meter value that indicates poor accuracy. The research results can provide guidance and basis for further research into the accurate control of compaction quality for roadbeds and pavements.展开更多
Frac-packing technology has been introduced to improve the development effect of weakly consolidated sandstone.It has double effects on increasing production and sand control.However,determining operation parameters o...Frac-packing technology has been introduced to improve the development effect of weakly consolidated sandstone.It has double effects on increasing production and sand control.However,determining operation parameters of frac-packing is the key factor due to the particularity of weakly consolidated sandstone.In order to study the mechanisms of hydraulic fracture propagation and reveal the effect of fracturing parameters on fracture morphology in weakly consolidated sandstone,finite element numerical model of fluid-solid coupling is established to carry out numerical simulation to analyze influences of mechanical characteristics,formation permeability,fracturing fluid injection rate and viscosity on fracture propagation.The result shows that lower elastic modulus is favorable for inducing short and wide fractures and controls the fracture length while Poisson ratio has almost no effect.Large injection rate and high viscosity of fracturing fluid are advantageous to fracture initiation and propagation.Suitable fractures are produced when the injection rate is approximate to3–4m3/min and fluid viscosity is over100mPa?s.The leak-off of fracturing fluid to formation is rising with the increase of formation permeability,which is adverse to fracture propagation.The work provides theoretical reference to determine the construction parameters for the frac-packing design in weakly consolidated reservoirs.展开更多
Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out...Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area was carried out with FLAC3D. The whole analytical process was divided into two steps, including the static analysis and the dynamic analysis which were used to simulate the influence of excavation process and blasting vibration respectively. The results show that the shape of right upper boundary is extremely irregular after excavation, and stress concentration occurs at many places and higher tensile stress appears. The maximum tensile stress is higher than the tensile strength of rock mass, and surrounding rock of right roof will be damaged with tension fracture. The maximum displacement of surrounding rock is 4.75 mm after excavation. However, the maximum displacement increases to 5.47 mm after the blasting dynamic load is applied. And the covering area of plastic zones expands obviously, especially at the foot of right upper slope. The analytical results are in basic accordance with the observed results on the whole. Damage and disturbance on surrounding rock to some degree are caused by excavation, while blasting dynamic load increases the possibility of occurrence of dynamic instability and destruction further. So the effective supporting and vibration reducing measures should be taken during mining.展开更多
Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an...Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an orthogonal design with two factors, three levels and nine runs was proposed, followed by three-dimensional numerical simulation using ANSYS and FLAC3~. Results show that surface settlement after excavation is concentrically ringed, and increases with the decrease of pillar width and distances to stope gobs. In the meantime, the ore-control fault at the ore-rock boundary and the fractured argillaceous dolomite with intercalated slate at the hanging wall deteriorate the roof settlement. Additionally, stope stability is challenged due to pillar rheological yield and stress concentration, and both are induced by redistribution of stress and plastic zones after mining. Following an objective function and a constraint function, room and pillar configuration with widths of 14 m and 16 m, respectively, is presented as the optimization for improving the ore recovery rate while maintaining a safe working environment.展开更多
Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated...Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method （FEM）. The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/（m2. K）. When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the： decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.展开更多
文摘During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.
基金supported by the National Natural Science Foundation of China(Grant Nos.52171318 and 12202329)Joint Foundation of the Ministry of Education(Grant No.8091B022105)。
文摘The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process greatly influences the evolution of confined blast loading and the subsequent structural response,which is crucial in confined blast scenarios.Given the complex nature of the reaction process,accurate analysis of the afterburning effect remains challenging.Previous studies have either overlooked the mechanisms of detonation product combustion or failed to provide experimental validation.This study introduces a three-dimensional model to effectively characterize the combustion of detonation products.The model integrates chemical reaction source terms into the governing equations to consider the combustion processes.Numerical simulations and experimental tests were conducted to analyze the combustion and energy release from the detonation products of fuel-rich explosives in confined spaces.Approximately 50%of the energy was released during the combustion of detonation products in a confined TNT explosion.Although the combustion of these products was much slower than the detonation process,it aligned with the dynamic response of the structure,which enhanced the explosive yield.Excluding afterburning from the analysis reduced the center-point deformation of the structure by 30%.Following the inclusion of afterburning,the simulated quasistatic pressure increased by approximately 45%.Subsequent comparisons highlighted the merits of the proposed approach over conventional methods.This approach eliminates the reliance on empirical parameters,such as the amount and rate of energy release during afterburning,thereby laying the foundation for understanding load evolution in more complex environments,such as ships,buildings,and underground tunnels.
文摘Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.
基金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.
基金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.
基金Project supported by Joint Innovative Center for Safe and Effective Mining Technology and Equipment of Coal Resources of Shandong Province,ChinaProject supported by Taishan Scholar Program of Shandong Province,China+2 种基金Project(2014JQJH106)supported by Shandong University of Science and Technology Research Fund,ChinaProject(LAK2012-1)supported by Science and Technology Development Program of Safety Production of Shandong Province,ChinaProject(2012ZHTD06)supported by Science Research Innovative Group of College of Mining and Safety Engineering of Shandong University of Science and Technology,China
文摘Coal mine belt fire develops very rapidly and is difficult to control. If not suppressed quickly, a belt fire could easily lead to airflow disorder and undermine the ventilation system. However, belt fire can be prevented effectively by establishing fire airflow control system. In this work, the 5th belt roadway of Kongzhuang coal mine was taken as the object of investigation, where geometrical models of this roadway were established firstly. Then, based on mathematical model of fire smoke flow, the CO volume fraction, smoke density distribution, air temperature and pollutant velocity vector in the roadway before and after taking airflow control measures were simulated by using Fluent software. It can be known from the simulation that with the normal ventilation status in 5th belt roadway, the countercurrent of smoke does not happen when a fire occurs; the roadway's section is almost filled with CO at 10 m downstream from the fire source, and with air velocity getting stable gradually, the CO concentration reaches about 15 %. After taking airflow control measures, the effect range of temperature field which are harmful to the miners decreases from 69 m to 30 m; and the distance of the roadway fully filled with CO is 5 m farther than that before taking measures. Finally, according to the numerical simulation results and the actual condition of the belt roadway, the warning and automatic remote airflow control system with short-circuit method for the 5th belt roadway was designed to guarantee the safety production.
基金Projects(51274099,51474106)supported by the National Natural Science Foundation of ChinaProject(10C0660)supported by the Scientific Research Fund of Hunan Provincial Education Department,China
文摘In order to research spread law and distribution law of temperature nearby fire sources on roadway in mine, according to combustion theory and other basic, the theory model of temperature attenuation was determined under unsteady heat-exchange between wind and roadway wall. The full-size roadway fire simulation experiments were carried out in Chongqing Research Institute of China Coal Technology & Engineering Group Corporation. The development processes of mine fire and flow pattern of high temperature gas were analyzed. Experimental roadway is seen as physical model, and through using CFD software, the processes of mine fire have been simulated on computer. The results show that, after fire occurs, if the wind speed is less than the minimum speed which can prevent smoke from rolling back, then the smaller wind speed can cause smoke to roll back easily. Hot plume will lead to secondary disasters in upwind side. Because of roadway wall, hot plume released from roadway fire zone has caused the occurrence of the ceiling jet, and the hot plume has been forced down. Whereas, owing to the higher temperature, buoyancy effect is more obvious. Therefore, smoke rises gradually along the roadway in the flow process, and the hierarchical interface appears wavy.Oxygen-enriched combustion and fuel-enriched combustion are the two kinds of combustion states of fire. The oxygen content of downwind side of fire is maintained at around 15% for oxygen-enriched combustion, and the oxygen content of downwind side of fire is maintained at around 2% for fuel-enriched combustion. Furthermore, fuel-enriched combustion can lead to secondary disasters easily.
文摘Considering the fact that the temperature distribution in furnace of a tangential fired pulverized coal boiler is difficult to be measured and monitored, two-stage numerical simulation method was put forward. First, multi-field coupling simulation in typical work conditions was carried out off-line with the software CFX-4.3, and then the expression of temperature profile varying with operating parameter was obtained. According to real-time operating parameters, the temperature at arbitrary point of the furnace can be calculated by using this expression. Thus the temperature profile can be shown on-line and monitoring for combustion state in the furnace is realized. The simul-(ation) model was checked by the parameters measured in an operating boiler, (DG130-9.8/540.) The maximum of relative error is less than 12% and the absolute error is less than 120℃, which shows that the proposed two-stage simulation method is reliable and able to satisfy the requirement of industrial application.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
基金Project(51404243)supported by the National Natural Science Foundation of ChinaProject(2014QNB26)supported by the Fundamental Research Funds for the Central Universities,China
文摘In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading effect of roof vibration on roadway surrounding rocks as well as the impact on stability. The results show that, on one hand, hard-thick roof will result in high stress concentration on mining surrounding rocks; on the other hand, the breaking of hard-thick roof will lead to mining seismicity, causing dynamic loading effect on coal and rock mass. High stress concentration and dynamic loading combination reaches to the mechanical conditions for the occurrence of rock burst, which will induce rock burst. The mining induced seismic events occurring in the roof breaking act on the mining surrounding rocks in the form of stress wave. The stress wave then has a reflection on the free surface of roadway and the tensile stress will be generated around the free surface. Horizontal vibration of roadway surrounding particles will cause instant changes of horizontal stress of roadway surrounding rocks; the horizontal displacement is directly related to the horizontal stress but is not significantly correlated with the vertical stress; the increase of horizontal stress of roadway near surface surrounding rocks and the release of elastic deformation energy of deep surrounding coal and rock mass are immanent causes that lead to the impact instability of roadway surrounding rocks. The most significant measures for rock burst prevention are controlling of horizontal stress and vibration strength.Key words
基金Project(42177143) supported by the National Natural Science Foundation of ChinaProject(2020JDJQ0011) supported by the Science Foundation for Distinguished Young Scholars of Sichuan Province,China。
文摘Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was usually observed near the excavation boundaries. The formation mechanism of the “11·28” rockburst, which was a typical rockburst and occurred in a drainage tunnel under a deep burial depth, high in-situ stress state and complex geological conditions, has been difficult to explain. Realistic failure process analysis(RFPA3D) software was adopted to numerically simulate the whole failure process of the surrounding rock mass around the tunnel subjected to excavation. The spatial distribution of acoustic emission derived from numerical simulation contributed to explaining the mechanical responses of the process. Analyses of the stress, safety reserve coefficient and damage degree were performed to reveal the effect of faults on the formation of rockbursts in the deep tunnel. The existence of faults results in the formation of stress anomaly areas between the tunnel and the fault. The surrounding rock mass failure propagates toward the fault from the initial failure, to different degrees. The relative positions and angles of faults play significant roles in the extent and development of surrounding rock mass failure, respectively. The increase in the lateral stress coefficient leads to the aggravation of the surrounding rock mass damage, especially in the roof and floor of the tunnel. Moreover, as the rock strength-stress ratio increases, the failure mode of the near-fault tunnel gradually changes from the stress-controlled type to the compound-controlled type. These findings were consistent with the microseismic monitoring results and field observations, which was helpful to understand the mechanical behavior of tunnel excavation affected by faults. The achievements of this study can provide some references for analysis of the failure mechanisms of similar deep tunnels.
基金Project(2006AA03Z523) supported by the National High-Tech Research and Development Program of ChinaProject(08C26224302178) supported by the Innovation Foundation of Central South University,China
文摘To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.
基金Project(51874016)supported by the National Natural Science Foundation of China。
文摘In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by using FLUENT software,and the numerical simulation analysis of tunnel drilling ventilation and dust removal parameters was carried out.The results show that it is difficult to control the dust pollution of the face by conventional ventilation,and the drilling dust is distributed in the range of 10 m from the face;after the introduction of the long pressure and short suction ventilation scheme,when the ratio of compressed air volume to exhaust air volume is 0.72,the height of the pressure fan is 2.5 m,the distance between the pressure fan and the palm face is 20 m,and the exhaust fan is 12 m away from the palm,the dust concentration control efficiency of the working face is increased by about 60%.Therefore,in the similar long-distance single head tunnel construction,it is appropriate to adopt the dust removal method of long-distance short suction and exhaust fan to ensure the working environment.
基金Project(51878164) supported by the National Natural Science Foundation of ChinaProjects(BK20161421, BK20140109) supported by the Natural Science Foundation of Jiangsu Province, China+4 种基金Project(141076) supported by the Huoyingdong Foundation of the Ministry of Education of ChinaProject(BZ2017011) supported by the Science and Technology Support Project of Jiangsu Province, ChinaProject(2242015R30027) supported by the Fundamental Research Funds for the Central Universities, ChinaProject(grant number KFJ170106) supported by the Changsha University of Science & Technology via Open Fund of National Engineering Laboratory of Highway Maintenance Technology, ChinaProject(2018B51) supported by the Science and Technology Support Project of Qilu Transportation Development Group, China。
文摘During the compaction of a road subgrade, the mechanical parameters of the soil mass change in real time, but current research assumes that these parameters remain unchanged. In order to address this discrepancy, this paper establishes a relationship between the degree of compaction K and strain ε. The relationship between the compaction degree K and the shear strength of soil(cohesion c and frictional angle φ) was clearly established through indoor experiments. The subroutine UMAT in ABAQUS finite element numerical software was developed to realize an accurate calculation of the subgrade soil compaction quality. This value was compared and analyzed against the assumed compaction value of the model, thereby verifying the accuracy of the intelligent compaction calculation results for subgrade soil. On this basis, orthogonal tests of the influential factors(frequency, amplitude, and quality) for the degree of compaction and sensitivity analysis were carried out. Finally, the ‘acceleration intelligent compaction value’, which is based on the acceleration signal, is proposed for a compaction meter value that indicates poor accuracy. The research results can provide guidance and basis for further research into the accurate control of compaction quality for roadbeds and pavements.
基金Project(2016ZX05058-002-006)supported by National Science and Technology Major Projects of ChinaProject(2018CXTD346)supported by Innovative Research Team Program of Natural Science Foundation of Hainan Province,China
文摘Frac-packing technology has been introduced to improve the development effect of weakly consolidated sandstone.It has double effects on increasing production and sand control.However,determining operation parameters of frac-packing is the key factor due to the particularity of weakly consolidated sandstone.In order to study the mechanisms of hydraulic fracture propagation and reveal the effect of fracturing parameters on fracture morphology in weakly consolidated sandstone,finite element numerical model of fluid-solid coupling is established to carry out numerical simulation to analyze influences of mechanical characteristics,formation permeability,fracturing fluid injection rate and viscosity on fracture propagation.The result shows that lower elastic modulus is favorable for inducing short and wide fractures and controls the fracture length while Poisson ratio has almost no effect.Large injection rate and high viscosity of fracturing fluid are advantageous to fracture initiation and propagation.Suitable fractures are produced when the injection rate is approximate to3–4m3/min and fluid viscosity is over100mPa?s.The leak-off of fracturing fluid to formation is rising with the increase of formation permeability,which is adverse to fracture propagation.The work provides theoretical reference to determine the construction parameters for the frac-packing design in weakly consolidated reservoirs.
基金Project (50490272) supported by the National Natural Science Foundation of China project(NCET-05-0687) supportedby Programfor New Century Excellent Talents project (040109) supported bythe Doctor Degree Paper Innovation Engineering of CentralSouth University
文摘Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area was carried out with FLAC3D. The whole analytical process was divided into two steps, including the static analysis and the dynamic analysis which were used to simulate the influence of excavation process and blasting vibration respectively. The results show that the shape of right upper boundary is extremely irregular after excavation, and stress concentration occurs at many places and higher tensile stress appears. The maximum tensile stress is higher than the tensile strength of rock mass, and surrounding rock of right roof will be damaged with tension fracture. The maximum displacement of surrounding rock is 4.75 mm after excavation. However, the maximum displacement increases to 5.47 mm after the blasting dynamic load is applied. And the covering area of plastic zones expands obviously, especially at the foot of right upper slope. The analytical results are in basic accordance with the observed results on the whole. Damage and disturbance on surrounding rock to some degree are caused by excavation, while blasting dynamic load increases the possibility of occurrence of dynamic instability and destruction further. So the effective supporting and vibration reducing measures should be taken during mining.
基金Projects(50934002,51074013,51104100)supported by the National Natural Science Foundation of ChinaProject(IRT0950)supported by the Program for Changjiang Scholars and Innovative Research Team in University of China
文摘Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an orthogonal design with two factors, three levels and nine runs was proposed, followed by three-dimensional numerical simulation using ANSYS and FLAC3~. Results show that surface settlement after excavation is concentrically ringed, and increases with the decrease of pillar width and distances to stope gobs. In the meantime, the ore-control fault at the ore-rock boundary and the fractured argillaceous dolomite with intercalated slate at the hanging wall deteriorate the roof settlement. Additionally, stope stability is challenged due to pillar rheological yield and stress concentration, and both are induced by redistribution of stress and plastic zones after mining. Following an objective function and a constraint function, room and pillar configuration with widths of 14 m and 16 m, respectively, is presented as the optimization for improving the ore recovery rate while maintaining a safe working environment.
基金Project(20080431380) supported by China Postdoctoral Science Foundation
文摘Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method （FEM）. The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/（m2. K）. When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the： decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.
