Engineering topics

Rigby, David L.

Publications and source records attributed to Rigby, David L..

At least 19 records

NASA NTRS2017

Viscous Three-Dimensional Simulation of Flow in an Axial Low Pressure Compressor at Engine Icing Operating Points

Viscous three-dimensional simulations of the Honeywell ALF502R-5 low pressure compressor (sometimes called a booster) using the NASA Glenn code GlennHT have been carried out. A total of ten simulations were produced. Five operating points are investigated, with each point run with two different wall thermal conditions. These operating points are at, or near, points where engine icing has been determined to be likely. In the future, the results of this study will be used for further analysis such as predicting collection efficiency of ice particles and ice growth rates at various locations in the compressor. A mixing plane boundary condition is used between each blade row, resulting in convergence to steady state within each blade row. The k-omega turbulence model of Wilcox, combined with viscous grid spacing near the wall on the order of one, is used to resolve the turbulent boundary layers. For each of the operating points, heat transfer coefficients are generated on the blades and walls. The heat transfer coefficients are produced by running the operating point with two different wall thermal conditions and then solving simultaneously for the heat transfer coefficient and adiabatic wall temperature at each point. Average Nusselt numbers are calculated for the most relevant surfaces. The values are seen to scale with Reynolds number to approximately a power of 0.7. Additionally, images of surface distribution of Nusselt number are presented. Qualitative comparison between the five operating points show that there is relatively little change in the character of the distribution. The dominant observed effect is that of an overall scaling, which is expected due to Reynolds number differences. One interesting aspect about the Nusselt number distribution is observed on the casing (outer diameter) downstream of the exit guide vanes (EGVs). The Nusselt number is relatively high between the pairs of EGVs, with two lower troughs downstream of each EGV trailing edge. This is of particular interest since rather complex ice shapes have been observed in that region.

compressor↗

NASA NTRS2015

LEWICE3D/GlennHT Particle Analysis of the Honeywell Al502 Low Pressure Compressor

A flow and ice particle trajectory analysis was performed for the booster of the Honeywell AL502 engine. The analysis focused on two closely related conditions one of which produced a rollback and another which did not rollback during testing in the Propulsion Systems Lab at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.56 ice accretion software. The flow and particle analysis used a 3D steady flow, mixing plane approach to model the transport of flow and particles through the engine. The inflow conditions for the rollback case were: airspeed, 145 ms; static pressure, 33,373 Pa; static temperature, 253.3 K. The inflow conditions for the non-roll-back case were: airspeed, 153 ms; static pressure, 34,252 Pa; static temperature, 260.1 K. Both cases were subjected to an ice particle cloud with a median volume diameter of 24 microns, an ice water content of 2.0 gm3 and a relative humidity of 100 percent. The most significant difference between the rollback and non-rollback conditions was the inflow static temperature which was 6.8 K higher for the non-rollback case.

Ice Accretion Codes↗

NASA NTRS2015

Ice Particle Analysis of the Honeywell AL502 Engine Booster

A flow and ice particle trajectory analysis was performed for the booster of the Honeywell ALF502 engine. The analysis focused on two closely related conditions one of which produced an icing event and another which did not during testing of the ALF502 engine in the Propulsion Systems Lab (PSL) at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.63 ice accretion software. The inflow conditions for the two conditions were similar with the main differences being that the condition that produced the icing event was 6.8 K colder than the non-icing event case and the inflow ice water content (IWC) for the non-icing event case was 50% less than for the icing event case. The particle analysis, which considered sublimation, evaporation and phase change, was generated for a 5 micron ice particle with a sticky impact model and for a 24 micron median volume diameter (MVD), 7 bin ice particle distribution with a supercooled large droplet (SLD) splash model used to simulate ice particle breakup. The particle analysis did not consider the effect of the runback and re-impingement of water resulting from the heated spinner and anti-icing system. The results from the analysis showed that the amount of impingement for the components were similar for the same particle size and impact model for the icing and non-icing event conditions. This was attributed to the similar aerodynamic conditions in the booster for the two cases. The particle temperature and melt fraction were higher at the same location and particle size for the non-icing event than for the icing event case due to the higher incoming inflow temperature for the non-event case. The 5 micron ice particle case produced higher impact temperatures and higher melt fractions on the components downstream of the fan than the 24 micron MVD case because the average particle size generated by the particle breakup was larger than 5 microns which yielded less warming and melting. The analysis also showed that the melt fraction and wet bulb temperature icing criterion developed during tests in the Research Altitude Test Facility (RATFac) at the National Research Council (NRC) of Canada were useful in predicting icing events in the ALF502 engine. The development of an ice particle impact model which includes the effects of particle breakup, phase change, and surface state is necessary to further improve the prediction of ice particle transport with phase change through turbomachinery.

Ice Accretion Codes↗

NASA NTRS2014

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Aircraft Icing↗

NASA NTRS2014

A Thermal Analysis of a Hot-Wire Probe for Icing Applications

This paper presents a steady-state thermal model of a hot-wire instrument applicable to atmospheric measurement of water content in clouds. In this application, the power required to maintain the wire at a given temperature is used to deduce the water content of the cloud. The model considers electrical resistive heating, axial conduction, convection to the flow, radiation to the surroundings, as well as energy loss due to the heating, melting, and evaporation of impinging liquid and or ice. All of these parameters can be varied axially along the wire. The model further introduces a parameter called the evaporation potential which locally gauges the maximum fraction of incoming water that evaporates. The primary outputs of the model are the steady-state power required to maintain a spatially-average constant temperature as well as the variation of that temperature and other parameters along the wire. The model is used to understand the sensitivity of the hot-wire performance to various flow and boundary conditions including a detailed comparison of dry air and wet (i.e. cloud-on) conditions. The steady-state power values are compared to experimental results from a Science Engineering Associates (SEA) Multi-Element probe, a commonly used water-content measurement instrument. The model results show good agreement with experiment for both dry and cloud-on conditions with liquid water content. For ice, the experimental measurements under read the actual water content due to incomplete evaporation and splashing. Model results, which account for incomplete evaporation, are still higher than experimental results where the discrepancy is attributed to splashing mass-loss which is not accounted in the model.

Thermal Analysis↗

NASA NTRS2014

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Water Content↗

NASA NTRS2013

Multifunctional Low-Pressure Turbine for Core Noise Reduction, Improved Efficiency, and Nitrogen Oxide (NOx) Reduction

This work studied the feasibility of using Helmholtz resonator cavities embedded in low-pressure-turbine (LPT) airfoils to (1) reduce core noise by damping acoustic modes; (2) use the synthetic jets produced by the liner hole acoustic oscillations to improve engine efficiency by maintaining turbulent attached flow in the LPT at low-Reynolds-number cruise conditions; and (3) reduce engine nitrogen oxide emissions by lining the internal cavities with materials capable of catalytic conversion. Flat plates with embedded Helmholtz resonators, designed to resonate at either 3000 or at 400 Hz, were simulated using computational fluid dynamics. The simulations were conducted for two inlet Mach numbers, 0.25 and 0.5, corresponding to Reynolds numbers of 90 000 and 164 000 based on the effective chordwise distance to the resonator orifice. The results of this study are (1) the region of acoustic treatment may be large enough to have a benefit; (2) the jets may not possess sufficient strength to reduce flow separation (based on prior work by researchers in the flow control area); and (3) the additional catalytic surface area is not exposed to a high velocity, so it probably does not have any benefit.

Miller, Christopher J.↗

NASA NTRS2008

Unsteady Analysis of Blade and Tip Heat Transfer as Influenced by the Upstream Momentum and Thermal Wakes

The effect of the upstream wake on the blade heat transfer has been numerically examined. The geometry and the flow conditions of the first stage turbine blade of GE s E3 engine with a tip clearance equal to 2 percent of the span was utilized. Based on numerical calculations of the vane, a set of wake boundary conditions were approximated, which were subsequently imposed upon the downstream blade. This set consisted of the momentum and thermal wakes as well as the variation in modeled turbulence quantities of turbulence intensity and the length scale. Using a one-blade periodic domain, the distributions of unsteady heat transfer rate on the turbine blade and its tip, as affected by the wake, were determined. Such heat transfer coefficient distribution was computed using the wall heat flux and the adiabatic wall temperature to desensitize the heat transfer coefficient to the wall temperature. For the determination of the wall heat flux and the adiabatic wall temperatures, two sets of computations were required. The results were used in a phase-locked manner to compute the unsteady or steady heat transfer coefficients. It has been found that the unsteady wake has some effect on the distribution of the time averaged heat transfer coefficient on the blade and that this distribution is different from the distribution that is obtainable from a steady computation. This difference was found to be as large as 20 percent of the average heat transfer on the blade surface. On the tip surface, this difference is comparatively smaller and can be as large as four percent of the average.

Ameri, Ali A.↗

NASA NTRS2007

Unsteady Turbine Blade and Tip Heat Transfer Due to Wake Passing

The geometry and the flow conditions of the first stage turbine blade of GE s E3 engine have been used to obtain the unsteady three-dimensional blade and tip heat transfer. The isothermal wall boundary condition was used. The effect of the upstream wake of the first stage vane was of interest and was simulated by provision of a gust type boundary condition upstream of the blades. A one blade periodic domain was used. The consequence of this choice was explored in a preliminary study which showed little difference in the time mean heat transfer between 1:1 and 2:3 vane/blade domains. The full three-dimensional computations are of the blade having a clearance gap of 2 percent the span. Comparison between the time averaged unsteady and steady heat transfer is provided. It is shown that there is a significant difference between the steady and time mean of unsteady blade heat transfer in localized regions. The differences on the suction side of the blade in the near hub and near tip regions were found to be rather significant. Steady analysis underestimated the blade heat transfer by as much as 20 percent as compared to the time average obtained from the unsteady analysis. As for the blade tip, the steady analysis and the unsteady analysis gave results to within 2 percent.

Ameri, Ali A.↗

NASA NTRS2007

Effects of Unsteadiness Due to Wake Passing on Rotor Blade Heat Transfer

14. ABSTRACT In a gas turbine engine, the turbine rotor blades are buffeted by the wakes of the vanes located upstream. There is a transient effect from the passing of wakes on the blade heat transfer. This transient effect has been computed for a representative rotor by introducing a wake upstream via an unsteady inlet flow boundary condition, or "gust" condition. Two cases of turbulent flow and laminar flow with Reynolds numbers of 385,000 and 385 respectively were considered. For the turbulent flow case a quasi-steady calculation was also performed. The variation in the unsteady heat transfer coefficient was found to be as high as 120 percent of the mean. For the turbulent flow case a quasisteady calculation was also performed. The time mean of the unsteady heat transfer, the mean of the quasi-steady variations and the steady results agree reasonably well on all blade locations except for the turbulent results which differ near the leading edge. The quasi-steady heat transfer results do not agree with the instantaneous unsteady results, although the time-mean values are similar.

Ameri, Ali A.↗

NASA NTRS2006

Code for Multiblock CFD and Heat-Transfer Computations

The NASA Glenn Research Center General Multi-Block Navier-Stokes Convective Heat Transfer Code, Glenn-HT, has been used extensively to predict heat transfer and fluid flow for a variety of steady gas turbine engine problems. Recently, the Glenn-HT code has been completely rewritten in Fortran 90/95, a more object-oriented language that allows programmers to create code that is more modular and makes more efficient use of data structures. The new implementation takes full advantage of the capabilities of the Fortran 90/95 programming language. As a result, the Glenn-HT code now provides dynamic memory allocation, modular design, and unsteady flow capability. This allows for the heat-transfer analysis of a full turbine stage. The code has been demonstrated for an unsteady inflow condition, and gridding efforts have been initiated for a full turbine stage unsteady calculation. This analysis will be the first to simultaneously include the effects of rotation, blade interaction, film cooling, and tip clearance with recessed tip on turbine heat transfer and cooling performance. Future plans call for the application of the new Glenn-HT code to a range of gas turbine engine problems of current interest to the heat-transfer community. The new unsteady flow capability will allow researchers to predict the effect of unsteady flow phenomena upon the convective heat transfer of turbine blades and vanes. Work will also continue on the development of conjugate heat-transfer capability in the code, where simultaneous solution of convective and conductive heat-transfer domains is accomplished. Finally, advanced turbulence and fluid flow models and automatic gridding techniques are being developed that will be applied to the Glenn-HT code and solution process.

Fabian, John C.↗

NASA NTRS2004

TopMaker: Technique Developed for Automatic Multiblock Topology Generation Using the Medial Axis

The TopMaker technique was developed in an effort to reduce the time required for grid generation in complex numerical studies. Topology generation accounts for much of the man-hours required for structured multiblock grids. With regard to structured multiblock grids, topology refers to how the blocks are arranged and connected. A two-dimensional multiblock topology generation technique has been developed at the NASA Glenn Research Center. Very general configurations can be addressed by the technique. A configuration is defined by a collection of non-intersecting closed curves, which will be referred to as loops. More than a single loop implies that holes exist in the domain, which poses no problem. This technique requires only the medial vertices and the touch points that define each vertex. From the information about the medial vertices, the connectivity between medial vertices is generated. The physical shape of the medial edge is not required. By applying a few simple rules to each medial edge, a multiblock topology can be generated without user intervention. The resulting topologies contain only the level of complexity dictated by the configurations. Grid lines remain attached to the boundary except at sharp concave turns, where a change in index family is introduced as would be desired. Keeping grid lines attached to the boundary is especially important in computational fluid dynamics, where highly clustered grids are used near no-slip boundaries. This technique is simple and robust and can easily be incorporated into the overall grid-generation process.

Rigby, David L.↗

NASA NTRS2004

TopMaker: A Technique for Automatic Multi-Block Topology Generation Using the Medial Axis

A two-dimensional multi-block topology generation technique has been developed. Very general configurations are addressable by the technique. A configuration is defined by a collection of non-intersecting closed curves, which will be referred to as loops. More than a single loop implies that holes exist in the domain, which poses no problem. This technique requires only the medial vertices and the touch points that define each vertex. From the information about the medial vertices, the connectivity between medial vertices is generated. The physical shape of the medial edge is not required. By applying a few simple rules to each medial edge, the multiblock topology is generated with no user intervention required. The resulting topologies contain only the level of complexity dictated by the configurations. Grid lines remain attached to the boundary except at sharp concave turns where a change in index family is introduced as would be desired. Keeping grid lines attached to the boundary is especially important in the area of computational fluid dynamics where highly clustered grids are used near no-slip boundaries. This technique is simple and robust and can easily be incorporated into the overall grid generation process.

Heidmann, James D.↗

NASA NTRS2003

Glenn-HT Code Validated for Complex Turbine Blade Cooling Passage

This work is motivated by the need to accurately predict heat transfer in turbomachinery. For efficient gas turbine operation, flow temperatures in the hot gas path exceed acceptable metal temperatures in many regions of the engine. So that the integrity of the parts can be maintained for an acceptable engine life, the parts must be cooled. Efficient cooling schemes require accurate heat transfer prediction to minimize regions that are overcooled and, even more importantly, to ensure adequate cooling in high-heat-flux regions.

Rigby, David L.↗

NASA NTRS2003

An Aerodynamic Simulation Process for Iced Lifting Surfaces and Associated Issues

This paper discusses technologies and software tools that are being implemented in a software toolkit currently under development at NASA Glenn Research Center. Its purpose is to help study the effects of icing on airfoil performance and assist with the aerodynamic simulation process which consists of characterization and modeling of ice geometry, application of block topology and grid generation, and flow simulation. Tools and technologies for each task have been carefully chosen based on their contribution to the overall process. For the geometry characterization and modeling, we have chosen an interactive rather than automatic process in order to handle numerous ice shapes. An Appendix presents features of a software toolkit developed to support the interactive process. Approaches taken for the generation of block topology and grids, and flow simulation, though not yet implemented in the software, are discussed with reasons for why particular methods are chosen. Some of the issues that need to be addressed and discussed by the icing community are also included.

Choo, Yung K.↗

NASA NTRS2002

Heat Transfer in a Complex Trailing Edge Passage for a High Pressure Turbine Blade: Simulation Results - Part 2:

A combined experimental and numerical study to investigate the heat transfer distribution in a complex blade trailing edge passage was conducted. The geometry consists of a two pass serpentine passage with taper toward the trailing edge, as well as from hub to tip. The upflow channel has an average aspect ratio of roughly 14:1, while the exit passage aspect ratio is about 5:1. The upflow channel is split in an interrupted way and is smooth on the trailing edge side of the split and turbulated on the other side. A turning vane is placed near the tip of the upflow channel. Reynolds numbers in the range of 31,000 to 61,000, based on inlet conditions, were simulated numerically. The simulation was performed using the Glenn-HT code, a full three-dimensional Navier-Stokes solver using the Wilcox k-omega turbulence model. A structured multi-block grid is used with approximately 4.5 million cells and average y+ values on the order of unity. Pressure and heat transfer distributions are presented with comparison to the experimental data. While there are some regions with discrepancies, in general the agreement is very good for both pressure and heat transfer.

Rigby, David L.↗

NASA NTRS2000

Heat Transfer in a Complex Trailing Edge Passage for a High Pressure Turbine Blade - Part 1: Experimental Measurements: Experimental Measurements - Part 1

A combined experimental and computational study has been performed to investigate the detailed heat transfer coefficient distributions within a complex blade trailing edge passage. The experimental measurements are made using a steady liquid crystal thermography technique applied to one major side of the passage. The geometry of the trailing edge passage is that of a two-pass serpentine circuit with a sharp 180-degree turning region at the tip. The upflow channel is split by interrupted ribs into two major subchannels, one of which is turbulated. This channel has an average aspect ratio of roughly 14:1. The spanwise extent of the channel geometry includes both area convergence from root to tip, as well as taper towards the trailing edge apex. The average section Reynolds numbers tested in this upflow channel range from 55,000 to 98,000. The tip section contains a turning vane near the extreme comer. The downflow channel has an aspect ratio of about 5:1, and also includes convergence and taper. Turbulators of varying sizes are included in this channel also. Both detailed heat transfer and pressure distribution measurements are presented. The pressure measurements are incorporated into a flow network model illustrating the major loss contributors.

Bunker, Ronald S.↗

NASA NTRS2000

Numerical Simulation of Non-rotating and Rotating Coolant Channel Flow Fields

Numerical simulations of the three-dimensional flow and heat transfer in a rectangular duct with a bend were performed. Results are presented for Reynolds numbers of 17,000 and 37,000 and for aspect ratios of 0.5 and 1.0. A k-(omega) turbulence model with no reference to distance to a wall is used. Direct comparison between single block and multiblock grid calculations are made. Heat transfer and velocity distributions are compared to available literature with good agreement. The multi-block grid system is seen to produce more accurate results compared to a single-block grid with the same number of cells.

Rigby, David L.↗