Fully Parametric High-Fidelity CFD Model for the Design Optimisation of the Cyclic Stagger Pattern of a Set of Fan Outlet Guide Vanes

Author(s):  
Andrea Milli ◽  
Olivier Bron

The present paper deals with the redesign of cyclic variation of a set of fan outlet guide vanes by means of high-fidelity full-annulus CFD. The necessity for the aerodynamic redesign originated from a change to the original project requirement, when the customer requested an increase in specific thrust above the original engine specification. The main objectives of this paper are: 1) make use of 3D CFD simulations to accurately model the flow field and identify high-loss regions; 2) elaborate an effective optimisation strategy using engineering judgement in order to define realistic objectives, constraints and design variables; 3) emphasise the importance of parametric geometry modelling and meshing for automatic design optimisation of complex turbomachinery configurations; 4) illustrate that the combination of advanced optimisation algorithms and aerodynamic expertise can lead to successful optimisations of complex turbomachinery components within practical time and costs constrains. The current design optimisation exercise was carried out using an in-house set of software tools to mesh, resolve, analyse and optimise turbomachinery components by means of Reynolds-averaged Navier-Stokes simulations. The original configuration was analysed using the 3D CFD model and thereafter assessed against experimental data and flow visualisations. The main objective of this phase was to acquire a deep insight of the aerodynamics and the loss mechanisms. This was important to appropriately limit the design scope and to drive the optimisation in the desirable direction with a limited number of design variables. A mesh sensitivity study was performed in order to minimise computational costs. Partially converged CFD solutions with restart and response surface models were used to speed up the optimisation loop. Finally, the single-point optimised circumferential stagger pattern was manually adjusted to increase the robustness of the design at other flight operating conditions. Overall, the optimisation resulted in a major loss reduction and increased operating range. Most important, it provided the project with an alternative and improved design within the time schedule requested and demonstrated that CFD tools can be used effectively not only for the analysis but also to provide new design solutions as a matter of routine even for very complex geometry configurations.

2015 ◽  
Vol 59 (04) ◽  
pp. 209-226
Author(s):  
Nitin Garg ◽  
Gaetan K.W. Kenway ◽  
Zhoujie Lyu ◽  
Joaquim R.R.A. Martins ◽  
Yin L. Young

With recent advances in high-performance computing, computational fluid dynamics (CFD) modeling has become an integral part in the engineering analysis and even in the design process of marine vessels and propulsors. In aircraft wing design, CFD has been integrated with numerical optimization and adjoint methods to enable high fidelity aerodynamic shape optimization with respect to large numbers of design variables. There is a potential to use some of these techniques for maritime applications, but there are new challenges that need to be addressed to realize that potential. This work presents a solution to some of those challenges by developing a CFD-based hydrodynamic shape optimization tool that considers cavitation and a wide range of operating conditions. A previously developed three-dimensional compressible Reynold saveraged Navier-Stokes (RANS) solver is extended to solve for nearly incompressible flows, using a low-speed preconditioner. An efficient gradient-based optimizer and the adjoint method are used to carry out the optimization. The modified CFD solver is validated and verified for a tapered NACA 0009 hydrofoil. The need for a large number of design variables is demonstrated by comparing the optimized solution obtained using different number of shape design variables. The results showed that at least 200 design variables are needed to get a converged optimal solution for the hydrofoil considered. The need for a high-fidelity hydrodynamic optimization tool is also demonstrated by comparing RANS-based optimization with Euler-based optimization. The results show that at high lift coefficient (CL) values, the Euler-based optimization leads to a geometry that cannot meet the required lift at the same angle of attack as the original foil due to inability of the Euler solver to predict viscous effects. Single-point optimization studies are conducted for various target CL values and compared with the geometry and performance of the original NACA 0009 hydrofoil, as well as with the results from a multipoint optimization study. A total of 210 design variables are used in the optimization studies. The optimized foil is found to have a much lower negative suction peak, and hence delayed cavitation inception, in addition to higher efficiency, compared to the original foil at the design CL value. The results show significantly different optimal geometry for each CL, which means an active morphing capability was needed to achieve the best possible performance for all conditions. For the single-point optimization, using the highest CL as the design point, the optimized foil yielded the best performance at the design point, but the performance degraded at the off-design CL points compared to the multipoint design. In particular, the foil optimized for the highest CL showed inferior performance even compared to the original foil at the lowest CL condition. On the other hand, the multipoint optimized hydrofoil was found to perform better than the original NACA 0009 hydrofoil over the entire operation profile, where the overall efficiency weighted by the probability of operation at each CL, is improved by 14.4%. For the multipoint optimized foil, the geometry remains fixed throughout the operation profile and the overall efficiency was only 1.5% lower than the hypothetical actively morphed foil with the optimal geometry at each CL. The new methodology presented herein has the potential to improve the design of hydrodynamic lifting surfaces such as propulsors, hydrofoils, and hulls.


Author(s):  
K. R. Mrinal ◽  
Md. Hamid Siddique ◽  
Abdus Samad

A progressive cavity pump (PCP) is a positive displacement pump and has been used as an artificial lift method in the oil and gas industry for pumping fluid with solid content and high viscosity. In a PCP, a single-lobe rotor rotates inside a double-lobe stator. Articles on computational works for flows through a PCP are limited because of transient behavior of flow, complex geometry and moving boundaries. In this paper, a 3D CFD model has been developed to predict the flow variables at different operating conditions. The flow is considered as incompressible, single phase, transient, and turbulent. The dynamic mesh model in Ansys-Fluent for the rotor mesh movement is used, and a user defined function (UDF) written in C language defines the rotor’s hypocycloid path. The mesh deformation is done with spring based smoothing and local remeshing technique. The computational results are compared with the experiment results available in the literature. Thepump gives maximum flowrate at zero differential pressure.


Author(s):  
Jin Yan ◽  
Andrew Mallner

Boiling water reactors (BWRs) are equipped with a standby liquid control system (SLCS). The SLCS is used to inject boron to shutdown the reactor from full power condition in the event that the control rods fail to insert. In order for the SLCS system to shut down the reactor, adequate mixing of the borated solution with the reactor coolant is necessary. In BWRs prior to BWR 5, the boron injection points are located in the lower plenum. The objective of this project is to evaluate the impact of the operating conditions on the boron injection based on the understanding of the behavior of multi-species flow in typical pre-BWR 5 reactors by using computational fluid dynamics (CFD). The project is divided into two phases. At the first phase, a CFD model based on the test configurations of GE 1/6 scale test program was established. The results were validated against measurements conducted by GE during the 1/6 scale test program performed in 1981. The validation shows that the CFD can give accurate predictions of the boron mixing. The technical approach employed in the CFD model was adequate to capture the boron mixing process in the BWR lower plenum. The second phase of the project is the sensitivity study based on the same technical approach developed in the first phase. However, a simplified BWR lower plenum model was used due to the time constrain. In the sensitivity study, the baseline case and four additional cases with different operating conditions were investigated using the same CFD approach.


2015 ◽  
Vol 137 (9) ◽  
Author(s):  
Hadrien Montanelli ◽  
Marc Montagnac ◽  
François Gallard

This paper presents the application of the gradient span analysis (GSA) method to the multipoint optimization of the two-dimensional LS89 turbine distributor. The cost function (total pressure loss) and the constraint (mass flow rate) are computed from the resolution of the Reynolds-averaged Navier–Stokes equations. The penalty method is used to replace the constrained optimization problem with an unconstrained problem. The optimization process is steered by a gradient-based quasi-Newton algorithm. The gradient of the cost function with respect to design variables is obtained with the discrete adjoint method, which ensures an efficient computation time independent of the number of design variables. The GSA method gives a minimal set of operating conditions to insert into the weighted sum model to solve the multipoint optimization problem. The weights associated to these conditions are computed with the utopia point method. The single-point optimization at the nominal condition and the multipoint optimization over a wide range of conditions of the LS89 blade are compared. The comparison shows the strong advantages of the multipoint optimization with the GSA method and utopia-point weighting over the traditional single-point optimization.


Author(s):  
Alexander O. Pugachev

The study deals with optimization of leakage and rotordynamic characteristics of liquid annular seals. A nonlinear constrained multi-objective optimization problem is considered. An objective function is a weighted sum of leakage and whirl-frequency ratio of the seal. Side constraints are imposed on design variables. The seal consists of two rings which shape can be either cylinder or converging taper or diverging taper. There are four design variables — seal length and diameters of the rings. A non-gradient-based method is used for solving the optimization problem. Analysis of the seal performance is based on computational fluid dynamics (CFD). A full 3D eccentric CFD model of the seal including upstream and downstream regions is constructed in ANSYS CFX. The solution procedures for prediction of rotordynamic coefficients are discussed and compared. The whirling rotor method under the assumption of centered circular orbit is used in optimization runs. The CFD model of the seal is validated against experimental data taken from the literature. A mesh independence study is carried out. An optimization environment includes automatic grid generation, parallel CFD calculations of the seal, and optimization algorithm. Two optimization runs corresponding to low-speed and high-speed cases are performed. Seals with improved characteristics include near-cylindrical and divergent-tapered rings. Performance of three seals from the Pareto set is calculated for different rotational speeds and inlet pressures. Generally, the rotordynamic performance degrades at other operating conditions. Additional study for the seals with enlarged clearance is carried out to model effect of wear.


Author(s):  
A. G. Sanvito ◽  
G. Persico ◽  
M. S. Campobasso

Abstract This study provides a novel contribution towards the establishment of a new high–fidelity simulation–based design methodology for stall–regulated horizontal axis wind turbines. The aerodynamic design of these machines is complex, due to the difficulty of reliably predicting stall onset and post–stall characteristics. Low–fidelity design methods, widely used in industry, are computationally efficient, but are often affected by significant uncertainty. Conversely, Navier–Stokes CFD can reduce such uncertainty, resulting in lower development costs by reducing the need of field testing of designs not fit for purpose. Here, the compressible CFD research code COSA is used to assess the performance of two alternative designs of a 13–meter stall–regulated rotor over a wide range of operating conditions. Validation of the numerical methodology is based on thorough comparisons of novel simulations and measured data of the NREL Phase VI turbine rotor, and one of the two industrial rotor designs. An excellent agreement is found in all cases. All simulations of the two industrial rotors are time–dependent, to capture the unsteadiness associated with stall which occurs at most wind speeds. The two designs are cross-compared, with emphasis on the different stall patterns resulting from particular design choices. The key novelty of this work is the CFD–based assessment of the correlation among turbine power, blade aerodynamics, and blade design variables (airfoil geometry, blade planform and twist) over most operational wind speeds.


Author(s):  
Vassilios Pachidis ◽  
Pericles Pilidis ◽  
Geoffrey Guindeuil ◽  
Anestis Kalfas ◽  
Ioannis Templalexis

This study focuses on a simulation strategy that will allow the performance characteristics of an isolated gas turbine engine component, resolved from a detailed, high-fidelity analysis, to be transferred to an engine system analysis carried out at a lower level of resolution. This work will enable component-level, complex physical processes to be captured and analyzed in the context of the whole engine performance, at an affordable computing resource and time. The technique described in this paper utilizes an object-oriented, zero-dimensional (0-D) gas turbine modeling and performance simulation system and a high-fidelity, three-dimensional (3-D) computational fluid dynamics (CFD) component model. The technique is called ‘partially integrated’ zooming, in that there is no automatic link between the 0-D engine cycle and the 3-D CFD model. It can be applied to all engine components and involves the generation of a component characteristic map via an iterative execution of the 0-D cycle and the 3-D CFD model. This work investigates relative changes in the simulated engine performance after integrating the CFD-generated component map into the 0-D engine analysis. This paper attempts to demonstrate the ‘partially integrated’ approach to component zooming by using a 3-D CFD intake model of a high by-pass ratio (HBR) turbofan as a case study. The CFD model is based on the geometry of the intake of the CFM56-5B2 engine. The CFD-generated performance map can fully define the characteristic of the intake at several operating conditions and is subsequently used to provide a more accurate, physics-based estimate of intake performance (i.e. pressure recovery) and hence, engine performance, replacing the default, empirical values within the 0-D cycle model. A detailed comparison between the baseline engine performance (empirical pressure recovery) and the engine performance obtained after using the CFD-generated map is presented in this paper. The analysis carried out by this study, demonstrates relative changes in the simulated engine performance larger than 1%.


2021 ◽  
Vol 11 (12) ◽  
pp. 5430
Author(s):  
Paolo Neri ◽  
Alessandro Paoli ◽  
Ciro Santus

Vibration measurements of turbomachinery components are of utmost importance to characterize the dynamic behavior of rotating machines, thus preventing undesired operating conditions. Local techniques such as strain gauges or laser Doppler vibrometers are usually adopted to collect vibration data. However, these approaches provide single-point and generally 1D measurements. The present work proposes an optical technique, which uses two low-speed cameras, a multimedia projector, and three-dimensional digital image correlation (3D-DIC) to provide full-field measurements of a bladed disk undergoing harmonic response analysis (i.e., pure sinusoidal excitation) in the kHz range. The proposed approach exploits a downsampling strategy to overcome the limitations introduced by low-speed cameras. The developed experimental setup was used to measure the response of a bladed disk subjected to an excitation frequency above 6 kHz, providing a deep insight in the deformed shapes, in terms of amplitude and phase distributions, which could not be feasible with single-point sensors. Results demonstrated the system’s effectiveness in measuring amplitudes of few microns, also evidencing blade mistuning effects. A deeper insight into the deformed shape analysis was provided by considering the phase maps on the entire blisk geometry, and phase variation lines were observed on the blades for high excitation frequency.


Water ◽  
2021 ◽  
Vol 13 (4) ◽  
pp. 463
Author(s):  
Gopinathan R. Abhijith ◽  
Leonid Kadinski ◽  
Avi Ostfeld

The formation of bacterial regrowth and disinfection by-products is ubiquitous in chlorinated water distribution systems (WDSs) operated with organic loads. A generic, easy-to-use mechanistic model describing the fundamental processes governing the interrelationship between chlorine, total organic carbon (TOC), and bacteria to analyze the spatiotemporal water quality variations in WDSs was developed using EPANET-MSX. The representation of multispecies reactions was simplified to minimize the interdependent model parameters. The physicochemical/biological processes that cannot be experimentally determined were neglected. The effects of source water characteristics and water residence time on controlling bacterial regrowth and Trihalomethane (THM) formation in two well-tested systems under chlorinated and non-chlorinated conditions were analyzed by applying the model. The results established that a 100% increase in the free chlorine concentration and a 50% reduction in the TOC at the source effectuated a 5.87 log scale decrement in the bacteriological activity at the expense of a 60% increase in THM formation. The sensitivity study showed the impact of the operating conditions and the network characteristics in determining parameter sensitivities to model outputs. The maximum specific growth rate constant for bulk phase bacteria was found to be the most sensitive parameter to the predicted bacterial regrowth.


Author(s):  
Zumrat Usmanova ◽  
Emin Sunbuloglu

Numerical simulation of automotive tires is still a challenging problem due to their complex geometry and structures, as well as the non-uniform loading and operating conditions. Hysteretic loss and rolling resistance are the most crucial features of tire design for engineers. A decoupled numerical model was proposed to predict hysteretic loss and temperature distribution in a tire, however temperature dependent material properties being utilized only during the heat generation analysis stage. Cyclic change of strain energy values was extracted from 3-D deformation analysis, which was further used in a thermal analysis as input to predict temperature distribution and thermal heat generation due to hysteretic loss. This method was compared with the decoupled model where temperature dependence was ignored in both deformation and thermal analysis stages. Deformation analysis results were compared with experimental data available. The proposed method of numerical modeling was quite accurate and results were found to be close to the actual tire behavior. It was shown that one-way-coupled method provides rolling resistance and peak temperature values that are in agreement with experimental values as well.


Sign in / Sign up

Export Citation Format

Share Document