Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
To study the effects of factors such as electromagnetic field and electron temperature on the ion extraction characteristics (extraction time, collision loss ratio) in atomic vapor laser isotope separation (AVLIS),...To study the effects of factors such as electromagnetic field and electron temperature on the ion extraction characteristics (extraction time, collision loss ratio) in atomic vapor laser isotope separation (AVLIS), 2 D electron equilibrium model was used to analyze and simulate the ion extraction process. The collisions between particles mainly considered charge exchanges between isotopic ions and atoms, which were treated by using cross section. The results show that the electric field and electron temperature have significant effects on the ion extraction characteristics: the stronger the electric field is, the higher the electron temperature is; the shorter the extraction time is, the less the collision loss ratio is, and moderate constant magnetic field has few effects on them. Key words atomic vapor laser isotope separation (AVLIS); ion extraction; charge exchange cross section; collision loss ratio展开更多
To improve the performance of the positive displacement blower, it is imperative to understand the detailed internal flow characteristics or enable a visualization of flow status. However, the existing two-dimensional...To improve the performance of the positive displacement blower, it is imperative to understand the detailed internal flow characteristics or enable a visualization of flow status. However, the existing two-dimensional unsteady, three-dimensional steady or quasi-unsteady numerical simulation and theoretical analysis cannot provide the detailed flow information, which is unfavorable to improve the performance of positive displacement blower. Therefore, the unsteady flow characteristics in a three-lobe positive displacement blower are numerically investigated by solving the three-dimensional, unsteady, compressible Navier-Stokes equations coupled with RNG k-e turbulent model. In the numerical simulation, the dynamic mesh technique and overset mesh updating method are adopted. Due to the air being compressed in the process of the rotors rotating, the variation of the temperature field in the positive displacement blower is considered. By comparing the experimental measurements and the numerical results on the variation of flow rate with the outlet pressure, the maximum relative error of the flow rate is less than 2.15% even at the maximum outlet pressure condition, which means that the calculation model and numerical computational method used are effective. The numerical results show that in the intake region, the fluctuations of the inlet flow are greatly affected by the direction of the velocity vectors. In the exhaust region, the temperature changes significantly, which leads to the increase of the airflow pulsation. Through analysis on the velocity, pressure and temperature fields obtained from the numerical simulations, three-dimensional unsteady flow characteristics in the positive displacement blower are revealed. The studied results will provide useful reference for improving the performance and empirical correction in the design of the positive displacement blower.展开更多
The application of small diameter arterial grafts is limited due to the fact of relatively poor long-time patency which is caused by thrombosis formation in the short term and intimal hyperplasia(IH) in the medium and...The application of small diameter arterial grafts is limited due to the fact of relatively poor long-time patency which is caused by thrombosis formation in the short term and intimal hyperplasia(IH) in the medium and long term.Thrombosis,obstructing the flow of blood展开更多
The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for...The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.展开更多
A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusio...A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.展开更多
The application of unsteady simulation to lightening urban storm waterlogged disaster involves six blocks: the main program, the database of drainage networks, the abstractions and losses block, the calculating inflow...The application of unsteady simulation to lightening urban storm waterlogged disaster involves six blocks: the main program, the database of drainage networks, the abstractions and losses block, the calculating inflow of sewer block, the drainage networks routing block and the waterlogged areas calculating block. This model can be used for simulating sewer surcharge and calculating area and water depth of waterlogged areas under storm conditions. The best design of rebuilding storm sewer system or combined sewer system can be provided by simulating existing state of built drainage pipeline. The model was used for rebuilding Shenyang drainage works last year, and the good results were obtained.展开更多
Numerical studies of the flow past a circular cylinder at Reynolds number 1.4 × 105 and NACA0021 airfoil at the angle of attack 60° have been carried out by scale-adaptive simulation (SAS) and detached edd...Numerical studies of the flow past a circular cylinder at Reynolds number 1.4 × 105 and NACA0021 airfoil at the angle of attack 60° have been carried out by scale-adaptive simulation (SAS) and detached eddy simu- lation (DES), in comparison with the existing experimental data. The new version of the model developed by Egorov and Menter is assessed, and advantages and disadvantages of the SAS simulation are analyzed in detail to provide guidance for industrial application in the future. Moreover, the mechanism of the scale-adaptive characteristics in separated regions is discussed, which is obscure in previous analyses. It is con- cluded that: the mean flow properties satisfactorily agree with the experimental results for the SAS simulation, although the prediction of the second order turbulent statistics in the near wake region is just reasonable. The SAS model can produce a larger magnitude of the turbulent kinetic energy in the recir- culation bubble, and, consequently, a smaller recirculation region and a more rapid recovery of the mean velocity out- side the recirculation region than the DES approach with the same grid resolution. The vortex shedding is slightly less irregular with the SAS model than with the DES approach, probably due to the higher dissipation of the SAS simulation under the condition of the coarse mesh.展开更多
The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller (HALLER, G. Exact theor...The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller (HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 (2004)). By analyzing the distribution of the finite-time Lyapunov exponent (FTLE) along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.展开更多
In this paper,a numerical simulation method is used to calculate a 1.5-stage axial transonic compressor to explore its unsteady flow mechanism.The performance curve is compared with the experimental data to verify the...In this paper,a numerical simulation method is used to calculate a 1.5-stage axial transonic compressor to explore its unsteady flow mechanism.The performance curve is compared with the experimental data to verify the calculation method with a high numerical accuracy,which shows that the unsteady calculation has good reliability.According to the analysis of the data from the monitoring points under the near-stall condition,the unsteady disturbances originate from the tip region of blade and perform the strongest at the blade pressure surface with a broadband characteristic.Further analysis is conducted by combining with the characteristics of the transient flow field at the tip of blade.The results show that the unsteady pressure fluctuations are caused by the migration of the new vortex cores.These new vortex cores are generated by the breakdown of leakage vortex in the downstream,which is induced by the leakage vortex and shock wave interference.Moreover,the relationship between the unsteady flow characteristics and the working conditions is also studied.The leakage vortex intensity and the shock wave strength gradually increase with the decrease of flow rate.When the combination of the leakage vortex intensity and shock wave strength reaches the first threshold,a single frequency of unsteady disturbances appears at the blade tip.When the combination of the leakage vortex intensity and shock wave strength reaches the second threshold,the frequency of unsteady disturbances changes to a broadband.展开更多
Stall phenomena increase the complexity of the internal flow in centrifugal pump impellers.In order to tackle this problem,in the present work,a large eddy simulation(LES)approach is applied to determine the character...Stall phenomena increase the complexity of the internal flow in centrifugal pump impellers.In order to tackle this problem,in the present work,a large eddy simulation(LES)approach is applied to determine the characteristics of these unstable flows.Moreover,a vorticity identification method is used to characterize quantitatively the vortex position inside the impeller and its influencing area.By comparing the outcomes of the numerical simulations and experimental results provided by a Particle Image Velocimetry(PIV)technique,it is shown that an apparent“alternating stall”phenomenon exists inside the impeller when relatively small flow rate conditions are considered.The stall is generated near the suction side of the blade inlet,grows towards the high-pressure side of the blade in the circumferential direction,and gradually attenuates.As the flow rate decreases,the number of stalls remains unchanged,while the related influencing area and strength gradually increase and the circumferential velocity increases.展开更多
For Francis hydraulic turbines, unsteady flow caused by vortex ropes in the draft tube leads to a problem of stability in operation. The unsteady flow field of a model Francis hydraulic turbine was simulated under par...For Francis hydraulic turbines, unsteady flow caused by vortex ropes in the draft tube leads to a problem of stability in operation. The unsteady flow field of a model Francis hydraulic turbine was simulated under part-load operation. A sliding mesh model was used to calculate a time-accurate solution for the strong rotor-stator interactions between the runner and guide vanes, and the draft tube. Based on three-dimensional incompressible Reynolds averaged Navier-Stokes equations and on a renormalization group k-?turbulence model, spatial discretization was obtained by using the finite volume method with unstructured grid elements, and a second order fully implicit scheme was applied for time. Pressure fluctuations in the draft tube were recorded and analyzed via a fast Fourier transform calculation. The results were compared with the experimental data, and show that the vortex rope in the draft tube and the induced pressure fluctuations are well simulated.展开更多
Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks th...Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.展开更多
This paper describes the numerical simulation of unsteady flows due to incoming wakes and/or varying back pressure,The solution method is based upon the one-step finite-volume TVD Lax-Wendroff scheme.Dual time-step ap...This paper describes the numerical simulation of unsteady flows due to incoming wakes and/or varying back pressure,The solution method is based upon the one-step finite-volume TVD Lax-Wendroff scheme.Dual time-step approach and multigrid algorithm are adopted to improve the computational efficiency of the baseline scheme.Numerical results for the transonic unsteady flow in a channel bump and the unsteady flow in a flat plate cascade and the VKI cascade are presented.展开更多
The vortex-in-cell method in the discrete vortex simulation for the separated flow around a bluff body is discussed,Some improvements are made.The separated flow around a circular cylinder in oscillatory flow is inves...The vortex-in-cell method in the discrete vortex simulation for the separated flow around a bluff body is discussed,Some improvements are made.The separated flow around a circular cylinder in oscillatory flow is investigated.展开更多
A CFD based numerical simulation of flow velocity of hydrocyclone was conducted with different structural and operational parameters to investigate its distribution characteristics and influencing mechanism. The resul...A CFD based numerical simulation of flow velocity of hydrocyclone was conducted with different structural and operational parameters to investigate its distribution characteristics and influencing mechanism. The results show there exist several unsymmetrical envelopes of equal vertical velocities in both upward inner flows and downward outer flows in the hydrocyclone, and the cone angle and apex diameter have remarkable influence on the vertical location of the cone bottom of the envelope of zero vertical velocity. It is also found that the tangential velocity isolines exist in the horizontal planes located in the effective separation region of hydrocyclone. The increase of feed pressure has almost no effect on the distribution characteristics of both vertical velocity and tangential velocity in hydrocyclone, but the magnitude and gradient of tangential velocity are increased obviously to make the motion velocity of high density particles to the wall increased and to make the cyclonic separation effect improved.展开更多
The operation methods of channel and the speed of gate regulation have great influence on the transformation of flow in water conveyance channels. Based on characteristics method, a 1-D unsteady flow numerical model f...The operation methods of channel and the speed of gate regulation have great influence on the transformation of flow in water conveyance channels. Based on characteristics method, a 1-D unsteady flow numerical model for gate regulation was established in this study. The process of water flow was simulated under different boundary conditions. The influence of gate regulation speed and channel operation methods on flow transition process was analyzed. The numerical results show that under the same conditions, with increasing regulation speed of the gate, the change rates of discharge and water level increase, while the response time of channel becomes shorter, and ultimately the discharge and water level will transit to the same equilibrium states. Moreover, the flow is easier to reach stable state, if the water level in front of the sluice is kept constant, instead of,behind the sluice. This study will be important to the scheme design of automatic operation control in water conveyance channels.展开更多
Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studie...Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.展开更多
The unsteady flow around a ship advancing at a constant speed was studied in time domain by using a higher-order boundary element method. The Rankine source was employed. The ship hull and free surface were discretize...The unsteady flow around a ship advancing at a constant speed was studied in time domain by using a higher-order boundary element method. The Rankine source was employed. The ship hull and free surface were discretized into 8-node isoparametric elements. A modified Green formula was used to overcome the numerical difficulty arising from singular or quasi-singular integral equation. In the case of τ>0.25, a rigid wall condition was applied to the upstream and a numerical beach was applied to the downstream to absorb the outgoing waves. The numerical results are satisfactory compared with the experimental results.展开更多
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
文摘To study the effects of factors such as electromagnetic field and electron temperature on the ion extraction characteristics (extraction time, collision loss ratio) in atomic vapor laser isotope separation (AVLIS), 2 D electron equilibrium model was used to analyze and simulate the ion extraction process. The collisions between particles mainly considered charge exchanges between isotopic ions and atoms, which were treated by using cross section. The results show that the electric field and electron temperature have significant effects on the ion extraction characteristics: the stronger the electric field is, the higher the electron temperature is; the shorter the extraction time is, the less the collision loss ratio is, and moderate constant magnetic field has few effects on them. Key words atomic vapor laser isotope separation (AVLIS); ion extraction; charge exchange cross section; collision loss ratio
基金Supported by Fundamental Research Funds for the Central Universities,China(Grant No.xjj20100073)Science and Technology Innovation Project of Shaanxi Province of China(Grant No.2011KTCL01-04)
文摘To improve the performance of the positive displacement blower, it is imperative to understand the detailed internal flow characteristics or enable a visualization of flow status. However, the existing two-dimensional unsteady, three-dimensional steady or quasi-unsteady numerical simulation and theoretical analysis cannot provide the detailed flow information, which is unfavorable to improve the performance of positive displacement blower. Therefore, the unsteady flow characteristics in a three-lobe positive displacement blower are numerically investigated by solving the three-dimensional, unsteady, compressible Navier-Stokes equations coupled with RNG k-e turbulent model. In the numerical simulation, the dynamic mesh technique and overset mesh updating method are adopted. Due to the air being compressed in the process of the rotors rotating, the variation of the temperature field in the positive displacement blower is considered. By comparing the experimental measurements and the numerical results on the variation of flow rate with the outlet pressure, the maximum relative error of the flow rate is less than 2.15% even at the maximum outlet pressure condition, which means that the calculation model and numerical computational method used are effective. The numerical results show that in the intake region, the fluctuations of the inlet flow are greatly affected by the direction of the velocity vectors. In the exhaust region, the temperature changes significantly, which leads to the increase of the airflow pulsation. Through analysis on the velocity, pressure and temperature fields obtained from the numerical simulations, three-dimensional unsteady flow characteristics in the positive displacement blower are revealed. The studied results will provide useful reference for improving the performance and empirical correction in the design of the positive displacement blower.
文摘The application of small diameter arterial grafts is limited due to the fact of relatively poor long-time patency which is caused by thrombosis formation in the short term and intimal hyperplasia(IH) in the medium and long term.Thrombosis,obstructing the flow of blood
基金The project is supported by the National Nature Science Foundation of China(NNSFC)
文摘The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.
文摘A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.
文摘The application of unsteady simulation to lightening urban storm waterlogged disaster involves six blocks: the main program, the database of drainage networks, the abstractions and losses block, the calculating inflow of sewer block, the drainage networks routing block and the waterlogged areas calculating block. This model can be used for simulating sewer surcharge and calculating area and water depth of waterlogged areas under storm conditions. The best design of rebuilding storm sewer system or combined sewer system can be provided by simulating existing state of built drainage pipeline. The model was used for rebuilding Shenyang drainage works last year, and the good results were obtained.
基金the National Basic Research Program of China ("973" Project) (Grant No. 2009CB724104)
文摘Numerical studies of the flow past a circular cylinder at Reynolds number 1.4 × 105 and NACA0021 airfoil at the angle of attack 60° have been carried out by scale-adaptive simulation (SAS) and detached eddy simu- lation (DES), in comparison with the existing experimental data. The new version of the model developed by Egorov and Menter is assessed, and advantages and disadvantages of the SAS simulation are analyzed in detail to provide guidance for industrial application in the future. Moreover, the mechanism of the scale-adaptive characteristics in separated regions is discussed, which is obscure in previous analyses. It is con- cluded that: the mean flow properties satisfactorily agree with the experimental results for the SAS simulation, although the prediction of the second order turbulent statistics in the near wake region is just reasonable. The SAS model can produce a larger magnitude of the turbulent kinetic energy in the recir- culation bubble, and, consequently, a smaller recirculation region and a more rapid recovery of the mean velocity out- side the recirculation region than the DES approach with the same grid resolution. The vortex shedding is slightly less irregular with the SAS model than with the DES approach, probably due to the higher dissipation of the SAS simulation under the condition of the coarse mesh.
基金supported by the National Natural Science Foundation of China(Nos.11372340 and 11732016)
文摘The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller (HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 (2004)). By analyzing the distribution of the finite-time Lyapunov exponent (FTLE) along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.
基金the support of the grants of Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA29050500)。
文摘In this paper,a numerical simulation method is used to calculate a 1.5-stage axial transonic compressor to explore its unsteady flow mechanism.The performance curve is compared with the experimental data to verify the calculation method with a high numerical accuracy,which shows that the unsteady calculation has good reliability.According to the analysis of the data from the monitoring points under the near-stall condition,the unsteady disturbances originate from the tip region of blade and perform the strongest at the blade pressure surface with a broadband characteristic.Further analysis is conducted by combining with the characteristics of the transient flow field at the tip of blade.The results show that the unsteady pressure fluctuations are caused by the migration of the new vortex cores.These new vortex cores are generated by the breakdown of leakage vortex in the downstream,which is induced by the leakage vortex and shock wave interference.Moreover,the relationship between the unsteady flow characteristics and the working conditions is also studied.The leakage vortex intensity and the shock wave strength gradually increase with the decrease of flow rate.When the combination of the leakage vortex intensity and shock wave strength reaches the first threshold,a single frequency of unsteady disturbances appears at the blade tip.When the combination of the leakage vortex intensity and shock wave strength reaches the second threshold,the frequency of unsteady disturbances changes to a broadband.
基金This research was funded by the Zhejiang Provincial Natural Science Foundation of China(Grant Nos.LGG21E090002,LY21E060004,LGG21E090003)National Natural Science Foundation of China(Grant No.51779226)the China Postdoctoral Science Foundation(Grant No.2021M691383).
文摘Stall phenomena increase the complexity of the internal flow in centrifugal pump impellers.In order to tackle this problem,in the present work,a large eddy simulation(LES)approach is applied to determine the characteristics of these unstable flows.Moreover,a vorticity identification method is used to characterize quantitatively the vortex position inside the impeller and its influencing area.By comparing the outcomes of the numerical simulations and experimental results provided by a Particle Image Velocimetry(PIV)technique,it is shown that an apparent“alternating stall”phenomenon exists inside the impeller when relatively small flow rate conditions are considered.The stall is generated near the suction side of the blade inlet,grows towards the high-pressure side of the blade in the circumferential direction,and gradually attenuates.As the flow rate decreases,the number of stalls remains unchanged,while the related influencing area and strength gradually increase and the circumferential velocity increases.
文摘For Francis hydraulic turbines, unsteady flow caused by vortex ropes in the draft tube leads to a problem of stability in operation. The unsteady flow field of a model Francis hydraulic turbine was simulated under part-load operation. A sliding mesh model was used to calculate a time-accurate solution for the strong rotor-stator interactions between the runner and guide vanes, and the draft tube. Based on three-dimensional incompressible Reynolds averaged Navier-Stokes equations and on a renormalization group k-?turbulence model, spatial discretization was obtained by using the finite volume method with unstructured grid elements, and a second order fully implicit scheme was applied for time. Pressure fluctuations in the draft tube were recorded and analyzed via a fast Fourier transform calculation. The results were compared with the experimental data, and show that the vortex rope in the draft tube and the induced pressure fluctuations are well simulated.
基金supported by National Natural Science Foundation(52204050)Sichuan Science and Technology Program(2021ZHCG0013,22ZDYF3009)。
文摘Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.
文摘This paper describes the numerical simulation of unsteady flows due to incoming wakes and/or varying back pressure,The solution method is based upon the one-step finite-volume TVD Lax-Wendroff scheme.Dual time-step approach and multigrid algorithm are adopted to improve the computational efficiency of the baseline scheme.Numerical results for the transonic unsteady flow in a channel bump and the unsteady flow in a flat plate cascade and the VKI cascade are presented.
文摘The vortex-in-cell method in the discrete vortex simulation for the separated flow around a bluff body is discussed,Some improvements are made.The separated flow around a circular cylinder in oscillatory flow is investigated.
基金Project (50974033) supported by the National Natural Science Foundation of ChinaProject (N100301002) supported by the Fundamental Research Funds for the Universities, China
文摘A CFD based numerical simulation of flow velocity of hydrocyclone was conducted with different structural and operational parameters to investigate its distribution characteristics and influencing mechanism. The results show there exist several unsymmetrical envelopes of equal vertical velocities in both upward inner flows and downward outer flows in the hydrocyclone, and the cone angle and apex diameter have remarkable influence on the vertical location of the cone bottom of the envelope of zero vertical velocity. It is also found that the tangential velocity isolines exist in the horizontal planes located in the effective separation region of hydrocyclone. The increase of feed pressure has almost no effect on the distribution characteristics of both vertical velocity and tangential velocity in hydrocyclone, but the magnitude and gradient of tangential velocity are increased obviously to make the motion velocity of high density particles to the wall increased and to make the cyclonic separation effect improved.
基金supported by the National Natural SciencFoundation of China (Grant No. 59879016)the National BasiResearch Program of China (973 Program, Grant No2003CB415200)
文摘The operation methods of channel and the speed of gate regulation have great influence on the transformation of flow in water conveyance channels. Based on characteristics method, a 1-D unsteady flow numerical model for gate regulation was established in this study. The process of water flow was simulated under different boundary conditions. The influence of gate regulation speed and channel operation methods on flow transition process was analyzed. The numerical results show that under the same conditions, with increasing regulation speed of the gate, the change rates of discharge and water level increase, while the response time of channel becomes shorter, and ultimately the discharge and water level will transit to the same equilibrium states. Moreover, the flow is easier to reach stable state, if the water level in front of the sluice is kept constant, instead of,behind the sluice. This study will be important to the scheme design of automatic operation control in water conveyance channels.
基金National Natural Science Foundation of China (51176013)Ph.D. Programs Foundation of Ministry of Education of China (20091101110014)National High-tech Research and Development Program of China (2007AA050502)
文摘Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.
基金the National Natural Science Foundation of China and the Chinese State Edu-cation Commission Grant.
文摘The unsteady flow around a ship advancing at a constant speed was studied in time domain by using a higher-order boundary element method. The Rankine source was employed. The ship hull and free surface were discretized into 8-node isoparametric elements. A modified Green formula was used to overcome the numerical difficulty arising from singular or quasi-singular integral equation. In the case of τ>0.25, a rigid wall condition was applied to the upstream and a numerical beach was applied to the downstream to absorb the outgoing waves. The numerical results are satisfactory compared with the experimental results.