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 Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resour...The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.展开更多
This study demonstrates an active flow control for deflecting a direction of wake vortex structures behind a NACA0012 airfoil using an active morphing flap. Two-dimensional direct numerical simulations are performed f...This study demonstrates an active flow control for deflecting a direction of wake vortex structures behind a NACA0012 airfoil using an active morphing flap. Two-dimensional direct numerical simulations are performed for flows at the chord Reynolds number of 10,000, and the vortex pattern in the controlled and noncontrolled wakes as well as the effect of an actuation frequency on the control ability are rigorously investigated. It is found that there is an optimum actuation-frequency regime at around <em>F <sup>+</sup></em> = 2.00 which is normalized by the chord length and freestream velocity. The wake vortex pattern of the well-controlled case is classified as the 2P wake pattern according to the Williamson’s categorization [<a href="#ref1">1</a>] [<a href="#ref2">2</a>], where the forced oscillation frequency corresponds to the natural vortex shedding frequency without control. The present classification of wake vortex patterns and finding of the optimum frequency regime in the wake deflection control can lead to a more robust design suitable for vortex-induced-vibration (VIV) related engineering systems.展开更多
A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iter...A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iterative errors generated within one time step is constructed.With an estimate of temporal and spatial error using an embedded RungeKutta scheme and a higher order spatial discretization,an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously infuencing the total error of the discretization.The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection,steady-state fow past a fat plate,Taylor-Green vortex and turbulent fow over a circular cylinder at Re=3900.The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efciency is obtained for unsteady and steady,well-resolved and under-resolved simulations.展开更多
The use of the specified time interval numerical schemes has been popular in applying the method of characteristics to unsteady oped-channel flow problem. Studies and analyses of several variants of those schemes have...The use of the specified time interval numerical schemes has been popular in applying the method of characteristics to unsteady oped-channel flow problem. Studies and analyses of several variants of those schemes have led to the derivation of a complete set of timeline-interpolation scheme in applying the method of characteristics to Saint-Venant equations. A computer model based upon this scheme has ben successfully developed for predicting free surface unsteady flow, and has been verified by field experiment. Timeline interpolations, where the characteristic lines are projected back before the current time step,demonstrate much less damping than the corresponding spatial interpolation scheme at the same discretization.展开更多
Background Superconducting cavity is usually needed to be gradually cooled from room temperature to the superconducting temperature zone(4.2 K and below)in the testing and sophisticated operation process of supercondu...Background Superconducting cavity is usually needed to be gradually cooled from room temperature to the superconducting temperature zone(4.2 K and below)in the testing and sophisticated operation process of superconducting cavity.Purpose The purpose of this paper is to study the cooling law on the helium cooldown process for the 650 MHz two-cell superconducting cavity with the unsteady numerical simulation.Method A three-dimensional coupled heat-flow model of 650 MHz two-cell superconducting cavity was established.The unsteady numerical simulation of different inlet temperatures,flow rates and pressure conditions was carried out.The equiva-lent convective heat transfer coefficient and temperature distribution of 650 MHz two-cell superconducting cavity during cooldown process were obtained.The effects of cooling time and entrance parameters on the cooldown process were analyzed.Results The temperature distribution of the lower intersection lines has a large drop in the initial stage of cooldown process(120 s),while the temperature near the flanges at the both ends is still higher(remaining at the initial temperature of 300 K).With the passage of time,the temperature of the upper and lower intersection lines decreases.The maximum temperature difference on the lower intersections is within 2 K in the final stage of cooldown process(3600 s).The maximum temperature difference increases by 180%,and the difference between the maximum temperature and the minimum temperature(dT)at the end of a cooldown stage increases by 130%after 1 h,respectively,when the inlet temperature drops from 290 to 270 K(under the condition of the initial temperature of 300 K).Conclusions The maximum temperature difference and the dT at the end of a cooldown stage increase with the decrease in the inlet temperature.The maximum temperature difference increases with the increase in the inlet flow rate,while the dT at the end of a cooldown stage decreases with the increase in the inlet flow rate.The effect of changing the inlet flow rate on the cooling rate is not as obvious as changing the inlet temperature.Once there is a certain flow rate,the advantage of further increasing the flow rate to reduce the temperature of the superconducting cavity is not so great.展开更多
An unsteady numerical analysis has been conducted to study the strong interaction between impeller blade and volute tongue of a centrifugal pump. The 3-D-URANS equations were solved with the shear stress transport tur...An unsteady numerical analysis has been conducted to study the strong interaction between impeller blade and volute tongue of a centrifugal pump. The 3-D-URANS equations were solved with the shear stress transport turbulence model for a wide range of flow rates. These unsteady interactions are mostly related to the unsteady radial force due to an imbalance in the pressure field at the impeller periphery. This force represents dynamic load that are one of the most important sources of vibration and hydraulic noise. Based on this phenomenon, this work analyzes and gives a more realistic prediction of the pressure fluctuation and the radial force during steady and unsteady calculation by considering the effect of the change in the pump operating point. Actually, the pressure fluctuations in the impeller and the volute were recorded by mounting nine monitoring points on the impeller and volute casing. The results of the existing analysis has proven that the pressure fluctuation is periodic due to the relative position of impeller blade to volute tongue. The characteristics of the time domain and frequency domain of the pressure pulsation were analyzed under different coupling conditions. Fast Fourier transform was performed to obtain the spectra of pressure pulsation. Besides, the steady and unsteady forces were calculated around the impeller periphery to fully characterize the pump behavior. The obtained pump performance curves were numerically compared with the experimental ones, and the outcome have shown an acceptable agreement between both curves.展开更多
This paper presents a new numerical method to simulate the high velocity turbulent flow with free surface by solving two-dimensional incompressible unsteady Navier-Stokes Eqs. , together with the k-ε turbulence model...This paper presents a new numerical method to simulate the high velocity turbulent flow with free surface by solving two-dimensional incompressible unsteady Navier-Stokes Eqs. , together with the k-ε turbulence model. In order to treat the non-rectangular boundary (or curvilinear boundary), orthogonal boundary-fitted grid is used and the Navier-Stokes Eqs. and k-ε turbulence model are rewritten and discreted in orthogonal curvilinear coordinates. Meanwhile, gas-liquid two-field model theory is introduced to treat the free-surface problem.展开更多
As it is almost impossible to carry out the prototype hydro-turbine experiment be- fore the power plant is built up, rational prediction of pressure fluctuations in the prototype turbine is very important at the desig...As it is almost impossible to carry out the prototype hydro-turbine experiment be- fore the power plant is built up, rational prediction of pressure fluctuations in the prototype turbine is very important at the design stage. From this viewpoint, we at first treated the unsteady turbulent flow computation based on the modified RNG k-ε turbulence model through the whole flow passage to simulate the pressure fluctuation in a model turbine. Since fair agreement was recognized between the numerical results and the experimental data, this numerical method was applied to simulate the pressure fluctuations in the prototype turbine. From the comparison of them with the model turbine results, it is seen that their qualitative trend of pres- sure fluctuations are similar, but an appreciable difference is observed between the amplitudes of pressure fluctuation of the prototype turbine and that of the model turbine. Though the present findings may be explained by the effect of Reynolds number, further studies are expected for quantitative interpretation. We paid atten- tion to the interaction between the fluid and turbine structure. Adopting a weak fluid-solid coupling method, we studied the pressure fluctuation in the prototype turbine to clarify how the elastic behavior of runner blades influenced the charac- teristics of pressure fluctuation.展开更多
The fluid flow in tundish is a non-isothermal process and the temperature variation of stream from teeming ladle dominates the fluid flow and thermal distribution in tundish. A numerical model was established to inves...The fluid flow in tundish is a non-isothermal process and the temperature variation of stream from teeming ladle dominates the fluid flow and thermal distribution in tundish. A numerical model was established to investigate the effect of inlet cooling rate on fluid flow and temperature distribution in tundish based on a FTSC (Flexible Thin Slab Casting) tundish. The inlet cooling rate varies from 0. 5 to 0. 25 ~C/rain. Under the present calculation conditions, the following conclusions were made. When the stream temperature from teeming ladle drops seriously (for inlet cooling rate of 0.5℃/min), there is a "backward flow" at the coming end of casting. The horizontal flow along the free surface turns to flow along the bottom of tundish. The bottom flow shortens the fluid flow route in tundish and deteriorates the removal effect of nonmetallic inclusions from molten steel. Nevertheless, when the inlet cooling rate decreases to 0.25℃/min, the horizontal flow is sustained during the whole casting period. The present research provides theoretical directions for temperature control in teeming ladle and continuous casting tundish during production of advanced steels.展开更多
A turbine based combined cycle(TBCC)propulsion system uses a turbine-based engine to accelerate the vehicle from takeoff to the mode transition flight condition,at which point,the propulsion system performs a“mode tr...A turbine based combined cycle(TBCC)propulsion system uses a turbine-based engine to accelerate the vehicle from takeoff to the mode transition flight condition,at which point,the propulsion system performs a“mode transition”from the turbine to ramjet engine.Smooth inlet mode transition is accomplished when flow is diverted from one flowpath to the other,without experiencing unstart or buzz.The smooth inlet mode transition is a complex unsteady process and it is one of the enabling technologies for combined cycle engine to become a functional reality.In order to unveil the unsteady process of inlet mode transition,the research of over/under TBCC inlet mode transition was conducted through a numerical simulation.It shows that during the mode transition the terminal shock oscillates in the inlet.During the process of inlet mode transition mass flow rate and Mach number of turbojet flowpath reduce with oscillation.While in ramjet flowpath the flow field is non-uniform at the beginning of inlet mode transition.The speed of mode transition and the operation states of the turbojet and ramjet engines will affect the motion of terminal shock.The result obtained in present paper can help us realize the unsteady flow characteristic during the mode transition and provide some suggestions for TBCC inlet mode transition based on the smooth transition of thrust.展开更多
Large-eddy simulation(LES) was originally proposed for simulating atmospheric flows in the 1960 s and has become one of the most promising and successful methodology for simulating turbulent flows with the improveme...Large-eddy simulation(LES) was originally proposed for simulating atmospheric flows in the 1960 s and has become one of the most promising and successful methodology for simulating turbulent flows with the improvement of computing power. It is now feasible to simulate complex engineering flows using LES. However, apart from the computing power, significant challenges still remain for LES to reach a level of maturity that brings this approach to the mainstream of engineering and industrial computations. This paper will describe briefly LES formalism first, present a quick glance at its history, review its current state focusing mainly on its applications in transitional flows and gas turbine combustor flows, discuss some major modelling and numerical challenges/issues that we are facing now and in the near future, and finish with the concluding remarks.展开更多
The internal flow in an axial flow rotating machinery is affected by the rotating characteristics, often accompanied by a strong rotating separation under small flow conditions. At present, the very large eddy simulat...The internal flow in an axial flow rotating machinery is affected by the rotating characteristics, often accompanied by a strong rotating separation under small flow conditions. At present, the very large eddy simulation (VLES) model commonly used for the separation flow simulation still has certain limitations in simulating such rotating separation flow: (1) The Reynolds stress level is overestimated in the near-wall region. (2) The influence of the rotating effect cannot be effectively considered. The above two limitations affect the simulation accuracy of the VLES model for the rotating separation flow under small flow conditions in the axial flow rotating machinery. The objective of this paper is to provide a new hybrid unsteady Reynolds average Navier-Stokes/large eddy simulation (URANS/LES) model suitable for the simulation of the rotating separation flow in an axial flow rotating machinery. Compared with the original VLES method, the modifications are as follows: (1) A Reynolds stress damping function in the near-wall region is introduced to reduce the overestimation of the Reynolds stress caused by the near-wall Reynolds average Navier-Stokes (RANS) behavior of the VLES model. (2) A control function driven by the vortex is introduced to reflect the influence of the rotating effect. Three typical cases are used to verify the calculation accuracy of the modified model. It is shown that the modified model can capture more turbulent vortices based on the URANS grids, and the prediction accuracy of the rotating separation flow is effectively improved. Compared with the original VLES model, the modified model can accurately predict the head change in the hump region of the axial flow pump.展开更多
Casing treatment is an effective technique in extending stall margin of axial and centrifugal compressor.However,its impacts on the stall behaviour of mixed-flow compressor are still not completely understood until no...Casing treatment is an effective technique in extending stall margin of axial and centrifugal compressor.However,its impacts on the stall behaviour of mixed-flow compressor are still not completely understood until now.To conquer this issue,unsteady full-annulus simulations were conducted to investigate the stall mechanism of a mixed-flow compressor with and without axial slot casing treatment(ASCT).The circumferential propagating speed of spike inception resolved by the numerical approach is 87.1%of the shaft speed,which is identical to the test data.The numerical results confirmed that the mixed-flow compressor fell into rotating stall via spike-type with and without ASCT.The flow structure of the spike inception was investigated at 50%design rotational speed.Instantaneous static pressure traces extracted upstream of the leading edge had shown a classic spiky wave.Furthermore,it was found that with and without ASCT,the mixed-flow compressor stalled through spike with the characteristic of tip leakage spillage at leading edge and tip leakage backflow from trailing edge,which is different from a fraction of the centrifugal compressor.The resultant phenomenon provides conoborating evidence for that unlike in axial-flow compressor,the addition of ASCT does not change the stall characteristics of the mixed-flow compressor.The flow structure that induced spike inception with ASCT is similar to the case with smooth casing.In the throttling process,tip leakage flow vortex had been involved in the formation of tornado vortices,with one end at the suction side,and the other end at the casing-side.The low-pressure region relevant to the downward spike is caused by leading-edge separation vortex or tornado vortex.The high-pressure region relevant to the upward spike is induced by blockage from the passage vortex.These results not only can provide guidance for the design of casing treatment in mixed-flow compressor,but also can pave the way for the stall waring in the highly-loaded compressors of next-generation aeroengines.展开更多
The variable geometry supersonic inlet tends to decrease the throat area to reduce the Mach number upstream of the terminal shock,so as to reduce the flow loss.However,excessive Internal Contraction Ratio(ICR)exposes ...The variable geometry supersonic inlet tends to decrease the throat area to reduce the Mach number upstream of the terminal shock,so as to reduce the flow loss.However,excessive Internal Contraction Ratio(ICR)exposes the inlet to a greater risk of unstart,which inevitably results in a process of increasing the throat area to aid the inlet restart.In the above throat regulation process,the inlet undergoes the start,unstart,and restart states in turn.In order to reveal the flow structure and mechanism of this process,a two-dimensional unsteady numerical simulation combined with a dynamic mesh technique were employed.The shock-on-lip Mach number of the studied inlet is 4.0 and the flight angle of attack is+6°.Analysis was focused on the state with a freestream Mach number of 3.0.The results clearly show that the flow response hysteresis appears,and restart is only realized when the throat area is obviously increased as compared to that of unstart due to the historical unstart flow structure.In addition,three typical flow fields were analyzed,and it is found that the separation ahead of the inlet was the key factor affecting the hysteresis.Finally,unstart and restart boundaries of the inlet were discussed,and the factors influencing its deviation from the typical boundaries of dual-solution area were analyzed.The newly predicted unstart and restart boundaries are much closer to the CFD results.展开更多
文摘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
基金supported by National Natural Science Foundation of China (Grant No.11172040)
文摘The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.
文摘This study demonstrates an active flow control for deflecting a direction of wake vortex structures behind a NACA0012 airfoil using an active morphing flap. Two-dimensional direct numerical simulations are performed for flows at the chord Reynolds number of 10,000, and the vortex pattern in the controlled and noncontrolled wakes as well as the effect of an actuation frequency on the control ability are rigorously investigated. It is found that there is an optimum actuation-frequency regime at around <em>F <sup>+</sup></em> = 2.00 which is normalized by the chord length and freestream velocity. The wake vortex pattern of the well-controlled case is classified as the 2P wake pattern according to the Williamson’s categorization [<a href="#ref1">1</a>] [<a href="#ref2">2</a>], where the forced oscillation frequency corresponds to the natural vortex shedding frequency without control. The present classification of wake vortex patterns and finding of the optimum frequency regime in the wake deflection control can lead to a more robust design suitable for vortex-induced-vibration (VIV) related engineering systems.
基金Zhen-Guo Yan acknowledges supports from the National Natural Science Foundation of China(Grant no.11902344)National Numerical Windtunnel Project.The development of the implicit solver in Nektar++has been supported by EPSRC grant(EP/R029423/1)UK Turbulence Consortium grant(EP/R029326/1).
文摘A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iterative errors generated within one time step is constructed.With an estimate of temporal and spatial error using an embedded RungeKutta scheme and a higher order spatial discretization,an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously infuencing the total error of the discretization.The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection,steady-state fow past a fat plate,Taylor-Green vortex and turbulent fow over a circular cylinder at Re=3900.The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efciency is obtained for unsteady and steady,well-resolved and under-resolved simulations.
文摘The use of the specified time interval numerical schemes has been popular in applying the method of characteristics to unsteady oped-channel flow problem. Studies and analyses of several variants of those schemes have led to the derivation of a complete set of timeline-interpolation scheme in applying the method of characteristics to Saint-Venant equations. A computer model based upon this scheme has ben successfully developed for predicting free surface unsteady flow, and has been verified by field experiment. Timeline interpolations, where the characteristic lines are projected back before the current time step,demonstrate much less damping than the corresponding spatial interpolation scheme at the same discretization.
文摘Background Superconducting cavity is usually needed to be gradually cooled from room temperature to the superconducting temperature zone(4.2 K and below)in the testing and sophisticated operation process of superconducting cavity.Purpose The purpose of this paper is to study the cooling law on the helium cooldown process for the 650 MHz two-cell superconducting cavity with the unsteady numerical simulation.Method A three-dimensional coupled heat-flow model of 650 MHz two-cell superconducting cavity was established.The unsteady numerical simulation of different inlet temperatures,flow rates and pressure conditions was carried out.The equiva-lent convective heat transfer coefficient and temperature distribution of 650 MHz two-cell superconducting cavity during cooldown process were obtained.The effects of cooling time and entrance parameters on the cooldown process were analyzed.Results The temperature distribution of the lower intersection lines has a large drop in the initial stage of cooldown process(120 s),while the temperature near the flanges at the both ends is still higher(remaining at the initial temperature of 300 K).With the passage of time,the temperature of the upper and lower intersection lines decreases.The maximum temperature difference on the lower intersections is within 2 K in the final stage of cooldown process(3600 s).The maximum temperature difference increases by 180%,and the difference between the maximum temperature and the minimum temperature(dT)at the end of a cooldown stage increases by 130%after 1 h,respectively,when the inlet temperature drops from 290 to 270 K(under the condition of the initial temperature of 300 K).Conclusions The maximum temperature difference and the dT at the end of a cooldown stage increase with the decrease in the inlet temperature.The maximum temperature difference increases with the increase in the inlet flow rate,while the dT at the end of a cooldown stage decreases with the increase in the inlet flow rate.The effect of changing the inlet flow rate on the cooling rate is not as obvious as changing the inlet temperature.Once there is a certain flow rate,the advantage of further increasing the flow rate to reduce the temperature of the superconducting cavity is not so great.
文摘An unsteady numerical analysis has been conducted to study the strong interaction between impeller blade and volute tongue of a centrifugal pump. The 3-D-URANS equations were solved with the shear stress transport turbulence model for a wide range of flow rates. These unsteady interactions are mostly related to the unsteady radial force due to an imbalance in the pressure field at the impeller periphery. This force represents dynamic load that are one of the most important sources of vibration and hydraulic noise. Based on this phenomenon, this work analyzes and gives a more realistic prediction of the pressure fluctuation and the radial force during steady and unsteady calculation by considering the effect of the change in the pump operating point. Actually, the pressure fluctuations in the impeller and the volute were recorded by mounting nine monitoring points on the impeller and volute casing. The results of the existing analysis has proven that the pressure fluctuation is periodic due to the relative position of impeller blade to volute tongue. The characteristics of the time domain and frequency domain of the pressure pulsation were analyzed under different coupling conditions. Fast Fourier transform was performed to obtain the spectra of pressure pulsation. Besides, the steady and unsteady forces were calculated around the impeller periphery to fully characterize the pump behavior. The obtained pump performance curves were numerically compared with the experimental ones, and the outcome have shown an acceptable agreement between both curves.
文摘This paper presents a new numerical method to simulate the high velocity turbulent flow with free surface by solving two-dimensional incompressible unsteady Navier-Stokes Eqs. , together with the k-ε turbulence model. In order to treat the non-rectangular boundary (or curvilinear boundary), orthogonal boundary-fitted grid is used and the Navier-Stokes Eqs. and k-ε turbulence model are rewritten and discreted in orthogonal curvilinear coordinates. Meanwhile, gas-liquid two-field model theory is introduced to treat the free-surface problem.
基金Supported by the National Natural Science Foundation of China(Grant No.90410019)
文摘As it is almost impossible to carry out the prototype hydro-turbine experiment be- fore the power plant is built up, rational prediction of pressure fluctuations in the prototype turbine is very important at the design stage. From this viewpoint, we at first treated the unsteady turbulent flow computation based on the modified RNG k-ε turbulence model through the whole flow passage to simulate the pressure fluctuation in a model turbine. Since fair agreement was recognized between the numerical results and the experimental data, this numerical method was applied to simulate the pressure fluctuations in the prototype turbine. From the comparison of them with the model turbine results, it is seen that their qualitative trend of pres- sure fluctuations are similar, but an appreciable difference is observed between the amplitudes of pressure fluctuation of the prototype turbine and that of the model turbine. Though the present findings may be explained by the effect of Reynolds number, further studies are expected for quantitative interpretation. We paid atten- tion to the interaction between the fluid and turbine structure. Adopting a weak fluid-solid coupling method, we studied the pressure fluctuation in the prototype turbine to clarify how the elastic behavior of runner blades influenced the charac- teristics of pressure fluctuation.
文摘The fluid flow in tundish is a non-isothermal process and the temperature variation of stream from teeming ladle dominates the fluid flow and thermal distribution in tundish. A numerical model was established to investigate the effect of inlet cooling rate on fluid flow and temperature distribution in tundish based on a FTSC (Flexible Thin Slab Casting) tundish. The inlet cooling rate varies from 0. 5 to 0. 25 ~C/rain. Under the present calculation conditions, the following conclusions were made. When the stream temperature from teeming ladle drops seriously (for inlet cooling rate of 0.5℃/min), there is a "backward flow" at the coming end of casting. The horizontal flow along the free surface turns to flow along the bottom of tundish. The bottom flow shortens the fluid flow route in tundish and deteriorates the removal effect of nonmetallic inclusions from molten steel. Nevertheless, when the inlet cooling rate decreases to 0.25℃/min, the horizontal flow is sustained during the whole casting period. The present research provides theoretical directions for temperature control in teeming ladle and continuous casting tundish during production of advanced steels.
基金The authors gratefully acknowledge the financial support received from the Aviation Foundation Project(Grant No.2012ZB52031)the Fundamental Research Funds for the Central Universities(Grant No.NJ20140021)for this project.
文摘A turbine based combined cycle(TBCC)propulsion system uses a turbine-based engine to accelerate the vehicle from takeoff to the mode transition flight condition,at which point,the propulsion system performs a“mode transition”from the turbine to ramjet engine.Smooth inlet mode transition is accomplished when flow is diverted from one flowpath to the other,without experiencing unstart or buzz.The smooth inlet mode transition is a complex unsteady process and it is one of the enabling technologies for combined cycle engine to become a functional reality.In order to unveil the unsteady process of inlet mode transition,the research of over/under TBCC inlet mode transition was conducted through a numerical simulation.It shows that during the mode transition the terminal shock oscillates in the inlet.During the process of inlet mode transition mass flow rate and Mach number of turbojet flowpath reduce with oscillation.While in ramjet flowpath the flow field is non-uniform at the beginning of inlet mode transition.The speed of mode transition and the operation states of the turbojet and ramjet engines will affect the motion of terminal shock.The result obtained in present paper can help us realize the unsteady flow characteristic during the mode transition and provide some suggestions for TBCC inlet mode transition based on the smooth transition of thrust.
基金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)
文摘Large-eddy simulation(LES) was originally proposed for simulating atmospheric flows in the 1960 s and has become one of the most promising and successful methodology for simulating turbulent flows with the improvement of computing power. It is now feasible to simulate complex engineering flows using LES. However, apart from the computing power, significant challenges still remain for LES to reach a level of maturity that brings this approach to the mainstream of engineering and industrial computations. This paper will describe briefly LES formalism first, present a quick glance at its history, review its current state focusing mainly on its applications in transitional flows and gas turbine combustor flows, discuss some major modelling and numerical challenges/issues that we are facing now and in the near future, and finish with the concluding remarks.
基金the National Natural Science Foundation of China(Grant Nos.51836010,51779258).
文摘The internal flow in an axial flow rotating machinery is affected by the rotating characteristics, often accompanied by a strong rotating separation under small flow conditions. At present, the very large eddy simulation (VLES) model commonly used for the separation flow simulation still has certain limitations in simulating such rotating separation flow: (1) The Reynolds stress level is overestimated in the near-wall region. (2) The influence of the rotating effect cannot be effectively considered. The above two limitations affect the simulation accuracy of the VLES model for the rotating separation flow under small flow conditions in the axial flow rotating machinery. The objective of this paper is to provide a new hybrid unsteady Reynolds average Navier-Stokes/large eddy simulation (URANS/LES) model suitable for the simulation of the rotating separation flow in an axial flow rotating machinery. Compared with the original VLES method, the modifications are as follows: (1) A Reynolds stress damping function in the near-wall region is introduced to reduce the overestimation of the Reynolds stress caused by the near-wall Reynolds average Navier-Stokes (RANS) behavior of the VLES model. (2) A control function driven by the vortex is introduced to reflect the influence of the rotating effect. Three typical cases are used to verify the calculation accuracy of the modified model. It is shown that the modified model can capture more turbulent vortices based on the URANS grids, and the prediction accuracy of the rotating separation flow is effectively improved. Compared with the original VLES model, the modified model can accurately predict the head change in the hump region of the axial flow pump.
基金The authors are grateful for the support of National Natural Science Foundation of China for the grant with Project No.51636001,No.51790510,and No.51922098the National Science and Technology Major Project(2017-11-0004-0017,2017-H-005-0018)The authors would also appreciate the kind help of Prof.Seume,and the German Research Foundation(DFG)for supporting the experiments as part of the Collaborative Research Centre 880(Sonderforschungsbereich SFB 880).
文摘Casing treatment is an effective technique in extending stall margin of axial and centrifugal compressor.However,its impacts on the stall behaviour of mixed-flow compressor are still not completely understood until now.To conquer this issue,unsteady full-annulus simulations were conducted to investigate the stall mechanism of a mixed-flow compressor with and without axial slot casing treatment(ASCT).The circumferential propagating speed of spike inception resolved by the numerical approach is 87.1%of the shaft speed,which is identical to the test data.The numerical results confirmed that the mixed-flow compressor fell into rotating stall via spike-type with and without ASCT.The flow structure of the spike inception was investigated at 50%design rotational speed.Instantaneous static pressure traces extracted upstream of the leading edge had shown a classic spiky wave.Furthermore,it was found that with and without ASCT,the mixed-flow compressor stalled through spike with the characteristic of tip leakage spillage at leading edge and tip leakage backflow from trailing edge,which is different from a fraction of the centrifugal compressor.The resultant phenomenon provides conoborating evidence for that unlike in axial-flow compressor,the addition of ASCT does not change the stall characteristics of the mixed-flow compressor.The flow structure that induced spike inception with ASCT is similar to the case with smooth casing.In the throttling process,tip leakage flow vortex had been involved in the formation of tornado vortices,with one end at the suction side,and the other end at the casing-side.The low-pressure region relevant to the downward spike is caused by leading-edge separation vortex or tornado vortex.The high-pressure region relevant to the upward spike is induced by blockage from the passage vortex.These results not only can provide guidance for the design of casing treatment in mixed-flow compressor,but also can pave the way for the stall waring in the highly-loaded compressors of next-generation aeroengines.
基金co-supported by the National Natural Science Foundation of China(Nos.U20A2070,12025202,11772156,51806102,and 51906104)。
文摘The variable geometry supersonic inlet tends to decrease the throat area to reduce the Mach number upstream of the terminal shock,so as to reduce the flow loss.However,excessive Internal Contraction Ratio(ICR)exposes the inlet to a greater risk of unstart,which inevitably results in a process of increasing the throat area to aid the inlet restart.In the above throat regulation process,the inlet undergoes the start,unstart,and restart states in turn.In order to reveal the flow structure and mechanism of this process,a two-dimensional unsteady numerical simulation combined with a dynamic mesh technique were employed.The shock-on-lip Mach number of the studied inlet is 4.0 and the flight angle of attack is+6°.Analysis was focused on the state with a freestream Mach number of 3.0.The results clearly show that the flow response hysteresis appears,and restart is only realized when the throat area is obviously increased as compared to that of unstart due to the historical unstart flow structure.In addition,three typical flow fields were analyzed,and it is found that the separation ahead of the inlet was the key factor affecting the hysteresis.Finally,unstart and restart boundaries of the inlet were discussed,and the factors influencing its deviation from the typical boundaries of dual-solution area were analyzed.The newly predicted unstart and restart boundaries are much closer to the CFD results.