Search Results

Testing of ducted fuel injection (DFI) in a single-cylinder engine with production-like hardware previously showed that adding a duct structure increased soot emissions at the full load, rated speed operating point [1]. The authors hypothesized that the DFI flame, which travels faster than a conventional diesel combustion (CDC) flame, and has a shorter distance to travel, was being re-entrained into the on-going fuel injection around the lift-off length (LOL), thus reducing air entrainment into the on-going injection. The engine operating condition and the engine combustion chamber geometry were duplicated in a constant pressure vessel. The experimental setup used a 3D piston section combined with a glass fire deck allowing for a comparison between a CDC flame and a DFI flame via high-speed imaging. CH* imaging of the 3D piston profile view clearly confirmed the re-entrainment hypothesis presented in the previous engine work.

A research program was conducted to evaluate the effectiveness of icing tunnel hybrid model design. A hybrid design is where the full-scale leading edge of a wing section is maintained only to a certain percentage of the local chord, while the aft section of the model is redesigned into a shortened or truncated planform. An initial study was conducted in 2020 where the ice shape geometries on a full-chord length version of the swept CRM65 wing model were compared to those from the hybrid version of CRM65 that were obtained in the NASA Icing Research Tunnel in 2015. The results were reported in a 2021 paper. For most test conditions, the overall size and shape of the ice shapes compared well. However, the ice shapes from the full-chord model were generally slightly smaller than those from the hybrid model.

Predicting the aerodynamic performance of an aircraft in icing conditions is critical as failures in an aircraft’s ice protection system can compromise flight safety. Aerodynamic effects of icing have typically relied on RANS modeling, which usually struggles to predict stall behavior, including those induced by surface roughness. Encouraged by recent studies using LES that demonstrate the ability to predict stall characteristics on full aircraft with smooth wings at an affordable cost [1], this study seeks to apply this methodology to icing conditions. Measurements of lift, drag, and pitching moments of a NACA23012 airfoil under clean and iced conditions are collected at Re = 1.8M. Using laser scanned, detailed representations of the icing geometries, LES calculations are conducted to compare integrated loads against experimental measurements in both clean and iced conditions at various angles of attack through the onset of stall [2].

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

The Icing Research Tunnel at NASA Glenn follows the recommended practice for calibration outlined in SAE’s ARP5905. The calibration team has followed the schedule of a full calibration every five years with a check calibration done every six months following. The liquid water content of the IRT has maintained stability within the stated specifications of variation within +/- 10% of the curve fit equation generated from calibration data. Using past measurements and data trends, IRT characterization engineers wanted to develop methods for the ability to know when data were not within variation. Trends can be observed in the liquid water content measurement process by constructing statistical process control charts. This paper describes data processing procedures for the Multi-Element Sensor in the IRT, including collision efficiency corrections, canonical correlation analysis, process for rejection of data, and construction of control charts.

When examining the literature on the adhesion strength of impact ice, there have been a wide range of methodologies tried to measure the required stresses to induce interfacial delamination. Utilizing the Icing Research Tunnel at the NASA Glenn Research Center to generate the impact ice required for this work, several different mechanical tests have been and are being developed to determine the stresses along the interface between ice and coupon. This set of tests includes the technical mature modified lap joint test which has been used to conduct ice adhesion studies through a wide sweep of icing conditions. To conduct in situ ice adhesion measurements inside of the Icing Research Tunnel, several new experiments are currently being developed to make ice adhesion measurements during and immediately after ice accretion.

A process of using machine learning to segment impact ice microstructure is presented and analyzed. The microstructure of impact ice has been shown to correlate with the adhesion strength of ice. Machine vision techniques are explored as a method of decreasing analysis time. The segmentation was conducted with the goal of obtaining average grain size estimations. The model was trained on a set of micrographs of impact ice grown at NASA Glenn’s Icing Research Tunnel. The model leveraged a model pre-trained on a large set of micrographs of various materials as a starting point. Post-processing of the segmented images was done to connect broken boundaries. An automatic method of determining grain size following an ASTM standard was implemented. Segmentation results using different training sets as well as different encoder and decoder pairs are presented. Calculated sizes are compared to manual grain size measurement methods.

This paper details results of a numerical and experimental investigation into the relative performance of vaned and vaneless mixed flow turbines for application to medium and heavy-duty diesel engines utilizing pulse exhaust systems. Previous investigations into the impact of nozzle vanes on turbine performance considered only open turbine housings, whereas a majority of medium and heavy-duty diesel engine applications are six-cylinder engines using pulse exhaust systems with divided turbines. The two turbine stages for this investigation were carefully designed to meet the constraints of engines with pulse exhaust systems and to control confounding factors that would undermine the vaned vs vaneless performance comparison. Detailed CFD analysis and turbine dynamometer test results confirm a significant efficiency advantage for the vaned turbine stage under both full and partial admission conditions.

Numerical solutions have been generated which simulate flow inside an aircraft engine flying at altitude through an ice crystal cloud. The geometry used for this study is the Honeywell Uncertified Research Engine (HURE) which was recently tested in the NASA Propulsion Systems Laboratory (PSL) in January 2018. The simulations were carried out at predicted operating points with a potential risk of ice accretion. The extent of the simulation is from upstream of the engine inlet to downstream past the strut in the core and bypass. The flow solution is produced using GlennHT, a NASA in-house code. A mixing plane approximation is used upstream and downstream of the fan. The use of the mixing plane allows for steady state solutions in the relative frame. The flow solution is then passed on to LEWICE3D for particle trajectory, impact and breakup prediction. The LEWICE3D code also uses a mixing plane approximation at the boundaries upstream and downstream of the fan.

NASA and the FAA conducted two flight campaigns to quantify onboard weather radar measurements with in-situ measurements of high concentrations of ice crystals found in deep convective storms. The ultimate goal of this research was to improve the understanding of high ice water content (HIWC) and develop onboard weather radar processing techniques to detect regions of HIWC ahead of an aircraft to enable tactical avoidance of the potentially hazardous conditions. Both HIWC RADAR campaigns utilized the NASA DC-8 Airborne Science Laboratory equipped with a Honeywell RDR-4000 weather radar and in-situ microphysical instruments to characterize the ice crystal clouds. The purpose of this paper is to summarize how these campaigns were conducted and highlight key results. The first campaign was conducted in August 2015 with a base of operations in Ft. Lauderdale, Florida.

Artificial ice shapes of various geometric fidelity were tested on a wing model based on the Common Research Model. Low Reynolds number tests were conducted at Wichita State University’s Walter H. Beech Memorial Wind Tunnel utilizing an 8.9% scale model, and high Reynolds number tests were conducted at ONERA’s F1 wind tunnel utilizing a 13.3% scale model. Several identical geometrically-scaled ice shapes were tested at both facilities, and the results were compared at overlapping Reynolds and Mach numbers. This was to ensure that the results and trends observed at low Reynolds number could be applied and continued to high, near-flight Reynolds number. The data from Wichita State University and ONERA F1 agreed well at matched Reynolds and Mach numbers. The lift and pitching moment curves agreed very well for most configurations.

Understanding the aerodynamic impact of swept-wing ice accretions is a crucial component of the design of modern aircraft. Computer-simulation tools are commonly used to approximate ice shapes, so the necessary level of detail or fidelity of those simulated ice shapes must be understood relative to high-fidelity representations of the ice. Previous tests were performed in the NASA Icing Research Tunnel to acquire high-fidelity ice shapes. From this database, full-span artificial ice shapes were designed and manufactured for both an 8.9%-scale and 13.3%-scale semispan wing model of the CRM65 which has been established as the full-scale baseline for this swept-wing project. These models were tested in the Walter H. Beech wind tunnel at Wichita State University and at the ONERA F1 facility, respectively. The data collected in the Wichita St.

The influence of freestream static temperature on roughness temporal evolution and spatial variation was investigated in the Icing Research Tunnel (IRT) at NASA Glenn Research Center. A 53.34 cm (21-in.) NACA 0012 airfoil model and a 152.4 cm (60-in.) HAARP-II business jet airfoil model were exposed to Appendix C clouds for fixed exposure times and thus fixed ice accumulation parameter. For the base conditions, the static temperature was varied to produce different stagnation point freezing fractions. The resulting ice shapes were then scanned using a ROMER Absolute Arm system and analyzed using the self-organizing map approach of McClain and Kreeger. The ice accretion prediction program LEWICE was further used to aid in interrogations of the ice accretion point clouds by using the predicted surface variations of local collection efficiency.

Topology of liner finish is critical to the performance of internal combustion engines. Proper liner finish simulation processes lead to efficient engine design and research. Fourier methods have been well studied to numerically generate liner topology. However, three major issues wait to be addressed to make the generation processes feasible and reliable. First, in order to simulate real plateau honed liners, approaches should be developed to calculate accurate liner geometric parameters. These parameters are served as the input of the generation algorithm. Material ratio curve, the common geometry calculation method, should be modified so that accurate root mean square of plateau height distribution could be obtained. Second, the set of geometric parameters used in generating liner finish (ISO 13565-2) is different from the set of parameters used in manufacturing industry (ISO 13565-3). Quantitative relations between these two sets should be studied.

The Society of Automotive Engineers Fatigue Design & Evaluation Committee [SAE FD&E] is actively working on a total life project for weldments, in which the welding residual stress is a key contributor to an accurate assessment of fatigue life. Physics-based welding process simulation and various types of residual stress measurements were pursued to provide a representation of the residual stress field at the failure location in the fatigue samples. A well-controlled and documented robotic welding process was used for all sample fabrications to provide accurate inputs for the welding simulations. One destructive (contour method) residual stress measurement and several non-destructive residual stress measurements-surface X-ray diffraction (XRD), energy dispersive X-ray diffraction (EDXRD), and neutron diffraction (ND)-were performed on the same or similarly welded samples.

This paper introduces a data driven model for predicting airport arrival capacity with 2-8 hour look-ahead forecast data. The model is suitable for air traffic flow management by explicitly investigating the impact of convective weather on airport arrival meter fix throughput. Estimation of the arrival airport capacity under arrival meter fix flow constraints due to severe weather is an important part of Air Traffic Management (ATM). Airport arrival capacity can be reduced if one or more airport arrival meter fixes are partially or completely blocked by convective weather. When the predicted airport arrival demands exceed the predicted available airport’s arrival capacity for a sustained period, Ground Delay Program (GDP) operations will be triggered by ATM system.

A new multi-scale curved beam based model was developed for piston-ring designs. This tool is able to characterize the behavior of a ring with any cross section design. This paper describes the conformability and ring static twist calculation. The conformability part model the static behavior of the ring in working engine conditions. The model employs the computation scheme that separates the meshing of the structure and local force generation. Additional to the conventional static ring-bore conformability analysis, the conformability model is designed to examine ring-bore and ring-groove interactions in a running engine under varying driving forces and localized lubrication conditions. We made Improvements on the way to handle the effects of the radial temperature gradient compared to the existing models. Examples are given on the effects of ring rotation on the interaction of the ring and a distorted bore as well as the change of local lubrication conditions.

A new multi-scale curved beam based model was developed for piston-ring designs. This paper describes the free-shape, force in circular bore and closed shape within a flexible band (ovality) related parts. Knowing any one of these distributions, this model determines the other two. This tool is useful in the sense that the characterization of the ring is carried out by measuring its closed shape within a flexible band which is more accurate than measuring its free shape or force distribution in circular bore. Thus, having a model that takes the closed shape within a flexible band as an input is more convenient and useful based on the experiments carried out to characterize the ring.

The blade crossing event of a coaxial counter-rotating rotor is a potential source of noise and impulsive blade loads. Blade crossings occur many times during each rotor revolution. In previous research by the authors, this phenomenon was analyzed by simulating two airfoils passing each other at specified speeds and vertical separation distances, using the compressible Navier-Stokes solver OVERFLOW. The simulations explored mutual aerodynamic interactions associated with thickness, circulation, and compressibility effects. Results revealed the complex nature of the aerodynamic impulses generated by upper/lower airfoil interactions. In this paper, the coaxial rotor system is simulated using two trains of airfoils, vertically offset, and traveling in opposite directions. The simulation represents multiple blade crossings in a rotor revolution by specifying horizontal distances between each airfoil in the train based on the circumferential distance between blade tips.

A new ring pack model has been developed based on the curved beam finite element method. This paper describes the first part of this model: simulating gas pressure in different regions above piston skirt and ring dynamic behavior of two compression rings and a twin-land oil control ring. The model allows separate grid divisions to resolve ring structure dynamics, local force/pressure generation, and gas pressure distribution. Doing so enables the model to capture both global and local processes at their proper length scales. The effects of bore distortion, piston secondary motion, and groove distortion are considered. Gas flows, gas pressure distribution in the ring pack, and ring structural dynamics are coupled with ring-groove and ring-liner interactions, and an implicit scheme is employed to ensure numerical stability. The model is applied to a passenger car engine to demonstrate its ability to predict global and local effects on ring dynamics and oil transport.