SEARCH · Engineering Papers

Results for “shock capturing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Grid adaptation using chimera composite overlapping meshes

The objective of this paper is to perform grid adaptation using composite overlapping meshes in regions of large gradient to accurately capture the salient features during computation. The chimera grid scheme, a multiple overset mesh technique, is used in combination with a Navier-Stokes solver. The numerical solution is first converged to a steady state based on an initial coarse mesh. Solution-adaptive enhancement is then performed by using a secondary fine grid system which oversets on top of the base grid in the high-gradient region, but without requiring the mesh boundaries to join in any special way. Communications through boundary interfaces between those separated grids are carried out using trilinear interpolation. Application to the Euler equations for shock reflections and to shock wave/boundary layer interaction problem are tested. With the present method, the salient features are well-resolved.

Kao, Kai-Hsiung↗

NASA NTRS2006

Aerothermodynamic Analyses of Towed Ballutes

A ballute (balloon-parachute) is an inflatable, aerodynamic drag device for application to planetary entry vehicles. Two challenging aspects of aerothermal simulation of towed ballutes are considered. The first challenge, simulation of a complete system including inflatable tethers and a trailing toroidal ballute, is addressed using the unstructured-grid, Navier-Stokes solver FUN3D. Auxiliary simulations of a semi-infinite cylinder using the rarefied flow, Direct Simulation Monte Carlo solver, DSV2, provide additional insight into limiting behavior of the aerothermal environment around tethers directly exposed to the free stream. Simulations reveal pressures higher than stagnation and corresponding large heating rates on the tether as it emerges from the spacecraft base flow and passes through the spacecraft bow shock. The footprint of the tether shock on the toroidal ballute is also subject to heating amplification. Design options to accommodate or reduce these environments are discussed. The second challenge addresses time-accurate simulation to detect the onset of unsteady flow interactions as a function of geometry and Reynolds number. Video of unsteady interactions measured in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel and CFD simulations using the structured grid, Navier-Stokes solver LAURA are compared for flow over a rigid spacecraft-sting-toroid system. The experimental data provides qualitative information on the amplitude and onset of unsteady motion which is captured in the numerical simulations. The presence of severe unsteady fluid - structure interactions is undesirable and numerical simulation must be able to predict the onset of such motion.

Gnoffo, Peter A.↗

NASA NTRS2009

Recovering Aerodynamic Side Loads on Rocket Nozzles using Quasi-Static Strain-Gage Measurements

During over-expanded operation of rocket nozzles, which is defined to be when the exit pressure is greater than internal pressure over some part of the nozzle, the nozzle will experience a transverse forcing function due to the pressure differential across the nozzle wall. Over-expansion occurs during the nozzle start-up and shutdown transient, even in high-altitude engines, because most test facilities cannot completely reproduce the near-vacuum pressures at those altitudes. During this transient, the pressure differential moves axially down the nozzle as it becomes pressurized, but this differential is never perfectly symmetric circumferentially. The character of the forcing function is highly complex and defined by a series of restricted and free shock separations. The subject of this paper is the determination of the magnitude of this loading during sub-scale testing via measurement of the structural dynamic response of the nozzle and its support structure. An initial attempt at back-calculating this load using the inverse of the transfer function was performed, but this attempt was shown to be highly susceptible to numerical error. The final method chosen was to use statically calibrated strain data and to filter out the system fundamental frequency such that the measured response yields close to the correct dynamic loading function. This method was shown to capture 93% of the pressure spectral energy using controlled load shaker testing. This method is one of the only practical ways for the inverse determination of the forcing function for non-stationary excitations, and, to the authors' knowledge, has not been described in the literature to date.

Brown, Andrew↗

NASA NTRS2020

SN 2018fif: The Explosion of a Large Red Supergiant Discovered in its Infancy by the Zwicky Transient Facility

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars' properties. In this paper, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of R = 744.0 (sup +183.0) (sub -128.0) R⨀ and an ejected mass of M(sub ej) = 9.3 (sup +0.4)(sub -5.8) M⨀. Our model also gives information on the explosion epoch, the progenitor's inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

Supernovae↗

NASA NTRS1991

Numerical performance estimates for a generic hypersonic forebody

The computation of a hypersonic flow past the forebody of a hypersonic vehicle is analyzed for two sets of flow conditions: the first represents flow conditions of an experimental test of the model in the Calspan 96-inch shock tunnel, and the second set is chosen to represent actual flight conditions. Solutions are derived for sharp- and blunt-nose versions of the geometry in order to understand the effects of the entropy layer on the forebody pressure and heat transfer rates. Some sensitivity of predicted heating rates to grid refinement is observed, but it is found to be small compared to the effects of bluntness. Real-gas effects are studied for the blunt-nose version at flight-like conditions, and these effects are found to have a significant effect on inlet-face performance measures such as mass capture and kinetic energy efficiency.

Lawrence, Scott L.↗

NASA NTRS1993

Grid adaption using Chimera composite overlapping meshes

The objective of this paper is to perform grid adaptation using composite over-lapping meshes in regions of large gradient to capture the salient features accurately during computation. The Chimera grid scheme, a multiple overset mesh technique, is used in combination with a Navier-Stokes solver. The numerical solution is first converged to a steady state based on an initial coarse mesh. Solution-adaptive enhancement is then performed by using a secondary fine grid system which oversets on top of the base grid in the high-gradient region, but without requiring the mesh boundaries to join in any special way. Communications through boundary interfaces between those separated grids are carried out using tri-linear interpolation. Applications to the Euler equations for shock reflections and to a shock wave/boundary layer interaction problem are tested. With the present method, the salient features are well resolved.

Kao, Kai-Hsiung↗

NASA NTRS1993

Grid adaptation using Chimera composite overlapping meshes

Kao, Kai-Hsiung↗

NASA NTRS1994

Carbonaceous Survivability on Impact

In order to gain knowledge about the potential contributions of comets and cosmic dust to the origin of life on Earth, we need to explore the survivability of their potential organic compounds on impact and the formation of secondary products that may have arisen from the chaotic events sustained by the carriers as they fell to Earth. We have performed a series of hypervelocity impact experiments using carbon-bearing impactors (diamond, graphite, kerogens, PAH crystals, and Murchison and Nogoya meteorites) into Al plate targets at velocities - 6 km/s. Estimated peak shock pressures probably did not exceed 120 GPa and peak shock temperatures were probably less than 4000 K for times of nano- to microsecs. Nominal crater dia. are less than one mm. The most significant results of these experiments are the preservation of the higher mass PAHs (e. g., pyrene relative to napthalene) and the formation of additional alkylated PAHs. We have also examined the residues of polystyrene projectiles impacted by a microparticle accelerator into targets at velocities up to 15 km/s. This talk will discuss the results of these experiments and their implications with respect to the survival of carbonaceous deliverables to early Earth. The prospects of survivability of organic molecules on "intact" capture of cosmic dust in space via soft: and hard cosmic dust collectors will also be discussed.

Bunch, T. E.↗

NASA NTRS2020

Simulating the Benesov Bolide Flowfield and Spectrum at Altitudes of 47 and 57 km

This paper develops a computational fluid dynamics capability for simulating the radiative emission from a meteor shock-layer and wake to an external observer. The developed capability includes the impact of radiation and ablation on the meteor flowfield, where coupled radiation accounts for the impact of radiation on the flowfield energy equations and coupled ablation accounts for the injection of ablation species from the meteoroid surface into the ow. This capability includes updated flowfield chemistry and nonequilibrium radiation models, which are compiled from the literature. To provide a level of validation for this capability, the Benesov spectral measurements are considered. Although the meteoroid was likely fragmented at the 47 and 57 km altitudes considered, and the measurement uncertainty is roughly ±50%, these measurements represent the best available spectral measurements for a relatively large (meter-class) bolide. To determine the equivalent meteoroid diameter to simulate, the continuum component of the spectrum in the 570 to 610 nm range is considered. This range is dominated by air emission from the high pressure and temperature shock-layer in front of the meteoroid, which is simulated with a relatively small uncertainty (because of the abundance of relevant shock tube measurements) and sensitive to the equivalent meteoroid diameter. This allows an equivalent diameter between 0.62 to 0.9 m to be determined, based on comparisons between measurements and simulations in the 570 to 610 nm range. The range of equivalent diameters is due to the various meteoroid geometries considered, where increasing the bluntness decreases the required diameter. Applying these diameters to simulations at 47 and 57 km results in simulated spectra that compare within 30% of the measured values, considering the wavelength integrated values between 400 and 650 nm. This agreement is well within the estimated measurement uncertainty of ±50%. All major spectral features are captured by the simulations. This agreement provides a level of validation for the developed model that has not been previously available.

Luminous efficiency↗

NASA NTRS1991

Flowfield analysis of modern helicopter rotors in hover by Navier-Stokes method

The viscous, three-dimensional, flowfields of UH60 and BERP rotors are calculated for lifting hover configurations using a Navier-Stokes computational fluid dynamics method with a view to understand the importance of planform effects on the airloads. In this method, the induced effects of the wake, including the interaction of tip vortices with successive blades, are captured as a part of the overall flowfield solution without prescribing any wake models. Numerical results in the form of surface pressures, hover performance parameters, surface skin friction and tip vortex patterns, and vortex wake trajectory are presented at two thrust conditions for UH60 and BERP rotors. Comparison of results for the UH60 model rotor show good agreement with experiments at moderate thrust conditions. Comparison of results with equivalent rectangular UH60 blade and BERP blade indicates that the BERP blade, with an unconventional planform, gives more thrust at the cost of more power and a reduced figure of merit. The high thrust conditions considered produce severe shock-induced flow separation for UH60 blade, while the BERP blade develops more thrust and minimal separation. The BERP blade produces a tighter tip vortex structure compared with the UH60 blade. These results and the discussion presented bring out the similarities and differences between the two rotors.

Srinivasan, G. R.↗

NASA NTRS2019

Acoustic Analysis of Sensor Ports

Acquiring data in rocket engines that is representative of the actual dynamic environment can often be difficult due to a multitude of influences. One source of contamination that is often not considered entirely is the response associated with the acoustic cavity created by a sensor offset. It is common to offset a sensor due to various reasons such as for mounting accessibility, thermal isolation, shock reduction, or prevention of debris impingement. While estimating the natural frequency of the acoustic cavity is straightforward, limited analysis has been described on the determination of the overall frequency response. The sensor port design approach usually attempts to ensure the port is short enough such that the acoustic response is negligible near the frequency-of-interest, but this requires knowledge of the frequency response and simple rules-of-thumb are not always guaranteed. Data correction and/or data interpretation is also often desired for an unsatisfactory response. The limited response analysis in the literature only offers approximations or neglects important contributions. A new approach is devised theoretically and computationally that captures the true acoustic response of a sensor port. This paper summarizes the acoustics background, the port response theoretical development, and provides comparisons of a port acoustic response using an analytical model and computational acoustics. The effects of nonlinear acoustics are also examined. Additionally, the paper summarizes the design of a specialized filter using the predicted sensor port response that can be applied to data for correction.

Casiano, M. J.↗

NASA NTRS1998

Particle Imaging Velocimetry Used in a Transonic Compressor Facility

Particle Imaging Velocimetry (PIV) is an optical technique whereby a pulsed laser sheet is used to illuminate particles entrained in a fluid across an extended planar cross section of a flow field. Electronic recording of the particle positions at two closely timed laser pulses permits the computation of the flow velocity. PIV captures the instantaneous flow field, permitting the study of unsteady flow phenomena. Mean flow statistics can be computed by acquiring several hundred images and averaging the results. The first-ever successful application of PIV to acquire measurements in a high-speed rotating turbomachinery blade row was completed in NASA Lewis Research Center's W-8 Single Stage Axial Compressor Facility. Measurements were acquired in a 20-in.-diameter transonic compressor rotor operating at 17,188 rpm. A custom-designed light-sheet generating probe was used to insert the high-energy, pulsed light-sheet illumination required for recording the unblurred images of particles entrained in the fluid. Measurements of the shock wave formed within the rotor blade passage and of unsteady structures within the blade wakes were acquired. These measurements provide insight into unsteady spatial structures in the flow field which cannot be measured with the more conventional laser anemometry technique. The PIV technique provides both instantaneous and average velocity data in a transonic compressor in an order of magnitude less time than required for other conventional optical diagnostic techniques.

Wernet, Mark P.↗

NASA NTRS2019

Field Evaluations of Sampling, Interviewing, and Flight Tracking of NASA's Low Boom Flight Demonstrator Aircraft

The first year’s effort identified sampling and interviewing as the principal risks to assessment of prompt reactions to overflights producing low-amplitude sonic booms. It also 1) established the utility of geo-information system-based route planning for LBFD flight missions, 2) developed and demonstrated a prototype of a geographically-distributed, Internet-enabled instrumentation system capable of wide-area tracking of LBFD aircraft in near-real time. The latter system permits synchronizing the conduct of interviews in multiple overflown communities with arrival times of shock waves at interviewing sites; and of measuring, archiving, and processing their acoustic signatures. Means were also recommended for constructing representative, telephone-based samples of eligible respondents living in households within carpet boom corridors adjacent to LBFD flight tracks, and for conducting interviews with cross-sectional (independent) samples of such respondents about their prompt reactions to exposure to low-amplitude sonic booms. A detailed study design was prepared and accepted by NASA for a set of single-contact attempt telephone interviews with a nationally representative sample of households. The study design focused on testing automated and live agent interview completion rates obtainable without callbacks. A minimal (two monitoring station) version of the aircraft tracking system was built and installed near a civil airport in a successful demonstration of the system’s ability to detect and track aircraft movements. The field exercise also demonstrated the ability of the system to capture the acoustic emissions of departing aircraft, and to serve aircraft position and sound level information to remote, geographically-distributed analysts in near-real time. Upon approval of OMB and IRB of the detailed study plan, a stratified, nationally representative sample of landline and wireless telephone-subscribing households was constructed. A total of 12,734 telephone interview contact attempts of the sort required by a straightforward cross-sectional study design were then made. These contact attempts demonstrated the impracticality of conducting a time-critical, cross-sectional study of prompt community response to low-amplitude sonic booms by means of “independent” (single contact attempt per respondent for each LBFD flight mission) telephone samples of respondents. The observed interview completion rates for these single telephone contact attempts were so low (~ 1% to 3% for automated and live agent interviews, respectively) that: 1) the representativeness of collected opinions would be susceptible to intuitive challenge as inadequate, even absent conclusive evidence of non-representativeness. Refuting challenges to representativeness would have to demonstrate that the composition of the actual sample did not differ from that of the target population, a task that is tantamount to proving a negative; 2) the information required to refute allegations of non-representativeness would require a questionnaire considerably lengthier than that required simply to determine the prevalence of boom-induced startle and annoyance. Such a questionnaire would have to inquire about potentially sensitive and intrusive matters, including respondents’ age, gender, education, employment, home ownership, income, ethnicity, family size, and other demographic factors; and 3) unreasonable numbers of attempts would be required to re-contact households with unsuccessful initial contact attempts, given the limited time available for doing so. For example, if about 500 completed interviews were desired in a supersonically overflown community, approximately 50,000 automated interview attempts would have to be made within ten to fifteen minutes of each LBFD overflight. Such large numbers of contact attempts could well exceed the numbers of households available for interview in areas of similar boom exposure levels in some communities near LBFD flight tracks. Such large numbers of interviews could be cost-effectively undertaken only by means of automated (i.e., outgoing interactive voice response) interviewing, a data collection method ill-suited for complex and sensitive questionnaire items. The infeasibility of independent sampling for evaluating prompt responses to LBFD overflights in a cross-sectional study is due in large part to simple non-response: that is, potential respondents – particularly those contacted on wireless telephones – refusing to answer calls with unfamiliar caller IDs. It is also due in part, however, to 1) the lack of time to attempt to contact the same respondent more than once within a few minutes after the arrival of a shock wave at the respondent’s location; and 2) the need to place calls during weekday/daytime hours, when response rates are notably lower than during evenings and weekends. Despite the poor interview completion rates achieved under the above constraints, cross sectional assessments of delayed reactions to LBFD overflights could still be feasible, if multiple attempts could be made to contact respondents during evening and weekend time periods, over extended time periods. Detailed plans for a longitudinal (panel) sample were developed as an alternative to a cross sectional sample design.

Fidell, Sanford↗

NASA NTRS2022

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Hydrodynamic instability↗

NASA NTRS2006

Very High-Speed Digital Video Capability for In-Flight Use

digital video camera system has been qualified for use in flight on the NASA supersonic F-15B Research Testbed aircraft. This system is capable of very-high-speed color digital imaging at flight speeds up to Mach 2. The components of this system have been ruggedized and shock-mounted in the aircraft to survive the severe pressure, temperature, and vibration of the flight environment. The system includes two synchronized camera subsystems installed in fuselage-mounted camera pods (see Figure 1). Each camera subsystem comprises a camera controller/recorder unit and a camera head. The two camera subsystems are synchronized by use of an MHub(TradeMark) synchronization unit. Each camera subsystem is capable of recording at a rate up to 10,000 pictures per second (pps). A state-of-the-art complementary metal oxide/semiconductor (CMOS) sensor in the camera head has a maximum resolution of 1,280 1,024 pixels at 1,000 pps. Exposure times of the electronic shutter of the camera range from 1/200,000 of a second to full open. The recorded images are captured in a dynamic random-access memory (DRAM) and can be downloaded directly to a personal computer or saved on a compact flash memory card. In addition to the high-rate recording of images, the system can display images in real time at 30 pps. Inter Range Instrumentation Group (IRIG) time code can be inserted into the individual camera controllers or into the M-Hub unit. The video data could also be used to obtain quantitative, three-dimensional trajectory information. The first use of this system was in support of the Space Shuttle Return to Flight effort. Data were needed to help in understanding how thermally insulating foam is shed from a space shuttle external fuel tank during launch. The cameras captured images of simulated external tank debris ejected from a fixture mounted under the centerline of the F-15B aircraft. Digital video was obtained at subsonic and supersonic flight conditions, including speeds up to Mach 2 and altitudes up to 50,000 ft (15.24 km). The digital video was used to determine the structural survivability of the debris in a real flight environment and quantify the aerodynamic trajectories of the debris.

Corda, Stephen↗

NASA NTRS1993

Laser Interferometer Skin-Friction measurements of crossing-shock wave/turbulent boundary-layer interactions

Wall shear stress measurements beneath crossingshock wave/turbulent boundary-layer interactions have been made for three interactions of different strengths. The interactions are generated by two sharp fins at symmetric angles of attack mounted on a flat plate. The shear stress measurements were made for fin angles of 7 and 11 degrees at Mach 3 and 15 degrees at Mach 4. The measurements were made using a Laser Interferometer Skin Friction (LISF) meter; a device which determines the wail shear by optically measuring the time rate of thinning of an oil film placed on the test model surface. Results of the measurements reveal high skin friction coefficients in the vicinity of the fin/plate junction and the presence of quasi-two-dimensional flow separation on the interaction centerline. Additionally, two Navier-Stokes computations, one using a Baldwin-Lomax turbulence model and one using a k-~ model, are compared to the experimental results for the Mach 4, 15 degree interaction case. While the k-~ model did a reasonable job of predicting the overall trend in portions of the skin friction distribution, neither computation fully captured the physics of the near surface flow in this complex interaction.

Garrison, T. J.↗

NASA NTRS1994

Two-Dimensional Scramjet Inlet Unstart Model: Wind-Tunnel Blockage and Actuation Systems Test

This supplement to NASA TM 109152 shows the Schlieren video (10 min. 52 sec., color, Beta and VHS) of the external flow field and a portion of the internal flow field of a two-dimensional scramjet inlet model in the NASA Langley 20-Inch Mach 6 Tunnel. The intent of the overall test program is to study (both experimentally and computationally) the dynamics of the inlet unstart; this (phase I) effort examines potential wind-tunnel blockage issues related to model sizing and the adequacy of the actuation systems in accomplishing the start and unstart. The model is equipped with both a moveable cowl and aft plug. Windows in the inlet sidewalls allow limited optical access to the internal shock structure. In the video, flow is from right to left, and the inlet is oriented inverted with respect to flight, i.e., with the cowl on top. The plug motion is obvious because the plug is visible in the aft window. The cowl motion, however, is not as obvious because the cowl is hidden from view by the inlet sidewall. The end of the cowl actuator arm, however, becomes visible above the inlet sidewalls between the windows when the cowl is up (see figure 1b of the primary document). The model is injected into the tunnel and observed though several actuation sequences with two plug configurations over a range of unit freestream Reynolds number at a nominal freestream Mach number of 6. The framing rate and shutter speed of the camera were too slow to fully capture the dynamics of the unstart but did prove sufficient to identify inlet start and unstart. This series of tests indicated that the model was appropriately sized for this facility and identified operability limits required first to allow the inlet to start and second to force the unstart.

Holland, Scott D.↗

NASA NTRS2001

Coupled-Flow Simulation of HP-LP Turbines Has Resulted in Significant Fuel Savings

Our objective was to create a high-fidelity Navier-Stokes computer simulation of the flow through the turbines of a modern high-bypass-ratio turbofan engine. The simulation would have to capture the aerodynamic interactions between closely coupled high- and low-pressure turbines. A computer simulation of the flow in the GE90 turbofan engine's high-pressure (HP) and low-pressure (LP) turbines was created at GE Aircraft Engines under contract with the NASA Glenn Research Center. The three-dimensional steady-state computer simulation was performed using Glenn's average-passage approach named APNASA. The areas upstream and downstream of each blade row mutually interact with each other during engine operation. The embedded blade row operating conditions are modeled since the average passage equations in APNASA actively include the effects of the adjacent blade rows. The turbine airfoils, platforms, and casing are actively cooled by compressor bleed air. Hot gas leaks around the tips of rotors through labyrinth seals. The flow exiting the high work HP turbines is partially transonic and, therefore, has a strong shock system in the transition region. The simulation was done using 121 processors of a Silicon Graphics Origin 2000 (NAS 02K) cluster at the NASA Ames Research Center, with a parallel efficiency of 87 percent in 15 hr. The typical average-passage analysis mesh size per blade row was 280 by 45 by 55, or approx.700,000 grid points. The total number of blade rows was 18 for a combined HP and LP turbine system including the struts in the transition duct and exit guide vane, which contain 12.6 million grid points. Design cycle turnaround time requirements ran typically from 24 to 48 hr of wall clock time. The number of iterations for convergence was 10,000 at 8.03x10(exp -5) sec/iteration/grid point (NAS O2K). Parallel processing by up to 40 processors is required to meet the design cycle time constraints. This is the first-ever flow simulation of an HP and LP turbine. In addition, it includes the struts in the transition duct and exit guide vanes.

Veres, Joseph P.↗