Numerical Resolution of Compressor Corner Separation Problem

Author(s):  
Wei Sun ◽  
Zhong-Nan Wang ◽  
Liping Xu

Abstract Corner separation is a common flow phenomenon within compressors that can significantly affect the compressor performance. The RANS turbulence closures, commonly used in the industrial CFD simulations, often struggle to predict corner separation with reasonable accuracy. In this paper, the results of two RANS-based modelling approaches are presented for the corner separation within a high-loaded Prescribed Velocity Distribution (PVD) compressor cascade. The flow characteristics are studied to facilitate understanding the causes of varying performance of RANS models. It is observed that mixing plays a crucial role in accurately predicting the type, location, and size of flow separation. The source terms that control the turbulence mixing in SA and SST models are identified, based on physical analyses. Both RANS models are modified to better model the mixing process. Based on the modified SST model, an improved RANS-LES blending function has been proposed for a hybrid RANS-LES model. This new blending function ensures reliable shielding of attached boundary layers by the RANS portion of the hybrid model. Finally, to gain further understanding of the endwall flow physics, the turbulence characteristics of the resolved corner separation flow are studied, in terms of the large-scale unsteadiness and loss generation mechanisms.

Author(s):  
Wei Ma ◽  
Feng Gao ◽  
Xavier Ottavy ◽  
Lipeng Lu ◽  
A. J. Wang

Recently bimodal phenomenon in corner separation has been found by Ma et al. (Experiments in Fluids, 2013, doi:10.1007/s00348-013-1546-y). Through detailed and accurate experimental results of the velocity flow field in a linear compressor cascade, they discovered two aperiodic modes exist in the corner separation of the compressor cascade. This phenomenon reflects the flow in corner separation is high intermittent, and large-scale coherent structures corresponding to two modes exist in the flow field of corner separation. However the generation mechanism of the bimodal phenomenon in corner separation is still unclear and thus needs to be studied further. In order to obtain instantaneous flow field with different unsteadiness and thus to analyse the mechanisms of bimodal phenomenon in corner separation, in this paper detached-eddy simulation (DES) is used to simulate the flow field in the linear compressor cascade where bimodal phenomenon has been found in previous experiment. DES in this paper successfully captures the bimodal phenomenon in the linear compressor cascade found in experiment, including the locations of bimodal points and the development of bimodal points along a line that normal to the blade suction side. We infer that the bimodal phenomenon in the corner separation is induced by the strong interaction between the following two facts. The first is the unsteady upstream flow nearby the leading edge whose angle and magnitude fluctuate simultaneously and significantly. The second is the high unsteady separation in the corner region.


Author(s):  
Jean-François Monier ◽  
Nicolas Poujol ◽  
Mathieu Laurent ◽  
Feng Gao ◽  
Jérôme Boudet ◽  
...  

The present study aims at analysing the Boussinesq constitutive relation validity in a corner separation flow of a compressor cascade. The Boussinesq constitutive relation is commonly used in Reynolds-averaged Navier-Stokes (RANS) simulations for turbomachinery design. It assumes an alignment between the Reynolds stress tensor and the zero-trace mean strain-rate tensor. An indicator that measures the alignment between these tensors is used to test the validity of this assumption in a high fidelity large-eddy simulation. Eddy-viscosities are also computed using the LES database and compared. A large-eddy simulation (LES) of a LMFA-NACA65 compressor cascade, in which a corner separation is present, is considered as reference. With LES, both the Reynolds stress tensor and the mean strain-rate tensor are known, which allows the construction of the indicator and the eddy-viscosities. Two constitutive relations are evaluated. The first one is the Boussinesq constitutive relation, while the second one is the quadratic constitutive relation (QCR), expected to render more anisotropy, thus to present a better alignment between the tensors. The Boussinesq constitutive relation is rarely valid, but the QCR tends to improve the alignment. The improvement is mainly present at the inlet, upstream of the corner separation. At the outlet, the correction is milder. The eddy-viscosity built with the LES results are of the same order of magnitude as those built as the ratio of the turbulent kinetic energy k and the turbulence specific dissipation rate ω. They also show that the main impact of the QCR is to rotate the mean strain-rate tensor in order to realign it with the Reynolds stress tensor, without dilating it.


2020 ◽  
Vol 142 (2) ◽  
Author(s):  
Yangwei Liu ◽  
Yumeng Tang ◽  
Ashley D. Scillitoe ◽  
Paul G. Tucker

Abstract Three-dimensional corner separation significantly affects compressor performance, but turbulence models struggle to predict it accurately. This paper assesses the capability of the original shear stress transport (SST) turbulence model to predict the corner separation in a linear highly loaded prescribed velocity distribution (PVD) compressor cascade. Modifications for streamline curvature, Menter’s production limiter, and the Kato-Launder production term are examined. Comparisons with experimental data show that the original SST model and the SST model with different modifications can predict the corner flow well at an incidence angle of −7 deg, where the corner separation is small. However, all the models overpredict the extent of the flow separation when the corner separation is larger, at an incidence angle of 0 deg. The SST model is then modified using the helicity to take account of the energy backscatter, which previous studies have shown to be important in the corner separation regions of compressors. A Reynolds stress model (RSM) is also used for comparison. By comparing the numerical results with experiments and RSM results, it can be concluded that sensitizing the SST model to helicity can greatly improve the predictive accuracy for simulating the corner separation flow. The accuracy is quite competitive with the RSM, whereas in terms of computational cost and robustness it is superior to the RSM.


Author(s):  
Tan Zheng ◽  
Xiaoqing Qiang ◽  
Jinfang Teng

This paper presents a numerical investigation to explore the effects of vortex generators on a high speed compressor cascade. Secondary flow effects like the corner separation vortex have an influence on the performance of a compressor cascade such as leading to increased losses. In order to control the corner separation vortex and reduce losses, an extensive study of vortex generators applied to a compressor cascade is conducted. A preliminary study by steady 3D RANS simulations is performed using the Spalart-Allmaras turbulence model. The aerodynamic performance as well as the behavior of the corner separation vortex is investigated in the compressor cascade without vortex generators. Then, a vortex generator is added to the cascade, which is numerically simulated. Various configurations are considered, which are decided by the height and installation angle of the vortex generator. Comparison of the performance attained by these configurations results in an optimum scheme that has minimum losses. Furthermore, unsteady 3D DES simulations are performed with the optimum configuration. This method that predicts the flow field more precisely could help verify the accuracy of the RANS results. Finally, by analyzing all the resulting aerodynamic performance and numerical flow phenomena, the mechanism of vortex-vortex interaction is presented and discussed, which could be a criterion to reduce the corner separation flow and enhance the performance of axial compressors.


2021 ◽  
pp. 1-30
Author(s):  
Jerome Boudet ◽  
Emmanuel Lévêque ◽  
Hatem Touil

Abstract Lattice-Boltzmann simulations of corner separation flow in a compressor cascade are presented. The lattice Boltzmann approach is rather new in the context of turbomachinery and the configuration is known to be particularly challenging for turbulence modelling. The present methodology is characterized by a quasi-autonomous meshing strategy and a limited computational cost (a net ratio of 5 compared to a previous finite-volume compressible Navier-Stokes simulation). The simulation of the reference case (4° incidence) shows a good agreement with the experimental data concerning the wall pressure distribution or the distribution of losses. A good description is also obtained when incidence angle is increased to 7°, with a span-wise development of the separation. Subsequently, the methodology is used to investigate the sensitivity of the flow to the end-wall boundary-layer thickness. A thinner boundary-layer results in a smaller corner separation, but not a complete elimination. Finally, the ingredients of the wall modelling are analysed in details. On the one hand, the curvature correction term promotes transition to turbulence on the blade suction side and avoids a spurious separation. On the other hand, the addition of the pressure-gradient correction term allows a wider and more realistic corner separation.


Author(s):  
Bingxiao Lu ◽  
Jinfang Teng ◽  
Mingmin Zhu ◽  
Xiaoqing Qiang ◽  
Wei Ma

Abstract Under the condition of large incidence angle in an axial compressor, corner separation will occur in corner region. When the blade loading increases, there may even be corner stall. This three dimensional complex flow structure is high intermittent and unsteady. In order to study the flow mechanisms, a hybrid RANS/LES turbulence model (SBES) was used to simulate the corner separation in the cascade. Firstly, the evolution of vortical structures in the corner separation region were analyzed. The boundary layer and the backflow in the corner separation region encounter, forming vortices with opposite rotation and developing downstream continuously. This cause the unsteady flow in corner separation. In addition, the shedding of detached eddies at the trailing edge is another sources of unsteadiness. Secondly, based on studying the turbulent characteristics, it can be seen that there is active turbulent fluctuation in the corner separation region. The turbulence is high anisotropic because the distribution and the values of Reynolds stresses in different directions are quite different. Thirdly, the mechanism of bimodal phenomenon was explained in this paper. Bimodal means there are two peaks in the velocity probability histogram, it is found in the region between separated flow and non-separated flow. This indicates that there are two non-periodic switching flow modes. Through the analysis of this paper, we find that it is related to the development of large scale coherent vortical structures.


1969 ◽  
Vol 9 (04) ◽  
pp. 403-411 ◽  
Author(s):  
B.K. Sinha ◽  
Harvey T. Kennedy

Abstract Recommendations are made for obtaining consistent and reproducible test data on drilling fluids having identical composition. Previously, such a procedure has been difficult to accomplish even when the fluids were mixed in similar equipment. A survey of work in this area indicates that previous methods have been unsatisfactory because previous methods have been unsatisfactory because (1) the muds are extremely sensitive to the duration and violence of agitation during a normal mixing routine, and (2) gelling of the muds occurs before the properties can reach constant values. This gelling is caused by water evaporation resulting from the increase in temperature associated with the agitation. The work shows that these problems largely can be overcome by (1) agitating the constituents of the drilling fluid more vigorously, (2) maintaining a fairly constant temperature, and(3) Protecting the fluid from evaporation. When these steps are followed, the fluid properties approach asymptotic values that do not change by prolonged or accelerated agitation or by aging for a month. The time required to reach asymptotic values or a stabilized state is from 2 to 6 hours and is a function of the mud composition. Introduction Preparation of drilling fluids in the laboratory to determine their suitability to meet specific drilling requirements or to serve as a base fluid to evaluate the effectiveness of thinners, dispersants or other additives normally begins with combining measured quantities of the constituents and stirring them for a short time in a low-speed mixer. This is done to obtain a uniform mixture and to hydrate clays. Then the fluid is further agitated in a higher-speed device (Hamilton Beach mixer or Waring blender) to disperse more thoroughly and clay particles The biggest obstacle in the laboratory investigation of drilling fluids has been the lack of a method of producing a mixture by which reproducible results of the measured properties could be obtained. Numerous investigators have encountered this difficulty. Prior to 1929, density was the only property of mud that customarily was measured. The use of Wyoming bentonite on a large scale after 1929 was mainly responsible for the development of more elaborate testing procedures and for the application of the principles of colloid chemistry to the drilling fluids. Ambrose and Loomis in 1931 were among the first to recognize the plastic flow characteristics of drilling fluids, although Bingham in 1916 had observed The same phenomenon with dilute clay suspensions. Marsh introduced the Marsh funnel for field testing in 1931. By this time, non-Newtonian characteristics of drilling fluids were established. The Stormer and MacMichael viscometers were used to study the rheological properties of the fluids. In the 1930's and early 1940's, the work conducted by several investigators contributed toward a better understanding of drilling fluids. In the mid 1930's, fluid-loss and the associated mud-cake-forming properties of drilling fluids were recognized as important to the behavior of these fluids. The other properties of drilling fluids, including gel strength, pH, and sand content soon were recognized. In 1937, API published its first recommended procedure for test methods. Since that time, these procedures have been revised periodically. The latest edition, RP-13B, was published in 1961 However, in spite of the recognized need for a method of mixing that provides drilling fluids with stabilized properties, no such method previously has been described. SPEJ P. 403


2019 ◽  
Vol 141 (4) ◽  
Author(s):  
H. D. Akolekar ◽  
J. Weatheritt ◽  
N. Hutchins ◽  
R. D. Sandberg ◽  
G. Laskowski ◽  
...  

Nonlinear turbulence closures were developed that improve the prediction accuracy of wake mixing in low-pressure turbine (LPT) flows. First, Reynolds-averaged Navier–Stokes (RANS) calculations using five linear turbulence closures were performed for the T106A LPT profile at isentropic exit Reynolds numbers 60,000 and 100,000. None of these RANS models were able to accurately reproduce wake loss profiles, a crucial parameter in LPT design, from direct numerical simulation (DNS) reference data. However, the recently proposed kv2¯ω transition model was found to produce the best agreement with DNS data in terms of blade loading and boundary layer behavior and thus was selected as baseline model for turbulence closure development. Analysis of the DNS data revealed that the linear stress–strain coupling constitutes one of the main model form errors. Hence, a gene-expression programming (GEP) based machine-learning technique was applied to the high-fidelity DNS data to train nonlinear explicit algebraic Reynolds stress models (EARSM), using different training regions. The trained models were first assessed in an a priori sense (without running any RANS calculations) and showed much improved alignment of the trained models in the region of training. Additional RANS calculations were then performed using the trained models. Importantly, to assess their robustness, the trained models were tested both on the cases they were trained for and on testing, i.e., previously not seen, cases with different flow features. The developed models improved prediction of the Reynolds stress, turbulent kinetic energy (TKE) production, wake-loss profiles, and wake maturity, across all cases.


2018 ◽  
Vol 611 ◽  
pp. A89 ◽  
Author(s):  
Yueh-Ning Lee ◽  
Patrick Hennebelle

Context. Understanding the origin of the initial mass function (IMF) of stars is a major problem for the star formation process and beyond. Aim. We investigate the dependence of the peak of the IMF on the physics of the so-called first Larson core, which corresponds to the point where the dust becomes opaque to its own radiation. Methods. We performed numerical simulations of collapsing clouds of 1000 M⊙ for various gas equations of state (eos), paying great attention to the numerical resolution and convergence. The initial conditions of these numerical experiments are varied in the companion paper. We also develop analytical models that we compare to our numerical results. Results. When an isothermal eos is used, we show that the peak of the IMF shifts to lower masses with improved numerical resolution. When an adiabatic eos is employed, numerical convergence is obtained. The peak position varies with the eos, and using an analytical model to infer the mass of the first Larson core, we find that the peak position is about ten times its value. By analyzing the stability of nonlinear density fluctuations in the vicinity of a point mass and then summing over a reasonable density distribution, we find that tidal forces exert a strong stabilizing effect and likely lead to a preferential mass several times higher than that of the first Larson core. Conclusions. We propose that in a sufficiently massive and cold cloud, the peak of the IMF is determined by the thermodynamics of the high-density adiabatic gas as well as the stabilizing influence of tidal forces. The resulting characteristic mass is about ten times the mass of the first Larson core, which altogether leads to a few tenths of solar masses. Since these processes are not related to the large-scale physical conditions and to the environment, our results suggest a possible explanation for the apparent universality of the peak of the IMF.


2021 ◽  
Vol 13 (23) ◽  
pp. 13206
Author(s):  
Luis Rodriguez-Lucas ◽  
Chen Ning ◽  
Marcelo Fajardo-Pruna ◽  
Yugui Yang

This paper presents a new concept called the urban vortex system (UVS). The UVS couples a vortex generator (V.G.) that produces updraft by artificial vortex and a vortex stability zone (VSZ) consisting of an assembly of four buildings acting as a chimney. Through this system, a stable, upward vortex flow can be generated. The Reynolds Averaged Navier–Stokes (RANS) simulation was carried out to investigate the flow field in the UVS. The Renormalized Group (RNG) k–ε turbulent model was selected to solve the complex turbulent flow. Validation of the numerical results was achieved by making a comparison with the large-size experimental model. The results reported that a steady-state vortex could be formed when a vapor-air mixture at 2 m/s and 450 K enters the vortex generator. This vortex presented a maximum negative central pressure of −6.81 Pa and a maximum velocity of 5.47 (m/s). Finally, the similarity method found four dimensionless parameters, which allowed all the flow characteristics to be transported on a large scale. The proposed large-scale UVS application is predicted to be capable, with have a maximum power of 2 M.W., a specific work of 3 kJ/kg, buildings 200-m high, and the ability to generate winds of 6.1 m/s (20 km/h) at 200 m up to winds of 1.5 m/s (5 km/h) at 400 m. These winds would cause the rupture of the gas capsule of the heat island phenomenon. Therefore, the city would balance its temperature with that of the surrounding rural areas.


Sign in / Sign up

Export Citation Format

Share Document