Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “highly deviated”

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

USM3D-ME Contributions to the 5th AIAA High Lift Prediction Workshop

This paper presents the results of Reynolds-averaged Navier-Stokes (RANS) simulations conducted by the NASA flow solver, mixed-element USM3D (USM3D-ME), for the 5th AIAA High-Lift Prediction Workshop. As part of the Fixed-Grid RANS Technology Focus Group (TFG), these simulations were performed to assess the accuracy and efficiency of the USM3D-ME solutions in predicting high-lift flows. The High-Lift Common Research Model (CRM-HL) served as the primary geometry. Several CRM-HL configurations were used for three case studies: a verification study (Case 1), a configuration buildup study (Case 2), and a Reynolds-number variation study (Case 3). Overall, USM3D-ME RANS results aligned with the solutions selected by the Fixed-Grid RANS TFG and available wind tunnel data. Simulations for Case 1 and Configuration 2.1 achieved machine-zero residual convergence, with aerodynamic coefficients converging to steady-state values. However, Configurations 2.2-2.4 and Case 3 encountered iterative- and grid-convergence challenges, particularly at high angles of attack. Compared with the experimental data available for Configurations 2.2-2.4, close agreement was demonstrated at low angles of attack. However, for angles of attack approaching the maximum lift conditions, the predicted lift coefficient and pitching moment deviated from experimental values. The drag-coefficient predictions were in relatively good agreement, however, slight overpredictions were observed at the highest angle of attack corresponding to the maximum-lift condition. Although iterative convergence for Configurations 2.2-2.4 at high angles of attack remains a persistent challenge, averaging aerodynamic coefficients over the last 5000 iterations yielded satisfactory agreement with the available wind tunnel experimental data. During the workshop, the lack of iterative convergence was attributed to the vortex structures emanating from the slat brackets. To investigate this issue further, post-workshop simulations were conducted on Configuration 2.2. In one study, RANS simulations were performed on a simplified geometry with the slat brackets removed. The second study focused on performing unsteady RANS (URANS) simulations on the original Configuration 2.2 geometry. Preliminary results from both studies are presented and compared with wind tunnel data for Configuration 2.2. Consistent with the findings of other participants in the Fixed-Grid RANS TFG, this study emphasizes the necessity for further exploration and advancement in RANS technology for predicting high-lift flows.

Aerodynamics↗
NASA NTRSDate not supplied

USM3D-ME Contributions to the 5th AIAA High Lift Prediction Workshop

This paper presents the results of Reynolds-averaged Navier-Stokes (RANS) simulations conducted by NASA’s flow solver, mixed-element USM3D (USM3D-ME), for the 5th AIAA High-Lift Prediction Workshop. As part of the Fixed-Grid RANS Technology Focus Group (TFG), these simulations were performed to assess the accuracy and efficiency of the USM3D-ME solutions in predicting high-lift flows. The High-Lift Common Research Model (CRM-HL) served as the primary geometry. Several CRM-HL configurations were used for three case studies: a verification study (Case 1), a configuration buildup study (Case 2), and a Reynolds-number variation study (Case 3). Overall, USM3D-ME RANS results aligned with the solutions selected by the Fixed-Grid RANS TFG and available wind tunnel data Simulations for Cases 1 and Configuration 2.1 achieved machine-zero residual convergence, with aerodynamic coefficients converging to steady-state values. However, Configurations 2.2-2.4 and Case 3 encountered iterative- and grid-convergence challenges, particularly at high angles of attack. Compared with the experimental data available for Configurations 2.2-2.4, close agreement was demonstrated at low angles of attack. However, for angles of attack approaching the maximum lift conditions, the predicted lift coefficient and pitching moment deviated from experimental values. The drag-coefficient predictions were in a relatively good agreement, however, slight overpredictions were observed at the highest angle of attack corresponding to the maximum-lift condition. Although iterative convergence for Configurations 2.2-2.4 at high angles of attack remains a persistent challenge, averaging aerodynamic coefficients over the last 5000 iterations yielded satisfactory agreement with the available wind tunnel experimental data. During the workshop, the lack of iterative convergence was attributed to the vortex structures emanating from the slat brackets. To investigate this issue further, post-workshop simulations were conducted on Configuration 2.2. In one study, RANS simulations were performed on a simplified geometry with the slat brackets removed. The second study focused on performing URANS simulations on the original Configuration 2.2 geometry. Preliminary results from both studies are presented and compared with wind tunnel data for Configuration 2.2. Consistent with the findings of other participants in the Fixed-Grid RANS TFG, this study emphasizes the necessity for further exploration and advancement in RANS technology for predicting high-lift flows.

CFD↗
NASA NTRS2000

High Order Discontinuous Gelerkin Methods for Convection Dominated Problems with Application to Aeroacoustics

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. On the analysis side, we have studied the efficient and stable discontinuous Galerkin framework for small second derivative terms, for example in Navier-Stokes equations, and also for related equations such as the Hamilton-Jacobi equations. This is a truly local discontinuous formulation where derivatives are considered as new variables. On the applied side, we have implemented and tested the efficiency of different approaches numerically. Related issues in high order ENO and WENO finite difference methods and spectral methods have also been investigated. Jointly with Hu, we have presented a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the RungeKutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method. Jointly with Hu, we have constructed third and fourth order WENO schemes on two dimensional unstructured meshes (triangles) in the finite volume formulation. The third order schemes are based on a combination of linear polynomials with nonlinear weights, and the fourth order schemes are based on combination of quadratic polynomials with nonlinear weights. We have addressed several difficult issues associated with high order WENO schemes on unstructured mesh, including the choice of linear and nonlinear weights, what to do with negative weights, etc. Numerical examples are shown to demonstrate the accuracies and robustness of the methods for shock calculations. Jointly with P. Montarnal, we have used a recently developed energy relaxation theory by Coquel and Perthame and high order weighted essentially non-oscillatory (WENO) schemes to simulate the Euler equations of real gas. The main idea is an energy decomposition under the form epsilon = epsilon(sub 1) + epsilon(sub 2), where epsilon(sub 1) is associated with a simpler pressure law (gamma)-law in this paper) and the nonlinear deviation epsilon(sub 2) is convected with the flow. A relaxation process is performed for each time step to ensure that the original pressure law is satisfied. The necessary characteristic decomposition for the high order WENO schemes is performed on the characteristic fields based on the epsilon(sub l) gamma-law. The algorithm only calls for the original pressure law once per grid point per time step, without the need to compute its derivatives or any Riemann solvers. Both one and two dimensional numerical examples are shown to illustrate the effectiveness of this approach.

Shu, Chi-Wang↗
NASA NTRS2021

System Identification for Propellers at High Incidence Angles

Propellers used for electric vertical takeoff and landing (eVTOL) aircraft propulsion systems experience a wide range of aerodynamic conditions, including large incidence angles relative to oncoming airflow. In oblique flow, propellers exhibit deviations in thrust and torque oriented along the propeller axis of rotation, as well as significant off-axis forces and moments. Although important for modeling eVTOL aircraft aerodynamics, sparse experimental data or mathematical models exist for propellers at incidence. This paper describes a propulsion system modeling methodology for the LA-8 tandem tilt-wing, eVTOL aircraft. System identification methods are applied to isolated propeller wind tunnel data gathered across the vehicle's flight envelope to develop a mathematical model of the propulsion system, including a static motor model, dynamic motor model, and propeller aerodynamic model. Modeling results validated against data withheld from the modeling process indicate good predictive capability and agree with theoretical expectations. The results are followed by a discussion of model implementation strategies into high-fidelity eVTOL aircraft simulations.

system identification↗
NASA NTRS2008

Compact, Highly Stable Ion Atomic Clock

A mercury-ion clock now at the breadboard stage of development (see figure) has a stability comparable to that of a hydrogen-maser clock: In tests, the clock exhibited an Allan deviation of between 2 x 10(exp -13) and 3 x 10(exp -13) at a measurement time of 1 second, averaging to about 10(exp -15) at 1 day. However, the clock occupies a volume of only about 2 liters . about a hundredth of the volume of a hydrogen-maser clock. The ion-handling parts of the apparatus are housed in a sealed vacuum tube, wherein only a getter pump is used to maintain the vacuum. Hence, this apparatus is a prototype of a generation of small, potentially portable high-precision clocks for diverse ground- and space-based navigation and radio science applications. Furthermore, this new ion-clock technology is about 100 times more stable and precise than the rubidium atomic clocks currently in use in the NAV STAR GPS Earth-orbiting satellites. In this clock, mercury ions are shuttled between a quadrupole and a 16-pole linear radio-frequency trap. In the quadrupole trap, the ions are tightly confined and optical state selection from a Hg-202 radio-frequency-discharge ultraviolet lamp is carried out. In the 16-pole trap, the ions are more loosely confined and atomic transitions resonant at frequency of about 40.507 GHz are interrogated by use of a microwave beam at that frequency. The trapping of ions effectively eliminates the frequency pulling caused by wall collisions inherent to gas-cell clocks. The shuttling of the ions between the two traps enables separation of the state-selection process from the clock microwave- resonance process, so that each of these processes can be optimized independently of the other. The basic ion-shuttling, two-trap scheme as described thus far is not new: it has been the basis of designs of prior larger clocks. The novelty of the present development lies in major redesigns of its physics package (the ion traps and the vacuum and optical subsystems) to effect the desired reduction of size to a volume of no more than a couple of liters. The redesign effort has included selection of materials for the vacuum tube, ion trap, and ultraviolet windows that withstand bakeout at a temperature of approx.450 C in preparation for sealing the tube to contain the vacuum. This part of the redesign effort follows the approach taken in the development of such other vacuum-tube electronic components as flight traveling- wave-tube amplifiers having operational and shelf lives as long as 15 years. The redesign effort has also included a thorough study of residual-gas-induced shifts of the ion-clock frequency and a study of alternative gases as candidates for use as a buffer gas within the sealed tube. It has been found that neon is more suitable than is helium, which has been traditionally used for this purpose, in that the pressure-induced frequency pulling by neon is between a third and a half of that of helium. In addition, because neon diffuses through solids much more slowly than does helium, the loss of neon by diffusion over the operational lifetime is expected to be negligible.

Prestage, John↗
NASA NTRSDate not supplied

System Identification for Propellers at High Incidence Angles

Propellers used for electric vertical takeoff and landing (eVTOL) aircraft propulsion systems experience a wide range of aerodynamic conditions, including large incidence angles relative to oncoming airflow. In oblique flow, propellers exhibit deviations in thrust and torque oriented along the propeller axis of rotation, as well as significant off-axis forces and moments. Although important for modeling eVTOL aircraft aerodynamics, sparse experimental data or mathematical models exist for propellers at incidence. This paper describes a propulsion system modeling methodology for the Langley Aerodrome No. 8 (LA-8) tandem tilt-wing, eVTOL aircraft. System identification methods are applied to isolated propeller wind tunnel data gathered across the vehicle's flight envelope to develop a mathematical model of the propulsion system, including a static motor model, dynamic motor model, and propeller aerodynamic model. Modeling results validated against data withheld from the modeling process indicate good predictive capability and agree with theoretical expectations. The results are followed by a discussion of model implementation strategies into high-fidelity eVTOL aircraft simulations.

system identification↗
NASA NTRS1981

Absolute S- and P-plane polarization efficiencies for high frequency holographic gratings in the VUV

High frequency plane gratings (3500 and 3600 gr/mm) have been holographically ruled and blazed for the VUV spectral region. All gratings were coated with 70 nm Al + 25 nm MgF2. Absolute unpolarized and S- and P-plane polarization efficiencies have been measured for the first and second orders in the 120- to 450-nm spectral region at 18.5 and 30 deg angles of deviation. For deep grooves, anomalous features are more pronounced for the P-plane polarization efficiency than for the S-plane polarization efficiency. Holographic gratings can be tailored to produce high polarization or low polarization in the VUV. For comparison, efficiencies and polarization of the best conventional high frequency gratings were also determined. Measurements show that scattered light is significantly lower for holographic gratings in the VUV when compared with the conventional gratings.

Caruso, A. J.↗
NASA NTRS1993

The FNAS Gamma Ray Observatory

Highly stressed and nonpotential magnetic fields in solar active region are believed to be the site and source of M or X class flares which release as much as 10(exp 32) ergs in a periods of 10 to 1000 s with intense emission detected over most of the electromagnetic spectrum. Such M or X class flares are the most likely gamma-ray emission candidates since sheared magnetic field is the means of energy storage for flare events. We have studied the correlation between 'nonpotentiality' parameters such as 'shear angle' and maximum magnetic field strength with the occurrences gamma-ray events. Also the angle of deviation of the magnetic field from the line of sight was studied since the bremsstrahlung gamma-ray beam produced by relativistic ions and electrons is highly directional in the forward direction and required non-vertical field-lines for the detection of the disk-center events.

Emslie, A. Gordon↗
NASA NTRS1983

The average solar wind in the inner heliosphere: Structures and slow variations

Measurements from the HELIOS solar probes indicated that apart from solar activity related disturbances there exist two states of the solar wind which might result from basic differences in the acceleration process: the fast solar wind (v 600 kms(-)1) emanating from magnetically open regions in the solar corona and the "slow" solar wind (v 400 kms(-)1) correlated with the more active regions and its mainly closed magnetic structures. In a comprehensive study using all HELIOS data taken between 1974 and 1982 the average behavior of the basic plasma parameters were analyzed as functions of the solar wind speed. The long term variations of the solar wind parameters along the solar cycle were also determined and numerical estimates given. These modulations appear to be distinct though only minor. In agreement with earlier studies it was concluded that the major modulations are in the number and size of high speed streams and in the number of interplanetary shock waves caused by coronal transients. The latter ones usually cause huge deviations from the averages of all parameters.

Schwenn, R.↗
NASA NTRS1984

Solar Terrestrial Influences on the D Region as Shown by the Level of Atmospheric Radio Noise

Measurements of the integrated atmospheric radio noise field strength at 27 kHz, used here, were made from 1965 to 1975 at Uppsala, Kuhlungsborn, and Prague-Panska Ves. The large scale meteorological situation was considered by comparing solar disturbed and undisturbed periods under similar weather situations. In order to show the effects of the precipitating high energy particle (HEP) flux and of the Forbush decrease on the noise level between pairs of stations were computed as deviations from the monthly median. Delta E (dB), day by day for all six periods was studied. The correlation coefficients for noon as well as for night values were computed. The correlation coefficients were compared with those for solar undisturbed periods.

Satori, G.↗
NASA NTRS1994

The evolution of synoptic ozone anomalies during the European Arctic Stratospheric Ozone Experiment in winter 1991/1992

The evolution of ozone anomalies over the middle and high latitudes of the Northern Hemisphere during the winter 1991-1992 is studied in this work. The largest monthly mean negative deviations in the middle latitudes of the Northern Hemisphere were about 10 percent in November and December, and up to 20 percent in January, February, and March over Eurasian territories, and much smaller over the Canadian sector. At the end of January, on individual days, total ozone values of 190-210 D.U. were observed over Eastern Europe and European part of Russia, that is 40-45 percent below normal. On the whole, the 1991-1992 winter was one of the most anomalous over all the period of ozone observations. Finally, an attempt is made to quantify the contribution of transport in the ozone layer changes over Europe during this period.

Zerefos, Christos S.↗
NASA NTRS2003

Comparison of Flow-Dependent and Static Error Correlation Models in the DAO Ozone Data Assimilation System

In a data assimilation system the forecast error covariance matrix governs the way in which the data information is spread throughout the model grid. Implementation of a correct method of assigning covariances is expected to have an impact on the analysis results. The simplest models assume that correlations are constant in time and isotropic or nearly isotropic. In such models the analysis depends on the dynamics only through assumed error standard deviations. In applications to atmospheric tracer data assimilation this may lead to inaccuracies, especially in regions with strong wind shears or high gradient of potential vorticity, as well as in areas where no data are available. In order to overcome this problem we have developed a flow-dependent covariance model that is based on short term evolution of error correlations. The presentation compares performance of a static and a flow-dependent model applied to a global three- dimensional ozone data assimilation system developed at NASA s Data Assimilation Office. We will present some results of validation against WMO balloon-borne sondes and the Polar Ozone and Aerosol Measurement (POAM) III instrument. Experiments show that allowing forecast error correlations to evolve with the flow results in positive impact on assimilated ozone within the regions where data were not assimilated, particularly at high latitudes in both hemispheres and in the troposphere. We will also discuss statistical characteristics of both models; in particular we will argue that including evolution of error correlations leads to stronger internal consistency of a data assimilation ,

Wargan, K.↗
NASA NTRS2005

A high performance frequency standard and distribution system for Cassini Ka-band experiment

This paper provides an overview and update of a specialized frequency reference system for the NASA Deep Space Network (DSN) to support Ka-band radio science experiments with the Cassini spacecraft, currently orbiting Saturn. Three major components, a Hydrogen Maser, Stabilized Fiber-optic Distribution Assembly (SFODA), and 10 Kelvin Cryocooled Sapphire Oscillator (10K CSO) and frequency-lock-loop, are integrated to achieve the very high performance, ground based frequency reference at a remote antenna site located 16 km from the hydrogen maser. Typical measured Allan Deviation is 1.6 -14 1 0a't 1 second and 1.7 x 10 -15 at 1000 seconds averaging intervals. Recently two 10K CSOs have been compared in situ while operating at the remote DSN site DSS-25. The CSO references were used operationally to downconvert the Ka band downlink received from the Cassini spacecraft in a series of occultation measurements performed over a 78 day period from March to June 2005.

sapphire↗
NASA NTRS2015

Scaling Relation for Occulter Manufacturing Errors

An external occulter is a spacecraft own along the line-of-sight of a space telescope to suppress starlight and enable high-contrast direct imaging of exoplanets. The shape of an external occulter must be specially designed to optimally suppress starlight and deviations from the ideal shape due to manufacturing errors can result loss of suppression in the shadow. Due to the long separation distances and large dimensions involved for a space occulter, laboratory testing is conducted with scaled versions of occulters etched on silicon wafers. Using numerical simulations for a flight Fresnel occulter design, we show how the suppression performance of an occulter mask scales with the available propagation distance for expected random manufacturing defects along the edge of the occulter petal. We derive an analytical model for predicting performance due to such manufacturing defects across the petal edges of an occulter mask and compare this with the numerical simulations. We discuss the scaling of an extended occulter test-bed.

Sirbu, Dan↗
NASA NTRS2020

Characterization of the Ryugu Surface By Means Of the Variability of the Near-Infrared Spectral Slope in NIRS3 Data

The Near-Earth Asteroid 162173 Ryugu (1999 JU3) was investigated by the JAXA Hayabusa2 mission from June 2018 to November 2019. The data acquired by NIRS3 spectrometer revealed a dark surface with a positive near-infrared spectral slope. In this work we investigated the spectral slope variations across the Ryugu surface, providing information about physical/chemical properties of the surface. We analysed the calibrated, thermally and photometrically corrected NIRS3 data and we estimated the mean value of spectral slope between 1.9 μm and 2.5 μm, corresponding to 0.163. Starting from the mean value of slope and moving in step of 1 standard deviation (0.022), we defined 9 “families of slope”, the Low-Red-Slope families (LR1, LR2 and LR3) and the High-Red-Sloped families (HR1, HR2, HR3, HR4, HR5, HR6). The mean values of some spectral parameters were estimated for each family, such as the reflectance factor at 1.9 μm, the spectral slope, the depth of bands at 2.7 μm and at 2.8 μm. A progressive spectral reddening, darkening and weakening/narrowing of OH bands is observed moving from the LR families to the HR families. We concluded that the spectral variability observed among families is the resulting contribution of the thermal metamorphism experienced by Ryugu after the catastrophic disruption of its parent body and space weathering processes that occurred on airless bodies as Ryugu, such as impact cratering and solar wind irradiation. As a consequence, the HR1, LR1, LR2 and LR3 families, corresponding to equatorial ridge and crater rims, are the less altered regions on Ryugu surface, which experienced the minor alteration and OH devolatilization; the HR2, HR3, HR4, HR5 families, coincident with floors and walls of impact craters, are the most altered areas, result of the three processes occurred on Ryugu. The strong reddening of the HR6 family (coincident with Ejima Saxum) is likely due to the fine-sized material covering the large boulder.

Ryugu↗
NASA NTRS1993

A study of the influence of satellite data on GLA analyses over the Atlantic Ocean during a period of blocking anticyclone development

This paper presents a statistical diagnosis of the differences between two Goddard Laboratory for Atmospheres analysis representations of temperature and wind speed, one containing satellite data influences (SAT analyses) and one without (NOSAT analyses), for a case of oceanic blocking anticyclone development. Results, obtained using area means and standard deviations of the two fields and correlation coefficients and root-mean-square differences between the two fields, indicate that the inclusion of satellite data can have significant impact on ocean-domain analyses. This is especially evident if one examines higher-order fields that represent gradients in the basic variable fields or contain covariances between two different variables (e.g., advections). More specifically, the inclusion of satellite data resulted in a cold (warm) temperature bias at low (high) levels and weaker temperature gradients, stronger winds, weaker vertical wind shears, and a warm-air advection bias at most levels. The SAT analyses also exhibited larger standard deviations than the NOSAT for wind speed, relative vorticity, vorticity advection, temperature advection, and 500-mb height tendencies, suggesting that for this case the spatial variability of the circulation features were enhanced by the inclusion of satellite data. Even so, since height tendencies forced by vorticity and temperature advection were less sensitive to the addition of satellite data than were the advection quantities themselves, the dynamics of the blocking system was apparently less influenced by the presence of satellite data than was the structure of the system.

Lamberty, Gregory L.↗
NASA NTRS1995

Rotational Energy Transfer Cross Sections in N2-N2 Collisions

Rotational inelastic transitions of N2 have been studied in the coupled state (CS, also called centrifugal sudden) and infinite-order-sudden (IOS) approximations, using the N2-N2 rigid-rotor potential of van der Avoird et al. For benchmarking purposes, close coupling (CC) calculations have also been carried out over a limited energy range and for even j - even j collisions only. Both the CC and CS cross sections have been obtained with and without exchange symmetry, whereas exchange is neglected in the IOS calculations. The CS results track the CC cross sections rather well. At total energies between 113 to 219 cm(exp -1) the average deviation is 14%. The deviation decrease with increasing energy, indicating that the CS approximation can be used as a substitute at higher energies when the CC calculations become impractical. Comparison between the CS and IOS cross sections at the high energy end of the CS calculation, 500 - 680 cm(exp-1), shows significant differences between the two. In addition, the IOS results exhibits sensitivity to the amount of inelasticity and the results for large DELTA J transitions are subjected to bigger errors. At total energy 113 cm(exp -1) and above, the average deviation between state-to-state cross sections calculated with even and odd exchange symmetries is 1.5%.

Huo, Winifred M.↗
NASA NTRS2003

Mixing and NO(x) Emission Calculations of Confined Reacting Jet Flows in a Cylindrical Duct

Rapid mixing of cold lateral jets with hot cross-stream flows in confined configurations is of practical interest in gas turbine combustors as it strongly affects combustor exit temperature quality, and gaseous emissions in for example rich-lean combustion. It is therefore important to further improve our fundamental understanding of the important processes of dilution jet mixing especially when the injected jet mass flow rate exceeds that of the cross-stream. The results reported in this report describe some of the main flow characteristics which develop in the mixing process in a cylindrical duct. A 3-dimensional tool has been used to predict the mixing flow field characteristics and NOx emission in a quench section of an RQL combustor, Eighteen configurations have been analyzed in a circular geometry in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying three parameters: 1) jet-to-mainstream momentum-flux ratio (J), 2) orifice shape or orifice aspect ratio, and 3) slot slant angle. The results indicate that the mixing flow field significantly varies with the value of the jet penetration and subsequently, slanting elongated slots generally improve the mixing uniformity at high J conditions. Round orifices produce more uniform mixing and low NO(x) emissions at low J due to the strong and adequate jet penetration. No significant correlation was found between the NO(x) production rates and the mixing deviation parameters, however, strong correlation was found between NO(x) formation and jet penetration. In the computational results, most of the NO(x) formation occurred behind the orifice starting at the orifice wake region. Additional NO(x) is formed upstream of the orifice in certain configurations with high J conditions due to the upstream recirculation.

Holdeman, James D.↗