Search NASA⌕ Search

SEARCH · Search NASA

Results for “Air Dispersion Modeling”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Baltimore Energy & Infrastructure: Assessing Urban Heat Vulnerability in Baltimore Neighborhoods to Inform Transportation Resiliency Planning Efforts

Across the world, climate change is altering cities and their local climate systems. Following global trends of rising mean temperatures, Baltimore, Maryland is projected to experience more frequent extreme heat events. Exacerbated by Urban Heat Island (UHI) effects, residents of the city face high temperatures from limited tree canopy, highly paved and impervious surfaces, limited air flow, and high concentrations of localized emissions. Urban heat is dispersed asymmetrically, with historically marginalized populations experiencing disproportionately severe heat. Recognizing the impacts of extreme urban heat, the Maryland Transit Administration (MTA) collaborated with us, the Spring 2024 Goddard Energy and Infrastructure NASA DEVELOP team to assess urban impacts of MTA users and assets and explore opportunities to enhance resiliency planning. We evaluated the feasibility of NASA Earth observations (EOs) in visualizing heat vulnerability in Baltimore City and County. We mapped UHI and extreme heat using remote sensing observations, tailored a Heat Priority Score (HPS) to analyze heat data in tandem with socioeconomic data from American Community Survey (ACS), and utilized SOlar and LongWave Environmental Irradiance Geometry (SOLWEIG) model to estimate MTA riders’ thermal comfort. NASA EOs served to quantify the distribution and severity of extreme heat at the block group level, which was found to correlate with a lack of vegetation and presence of densely built-up surroundings. The HPS revealed that extreme heat events most negatively impact communities that exhibit social vulnerability through indicators such as age, race, and income. These findings will inform neighborhood resiliency plans for the MTA that can be incorporated into their Adaptation Resiliency Toolbox (ART).

Baltimore,↗
DOE OSTI2024

Overview of the Alaskan Layered Pollution and Chemical Analysis (ALPACA) Field Experiment

The Alaskan Layered Pollution And Chemical Analysis (ALPACA) field experiment was a collaborative study designed to improve understanding of pollution sources and chemical processes during winter (cold climate and low-photochemical activity), to investigate indoor pollution, and to study dispersion of pollution as affected by frequent temperature inversions. A number of the research goals were motivated by questions raised by residents of Fairbanks, Alaska, where the study was held. This paper describes the measurement strategies and the conditions encountered during the January and February 2022 field experiment, and reports early examples of how the measurements addressed research goals, particularly those of interest to the residents. Outdoor air measurements showed high concentrations of particulate matter and pollutant gases including volatile organic carbon species. During pollution events, low winds and extremely stable atmospheric conditions trapped pollution below 73 m, an extremely shallow vertical scale. Tethered-balloon-based measurements intercepted plumes aloft, which were associated with power plant point sources through transport modeling. Because cold climate residents spend much of their time indoors, the study included an indoor air quality component, where measurements were made inside and outside a house to study infiltration and indoor sources. In the absence of indoor activities such as cooking and/or heating with a pellet stove, indoor particulate matter concentrations were lower than outdoors; however, cooking and pellet stove burns often caused higher indoor particulate matter concentrations than outdoors. The mass-normalized particulate matter oxidative potential, a health-relevant property measured here by the reactivity with dithiothreiol, of indoor particles varied by source, with cooking particles having less oxidative potential per mass than pellet stove particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Atomization and vaporization characteristics of airblast fuel injection inside a venturi tube

This paper describes the experimental and numerical characterization of the capillary fuel injection, atomization, dispersion, and vaporization of liquid fuel in a coflowing air stream inside a single venturi tube. The experimental techniques used are all laser-based. Phase Doppler analyzer was used to characterize the atomization and vaporization process. Planar laser-induced fluorescence visualizations give good qualitative picture of the fuel droplet and vapor distribution. Limited quantitative capabilities of the technique are also demonstrated. A modified version of the KIVA-II was used to simulate the entire spray process, including breakup and vaporization. The advantage of venturi nozzle is demonstrated in terms of better atomization, more uniform F/A distribution, and less pressure drop. Multidimensional spray calculations can be used as a design tool only if care is taken for the proper breakup model, and wall impingement process.

Sun, H.↗
NASA NTRS2003

Control of Thermal Convection in Layered Fluids Using Magnetic fields

Immiscible fluid layers are found in a host of applications ranging from materials processing, for example the use of encapsulants in float zone crystal growth technique and a buffer layer in industrial Czochralski growth of crystals to prevent Marangoni convection, to heat transfer phenomena in day-to-day processes like the presence of air pockets in heat exchangers. In the microgravity and space processing realm, the exploration of other planets requires the development of enabling technologies in several fronts. The reduction in the gravity level poses unique challenges for fluid handling and heat transfer applications. The present work investigates the efficacy of controlling thermal convective flow using magnetic fluids and magnetic fields. The setup is a two-layer immiscible liquids system with one of the fluids being a diluted ferrofluid (super paramagnetic nano particles dispersed in carrier fluid). Using an external magnetic field one can essentially dial in a volumetric force - gravity level, on the magnetic fluid and thereby affect the system thermo-fluid behavior. The paper will describe the experimental and numerical modeling approach to the problem and discuss results obtained to date.

Ramachandran, N.↗
NASA NTRS1994

Two space scatterer formalism calculation of bulk parameters of thunderclouds

In a previous study, we used a modified two-space scatterer formalism of Twersky to establish for a cloud modeled as a statistically homogeneous distribution of spherical water droplets, the dispersion relations that determine its bulk propagation numbers and bulk indexes of refraction in terms of the vector equivalent scattering amplitude and the dyadic scattering amplitude of the single water droplet in isolation. The results were specialized to the forward direction of scattering while demanding that the scatterers preserve the incident polarization. We apply this approach to obtain specific numerical values for the macroscopic parameters of the cloud. We work with a cloud of density rho = 100 cm(exp -3), a wavelength lambda = 0.7774 microns, and with spherical water droplets of common radius alpha = 10 microns. In addition, the scattering medium is divided into three parts, the medium outside the cloud, moist air (the medium inside the cloud but outside the droplets), and the medium inside the spherical water droplets. The results of this report are applicable to a cloud of any geometry since the boundary does not interfere with the calculations. Also, it is important to notice the plane wave nature of the incidence wave in the moist atmosphere.

Phanord, Dieudonne D.↗
NASA NTRS1990

Oxidation characteristics of Ti-14Al-21Nb ingot alloy

Static oxidation kinetics of Ti14Al21Nb (wt pct) ingot alloy were studied in air over the temperature interval of 649 to 1093 C in a thermogravimetric apparatus. The oxidation products were characterized by x ray diffraction, electron microprobe analysis, energy dispersive x ray analysis, and Auger electron spectroscopy. Cross-sections of the oxidized samples were also examined using light and scanning electron microscopy. The oxidation rate was substantially lower than the conventional alloys of titanium, but the kinetics displayed a complex behavior involving two or more oxidation rates depending on the temperature and duration of exposure. The primary oxide formed was TiO2, but this oxide was doped with Nb. Small amounts of Al2O3 and TiN were also present in the scale. Diffusion of oxygen into the alloy was observed and the diffusivity seemed to be dependent on the microstructure of the metal. A model was presented to explain the oxidation behavior of the alloy in terms of the reduction in the oxygen diffusivity in the oxide caused by the modification of the defect structure of TiO2 by Nb ions.

Sankaran, Sankara N.↗
DOE OSTI1995

Nuclear Safety [Vol. 36, No. 1, January-June 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 1 The Chornobyl Accident Revisited, Part II: The State of the Nuclear Fuel Located Within the Chornobyl Sarcophagus, A A. Borovoi and A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 33 Nuclear Power Safety in Central and Eastern Europe, R. Wilson; 46 Safety of Nuclear Power Reactors in the Former Eastern European Countries, S. Chakraborty; 53 Technical Note: On the Definition of Common-Cause Failures, H. Paula; ACCIDENT ANALYSIS: 58 Modeling and Analysis of Core-Debris Recriticality During Hypothetical Severe Accidents in the Advanced Neutron Source Reactor, S.-H. Kim, V. Georgevich, D. B. Simpson, C. O. Slater, and R. P. Taleyarkhan; 68 Ignitability of Hydrogen/Oxygen/Diluent Mixtures in the Presence of Hot Surfaces, R. K. Kumar and G. W. Koroll; 94 Coupled RELAP5 and CONTAIN Accident Analysis Using PVM, K. A. Smith, A. J. Baratta, and G. E. Robinson; CONTROL AND INSTRUMENTATION: 109 Application of Fuzzy Logic in Nuclear Reactor Control Part I: An Assessment of State-of-the-Art, A. S. Heger, N. K. Alang-Rashid, and M. Jamshidi; DESIGN FEATURES: 122 Twenty-Third DOE/NRC Nuclear Air-Cleaning and Treatment Conference, R. R. Bellamy, J. J. Hayes, and M. W. First; ENVIRONMENTAL EFFECTS: 135 Atmospheric Dispersion and the Radiological Consequences of Normal Airborne Effluents from a Nuclear Power Plant, D. Fang, C. Z. Sun, and L. Yang; 142 Calculation of Distribution Coefficients for Radionuclides in Soils and Sediments, I. Puigdomenech and U. Bergstrom: OPERATING EXPERIENCES: 155 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 158 Operating Experience Feedback Report—Reliability of Safety-Related Steam Turbine-Driven Standby Pumps Used in U.S. Commercial Nuclear Power Plants, J. R. Boardman; 166 Turbine Building Hazards, H. L Ornstein; RECENT DEVELOPMENTS: 169 Reports, Standards, and Safety Guides, D. S. Queener; 175 Proposed Rule Changes as of Dec. 31,1994; ANNOUNCEMENTS: 32 Harvard School of Public Health In-Place Filter Testing Workshop; 134 International Conference on Advances in the Operational Safety of Nuclear Power Plants; 193 30th Tennessee Industries Week; 193 DOE Technical Standards Program 1995 Workshop; 194 Multiphase Flow Experiments and Instrumentation; 180 The Authors; 185 Indexes to Nuclear Safety, Volumes 34 and 35.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

Kinetic Model Development of the Carbon Vapor Deposition Reactor to Predict Performance Versus Design

To enable complete closure of a Sabatier-based atmosphere revitalization system where oxygen recovery is predicated by carbon dioxide reduction to methane and water, there must be a means of recovering the hydrogen from the byproduct methane. Honeywell’s carbon vapor deposition (CVD) reactor is a potential technology that provides a thermochemical means by which methane decomposes at high temperatures (> 1,000 °C) into hydrogen, solid carbon, and a small percentage of hydrocarbons. This advanced methane pyrolysis reactor utilizes a carbon fiber substrate to capture and immobilize the solid pyrolytic carbon rather than allowing continuation toward undesirable soot. The CVD reactor is a highly dynamic system whose characteristic geometry changes while operating. The reactor model was developed to provide predictive performance assessment capabilities at different operating conditions, geometries, and configurations to aid in technology development. Modeling of the CVD reactor considered the multitude of simultaneous and interdependent physicochemical processes that are entailed, including: (1) the homogeneous gas-phase pyrolysis reactions, (2) the heterogeneous deposition reactions, and (3) the transport of gaseous species in the reactor domains. A 1-D CVD reactor model was developed based on a general packed-bed, axial-dispersion plug flow reactor model with a solid substrate balance to account for carbon deposition. The model includes 59 species and 243 reactions constituting the methane pyrolysis reaction network, carbon deposition, and interconversions of oxygenate species. The reactor model was compared to brassboard test data, where it was shown to be able to capture the pertinent physicochemical phenomena as evidenced by its agreement with short-duration run data (i.e., looking at different temperature and flow rate conditions), as well as trend agreement with the brassboard reactor’s long-duration test data. Recommendations are provided on further model development steps as well as its use cases.

Air Revitalization↗
DOE OSTI2025

Development of the United States GReenhouse Gas and Air Pollutants Emissions System (GRA 2 PES)

In the U.S., emissions of greenhouse gases and air pollutants are often developed independently. Here, we describe the GReenhouse gas And Air Pollutants Emissions System (GRA 2 PES), which provides gridded emissions of fossil-fuel carbon dioxide (ffCO 2 ) and 93 air quality (AQ) species for 17 combustion and non-combustion sectors at 4 km × 4 km spatial resolution across the contiguous US. We find that the AQ emissions most spatially correlated with ffCO 2 are nitrogen oxides (NO x , ρ = 0.67), followed by sulfur dioxide (SO 2 , ρ = 0.51), carbon monoxide (CO, ρ = 0.44), and fine particulate matter (PM 2.5 , ρ = 0.38). We evaluate GRA 2 PES ffCO 2 emissions with an ensemble of publicly available regional and global inventories at national (Normalized Mean Bias (NMB) = +1.4%), state (NMB = +1.5%, R 2 = 0.98), and urban (NMB = +11.5%, R 2 = 0.97) scales. Nationally, the differences of publicly available inventories from the ensemble average range from −10.0% to +5.7%, and consistency diverges at state and urban scales. We simulate GRA 2 PES ffCO 2 in a particle dispersion model and compare to measurements of radiocarbon ( 14 C)-derived ffCO 2 collected in Los Angeles (August 2021), with results suggesting that GRA 2 PES ffCO 2 may be low by 19% for this city, but well within model-observation differences for other publicly available inventories (−43% to +94%). GRA 2 PES AQ/ffCO 2 ratios converted to concentration space generally agree with field observations (NMB = +4%, log R 2 = 0.90). Lastly, we present a method by which to utilize GRA 2 PES to derive AQ emission fluxes from ffCO 2 emissions.

Lyu, Congmeng [National Oceanic and Atmospheric Ad↗
NASA NTRS1974

Investigation of the oxidation behavior of dispersion stabilized alloys when exposed to a dynamic high temperature environment

The oxidation behavior of TD-NiCr and TD-NiCrAlY alloys have been studied at 2000 and 2200 F in static and high speed flowing air environments. The TD-NiCrAlY alloys preoxidized to produce an Al2O3 scale on the surface showed good oxidation resistance in both types of environments. The TD-NiCr alloy which had a Cr2O3 oxide scale after preoxidation was found to oxidize more than an order of magnitude faster under the dynamic test conditions than at comparable static test conditions. Although Cr2O3 normally provides good oxidation protection, it was rapidly lost due to formation of volatile CrO3 when exposed to the high speed air stream. The preferred oxide arrangement for the dynamic test consisted of an external layer of NiO with a porous mushroom type morphology, an intermediate duplex layer of NiO and Cr2O3, and a continuous inner layer of Cr2O3 in contact with the alloy substrate. An oxidation model has been developed to explain the observed microstructure and overall oxidation behavior of all alloys.

Tenney, D. R.↗
DOE OSTI2022

Modeling multiple scattering transient of an ultrashort laser pulse by spherical particles

The multiple scattering of an ultrashort laser pulse by a turbid dispersive medium (namely a cloud of bubbles in water) is investigated by means of Monte Carlo simulations. The theory of Gouesbet and Gréhan (2000) is used to derive an energetic model of the scattering transient. It is shown that the spreading and extinction of the pulse can be decoupled from the transient of scattering, which allows to describe each phenomenon individually. The transient of scattering is modeled with the Lorenz-Mie Theory and thus is also valid for a relative refractive index lower than one, contrary to the Debye series expansion which does not converge close to the critical angle. This is made possible after the introduction of a new physical object, the Scattering Impulse Response Function (SIRF) which allows to detect the different modes of scattering transient, in time and direction. Here the present approach is more generic, as it enables to simulate clouds of air bubbles in water, which was not possible previously. Two different approaches are proposed within the Monte Carlo framework. The first is a pure Monte Carlo approach where the delay due to the scattering is randomly drawn at each event, while the second is based on the transport of the whole scattering signal. They are both embedded in the Monte Carlo code Scatter3D. Both models produce equivalent trends and are validated against published numerical results. They are then applied to the multiple scattering of ultra short pulse by a cloud of bubble in water in the forward direction. The pulse spread due to the propagation in water is computed for a wide range of traveled distances and pulse durations, and the optimal pulse duration is given to minimize the pulse spread at a given distance. The main result is that the scattered photons exit the turbid medium earlier than the ballistic photons and produce a double peak related to the refraction in the bubble. This demonstrates the possibility to develop new diagnostics to characterize dynamic bubbly flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

This is a new project which started in May 1996. The main objective of the experimental/numerical study is to improve the understanding of the physics of two-way coupling between the dispersed phase and turbulence in a prototypical turbulent shear flow - homogeneous shear, laden with small liquid droplets (in gas) or gaseous bubbles (in liquid). The method of direct numerical simulation (DNS) is used to solve the full three-dimensional, time-dependent Navier-Stokes equations including the terms describing the two-way coupling between the dispersed phase and the carrier flow. The results include the temporal evolution of the three-dimensional energy and dissipation spectra and the rate of energy transfer across the energy spectrum to understand the fundamental physics of turbulence modulation, especially the effects of varying the magnitude of gravitational acceleration. The mean-square displacement and diffusivity of the droplets (or bubbles) of a given size and the preferential accumulation of droplets in low vorticity regions and bubbles in high vorticity regions will be examined in detail for different magnitudes of gravitational acceleration. These numerical results which will be compared with their corresponding measured data will provide a data base from which a subgrid-scale (SGS) model can be developed and validated for use in large-eddy simulation (LES) of particle-laden shear flows. Two parallel sets of experiments will be conducted: bubbles in an immiscible liquid and droplets in air. In both experiments homogeneous shear will be imposed on the turbulent carrier flow. The instantaneous velocities of the fluid and polydispersed-size particles (droplets or bubbles) will be measured simultaneously using a two-component Phase-Doppler Particle Analyzer (PDPA). Also, the velocity statistics and energy spectra for the carrier flow will be measured.

Elghobashi, Said↗
NASA NTRS1984

Structure of nonevaporating sprays - Measurements and predictions

Structure measurements were completed within the dilute portion of axisymmetric nonevaporating sprays (SMD of 30 and 87 microns) injected into a still air environment, including: mean and fluctuating gas velocities and Reynolds stress using laser-Doppler anemometry; mean liquid fluxes using isokinetic sampling; drop sizes using slide impaction; and drop sizes and velocities using multiflash photography. The new measurements were used to evaluate three representative models of sprays: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop interaction with turbulent fluctuations were ignored; and (3) a stochastic separated flow (SSF) model, where effects of both interphase slip and turbulent fluctuations were considered using random sampling for turbulence properties in conjunction with random-walk computations for drop motion. The LHF and DSF models were unsatisfactory for present test conditions-both underestimating flow widths and the rate of spread of drops. In contrast, the SSF model provided reasonably accurate predictions, including effects of enhanced spreading rates of sprays due to drop dispersion by turbulence, with all empirical parameters fixed from earlier work.

Solomon, A. S. P.↗
NASA NTRS1988

Wave/current interaction model

The wave-current interaction for the application to remote sensing data via numerical simulations and data comparison is modelled. Using the field data of surface current shear, wind condition and ambient wave spectrum, the numerical simulations of directional wave spectrum evolution were used to interpret and to compare with the aircraft data from Radar Ocean Wave Spectrometer (ROWS) and Surface Contour Radar (SCR) across the front during Frontal Air Sea Interaction Experiment (FASINEX). The wave-ice interaction was inspired by the observation of large amplitude waves hundreds of kms inside the ice pack in the Weddell Sea, resulting in breakup of the ice pack. The developed analysis of processes includes the refraction of waves at the pack edge, the effects of pack compression on wave propagation, wave train stability and buckling stability in the ice pack. Sources of pack compression and interaction between wave momentum and pack compression are investigated. Viscous camping of propagating waves in the marginal ice zone are also studied. The analysis suggests an explanation for the change in wave dispersion observed from the ship and the sequence of processes that cause ice pack breakup, pressure ridge formation and the formation of open bands of water.

Liu, A. K.↗
DOE OSTI2021

Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times

Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.

42 ENGINEERING↗
DOE OSTI2023

Quantifying spatiotemporal variability in occupant exposure to an indoor airborne contaminant with an uncertain source location

Well-mixed zone models are often employed to compute indoor air quality and occupant exposures. While effective, a potential downside to assuming instantaneous, perfect mixing is underpredicting exposures to high intermittent concentrations within a room. When such cases are of concern, more spatially resolved models, like computational-fluid dynamics methods, are used for some or all of the zones. But, these models have higher computational costs and require more input information. A preferred compromise would be to continue with a multi-zone modeling approach for all rooms, but with a better assessment of the spatial variability within a room. Here to do so, we present a quantitative method for estimating a room's spatiotemporal variability, based on influential room parameters. Our proposed method disaggregates variability into the variability in a room's average concentration, and the spatial variability within the room relative to that average. This enables a detailed assessment of how variability in particular room parameters impacts the uncertain occupant exposures. To demonstrate the utility of this method, we simulate contaminant dispersion for a variety of possible source locations. We compute breathing-zone exposure during the releasing (source is active) and decaying (source is removed) periods. Using CFD methods, we found after a 30 minutes release the average standard deviation in the spatial distribution of exposure was approximately 28% of the source average exposure, whereas variability in the different average exposures was lower, only 10% of the total average. We also find that although uncertainty in the source location leads to variability in the average magnitude of transient exposure, it does not have a particularly large influence on the spatial distribution during the decaying period, or on the average contaminant removal rate. By systematically characterizing a room's average concentration, its variability, and the spatial variability within the room important insights can be gained as to how much uncertainty is introduced into occupant exposure predictions by assuming a uniform in-room contaminant concentration. We discuss how the results of these characterizations can improve our understanding of the uncertainty in occupant exposures relative to well-mixed models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRSDate not supplied

Nonlinear Acoustic Propagation in Low-Density Media

Nonlinear acoustic propagation in the atmosphere is usually modeled using an augmented Burgers equation accounting for atmospheric absorption and weak nonlinearity. However, the weak-nonlinearity assumption may not apply to acoustic signals propagating from the lower atmosphere that are subsequently refracted downward from the upper atmosphere (e.g., stratosphere and thermosphere) due to the decreasing air density with increasing altitude [Lonzaga, et al., Geophysical Journal International, 200(3), pp.1347-1361]. Consequently, this paper discusses the effects of a strong nonlinearity that lead to an amplitude-dependent increase in signal propagation speed. This propagation speed is obtained using a perturbation expansion where the leading-order term is simply the small-amplitude sound speed while the first-order term is the existing expression that gives rise to waveform steepening and stretching. Furthermore, the second-order term is proportional to the spectral density of the acoustic signal and is inversely proportional to the local density of the medium. Consequently, for an impulsive signal such as a sonic boom or an infrasound, the second-order effect causes a dispersion of the signal similar to the observed dispersion of acoustic signals from supersonic Concorde as well as from large explosions. This paper also discusses the extension of these results in general low-density propagation media.

sonic booms↗
NASA NTRS2011

In Situ Observations and Sampling of Volcanic Emissions with Unmanned Aircraft: A NASA/UCR Case Study at Turrialba Volcano, Costa Rica

Burgeoning new technology in the design and development of robotic aircraft-unmanned aerial vehicles (UAVs)-presents unprecedented opportunities for the volcanology community to observe, measure, and sample eruption plumes and drifting volcanic clouds in situ. While manned aircraft can sample dilute parts of such emissions, demonstrated hazards to air breathing, and most particularly turbine, engines preclude penetration of the zones of highest ash concentrations. Such areas within plumes are often of highest interest with respect to boundary conditions of applicable mass-loading retrieval models, as well as Lagrangian, Eulerian, and hybrid transport models used by hazard responders to predict plume trajectories, particularly in the context of airborne hazards. Before the 2010 Ejyafyallajokull eruption in Iceland, ICAO zero-ash-tolerance rules were typically followed, particularly for relatively uncrowded Pacific Rim airspace, and over North and South America, where often diversion of aircraft around ash plumes and clouds was practical. The 2010 eruption in Iceland radically changed the paradigm, in that critical airspace over continental Europe and the United Kingdom were summarily shut by local civil aviation authorities and EURO CONTROL. A strong desire emerged for better real-time knowledge of ash cloud characteristics, particularly ash concentrations, and especially for validation of orbital multispectral imaging. UAV platforms appear to provide a viable adjunct, if not a primary source, of such in situ data for volcanic plumes and drifting volcanic clouds from explosive eruptions, with prompt and comprehensive application to aviation safety and to the basic science of volcanology. Current work is underway in Costa Rica at Turrialba volcano by the authors, with the goal of developing and testing new small, economical UAV platforms, with miniaturized instrument payloads, within a volcanic plume. We are underway with bi-monthly deployments of tethered SO2-sondes and are in the planning stages for the deployment of the SIERRA UAV to our site in March 2013. We will be conducting in situ observations simultaneously with ASTER orbital multispectral TIR data acquisitions, in order to compare in situ measurements with estimates of SO2 mass loading and dispersion derived from ASTER data. Though small UAVs are now being considered for use in active volcanic areas for in situ sampling of emissions (e.g., efforts by our group, and by our colleagues at the INGV in Italy and the Applied Science University in Dusseldorf, Germany, and others in the United Kingdom and Iceland), and also for remote sensing, much more needs to be done in the way of instrument development, and in developing small UAVs for both low altitude (tropospheric) and high altitude (stratospheric) applications. In particular, the development of all weather and day/night operational flight capabilities in close proximity to hazardous topography is crucial to a truly responsive volcano in situ measurement system. Finally, it is imperative that national civil aviation authorities recognize the unique benefits of such platforms. It is important that authorities understand that severely restricting or not deploying such tools in airspace over restless volcanoes or within eruption plumes, ostensibly because of the perceived (small) risk that such unmanned aircraft pose to manned air operations, itself poses a bigger transcendental risk to proximal populations and particularly to the aviation community, itself.

Volcanic Ash Advisory Centers (VAACs)↗