SEARCH · Search NASA

Results for “Central Pattern Generators”

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.

126 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Scattering and the Point Spread Function of the New Generation Space Telescope

Preliminary design work on the New Generation Space Telescope (NGST) is currently under way. This telescope is envisioned as a lightweight, deployable Cassegrain reflector with an aperture of 8 meters, and an effective focal length of 80 meters. It is to be folded into a small-diameter package for launch by an Atlas booster, and unfolded in orbit. The primary is to consist of an octagon with a hole at the center, and with eight segments arranged in a flower petal configuration about the octagon. The comers of the petal-shaped segments are to be trimmed so that the package will fit atop the Atlas booster. This mirror, along with its secondary will focus the light from a point source into an image which is spread from a point by diffraction effects, figure errors, and scattering of light from the surface. The distribution of light in the image of a point source is called a point spread function (PSF). The obstruction of the incident light by the secondary mirror and its support structure, the trimmed corners of the petals, and the grooves between the segments all cause the diffraction pattern characterizing an ideal point spread function to be changed, with the trimmed comers causing the rings of the Airy pattern to become broken up, and the linear grooves causing diffraction spikes running radially away from the central spot, or Airy disk. Any figure errors the mirror segments may have, or any errors in aligning the petals with the central octagon will also spread the light out from the ideal point spread function. A point spread function for a mirror the size of the NGST and having an incident wavelength of 900 nm is considered. Most of the light is confined in a circle with a diameter of 0.05 arc seconds. The ring pattern ranges in intensity from 10(exp -2) near the center to 10(exp -6) near the edge of the plotted field, and can be clearly discerned in a log plot of the intensity. The total fraction of the light scattered from this point spread function is called the total integrated scattering (TIS), and the fraction remaining is called the Strehl ratio. The angular distribution of the scattered light is called the angle resolved scattering (ARS), and it shows a strong spike centered on a scattering angle of zero, and a broad , less intense distribution at larger angles. It is this scattered light, and its effect on the point spread function which is the focus of this study.

Schreur, Julian J.↗

DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

NASA NTRS2019

Overexpression of Catalase in Mitochondria Mitigates the Effects of Simulated Microgravity and Social Isolation on Cytokine Expression in Mouse Hippocampus

Space flight missions are becoming longer and more common and evidence points to the physiological toll the missions have on the human body. Aging, sedentary lifestyle, and spaceflight have similar degenerative effects on almost every part of our body; both exposure to the space environment and aging result in cardiovascular deconditioning, bone loss, muscle atrophy, brain changes, and immune response impairment. We hypothesize that exposure to the space environment generates excessive Reactive Oxygen Species (ROS), which results in neuroinflammation and aging-like degenerative symptoms in the brain. We used the hindlimb unloading (HU) model to mimic microgravity with either paired or single housed animals (social isolation). Responses to 30d (30 days) of HU were compared in wildtype or transgenic MCAT (Mitochondrial-targeted CATalase) mice, in which mitochondrial ROS is quenched by over-expression of human catalase. Expression of 4-Hydroxynonenal (4HNE) and Park7 (redox-sensitive chaperone and sensor of oxidative stress) were measured by ELISA (Enzyme-Linked ImmunoSorbent Assay), a protein array quantified from the hippocampal cytokines and 8-hydroxy-2'-deoxyguanosine in serum was measured by ELISA to assess oxidative DNA damage. Preliminary analysis of cage behavior patterns from video collected at the end of the study showed that MCAT HU mice (socially housed) were more active and conducted more exploratory activities compared to NL (Normally Loaded). Our biochemical results showed simulated microgravity and/or social isolation caused changes in levels of cytokines related to immune responses. Two-way ANOVA (Analysis of Variance) revealed significant interaction effects of HU and genotype in expression levels of five cytokines (out of 35) in socially-housed animals. Elevation of these generally pro-inflammatory cytokines by HU in WT (Wild Type) mice was mitigated in MCAT mice, suggesting a role for mitochondrial ROS signaling in inflammatory CNS (Central Nervous System) responses to microgravity. Interestingly, some of these cytokines in the hippocampus displayed strong correlations to the 4HNE levels. We also found substantive cytokine responses to social isolation in the hippocampus; housing and genotype interaction effects were significant (by 2-Factor ANOVA) for 15 cytokines, most of which were mitigated in MCAT mice.Taken together, our results showed that both simulated microgravity and social isolation influenced cytokine levels in the hippocampus and MCAT mice were at least partially protected from these changes. These findings implicate a potentially important role for mitochondrial ROS in CNS responses to the challenges posed by long duration spaceflight.

MCAT↗

DOE OSTI2022

Turbulent transport and reactions of plant-emitted hydrocarbons in an Amazonian rain forest

The processes governing the temporal and spatial patterns of isoprene and monoterpenes emitted by a rainforest in the central Amazon region of Brazil is investigated using a combination of field experiments and numerical simulations. Specifically, Large Eddy Simulations (LES) were used to resolve emissions of isoprene and monoterpenes, turbulent transport, and air chemistry. The coupled chemistry-transport LES included the effects of isoprene and monoterpenes reactivity due to reactions with hydroxyl radical and ozone. The LES results are used to compute vertically resolved budgets of isoprene and monoterpenes in the rainforest canopy in response to emissions, turbulent transport, surface deposition, and air chemistry. Our results indicated that emission and dispersion dominated the isoprene budget as the gases were transported out of the canopy space. In a region limited by nitrogen oxides (with prevailing nitric oxide levels of < 0.5 parts per billion), the in-canopy chemical destruction removed approximately 10% of locally emitted monoterpenes. Hydroxyl radical production rates from the ozonolysis of monoterpenes amounted to ≈ 2 x 10 6 radicals cm -3 s -1 and had similar magnitude to the light-dependent hydroxyl radical formation. One key conclusion was that the Amazonia rainforest abundantly emitted monoterpenes whose in-canopy ozonolysis yielded hydroxyl radicals in amounts similar to the magnitude of light-dependent formation. Reactions of monoterpenes and isoprene with hydroxyl radical and ozone were necessary for the maintenance of the Amazon rainforest canopy as a photochemically active environment suitable to generate oxidants and secondary organic aerosols.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Megadrought: A Series of Unfortunate La Niña Events?

Megadroughts are multidecadal periods of aridity more persistent than most droughts during the instrumental period. Paleoclimate evidence suggests that megadroughts occur in many parts of the world, including North America, Central America, western Europe, eastern Asia, and northern Africa. It remains unclear to what extent such megadroughts require external forcing or whether they can arise from internal climate variability alone. A novel statistical–dynamical approach is used to evaluate the possibility that such events arise solely as a function of interannual tropical sea surface temperature (SST) variations. A statistical emulator of tropical SST variations is constructed by using an empirical moving-blocks bootstrap approach that randomly samples multiyear sequences of the observational SST record. This approach preserves the power spectrum, seasonal cycle, and spatial pattern of El Niño-Southern Oscillation (ENSO) but removes longer timescale fluctuations embedded in the observational record. These resampled SST anomalies are then used to force an atmospheric model (the Community Atmosphere Model Version 5). As megadroughts emerge in this run, they should, therefore, be solely a consequence of La Niña sequences combined with internal atmospheric variability and persistence driven by soil moisture storage and other land-surface processes. We indeed find that megadroughts in this simulation have an amplitude-duration rate that is generally indistinguishable from the rate documented in paleoclimate records of the western United States. Our findings support the idea that megadroughts may occur randomly when the unforced climate system evolves freely over a sufficiently long period of time, implying that an unforced unusual but statistically plausible series of La Niña events may be sufficient to generate megadrought.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Development of Thin Gap GEM-µRWELL Hybrid Detectors at Jefferson Lab

Over the past few decades, Micro Pattern Gaseous Detector (MPGD) technologies have been increasingly adopted as tracking detector options in High Energy and Nuclear Physics experiments thanks to their good spatial resolution, high-rate capability, stability and more importantly their ability for large area coverage at a relatively low cost compared to the alternative. The thin gap GEM-µRWELL hybrid detector is the latest addition to the MPGD family, that was introduced to vastly improve the spatial resolution capability of gaseous trackers when deployed in the barrel region to cover large angular acceptance of the central tracker in a collider experiment. In this talk, I will re-introduce the concept and motivation for the development of thin gap GEM-µRWELL hybrid technology with an emphasis on the initial studies that establish the proof-of-concept of the technology. I will then discuss the more recent results from latest beam test campaign at Jefferson Lab in May 2025 to study detector efficiency performance with various gas mixtures. I will also briefly present the ongoing activities to develop large area thin gap GEM-µRWELL tracking detectors for the ePIC experiment of the future Electron Ion Collider as well as the exploration of the technology to provide large area tracking options to the muon system of experiments at a future Higgs Factory Collider such as the FCC-ee for example. Finally, I will conclude with some perspectives on new ideas under exploration to develop the next generation of thin gap MPGD technologies with enhanced timing and spatial resolution capabilities

Gnanvo, Kondo [Thomas Jefferson National Accelerat↗

DOE OSTI2021

Shape-Selective Synthesis of Alkylcyclopentenyl Cations in Zeolites and Spectroscopic Distinction of Constitutional Isomers

Alkylcyclopentenyl cations belong to the long-lived intermediates that make up the “hydrocarbon pool” during the catalytic conversion of methanol on zeolites, and recent works show that such cations contribute to olefin and aromatics formation. From liquid phase chemistry, two types of alkylcyclopentenyl cations are known and distinguished by the substituent at the central carbon (C-2) of the allylic system: a more stable type with a methyl group at the C-2 and a less stable type with hydrogen at the C-2. The following three linked objectives are pursued: (i) IR spectroscopic distinction of the different substitution patterns at the allylic system of alkylcyclopentenyl cations, (ii) the role of the zeolite framework in determining the substitution pattern, and (iii) identification of alkylcyclopentenyl cations from precursors relevant to methanol-to-olefins conversion. UV–vis and IR spectroscopy are applied in situ to characterize alkylcyclopentenyl cations produced by adsorption of a pentaalkylcyclopentadiene or by adsorption and thermally induced cyclization of 2,6-dimethyl-2,4,6-octatriene. Prior knowledge of electronic spectra is combined with DFT-computed and experimental IR spectra to establish the frequency of the C–H vibration of the C-2 hydrogen-substituted type, and a characteristic red shift of the asymmetric allylic stretching vibration (Δν ≈ -20 to -30 cm –1 ) after replacing hydrogen by methyl at the central carbon of the allylic system. Although DFT demonstrates that both types of ions fit into medium- and large-pore zeolites and that the C-2 methyl-substituted type is thermodynamically favored even in the pores of the considered zeolites, formation of alkylcyclopentenyl cations is found by UV–vis and IR spectra to be shape-selective. The bulkier C-2 methyl-substituted type is detected in large-pore zeolites (MOR, BEA) and in the intersections of medium-pore zeolites (MFI), whereas in channels of medium size (TON), the less bulky C-2 hydrogen-substituted type is exclusively formed. Finally, the type of ion formed and its overall size are dictated by the zeolite framework and are independent of the precursor; the same type of alkylcyclopentenyl cation as found through cyclization of dimethyloctatriene could be generated from ethene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1997

Active Microwave Remote Sensing Observations of Weddell Sea Ice

Since July 1991, the European Space Agency's ERS-1 and ERS-2 satellites have acquired radar data of the Weddell Sea, Antarctica. The Active Microwave Instrument on board ERS has two modes; SAR and Scatterometer. Two receiving stations enable direct downlink and recording of high bit-rate, high resolution SAR image data of this region. When not in an imaging mode, when direct SAR downlink is not possible, or when a receiving station is inoperable, the latter mode allows normalized radar cross-section data to be acquired. These low bit-rate ERS scatterometer data are tape recorded, downlinked and processed off-line. Recent advances in image generation from Scatterometer backscatter measurements enable complementary medium-scale resolution images to be made during periods when SAR images cannot be acquired. Together, these combined C-band microwave image data have for the first time enabled uninterrupted night and day coverage of the Weddell Sea region at both high (25 m) and medium-scale (-20 km) resolutions. C-band ERS-1 radar data are analyzed in conjunction with field data from two simultaneous field experiments in 1992. Satellite radar signature data are compared with shipborne radar data to extract a regional and seasonal signature database for recognition of ice types in the images. Performance of automated sea-ice tracking algorithms is tested on Antarctic data to evaluate their success. Examples demonstrate that both winter and summer ice can be effectively tracked. The kinematics of the main ice zones within the Weddell Sea are illustrated, together with the complementary time-dependencies in their radar signatures. Time-series of satellite images are used to illustrate the development of the Weddell Sea ice cover from its austral summer minimum (February) to its winter maximum (September). The combination of time-dependent microwave signatures and ice dynamics tracking enable various drift regimes to be defined which relate closely to the circulation of the sea ice in response to current and wind forcing and iceberg barriers. These are closely related to continental-shelf or central basin regimes, in which tidal forcing or barotropic circulation patterns appear to influence the sea-ice motion, respectively. These regimes provide valuable information about the regions of most prolific ice growth and influence of ice conditions upon air-sea-ice exchange processes in the Weddell Sea.

Drinkwater, Mark R.↗

DOE OSTI2023

Inclusive jet and hadron suppression in a multistage approach

We present a new study of jet interactions in the quark-gluon plasma created in high-energy heavy-ion collisions, using a multistage event generator within the jetscape framework. We focus on medium-induced modifications in the rate of inclusive jets and high transverse momentum (high-p T ) hadrons. Scattering-induced jet energy loss is calculated in two stages: a high virtuality stage based on the matter model, in which scattering of highly virtual partons modifies the vacuum radiation pattern, and a second stage at lower jet virtuality based on the LBT model, in which leading partons gain and lose virtuality by scattering and radiation. Coherence effects that reduce the medium-induced emission rate in the matter phase are also included. The TRENTO model is used for initial conditions, and the (2 + 1)dimensional VISHNU model is used for viscous hydrodynamic evolution. Jet interactions with the medium are modeled via 2-to-2 scattering with Debye screened potentials, in which the recoiling partons are tracked, hadronized, and included in the jet clustering. Holes left in the medium are also tracked and subtracted to conserve transverse momentum. Calculations of the nuclear modification factor (R AA ) for inclusive jets and high-p T hadrons are compared to experimental measurements at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC). Within this framework, we find that with one extra parameter which codifies the transition between stages of jet modification—along with the typical parameters such as the coupling in the medium, the start and stop criteria, etc.—we can describe these data at all energies for central and semicentral collisions without a rescaling of the jet transport coefficient q.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Mitigation of ICRF—edge plasma interaction in Alcator C-Mod

The parasitic interaction between ion cyclotron range of frequencies (ICRF) waves and the edge plasma has long limited high-power ICRF utilization in magnetically confined fusion plasmas. In this paper, we report on the successful mitigation of this interaction by power tapering a four-strap field-aligned (FA) antenna in the Alcator C-Mod tokamak. With the antenna operated in dipole phasing ($0\pi0\pi$), the ratio of the power coupled by the central two straps ($P_\mathrm {cent}$) to the power coupled by the outer two straps was varied at fixed total coupled ICRF power ($P_\mathrm {tot}$). With $P_\mathrm {cent}$/$P_\mathrm {tot} \approx 0.8-0.9$, no enhancement of the plasma potential was measured despite 1 MW of coupled ICRF power. When the antenna is powered, energy deposition increases on its four corners and on neighboring limiter tiles magnetically connected to some of those corners. The energy deposited on the antenna and limiter is minimized together with near-field plasma potentials. Additionally, we show that deviation from $P_\mathrm {cent}$/$P_\mathrm {tot} \approx 0.8$ enhances wave-plasma interactions in places not magnetically connected to the antenna. These far-field effects are correlated with changes in the antenna spectrum. Hence, on C-Mod, near-field and far-field interactions were minimized for about the same range of $P_\mathrm {cent}$/$P_\mathrm {tot}$. Optimal power tapering was also found to improve ICRF heating efficiency and facilitate H-mode access. Finally, we show that the heat flux pattern on the toroidally aligned antennas reverses with the direction of the RF-induced E×B flow in front of the antenna, a correlation not observed for the FA antenna. We argue that the RF-induced E×B flow, which carries high-density plasma from the near-SOL, is aligned with the FA antenna and does not intercept its corners, where the plasma potential is large. This could explain the previously reported smaller impurity generation by the FA antenna (Wukitch et al 2013 Phys. Plasmas 20 056117).

Alcator C-Mod↗

NASA NTRS1995

ROSAT observations of NGC 2146: Evidence for a starburst-driven superwind

We have imaged the edge-on starburst galaxy NGC 2146 with the Position Sensitive Proportional Counter (PSPC) and the High Resolution Imager (HRI) on board ROSAT and have compared these data to optical images and long-slit spectra. NGC 2146 possesses a very large X-ray nebula with a half-light radius of 1 min (4 kpc) and a maximum diameter of approximately 4 min, or 17 kpc. The X-ray emission is resolved by the PSPC and preferentially oriented along the minor axis, with a total flux of 1.1 x 10(exp -12) ergs/sq cm/s over 0.2 - 2.4 keV and a luminosity of approximately 3 x 10(exp 40) ergs/s. The inner X-ray nebula is resolved by the HRI into at least four bright knots together with strong diffuse emission responsible for at least 50% of the flux within a radius of 0.5 min (approximately 2 kpc). The brightest knot has a luminosity of (2 - 3) x 10(exp 39) ergs/s. The X-ray nebula has a spatial extent much larger than the starburst ridge seen at centimeter wavelengths by Kronberg & Biermann (1981) and is oriented in a `X-like' pattern along the galaxy minor axis at a position angle of approximately 30 degrees. This minor-axis X-ray emission is associated with a region of H alpha and dust filaments seen in optical images. Optical spectra show that the emission-line gas along the minor axis is characterized by relatively broad lines (approximately 250 km/s full width half-maximum (FWHM)) and by `shocklike' emission-line flux ratios. Together with the blue-asymmetric nuclear emission-line and NaD interstellar absorption-line profiles, these optical data strongly suggest the presence of a starburst-driven superwind. The X-ray spectrum extracted from the central 5 min contains a strong Fe L emission-line complex at 0.6 - 1.0 keV and a hard excess above 1.0 keV. The spectrum is best described with a two-component model, containing a soft (kT approximately 400 - 500 eV) Raymond-Smith thermal plasma together with either a Gamma = 1.7 power-law or a kT greater than 2.2 keV bremsstrahlung component. The soft thermal component provides approximately 30% of the total luminosity over 0.2 - 2.4 keV, or approximately 10(exp 40) ergs/s. The pressure derived from the soft component of the X-ray spectrum is consistent with that predicted from a starburst-driven superwind if the filling factor of the warm gas is approximately 1% - 10 %. If the hard X-ray component is thermal gas associated with the galactic outflow, the filling factor must be close to unity. Predictions of the luminosity, temperature, and size of an adiabatic starburst-generated windblown bubble are consistent with those measured for the soft thermal X-ray emission in NGC 2146. The hard X-ray component, however, has a luminosity much larger than predicted by the superwind model if this component is thermal emission from gas heated by an internal shock in the expanding bubble. We briefly review various possibilities as to the nature of the hard X-ray component in NGC 2146.

Armus, L.↗

DOE OSTI2019

Understanding the Response of Photosynthetic Metabolism in Tropical Forests to Seasonal Climate Variations

We propose a study focused on the basic question: What controls the response of photosynthesis in Amazonian forests to seasonal variations in climate? This question, despite its apparent simplicity, remains difficult for modern earth system models to answer, and is also the subject of continuing controversy in the remote-sensing literature. For example, in the modeling arena, four cutting-edge earth system models (ESM’s) show significant divergence in their seasonal patterns of photosynthesis from observed whole-system photosynthetic fluxes at two sites in the central Amazon. The overall objective of the project is to guide improvements in earth system models of tropical forest photosynthesis by collecting and integrating a suite of observations to 1) test several hypotheses (three core, conceptual hypotheses for explaining observed photosynthetic seasonality and a methodological hypotheses for scaling from leaves to canopy with hyperspectral cameras), and 2) perform a synthesis activity that applies our empirical work to earth system models of terrestrial carbon cycling. The project will provide an extensive suite of new and unique datasets that enable us to fill, through advanced modeling techniques and analysis, critical knowledge gaps in current understanding of what controls the response of canopy photosynthesis and related functions in Amazonian forests to seasonal variation in climate. Three major types of datasets and data products will be delivered: (1) in situ leaf and tree-scale measurements from intensive ecophysiological and ecohydological field campaigns, (2) time-series observations of leaf-to-crown scale forest reflectance properties and atmospheric radiation from two innovative, ground-based imaging sensors (respectively, the Hyperspectral Vegetation Imaging System and the High Dynamic Range All-Sky Imaging System), and (3) results from state-of-the-art models of 3-dimensional canopy processes for radiative transfer and photosynthesis that integrate and link our observations to tropical forest processes. These data products and the improved knowledge we achieve with them will contribute to testing and improving the treatment of tropical forest processes in ESMs. They will contribute data to and leverage related data from of the GOAmazon campaign, and make significant contributions to support the overall goals of GOAmazon. This work will also help establish a foundation for the Next Generation Ecosystem Experiments (NGEE) in the Tropics.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1984

Case studies of NOAA 6/TIROS N data impact on numerical weather forecasts

The impact of satellite temperatures from systems which predate the launching of the third generation of vertical sounding instruments aboard TIROS N (13 Oct 1978) and NOAA 6 (27 June 1979) is reported. The first evaluation of soundings from TIROS N found that oceanic, cloudy retrievals over NH mid latitudes show a cold bias in winter. It is confirmed for both satellite systems using a larger data base. It is shown that RMS differences between retrievals and colocated radiosonde observations within the swath 30-60N during the 1979-80 winter were generally 2-3K in clear air and higher for cloudy columns. A positive impact of TIROS N temperatures on the analysis of synoptic weather systems is shown. Analyses prepared from only satellite temperatures seemed to give a better definition to weather systems' thermal structure than that provided by corresponding NMC analyses without satellite data. The results of a set of 14 numerical forecast experiments performed with the PE model of the Israel Meteorological Service (IMS) are summarized; these were designed to test the impact of TIROS N and NOAA 6 temperatures within the IMS analysis and forecast cycle. The satellite data coverage over the NH, the mean area/period S1 and RMS verification scores and the spatial distribution of SAT versus NO SAT forecast differences are discussed and it is concluded that positive forecast impact occurs over ocean areas where the extra data improve the specification which is otherwise available from conventional observations. The forecast impact for three cases from the same set of experiments was examined and it is found that satellite temperatures, observed over the Atlantic Ocean contribute to better forecasts over Iceland and central Europe although a worse result was verified over Spain. It is also shown that the better scores of a forecast based also on satellite data and verified over North America actually represent a mixed impact on the forecast synoptic patterns. A superior 48 hr 500 mb forecast over the western US due to the better initial specification afforded by satellite observed temperatures over the North Pacific Ocean is shown.

Druyan, L. M.↗

NASA NTRS2001

Large Aperture Scanning Lidar Based on Holographic Optical Elements

Lidar remote sensing instruments can make a significant contribution to satisfying many of the required measurements of atmospheric and surface parameters for future spaceborne platforms, including topographic altimeters, atmospheric profiles of, wind, humidity, temperature, trace molecules, aerosols, and clouds. It is highly desirable to have wide measurement swaths for rapid coverage rather than just the narrow ribbon of data that is obtained with a nadir only observation. For most applications global coverage is required, and for wind measurements scanning or pointing is required in order to retrieve the full 3-D wind vector from multiple line-of-sight Doppler measurements. Conventional lidar receivers make up a substantial portion of the instrument's size and weight. Wide angle scanning typically requires a large scanning mirror in front of the receiver telescope, or pointing the entire telescope and aft optics assembly, Either of these methods entails the use of large bearings, motors, gearing and their associated electronics. Spaceborne instruments also need reaction wheels to counter the torque applied to the spacecraft by these motions. NASA has developed simplified conical scanning telescopes using Holographic Optical Elements (HOEs) to reduce the size, mass, angular momentum, and cost of scanning lidar systems. NASA has developed two operating lidar systems based on 40 cm diameter HOEs. The first such system, named Prototype Holographic Atmospheric Scanner for Environmental Remote Sensing (PHASERS) was a joint development between NASA Goddard Space Flight Center (GSFC) and the University of Maryland College Park. PHASERS is based on a reflection HOE for use at the doubled Nd:YAG laser wavelength of 532 nm and has recently undergone a number of design changes in a collaborative effort between GSFC and Saint Anselm College in New Hampshire. The next step was to develop IR transmission HOEs for use with the Nd:YAG fundamental in the Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE). The HOE spins like a compact disk in a large ring ball bearing. In an aircraft the HOE faces down, looking out through a window at an angle of 45 degrees off-nadir. The HOE diffracts 85% of the incident 532 nm light into a 160 micron spot at a focal length of 1 meter. HARLIE is a field deployable lidar measuring aerosol, cloud, and boundary layer backscatter for atmospheric research. It has flown several times and is also used from a ground-based trailer in an upward-looking mode. The HOE generates a 45 degree conical scan pattern by rotating at speeds up to 30 rpm. Like PHASERS, the HOE in HARLIE serves both as the laser collimating lens as well as the receiver telescope primary optic. The telescope is coupled to the receiver package via fiber optic. The transmitter is a diode pumped Nd:YAG laser operating at 1064 nm, delivering 1 mJ pulses at a 5 KHz rep-rate. The receiver has a 200 microradian field-of-view and a 0.5 nm optical bandpass. The photon counting data system utilizes a single Geiger-mode silicon avalanche photodiode detector, This new technology has also presented us with new data visualization challenges as well as new measurement techniques. The backscatter data obtained from a stationary (i.e. ground-based) scanning HOE lidar is on the surface of a cone, which when viewed over many consecutive scans can reveal atmospheric motions on this surface over time as the atmosphere advects over the site. In a moving platform such as an airplane or satellite, the data from consecutive scans cover different areas under the flight path, revealing atmospheric structure in 3-dimensions. An example of a visualization of HARLIE ground-based data is presented, showing aerosol backscatter on a 90 degree conical surface generated from one 360 degree scan of the lidar during the HOLO-1 field campaign on the afternoon of 10 March 1999. Higher backscatter levels are rendered as lighter signal against a dark background. Breaking Kelvin-Helmholtz waves are evident on the north side of the scan at an altitude of 10-11 km. Time series of successive scans made at regular intervals render unique views of atmospheric motions, from which vertical profiles of atmospheric wind vectors can be obtained using a unique data analysis approach. Wind vectors obtained from the lidar were compared with co-located radiosonde wind profiles during an intensive operating period in September-October 2000 at the Atmospheric Radiation Measurement Program's Southern Great Plains Central Facility.

Schwemmer, Geary K.↗

DOE OSTI2020

Climate-Related Trends of Within-Storm Intensities Using Dimensionless Temporal-Storm Distributions

Huff curves are probabilistic time distributions of rainfall expressed as dimensionless cumulative percentages of storm depth and duration. Previous studies have documented development factors, spatial robustness, and the utility of Huff curves in practical applications. However, the effects of trending rainfall on Huff curve intensity patterns have not yet been studied. As such, the goal of this paper is to fill this gap by studying Huff curve patterns in a watershed with demonstrated increasing trends of temperature and precipitation, with the intention that it can be generalized to other areas in the US and the world. To achieve this goal, the high temporal resolution precipitation data collected from a high spatial density, 72-year precipitation-gauge network on the 4.25-km 2 North Appalachian Experimental Watershed in east-central Ohio were used. Seasonal storm pattern trends from 1939 to 2010 were investigated using dimensionless depth (with the frequency of 50%, $d_{50}$) and the curve variability $(V = d_{80} - d_{20})$ at three dimensionless within-storms time periods (three verticals). The Spearman rank correlation procedure (correlation coefficient, ρ and significance probability, p) was used to statistically determine trends over time using 8 periods of 4-season sets of Huff curves over the 72 years. Two averages of ρ and p were computed: (1) by averaging the individual ρ and p obtained from 10 gauges (AvgI); and (2) by grouped averaging of all individual gauge values of $d_{50}$ and V and then computing ρ and p (AvgG). The test results of individual gauges showed that 4 cases for $d_{50}$ and 23 cases for V were significant for all seasons and verticals (total of 120 cases for each variable). The test results of AvgI for $d_{50}$ and V and AvgG for $d_{50}$ showed no significant trends in all seasons and verticals. Only the AvgG for V led to a significant trend for V in spring and fall at different times within storm patterns. The data do not provide sufficient evidence at the p=0.05 significance level to reject the null hypothesis of unchanging position of the dimensionless depth of the 50% Huff curves for individual or averages for all seasons and verticals. Also, there is insufficient evidence to reject the null hypothesis of unchanging variability, V, using the average ρ of individual gauges (AvgI for V) for all seasons and verticals. AvgG results showed a significant trend in V; however, this analysis may be affected by the nonindependence of storm data. The results suggest that it is likely that there is little if any effect of trending climate over approximately 70 years on Huff curve patterns. The results of this study add to the robustness characteristics of Huff curves and to their potential use in hydrological practice as design storms, as the foundation of stochastic storm generation, and for storm disaggregation, and they deserve further investigation. These different forms of inputs to watershed models have the potential to improve runoff estimation. The results suggest that, if verified in other studies, they have applicability to provide useful stationary precipitation patterns across the US and other areas of the world in areas of nonstationary climate. Also, individual rain gauge data may not be representative of trends even over small areas, and seasonal differences were noticeable as found in other studies. Recommendations are provided.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2011

Contributions of GRACE to Climate Monitoring

The NASA/German Gravity Recovery and Climate Experiment (GRACE) was launched in March 2002. Rather than looking downward, GRACE continuously monitors the locations of and precise distance between twin satellites which orbit in tandem about 200 km apart. Variations in mass near Earth's surface cause heterogeneities in its gravity field, which in turn affect the orbits of satellites. Thus scientists can use GRACE data to map Earth's gravity field with enough accuracy to discern month to month changes caused by ocean circulation and redistribution of water stored on and in the land. Other gravitational influences, such as atmospheric circulation, post-glacial rebound, and solid earth movements are either independently determined and removed or are negligible on a monthly to sub-decadal timescale. Despite its coarse spatial (>150,000 sq km at mid-latitudes) and temporal (approx monthly) resolutions, GRACE has enabled significant advancements in the oceanic, hydrologic, and cryospheric science, and has great potential for climate monitoring, because it is the only global observing system able to measure ocean bottom pressures, total terrestrial water storage, and ice mass changes. The best known GRACE results are estimates of Greenland and Antarctic ice sheet loss rates. Previously, scientists had estimated ice mass losses using ground and satellite based altimetry and surface mass balance estimates based on snowfall accumulation and glacier discharge. While such measurements are still very useful for their spatial detail, they are imperfectly correlated with large-scale ice mass changes, due to snow and ice compaction and incomplete spatial coverage. GRACE enables scientists to generate monthly time series of Greenland and Antarctic ice mass, which have confirmed the shrinking of the polar ice sheets, one of the most obvious and indisputable manifestations of climate change. Further, GRACE has located and quantified hot spots of ice loss in southeastern Greenland and western Antarctica. For 2002 to present, the rate of ice mass loss has been 200 to 300 GT/yr in Greenland and 70 to 210 GT/yr in Antarctica, and some scientists are suggesting that the rates are accelerating. Similarly, GRACE has been used to monitor mass changes in alpine glaciers. Tamisiea et al. first characterized glacier melt along the southern coast of Alaska, more recently estimated to be occurring at a rate of 84 GT/yr. Chen et al. estimated that Patagonian glaciers are melting at a rate of 28 GT/yr, and estimated that the high mountains of central Asia lose ice at a rate of 47 GT/yr. Tapley et al. and Wahr et al. presented the first GRACE based estimates of changes in column-integrated terrestrial water storage (TWS; the sum of ground-water, soil moisture, surface waters, snow, ice, and water stored in vegetation) at continental scales. Since then, dozens of studies have shown that GRACE based estimates of regional to continental scale TWS variations agree with independent information, and some innovative uses of GRACE data have been developed. Rodell et al. (2004) and Swenson and Wahr (2006) demonstrated that by combining GRACE derived terrestrial water storage changes with observations of precipitation and runoff in a river basin scale water budget, it was possible to produce new estimates of evapotranspiration and atmospheric moisture convergence, essential climate variables that are difficult to estimate accurately. Similarly, GRACE has been used to constrain estimates of global river discharge and the contribution of changes in TWS to sea level rise. Crowley et al. observed a negative correlation between interannual TWS anomalies in the Amazon and the Congo River basin. Yeh et al. and Rodell et al. estimated regionally averaged groundwater storage variations based on GRACE and auxiliary observations. Rodell et al. and Tiwari et al. applied that method to quantify massive groundwater depletion in northern India caused by over reliance on aquifers for irration, and Famiglietti et al. found a similar situation in California's Central Valley. Zaitchik et al. and Lo et al. described approaches to use GRACE to constrain hydrological models, enabling integration of GRACE data with other observations and achieving much higher spatial and temporal resolutions than GRACE alone. Such approaches are now supporting applications including drought and water resources monitoring. Oceanography has likewise benefitted from the independent nature of GRACE observations. One application is measurement of the mass component of sea level rise, which complements radar altimetry and in situ measurements. GRACE also measures ocean bottom pressures (OBP), which help to refine understanding and modeling of ocean circulation and the ocean's fresh water budget, among other things. For example, Hayakawa et al. showed that GRACE observes OBP patterns absent from the background models of oceanic variability. Morison et al. used GRACE to describe important decadal scale shifts in circulation and an ongoing trend of freshening of the western Arctic, important indicators of climate variability. The research of Song and Zlotnicki and Chambers and Willis on GRACE-derived ocean bottom pressures in the sub-polar gyre led to the discovery of an ENSO teleconnection and a long-term change in OBP in the North Pacific sub-polar gyre that was not predicted by an ocean model. Further, Chambers and Willis were able to identify an internal redistribution of mass between Atlantic and Pacific Oceans lasting at least six years, which was not predicted by ocean models and was the first direct evidence of sustained mass transport from one ocean basin to another on periods longer than a year. Boening et al. observed a record increase in OBP over part of the southeastern Pacific in late 2009 and early 2010, primarily caused by wind stress curl associated with a strong and persistent anticyclone and likely related to the concurrent Central Pacific El Nino. GRACE has far surpassed its 5-year design lifetime, but it will likely succumb to the aging of batteries and instrument systems sometime in the next few years. NASA has begun initial development of a follow-on to GRACE with very similar design, which could launch as soon as 2016 and would provide continuity in the data record while improving resolution slightly. Higher resolution time variable gravity missions are also on the drawing board .

Rodell, Matthew↗

NASA NTRS2004

Exhaust Nozzle for a Multitube Detonative Combustion Engine

An improved type of exhaust nozzle has been invented to help optimize the performances of multitube detonative combustion engines. The invention is applicable to both air-breathing and rocket engines used to propel some aircraft and spacecraft, respectively. In a detonative combustion engine, thrust is generated through the expulsion of combustion products from a detonation process in which combustion takes place in a reaction zone coupled to a shock wave. The combustion releases energy to sustain the shock wave, while the shock wave enhances the combustion in the reaction zone. The coupled shockwave/reaction zone, commonly referred to as a detonation, propagates through the reactants at very high speed . typically of the order of several thousands of feet per second (of the order of 1 km/s). The very high speed of the detonation forces combustion to occur very rapidly, thereby contributing to high thermodynamic efficiency. A detonative combustion engine of the type to which the present invention applies includes multiple parallel cylindrical combustion tubes, each closed at the front end and open at the rear end. Each tube is filled with a fuel/oxidizer mixture, and then a detonation wave is initiated at the closed end. The wave propagates rapidly through the fuel/oxidizer mixture, producing very high pressure due to the rapid combustion. The high pressure acting on the closed end of the tube contributes to forward thrust. When the detonation wave reaches the open end of the tube, it produces a blast wave, behind which the high-pressure combustion products are expelled from the tube. The process of filling each combustion tube with a detonable fuel/oxidizer mixture and then producing a detonation repeated rapidly to obtain repeated pulses of thrust. Moreover, the multiple combustion tubes are filled and fired in a repeating sequence. Hence, the pressure at the outlet of each combustion tube varies cyclically. A nozzle of the present invention channels the expansion of the pulsed combustion gases from the multiple combustion tubes into a common exhaust stream, in such a manner as to enhance performance in two ways: (1) It reduces the cyclic variations of pressure at the outlets of the combustion tubes so as to keep the pressure approximately constant near the optimum level needed for filling the tubes, regardless of atmospheric pressure at the altitude of operation; and (2) It maximizes the transfer of momentum from the exhaust gas to the engine, thereby maximizing thrust. The figure depicts a typical engine equipped with a nozzle according to the invention. The nozzle includes an interface section comprising multiple intake ports that couple the outlets of the combustion tubes to a common plenum. Proceeding from its upstream to its downstream end, the interface section tapers to a larger cross-sectional area for flow. This taper fosters expansion of the exhaust gases flowing from the outlets of the combustion tubes and contributes to the desired equalization of exhaust combustion pressure. The cross-sectional area for flow in the common plenum is greater than, or at least equal to, the combined cross-sectional flow areas of the combustor tubes. In the common plenum, the exhaust streams from the individual combustion tubes mix to form a single compound subsonic exhaust stream. Downstream of the common plenum is the throat that tapers to a smaller flow cross section. In this throat, the exhaust gases become compressed to form a compound sonic gas stream. Downstream of the throat is an expansion section, which typically has a bell or a conical shape. (The expansion section can be truncated or even eliminated in the case of an air-breathing engine.) After entering the expansion section, the exhaust gases expand rapidly from compound sonic to compound supersonic speeds and are then vented to the environment. The basic invention admits of numerous variations. For example, the combustion tubes can be arranged around the central axin a symmetrical or asymmetrical pattern other than the one shown in the figure. For another example, the flow cross-sectional area(s) of one or more of the intake ports in the interface section, of the common plenum, the throat, and/or the expansion section can be varied, either symmetrically or asymmetrically, to adjust dynamics of the exhaust stream or to direct the thrust vector away from the central axis.

Bratkovich, Thomas E.↗