Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RESPONSE FUNCTIONS”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

OLTARIS: On-Line Tool for the Assessment of Radiation in Space

The On-Line Tool for the Assessment of Radiation In Space (OLTARIS) is a World Wide Web based tool that assesses the effects of space radiation to humans in items such as spacecraft, habitats, rovers, and spacesuits. This document explains the basis behind the interface and framework used to input the data, perform the assessment, and output the results to the user as well as the physics, engineering, and computer science used to develop OLTARIS. The physics is based on the HZETRN2005 and NUCFRG2 research codes. The OLTARIS website is the successor to the SIREST website from the early 2000 s. Modifications have been made to the code to enable easy maintenance, additions, and configuration management along with a more modern web interface. Over all, the code has been verified, tested, and modified to enable faster and more accurate assessments. The next major areas of modification are more accurate transport algorithms, better uncertainty estimates, and electronic response functions. Improvements in the existing algorithms and data occur continuously and are logged in the change log section of the website.

Singleterry, Robert C., Jr.↗
NASA NTRS2009

On Calibrations Using the Crab Nebula as a Standard Candle

Inspired by a recent paper (Kirsch et al. 2005) on possible use of the Crab Nebula as a standard candle for calibrating X-ray response func tions, we examine possible consequences of intrinsic departures from a single (absorbed) power law upon such calibrations. We limited our analyses to three more modern X-ray instruments -- the ROSAT/PSPC, th e RXTE/PCA, and the XMM-Newton/EPIC-pn. The results are unexpected an d indicate a need to refine two of the three response functions studi ed. The implications for Chandra will be discussed.

Weisskopf, Martin↗
NASA NTRS2010

CHOICE - Directed Study: Consequences of longterm- Confinement and Hypobaric HypOxia on Immunity in the Antarctic Concordia Environment

Concerning ground-based space physiological research, the choice of analog must carefully match the system of interest. For spaceflight-associated immune dysregulation (SAID), Antarctica winter-over has emerged as potentially the best terrestrial analog. The prolonged mission durations, extreme/dangerous environment, station-based lifestyle, isolation from outside world, disrupted circadian rhythms, and other psychological aspects make this analog extremely high fidelity for exploration-class space missions (long duration lunar, Mars). NASA, ESA and RSA are currently investigating SAID, with NASA currently operating the Integrated Immune flight study. It is desirable to have a ground analog for SAID validated, so that potential countermeasures might be validated terrestrially prior to during flight. For this presentation, NASA data collected on the winterover 2009 crewmembers, baseline through early deployment will be presented. Through early deployment (approximately 2-3 weeks at Concordia), phenotypic alterations included increased levels of memory T cells, shifts among the CD8+ T cell compartment to a more mature phenotype, and increases in constitutively activated T cells. CD8+/IFNg+ T cell percentages, and T cell blastogenesis functional responses were depressed early deployment as compared to healthy controls. In four compatible subjects, secreted T cell Th1/Th2 cytokines were measured following culture stimulation, and a Th2 shift was observed as compared to controls. Post-winter over frozen sample return will be required to determine if this shift persisted during the winter over period. Additionally, circadian rhythms remained altered compared to baseline, as determined through 5x daily cortisol measurements. Latent viral reactivation will not be determined until frozen sample return occurs.

Sams, Clarence↗
NASA NTRS2010

NASA's Advancements in Space-Based Spectrometry Lead to Improvements in Weather Prediction and Understanding of Climate Processes

AIRS (Atmospheric Infra-Red Sounder), was launched, in conjunction with AMSU-A (Advanced Microwave Sounding Unit-A) on the NASA polar orbiting research satellite EOS (Earth Observing System) Aqua satellite in May 2002 as a next generation atmospheric sounding system. Atmospheric sounders provide information primarily about the vertical distribution of atmospheric temperature and water vapor distribution. This is achieved by measuring outgoing radiation in discrete channels (spectral intervals) which are sensitive primarily to variations of these geophysical parameters. The primary objectives of AIRS/AMSU were to utilize such information in order to improve the skill of numerical weather prediction as well as to measure climate variability and trends. AIRS is a multi-detector array grating spectrometer with 2378 channels covering the spectral range 650/cm (15 microns) to 2660/cm (3.6 microns) with a resolving power (i/a i) of roughly 1200 where a i is the spectral channel bandpass. Atmospheric temperature profile can be determined from channel observations taken within the 15 micron (the long-wave CO2 absorption band) and within the 4.2 micron (the short-wave CO2 absorption band). Radiances in these (and all other) spectral intervals in the infrared are also sensitive to the presence of clouds in the instrument?s field of view (FOV), which are present about 95% of the time. AIRS was designed so as to allow for the ability to produce accurate Quality Controlled atmospheric soundings under most cloud conditions. This was achieved by having 1) extremely low channel noise values in the shortwave portion of the spectrum and 2) a very flat spatial response function within a channel?s FOV. IASI, the high spectral resolution IR interferometer flying on the European METOP satellite, does not contain either of these important characteristics. The AIRS instrument was also designed to be extremely stabile with regard to its spectral radiometric characteristics, which is critical with regard to the ability to measure accurate long term trends.

Susskind, Joel↗
NASA NTRS2009

Elemental Composition of 433 Eros: New Calibration of the NEAR-Shoemaker XRS Data

We present a new calibration of the elemental-abundance data for Asteroid 433 Fros taken by the X-ray spectrometer (XRS) aboard the NEAR-Shoemaker spacecraft. (NEAR is an acronym for "Near-Earth Asteroid Rendezvous,") Quintification of the asteroid surface elemental abundance ratios depends critically on accurate knowledge of the incident solar X-ray spectrum, which was monitored simultaneously with asteroid observations. Previously published results suffered from incompletely characterized systematic uncertainties due to an imperfect ground calibrations of the NEAR gas solar monitor. The solar monitor response function and associated uncertainties have now been characterized by cross-calibration of a large sample of NEAR solar monitor flight data against. contemporary broadband solar X-ray data from the Earth-orbiting GOES-8 (Geostationary Operational Environmental Satellite). The results have been used to analyze XRS spectra acquired from Eros during eight major solar flares (including three that have not previously been reported). The end product of this analysis is a revised set of Eros surface elemental abundance ratios with new error estimates that more accurately reflect the remaining uncertainties in the solar flare spectra: Mg/Si=.753 +0.078/-0.055, Al/Si=0.069 +/-0.055, S/Si=0.005+/-0.008, Ca/Si=0.060+0.023/-0.024, and Fe/Si= 1.578+0.338/-0.320. These revised abundance ratios are consitent within cited uncertainties with the results of Nittler et al. [Nittler, L.R., and 14 colleagues, 2001. Meteorit Planet. Sci 36, 1673-1695] and thus support the prior conclusions that 433 Eros has major-element composition simular to ordinary chondrites with the exception of a stong depletoin in sulfur, most likely caused by space weathering.

Lim, Lucy F.↗
NASA NTRS2010

Dst Index in the 2008 GEM Modeling Challenge - Model Performance for Moderate and Strong Magnetic Storms

The GEM 2008 modeling challenge efforts are expanding beyond comparing in-situ measurements in the magnetosphere and ionosphere to include the computation of indices to be compared. The Dst index measures the largest deviations of the horizontal magnetic field at 4 equatorial magnetometers from the quiet-time background field and is commonly used to track the strength of the magnetic disturbance of the magnetosphere during storms. Models can calculate a proxy Dst index in various ways, including using the Dessler-Parker Sckopke relation and the energy of the ring current and Biot-Savart integration of electric currents in the magnetosphere. The GEM modeling challenge investigates 4 space weather events and we compare models available at CCMC against each other and the observed values of Ost. Models used include SWMF/BATSRUS, OpenGGCM, LFM, GUMICS (3D magnetosphere MHD models), Fok-RC, CRCM, RAM-SCB (kinetic drift models of the ring current), WINDMI (magnetosphere-ionosphere electric circuit model), and predictions based on an impulse response function (IRF) model and analytic coupling functions with inputs of solar wind data. In addition to the analysis of model-observation comparisons we look at the way Dst is computed in global magnetosphere models. The default value of Dst computed by the SWMF model is for Bz the Earth's center. In addition to this, we present results obtained at different locations on the Earth's surface. We choose equatorial locations at local noon, dusk (18:00 hours), midnight and dawn (6:00 hours). The different virtual observatory locations reveal the variation around the earth-centered Dst value resulting from the distribution of electric currents in the magnetosphere during different phases of a storm.

Rastaetter, Lutz↗
NASA NTRS2011

Development of the Orion Crew Module Static Aerodynamic Database: Hypersonic - Part 1

The Orion aerodynamic database provides force and moment coefficients given the velocity, attitude, configuration, etc. of the Crew Exploration Vehicle (CEV). The database is developed and maintained by the NASA CEV Aerosciences Project team from computational and experimental aerodynamic simulations. The database is used primarily by the Guidance, Navigation, and Control (GNC) team to design vehicle trajectories and assess flight performance. The initial hypersonic re-entry portion of the Crew Module (CM) database was developed in 2006. Updates incorporating additional data and improvements to the database formulation and uncertainty methodologies have been made since then. This paper details the process used to develop the CM database, including nominal values and uncertainties, for Mach numbers greater than 8 and angles of attack between 140deg and 180deg. The primary available data are more than 1000 viscous, reacting gas chemistry computational simulations using both the Laura and Dplr codes, over a range of Mach numbers from 2 to 37 and a range of angles of attack from 147deg to 172deg. Uncertainties were based on grid convergence, laminar-turbulent solution variations, combined altitude and code-to-code variations, and expected heatshield asymmetry. A radial basis function response surface tool, NEAR-RS, was used to fit the coefficient data smoothly in a velocity-angle-of-attack space. The resulting database is presented and includes some data comparisons and a discussion of the predicted variation of trim angle of attack and lift-to-drag ratio. The database provides a variation in trim angle of attack on the order of +/-2deg, and a range in lift-to-drag ratio of +/-0.035 for typical vehicle flight conditions.

Bibb, Karen L.↗
NASA NTRS2011

Tracking Climate Effects on Plant-Pollinator Interaction Phenology with Satellites and Honey Bee Hives

Background/Question/Methods: The complexity of plant-pollinator interactions, the large number of species involved, and the lack of species response functions present challenges to understanding how these critical interactions may be impacted by climate and land cover change on large scales. Given the importance of this interaction for terrestrial ecosystems, it is desirable to develop new approaches. We monitor the daily weight change of honey bee (Apis mellifera) colonies to record the phenology of the Honey Bee Nectar Flow (HBNF) in a volunteer network (honeybeenet.gsfc.nasa.gov). The records document the successful interaction of a generalist pollinator with a variety of plant resources. We extract useful HBNF phenology metrics for three seasons. Sites currently exist in 35 states/provinces in North America, with a concentration in the Mid-Atlantic region. HBNF metrics are compared to standard phenology metrics derived from remotely sensed vegetation indices from NASA's MODIS sensor and published results from NOAA's A VHRR. At any given time the percentage of plants producing nectar is usually a sma11 fraction of the total satellite sensor signal. We are interested in determining how well the 'bulk' satellite vegetation parameters relate to the phenology of the HBNF, and how it varies spatially on landscape to continental scales. Results/Conclusions: We found the median and peak seasonal HBNF dates to be robust, with variation between replicate scale hives of only a few days. We developed quality assessment protocols to identify abnormal colony artifacts. Temporally, the peak and median of the HBNF in the Mid-Atlantic show a significant advance of 0.58 d/y beginning about 1970, very similar to that observed by the A VHRR since 1982 (0.57 d/y). Spatially, the HBNF metrics are highly correlated with elevation and winter minimum temperature distribution, and exhibit significant but regionally coherent inter-annual variation. The relationship between median of the spring HBNF with the "Green-up" metric from the 500 meter MODIS NDVI phenology product, for sites throughout the Eastern US 2000-2009, is well described by a single linear fit (r(exp 2) = 0.72). We conclude.that for the tree-dominated areas of the Eastern US at least the spring HBNF can be tracked very well by MODIS phenology. Analysis of other regions and seasons is presently underway but with more limited data. Spatial patterns in the eastern US and management implications will be presented and discussed.

Esaias, Wayne E.↗
NASA NTRS2005

Simulation of Range Safety for the NASA Space Shuttle

This paper describes a simulation environment that seamlessly combines a number of safety and environmental models for the launch phase of a NASA Space Shuttle mission. The components of this simulation environment represent the different systems that must interact in order to determine the Expectation of casualties (E(sub c)) resulting from the toxic effects of the gas dispersion that occurs after a disaster affecting a Space Shuttle within 120 seconds of lift-off. The utilization of the Space Shuttle reliability models, trajectory models, weather dissemination systems, population models, amount and type of toxicants, gas dispersion models, human response functions to toxicants, and a geographical information system are all integrated to create this environment. This simulation environment can help safety managers estimate the population at risk in order to plan evacuation, make sheltering decisions, determine the resources required to provide aid and comfort, and mitigate damages in case of a disaster. This simulation environment may also be modified and used for the landing phase of a space vehicle but will not be discussed in this paper.

Rabelo, Luis↗
NASA NTRS2005

Range Systems Simulation for the NASA Shuttle: Emphasis on Disaster and Prevention Management During Lift-Off

This article describes a decision-making system composed of a number of safety and environmental models for the launch phase of a NASA Space Shuttle mission. The components of this distributed simulation environment represent the different systems that must collaborate to establish the Expectation of Casualties (E(sub c)) caused by a failed Space Shuttle launch and subsequent explosion (accidental or instructed) of the spacecraft shortly after liftoff. This decision-making tool employs Space Shuttle reliability models, trajectory models, a blast model, weather dissemination systems, population models, amount and type of toxicants, gas dispersion models, human response functions to toxicants, and a geographical information system. Since one of the important features of this proposed simulation environment is to measure blast, toxic, and debris effects, the clear benefits is that it can help safety managers not only estimate the population at risk, but also to help plan evacuations, make sheltering decisions, establish the resources required to provide aid and comfort, and mitigate damages in case of a disaster.

Rabelo, Lisa↗
NASA NTRS2011

Analysis of Pulsed Airborne Lidar Measurements of Atmospheric CO2 Column Absorption from 3-13 km Altitudes

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS space mission [1]. It uses two pulsed laser transmitters allowing simultaneous measurement of a CO2 absorption line in the 1575 nm band, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are precisely stepped in wavelength across the CO2 line and an O2 line region during the measurement. The direct detection receiver measures the energies of the laser echoes from the surface along with the range profile of scattering in the path. The column densities for the CO2 and O2 gases are estimated from the ratio of the on- and off-line signals via the integrated path differential absorption (IPDA) technique. The time of flight of the laser pulses is used to estimate the height of the scattering surface and to reject laser photons scattered in the atmosphere. We developed an airborne lidar to demonstrate an early version of the CO2 measurement from the NASA Glenn Lear-25 aircraft. The airborne lidar stepped the pulsed laser's wavelength across the selected CO2 line with 20 wavelength steps per scan. The line scan rate is 450 Hz, the laser pulse widths are 1 usec, and laser pulse energy is 24 uJ. The time resolved laser backscatter is collected by a 20 cm telescope, detected by a NIR photomultiplier and is recorded on every other reading by a photon counting system [2]. During August 2009 we made a series of 2.5 hour long flights and measured the atmospheric CO2 absorption and line shapes using the 1572.33 nm CO2 line. Measurements were made at stepped altitudes from 3-13 km over locations in the US, including the SGP ARM site in Oklahoma, central Illinois, north-eastern North Carolina, and over the Chesapeake Bay and the eastern shore of Virginia. Although the received signal energies were weaker than expected for ASCENDS, clear CO2 line shapes were observed at all altitudes, and some measurements were made through thin clouds. The Oklahoma and east coast flights were coordinated with a LaRC/ITT CO2 lidar on the LaRC UC-12 aircraft, and in-situ measurements were made using its CO2 sensor and radiosondes. We have conducted an analysis of the ranging and IPDA lidar measurements from these four flights. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We used a cross-correlation approach to process the laser echo records. This was used to estimate the range to the scattering surface, to define the edges of the laser pulses and to determine echo pulse energy at each wavelength. We used a minimum mean square approach to fit an instrument response function and to solve for the best-fit CO2 absorption line shape. We then calculated the differential optical depth (DOD) of the fitted CO2 line. We computed its statistics at the various altitude steps, and compare them to the DODs calculated from spectroscopy based on HITRAN 2008 and the column conditions calculated from the airborne in-situ readings. The results show the lidar and in-situ measurements have very similar DOD change with altitude and greater than 10 segments per flight where the scatter in the lidar measurements are less than or equal to 1ppm. We also present the results from subsequent CO2 column absorption measurements, which were made with stronger detected signals during three flights on the NASA DC-8 over the southwestern US in during July 2010.

Abshire, James B.↗
NASA NTRS2012

System Identification and Uncertainty Quantification Using Orthogonal Excitations and the Semi-span Super Sonic Transport (S4T) Model

Orthogonal harmonic multisine excitations were utilized in a wind tunnel test and in simulation of the SemiSpan Supersonic Transport model to assess aeroservoelastic characteristics. Fundamental issues associated with analyzing sinusoidal signals were examined, including spectral leakage, excitation truncation, and uncertainties on frequency response functions and mean-square coherence. Simulation allowed for evaluation of these issues relative to a truth model, while wind tunnel data introduced real-world implementation issues.

Heeg, Jennifer↗
NASA NTRS2012

An Alternative Inter-Satellite Calibration of the UMD HIRS OLR Retrievals

Outgoing Longwave Radiation (OLR) at the top-of-atmosphere (TOA) is a fundamental component of Earth's energy balance and represents the heat energy in the thermal bands rejected to space by the planet. Determination of OLR from satellites has a long and storied history, but the observational record remains largely fragmented with gaps in satellite measurements over the past three decades. Perhaps the most semi-continuous set of retrievals comes from the University of Maryland (UMD) algorithm that uses four HIRS (High Resolution Infrared Sounder) channels on the NOAA polar orbiting satellites to estimate OLR. This data set shows great promise in helping to bridge the discontinuous ERBS (Earth Radiation Budget Satellite) and CERES (Clouds and the Earth s Radiant Energy System) measurements. However, significant satellite inter-calibration biases persist with the present UMD data, principally outside the tropics. Difficulties relate to the combination of drift of the satellite equator crossing time through the diurnal cycle and changes in HIRS channel response function design. Here we show how an ad hoc recalibration of the UMD retrievals among the different satellites removes much of the remaining uncertainty due to diurnal drift of the satellite orbit. The adjusted HIRS data (using no other external information) show much better agreement with OLR from the European Center Interim Reanalysis (EC-Int), longer-term signals in the Global Energy and Water Cycle Experiment / Surface Radiation Budget (GEWEX/SRB) retrievals, and also agree well with ERBS and CERES OLR measurements. These results augur well for narrowing the uncertainties in multi-decadal estimates of this important climate variable.

Robertson, Franklin R.↗
NASA NTRS2013

FUN3D Analyses in Support of the First Aeroelastic Prediction Workshop

This paper presents the computational aeroelastic results generated in support of the first Aeroelastic Prediction Workshop for the Benchmark Supercritical Wing (BSCW) and the HIgh REynolds Number AeroStructural Dynamics (HIRENASD) configurations and compares them to the experimental data. The computational results are obtained using FUN3D, an unstructured grid Reynolds-averaged Navier-Stokes solver developed at NASA Langley Research Center. The analysis results for both configurations include aerodynamic coefficients and surface pressures obtained for steady-state or static aeroelastic equilibrium (BSCW and HIRENASD, respectively) and for unsteady flow due to a pitching wing (BSCW) or modally-excited wing (HIRENASD). Frequency response functions of the pressure coefficients with respect to displacement are computed and compared with the experimental data. For the BSCW, the shock location is computed aft of the experimentally-located shock position. The pressure distribution upstream of this shock is in excellent agreement with the experimental data, but the pressure downstream of the shock in the separated flow region does not match as well. For HIRENASD, very good agreement between the numerical results and the experimental data is observed at the mid-span wing locations.

Chwalowski, Pawel↗
NASA NTRS2011

Estimation of Time-Varying Pilot Model Parameters

Human control behavior is rarely completely stationary over time due to fatigue or loss of attention. In addition, there are many control tasks for which human operators need to adapt their control strategy to vehicle dynamics that vary in time. In previous studies on the identification of time-varying pilot control behavior wavelets were used to estimate the time-varying frequency response functions. However, the estimation of time-varying pilot model parameters was not considered. Estimating these parameters can be a valuable tool for the quantification of different aspects of human time-varying manual control. This paper presents two methods for the estimation of time-varying pilot model parameters, a two-step method using wavelets and a windowed maximum likelihood estimation method. The methods are evaluated using simulations of a closed-loop control task with time-varying pilot equalization and vehicle dynamics. Simulations are performed with and without remnant. Both methods give accurate results when no pilot remnant is present. The wavelet transform is very sensitive to measurement noise, resulting in inaccurate parameter estimates when considerable pilot remnant is present. Maximum likelihood estimation is less sensitive to pilot remnant, but cannot detect fast changes in pilot control behavior.

Zaal, Peter M. T.↗
NASA NTRS2013

Understanding How Kurtosis Is Transferred from Input Acceleration to Stress Response and Its Influence on Fatigue Llife

High cycle fatigue of metals typically occurs through long term exposure to time varying loads which, although modest in amplitude, give rise to microscopic cracks that can ultimately propagate to failure. The fatigue life of a component is primarily dependent on the stress amplitude response at critical failure locations. For most vibration tests, it is common to assume a Gaussian distribution of both the input acceleration and stress response. In real life, however, it is common to experience non-Gaussian acceleration input, and this can cause the response to be non-Gaussian. Examples of non-Gaussian loads include road irregularities such as potholes in the automotive world or turbulent boundary layer pressure fluctuations for the aerospace sector or more generally wind, wave or high amplitude acoustic loads. The paper first reviews some of the methods used to generate non-Gaussian excitation signals with a given power spectral density and kurtosis. The kurtosis of the response is examined once the signal is passed through a linear time invariant system. Finally an algorithm is presented that determines the output kurtosis based upon the input kurtosis, the input power spectral density and the frequency response function of the system. The algorithm is validated using numerical simulations. Direct applications of these results include improved fatigue life estimations and a method to accelerate shaker tests by generating high kurtosis, non-Gaussian drive signals.

Kihm, Frederic↗
NASA NTRS2013

Anti­-parallel Filament Flows and Bright Dots Observed in the EUV with Hi-­C

Hi-C obtained the highest spatial and temporal resolution observations ever taken in the solar EUV corona. Hi-C reveals dynamics and structure at the limit of its temporal and spatial resolution. Hi-C observed various fine-scale features that SDO/AIA could not pick out. For the first time in the corona, Hi-C revealed magnetic braiding and component reconnection consistent with coronal heating. Hi-C shows evidence of reconnection and heating in several different regions and magnetic configurations with plasma being heated to 0.3 - 8 x 10(exp 6) K temperatures. Surprisingly, many of the first results highlight plasma at temperatures that are not at the peak of the response functions.

Alexander, Caroline E.↗
NASA NTRS2014

Combustion Stability Assessments of the Black Brant Solid Rocket Motor

The Black Brant variation of the Standard Brant developed in the 1960's has been a workhorse motor of the NASA Sounding Rocket Project Office (SRPO) since the 1970's. In March 2012, the Black Brant Mk1 used on mission 36.277 experienced combustion instability during a flight at White Sands Missile Range, the third event in the last four years, the first occurring in November, 2009, the second in April 2010. After the 2010 event the program has been increasing the motor's throat diameter post-delivery with the goal of lowering the chamber pressure and increasing the margin against combustion instability. During the most recent combustion instability event, the vibrations exceeded the qualification levels for the Flight Termination System. The present study utilizes data generated from T-burner testing of multiple Black Brant propellants at the Naval Air Warfare Center at China Lake, to improve the combustion stability predictions for the Black Brant Mk1 and to generate new predictions for the Mk2. Three unique one dimensional (1-D) stability models were generated, representing distinct Black Brant flights, two of which experienced instabilities. The individual models allowed for comparison of stability characteristics between various nozzle configurations. A long standing "rule of thumb" states that increased stability margin is gained by increasing the throat diameter. In contradiction to this experience based rule, the analysis shows that little or no margin is gained from a larger throat diameter. The present analysis demonstrates competing effects resulting from an increased throat diameter accompanying a large response function. As is expected, more acoustic energy was expelled through the nozzle, but conversely more acoustic energy was generated due to larger gas velocities near the propellant surfaces.

Fischbach, Sean↗