Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
The turbulent characteristics of the top-blown Laval nozzle and the influence of pressure and Mach number were studied through numerical simulation.With 2.72%error between the results and the empirical formula,the res...The turbulent characteristics of the top-blown Laval nozzle and the influence of pressure and Mach number were studied through numerical simulation.With 2.72%error between the results and the empirical formula,the results are reliable.Nozzle fluid is influenced by pipe structure,causing pressure and density to drop as speed increases.Differences in pressure and velocity between the jet and surrounding gas lead to jet velocity attenuation,flow expansion,deflection,and eddy currents.The optimal top blowing pressure is 0.6 MPa,and the center velocity and width of the jet are 345 m/s and 0.124 m,respectively,at 20De(De is the nozzle exit diameter).It achieves a maximum jet velocity of 456 m/s.The optimal nozzle Mach number is 1.75,with a maximum jet velocity of 451 m/s.At 20D_(e),the jet center velocity is 338 m/s,with a width of 0.12 m.展开更多
In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components ...In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components in aircraft.The interfacial plastic flow behavior and bonding mechanism of this process were investigated by a developed coupling EulerianLagrangian numerical model using software ABAQUS and a novel thermo-physical simulation method with designed embedded hot compression specimen.In addition,the formation mechanism and control method of welding defects caused by uneven plastic flow were discussed.The results reveal that the plastic flow along oscillating direction of this process is even and sufficient.In the direction perpendicular to oscillation,thermo-plastic metals mainly flow downward along welding interface under coupling of shear stress and interfacial pressure,resulting in the interfacial plastic zone shown as an inverted“V”shape.The upward plastic flow in this direction is relatively weak,and only a small amount of flash is extruded from top of joint.Moreover,the wedge block and welding components at top of joint are always in un-steady friction stage,leading to nonuniform temperature field distribution and un-welded defects.According to the results of numerical simulation,high oscillating frequency combined with low pressure and small amplitude is considered as appropriate parameter selection scheme to improve the upward interfacial plastic flow at top of joint and suppress the un-welded defects.The results of thermo-physical simulation illustrate that continuous dynamic recrystallization(CDRX)induces the bonding of interface,accompanying by intense dislocation movement and creation of many low-angle grain boundaries.In the interfacial bonding area,grain orientation is random with relatively low texture density(5.0 mud)owing to CDRX.展开更多
Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphas...Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphaseflow within the furnace.Understanding the flow structure and temperature distribution in this setup is crucial foroptimizing production.In this study,gas-liquid interactions,and temperature profiles under varying air-injectionconditions are examined by means of numerical simulation for a 3.2 m×20 m furnace.The results indicate that thehigh-velocity regions are essentially distributed near the lance within the reaction region and the flue gas outlet,while low-velocity regions are located close to the furnace walls on both side of the reaction region.Dead regionsappear in the sedimentation region,with gas velocities surpassing those of the molten phase.As the injection rateincreases from 0.50 to 0.80 Nm3/s,the stabilization time of the average liquid surface velocity decreases from 2.6 sto 1.9 s,exhibiting a similar trend to the gas holdup.During stabilization,the average liquid surface velocity risesfrom 0.505 to 0.702 m/s.The average turbulent kinetic energy(TKE)of the fluid in the molten bath increases from0.095 to 0.162 m^(2)/s^(2).The proportion of the area distribution with TKE greater than 0.10 m^(2)/s^(2) and the gas holdupat steady state both rise with an increase in the injection quantity.The maximum splashing height of the melt growsfrom approximately 0.756 to 1.154 m,with the affected area expanding from 14.239 to 20.498 m^(2).Under differentworking conditions with varying injection quantities,the average temperature changes in melt zone and flue gaszone of the furnace are small.The temperature in the melt and in the flue-gas zone spans the interval 1200℃–1257℃,and 1073℃–1121℃,respectively.The temperature distribution of the melt and flue gas reveals a patterncharacterized by elevated temperatures in the reaction zone,gradually transitioning to lower temperatures in thesedimentation region.展开更多
To mitigate the sand burial of highways in sandy regions,a separated subgrade design was widely adopted in the embankments of high-grade highways,but the problem of sand deposition on subgrade slopes and pavements sti...To mitigate the sand burial of highways in sandy regions,a separated subgrade design was widely adopted in the embankments of high-grade highways,but the problem of sand deposition on subgrade slopes and pavements still happens frequently.Based on the theory of wind-sand two-phase flow,this paper constructed a three-dimensional model of the separated subgrade,the wind-sand flow transport law around the subgrade with varying median strip widths and concave depths was simulated by Fluent software.After comparison and analysis of seven subgrade models,the flow field distribution,wind speed horizontal variation,and erosion-deposition characteristics were investigated.The findings are as follows:(1)The width of the median strip in the separated subgrade had significant influences on the wind-sand flow.The smooth passage of wind-sand flow over the road surface was facilitated with the increase of the median strip width.However,sand deposition in the median strip happened.It can lead to secondary sand damage of downwind subgrade and increase the work load of road sand removal for subsequent maintenance.(2)The obstruction to airflow and sand accumulation was aggravated with greater concave depth of the median strip.Therefore,it is advisable to minimize the concave depth of the median strip in case of more sand damage.(3)A median strip width exceeding 12 m(possibly without guardrails)for an integral embankment without enough road land is recommended.Conversely,median strip width of over 40 m for separate subgrade with unrestricted land is suggested.(4)In the case of sand deposition in the existing separated subgrade,the median strip can be filled by sand deposition or other materials,then was covered with gravel to form a flat ground like Gobi smooth surface,which can let the wind-blown sand flow pass through the subgrade section without sand deposition.展开更多
A new three-dimensional multiphase numerical model was built.The volume of fluid and k-ε turbulence models were used to investigate the hot metal ladle pouring process.During the pouring process,issues such as iron s...A new three-dimensional multiphase numerical model was built.The volume of fluid and k-ε turbulence models were used to investigate the hot metal ladle pouring process.During the pouring process,issues such as iron splashing,overflow,and significant heat loss are prevalent.To realize efficient and stable pouring,the effects of ladle tilting velocity,flow rate,and converter tilting angle on the pouring process were examined.The model was verified by comparing the actual pouring time with the numerical results.It is shown that there is a nonlinear relationship between pouring velocity and hot metal flow rate at the ladle mouth.As the mass flow increased and the converter tilting angle decreased,the impact point of the hot metal into the converter pool shifted from the side wall to the bottom,and the impact force increased accordingly.The pouring velocity curve was optimized by the volume difference of the ladle at different angles,and an empirical formula was derived.After the optimization of pouring speed,the flow rate was stabilized between 4000 and 6000 kg/s,and the pouring time was reduced by approximately 30 s.After applying this model in actual production,the hot metal temperature inside the converter increased by approximately 5℃statistically.This model is potential to enhance the production efficiency,stability,and safety of the pouring process between open containers.展开更多
A CFD based numerical simulation of flow velocity of hydrocyclone was conducted with different structural and operational parameters to investigate its distribution characteristics and influencing mechanism. The resul...A CFD based numerical simulation of flow velocity of hydrocyclone was conducted with different structural and operational parameters to investigate its distribution characteristics and influencing mechanism. The results show there exist several unsymmetrical envelopes of equal vertical velocities in both upward inner flows and downward outer flows in the hydrocyclone, and the cone angle and apex diameter have remarkable influence on the vertical location of the cone bottom of the envelope of zero vertical velocity. It is also found that the tangential velocity isolines exist in the horizontal planes located in the effective separation region of hydrocyclone. The increase of feed pressure has almost no effect on the distribution characteristics of both vertical velocity and tangential velocity in hydrocyclone, but the magnitude and gradient of tangential velocity are increased obviously to make the motion velocity of high density particles to the wall increased and to make the cyclonic separation effect improved.展开更多
The friction and wear tests were performed using Nitinol 60 alloy pin sliding over GCr15 steel disk in a pin-on-disk tribometer system under PAO oil lubrication conditions. It was found that Nitinol 60 alloy can be lu...The friction and wear tests were performed using Nitinol 60 alloy pin sliding over GCr15 steel disk in a pin-on-disk tribometer system under PAO oil lubrication conditions. It was found that Nitinol 60 alloy can be lubricated well and has shown remarkable tribological performance. Average coefficient of friction （COF） of Nitinol 60 is 0.6 under dry friction; however, average COF decreases to 0.1 under PAO oil lubrication. SEM image of the worn surface shows that Nitinol 60 exhibits excellent wear resistance and the wear mechanism is mainly adhesive wear. Flow pattern of oil-air flow in oil pipe was simulated by FLUENT software with VOF model for acquiring working performance of oil-air lubrication. The optimum velocity of oil and air at the inlet was achieved, which provides the great proposal for the design of experiment of oil-air lubrication of Nitinol 60 alloy. The simulation results showed that the optimum annular flow of flow pattern was obtained when air velocity is 10 m/s and oil velocity is 0.05 m/s. The formation mechanism of annular flow was also discussed in the present study.展开更多
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.展开更多
Investigating the wind-sand flow response regularity in the longitudinal slope sections of desert highways provides a scientific basis for selecting the slope of desert roads.This study uses the Tengger Desert section...Investigating the wind-sand flow response regularity in the longitudinal slope sections of desert highways provides a scientific basis for selecting the slope of desert roads.This study uses the Tengger Desert section of the Wuhai-Maqin Expressway as a case study,employing CFD numerical simulation methods to calculate and analyze the wind-sand flow field distribution characteristics in different longitudinal slope sections.The results show that:(1)Along with the direction of the incoming flow,the windward and leeward slope toes of the embankment are low-wind-speed zones,with the wind speed at the leeward slope toe being even lower.The higher the embankment,the larger the low-wind-speed zone at the windward and leeward slope toes.As the longitudinal slope increases,the extent of the lowwind-speed zone at the same location along the route also increases.(2)Along the route direction,the wind speed at the windward and leeward slope toes decreases as embankment height increases.At the embankment toe,sand particles are transported from the top to the bottom of the longitudinal slope,and the greater the longitudinal slope,the stronger the transport effect.(3)Along the route direction,the sand accumulation around the embankment gradually gathers toward the bottom of the longitudinal slope as the slope increases.When the longitudinal slope is 3%and 4%,the trend of sand accumulation moving from the windward side at the end of the route to the leeward side at the start of the route is more significant.When the longitudinal slope is less than or equal to 3%,severe sand accumulation within the embankment range is reduced by 86.4%or more compared to when the slope is 4%.(4)Under the same longitudinal slope,the higher the embankment height,the smaller its transport rate.When the embankment height is the same,the greater the longitudinal slope,the greater the embankment transport rate.展开更多
Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On ...Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On the other hand,the entry of ice crystal particles into the combustion chamber can cause a decrease in temperature or even flameout,leading to engine surge or shutdown.Therefore,it is necessary to conduct multiphase flow tests on ice crystals for aircraft components such as aircraft engines.Conducting ice crystal multiphase flow tests on aircraft is an effective research method,but it requires the construction of an ice crystal multiphase flow test platform that meets relevant technical requirements.The paper focuses on the relevant experimental requirements and combines wind tunnel test structures to conduct multiphase flow numerical simulations on various forms of jet pipelines,obtaining particle motion distribution results.After comparison,the optimal form of jet structure is obtained,providing the best selection scheme for the design of relevant wind tunnel structures.展开更多
A CFD-based numerical model was employed to quantitatively analyze the flow characteristics of double-side-blown gas−liquid flow.Key parameters were extracted,and Spearman correlation analysis was used to quantify the...A CFD-based numerical model was employed to quantitatively analyze the flow characteristics of double-side-blown gas−liquid flow.Key parameters were extracted,and Spearman correlation analysis was used to quantify the relationships among bubble behavior,circulating flow,and liquid oscillations.The results show that periodic bubble behavior under steady injection drives the circulating flow of the liquid on both sides.The asynchronism of bubble behavior on both sides results in the alternation of circulating intensity,which significantly enhances gas−liquid mixing efficiency at certain liquid levels of 200 and 220 mm.Flow patterns of the double-side-blown process are classified into weak circulation,strong−weak alternating circulation,and strong circulation modes based on the influence of circulating flows on the penetration depth.The penetration depth in the strong−weak alternating circulation mode is generally greater than that in the single-side-blown process.The imbalance of circulating intensities on both sides primarily leads to the stable fluctuation in the injecting direction,which reveals the appearance of periodic oscillations in the molten bath.The effect of control parameters such as liquid level and gas flow rate on the liquid oscillations were discussed.展开更多
The physical and mathematical model of an operating electrowinning cell was established, and the flow of electrolyte was numerically simulated by the commercial software Fluent. The results indicate that there are two...The physical and mathematical model of an operating electrowinning cell was established, and the flow of electrolyte was numerically simulated by the commercial software Fluent. The results indicate that there are two circulations at the surface flow where part of electrolyte backflows to the inlet from the side of cell, and the rest flows directly to the outlet, and the separation of two circulations with opposite direction occurs at the 20th pair of anode-cathode. This phenomenon was observed in the real operation. The electrolyte flows into the space between anode and cathode from the side portion of the cell. Meanwhile, the interelectrode effective flow rate （IEFR） is put forward to describe quantitively the flow field characteristics and is defined as the ratio of electrolyte flow between the anode and cathode to the total flow area. The influences of structure parameters and operating conditions on IEFR, such as the inlet angle, the volumetric flow rate, the inlet position and the height of steel baffles were simulated. The inlet position has a significant influence on the IEFR and its optimal value is 0.9 m below free surface. The inlet angle should be in the range from -10° to 10°. IEFR is in linear proportion with the volumetric flow rate, and the height of the steel baffle has little influence on the flow field.展开更多
Based on K ε two equation turbulence model,the air distribution and contamination field under different conditions was numerically simulated in a vector flow cleanroom.Special mesh system was introduced to deal wit...Based on K ε two equation turbulence model,the air distribution and contamination field under different conditions was numerically simulated in a vector flow cleanroom.Special mesh system was introduced to deal with the quarter circle shaped inlets.Model experiments were also made.By analysis of numerical as well as experimental results,we made some predictions about flow characteristics,contaminant control effect and ventilation performance of this energy saving cleanroom.It has been proved that the vector flow can meet class 100 of Fed. St.with 13 to 12 conventional air change per hour,and reduce layer height and simplify air clean system.展开更多
Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of th...Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.展开更多
This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-laye...This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.展开更多
Dynamic modeling and numerical simulation of hydrate slurry flow behavior are of great importance to offshore hydrate management.For this purpose, a dynamic model of hydrate agglomeration was proposed in this paper.Ba...Dynamic modeling and numerical simulation of hydrate slurry flow behavior are of great importance to offshore hydrate management.For this purpose, a dynamic model of hydrate agglomeration was proposed in this paper.Based on population balance equation, the frame of the dynamic model was established first, which took both hydrate agglomeration and hydrate breakage into consideration.Then, the calculating methods of four key parameters involved in the dynamic model were given according to hydrate agglomeration dynamics.The four key parameters are collision frequency, agglomeration efficiency, breakage frequency and the size distribution of sub particles resulting from particle breakage.After the whole dynamic model was built, it was combined with several traditional solid–liquid flow models and then together solved by the CFD software FLUENT 14.5.Finally, using this method, the influences of flow rate and hydrate volume fraction on hydrate particle size distribution, hydrate volume concentration distribution and pipeline pressure drop were simulated and analyzed.展开更多
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
基金supported by the National Key Research and Development Project of China(No.2022YFC3902001)the National Natural Science Foundation of China(No.52004340)the Guangxi Innovation-driven Development Project,China(Nos.Gui 2021AA12006 and 2021AB26024)。
文摘The turbulent characteristics of the top-blown Laval nozzle and the influence of pressure and Mach number were studied through numerical simulation.With 2.72%error between the results and the empirical formula,the results are reliable.Nozzle fluid is influenced by pipe structure,causing pressure and density to drop as speed increases.Differences in pressure and velocity between the jet and surrounding gas lead to jet velocity attenuation,flow expansion,deflection,and eddy currents.The optimal top blowing pressure is 0.6 MPa,and the center velocity and width of the jet are 345 m/s and 0.124 m,respectively,at 20De(De is the nozzle exit diameter).It achieves a maximum jet velocity of 456 m/s.The optimal nozzle Mach number is 1.75,with a maximum jet velocity of 451 m/s.At 20D_(e),the jet center velocity is 338 m/s,with a width of 0.12 m.
基金co-supported by the National Natural Science Foundation of China(Nos.52105411,52105400and 52305420)the China Postdoctoral Science Foundation(No.2023M742830)Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University(No.CX2023008).
文摘In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components in aircraft.The interfacial plastic flow behavior and bonding mechanism of this process were investigated by a developed coupling EulerianLagrangian numerical model using software ABAQUS and a novel thermo-physical simulation method with designed embedded hot compression specimen.In addition,the formation mechanism and control method of welding defects caused by uneven plastic flow were discussed.The results reveal that the plastic flow along oscillating direction of this process is even and sufficient.In the direction perpendicular to oscillation,thermo-plastic metals mainly flow downward along welding interface under coupling of shear stress and interfacial pressure,resulting in the interfacial plastic zone shown as an inverted“V”shape.The upward plastic flow in this direction is relatively weak,and only a small amount of flash is extruded from top of joint.Moreover,the wedge block and welding components at top of joint are always in un-steady friction stage,leading to nonuniform temperature field distribution and un-welded defects.According to the results of numerical simulation,high oscillating frequency combined with low pressure and small amplitude is considered as appropriate parameter selection scheme to improve the upward interfacial plastic flow at top of joint and suppress the un-welded defects.The results of thermo-physical simulation illustrate that continuous dynamic recrystallization(CDRX)induces the bonding of interface,accompanying by intense dislocation movement and creation of many low-angle grain boundaries.In the interfacial bonding area,grain orientation is random with relatively low texture density(5.0 mud)owing to CDRX.
基金Supported by Yunnan Fundamental Research Projects(Nos.202301AT070469,202301AT070275)Supported by Yunnan Major Scientific and Technological Projects(No.202202AG050002).
文摘Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphaseflow within the furnace.Understanding the flow structure and temperature distribution in this setup is crucial foroptimizing production.In this study,gas-liquid interactions,and temperature profiles under varying air-injectionconditions are examined by means of numerical simulation for a 3.2 m×20 m furnace.The results indicate that thehigh-velocity regions are essentially distributed near the lance within the reaction region and the flue gas outlet,while low-velocity regions are located close to the furnace walls on both side of the reaction region.Dead regionsappear in the sedimentation region,with gas velocities surpassing those of the molten phase.As the injection rateincreases from 0.50 to 0.80 Nm3/s,the stabilization time of the average liquid surface velocity decreases from 2.6 sto 1.9 s,exhibiting a similar trend to the gas holdup.During stabilization,the average liquid surface velocity risesfrom 0.505 to 0.702 m/s.The average turbulent kinetic energy(TKE)of the fluid in the molten bath increases from0.095 to 0.162 m^(2)/s^(2).The proportion of the area distribution with TKE greater than 0.10 m^(2)/s^(2) and the gas holdupat steady state both rise with an increase in the injection quantity.The maximum splashing height of the melt growsfrom approximately 0.756 to 1.154 m,with the affected area expanding from 14.239 to 20.498 m^(2).Under differentworking conditions with varying injection quantities,the average temperature changes in melt zone and flue gaszone of the furnace are small.The temperature in the melt and in the flue-gas zone spans the interval 1200℃–1257℃,and 1073℃–1121℃,respectively.The temperature distribution of the melt and flue gas reveals a patterncharacterized by elevated temperatures in the reaction zone,gradually transitioning to lower temperatures in thesedimentation region.
基金supported by the Third Xinjiang Scientific Expedition and Research Program-Investigation and Risk Assessment of Drought and Aeolian Disasters in Tarim River Basin(No.2021xjkk0300)the National Natural Science Foundation of China(No.62466056)the subject of'the technical scheme and application demonstration of sand disaster prevention and control of Xinjiang expressway to engineering practice,Xinjiang Transportation Investment(Group)Co.,Ltd.(No.XJJTZKX-FWCG-202401-0043).
文摘To mitigate the sand burial of highways in sandy regions,a separated subgrade design was widely adopted in the embankments of high-grade highways,but the problem of sand deposition on subgrade slopes and pavements still happens frequently.Based on the theory of wind-sand two-phase flow,this paper constructed a three-dimensional model of the separated subgrade,the wind-sand flow transport law around the subgrade with varying median strip widths and concave depths was simulated by Fluent software.After comparison and analysis of seven subgrade models,the flow field distribution,wind speed horizontal variation,and erosion-deposition characteristics were investigated.The findings are as follows:(1)The width of the median strip in the separated subgrade had significant influences on the wind-sand flow.The smooth passage of wind-sand flow over the road surface was facilitated with the increase of the median strip width.However,sand deposition in the median strip happened.It can lead to secondary sand damage of downwind subgrade and increase the work load of road sand removal for subsequent maintenance.(2)The obstruction to airflow and sand accumulation was aggravated with greater concave depth of the median strip.Therefore,it is advisable to minimize the concave depth of the median strip in case of more sand damage.(3)A median strip width exceeding 12 m(possibly without guardrails)for an integral embankment without enough road land is recommended.Conversely,median strip width of over 40 m for separate subgrade with unrestricted land is suggested.(4)In the case of sand deposition in the existing separated subgrade,the median strip can be filled by sand deposition or other materials,then was covered with gravel to form a flat ground like Gobi smooth surface,which can let the wind-blown sand flow pass through the subgrade section without sand deposition.
基金supported by the National Natural Science Foundation of China(Grant Nos.52104318 and 52374322).
文摘A new three-dimensional multiphase numerical model was built.The volume of fluid and k-ε turbulence models were used to investigate the hot metal ladle pouring process.During the pouring process,issues such as iron splashing,overflow,and significant heat loss are prevalent.To realize efficient and stable pouring,the effects of ladle tilting velocity,flow rate,and converter tilting angle on the pouring process were examined.The model was verified by comparing the actual pouring time with the numerical results.It is shown that there is a nonlinear relationship between pouring velocity and hot metal flow rate at the ladle mouth.As the mass flow increased and the converter tilting angle decreased,the impact point of the hot metal into the converter pool shifted from the side wall to the bottom,and the impact force increased accordingly.The pouring velocity curve was optimized by the volume difference of the ladle at different angles,and an empirical formula was derived.After the optimization of pouring speed,the flow rate was stabilized between 4000 and 6000 kg/s,and the pouring time was reduced by approximately 30 s.After applying this model in actual production,the hot metal temperature inside the converter increased by approximately 5℃statistically.This model is potential to enhance the production efficiency,stability,and safety of the pouring process between open containers.
基金Project (50974033) supported by the National Natural Science Foundation of ChinaProject (N100301002) supported by the Fundamental Research Funds for the Universities, China
文摘A CFD based numerical simulation of flow velocity of hydrocyclone was conducted with different structural and operational parameters to investigate its distribution characteristics and influencing mechanism. The results show there exist several unsymmetrical envelopes of equal vertical velocities in both upward inner flows and downward outer flows in the hydrocyclone, and the cone angle and apex diameter have remarkable influence on the vertical location of the cone bottom of the envelope of zero vertical velocity. It is also found that the tangential velocity isolines exist in the horizontal planes located in the effective separation region of hydrocyclone. The increase of feed pressure has almost no effect on the distribution characteristics of both vertical velocity and tangential velocity in hydrocyclone, but the magnitude and gradient of tangential velocity are increased obviously to make the motion velocity of high density particles to the wall increased and to make the cyclonic separation effect improved.
基金Project (2012M511993) supported by China Postdoctoral Science FoundationProject (TPL1202) supported by the Open Fund Program of the State Key Laboratory of Traction Power, Southwest Jiaotong University, China
文摘The friction and wear tests were performed using Nitinol 60 alloy pin sliding over GCr15 steel disk in a pin-on-disk tribometer system under PAO oil lubrication conditions. It was found that Nitinol 60 alloy can be lubricated well and has shown remarkable tribological performance. Average coefficient of friction （COF） of Nitinol 60 is 0.6 under dry friction; however, average COF decreases to 0.1 under PAO oil lubrication. SEM image of the worn surface shows that Nitinol 60 exhibits excellent wear resistance and the wear mechanism is mainly adhesive wear. Flow pattern of oil-air flow in oil pipe was simulated by FLUENT software with VOF model for acquiring working performance of oil-air lubrication. The optimum velocity of oil and air at the inlet was achieved, which provides the great proposal for the design of experiment of oil-air lubrication of Nitinol 60 alloy. The simulation results showed that the optimum annular flow of flow pattern was obtained when air velocity is 10 m/s and oil velocity is 0.05 m/s. The formation mechanism of annular flow was also discussed in the present study.
基金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.
基金financially supported by Youth Science“Research on Failure Mechanism and Evaluation Method of Sand Control Measures for Railway Machinery in Sandy Area”(12302511)Ningxia Transportation Department Science and Technology Project(20200173)The Central Guidance on Local Science and Technology Development Funds(22ZY1QA005)。
文摘Investigating the wind-sand flow response regularity in the longitudinal slope sections of desert highways provides a scientific basis for selecting the slope of desert roads.This study uses the Tengger Desert section of the Wuhai-Maqin Expressway as a case study,employing CFD numerical simulation methods to calculate and analyze the wind-sand flow field distribution characteristics in different longitudinal slope sections.The results show that:(1)Along with the direction of the incoming flow,the windward and leeward slope toes of the embankment are low-wind-speed zones,with the wind speed at the leeward slope toe being even lower.The higher the embankment,the larger the low-wind-speed zone at the windward and leeward slope toes.As the longitudinal slope increases,the extent of the lowwind-speed zone at the same location along the route also increases.(2)Along the route direction,the wind speed at the windward and leeward slope toes decreases as embankment height increases.At the embankment toe,sand particles are transported from the top to the bottom of the longitudinal slope,and the greater the longitudinal slope,the stronger the transport effect.(3)Along the route direction,the sand accumulation around the embankment gradually gathers toward the bottom of the longitudinal slope as the slope increases.When the longitudinal slope is 3%and 4%,the trend of sand accumulation moving from the windward side at the end of the route to the leeward side at the start of the route is more significant.When the longitudinal slope is less than or equal to 3%,severe sand accumulation within the embankment range is reduced by 86.4%or more compared to when the slope is 4%.(4)Under the same longitudinal slope,the higher the embankment height,the smaller its transport rate.When the embankment height is the same,the greater the longitudinal slope,the greater the embankment transport rate.
文摘Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On the other hand,the entry of ice crystal particles into the combustion chamber can cause a decrease in temperature or even flameout,leading to engine surge or shutdown.Therefore,it is necessary to conduct multiphase flow tests on ice crystals for aircraft components such as aircraft engines.Conducting ice crystal multiphase flow tests on aircraft is an effective research method,but it requires the construction of an ice crystal multiphase flow test platform that meets relevant technical requirements.The paper focuses on the relevant experimental requirements and combines wind tunnel test structures to conduct multiphase flow numerical simulations on various forms of jet pipelines,obtaining particle motion distribution results.After comparison,the optimal form of jet structure is obtained,providing the best selection scheme for the design of relevant wind tunnel structures.
基金financial support for this research work from the National Key Research and Development Program of China(No.2022YFB3304901)。
文摘A CFD-based numerical model was employed to quantitatively analyze the flow characteristics of double-side-blown gas−liquid flow.Key parameters were extracted,and Spearman correlation analysis was used to quantify the relationships among bubble behavior,circulating flow,and liquid oscillations.The results show that periodic bubble behavior under steady injection drives the circulating flow of the liquid on both sides.The asynchronism of bubble behavior on both sides results in the alternation of circulating intensity,which significantly enhances gas−liquid mixing efficiency at certain liquid levels of 200 and 220 mm.Flow patterns of the double-side-blown process are classified into weak circulation,strong−weak alternating circulation,and strong circulation modes based on the influence of circulating flows on the penetration depth.The penetration depth in the strong−weak alternating circulation mode is generally greater than that in the single-side-blown process.The imbalance of circulating intensities on both sides primarily leads to the stable fluctuation in the injecting direction,which reveals the appearance of periodic oscillations in the molten bath.The effect of control parameters such as liquid level and gas flow rate on the liquid oscillations were discussed.
基金Project(2011AA061003)supported by the National High-Tech Research and Development Program of China
文摘The physical and mathematical model of an operating electrowinning cell was established, and the flow of electrolyte was numerically simulated by the commercial software Fluent. The results indicate that there are two circulations at the surface flow where part of electrolyte backflows to the inlet from the side of cell, and the rest flows directly to the outlet, and the separation of two circulations with opposite direction occurs at the 20th pair of anode-cathode. This phenomenon was observed in the real operation. The electrolyte flows into the space between anode and cathode from the side portion of the cell. Meanwhile, the interelectrode effective flow rate （IEFR） is put forward to describe quantitively the flow field characteristics and is defined as the ratio of electrolyte flow between the anode and cathode to the total flow area. The influences of structure parameters and operating conditions on IEFR, such as the inlet angle, the volumetric flow rate, the inlet position and the height of steel baffles were simulated. The inlet position has a significant influence on the IEFR and its optimal value is 0.9 m below free surface. The inlet angle should be in the range from -10° to 10°. IEFR is in linear proportion with the volumetric flow rate, and the height of the steel baffle has little influence on the flow field.
文摘Based on K ε two equation turbulence model,the air distribution and contamination field under different conditions was numerically simulated in a vector flow cleanroom.Special mesh system was introduced to deal with the quarter circle shaped inlets.Model experiments were also made.By analysis of numerical as well as experimental results,we made some predictions about flow characteristics,contaminant control effect and ventilation performance of this energy saving cleanroom.It has been proved that the vector flow can meet class 100 of Fed. St.with 13 to 12 conventional air change per hour,and reduce layer height and simplify air clean system.
基金China Postdoctoral Science Foundation (20100481368)National Key Laboratory Foundation of China
文摘Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.
基金supported by the National Basic Research Program of China(2009CB724100)the National Natural Science Foundation of China(10632050,10872205,11072248).
文摘This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.
基金Supported by Shandong Provincial Natural Science Foundation,China(ZR2017MEE057)the Fundamental Research Funds for the Central Universities(14CX02207A,17CX05006,17CX06017)the Graduate Innovation Project of China University of Petroleum(East China)(YCX2017062)
文摘Dynamic modeling and numerical simulation of hydrate slurry flow behavior are of great importance to offshore hydrate management.For this purpose, a dynamic model of hydrate agglomeration was proposed in this paper.Based on population balance equation, the frame of the dynamic model was established first, which took both hydrate agglomeration and hydrate breakage into consideration.Then, the calculating methods of four key parameters involved in the dynamic model were given according to hydrate agglomeration dynamics.The four key parameters are collision frequency, agglomeration efficiency, breakage frequency and the size distribution of sub particles resulting from particle breakage.After the whole dynamic model was built, it was combined with several traditional solid–liquid flow models and then together solved by the CFD software FLUENT 14.5.Finally, using this method, the influences of flow rate and hydrate volume fraction on hydrate particle size distribution, hydrate volume concentration distribution and pipeline pressure drop were simulated and analyzed.
