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.展开更多
Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical s...Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical simulation study on the dynamics of blood flow in a stenosed artery,focusing on the effects of copper and alumina nanoparticles,is conducted.The study employs a 2-dimensional Newtonian blood flow model infused with copper and alumina nanoparticles,considering the influence of a magnetic field,thermal radiation,and various flow parameters.The governing differential equations are first non-dimensionalized to facilitate analysis and subsequently solved using the 4th order collocation method,bvp4c module in MATLAB.This approach obtains velocity and temperature profiles,revealing the impact of relevant parameters crucial in the biomedical field.The findings of this study underscore the significance of understanding blood flow dynamics in stenosed arteries and the potential benefits of utilizing copper and alumina nanoparticles in treatment strategies.The incorporation of nanoparticles introduces novel avenues for enhancing therapeutic interventions,particularly in mitigating the effects of stenosis.The elucidation of velocity and temperature profiles provides valuable insights into the behavior of blood flow under different conditions,thereby informing the development of targeted biomedical applications.The arterial curvature flow parameter influences temperature profiles,with increased parameters promoting more efficient heat dissipation.The elevated values of Prandtl number and thermal radiation parameter showcase the diminished temperature profiles,indicating stronger dominance of momentum diffusion over thermal diffusion and radiative heat transfer mechanism.Sensitivity analysis of the pertinent physical parameters reveals that the Prandtl number has the most significant impact on blood flow dynamics.A statistical analysis of the present results and existing literature has also been included in the study.Overall,this research contributes to advancing our understanding of vascular health and lays the groundwork for innovative approaches in stenosis treatment and related biomedical fields.展开更多
A mild gasification process has been implemented to provide an alternative form of clean coal technology called the Integrated Mild Gasification Combined Cycle (IMGCC), which can be utilized to build a new, highly eff...A mild gasification process has been implemented to provide an alternative form of clean coal technology called the Integrated Mild Gasification Combined Cycle (IMGCC), which can be utilized to build a new, highly efficient, and compact power plant or to retrofit an existing coal-fired power plant in order to achieve lower emissions and significantly improved thermal efficiency. The core technology of the mild gasification power plant lies on the design of a compact and effective mild gasifier that can produce synthesis gases with high energy volatiles through a hybrid system: utilizing the features of both entrained-flow and fluidized bed gasifiers. To aid in the design of the mild gasifier, a computational model has been implemented to investigate the thermal-flow and gasification process inside this mild gasifier using the commercial CFD (Computational Fluid Dynamics) solver ANSYS/FLUENT. The Eulerian-Eulerian method is employed to model both the primary phase (air) and the secondary phase (coal particles). However, the Eulerian-Eulerian model used in the software does not facilitate any built-in devolatilization model. The objective of this study is therefore to implement a devolatilization model (along with demoisturization) and incorporate it into the existing code. The Navier-Stokes equations and seven species transport equations are solved with three heterogeneous (gas-solid) and two homogeneous (gas-gas) global gasification reactions. Implementation of the complete model starts from adding demoisturization first, then devolatilization, and then adding one chemical equation at a time until finally all reactions are included in the multiphase flow. The result shows that the demoisturization and devolatilization models are successfully incorporated and a large amount of volatiles are preserved as high-energy fuels in the syngas stream without being further cracked or reacted into lighter gases. The overall results are encouraging but require future experimental data for verification.展开更多
Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy ...Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.展开更多
Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the eff...Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the effectiveness of the debris flow control project and evaluated the debris flow hazards.Through field investigation and numerical simulation methods,the indicators of flow intensity reduction rate and storage capacity fullness were proposed to quantify the effectiveness of the engineering measures in the debris flow event.The simulation results show that the debris flow control project reduced the flow intensity by41.05%to 64.61%.The storage capacity of the dam decreases gradually from upstream to the mouth of the gully,thus effectively intercepting and controlling the debris flow.By evaluating the debris flow of different recurrence intervals,further measures are recommended for managing debris flow events.展开更多
A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten st...A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.展开更多
A coupled Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)approach is used to calculate the interaction of a flexible rag transported by a fluid current with a fixed solid cylinder.More specifically a hyb...A coupled Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)approach is used to calculate the interaction of a flexible rag transported by a fluid current with a fixed solid cylinder.More specifically a hybrid Eulerian-Lagrangian approach is used with the rag being modeled as a set of interconnected particles.The influence of various parameters is considered,namely the inlet velocity(1.5,2.0,and 2.5 m/s,respectively),the angle formed by the initially straight rag with the flow direction(45°,60°and 90°,respectively),and the inlet position(90,100,and 110 mm,respectively).The results show that the flow rate has a significant impact on the permeability of the rag.The higher the flow rate,the higher the permeability and the rag speed difference.The angle has a minor effect on rag permeability,with 45°being the most favorable angle for permeability.The inlet position has a small impact on rag permeability,while reducing the initial distance between the rag an the cylinder makes it easier for rags to pass through.展开更多
Casting speed,casting temperature and secondary cooling water flow rate are the main process parameters affecting the DC casting process.These parameters significantly influence the flow and temperature fields during ...Casting speed,casting temperature and secondary cooling water flow rate are the main process parameters affecting the DC casting process.These parameters significantly influence the flow and temperature fields during casting,which are crucial for the quality of the ingot and can determine the success or failure of the casting operation.Numerical simulation,with the advantages of low cost,rapid execution,and visualized results,is an important method to study and optimize the DC casting process.In the present work,a simulation model of DC casting 2024 aluminum alloy was established,and the reliability of the model was verified.Then,the influence of casting parameters on flow field and temperature field was studied in detail by numerical simulation method.Results show that with the increase of casting speed,the melt flow becomes faster,the depths of slurry zone and mushy zone increase,and the variation of slurry zone depth is greater than that of mushy zone.With an increase in casting temperature,the melt flow rate increases,the depth of the slurry zone becomes shallower,and the depth of the mushy zone experiences only minor changes.The simulation results further indicate that the increase of the flow rate of the secondary cooling water slightly reduces the depths of both slurry and mushy zone.展开更多
Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flo...Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flow field around roadbeds, overlooking detailed examinations of sand particle erosion and deposition patterns near bridges and roadbeds. This study employs numerical simulations to analyze the influence of varying heights and wind speeds on sand deposition and erosion characteristics at different locations: the bridge-road transition section(side piers), middle piers, and roadbeds. The results show that the side piers, experience greater accumulation than the middle piers. Similarly, the leeward side of the roadbed witnesses more deposition compared to the windward side. Another finding reveals a reduced sand deposition length as the vertical profile, in alignment with the wind direction, moves further from the bridge abutments at the same clearance height. As wind speeds rise, there’s a decline in sand deposition and a marked increase in erosion around the side piers, middle piers and roadbeds. In conclusion, a bridge clearance that’s too low can cause intense sand damage near the side piers, while an extremely high roadbed may lead to extensive surface sand deposition. Hence, railway bridges in areas prone to sandy winds should strike a balance in clearance height. This research provides valuable guidelines for determining the most suitable bridge and roadbed heights in regions affected by wind and sand.展开更多
It is of great practical value to explore the correlation between the vertical curve radius of desert highway and the increase of sand accumulation in local lines,and to select the appropriate vertical curve radius fo...It is of great practical value to explore the correlation between the vertical curve radius of desert highway and the increase of sand accumulation in local lines,and to select the appropriate vertical curve radius for reducing the risk of sand accumulation.In this study,three-dimensional models of desert highway embankments with different vertical curve radii were constructed,and Fluent software was used to simulate the wind-sand flow field and sand accumulation distribution of vertical curve embankments.The results show that:(1)Along the direction of the road,the concave and the convex vertical curve embankments have the effect of collecting and diverging the wind-sand flow,respectively.When the radius of the concave vertical curve is 3000 m,5000 m,8000 m,10000 m and 20000 m,the wind velocity in the middle of the vertical curve is 31.76%,22.58%,10.78%,10.53%and 10.44%,higher than that at both ends.When the radius of the convex vertical curve is 6500 m,8000 m,10000 m,20000 m and 30000 m,the wind velocity at both ends of the vertical curve is 14.06%,9.99%,6.14%,3.22%and 2.41%,higher than that in the middle.The diversion effect also decreases with the increase of the radius.(2)The conductivity of the concave and convex vertical curve embankments with different radii is greater than 1,which is the sediment transport roadbed.The conductivity increases with the increase of radius and gradually tends to be stable.When the radius of the concave and convex vertical curves reaches 8000 m and 20000 m respectively,the phenomenon of sand accumulation is no longer serious.Under the same radius condition,the concave vertical curve embankment is more prone to sand accumulation than the convex one.(3)Considering the strength of the collection and diversion of the vertical curve embankment with different radii,and the sand accumulation of the vertical curve embankment in the desert section of Wuma Expressway,the radius of the concave vertical curve is not less than 8000 m,and the radius of the convex vertical curve is not less than 20000 m,which can effectively reduce the sand accumulation of the vertical curve embankment.In the desert highway area,the research results of this paper can provide reference for the design of vertical curve to ensure the safe operation of desert highway.展开更多
Under the pressure of carbon neutrality,many carbon capture,utilization and storage technologies have witnessed rapid development in the recent years,including oxy-fuel combustion(OFC)technology.However,the convention...Under the pressure of carbon neutrality,many carbon capture,utilization and storage technologies have witnessed rapid development in the recent years,including oxy-fuel combustion(OFC)technology.However,the conventional OFC technology usually depends on the flue gas recirculation system,which faces significant investment,high energy consumption,and potential low-temperature corrosion problem.Considering these deficiencies,the direct utilization of pure oxygen to achieve particle fluidization and fuel combustion may reduce the overall energy consumption and CO_(2)-capture costs.In this paper,the fundamental structure of a self-designed 130 t·h^(-1) pure-oxygen combustion circulating fluidized bed(CFB)boiler was provided,and the computational particle fluid dynamics method was used to analyze the gas-solid flow characteristics of this new-concept boiler under different working conditions.The results indicate that through the careful selection of design or operational parameters,such as average bed-material size and fluidization velocity,the pure-oxygen combustion CFB system can maintain the ideal fluidization state,namely significant internal and external particle circulation.Besides,the contraction section of the boiler leads to the particle backflow in the lower furnace,resulting in the particle suspension concentration near the wall region being higher than that in the center region.Conversely,the upper furnace still retains the classic core-annulus flow structure.In addition to increasing solid circulation rate by reducing the average bed-material size,altering primary gas ratio and bed inventory can also exert varying degrees of influence on the gas-solid flow characteristics of the pure-oxygen combustion CFB boiler.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
The flow field and flow state of thin-film evaporators are complex,and it is significant to effectively divide and quantify the flow field and flow state,as well as to study the internal flow field distribution and ma...The flow field and flow state of thin-film evaporators are complex,and it is significant to effectively divide and quantify the flow field and flow state,as well as to study the internal flow field distribution and material mixing characteristics to improve the efficiency of thin-film evaporators.By using computational fluid dynamics(CFD)numerical simulation,the distribution pattern of the high-viscosity fluid flow field in the thin-film evaporators was obtained.It was found that the staggered interrupted blades could greatly promote material mixing and transportation,and impact the film formation of high-viscosity materials on the evaporator wall.Furthermore,a flow field state recognition method based on radial volume fraction statistics was proposed,and could quantitatively describe the internal flow field of thin-film evaporators.The method divides the high-viscosity materials in the thin-film evaporators into three flow states,the liquid film state,the exchange state and the liquid mass state.The three states of materials could be quantitatively described.The results show that the materials in the exchange state can connect the liquid film and the liquid mass,complete the material mixing and exchange,renew the liquid film,and maintain continuous and efficient liquid film evaporation.展开更多
As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary ...As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.展开更多
In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering comp...In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.展开更多
Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes...Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.展开更多
The fluid of casting process is a typical kind of multi-phase flow. Actually, many casting phenomena have close relationship with the multi-phase flow, such as molten metal filling process, air entrapment, slag moveme...The fluid of casting process is a typical kind of multi-phase flow. Actually, many casting phenomena have close relationship with the multi-phase flow, such as molten metal filling process, air entrapment, slag movement, venting process of die casting, gas escaping of lost foam casting and so on. Obviously, in order to analyze these phenomena accurately, numerical simulation of the multi-phase fluid is necessary. Unfortunately, so far, most of the commercial casting simulation systems do not have the ability of multi-phase flow modeling due to the difficulty in the multi-phase flow calculation. In the paper, Finite Different Method (FDM) technique was adopt to solve the multi-phase fluid model. And a simple object of the muiti-phase fluid was analyzed to obtain the fluid rates of the liquid phase and the entrapped air phase.展开更多
Based on a lattice Boltzmann method and general principles of porous flow, a numerical technique is presented for analysing the separation of multi-phase immiscible fluids in porous media. The total body force acting ...Based on a lattice Boltzmann method and general principles of porous flow, a numerical technique is presented for analysing the separation of multi-phase immiscible fluids in porous media. The total body force acting on fluid particles is modified by axiding relative permeability in Nithiarasu's expression with an axiditional surface tension term. As a test of this model, we simulate the phase separation for the case of two immiscible fluids. The numerical results show that the two coupling relative permeability coefficients K12 and K21 have the same magnitude, so the linear flux-forcing relationships satisfy Onsager reciprocity. Phase separation phenomenon is shown with the time evolution of density distribution and bears a strong similarity to the results obtained from other numerical models and the flows in sands. At the same time, the dynamical rules in this model are local, therefore it can be run on massively parallel computers with well computational efficiency.展开更多
基金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.
基金funded by Universiti Teknikal Malaysia Melaka and Ministry of Higher Education(MoHE)Malaysia,grant number FRGS/1/2024/FTKM/F00586.
文摘Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical simulation study on the dynamics of blood flow in a stenosed artery,focusing on the effects of copper and alumina nanoparticles,is conducted.The study employs a 2-dimensional Newtonian blood flow model infused with copper and alumina nanoparticles,considering the influence of a magnetic field,thermal radiation,and various flow parameters.The governing differential equations are first non-dimensionalized to facilitate analysis and subsequently solved using the 4th order collocation method,bvp4c module in MATLAB.This approach obtains velocity and temperature profiles,revealing the impact of relevant parameters crucial in the biomedical field.The findings of this study underscore the significance of understanding blood flow dynamics in stenosed arteries and the potential benefits of utilizing copper and alumina nanoparticles in treatment strategies.The incorporation of nanoparticles introduces novel avenues for enhancing therapeutic interventions,particularly in mitigating the effects of stenosis.The elucidation of velocity and temperature profiles provides valuable insights into the behavior of blood flow under different conditions,thereby informing the development of targeted biomedical applications.The arterial curvature flow parameter influences temperature profiles,with increased parameters promoting more efficient heat dissipation.The elevated values of Prandtl number and thermal radiation parameter showcase the diminished temperature profiles,indicating stronger dominance of momentum diffusion over thermal diffusion and radiative heat transfer mechanism.Sensitivity analysis of the pertinent physical parameters reveals that the Prandtl number has the most significant impact on blood flow dynamics.A statistical analysis of the present results and existing literature has also been included in the study.Overall,this research contributes to advancing our understanding of vascular health and lays the groundwork for innovative approaches in stenosis treatment and related biomedical fields.
文摘A mild gasification process has been implemented to provide an alternative form of clean coal technology called the Integrated Mild Gasification Combined Cycle (IMGCC), which can be utilized to build a new, highly efficient, and compact power plant or to retrofit an existing coal-fired power plant in order to achieve lower emissions and significantly improved thermal efficiency. The core technology of the mild gasification power plant lies on the design of a compact and effective mild gasifier that can produce synthesis gases with high energy volatiles through a hybrid system: utilizing the features of both entrained-flow and fluidized bed gasifiers. To aid in the design of the mild gasifier, a computational model has been implemented to investigate the thermal-flow and gasification process inside this mild gasifier using the commercial CFD (Computational Fluid Dynamics) solver ANSYS/FLUENT. The Eulerian-Eulerian method is employed to model both the primary phase (air) and the secondary phase (coal particles). However, the Eulerian-Eulerian model used in the software does not facilitate any built-in devolatilization model. The objective of this study is therefore to implement a devolatilization model (along with demoisturization) and incorporate it into the existing code. The Navier-Stokes equations and seven species transport equations are solved with three heterogeneous (gas-solid) and two homogeneous (gas-gas) global gasification reactions. Implementation of the complete model starts from adding demoisturization first, then devolatilization, and then adding one chemical equation at a time until finally all reactions are included in the multiphase flow. The result shows that the demoisturization and devolatilization models are successfully incorporated and a large amount of volatiles are preserved as high-energy fuels in the syngas stream without being further cracked or reacted into lighter gases. The overall results are encouraging but require future experimental data for verification.
基金financially supported by the Chang Jiang Scholar and Innovation Team Development Plan of China (IRT_15R29)the Basic Research Innovation Group Project of Gansu Province, China (21JR7RA347)the Natural Science Foundation of Gansu Province, China (20JR10RA231)。
文摘Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.
基金supported by the project of the China Geological Survey(No.DD20221746)the National Natural Science Foundation of China(Grant Nos.41101086)。
文摘Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the effectiveness of the debris flow control project and evaluated the debris flow hazards.Through field investigation and numerical simulation methods,the indicators of flow intensity reduction rate and storage capacity fullness were proposed to quantify the effectiveness of the engineering measures in the debris flow event.The simulation results show that the debris flow control project reduced the flow intensity by41.05%to 64.61%.The storage capacity of the dam decreases gradually from upstream to the mouth of the gully,thus effectively intercepting and controlling the debris flow.By evaluating the debris flow of different recurrence intervals,further measures are recommended for managing debris flow events.
基金the National Key R&D Program(No.2023YFB3709900)the National Natural Science Foundation of China(Nos.U22A20171 and 52104343)the High Steel Central(HSC)at North China University of Science and Technology and Yanshan Univ ersity,China。
文摘A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.
基金funded by the Zhejiang Provincial Natural Science Foundation of China(Grant Nos.LY21E060004,LGG22E060011)National Natural Science Foundation of China(Grant No.51976193).
文摘A coupled Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)approach is used to calculate the interaction of a flexible rag transported by a fluid current with a fixed solid cylinder.More specifically a hybrid Eulerian-Lagrangian approach is used with the rag being modeled as a set of interconnected particles.The influence of various parameters is considered,namely the inlet velocity(1.5,2.0,and 2.5 m/s,respectively),the angle formed by the initially straight rag with the flow direction(45°,60°and 90°,respectively),and the inlet position(90,100,and 110 mm,respectively).The results show that the flow rate has a significant impact on the permeability of the rag.The higher the flow rate,the higher the permeability and the rag speed difference.The angle has a minor effect on rag permeability,with 45°being the most favorable angle for permeability.The inlet position has a small impact on rag permeability,while reducing the initial distance between the rag an the cylinder makes it easier for rags to pass through.
基金financially supported by the National Natural Science Foundation of China(No.51674078)。
文摘Casting speed,casting temperature and secondary cooling water flow rate are the main process parameters affecting the DC casting process.These parameters significantly influence the flow and temperature fields during casting,which are crucial for the quality of the ingot and can determine the success or failure of the casting operation.Numerical simulation,with the advantages of low cost,rapid execution,and visualized results,is an important method to study and optimize the DC casting process.In the present work,a simulation model of DC casting 2024 aluminum alloy was established,and the reliability of the model was verified.Then,the influence of casting parameters on flow field and temperature field was studied in detail by numerical simulation method.Results show that with the increase of casting speed,the melt flow becomes faster,the depths of slurry zone and mushy zone increase,and the variation of slurry zone depth is greater than that of mushy zone.With an increase in casting temperature,the melt flow rate increases,the depth of the slurry zone becomes shallower,and the depth of the mushy zone experiences only minor changes.The simulation results further indicate that the increase of the flow rate of the secondary cooling water slightly reduces the depths of both slurry and mushy zone.
基金financially supported by the fellowship of the China Postdoctoral Science Foundation (2021M703466)the Natural Science Foundation of Gansu Province, China (20JR10RA231)the Natural Science Foundation of Gansu Province, China (22JR5RA050)。
文摘Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flow field around roadbeds, overlooking detailed examinations of sand particle erosion and deposition patterns near bridges and roadbeds. This study employs numerical simulations to analyze the influence of varying heights and wind speeds on sand deposition and erosion characteristics at different locations: the bridge-road transition section(side piers), middle piers, and roadbeds. The results show that the side piers, experience greater accumulation than the middle piers. Similarly, the leeward side of the roadbed witnesses more deposition compared to the windward side. Another finding reveals a reduced sand deposition length as the vertical profile, in alignment with the wind direction, moves further from the bridge abutments at the same clearance height. As wind speeds rise, there’s a decline in sand deposition and a marked increase in erosion around the side piers, middle piers and roadbeds. In conclusion, a bridge clearance that’s too low can cause intense sand damage near the side piers, while an extremely high roadbed may lead to extensive surface sand deposition. Hence, railway bridges in areas prone to sandy winds should strike a balance in clearance height. This research provides valuable guidelines for determining the most suitable bridge and roadbed heights in regions affected by wind and sand.
基金The research described in this paper was financially supported by Youth Science Foundation Project’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)Central guide local science and technology development funds(22ZY1QA005)。
文摘It is of great practical value to explore the correlation between the vertical curve radius of desert highway and the increase of sand accumulation in local lines,and to select the appropriate vertical curve radius for reducing the risk of sand accumulation.In this study,three-dimensional models of desert highway embankments with different vertical curve radii were constructed,and Fluent software was used to simulate the wind-sand flow field and sand accumulation distribution of vertical curve embankments.The results show that:(1)Along the direction of the road,the concave and the convex vertical curve embankments have the effect of collecting and diverging the wind-sand flow,respectively.When the radius of the concave vertical curve is 3000 m,5000 m,8000 m,10000 m and 20000 m,the wind velocity in the middle of the vertical curve is 31.76%,22.58%,10.78%,10.53%and 10.44%,higher than that at both ends.When the radius of the convex vertical curve is 6500 m,8000 m,10000 m,20000 m and 30000 m,the wind velocity at both ends of the vertical curve is 14.06%,9.99%,6.14%,3.22%and 2.41%,higher than that in the middle.The diversion effect also decreases with the increase of the radius.(2)The conductivity of the concave and convex vertical curve embankments with different radii is greater than 1,which is the sediment transport roadbed.The conductivity increases with the increase of radius and gradually tends to be stable.When the radius of the concave and convex vertical curves reaches 8000 m and 20000 m respectively,the phenomenon of sand accumulation is no longer serious.Under the same radius condition,the concave vertical curve embankment is more prone to sand accumulation than the convex one.(3)Considering the strength of the collection and diversion of the vertical curve embankment with different radii,and the sand accumulation of the vertical curve embankment in the desert section of Wuma Expressway,the radius of the concave vertical curve is not less than 8000 m,and the radius of the convex vertical curve is not less than 20000 m,which can effectively reduce the sand accumulation of the vertical curve embankment.In the desert highway area,the research results of this paper can provide reference for the design of vertical curve to ensure the safe operation of desert highway.
基金supported by the National Key Research and Development Program of China(2022YFB4100305).
文摘Under the pressure of carbon neutrality,many carbon capture,utilization and storage technologies have witnessed rapid development in the recent years,including oxy-fuel combustion(OFC)technology.However,the conventional OFC technology usually depends on the flue gas recirculation system,which faces significant investment,high energy consumption,and potential low-temperature corrosion problem.Considering these deficiencies,the direct utilization of pure oxygen to achieve particle fluidization and fuel combustion may reduce the overall energy consumption and CO_(2)-capture costs.In this paper,the fundamental structure of a self-designed 130 t·h^(-1) pure-oxygen combustion circulating fluidized bed(CFB)boiler was provided,and the computational particle fluid dynamics method was used to analyze the gas-solid flow characteristics of this new-concept boiler under different working conditions.The results indicate that through the careful selection of design or operational parameters,such as average bed-material size and fluidization velocity,the pure-oxygen combustion CFB system can maintain the ideal fluidization state,namely significant internal and external particle circulation.Besides,the contraction section of the boiler leads to the particle backflow in the lower furnace,resulting in the particle suspension concentration near the wall region being higher than that in the center region.Conversely,the upper furnace still retains the classic core-annulus flow structure.In addition to increasing solid circulation rate by reducing the average bed-material size,altering primary gas ratio and bed inventory can also exert varying degrees of influence on the gas-solid flow characteristics of the pure-oxygen combustion CFB boiler.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
基金National Natural Science Foundation of China(Nos.51905089 and 52075093)Special Fund for Basic Research and Operating Costs of Central Colleges and Universities,China(No.22320D-31)Open Fund for National Key Laboratory of Tribology of Tsinghua University,China(No.SKLTKF20B05)。
文摘The flow field and flow state of thin-film evaporators are complex,and it is significant to effectively divide and quantify the flow field and flow state,as well as to study the internal flow field distribution and material mixing characteristics to improve the efficiency of thin-film evaporators.By using computational fluid dynamics(CFD)numerical simulation,the distribution pattern of the high-viscosity fluid flow field in the thin-film evaporators was obtained.It was found that the staggered interrupted blades could greatly promote material mixing and transportation,and impact the film formation of high-viscosity materials on the evaporator wall.Furthermore,a flow field state recognition method based on radial volume fraction statistics was proposed,and could quantitatively describe the internal flow field of thin-film evaporators.The method divides the high-viscosity materials in the thin-film evaporators into three flow states,the liquid film state,the exchange state and the liquid mass state.The three states of materials could be quantitatively described.The results show that the materials in the exchange state can connect the liquid film and the liquid mass,complete the material mixing and exchange,renew the liquid film,and maintain continuous and efficient liquid film evaporation.
基金supported by State Grid Ningxia Electric Power Co.,Ltd.under Grant 5229CG220006Natural Science Foundation of Ningxia Province under Grant 2022AAC03629.
文摘As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.
文摘In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.
文摘Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.
文摘The fluid of casting process is a typical kind of multi-phase flow. Actually, many casting phenomena have close relationship with the multi-phase flow, such as molten metal filling process, air entrapment, slag movement, venting process of die casting, gas escaping of lost foam casting and so on. Obviously, in order to analyze these phenomena accurately, numerical simulation of the multi-phase fluid is necessary. Unfortunately, so far, most of the commercial casting simulation systems do not have the ability of multi-phase flow modeling due to the difficulty in the multi-phase flow calculation. In the paper, Finite Different Method (FDM) technique was adopt to solve the multi-phase fluid model. And a simple object of the muiti-phase fluid was analyzed to obtain the fluid rates of the liquid phase and the entrapped air phase.
基金Project supported by the National Natural Science Foundation of China (Grant No 10302018), the Research Grants Council of the Government of the HKSAR, China (Grant No PolyU5172/020), and the Natural Science Foundation of Zhejiang Province, China (Grant No M103082).
文摘Based on a lattice Boltzmann method and general principles of porous flow, a numerical technique is presented for analysing the separation of multi-phase immiscible fluids in porous media. The total body force acting on fluid particles is modified by axiding relative permeability in Nithiarasu's expression with an axiditional surface tension term. As a test of this model, we simulate the phase separation for the case of two immiscible fluids. The numerical results show that the two coupling relative permeability coefficients K12 and K21 have the same magnitude, so the linear flux-forcing relationships satisfy Onsager reciprocity. Phase separation phenomenon is shown with the time evolution of density distribution and bears a strong similarity to the results obtained from other numerical models and the flows in sands. At the same time, the dynamical rules in this model are local, therefore it can be run on massively parallel computers with well computational efficiency.