This study aims to analyze mixed convection in a square cavity with two moving vertical walls by finite volume method.The cavity filled with Non-Newtonian fluid of Bingham model is heated from below and cooled by the ...This study aims to analyze mixed convection in a square cavity with two moving vertical walls by finite volume method.The cavity filled with Non-Newtonian fluid of Bingham model is heated from below and cooled by the other walls.This study has been conducted for certain parameters of Reynolds number(Re=1-100),Richardson number(Ri=1-20),Prandtl number(Pr=1-500),and Bingham number has been studied from 0 to 10.The results indicate that the increase in yield stress drops the heat transfer and the flow become flatter,while increasing Reynolds number augments it.The convective transport is dominant when increasing Richardson number which leads to enhance heat transfer in the cavity for both Newtonian and Non-Newtonian fluid.A correlation of Nusselt number is given in function of different parameters.展开更多
The increase of wave energy in electricity production is an objective shared by many countries to meet growing demand and global warming. To analyze devices capable of converting the energy of sea waves into electrica...The increase of wave energy in electricity production is an objective shared by many countries to meet growing demand and global warming. To analyze devices capable of converting the energy of sea waves into electrical energy, it is important to master the various theories of gravity waves and generation. We will in our work consider a numerical waves tank for an amplitude A=0.5, a wavelength λ=0.25 , an average height H<sub>e</sub>=10 and a Froude number fixed at 1 × 10<sup>5</sup>. Numerical wave channel analysis is used to reproduce the natural phenomenon of wave propagation in an experimental model. Wave makers are usually used to generate waves in the channel. In theory, the influence of an incident wave can be considered, as in the case of our study. In this study, the evolution of the hydrodynamic parameters and the energy transported in one wavelength can be determined by calculation. A change of variable will be done in this work to facilitate the writing of the boundary conditions at the free surface and at the bottom. The nonlinear Stokes theory will be studied in this case in order to provide hydrodynamic solutions through the Navier-Stokes equations to finally deduce the energetic results. To do this, the finite difference method will be used for the hydrodynamic results such as the velocity potential and the free surface elevation and the trapezium method of Newton for the energetic results. Thus, we will determine the energetic potential according to the decrease in the slope of the tank. To do this, we will take as values of beta representing the inverse of the slope of the tank, β=100, β=105, β=110 and β=105. .展开更多
It is widely accepted that in a turbulent boundary layer (TBL) with adverse pressure gradient (APG) an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this ou...It is widely accepted that in a turbulent boundary layer (TBL) with adverse pressure gradient (APG) an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this outer peak is not clearly understood. In this paper, the effect of APG is analysed using the numerical and experimental results in the literature. Because the effect of upstream flow is inherent in the TBL, we first analyse this effect in TBLs with zero pressure gradient on flat plates. Under the individual effect of upstream flow, an outer peak already appears in the profile of streamwise Reynolds stress when the TBL continues developing in the streamwise direction. The APG accelerates the appearance of the outer peak, instead of being a trigger.展开更多
Natural convection heat transfer in open or closed cavities takes place in different engineering areas. The hemispherical cavity is a part of basic geometries although it is not widely studied. The present paper repor...Natural convection heat transfer in open or closed cavities takes place in different engineering areas. The hemispherical cavity is a part of basic geometries although it is not widely studied. The present paper reports the numerical study of natural convection in a closed hemispherical annulus delimited by two vertically eccentric hemispheres filled with Newtonian fluid (air in this case with <em>Pr</em> = 0.7) is conducted. The inner hemisphere is heated by a heat flux of constant density and the outer one is maintained isothermal. Based on the Boussinesq assumptions, the governing equations are numerically studied using unsteady natural convection formulated with vorticity and stream-function variables. These equations are written by using bispherical coordinates system and solved by using a finite difference method. The effect of the control parameters such as the Rayleigh number (<span style="white-space:nowrap;">10<sup>3</sup> ≤ <em>Ra</em> ≤ 10<sup>6</sup></span>) or the eccentricity (<em>e</em> = ±0.2, ±0.5, 0) in the dynamic and thermal behaviours of the fluid is investigated.展开更多
Contrary to natural cavitation,ventilated cavitation is controllable and is not harmful.It is particularly used to reduce the drag of the hydraulic vehicles.The ventilated cavitation is characterized by various gas re...Contrary to natural cavitation,ventilated cavitation is controllable and is not harmful.It is particularly used to reduce the drag of the hydraulic vehicles.The ventilated cavitation is characterized by various gas regimes.The mechanisms of ventilated cavitation are investigated in the present work with CFD based on a 2D solver.The attention is especially focused on the transition between the reentrant jet and twin vortex regimes.The results confirmthat the product of ventilated cavitation number and Froude number is lower than 1(σcFr<1)in the twin vortex regime,while it is higher than 1(σcFr>1)in the reentrant jet regime,as reported in the literature.Further analysis shows that ventilated cavitation is significantly influenced by the natural cavitation number.展开更多
Abstract: This work studies the active control of chemical oscillations governed by a forced modified Van der Pol-Duffing oscillator. We considered the dynamics of nonlinear chemical systems subjected to an external s...Abstract: This work studies the active control of chemical oscillations governed by a forced modified Van der Pol-Duffing oscillator. We considered the dynamics of nonlinear chemical systems subjected to an external sinusoidal excitation. The approximative solution to the first order of the modified Van der Pol-Duffing oscillator is found using the Lindstedt’s perturbation method. The harmonic balance method is used to find the amplitudes of the oscillatory states of the system under control. The effects of the constraint parameter and the control parameter of the model on the amplitude of oscillations are presented. The effects of the active control on the behaviors of the model are analyzed and it appears that with the appropriate selection of the coupling parameter, the chaotic behavior of the model has given way to periodic movements. Numerical simulations are used to validate and complete the analytical results obtained.展开更多
In our study, we investigate the differences between the combustion of different hydrocarbon fuels CH4, C3H8, C4H10. A numerical simulation of an impinging jet diffusion flames is used. The jet injector has a 10 mm in...In our study, we investigate the differences between the combustion of different hydrocarbon fuels CH4, C3H8, C4H10. A numerical simulation of an impinging jet diffusion flames is used. The jet injector has a 10 mm in diameter and the distance between the jet flame and the vertical wall is 2 time half diameter. The fuel jet velocity was fixed for 11.8 m/s, corresponding to a Reynolds number of 6881. The flame characteristics varied from hydrocarbon to another for the same Reynolds number. The combustion products of CO, CO2, NO, OH, are depending on the methane and propane and butane flames for the same conditions. The temperature of the flame was varied from hydrocarbon to another the same as for the chemical species production rate. The concentration of the thermal and prompt NO pollutant depends on the temperature flow field and on the thermochemical characteristics of the hydrocarbon fuels.展开更多
The subject of natural convection heat transfer is motivated by a wide range of applications in engineering technology. The hemispherical cavity is a part of basic geometries although it is not widely studied. The eff...The subject of natural convection heat transfer is motivated by a wide range of applications in engineering technology. The hemispherical cavity is a part of basic geometries although it is not widely studied. The effect of inclinaison on natural convection fluid motions in the gap between two eccentric hemispheres is numerically studied. The inner hemisphere is subjected to a heat flux of a constant density and the outer one is maintened isothermal. The walls separating the two hemispheres are thermally adiabatic. Equations are formulated with vorticity and stream-functions variables. It is also assumed the fluid incompressible and obeys the approximation of Boussinesq. These equations are written by using bispherical coordinates system and solved by using a finite difference method. The results show the topology of flow is strongly dependent on the inclinaison because the flow can change from a unicellular regime to a multicellular regime by varying the inclination from 0 to π. By increasing the Rayleigh number (10<sup>3</sup><<i>Ra</i><10<sup>7</sup>), the flow intensifies. T<span style="letter-spacing:-0.05pt;">he results are shown in terms of streamlines and isotherms during th</span>eir transient evolution.展开更多
The so-called blisks,i.e.integrally bladed disks,are characterized by very low viscous material damping and make the flutter prediction much more critical.In that framework,a two-dimensional numerical study of a space...The so-called blisks,i.e.integrally bladed disks,are characterized by very low viscous material damping and make the flutter prediction much more critical.In that framework,a two-dimensional numerical study of a space turbine blisk featuring complex deformation of blades and high eigenfrequency(>40kHz)is performed.The simulations are based on unsteady Reynolds Averaged Navier Stokes computations linearized in the frequency domain and consist in the superposition of an unsteady linear(in time)pressure field,generated by a harmonic perturbation,upon a steady nonlinear(in space)flow.The aerodynamic damping coefficient is calculated over a range of nodal diameters,and the blades are predicted aeroelastically stable.However,violent changes occur and are rather critical since sudden and large deviations in stability appear.In that context,the nature of the waves propagating from the cascade are evaluated.Such an approach provides fundamental knowledge about the perturbations which can either propagate to the far-field(cut-on mode)or decay(cut-off mode).It is expected that the ability of the flow to damp or to amplify the blade motion is strongly affected by the way unsteady perturbations are transferred from the cascade to the far-field.The nature of the waves are first assessed from the aforementioned linearized results,then they are evaluated analytically and finally compared.A good agreement is found despite the strong assumptions of the analytical model.The results show a clear correlation between the cut-on/cut-off conditions and stability.The least stable configuration corresponds to cut-off mode at the inlet and no wave at the outlet.Without outgoing waves from the cascade,the blade is prone to be less stable:the energy from the blades vibration is necessarily dissipated or sent out by the cascade.展开更多
The temporal behaviour of a flow separation in the hub-suction side corner of a transonic diffuser is studied thanks to unsteady numerical simulations based on the phase-lagged approach. The validity of the numerical ...The temporal behaviour of a flow separation in the hub-suction side corner of a transonic diffuser is studied thanks to unsteady numerical simulations based on the phase-lagged approach. The validity of the numerical results is confirmed by comparison with experimental unsteady pressure measurements. An analysis of the instantaneous skin-friction pattern and particles trajectories is presented. It highlights the topology of the separation and its temporal behaviour. The major result is that, despite of a highly time-dependent core flow, the separation is found to be a "fixed unsteady separation" characterized by a fixed location of the main saddle of the separation but an extent of the stall region modulated by the pressure waves induced by the impeller-diffuser interaction.展开更多
A modal analysis method of the rotor-stator interactions in multistage compressors has been developed by LMFA.This method,based on a double modal decomposition of the flow over space and time,has been applied to nu-me...A modal analysis method of the rotor-stator interactions in multistage compressors has been developed by LMFA.This method,based on a double modal decomposition of the flow over space and time,has been applied to nu-merical and experimental results of the high-speed 3?-stage compressor CREATE based at LMFA,Lyon-France.It reveals the presence of a very strong rotor-stator interaction which completely drives the flow at casing behind all the rotors.This modal analysis method applied to an unsteady RANS simulation permits to calculate the en-ergy of the rotor-stator interactions and to plot energetic meridian maps to explain experimental results and to analyze the interaction in the whole machine.展开更多
Full annulus simulations of the flow which develops in a transonic centrifugal compressor are performed at two stable operating points(peak efficiency and near surge) and during the path to surge. At stable conditions...Full annulus simulations of the flow which develops in a transonic centrifugal compressor are performed at two stable operating points(peak efficiency and near surge) and during the path to surge. At stable conditions, the flow field properties are analyzed by comparisons with experimental data and numerical simulations using a phase lagged approach previously carried out. Regarding the stage overall performance, an excellent agreement is obtained between the numerical results(both with time lagged approach and full-annulus calculation) and the experiments. From the full-annulus simulations, the change in flow pattern from peak efficiency to surge is found to be perfectly similar to that obtained from the simulations using the time lagged approach. In particular, provided that the operating point is stable, the flow proves to be chorochronic. The full-annulus simulations were continued after a unique small change in the throttle law applied at the exit of the numerical domain. The mass flow, pressure ratio and efficiency then significantly drop all the more the time progresses. The simulation becomes unstable and the surge inception well underway. The path to surge is found to be due to the enlargement of the boundary layer separation on the suction side of the diffuser vanes in accordance with the conclusions drawn from the chorochronic simulations and experiments. But as the time progresses, the flow loses its chorochronic character. Stall cells rotating at around 7% of the rotor speed develop and lead to surge in around 5 revolutions.展开更多
This paper presents a numerical study of the flow topologies of three-dimensional (3D) flows in a high pressure compressor stator blade row without and with boundary layer aspiration on the hub wall. The stator blad...This paper presents a numerical study of the flow topologies of three-dimensional (3D) flows in a high pressure compressor stator blade row without and with boundary layer aspiration on the hub wall. The stator blade is representative of the first stage operating under transonic inlet conditions and the blade design encourages development of highly complex 3D flows. The blade has a small tip clearance. The computational fluid dynamics (CFD) studies show progressive increase of hub corner stall with the increase in incidence. Aspiration is implemented on the hub wall via a slot in the comer between the hub wall and the suction surface. The CFD studies show aspiration to be sensitive to the suction flow rate; lower rate leads to very complex flow struc- tures and increased level of losses whereas higher rate renders aspiration effective for control of hub comer separation. The flow topologies are studied by trace of skin friction lines on the walls. The nature of flow can be explained by the topological rules of closed separation. Furthermore, a deeper analysis is done for a particular case with advanced criterion to test the non-degeneracy of critical points in the flow field.展开更多
Numerical and experimental investigations were conducted in a transonic centrifugal compressor stage composed of a backswept splittered unshrouded impeller and a vaned diffuser. A detailed analysis of the flow in the ...Numerical and experimental investigations were conducted in a transonic centrifugal compressor stage composed of a backswept splittered unshrouded impeller and a vaned diffuser. A detailed analysis of the flow in the inducer (i.e. the entry zone of the impeller between the main blade leading edge and the splitter blade leading edge) is proposed from choke to surge. Steady and unsteady simulations were performed using the code elsA, which uses a multi-domain approach on structured meshes and solves the compressible RANS equations, associated with a two-equation turbulence model k-l in the rotating frame of reference. The 1MW LMFA-ECL test rig was used for carrying out the tests in the compressor stage. Unsteady pressure measurements up to 150 kHz and Laser Doppler Anemometry measurements were performed in the inducer. A good agreement is obtained between the experimental and numerical data even if an over dissipation is noticed in the numerical results. The change in flow pattern from choke to surge is mainly due to a change in the tip leakage flow trajectory which straightens, leading to a flow blockage of an individual passage near shroud. A spectral analysis shows that only the blade passing frequency and its harmonics compose the various spectra obtained from choke to surge.展开更多
This paper presents the unsteady data acquisition system used to measure the pressure field in high speed compressors.Details and electronic sketches are given for the conditioners developed in-house that have been us...This paper presents the unsteady data acquisition system used to measure the pressure field in high speed compressors.Details and electronic sketches are given for the conditioners developed in-house that have been used to amplify and to filter the pressure signal with the aim of acquiring data up to 150 kHz.A discussion of the experimental results carried out in a centrifugal compressor is proposed.Through different processing of the pressure signals and a comparison with URANS simulations,the excitation of the pressure transducers by the pressure waves generated by shock waves that occur between the impeller and the diffuser is highlighted.The levels of pressure fluctuations measured when entering into surge are also presented and reveal very repetitive behaviour of the flow instabilities.展开更多
Doppler global velocimetry (DGV) has already been shown to be an interesting technique capable of measuringthe three components of velocity in a plane. A 1-component DGV system is currently under development at theLMF...Doppler global velocimetry (DGV) has already been shown to be an interesting technique capable of measuringthe three components of velocity in a plane. A 1-component DGV system is currently under development at theLMFA, using a stabilized continuous wave (CW) argon ion laser for emission. The receiver features only onecamera for both signal and reference images and incorporates a DEFI system to adjust the incident laser lightfrequency and its transmission coefficient on the iodine cell absorption line. A description of the whole system ispresented and a validation with measurements of axial velocities at several positions in a round free jet is proposed.展开更多
In the case of three-dimensional flows, the separation can be defined in more than one way. Discussions about three-dimensional boundary layer separation in the literature have found a rational extension of the zero s...In the case of three-dimensional flows, the separation can be defined in more than one way. Discussions about three-dimensional boundary layer separation in the literature have found a rational extension of the zero skin friction. Attempts have been made to establish the identity of “separation lines”. Among definitions, these may be found: (1) envelopes of limiting streamlines, (2) lines dividing flow which has come from different regions, (3) lines of singularities (problems of topology), (4) lines on which some component of the skin friction vanishes. Each of these is valid under certain conditions, but none is universally valid. In the present work, we use the definition (4), i.e. at low incidence of a blunt body the separation line is identified as the zero of the meridian skin friction component (Wang 1975). So the separation line on a flattened spheroid (6:3:1) at 6° of incidence is calculated, as well as experimentally determined by using the electrochemical method, which allows to follow the evolution of the parietal velocity gradient.展开更多
文摘This study aims to analyze mixed convection in a square cavity with two moving vertical walls by finite volume method.The cavity filled with Non-Newtonian fluid of Bingham model is heated from below and cooled by the other walls.This study has been conducted for certain parameters of Reynolds number(Re=1-100),Richardson number(Ri=1-20),Prandtl number(Pr=1-500),and Bingham number has been studied from 0 to 10.The results indicate that the increase in yield stress drops the heat transfer and the flow become flatter,while increasing Reynolds number augments it.The convective transport is dominant when increasing Richardson number which leads to enhance heat transfer in the cavity for both Newtonian and Non-Newtonian fluid.A correlation of Nusselt number is given in function of different parameters.
文摘The increase of wave energy in electricity production is an objective shared by many countries to meet growing demand and global warming. To analyze devices capable of converting the energy of sea waves into electrical energy, it is important to master the various theories of gravity waves and generation. We will in our work consider a numerical waves tank for an amplitude A=0.5, a wavelength λ=0.25 , an average height H<sub>e</sub>=10 and a Froude number fixed at 1 × 10<sup>5</sup>. Numerical wave channel analysis is used to reproduce the natural phenomenon of wave propagation in an experimental model. Wave makers are usually used to generate waves in the channel. In theory, the influence of an incident wave can be considered, as in the case of our study. In this study, the evolution of the hydrodynamic parameters and the energy transported in one wavelength can be determined by calculation. A change of variable will be done in this work to facilitate the writing of the boundary conditions at the free surface and at the bottom. The nonlinear Stokes theory will be studied in this case in order to provide hydrodynamic solutions through the Navier-Stokes equations to finally deduce the energetic results. To do this, the finite difference method will be used for the hydrodynamic results such as the velocity potential and the free surface elevation and the trapezium method of Newton for the energetic results. Thus, we will determine the energetic potential according to the decrease in the slope of the tank. To do this, we will take as values of beta representing the inverse of the slope of the tank, β=100, β=105, β=110 and β=105. .
基金supported by the Sino-French Project AX-IOOM (Advanced Experiments and Simulations of Complex Flows in Turbomachines)the National Natural Science Foundation of China (51136003, 50976010)the National Basic Research Program of China (2012CB720205)
文摘It is widely accepted that in a turbulent boundary layer (TBL) with adverse pressure gradient (APG) an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this outer peak is not clearly understood. In this paper, the effect of APG is analysed using the numerical and experimental results in the literature. Because the effect of upstream flow is inherent in the TBL, we first analyse this effect in TBLs with zero pressure gradient on flat plates. Under the individual effect of upstream flow, an outer peak already appears in the profile of streamwise Reynolds stress when the TBL continues developing in the streamwise direction. The APG accelerates the appearance of the outer peak, instead of being a trigger.
文摘Natural convection heat transfer in open or closed cavities takes place in different engineering areas. The hemispherical cavity is a part of basic geometries although it is not widely studied. The present paper reports the numerical study of natural convection in a closed hemispherical annulus delimited by two vertically eccentric hemispheres filled with Newtonian fluid (air in this case with <em>Pr</em> = 0.7) is conducted. The inner hemisphere is heated by a heat flux of constant density and the outer one is maintained isothermal. Based on the Boussinesq assumptions, the governing equations are numerically studied using unsteady natural convection formulated with vorticity and stream-function variables. These equations are written by using bispherical coordinates system and solved by using a finite difference method. The effect of the control parameters such as the Rayleigh number (<span style="white-space:nowrap;">10<sup>3</sup> ≤ <em>Ra</em> ≤ 10<sup>6</sup></span>) or the eccentricity (<em>e</em> = ±0.2, ±0.5, 0) in the dynamic and thermal behaviours of the fluid is investigated.
基金performed in the scope of project ANR-12-ASTR-0017-03 "BF-DRAINH" in collaboration with the IRENav Laboratory (French Naval Academy,Brest,France) and the IMFT laboratory (Toulouse,France)
文摘Contrary to natural cavitation,ventilated cavitation is controllable and is not harmful.It is particularly used to reduce the drag of the hydraulic vehicles.The ventilated cavitation is characterized by various gas regimes.The mechanisms of ventilated cavitation are investigated in the present work with CFD based on a 2D solver.The attention is especially focused on the transition between the reentrant jet and twin vortex regimes.The results confirmthat the product of ventilated cavitation number and Froude number is lower than 1(σcFr<1)in the twin vortex regime,while it is higher than 1(σcFr>1)in the reentrant jet regime,as reported in the literature.Further analysis shows that ventilated cavitation is significantly influenced by the natural cavitation number.
文摘Abstract: This work studies the active control of chemical oscillations governed by a forced modified Van der Pol-Duffing oscillator. We considered the dynamics of nonlinear chemical systems subjected to an external sinusoidal excitation. The approximative solution to the first order of the modified Van der Pol-Duffing oscillator is found using the Lindstedt’s perturbation method. The harmonic balance method is used to find the amplitudes of the oscillatory states of the system under control. The effects of the constraint parameter and the control parameter of the model on the amplitude of oscillations are presented. The effects of the active control on the behaviors of the model are analyzed and it appears that with the appropriate selection of the coupling parameter, the chaotic behavior of the model has given way to periodic movements. Numerical simulations are used to validate and complete the analytical results obtained.
文摘In our study, we investigate the differences between the combustion of different hydrocarbon fuels CH4, C3H8, C4H10. A numerical simulation of an impinging jet diffusion flames is used. The jet injector has a 10 mm in diameter and the distance between the jet flame and the vertical wall is 2 time half diameter. The fuel jet velocity was fixed for 11.8 m/s, corresponding to a Reynolds number of 6881. The flame characteristics varied from hydrocarbon to another for the same Reynolds number. The combustion products of CO, CO2, NO, OH, are depending on the methane and propane and butane flames for the same conditions. The temperature of the flame was varied from hydrocarbon to another the same as for the chemical species production rate. The concentration of the thermal and prompt NO pollutant depends on the temperature flow field and on the thermochemical characteristics of the hydrocarbon fuels.
文摘The subject of natural convection heat transfer is motivated by a wide range of applications in engineering technology. The hemispherical cavity is a part of basic geometries although it is not widely studied. The effect of inclinaison on natural convection fluid motions in the gap between two eccentric hemispheres is numerically studied. The inner hemisphere is subjected to a heat flux of a constant density and the outer one is maintened isothermal. The walls separating the two hemispheres are thermally adiabatic. Equations are formulated with vorticity and stream-functions variables. It is also assumed the fluid incompressible and obeys the approximation of Boussinesq. These equations are written by using bispherical coordinates system and solved by using a finite difference method. The results show the topology of flow is strongly dependent on the inclinaison because the flow can change from a unicellular regime to a multicellular regime by varying the inclination from 0 to π. By increasing the Rayleigh number (10<sup>3</sup><<i>Ra</i><10<sup>7</sup>), the flow intensifies. T<span style="letter-spacing:-0.05pt;">he results are shown in terms of streamlines and isotherms during th</span>eir transient evolution.
基金the Centre National d'Etudes Spatiales (CNES) and Snecma for their financial supportthe Centre Informatique National de l'Enseignement Supérieur (CINES) for the computational resources,and the Agence Nationale de la Recherche(ANR) for sponsoring the project ANR-08-2009 CapCAO (parametrization-aided optimized aeroelastic design)
文摘The so-called blisks,i.e.integrally bladed disks,are characterized by very low viscous material damping and make the flutter prediction much more critical.In that framework,a two-dimensional numerical study of a space turbine blisk featuring complex deformation of blades and high eigenfrequency(>40kHz)is performed.The simulations are based on unsteady Reynolds Averaged Navier Stokes computations linearized in the frequency domain and consist in the superposition of an unsteady linear(in time)pressure field,generated by a harmonic perturbation,upon a steady nonlinear(in space)flow.The aerodynamic damping coefficient is calculated over a range of nodal diameters,and the blades are predicted aeroelastically stable.However,violent changes occur and are rather critical since sudden and large deviations in stability appear.In that context,the nature of the waves propagating from the cascade are evaluated.Such an approach provides fundamental knowledge about the perturbations which can either propagate to the far-field(cut-on mode)or decay(cut-off mode).It is expected that the ability of the flow to damp or to amplify the blade motion is strongly affected by the way unsteady perturbations are transferred from the cascade to the far-field.The nature of the waves are first assessed from the aforementioned linearized results,then they are evaluated analytically and finally compared.A good agreement is found despite the strong assumptions of the analytical model.The results show a clear correlation between the cut-on/cut-off conditions and stability.The least stable configuration corresponds to cut-off mode at the inlet and no wave at the outlet.Without outgoing waves from the cascade,the blade is prone to be less stable:the energy from the blades vibration is necessarily dissipated or sent out by the cascade.
文摘The temporal behaviour of a flow separation in the hub-suction side corner of a transonic diffuser is studied thanks to unsteady numerical simulations based on the phase-lagged approach. The validity of the numerical results is confirmed by comparison with experimental unsteady pressure measurements. An analysis of the instantaneous skin-friction pattern and particles trajectories is presented. It highlights the topology of the separation and its temporal behaviour. The major result is that, despite of a highly time-dependent core flow, the separation is found to be a "fixed unsteady separation" characterized by a fixed location of the main saddle of the separation but an extent of the stall region modulated by the pressure waves induced by the impeller-diffuser interaction.
基金the CNRS and the company Snecma (SAFRAN) which support the compressor CREATE research program
文摘A modal analysis method of the rotor-stator interactions in multistage compressors has been developed by LMFA.This method,based on a double modal decomposition of the flow over space and time,has been applied to nu-merical and experimental results of the high-speed 3?-stage compressor CREATE based at LMFA,Lyon-France.It reveals the presence of a very strong rotor-stator interaction which completely drives the flow at casing behind all the rotors.This modal analysis method applied to an unsteady RANS simulation permits to calculate the en-ergy of the rotor-stator interactions and to plot energetic meridian maps to explain experimental results and to analyze the interaction in the whole machine.
基金the HPC resources of CINES under the allocation 2012- 2a6356 and 2013-2a6356
文摘Full annulus simulations of the flow which develops in a transonic centrifugal compressor are performed at two stable operating points(peak efficiency and near surge) and during the path to surge. At stable conditions, the flow field properties are analyzed by comparisons with experimental data and numerical simulations using a phase lagged approach previously carried out. Regarding the stage overall performance, an excellent agreement is obtained between the numerical results(both with time lagged approach and full-annulus calculation) and the experiments. From the full-annulus simulations, the change in flow pattern from peak efficiency to surge is found to be perfectly similar to that obtained from the simulations using the time lagged approach. In particular, provided that the operating point is stable, the flow proves to be chorochronic. The full-annulus simulations were continued after a unique small change in the throttle law applied at the exit of the numerical domain. The mass flow, pressure ratio and efficiency then significantly drop all the more the time progresses. The simulation becomes unstable and the surge inception well underway. The path to surge is found to be due to the enlargement of the boundary layer separation on the suction side of the diffuser vanes in accordance with the conclusions drawn from the chorochronic simulations and experiments. But as the time progresses, the flow loses its chorochronic character. Stall cells rotating at around 7% of the rotor speed develop and lead to surge in around 5 revolutions.
文摘This paper presents a numerical study of the flow topologies of three-dimensional (3D) flows in a high pressure compressor stator blade row without and with boundary layer aspiration on the hub wall. The stator blade is representative of the first stage operating under transonic inlet conditions and the blade design encourages development of highly complex 3D flows. The blade has a small tip clearance. The computational fluid dynamics (CFD) studies show progressive increase of hub corner stall with the increase in incidence. Aspiration is implemented on the hub wall via a slot in the comer between the hub wall and the suction surface. The CFD studies show aspiration to be sensitive to the suction flow rate; lower rate leads to very complex flow struc- tures and increased level of losses whereas higher rate renders aspiration effective for control of hub comer separation. The flow topologies are studied by trace of skin friction lines on the walls. The nature of flow can be explained by the topological rules of closed separation. Furthermore, a deeper analysis is done for a particular case with advanced criterion to test the non-degeneracy of critical points in the flow field.
基金performed using HPC resources from GENCI-[CCRT/CINES/IDRIS] (Grant2010-i2010026356 and 2011-x2011026356)
文摘Numerical and experimental investigations were conducted in a transonic centrifugal compressor stage composed of a backswept splittered unshrouded impeller and a vaned diffuser. A detailed analysis of the flow in the inducer (i.e. the entry zone of the impeller between the main blade leading edge and the splitter blade leading edge) is proposed from choke to surge. Steady and unsteady simulations were performed using the code elsA, which uses a multi-domain approach on structured meshes and solves the compressible RANS equations, associated with a two-equation turbulence model k-l in the rotating frame of reference. The 1MW LMFA-ECL test rig was used for carrying out the tests in the compressor stage. Unsteady pressure measurements up to 150 kHz and Laser Doppler Anemometry measurements were performed in the inducer. A good agreement is obtained between the experimental and numerical data even if an over dissipation is noticed in the numerical results. The change in flow pattern from choke to surge is mainly due to a change in the tip leakage flow trajectory which straightens, leading to a flow blockage of an individual passage near shroud. A spectral analysis shows that only the blade passing frequency and its harmonics compose the various spectra obtained from choke to surge.
基金the company Turbomeca which supports the centrifugal compressor research program
文摘This paper presents the unsteady data acquisition system used to measure the pressure field in high speed compressors.Details and electronic sketches are given for the conditioners developed in-house that have been used to amplify and to filter the pressure signal with the aim of acquiring data up to 150 kHz.A discussion of the experimental results carried out in a centrifugal compressor is proposed.Through different processing of the pressure signals and a comparison with URANS simulations,the excitation of the pressure transducers by the pressure waves generated by shock waves that occur between the impeller and the diffuser is highlighted.The levels of pressure fluctuations measured when entering into surge are also presented and reveal very repetitive behaviour of the flow instabilities.
基金the CIRT (Consortium Indus trie Recherche en Turbomachines) which is supporting this workmore particularly at the LMFA, Alain Ef fernelli for building many elements of the receiver, Gil bert Halter for helping in the collimation processAndre Vouillarmet for his advice in this project.
文摘Doppler global velocimetry (DGV) has already been shown to be an interesting technique capable of measuringthe three components of velocity in a plane. A 1-component DGV system is currently under development at theLMFA, using a stabilized continuous wave (CW) argon ion laser for emission. The receiver features only onecamera for both signal and reference images and incorporates a DEFI system to adjust the incident laser lightfrequency and its transmission coefficient on the iodine cell absorption line. A description of the whole system ispresented and a validation with measurements of axial velocities at several positions in a round free jet is proposed.
文摘In the case of three-dimensional flows, the separation can be defined in more than one way. Discussions about three-dimensional boundary layer separation in the literature have found a rational extension of the zero skin friction. Attempts have been made to establish the identity of “separation lines”. Among definitions, these may be found: (1) envelopes of limiting streamlines, (2) lines dividing flow which has come from different regions, (3) lines of singularities (problems of topology), (4) lines on which some component of the skin friction vanishes. Each of these is valid under certain conditions, but none is universally valid. In the present work, we use the definition (4), i.e. at low incidence of a blunt body the separation line is identified as the zero of the meridian skin friction component (Wang 1975). So the separation line on a flattened spheroid (6:3:1) at 6° of incidence is calculated, as well as experimentally determined by using the electrochemical method, which allows to follow the evolution of the parietal velocity gradient.