According to fluid dynamics analysis during the fire, the criteria k-ε two-equation model for solving three-dimensional turbulence was determined, the pollutants generated in the fire disaster were set by adopting Mi...According to fluid dynamics analysis during the fire, the criteria k-ε two-equation model for solving three-dimensional turbulence was determined, the pollutants generated in the fire disaster were set by adopting Mixture multiphase flow, and the SIMPLE algorithm was used for solving on the basis of comprehensive consideration on the heat radia- tion and components transmission during fire. By simulating the airflow flowing state inside the tunnel during fire disaster of downward ventilation, drift ventilation, and ascensional ventilation, respectively, with regard to the actual situation of No.l, No.3, and No.5 belt roadway in Kongzhuang Coal Mine, the velocity vector distributions of pollutants under different inlet air volumes were obtained, and the damage degree and influential factors of disaster were also clear, which is helpful to control and avoid disaster during belt roadway fire.展开更多
A numerical analysis of the log-law behavior for the turbulent boundary layer of a wall-bounded flow is performed over a flat plate immersed in three nanofluids(Zn O-water,SiO_(2)-water,TiO_(2)-water).Numerical simula...A numerical analysis of the log-law behavior for the turbulent boundary layer of a wall-bounded flow is performed over a flat plate immersed in three nanofluids(Zn O-water,SiO_(2)-water,TiO_(2)-water).Numerical simulations using CFD code are employed to investigate the boundary layer and the hydrodynamic flow.To validate the current numerical model,measurement points from published works were used,and the compared results were in good compliance.Simulations were carried out for the velocity series of 0.04,0.4 and 4 m/s and nanoparticle concentrations0.1% and 5%.The influence of nanoparticles’ concentration on velocity,temperature profiles,wall shear stress,and turbulent intensity was investigated.The obtained results showed that the viscous sub-layer,the buffer layer,and the loglaw layer along the potential-flow layer could be analyzed based on their curving quality in the regions which have just a single wall distance.It was seen that the viscous sub-layer is the biggest area in comparison with other areas.Alternatively,the section where the temperature changes considerably correspond to the thermal boundary layer’s thickness goes a downward trend when the velocity decreases.The thermal boundary layer gets deep away from the leading edge.However,a rise in the volume fraction of nanoparticles indicated a minor impact on the shear stress developed in the wall.In all cases,the thickness of the boundary layer undergoes a downward trend as the velocity increases,whereas increasing the nanoparticle concentrations would enhance the thickness.More precisely,the log layer is closed with log law,and it is minimal between Y^(+)=50 and Y^(+)=95.The temperature for nanoparticle concentration φ=5%is higher than that for φ=0.1%,in boundary layers,for all studied nanofluids.However,it is established that the behavior is inverted from the value of Y^(+)=1 and the temperature for φ =0.1% is more important than the case of φ =5%.For turbulence intensity peak,this peak exists at Y^(+)=100 for v=4 m/s,Y^(+)=10 for v=0.4 m/s and Y^(+)=8 for v=0.04 m/s.展开更多
Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numeri...Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture muitiphase model to simulate the bubbly wake around the KCS hull. The realizable k-e two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.展开更多
In this study,a suitable CFD(computational fluid dynamics)model has been developed to investigate the influence of liquid height on the discharge coefficient of the orifice-type liquid distributors.The orifice flow in...In this study,a suitable CFD(computational fluid dynamics)model has been developed to investigate the influence of liquid height on the discharge coefficient of the orifice-type liquid distributors.The orifice flow in different diameters and liquid heights has been realized using the shear stress transport(SST)turbulence model and the Gamma Theta transition(GTT)model.In the ANSYS CFX software,two models are used in conjunction with an automatic wall treatment which allows for a smooth shift from a wall function(WF)to a low turbulent-Re near wall formulation(LTRW).The results of the models coupled with LTRW are closer to the experimental results compared with the models with WF,indicating that LTRW is more appropriate for the prediction of boundary layer characteristics of orifice flow.Simulation results show that the flow conditions of orifices change with the variation of liquid height.With respect to the turbulence in orifice,the SST model coupled with LTRW is recommended.However,with respect to the transition to turbulence in orifice with an increase in liquid height,the predictions of GTT model coupled with LTRW are superior to those obtained using other models.展开更多
The kinematical characteristics and thermophoretic deposition of inhalable particles with the diameters of 0-2.5μm (hereafter referred to as PM2.5) suspended in turbulent air flow in a rectangular duct with tempera...The kinematical characteristics and thermophoretic deposition of inhalable particles with the diameters of 0-2.5μm (hereafter referred to as PM2.5) suspended in turbulent air flow in a rectangular duct with temperature distribution were experimentally studied. Particle dynamics analyzer (PDA) was used for the on-line measurement of particle motion and particle concentration distribution in the cross-sections of the duct. The influences of the parameters such as the ratio of the bulk air temperature to the cold wall temperature and the air flow rate in the duct on the kinematical characteristics and the deposition efficiencies of PM2.5 were investigated. The experimental re- sults show that the deposition efficiencies of PM2.5 mainly depend on the temperature difference between the air and the cold wail, wffile the air flow rate and the particlecon^centration almost affect hardly tile clep0si-tion-effi ciency. The radial force thermophoresis to push PM2.5 to the cold wail is found the key factor for PM2.5 deposition.Based on the experimental results, an empirical modified Romay correlation for the calculation of thermophoretic deposition efficiency of PM2.5 is presenlext. The empirical correlation agrees reasonably well with the experimental data.展开更多
Drop breakage and coalescence influence the particle formation in liquid-liquid dispersions. In order to reduce the influencing factors of the whole dispersion process, single drops where coalescence processes can be ...Drop breakage and coalescence influence the particle formation in liquid-liquid dispersions. In order to reduce the influencing factors of the whole dispersion process, single drops where coalescence processes can be neglected were analyzed in this work. Drops passing the turbulent vicinity of a single stirrer blade were investi- gated by high-speed imaging. In order to gain a statistically relevant amount of drops passing the area of interest and corresponding breakage events, at least 1600 droplets were considered for each parameter set of this work. A specially developed fully automatic image analysis based on Matlab was used for the evaluation of the resulting high amount of image data. This allowed the elimination of the time-consuming manual analysis and further- more, allowed the objective evaluation of the drops' behavior. Different deformation parameters were consid- ered in order to describe the drop deformation dynamics properly. Regarding the ratio of both main particle axes (0axes), which was therefore approximated through an ellipse, allowed the determination of very small de- viations from the spherical shape. The perimeter of the particle (0peri) was used for the description of highly de- formed shapes. In this work the results of a higher viscosity paraffin oil (ηd =127 mPa. s) and a low viscosity solvent (petroleum, ηd = 1.7 mPa-s) are presented with and without the addition of SDS to the continuous water phase. All results show that the experimentally determined oscillation but also deformation times underlie a wide spreading. Drop deformations significantly increased not only with increasing droplet viscosity, but also with decreasing interfacial tension. Highly deformed particles of one droplet species were more likely to break than more or less spherical particles. As droplet fragmentation results from a variety of different macro-scale de- formed particles, it is not assumed that a critical deformation value must be reached for the fragmentation pro- cess to occur. Especially for highly deformed particles thin particle filaments are assumed to induce the breakage process and, therefore, be responsible for the separation of drops.展开更多
Wind force coefficients for designing porous canopy roofs have been investigated based on a series of wind tunnel experiments. Gable, troughed and mono-sloped roofs were tested. The roof models were made of 0.5 mm thi...Wind force coefficients for designing porous canopy roofs have been investigated based on a series of wind tunnel experiments. Gable, troughed and mono-sloped roofs were tested. The roof models were made of 0.5 mm thick perforated duralumin plates, the porosity of which was changed from 0 to about 0.4. Overall aerodynamic forces and moments acting on the roof model were measured in a turbulent boundary layer with a six-component force balance for various wind directions. The results indicate that the wind loads on canopy roofs generally decrease with an increase in porosity of the roof. Assuming that the roof is rigid and supported by the four corner columns with no walls, the axial forces induced in the columns are regarded as the most important load effect for discussing the design wind loads. Two loading patterns causing the maximum tension and compression in the columns are considered. Based on a combination of the lift and moment coefficients, the design wind force coefficients on the windward and leeward halves of the roof are presented for the two loading patterns as a function of the roof pitch and porosity. The effect of porosity is taken into account as a reduction factor of the wind loads.展开更多
According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in...According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.展开更多
In the radiant section of cracking furnace,the thermal cracking process is highly coupled with turbulent flow,heat transfer and mass transfer.In this paper,a three-dimensional simulation of propane pyrolysis reactor t...In the radiant section of cracking furnace,the thermal cracking process is highly coupled with turbulent flow,heat transfer and mass transfer.In this paper,a three-dimensional simulation of propane pyrolysis reactor tube is performed based on a detailed kinetic radical cracking scheme,combined with a comprehensive rigorous computational fluid dynamics(CFD)model.The eddy-dissipation-concept(EDC)model is introduced to deal with turbulence-chemistry interaction of cracking gas,especially for the multi-step radical kinetics.Considering the high aspect ratio and severe gradient phenomenon,numerical strategies such as grid resolution and refinement,stepping method and relaxation technique at different levels are employed to accelerate convergence.Large scale of radial nonuniformity in the vicinity of the tube wall is investigated.Spatial distributions of each radical reaction rate are first studied,and made it possible to identify the dominant elementary reactions.Additionally,a series of operating conditions including the feedstock feed rate,wall temperature profile and heat flux profile towards the reactor tubes are investigated.The obtained results can be used as scientific guide for further technical retrofit and operation optimization aiming at high conversion and selectivity of pyrolysis process.展开更多
An explicit expression for local, instantaneous NO production rate model was proposed to simulate NO formation in turbulent methane-air combustion. The average production rates of mixture fraction and scalar dissipati...An explicit expression for local, instantaneous NO production rate model was proposed to simulate NO formation in turbulent methane-air combustion. The average production rates of mixture fraction and scalar dissipation were obtained from asymptotes through approximation of two single-variable probability-density function. The theory predicted significant contributions from the Zeldovich mechanism, but negligible contributions from the nitrous-oxide mechanism in the oxygenconsumption zone. The proposed model was used to simulate NO formation in the pilot methane-air jet diffusion combustion. The simulation results were compared with those obtained by the CFD software FLUENT module. Validation of predictions with the experimental data given by Sandia National Laboratory of the USA indicates that the proposed model yields better results than other models, and the deviation is under 5%. And in some complete reaction zones, the simulation results are even the same as the experimental data. Realizable κ-ε model, Reynold stress model and standard κ-ε model were also investigated to predict the turbulent combustion reaction, which shows that the simulation results of velocities, temperatures, and concentrations of combustion productions by standard κ-ε model are in accordance with the experimental data.展开更多
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.展开更多
We derive higher order magneto-hydrodynamic (MHD) equations from a microscopic picture using pro-jection and perturbation formalism. In an application to Hartmann flow we find velocity profiles flattening towards th...We derive higher order magneto-hydrodynamic (MHD) equations from a microscopic picture using pro-jection and perturbation formalism. In an application to Hartmann flow we find velocity profiles flattening towards the center at the onset of turbulence in hydrodynamic limit. Comparison with the system under the effect of a uniform magnetic field yields difference in the onset of turbulence consistent with observations, showing that the presence of magnetic field inhibits onset of instability or turbulence. The laminar-turbulent transition is demonstrated in a phase transition plot of the development in time of the relative average velocities vs. Reynolds number showing a sharp increase of the relative average velocity at the transition point as determined by the critical Reynolds number.展开更多
In this paper, the dynamics behaviors on fo-δ parameter surface is investigated for Gledzer-Ohkitani- Yamada model We indicate the type of intermittency chaos transitions is saddle node bifurcation. We plot phase dia...In this paper, the dynamics behaviors on fo-δ parameter surface is investigated for Gledzer-Ohkitani- Yamada model We indicate the type of intermittency chaos transitions is saddle node bifurcation. We plot phase diagram on fo-δ parameter surface, which is divided into periodic, quasi-periodic, and intermittent chaos areas. By means of varying Taylor-microscale Reynolds number, we calculate the extended self-similarity of velocity structure function.展开更多
The flow field of liquid phase (water) of agitated extraction columns is simulated with the help of computational fluid dynamics (CFD). Four kinds of Reynolds-averaged turbulence models, i.e. the standard k-ε model, ...The flow field of liquid phase (water) of agitated extraction columns is simulated with the help of computational fluid dynamics (CFD). Four kinds of Reynolds-averaged turbulence models, i.e. the standard k-ε model, the RNG (renormalization group) k-s model, the realizable k-ε model and the Reynolds stress model, are compared in detail in order to judge which is the best model in terms of the accuracy, less CPU time and memory required. The performance of the realizable k-s model is obviously improved by reducing the model constant from C2 = 1.90 to C2 = 1.61. It is concluded that the improved realizable k-e model is the optimal model.展开更多
基金Supported by the International Science and Technology Cooperation Projects(2009DFA71840)Basic Research Business Projects of China Academy of Safety Science and Technology(2009JBKY07)
文摘According to fluid dynamics analysis during the fire, the criteria k-ε two-equation model for solving three-dimensional turbulence was determined, the pollutants generated in the fire disaster were set by adopting Mixture multiphase flow, and the SIMPLE algorithm was used for solving on the basis of comprehensive consideration on the heat radia- tion and components transmission during fire. By simulating the airflow flowing state inside the tunnel during fire disaster of downward ventilation, drift ventilation, and ascensional ventilation, respectively, with regard to the actual situation of No.l, No.3, and No.5 belt roadway in Kongzhuang Coal Mine, the velocity vector distributions of pollutants under different inlet air volumes were obtained, and the damage degree and influential factors of disaster were also clear, which is helpful to control and avoid disaster during belt roadway fire.
基金support he received through General Research Project under the grant number (R.G.P.2/138/42)。
文摘A numerical analysis of the log-law behavior for the turbulent boundary layer of a wall-bounded flow is performed over a flat plate immersed in three nanofluids(Zn O-water,SiO_(2)-water,TiO_(2)-water).Numerical simulations using CFD code are employed to investigate the boundary layer and the hydrodynamic flow.To validate the current numerical model,measurement points from published works were used,and the compared results were in good compliance.Simulations were carried out for the velocity series of 0.04,0.4 and 4 m/s and nanoparticle concentrations0.1% and 5%.The influence of nanoparticles’ concentration on velocity,temperature profiles,wall shear stress,and turbulent intensity was investigated.The obtained results showed that the viscous sub-layer,the buffer layer,and the loglaw layer along the potential-flow layer could be analyzed based on their curving quality in the regions which have just a single wall distance.It was seen that the viscous sub-layer is the biggest area in comparison with other areas.Alternatively,the section where the temperature changes considerably correspond to the thermal boundary layer’s thickness goes a downward trend when the velocity decreases.The thermal boundary layer gets deep away from the leading edge.However,a rise in the volume fraction of nanoparticles indicated a minor impact on the shear stress developed in the wall.In all cases,the thickness of the boundary layer undergoes a downward trend as the velocity increases,whereas increasing the nanoparticle concentrations would enhance the thickness.More precisely,the log layer is closed with log law,and it is minimal between Y^(+)=50 and Y^(+)=95.The temperature for nanoparticle concentration φ=5%is higher than that for φ=0.1%,in boundary layers,for all studied nanofluids.However,it is established that the behavior is inverted from the value of Y^(+)=1 and the temperature for φ =0.1% is more important than the case of φ =5%.For turbulence intensity peak,this peak exists at Y^(+)=100 for v=4 m/s,Y^(+)=10 for v=0.4 m/s and Y^(+)=8 for v=0.04 m/s.
文摘Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture muitiphase model to simulate the bubbly wake around the KCS hull. The realizable k-e two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.
基金the financial support from the National Basic Research Program of China(No.2009CB219905)the Program for Changjiang Scholars and Innovative Research Team in University(No.IRT0936)the National Natural Science Foundation of China(No.21176172)
文摘In this study,a suitable CFD(computational fluid dynamics)model has been developed to investigate the influence of liquid height on the discharge coefficient of the orifice-type liquid distributors.The orifice flow in different diameters and liquid heights has been realized using the shear stress transport(SST)turbulence model and the Gamma Theta transition(GTT)model.In the ANSYS CFX software,two models are used in conjunction with an automatic wall treatment which allows for a smooth shift from a wall function(WF)to a low turbulent-Re near wall formulation(LTRW).The results of the models coupled with LTRW are closer to the experimental results compared with the models with WF,indicating that LTRW is more appropriate for the prediction of boundary layer characteristics of orifice flow.Simulation results show that the flow conditions of orifices change with the variation of liquid height.With respect to the turbulence in orifice,the SST model coupled with LTRW is recommended.However,with respect to the transition to turbulence in orifice with an increase in liquid height,the predictions of GTT model coupled with LTRW are superior to those obtained using other models.
基金the Special Funds for Major State Basic Research Project of China (2002CB211604)
文摘The kinematical characteristics and thermophoretic deposition of inhalable particles with the diameters of 0-2.5μm (hereafter referred to as PM2.5) suspended in turbulent air flow in a rectangular duct with temperature distribution were experimentally studied. Particle dynamics analyzer (PDA) was used for the on-line measurement of particle motion and particle concentration distribution in the cross-sections of the duct. The influences of the parameters such as the ratio of the bulk air temperature to the cold wall temperature and the air flow rate in the duct on the kinematical characteristics and the deposition efficiencies of PM2.5 were investigated. The experimental re- sults show that the deposition efficiencies of PM2.5 mainly depend on the temperature difference between the air and the cold wail, wffile the air flow rate and the particlecon^centration almost affect hardly tile clep0si-tion-effi ciency. The radial force thermophoresis to push PM2.5 to the cold wail is found the key factor for PM2.5 deposition.Based on the experimental results, an empirical modified Romay correlation for the calculation of thermophoretic deposition efficiency of PM2.5 is presenlext. The empirical correlation agrees reasonably well with the experimental data.
基金supported by the German Research Foundation (DFG) within the project "Modelling,Simulation,and Control of Drop Size Distributions in Stirred Liquid/liquid Systems - KR1639/15-1"the "Max-Buchner-Forschungsstiftung"
文摘Drop breakage and coalescence influence the particle formation in liquid-liquid dispersions. In order to reduce the influencing factors of the whole dispersion process, single drops where coalescence processes can be neglected were analyzed in this work. Drops passing the turbulent vicinity of a single stirrer blade were investi- gated by high-speed imaging. In order to gain a statistically relevant amount of drops passing the area of interest and corresponding breakage events, at least 1600 droplets were considered for each parameter set of this work. A specially developed fully automatic image analysis based on Matlab was used for the evaluation of the resulting high amount of image data. This allowed the elimination of the time-consuming manual analysis and further- more, allowed the objective evaluation of the drops' behavior. Different deformation parameters were consid- ered in order to describe the drop deformation dynamics properly. Regarding the ratio of both main particle axes (0axes), which was therefore approximated through an ellipse, allowed the determination of very small de- viations from the spherical shape. The perimeter of the particle (0peri) was used for the description of highly de- formed shapes. In this work the results of a higher viscosity paraffin oil (ηd =127 mPa. s) and a low viscosity solvent (petroleum, ηd = 1.7 mPa-s) are presented with and without the addition of SDS to the continuous water phase. All results show that the experimentally determined oscillation but also deformation times underlie a wide spreading. Drop deformations significantly increased not only with increasing droplet viscosity, but also with decreasing interfacial tension. Highly deformed particles of one droplet species were more likely to break than more or less spherical particles. As droplet fragmentation results from a variety of different macro-scale de- formed particles, it is not assumed that a critical deformation value must be reached for the fragmentation pro- cess to occur. Especially for highly deformed particles thin particle filaments are assumed to induce the breakage process and, therefore, be responsible for the separation of drops.
文摘Wind force coefficients for designing porous canopy roofs have been investigated based on a series of wind tunnel experiments. Gable, troughed and mono-sloped roofs were tested. The roof models were made of 0.5 mm thick perforated duralumin plates, the porosity of which was changed from 0 to about 0.4. Overall aerodynamic forces and moments acting on the roof model were measured in a turbulent boundary layer with a six-component force balance for various wind directions. The results indicate that the wind loads on canopy roofs generally decrease with an increase in porosity of the roof. Assuming that the roof is rigid and supported by the four corner columns with no walls, the axial forces induced in the columns are regarded as the most important load effect for discussing the design wind loads. Two loading patterns causing the maximum tension and compression in the columns are considered. Based on a combination of the lift and moment coefficients, the design wind force coefficients on the windward and leeward halves of the roof are presented for the two loading patterns as a function of the roof pitch and porosity. The effect of porosity is taken into account as a reduction factor of the wind loads.
基金Project(51074027)supported by the National Natural Science Foundation of China
文摘According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.
基金Supported by the National Science&Technology Supporting Plan(2012BAF05B00)the National Basic Research Program(2012CB720500)
文摘In the radiant section of cracking furnace,the thermal cracking process is highly coupled with turbulent flow,heat transfer and mass transfer.In this paper,a three-dimensional simulation of propane pyrolysis reactor tube is performed based on a detailed kinetic radical cracking scheme,combined with a comprehensive rigorous computational fluid dynamics(CFD)model.The eddy-dissipation-concept(EDC)model is introduced to deal with turbulence-chemistry interaction of cracking gas,especially for the multi-step radical kinetics.Considering the high aspect ratio and severe gradient phenomenon,numerical strategies such as grid resolution and refinement,stepping method and relaxation technique at different levels are employed to accelerate convergence.Large scale of radial nonuniformity in the vicinity of the tube wall is investigated.Spatial distributions of each radical reaction rate are first studied,and made it possible to identify the dominant elementary reactions.Additionally,a series of operating conditions including the feedstock feed rate,wall temperature profile and heat flux profile towards the reactor tubes are investigated.The obtained results can be used as scientific guide for further technical retrofit and operation optimization aiming at high conversion and selectivity of pyrolysis process.
文摘An explicit expression for local, instantaneous NO production rate model was proposed to simulate NO formation in turbulent methane-air combustion. The average production rates of mixture fraction and scalar dissipation were obtained from asymptotes through approximation of two single-variable probability-density function. The theory predicted significant contributions from the Zeldovich mechanism, but negligible contributions from the nitrous-oxide mechanism in the oxygenconsumption zone. The proposed model was used to simulate NO formation in the pilot methane-air jet diffusion combustion. The simulation results were compared with those obtained by the CFD software FLUENT module. Validation of predictions with the experimental data given by Sandia National Laboratory of the USA indicates that the proposed model yields better results than other models, and the deviation is under 5%. And in some complete reaction zones, the simulation results are even the same as the experimental data. Realizable κ-ε model, Reynold stress model and standard κ-ε model were also investigated to predict the turbulent combustion reaction, which shows that the simulation results of velocities, temperatures, and concentrations of combustion productions by standard κ-ε model are in accordance with the experimental data.
基金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.
文摘We derive higher order magneto-hydrodynamic (MHD) equations from a microscopic picture using pro-jection and perturbation formalism. In an application to Hartmann flow we find velocity profiles flattening towards the center at the onset of turbulence in hydrodynamic limit. Comparison with the system under the effect of a uniform magnetic field yields difference in the onset of turbulence consistent with observations, showing that the presence of magnetic field inhibits onset of instability or turbulence. The laminar-turbulent transition is demonstrated in a phase transition plot of the development in time of the relative average velocities vs. Reynolds number showing a sharp increase of the relative average velocity at the transition point as determined by the critical Reynolds number.
基金supported by National Natural Science Foundation-the Science Foundation of China Academy of Engineering Physics (NSAF) under Grant No. 10576076the Major Projects of National Natural Science Foundation of China under Grant No. 10335010
文摘In this paper, the dynamics behaviors on fo-δ parameter surface is investigated for Gledzer-Ohkitani- Yamada model We indicate the type of intermittency chaos transitions is saddle node bifurcation. We plot phase diagram on fo-δ parameter surface, which is divided into periodic, quasi-periodic, and intermittent chaos areas. By means of varying Taylor-microscale Reynolds number, we calculate the extended self-similarity of velocity structure function.
文摘The flow field of liquid phase (water) of agitated extraction columns is simulated with the help of computational fluid dynamics (CFD). Four kinds of Reynolds-averaged turbulence models, i.e. the standard k-ε model, the RNG (renormalization group) k-s model, the realizable k-ε model and the Reynolds stress model, are compared in detail in order to judge which is the best model in terms of the accuracy, less CPU time and memory required. The performance of the realizable k-s model is obviously improved by reducing the model constant from C2 = 1.90 to C2 = 1.61. It is concluded that the improved realizable k-e model is the optimal model.