A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient o...A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient of restitution due to non-ideal particle collisions on the simulated results were tested. It is demonstrated that the simulated result is strongly affected by the coefficient of restitution. Comparison of simulations with experiments in a small spout-fluid bed showed that an appropriate coefficient of restitution of 0.93 was necessary to simulate the flow characteristics in an underdesigned large size of spout-fluid bed coal gasifier with diameter of lm and height of 6m. The internal jet and gas/solid flow patterns at different operating conditions were obtained. The simulations show that an optimal gas/solid flow pattern for coal gasification is found when the spouting gas flow rate is equal to the fluidizing gas flow rate and the total of them is two and a half times the minimum fluidizing gas flow rate. Besides, the radial distributions of particle velocity and gas velocity show similar tendencies; the radial distributions of particle phase pressure due to particle collisions and the particle pseudo-temperature corresponding to the macroscopic kinetic energy of the random particle motion also show similar tendencies. These indicate that both gas drag force and particle collisions dominate the movement of particles.展开更多
To get fiber motion in condensing zone of compact spinning,velocity of this area is achieved by simulation,and then a bead-elastic rod fiber model is established.Based on simulation and dynamic analysis on this zone,g...To get fiber motion in condensing zone of compact spinning,velocity of this area is achieved by simulation,and then a bead-elastic rod fiber model is established.Based on simulation and dynamic analysis on this zone,governing equations are constructed and Runge-Kutta approach is used.Lastly,trajectories of fibers are calculated by specially designed Matlab procedure according to the principles mentioned above.Results show that fiber motions at different initial positions are different;X-axis velocity component makes fibers gathering on sides of suction slot;Y-axis airflow gets fibers gradually close and then stick to the surface of lattice apron.Fiber motions also reflect that the compact spinning process in condensing zone can be divided into three parts:fast convergence zone,adjustment convergence zone,and steady convergence zone.展开更多
A computational study on the flow behavior of a gas-solid injector by Eulerian approach was carried out. The gas phase was modeled with k-ε turbulent model and the particle phase was modeled with kinetic theory of gr...A computational study on the flow behavior of a gas-solid injector by Eulerian approach was carried out. The gas phase was modeled with k-ε turbulent model and the particle phase was modeled with kinetic theory of granular flow. The simulations by Eulerian two-fluid model (TFM) were compared with the corresponding results by discrete element method (DEM) and experiments. It was showed that TFM simulated results were in reasonable agreement with the experimental and DEM simulated results. Based on TFM simulations, gas-solid flow pattern, gas velocity, particle velocity and the static pressure under different driving jet velocity, backpressure and convergent section angle were obtained. The results showed that the time average axial gas velocity sharply decreased and then slightly increased to a constant value in the horizontal conveying pipe. The time average axial particle velocity increased initially and then decreased, but in the outlet region of the convergent section the particle velocity remarkably increased once more to the maximal value. As a whole, the static pressure distribution change trends were found to be independent on driving gas velocity, backpressure and convergent section angle. However, the static pressure increased with increase of convergent section angle and gas jet velocities. The difference of static pressure to backpressure increased with increasing backpressure.展开更多
In accordance to the anisotropic feature of turbulent flow, ananisotropic algebraic stress model is adopted to predict theturbulent flow field and turbulent characteristics generated by aRushton disc turbine with the ...In accordance to the anisotropic feature of turbulent flow, ananisotropic algebraic stress model is adopted to predict theturbulent flow field and turbulent characteristics generated by aRushton disc turbine with the improved inner-outer iterativeprocedure. The predicted turbulent flow is compared with experimentaldata and the simulation by the standard k-ε turbulence model. Theanisotropic algebraic stress model is found to give better predictionthan the standard k-ε turbulence model. The predicted turbulent flowfield is in accordance to experimental data and the trend of theturbulence intensity can be effectively reflected in the simulation.展开更多
The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these th...The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the moving cylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/ D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D = 0.5 and fv= 0.4 (where A is the forced amplitude, D is the cylinder diameter, and fv stands for the vibrating frequency, respectively). The calculated results not only indicate that the spacing ratios T/D (T is the center-to-center spacing between the two upstream cylinders) have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D (where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder).展开更多
A 3D numerical investigation has been carried out to examine periodic laminar flow and heat transfer characteristics in a circular tube with 45°V-baffles with isothermal wall.The computations are based on the fin...A 3D numerical investigation has been carried out to examine periodic laminar flow and heat transfer characteristics in a circular tube with 45°V-baffles with isothermal wall.The computations are based on the finite volume method(FVM),and the SIMPLE algorithm has been implemented.The fluid flow and heat transfer characteristics are presented for Reynolds numbers ranging from 100 to 2000.To generate main longitudinal vortex flows through the tested section,V-baffles with an attack angle of 45°are mounted in tandem and in-line arrangement on the opposite positions of the circular tube.Effects of tube blockage ratio,flow direction on heat transfer and pressure drop in the tube are studied.It is apparent that a pair of longitudinal twisted vortices(P-vortex)created by a V-baffle can induce impingement on a wall of the inter-baffle cavity and lead a drastic increase in heat transfer rate at tube wall.In addition,the larger blockage ratio results in the higher Nusselt number and friction factor values.The computational results show that the optimum thermal enhancement factor is around 3.20 at baffle height of B=0.20 and B=0.25 times of the tube diameter for the V-upstream and V-downstream,respectively.展开更多
The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of eros...The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of erosion and deposition processes during debris flow movement, no entrainment assumption is unreasonable. The numerical model considering these processes is used for simulating debris flow. Raster grid networks of a digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the numerical model is solved numerically using the Leap-frog finite difference method. Finally, the simulation results can be displayed by GIS easily and used to debris flow evaluation. To illustrate this approach, the proposed methodology is applied to the Yohutagawa debris flow that occurred on 2oth October 2010, in Amami- Oshima area, Japan. The simulation results that reproduced the movement, erosion and deposition are in good agreement with the field investigation. The effectiveness of the dam in this real-ease is also verified by this approach. Comparison with the results were simulated by other models, shows that the present coupled model is more rational and effective.展开更多
A three-dimensional transient numerical simulation was conducted to study the pressure fluctuations in low-specific-speed centrifugal pumps. The characteristics of the inner flow were investigated using the SST k-ω t...A three-dimensional transient numerical simulation was conducted to study the pressure fluctuations in low-specific-speed centrifugal pumps. The characteristics of the inner flow were investigated using the SST k-ω turbulence model. The distributions of pressure fluctuations in the impeller and the volute were recorded, and the pressure fluctuation intensity was analyzed comprehensively, at the design condition, using statistical methods. The results show that the pressure fluctuation intensity increases along the impeller streamline from the leading edge to the trailing edge. In the impeller passage, the intensity near the shroud is much higher than that near the hub at thc inlet. However, the intensity at the middle passage is almost equal to the intensity at the outlet. The pressure fluctuation intensity is the highest at the trailing edge on the pressure side and near the tongue because of the rotor-stator interaction. The distribution of pressure fluctuation intensity is symmetrical in the axial cross sections of the volute channel. However, this intensity decreases with increasing radial distance. Hence, the pressure fluctuation intensity can be reduced by modifying the geometry of the leading edge in the impeller and the tongue in the volute.展开更多
An attempt was made to extend mild combustion to forward flow furnace, such as the refinery and petrochemical tube furnace. Three dimensional numerical simulation was carried out to study the performance of this furna...An attempt was made to extend mild combustion to forward flow furnace, such as the refinery and petrochemical tube furnace. Three dimensional numerical simulation was carried out to study the performance of this furnace. The Eddy Dissipation Concept(EDC) model coupled with the reaction mechanism DRM-19 was used. The prediction showed a good agreement with the measurement. The effect of air nozzle circle(D), air nozzle diameter(d), air nozzle number(N), and air preheating temperature(Tair) on the flow, temperature and species fields, and the CO and NO emissions was investigated. The results indicate that there are four zones in the furnace, viz.: a central jet zone, an ignition zone, a combustion reaction zone, and a flue gas zone, according to the distribution profiles of H_2 CO and OH. The central jet entrains more flue gas in the furnace upstream with an increasing D while the effect of D is negligible in the downstream. The air jet momentum increases with a decreasing d or an increasing Tair, and entrains more flue gas. The effect of N is mainly identified near the burner exit. More heat is absorbed in the radiant section and less heat is discharged to the atmosphere with a decreasing d and an increasing N as evidenced by the flue gas temperature. The CO and NO emissions are less than 50 μL/L and 10 μL/L, respectively, in most of conditions.展开更多
As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics (CFD). A numerical study of one transitional flow in two-d...As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics (CFD). A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation (RANS) in this paper. Turbulence model plays a significant role in the complex flows' simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox (2006) k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition's behavior.展开更多
Heat and mass transfer between porous media and fluid is a complex coupling process, which is widely used in various fields of engineering applications, especially for natural and artificial fractures in oil and gas e...Heat and mass transfer between porous media and fluid is a complex coupling process, which is widely used in various fields of engineering applications, especially for natural and artificial fractures in oil and gas extraction. In this study, a new method is proposed to deal with the flow and heat transfer problem of steady flow in a fracture. The fluid flow in a fracture was described using the same method as Mohais, who considered a fracture as a channel with porous wall, and the perturbation method was used to solve the mathematical model. Unlike previous studies, the shear jump boundary condition proposed by Ochoa-Tapia and Whitaker was used at the interface between the fluid and porous media. The main methods were perturbation analysis and the application of shear jump boundary conditions. The influence of permeability, channel width, shear jump degree and effective dynamic viscosity on the flow and heat transfer in the channel was studied by analysing the analytical solution. The distribution of axial velocity in the channel with the change of the typical parameters and the sensitivity of the heat transfer was obtained.展开更多
This paper is concerned with experimental and numerical research on 3D flow past prismatic turbine cascade SE1050 (known in QNET network as open test case SE1050). The primary goal was to assess the influence of the...This paper is concerned with experimental and numerical research on 3D flow past prismatic turbine cascade SE1050 (known in QNET network as open test case SE1050). The primary goal was to assess the influence of the inlet velocity profile on the flow structures in the interblade channel and on the flow field parameters at the cascade exit and to compare these findings to results of numerical simulations. Investigations of 3D flow past the cascade with non-uniform inlet velocity profile were carried out both experimentally and numerically at subsonic (M2 = 0.8) and at transonic (M2 = 1.2) regime at design angle of incidence. Experimental data was obtained using a traversing device with a five-hole conical probe. Numerically, the 3D flow was simulated by open source code OpenFOAM and in-house code. Analyses of experimental data and CFD simulations have revealed the development of distinctive vortex structures resulting from non-uniform inlet velocity profile. Origin of these structures results in increased loss of kinetic energy and spanwise shift of kinetic energy loss coefficient distribution. Differences found between the subsonic and the transonic case confirm earlier findings available in the literature. Results of CFD and experiments agree reasonably well.展开更多
The stability of the centrifugal pump has not been well revealed because of the complexity of internal flow. To analyze the flow characteristics of a centrifugal pump operating at low capacity, methods of numerical si...The stability of the centrifugal pump has not been well revealed because of the complexity of internal flow. To analyze the flow characteristics of a centrifugal pump operating at low capacity, methods of numerical simulation and experimental research were adopted in this paper. Characteristics of the inner flow were obtained. Standard k-s turbulence models were used to calculate the inner flow of the pump under off-design conditions. The distri- bution of the energy gradient function K was obtained by three-dimensional numerical simulation at different flow rates. The relative velocity component was acquired from the absolute velocity obtained in particle image velocimetry. By comparing with experimental results, it was found that flow instability occurs at the position of maximum K. The flow stability reduces with an increasing flow rate. The research results provide a theoretical basis for the optimization design of a centrifugal pump.展开更多
The paper deals with development and application the numerical model for solution of processes at combustion chamber of the thermal power plant boiler. Mathematical simulation is based on solution of physical and chem...The paper deals with development and application the numerical model for solution of processes at combustion chamber of the thermal power plant boiler. Mathematical simulation is based on solution of physical and chemical processes occuring at burning pulverized coal in the furnace model. Three-dimensional flows, heat and mass transfer, chemical kinetics of the processes, effects of thermal radiation are considered. Obtained results give quantitative information on velocity distributions, temperature and concentration profiles of the components, the amount of combustion products including harmful substances. The numerical model becomes a tool for investigation and design of combustion chambers with high-efficiency and reliable operation of boiler at thermal power plants.展开更多
The blade tip clearance flow in axial-flow pump is simulated based on three-dimensional N-S equations, RNG k -e turbulence model, and SIMPLEC algorithm. It shows that numerical results agree well with experiment data ...The blade tip clearance flow in axial-flow pump is simulated based on three-dimensional N-S equations, RNG k -e turbulence model, and SIMPLEC algorithm. It shows that numerical results agree well with experiment data measured by 5-hole probe through validation. Flow fields at the blade tip and velocity distribution at the exit of rotor are analyzed in detail. The numerical results show that the increase in tip clearance reduces hydro-head, especially at small flow rate. Experiment equipment is also introduced.展开更多
A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physica...A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physical domain, and the particle streaming process is carried out in the corresponding computational domain, which is transferred from the physical domain using interpolation method. For the interpolation calculation, a second-order upwind interpolation method is adopted on internal lattice nodes in flow fields while a second-order central interpolation algorithm is employed at neighbor-boundary lattice nodes. Then the above-mentioned model and algorithms are used to numerically simulate the 3D flows in the lid-driven cavity at Reynolds numbers of 100, 400 and 1000 on non-uniform meshes. Various vortices on the x-y, y-z and x-z symmetrical planes are successfully predicted, and their changes in position with the Reynolds number increasing are obtained. The velocity profiles of u component along the vertical centerline and w component along the horizontal centerline are both in good agreement with the data in literature and the calculated results on uniform meshes. Besides, the velocity vector distributions on various cross sections in lid-driven cavity predicted on non-uniform meshes are compared with those simulated on uniform meshes and those in the literature. All the comparisons and validations show that the 3D lattice Boltzmann model and all the numerical algorithms on non-uniform meshes are accurate and reliable to predict effectively flow fields.展开更多
How control in turbomachinery is very difficult because of the complexity of its fully 3-D flow structure. The authors propose to introduce streamwise vortices into the control of internal flows. A simple configuratio...How control in turbomachinery is very difficult because of the complexity of its fully 3-D flow structure. The authors propose to introduce streamwise vortices into the control of internal flows. A simple configuration of vortices was investigated in order to better understand the flow control methods by means of streamwise vortices. The research presented here concerns streamwise vortex interaction with a horseshoe vortex. The effects of such an interaction are significantly dependent on the relative location of the streamwise vortex in respect to the leading edge of the profile. The streamwise vortex is induced by an air jet. The horseshoe vortex is generated by the leading edge of a symmetric profile. Such a configuration gives possibility to investigate the interaction of these two vortices alone. The presented analysis is based on numerical simulations by means of N-S compressible solver with a two-equation turbulence model.展开更多
基金Supported by the National Key Program of Basic Research in China (No.2004CB217702, No.2005CB221202, No.2006CB20030201) and the National Natural Science Foundation of China (No.20590367, No.50676021, No.50606006).
文摘A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient of restitution due to non-ideal particle collisions on the simulated results were tested. It is demonstrated that the simulated result is strongly affected by the coefficient of restitution. Comparison of simulations with experiments in a small spout-fluid bed showed that an appropriate coefficient of restitution of 0.93 was necessary to simulate the flow characteristics in an underdesigned large size of spout-fluid bed coal gasifier with diameter of lm and height of 6m. The internal jet and gas/solid flow patterns at different operating conditions were obtained. The simulations show that an optimal gas/solid flow pattern for coal gasification is found when the spouting gas flow rate is equal to the fluidizing gas flow rate and the total of them is two and a half times the minimum fluidizing gas flow rate. Besides, the radial distributions of particle velocity and gas velocity show similar tendencies; the radial distributions of particle phase pressure due to particle collisions and the particle pseudo-temperature corresponding to the macroscopic kinetic energy of the random particle motion also show similar tendencies. These indicate that both gas drag force and particle collisions dominate the movement of particles.
文摘To get fiber motion in condensing zone of compact spinning,velocity of this area is achieved by simulation,and then a bead-elastic rod fiber model is established.Based on simulation and dynamic analysis on this zone,governing equations are constructed and Runge-Kutta approach is used.Lastly,trajectories of fibers are calculated by specially designed Matlab procedure according to the principles mentioned above.Results show that fiber motions at different initial positions are different;X-axis velocity component makes fibers gathering on sides of suction slot;Y-axis airflow gets fibers gradually close and then stick to the surface of lattice apron.Fiber motions also reflect that the compact spinning process in condensing zone can be divided into three parts:fast convergence zone,adjustment convergence zone,and steady convergence zone.
基金Supported by the National High Technology Research and Development Program of China (2006AA05A103), the National Natural Science Foundation of China (50706007), Foundation of Graduate Creative Program of Jiangsu (CX08B-060Z), and the Foundation for Excellent Ph.D. Thesis of Southeast University. ACKNOWLEDGEMENTS The authors also expressed sincere gratitude to Professors M. Horio, B. Leckner, A. Kane and E.J. Anthony for constructive advice during their visiting period in Southeast University, which contributed to our research.
文摘A computational study on the flow behavior of a gas-solid injector by Eulerian approach was carried out. The gas phase was modeled with k-ε turbulent model and the particle phase was modeled with kinetic theory of granular flow. The simulations by Eulerian two-fluid model (TFM) were compared with the corresponding results by discrete element method (DEM) and experiments. It was showed that TFM simulated results were in reasonable agreement with the experimental and DEM simulated results. Based on TFM simulations, gas-solid flow pattern, gas velocity, particle velocity and the static pressure under different driving jet velocity, backpressure and convergent section angle were obtained. The results showed that the time average axial gas velocity sharply decreased and then slightly increased to a constant value in the horizontal conveying pipe. The time average axial particle velocity increased initially and then decreased, but in the outlet region of the convergent section the particle velocity remarkably increased once more to the maximal value. As a whole, the static pressure distribution change trends were found to be independent on driving gas velocity, backpressure and convergent section angle. However, the static pressure increased with increase of convergent section angle and gas jet velocities. The difference of static pressure to backpressure increased with increasing backpressure.
基金the National Natural Science Foundation of China (No. 29792074).
文摘In accordance to the anisotropic feature of turbulent flow, ananisotropic algebraic stress model is adopted to predict theturbulent flow field and turbulent characteristics generated by aRushton disc turbine with the improved inner-outer iterativeprocedure. The predicted turbulent flow is compared with experimentaldata and the simulation by the standard k-ε turbulence model. Theanisotropic algebraic stress model is found to give better predictionthan the standard k-ε turbulence model. The predicted turbulent flowfield is in accordance to experimental data and the trend of theturbulence intensity can be effectively reflected in the simulation.
基金Support by the National Natural Science Foundation of China under Grant Nos.10932010 and 11072220the Natural Science Foundation of Zhejiang Province under Grant Nos.Y607425,Z6090556the Foundation Project for Youths of Zhijiang Normal University under Grant No.KJ20090102
文摘The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the moving cylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/ D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D = 0.5 and fv= 0.4 (where A is the forced amplitude, D is the cylinder diameter, and fv stands for the vibrating frequency, respectively). The calculated results not only indicate that the spacing ratios T/D (T is the center-to-center spacing between the two upstream cylinders) have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D (where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder).
基金Supported by the King Mongkut's Institute of Technology Ladkrabang research fund,Thailand(KREF015611)
文摘A 3D numerical investigation has been carried out to examine periodic laminar flow and heat transfer characteristics in a circular tube with 45°V-baffles with isothermal wall.The computations are based on the finite volume method(FVM),and the SIMPLE algorithm has been implemented.The fluid flow and heat transfer characteristics are presented for Reynolds numbers ranging from 100 to 2000.To generate main longitudinal vortex flows through the tested section,V-baffles with an attack angle of 45°are mounted in tandem and in-line arrangement on the opposite positions of the circular tube.Effects of tube blockage ratio,flow direction on heat transfer and pressure drop in the tube are studied.It is apparent that a pair of longitudinal twisted vortices(P-vortex)created by a V-baffle can induce impingement on a wall of the inter-baffle cavity and lead a drastic increase in heat transfer rate at tube wall.In addition,the larger blockage ratio results in the higher Nusselt number and friction factor values.The computational results show that the optimum thermal enhancement factor is around 3.20 at baffle height of B=0.20 and B=0.25 times of the tube diameter for the V-upstream and V-downstream,respectively.
基金finanicial support from the Global Environment Research Fund of Japan(S-8)from Grants-in-Aid for Scientific Research(Scientific Research(B),22310113,G.Chen)from Japan Society for the Promotion of Science
文摘The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of erosion and deposition processes during debris flow movement, no entrainment assumption is unreasonable. The numerical model considering these processes is used for simulating debris flow. Raster grid networks of a digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the numerical model is solved numerically using the Leap-frog finite difference method. Finally, the simulation results can be displayed by GIS easily and used to debris flow evaluation. To illustrate this approach, the proposed methodology is applied to the Yohutagawa debris flow that occurred on 2oth October 2010, in Amami- Oshima area, Japan. The simulation results that reproduced the movement, erosion and deposition are in good agreement with the field investigation. The effectiveness of the dam in this real-ease is also verified by this approach. Comparison with the results were simulated by other models, shows that the present coupled model is more rational and effective.
基金Projects(51239005,51009072) supported by the National Natural Science Foundation of ChinaProject(2011BAF14B04) supported by the National Science&Technology Pillar Program of ChinaProject(13JDG084) supported by the Research Foundation for Advanced Talents of Jiansu University,China
文摘A three-dimensional transient numerical simulation was conducted to study the pressure fluctuations in low-specific-speed centrifugal pumps. The characteristics of the inner flow were investigated using the SST k-ω turbulence model. The distributions of pressure fluctuations in the impeller and the volute were recorded, and the pressure fluctuation intensity was analyzed comprehensively, at the design condition, using statistical methods. The results show that the pressure fluctuation intensity increases along the impeller streamline from the leading edge to the trailing edge. In the impeller passage, the intensity near the shroud is much higher than that near the hub at thc inlet. However, the intensity at the middle passage is almost equal to the intensity at the outlet. The pressure fluctuation intensity is the highest at the trailing edge on the pressure side and near the tongue because of the rotor-stator interaction. The distribution of pressure fluctuation intensity is symmetrical in the axial cross sections of the volute channel. However, this intensity decreases with increasing radial distance. Hence, the pressure fluctuation intensity can be reduced by modifying the geometry of the leading edge in the impeller and the tongue in the volute.
基金supported by the technology development fund of China Petroleum & Chemical Corporation (Sinopec 312016 and 314054)
文摘An attempt was made to extend mild combustion to forward flow furnace, such as the refinery and petrochemical tube furnace. Three dimensional numerical simulation was carried out to study the performance of this furnace. The Eddy Dissipation Concept(EDC) model coupled with the reaction mechanism DRM-19 was used. The prediction showed a good agreement with the measurement. The effect of air nozzle circle(D), air nozzle diameter(d), air nozzle number(N), and air preheating temperature(Tair) on the flow, temperature and species fields, and the CO and NO emissions was investigated. The results indicate that there are four zones in the furnace, viz.: a central jet zone, an ignition zone, a combustion reaction zone, and a flue gas zone, according to the distribution profiles of H_2 CO and OH. The central jet entrains more flue gas in the furnace upstream with an increasing D while the effect of D is negligible in the downstream. The air jet momentum increases with a decreasing d or an increasing Tair, and entrains more flue gas. The effect of N is mainly identified near the burner exit. More heat is absorbed in the radiant section and less heat is discharged to the atmosphere with a decreasing d and an increasing N as evidenced by the flue gas temperature. The CO and NO emissions are less than 50 μL/L and 10 μL/L, respectively, in most of conditions.
基金Foundation item: Supported by the National Natural Science Foundation of China (Nos. 51309040, 51379025), and the Fundamental Research Funds for the Central Universities (Nos. 3132014224, 3132014318).
文摘As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics (CFD). A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation (RANS) in this paper. Turbulence model plays a significant role in the complex flows' simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox (2006) k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition's behavior.
基金financially supported by National Natural Science Foundation of China(Grant No.51305238)
文摘Heat and mass transfer between porous media and fluid is a complex coupling process, which is widely used in various fields of engineering applications, especially for natural and artificial fractures in oil and gas extraction. In this study, a new method is proposed to deal with the flow and heat transfer problem of steady flow in a fracture. The fluid flow in a fracture was described using the same method as Mohais, who considered a fracture as a channel with porous wall, and the perturbation method was used to solve the mathematical model. Unlike previous studies, the shear jump boundary condition proposed by Ochoa-Tapia and Whitaker was used at the interface between the fluid and porous media. The main methods were perturbation analysis and the application of shear jump boundary conditions. The influence of permeability, channel width, shear jump degree and effective dynamic viscosity on the flow and heat transfer in the channel was studied by analysing the analytical solution. The distribution of axial velocity in the channel with the change of the typical parameters and the sensitivity of the heat transfer was obtained.
基金supported by the Czech Science Foundation under the grant No. GAP 101/10/1329
文摘This paper is concerned with experimental and numerical research on 3D flow past prismatic turbine cascade SE1050 (known in QNET network as open test case SE1050). The primary goal was to assess the influence of the inlet velocity profile on the flow structures in the interblade channel and on the flow field parameters at the cascade exit and to compare these findings to results of numerical simulations. Investigations of 3D flow past the cascade with non-uniform inlet velocity profile were carried out both experimentally and numerically at subsonic (M2 = 0.8) and at transonic (M2 = 1.2) regime at design angle of incidence. Experimental data was obtained using a traversing device with a five-hole conical probe. Numerically, the 3D flow was simulated by open source code OpenFOAM and in-house code. Analyses of experimental data and CFD simulations have revealed the development of distinctive vortex structures resulting from non-uniform inlet velocity profile. Origin of these structures results in increased loss of kinetic energy and spanwise shift of kinetic energy loss coefficient distribution. Differences found between the subsonic and the transonic case confirm earlier findings available in the literature. Results of CFD and experiments agree reasonably well.
基金supported by Key Projects in the National Science & Technology Pillar Program (No. 2013 BAF05B01)Public Projects of Zhejiang Province (No. 2014C31116)+1 种基金521 Talent Project of Zhejiang Sci-Tech University, Fluid Engineering Innovation Team of Zhejiang Sci-Tech University (No. 11132932611309)National Natural Science Foundation of China (No. 51409233)
文摘The stability of the centrifugal pump has not been well revealed because of the complexity of internal flow. To analyze the flow characteristics of a centrifugal pump operating at low capacity, methods of numerical simulation and experimental research were adopted in this paper. Characteristics of the inner flow were obtained. Standard k-s turbulence models were used to calculate the inner flow of the pump under off-design conditions. The distri- bution of the energy gradient function K was obtained by three-dimensional numerical simulation at different flow rates. The relative velocity component was acquired from the absolute velocity obtained in particle image velocimetry. By comparing with experimental results, it was found that flow instability occurs at the position of maximum K. The flow stability reduces with an increasing flow rate. The research results provide a theoretical basis for the optimization design of a centrifugal pump.
基金funded by the Ministry of Education and Science of Kazakhstan Republic,№0112РК01095support from the Technology Agency of the Czech Republic in the frame of the Competence Centre Advanced Technology of Heat and Electricity Output,No.TE01020036
文摘The paper deals with development and application the numerical model for solution of processes at combustion chamber of the thermal power plant boiler. Mathematical simulation is based on solution of physical and chemical processes occuring at burning pulverized coal in the furnace model. Three-dimensional flows, heat and mass transfer, chemical kinetics of the processes, effects of thermal radiation are considered. Obtained results give quantitative information on velocity distributions, temperature and concentration profiles of the components, the amount of combustion products including harmful substances. The numerical model becomes a tool for investigation and design of combustion chambers with high-efficiency and reliable operation of boiler at thermal power plants.
文摘The blade tip clearance flow in axial-flow pump is simulated based on three-dimensional N-S equations, RNG k -e turbulence model, and SIMPLEC algorithm. It shows that numerical results agree well with experiment data measured by 5-hole probe through validation. Flow fields at the blade tip and velocity distribution at the exit of rotor are analyzed in detail. The numerical results show that the increase in tip clearance reduces hydro-head, especially at small flow rate. Experiment equipment is also introduced.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51179192, 50779069, 51139007)the Program for New Century Excellent Talents in University (NCET) (Grant No. NETC-10-0784)+1 种基金the National Hi-Tech Research and Development Program of China ("863" Project) (Grant No. 2011AA100505)the Chinese Universities Scientific Fund (Grant No. 2013RC045)
文摘A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physical domain, and the particle streaming process is carried out in the corresponding computational domain, which is transferred from the physical domain using interpolation method. For the interpolation calculation, a second-order upwind interpolation method is adopted on internal lattice nodes in flow fields while a second-order central interpolation algorithm is employed at neighbor-boundary lattice nodes. Then the above-mentioned model and algorithms are used to numerically simulate the 3D flows in the lid-driven cavity at Reynolds numbers of 100, 400 and 1000 on non-uniform meshes. Various vortices on the x-y, y-z and x-z symmetrical planes are successfully predicted, and their changes in position with the Reynolds number increasing are obtained. The velocity profiles of u component along the vertical centerline and w component along the horizontal centerline are both in good agreement with the data in literature and the calculated results on uniform meshes. Besides, the velocity vector distributions on various cross sections in lid-driven cavity predicted on non-uniform meshes are compared with those simulated on uniform meshes and those in the literature. All the comparisons and validations show that the 3D lattice Boltzmann model and all the numerical algorithms on non-uniform meshes are accurate and reliable to predict effectively flow fields.
文摘How control in turbomachinery is very difficult because of the complexity of its fully 3-D flow structure. The authors propose to introduce streamwise vortices into the control of internal flows. A simple configuration of vortices was investigated in order to better understand the flow control methods by means of streamwise vortices. The research presented here concerns streamwise vortex interaction with a horseshoe vortex. The effects of such an interaction are significantly dependent on the relative location of the streamwise vortex in respect to the leading edge of the profile. The streamwise vortex is induced by an air jet. The horseshoe vortex is generated by the leading edge of a symmetric profile. Such a configuration gives possibility to investigate the interaction of these two vortices alone. The presented analysis is based on numerical simulations by means of N-S compressible solver with a two-equation turbulence model.