Effects of Inlet Disturbances on Fan Stability

Author(s):  
Kuen-Bae Lee ◽  
Mark Wilson ◽  
Mehdi Vahdati

This paper investigates the effects of inlet disturbance caused by crosswind on a fan blade operation and addresses the possible aerodynamic instabilities, which can arise for such a fan-intake system. The work is carried out by using a three-dimensional unsteady computational fluid dynamics (CFD) model (AU3D), and for a modern low-speed fan rig for which extensive measured data is available. The computational domain includes the fan with outlet guide vanes for the bypass flow, engine-section stators for the core, and a symmetric intake upstream of the fan (a whole low-pressure domain). The unsteady full annulus simulations under crosswind are performed to analyze the effects of inlet disturbances on the operation of this blade. It was observed that, for sufficiently high amplitudes of crosswind, the intake lip separates, and results in a significant loss of stall margin. Moreover, even in the absence of lip separation, the blade can still stall prematurely due to nonhomogeneous flow caused by the two contra-rotating trailing vortices. In the second phase of this study, the effects of fan loading on the suppression of flow separation in the intake, and the consequent stall margin of the fan blade, were explored. The results indicated that, as the fan speed increases, it becomes more capable of reducing the inlet distortion levels, and consequently, the loss in stall margin decreases.

Author(s):  
Kuen-Bae Lee ◽  
Mark Wilson ◽  
Mehdi Vahdati

This paper investigates the effects of inlet disturbance caused by crosswind on a fan blade operation and addresses the possible aerodynamic instabilities which can arise for such a fan-intake system. The work is carried out by using a three-dimensional unsteady CFD model (AU3D), and for a modern low-speed fan rig for which extensive measured data is available. The computational domain includes the fan, with outlet guide vanes (OGV) for the bypass flow and engine-section stators (ESS) for the core and a symmetric intake upstream of the fan (a whole low-pressure domain). The unsteady full annulus simulations under crosswind are performed to analyze the effects of inlet disturbances on the operation of this blade. It was observed that, for sufficiently high amplitudes of crosswind, the intake lip separates, and results in a significant loss of stall margin. Moreover, even in the absence of lip separation, the blade can still stall prematurely due to non-homogeneous flow caused by the two contra-rotating trailing vortices. In the 2nd phase of this study, the effects of fan loading on the suppression of flow separation in the intake, and the consequent stall margin of the fan blade was explored. The results indicated that, as the fan speed increases it becomes more capable of reducing the inlet distortion levels, and consequently it can extend the stall margin of the fan blade.


2004 ◽  
Vol 127 (2) ◽  
pp. 349-351 ◽  
Author(s):  
M. Vahdati ◽  
A. I. Sayma ◽  
C. Freeman ◽  
M. Imregun

This paper describes a novel way of prescribing computational fluid dynamics (CFD) boundary conditions for axial-flow compressors. The approach is based on extending the standard single passage computational domain by adding an intake upstream and a variable nozzle downstream. Such a route allows us to consider any point on a given speed characteristic by simply modifying the nozzle area, the actual boundary conditions being set to atmospheric ones in all cases. Using a fan blade, it is shown that the method not only allows going past the stall point but also captures the typical hysteresis loop behavior of compressors.


2011 ◽  
Vol 133 (12) ◽  
Author(s):  
Feng Zhou ◽  
Nicholas E. Hansen ◽  
David J. Geb ◽  
Ivan Catton

Modeling of fin-and-tube heat exchangers based on the volume averaging theory (VAT) requires proper closure of the VAT based governing equations. Closure can be obtained from reasonable lower scale solutions of a computational fluid dynamics (CFD) code, which means the tube row number chosen should be large enough, so that the closure can be evaluated for a representative elementary volume (REV) that is, not affected by the entrance or recirculation at the outlet of the fin gap. To determine the number of tube rows, three-dimensional numerical simulations for plate fin-and-tube heat exchangers were performed, with the Reynolds number varying from 500 to 6000 and the number of tube rows varying from 1 to 9. A clear perspective of the variations of both overall and local fiction factor and the Nusselt number as the tube row number increases are presented. These variation trends are explained from the view point of the field synergy principle (FSP). Our investigation shows that 4 + 1 + 1 tube rows is the minimum number to get reasonable lower scale solutions. A computational domain including 5 + 2 + 2 tube rows is recommended, so that the closure formulas for drag resistance coefficient and heat transfer coefficient could be evaluated for the sixth and seventh elementary volumes to close the VAT based model.


Author(s):  
Jonathan Timo Lipkowicz ◽  
Damien Nativel ◽  
Sean Cooper ◽  
Irenäus Wlokas ◽  
Mustapha Fikri ◽  
...  

Abstract Highly resolved two- and three-dimensional computational fluid dynamics (CFD) simulations are presented for shock-tube experiments containing hydrogen/oxygen (H2/O2) mixtures, to investigate mechanisms leading to remote ignition. The results of the reactive cases are compared against experimental results from Meyer and Oppenheim (Proc Combust Inst 13(1): 1153–1164, 1971. 10.1016/s0082-0784(71)80112-1) and Hanson et al. (Combust Flame 160(9): 1550–1558, 2013. 10.1016/j.combustflame.2013.03.026). The results of the non-reactive case are compared against shock tube experiments, recently carried out in Duisburg and Texas. The computational domain covers the end-wall region of the shock tube and applies high order numerics featuring an all-speed approximate Riemann scheme, combined with a 5th order interpolation scheme. Direct chemistry is employed using detailed reaction mechanisms with 11 species and up to 40 reactions, on a grid with up to 2.2 billion cells. Additional two-dimensional simulations are performed for non-reactive conditions to validate the treatment of boundary-layer effects at the inlet of the computational domain. The computational domain covers a region at the end part of the shock tube. The ignition process is analyzed by fields of localized, expected ignition times. Instantaneous fields of temperature, pressure, entropy, and dissipation rate are presented to explain the flow dynamics, specifically in the case of a bifurcated reflected shock. In all cases regions with locally increased temperatures were observed, reducing the local ignition-delay time in areas away from the end wall significantly, thus compensating for the late compression by the reflected shock and therefore leading for first ignition at a remote location, i.e., away from the end wall where the ignition would occur under ideal conditions. In cases without a bifurcated reflected shock, the temperature increase results from shock attenuation. In cases with a bifurcated reflected shock, the formation of a second normal shock and shear near the slip line is found to be crucial for the remote ignition to take place. Overall, the two- and three-dimensional simulations were found to qualitatively explain the occurrence of remote ignition and to be quantitatively correct, implying that they include the correct physics.


Author(s):  
Kuen-Bae Lee ◽  
Mark Wilson ◽  
Mehdi Vahdati

Over recent years engine designs have moved increasingly toward low specific thrust cycles to deliver significant specific fuel consumption (SFC) improvements. Such fan blades are more loaded than conventional fan blades and therefore can be more prone to aerodynamic and aeroelastic instabilities. The aim of this paper is to analyse the flutter stability of a low speed/low pressure ratio fan blade. By using a validated CFD model (AU3D), three dimensional unsteady simulations are performed for a modern low speed fan rig for which extensive measured data are available. The computational domain contains a complete fan assembly with an intake duct and the downstream OGVs (whole LP domain). Flutter simulations are conducted over a range of speeds to understand flutter characteristics of this blade. Only the 1F (first flap) mode is considered in this work. Measured rig data obtained by using the same fan set but with two different intakes showed a significant difference in the flutter boundary for the two intakes. AU3D computations were performed for both intakes and were used to explain this difference between the two intakes, and showed that intake reflections play an important role in flutter of this blade. In the next phase of this work, two possible modifications for increasing the flutter margin of the fan blade were explored: 1. Changing the mode shape of the blade 2. Using acoustic liners in the casing The results show that it is possible to increase the flutter margin of the blade by either decreasing the ratio of the twisting to plunging motion in 1F mode or by introducing deep acoustic liners in the intake. The liners have to be deep enough to attenuate the flutter pressure waves and hence influence the stability. The results indicate the importance of reflection in flutter stability of the fan blade, and clearly show that intake duct needs to be included in flutter study of any fan blade.


Author(s):  
John P. Borg

In this paper the numerically determined pressure increase (ΔP) versus volumetric flow rate (CFM) curve at a given fan speed for a plenum fan is compared to experimental data. The simulations were carried out using a rotating fan blade section and a fixed inlet and outlet plenum with a three-dimensional tetrahedral computational mesh. The objective of this work is to assess the feasibility of approximating a 36” plenum ΔP-CFM fan curve using a low order computational approach and thereby assessing the effectiveness of using such an approach as a real time design tool. The measures of success of this work include demonstrating the ability to capture pertinent characteristics of the fan curve such as slope and roll-off of the ΔP-CFM curve. It was found that a fairly high resolution was required near the fan blade section in order to better approximate the ΔP-CFM curve. This higher resolution greatly increased the runtime. In addition, including a k-e turbulent model improved the pressure drop characteristics as compared having no turbulence model. It is hoped that an approach such as this will be adopted in the real time design and manufacture of plenum fans.


Author(s):  
Tamas Szili-Torok ◽  
Jens Rump ◽  
Torsten Luther ◽  
Sing-Chien Yap

Abstract Better understanding of the lead curvature, movement and their spatial distribution may be beneficial in developing lead testing methods, guiding implantations and improving life expectancy of implanted leads. Objective The aim of this two-phase study was to develop and test a novel biplane cine-fluoroscopy-based method to evaluate input parameters for bending stress in leads based on their in vivo 3D motion using precisely determined spatial distributions of lead curvatures. Potential tensile, compressive or torque forces were not subjects of this study. Methods A method to measure lead curvature and curvature evolution was initially tested in a phantom study. In the second phase using this model 51 patients with implanted ICD leads were included. A biplane cine-fluoroscopy recording of the intracardiac region of the lead was performed. The lead centerline and its motion were reconstructed in 3D and used to define lead curvature and curvature changes. The maximum absolute curvature Cmax during a cardiac cycle, the maximum curvature amplitude Camp and the maximum curvature Cmax@amp at the location of Camp were calculated. These parameters can be used to characterize fatigue stress in a lead under cyclical bending. Results The medians of Camp and Cmax@amp were 0.18 cm−1 and 0.42 cm−1, respectively. The median location of Cmax was in the atrium whereas the median location of Camp occurred close to where the transit through the tricuspid valve can be assumed. Increased curvatures were found for higher slack grades. Conclusion Our results suggest that reconstruction of 3D ICD lead motion is feasible using biplane cine-fluoroscopy. Lead curvatures can be computed with high accuracy and the results can be implemented to improve lead design and testing.


Vibration ◽  
2020 ◽  
Vol 4 (1) ◽  
pp. 49-63
Author(s):  
Waad Subber ◽  
Sayan Ghosh ◽  
Piyush Pandita ◽  
Yiming Zhang ◽  
Liping Wang

Industrial dynamical systems often exhibit multi-scale responses due to material heterogeneity and complex operation conditions. The smallest length-scale of the systems dynamics controls the numerical resolution required to resolve the embedded physics. In practice however, high numerical resolution is only required in a confined region of the domain where fast dynamics or localized material variability is exhibited, whereas a coarser discretization can be sufficient in the rest majority of the domain. Partitioning the complex dynamical system into smaller easier-to-solve problems based on the localized dynamics and material variability can reduce the overall computational cost. The region of interest can be specified based on the localized features of the solution, user interest, and correlation length of the material properties. For problems where a region of interest is not evident, Bayesian inference can provide a feasible solution. In this work, we employ a Bayesian framework to update the prior knowledge of the localized region of interest using measurements of the system response. Once, the region of interest is identified, the localized uncertainty is propagate forward through the computational domain. We demonstrate our framework using numerical experiments on a three-dimensional elastodynamic problem.


2019 ◽  
Vol 2019 ◽  
pp. 1-14
Author(s):  
Gang Xu ◽  
Guangwei Zhao ◽  
Jing Chen ◽  
Shuqi Wang ◽  
Weichao Shi

The value of the tangential velocity on the Boundary Value Problem (BVP) is inaccurate when comparing the results with analytical solutions by Indirect Boundary Element Method (IBEM), especially at the intersection region where the normal vector is changing rapidly (named nonsmooth boundary). In this study, the singularity of the BVP, which is directly arranged in the center of the surface of the fluid computing domain, is moved outside the computational domain by using the Desingularized Boundary Integral Equation Method (DBIEM). In order to analyze the accuracy of the IBEM/DBIEM and validate the above-mentioned problem, three-dimensional uniform flow over a sphere has been presented. The convergent study of the presented model has been investigated, including desingularized distance in the DBIEM. Then, the numerical results were compared with the analytical solution. It was found that the accuracy of velocity distribution in the flow field has been greatly improved at the intersection region, which has suddenly changed the boundary surface shape of the fluid domain. The conclusions can guide the study on the flow over nonsmooth boundaries by using boundary value method.


Author(s):  
R. C. Schlaps ◽  
S. Shahpar ◽  
V. Gümmer

In order to increase the performance of a modern gas turbine, compressors are required to provide higher pressure ratio and avoid incurring higher losses. The tandem aerofoil has the potential to achieve a higher blade loading in combination with lower losses compared to single vanes. The main reason for this is due to the fact that a new boundary layer is generated on the second blade surface and the turning can be achieved with smaller separation occurring. The lift split between the two vanes with respect to the overall turning is an important design choice. In this paper an automated three-dimensional optimisation of a highly loaded compressor stator is presented. For optimisation a novel methodology based on the Multipoint Approximation Method (MAM) is used. MAM makes use of an automatic design of experiments, response surface modelling and a trust region to represent the design space. The CFD solutions are obtained with the high-fidelity 3D Navier-Stokes solver HYDRA. In order to increase the stage performance the 3D shape of the tandem vane is modified changing both the front and rear aerofoils. Moreover the relative location of the two aerofoils is controlled modifying the axial and tangential relative positions. It is shown that the novel optimisation methodology is able to cope with a large number of design parameters and produce designs which performs better than its single vane counterpart in terms of efficiency and numerical stall margin. One of the key challenges in producing an automatic optimisation process has been the automatic generation of high-fidelity computational meshes. The multi block-structured, high-fidelity meshing tool PADRAM is enhanced to cope with the tandem blade topologies. The wakes of each aerofoil is properly resolved and the interaction and the mixing of the front aerofoil wake and the second tandem vane are adequately resolved.


Sign in / Sign up

Export Citation Format

Share Document