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 109 records · Page 6

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.↗

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↗

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.↗

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.↗

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↗

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)↗

Commercial Non-Dispersive Infrared Spectroscopy Sensors for Sub-Ambient Carbon Dioxide Detection

Carbon dioxide produced through respiration can accumulate rapidly within closed spaces. If not managed, a crew's respiratory rate increases, headaches and hyperventilation occur, vision and hearing are affected, and cognitive abilities decrease. Consequently, development continues on a number of CO2 removal technologies for human spacecraft and spacesuits. Terrestrially, technology development requires precise performance characterization to qualify promising air revitalization equipment. On-orbit, instrumentation is required to identify and eliminate unsafe conditions. This necessitates accurate in situ CO2 detection. Recursive compensation algorithms were developed for sub-ambient detection of CO2 with commercial off-the-shelf (COTS) non-dispersive infrared (NDIR) sensors. In addition, the source of the exponential loss in accuracy is developed theoretically. The basis of the loss can be explained through thermal, Doppler, and Lorentz broadening effects that arise as a result of the temperature, pressure, and composition of the gas mixture under analysis. The objective was to develop a mathematical routine to compensate COTS CO2 sensors relying on NDIR over pressures, temperatures, and compositions far from calibration conditions. The routine relies on a power-law relationship for the pressure dependency of the sensors along with an equivalent pressure to account for the composition dependency. A Newton-Raphson iterative technique solves for actual carbon dioxide concentration based on the reported concentration. Moreover, first principles routines were established to predict mixed-gas spectra based on sensor specifications (e.g., optical path length). The first principles model can be used to parametrically optimize sensors or sensor arrays across a wide variety of pressures/temperatures/ compositions. In this work, heuristic scaling arguments were utilized to develop reasonable compensation techniques. Experimental results confirmed this approach and provided evidence that composition broadening significantly alters spectra when pressure is reduced. Consequently, a recursive compensation technique was developed with the Newton-Raphson method, which was subsequently verified through experimentation.

Swickrath, Michael J.↗

Absolute Position of Targets Measured Through a Chamber Window Using Lidar Metrology Systems

Lidar is a useful tool for taking metrology measurements without the need for physical contact with the parts under test. Lidar instruments are aimed at a target using azimuth and elevation stages, then focus a beam of coherent, frequency modulated laser energy onto the target, such as the surface of a mechanical structure. Energy from the reflected beam is mixed with an optical reference signal that travels in a fiber path internal to the instrument, and the range to the target is calculated based on the difference in the frequency of the returned and reference signals. In cases when the parts are in extreme environments, additional steps need to be taken to separate the operator and lidar from that environment. A model has been developed that accurately reduces the lidar data to an absolute position and accounts for the three media in the testbed air, fused silica, and vacuum but the approach can be adapted for any environment or material. The accuracy of laser metrology measurements depends upon knowing the parameters of the media through which the measurement beam travels. Under normal conditions, this means knowledge of the temperature, pressure, and humidity of the air in the measurement volume. In the past, chamber windows have been used to separate the measuring device from the extreme environment within the chamber and still permit optical measurement, but, so far, only relative changes have been diagnosed. The ability to make accurate measurements through a window presents a challenge as there are a number of factors to consider. In the case of the lidar, the window will increase the time-of-flight of the laser beam causing a ranging error, and refract the direction of the beam causing angular positioning errors. In addition, differences in pressure, temperature, and humidity on each side of the window will cause slight atmospheric index changes and induce deformation and a refractive index gradient within the window. Also, since the window is a dispersive media, the effect of both phase and group indices have to be considered. Taking all these factors into account, a method was developed to measure targets through multiple regions of different materials and produce results that are absolute measurements of target position in three-dimensional space, rather than simply relative position. The environment in which the lidar measurements are taken must be broken down into separate regions of interest and each region solved for separately. In this case, there were three regions of interest: air, fused silica, and vacuum. The angular position of the target inside the chamber is solved using only phase index and phase velocity, while the ranging effects due to travel from air to glass to vacuum/air are solved with group index and group velocity. When all parameters are solved simultaneously, an absolute knowledge of the position of each target within an environmental chamber can be derived. Novel features of this innovation include measuring absolute position of targets through multiple dispersive and non-dispersive media, deconstruction of lidar raw data from a commercial off-the-shelf unit into reworkable parameters, and use of group velocities to reduce range data. Measurement of structures within a vacuum chamber or other harsh environment, such as a furnace, may now be measured as easily as if they were in an ambient laboratory. This analysis permits transformation of the raw data into absolute spatial units (e.g., mm). This technique has also been extended to laser tracker, theodolite, and cathetometer measurements through refractive media.

Kubalak, David↗

AG Channel Measurement and Modeling Results for Over-Sea Conditions

This report describes results from flight tests conducted in an over-sea environment, for the purpose of characterizing the air-to-ground (AG) channel, for future unmanned aircraft system (UAS) communication system analysis and design. These results are for the first of a set of several flight tests conducted in different ground site (GS) environments. An ultimate aim of all these tests is the development of models for the AG channel that can be used in communication system evaluation. In this report we provide measured results for propagation path loss, root-mean square delay spread (RMS-DS), and the correlation coefficient of the primary received signal components on the four antennas (two antennas for C-band, two for L-band). For path loss, the curved-earth two-ray model provides a reasonable fit to the measured data, altered by several dB at the shortest link distances by aircraft antenna pattern effects. This two-ray model also accounts for the majority of measured RMS-DS results of a few tens of nanoseconds, except for the occasional intermittent reflections from surface objects. These intermittent reflections yield RMS-DS values up to several hundred nanoseconds. For portions of the flight path that were over a harbor area highly populated with boats, the channel was found to be more "continuously dispersive," with RMS-DS reaching approximately 250 ns. A separate model will be developed for this over-harbor setting. The correlation coefficient results are still undergoing analysis; preliminary observations are that correlation between signals on the same-band antennas is generally large (>0.6) for the C-band straight flight paths, whereas for the L-band signals and for the oval-shaped flight paths the correlation is generally small (below 0.4). Inter-band correlations are typically very small, and are well modeled as zero-mean Gaussian in distribution, with a standard deviation less than 0.2. Hence the over-sea channel effects in the two bands can be considered uncorrelated, which will allow for good diversity gains in dual-band systems. We describe initial modeling approaches for the over-sea channel; complete models for this and the over-harbor setting will appear in a subsequent report.

Aircraft Communications↗

The INSTEP Monitoring Network: Merging High-and-Low Cost Measurements to Characterize California Wildfires

Despite challenges with data quality and scope, low-cost sensor networks have skyrocketed in popularity over the last 15 years, making air quality data available on refined spatial scales. More recently, studies have leveraged both high and low-quality instruments to create stronger “hybrid” models, with most studies focusing on particulate matter. Low-cost measurements typically represent ground-level emissions only, providing context for human health issues from climate change-driven events such as wildfires. Since low-cost sensors’ capabilities are localized, daily events and microclimates tend to dominate the data rather than larger regional or atmospheric trends. Likewise, their low cost explains their high uncertainty. In contrast, some regulatory-grade instruments produce column measurements as well, providing reliable information on a broader scope. To bridge this gap while expanding into gas-phase measurements, we deployed 12 air quality sensor packages in California, USA during the 2022 wildfire season. These INSTEP (Inexpensive Network Sensor Technology Exploring Pollution) monitors measure carbon monoxide (CO), carbon dioxide (CO2), ozone (O3), nitrogen dioxide (NO2), and several hydrocarbons including methane (CH4) and formaldehyde (HCHO). Half of the monitors were co-located with remote sensing spectrometers: NASA Pandora and Total Column Carbon Observing Network (TCCON). The overlap in pollutants includes NO2, O3, and HCHO between the INSTEP monitors and the Pandora column measurements. TCCON covers column CO, CO2, and CH4, rounding out our comparison. Most of the monitors were distributed throughout the San Francisco Bay area, and an additional three were located within 100 km of Los Angeles. The sites ranged in geographic and population characteristics, including desert, mountainous, coastal, and urban locations. Since varying environmental conditions such as temperature and pressure are known to challenge sensor performance, we will apply newer sensor “calibration” techniques meant to combat this. We will normalize our sensor signals by z-scoring them prior to applying a single calibration model in the form of multivariate linear regression or an artificial neural network. While this technique has been validated for the hydrocarbon and ozone sensor types (metal oxide), it has not yet been tested on electrochemical and non-dispersive infrared sensors, which are also used in the INSTEP monitors. This will serve as a test to see if this normalization technique – or another – is most effective in accounting for environmental differences among sensors. Related data analysis efforts have found success with a variety of geospatial analysis techniques, including weighted network models in which high-quality instruments are given higher weights than their low-cost counterparts. Our preliminary analysis will focus on kriging, which uses a Gaussian algorithm to assign weights, providing estimated pollution levels at locations between monitors. Smoke trajectory and evolution will also be considered using both measurement types. We also aim to baseline subtract our emission estimates from each region to determine which portion of emissions are regional and local, further characterizing burn differences in northern and southern California fires. Future directions include using INSTEP jointly with TEMPO satellite data, and mobile deployments on aircraft and uncrewed aerial vehicles (UAV).

Low-cost sensors↗

An Overview of Experiments and Modeling of Polysiloxane-Coated Thermal Protection Systems for Missions to Mars, Titan, and Beyond.

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the historic Stardust mission, where it successfully returned samples from a comet’s tail and has since been instrumental in delivering payloads to the surface of Mars [1-3]. Most recently, PICA enabled the safe return of samples collected from asteroid Bennu as part of the OSIRIS-REx mission. This rich legacy underscores PICA’s critical role in allowing NASA’s most ambitious exploration missions. However, the friable nature of its phenolic phase presents challenges during handling and pre-launch activities. To mitigate this issue, PICA is coated with a polysiloxane resin system, which serves to suppress particulate dispersion and thereby safeguard spacecraft components. A comprehensive understanding of the polysiloxane resin’s behavior is imperative, as it profoundly shapes the material response of PICA during atmospheric entry by influencing its thermal and oxidative stability. This influence extends to thermocouple plugs embedded within thermal protection systems. These plugs have demonstrated their significance in missions such as Mars Science Laboratory (MSL) and Mars 2020, where the MEDLI and MEDLI2 instrumentation suites delivered in-valuable insights into the performance of thermal protection systems during entry into the Martian atmosphere [4]. Looking ahead, missions such as Dragonfly, set to descend into Titan’s dense atmosphere, aim to leverage advanced sensor technologies to further refine our understanding of thermal protection response [5]. Moreover, thermocouple plugs play an essential role in validating cutting-edge material response models, such as those pioneered under NASA’s Entry Systems Modeling Project (ESM), designed, in-part, to predict the operational integrity of thermal protection systems under the extreme stresses of atmospheric entry. To achieve these modeling goals, ground-based experiments are crucial to provide the foundational data necessary for developing and refining these predictive tools. To this end, an extensive test campaign was conducted at the Hypersonic Materials Environmental Test System (HyMETS) to investigate the high-temperature behavior of the polysiloxane resin in an air environment [6]. These experiments revealed critical phenomena, including the formation of a silicon oxycarbide layer that enhances oxidation resistance, moderates surface temperatures, and alters in-depth thermal response. Building on these findings, subsequent tests were designed to simulate atmospheric entry conditions in reactive gases, such as CO2 and N2, to mimic the environments of Mars and Titan, respectively, as well as non-reactive gases representing the atmospheres of the Ice Giants (Neptune and Uranus). A heating rate dependent decomposition mechanism has been identified for the polysiloxane resin under oxidizing conditions (Fig. 1). In the initial stage, the resin and the underlying thermal protection system undergo pyrolysis, rapidly generating a thin amorphous silicon oxycarbide interwoven with carbonaceous char and residual fibers from PICA. During the second stage, the nascent oxide layer establishes a robust, oxidation-resistant thermal barrier coating, which significantly impedes heat transfer to the underlying carbonaceous char, resulting in a stagnation of the surface temperature. A key factor contributing to this thermal resistance is the low recombination efficiency of atomic oxygen (γ), which further diminishes the heat load on the material’s interior layers [7]. Moreover, as the surface temperature stagnates, the silicon oxycarbide phase separates into distinct regions of silica and free graphite. Ultimately, when the heat flux reaches a critical threshold, a third stage is triggered, leading to the breakdown of the coating through carbothermal reduction, exposing the underlying char layer. This exposure leads to a dramatic surface temperature spike, driven by highly exothermic reactions between atomic oxygen and the char layer, further accelerating material degradation. A detailed mass and heat transfer model of PICA coated with polysiloxane resin was implemented in the Porous material Analysis Toolbox based on OpenFOAM, PATO [8]. The initial stage was considered negligible in this model because the resin decomposition occurs rapidly within a thin surface layer. Instead, the coating was directly considered as an oxygen-resistant thermal barrier coating. For the second stage, the thin amorphous silicon oxycarbide was treated as a pure silica surface to simplify the thermochemical behavior. The model ac-counts for surface equilibrium processes using representative elements of the coating-environment system. For the third stage, specific boundary conditions were developed to estimate the onset and progression of the coating removal. Two-dimensional material response simulations were conducted to compare uncoated and coated PICA using boundary conditions calibrated with HyMETS data. Fig. 2 illustrates that the simulations closely align with experimental data, successfully reproducing measured temperature profiles. This work will include the latest advancements in the coating model, including the calibration of recombination of atomic oxygen at the surface during the second phase. These simulated results will be further validated against additional CO2 data points from HyMETS, reinforcing the models’ predictive capabilities. These mechanisms and their effects on thermal protection systems, including thermochemical behavior and thermocouple probe performance in extreme environments, provide crucial insights for optimizing spacecraft designs that safeguard scientific payload and ensure mission success in future planetary exploration endeavors.

Active Oxidation↗

Cyclic Oxidation of High-Temperature Alloy Wires in Air

High-temperature alloy wires are proposed for use in seal applications for future re-useable space vehicles. These alloys offer the potential for improved wear resistance of the seals. The wires must withstand the high temperature environments the seals are subjected to as well as maintain their oxidation resistance during the heating and cooling cycles of vehicle re-entry. To model this, the wires were subjected to cyclic oxidation in stagnant air. of this layer formation is dependent on temperature. Slow growing oxides such as chromia and alumina are desirable. Once the oxide is formed it can prevent the metal from further reacting with its environment. Cyclic oxidation models the changes in temperature these wires will undergo in application. Cycling the temperature introduces thermal stresses which can cause the oxide layer to break off. Re-growth of the oxide layer consumes more metal and therefore reduces the properties and durability of the material. were used for cyclic oxidation testing. The baseline material, Haynes 188, has a Co base and is a chromia former while the other two alloys, Kanthal A1 and PM2000, both have a Fe base and are alumina formers. Haynes 188 and Kanthal A1 wires are 250 pm in diameter and PM2000 wires are 150 pm in diameter. The coiled wire has a total surface area of 3 to 5 sq cm. The wires were oxidized for 11 cycles at 1204 C, each cycle containing a 1 hour heating time and a minimum 20 minute cooling time. Weights were taken between cycles. After 11 cycles, one wire of each composition was removed for analysis. The other wire continued testing for 70 cycles. Post-test analysis includes X-Ray Diffraction (XRD), Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for phase identification and morphology.

Reigel, Marissa M.↗

Characterizing the wake vortex signature for an active line of sight remote sensor

A recurring phenomenon, described as a wake vortex, develops as an aircraft approaches the runway to land. As the aircraft moves along the runway, each of the wing tips generates a spiraling and expanding cone of air. During the lifetime of this turbulent event, conditions exist over the runway which can be hazardous to following aircraft, particularly when a small aircraft is following a large aircraft. Left to themselves, these twin vortex patterns will converge toward each other near the center of the runway, harmlessly dissipating through interaction with each other or by contact with the ground. Unfortunately, the time necessary to disperse the vortex is often not predictable, and at busy airports can severely impact terminal area productivity. Rudimentary methods of avoidance are in place. Generally, time delays between landing aircraft are based on what is required to protect a small aircraft. Existing ambient wind conditions can complicate the situation. Reliable detection and tracking of a wake vortex hazard is a major technical problem which can significantly impact runway productivity. Landing minimums could be determined on the basis of the actual hazard rather than imposed on the basis of a worst case scenario. This work focuses on using a windfield description of a wake vortex to generate line-of-sight Doppler velocity truth data appropriate to an arbitrarily located active sensor such as a high resolution radar or lidar. The goal is to isolate a range Doppler signature of the vortex phenomenon that can be used to improve detection. Results are presented based on use of a simplified model of a wake vortex pattern. However, it is important to note that the method of analysis can easily be applied to any vortex model used to generate a windfield snapshot. Results involving several scan strategies are shown for a point sensor with a range resolution of 1 to 4 meters. Vortex signatures presented appear to offer potential for detection and tracking.

Heil, Robert Milton↗

An exact solution of an augmented Burgers equation and amplitude-dependent acoustic propagation speed

Nonlinear sound propagation in the atmosphere is usually modeled using an augmented Burgers equation accounting for a weak nonlinearity and atmospheric absorption. Because the absorption includes the molecular vibrational relaxation, such a Burgers equation is more complex than the regular Burgers equation that only accounts for the thermoviscous dissipation in the absorption. Although an exact solution of the regular Burgers equation has long been derived using the Cole- Hopf transform, an exact solution of the augmented Burgers equation has not been derived previously. Thus, this paper presents an exact solution of the augmented Burgers equation. This novel solution is shown to be equivalent to the solution using the Cole-Hopf transform when the absorption only involves thermoviscous dissipation. It can also be reduced to the known solution of an N-wave when the absorption is ignored. The augmented Burgers equation is an approximation valid for weak nonlinearity. However, this assumption may not be accurate for acoustic signals propagating from the lower atmosphere and which 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, the current paper also discusses the effects of a strong nonlinearity that lead to an amplitude-dependent increase in signal propagation speed. For an impulsive signal such as a sonic boom, these effects cause a dispersion of the signal similar to the observed dispersion of acoustic signals from supersonic Concorde as well as from large explosions.

Joel B Lonzaga↗

Modeling the Frozen-In Anticyclone in the 2005 Arctic Summer Stratosphere

Immediately following the breakup of the 2005 Arctic spring stratospheric vortex, a tropical air mass, characterized by low potential vorticity (PV) and high nitrous oxide (N2O), was advected poleward and became trapped in the easterly summer polar vortex. This feature, known as a "Frozen-In Anticyclone (FrIAC)", was observed in Earth Observing System (EOS) Aura Microwave Limb Sounder (MLS) data to span the potential temperature range from approximately 580 to 1100 K (approximately 25 to 40 km altitude) and to persist from late March to late August 2005. This study compares MLS N2O observations with simulations from the Global Modeling Initiative (GMI) chemistry and transport model, the GEOS-5/MERRA Replay model, and the VanLeer Icosahedral Triangular Advection isentropic transport model to elucidate the processes involved in the lifecycle of the FrIAC which is here divided into three distinct phases. During the "spin-up phase" (March to early April), strong poleward flow resulted in a tight isolated anticyclonic vortex at approximately 70-90 deg N, marked with elevated N2O. GMI, Replay, and VITA all reliably simulted the spin-up of the FrIAC, although the GMI and Replay peak N2O values were too low. The FrIAC became trapped in the developing summer easterly flow and circulated around the polar region during the "anticyclonic phase" (early April to the end of May). During this phase, the FrIAC crossed directly over the pole between the 7th and 14th of April. The VITA and Replay simulations transported the N2O anomaly intact during this crossing, in agreement with MLS, but unrealistic dispersion of the anomaly occurred in the GMI simulation due to excessive numerical mixing of the polar cap. The vortex associated with the FrIAC was apparently resistant to the weak vertical hear during the anticyclonic phase, and it thereby protected the embedded N20 anomaly from stretching. The vortex decayed in late May due to diabatic processes, leaving the N2O anomaly exposed to horizontal and vertical wind shears during the "shearing phase" (June to August). The observed lifetime of the FrIAC during this phase is consistent with time-scales calculated from the ambient horizontal and vertical wind shear. Replay maintained the horizontal structure of the N2O anomaly similar to NILS well into August. The VITA simulation also captured the horizontal structure of the FrIAC during this phase, but VITA eventually developed fine-scale N2O structure not observed in MLS data.

Allen, D. R.↗

Turbulent dispersion of the icing cloud from spray nozzles used in icing tunnels

To correctly simulate flight in natural icing conditions, the turbulence in an icing simulator must be as low as possible. But some turbulence is required to mix the droplets from the spray nozzles and achieve an icing cloud of uniform liquid water content. The goal for any spray system is to obtain the widest possible spray cloud with the lowest possible turbulence in the test section of a icing tunnel. This investigation reports the measurement of turbulence and the three-dimensional spread of the cloud from a single spray nozzle. The task was to determine how the air turbulence and cloud width are affected by spray bars of quite different drag coefficients, by changes in the turbulence upstream of the spray, the droplet size, and the atomizing air. An ice accretion grid, located 6.3 m downstream of the single spray nozzle, was used to measure cloud spread. Both the spray bar and the grid were located in the constant velocity test section. Three spray bar shapes were tested: the short blunt spray bar used in the NASA Lewis Icing Research Tunnel, a thin 14.6 cm chord airfoil, and a 53 cm chord NACA 0012 airfoil. At the low airspeed (56 km/hr) the ice accretion pattern was axisymmetric and was not affected by the shape of the spray bar. At the high airspeed (169 km/hr) the spread was 30 percent smaller than at the low airspeed. For the widest cloud the spray bars should be located as far upstream in the low velocity plenum of the icing tunnel. Good comparison is obtained between the cloud spread data and predicitons from a two-dimensional cloud mixing computer code using the two equation turbulence (k epsilon g) model.

Marek, C. J.↗

Particle-laden weakly swirling free jets: Measurements and predictions

A theoretical and experimental investigation of particle-laden, weakly swirling, turbulent free jets was conducted. Glass particles, having a Sauter mean diameter of 39 microns, with a standard deviation of 15 microns, were used. A single loading ratio (the mass flow rate of particles per unit mass flow rate of air) of 0.2 was used in the experiments. Measurements are reported for three swirl numbers, ranging from 0 to 0.33. The measurements included mean and fluctuating velocities of both phases, and particle mass flux distributions. Measurements were also completed for single-phase non-swirling and swirling jets, as baselines. Measurements were compared with predictions from three types of multiphase flow analysis, as follows: (1) locally homogeneous flow (LHF) where slip between the phases was neglected; (2) deterministic separated flow (DSF), where slip was considered but effects of turbulence/particle interactions were neglected; and (3) stochastic separated flow (SSF), where effects of both interphase slip and turbulence/particle interactions were considered using random sampling for turbulence properties in conjunction with random-walk computations for particle motion. Single-phase weakly swirling jets were considered first. Predictions using a standard k-epsilon turbulence model, as well as two versions modified to account for effects of streamline curvature, were compared with measurements. Predictions using a streamline curvature modification based on the flux Richardson number gave better agreement with measurements for the single-phase swirling jets than the standard k-epsilon model. For the particle-laden jets, the LHF and DSF models did not provide very satisfactory predictions. The LHF model generally overestimated the rate of decay of particle mean axial and angular velocities with streamwise distance, and predicted particle mass fluxes also showed poor agreement with measurements, due to the assumption of no-slip between phases. The DSF model also performed quite poorly for predictions of particle mass flux because turbulent dispersion of the particles was neglected. The SSF model, which accounts for both particle inertia and turbulent dispersion of the particles, yielded reasonably good predictions throughout the flow field for the particle-laden jets.

Bulzan, Daniel L.↗

The Use of Satellite-Measured Aerosol Optical Depth to Constrain Biomass Burning Emissions Source Strength in a Global Model GOCART

Small particles in the atmosphere, called "atmospheric aerosol" have a direct effect on Earth climate through scattering and absorbing sunlight, and also an indirect effect by changing the properties of clouds, as they interact with solar radiation as well. Aerosol typically stays in the atmosphere for several days, and can be transported long distances, affecting air quality, visibility, and human health not only near the source, but also far downwind. Smoke from vegetation fires is one of the main sources of atmospheric aerosol; other sources include anthropogenic pollution, dust, and sea salt. Chemistry transport models (CTMs) are among the major tools for studying the atmospheric and climate effects of aerosol. Due to the considerable variation of aerosol concentrations and properties on many temporal and spatial scales, and the complexity of the processes involved, the uncertainties in aerosol effects on climate are large, as is featured in the latest report of Intergovernmental Panel on Climate Change (IPCC) in 2007. Reducing this uncertainty in the models is very important both for predicting future climate scenarios and for regional air quality forecasting and mitigation. During vegetation fires, also called biomass burning (BB) events, complex mixture of gases and particles is emitted. The amount of BB emissions is usually estimated taking into account the intensity and size of the fire and the properties of burning vegetation. These estimates are input into CTMs to simulate BB aerosol. Unfortunately, due to large variability of fire and vegetation properties, the quantity of BB emissions is very difficult to estimate and BB emission inventories provide numbers that can differ by up to the order of magnitude in some regions. Larger uncertainties in data input make uncertainties in model output larger as well. A powerful way to narrow the range of possible model estimates is to compare model output to observations. We use satellite observations of aerosol properties, specifically aerosol optical depth, which is directly proportional to the amount of aerosol in the atmosphere, and compare it to the model output. Assuming the model represents aerosol transport and particle properties correctly, the amount of BB emissions determines the simulated aerosol optical depth. In this study, we explore the regional performance of 13 commonly used emission estimates. These are each input to global Goddard Chemistry Aerosol Radiation and Transport (GOCART) model. We then evaluate how well each emission estimate reproduces the smoke aerosol optical depth measured by the MODIS instrument. We compared GOCART-simulate aerosol optical depth with that measured from the satellite for 124 fire cases around the world during 2006 and 2007. We summarize the regional performance of each emission inventory and discuss reasons for their differences by considering the assumptions made during their development. We also show that because stronger wind disperses smoke plumes more readily, in cases with stronger wind, a larger increase in emission amount is needed to increase aerosol optical depth. In quiet, low-wind-speed environments, BB emissions produce a more significant increase in aerosol optical depth, other things being equal. Using the region-specific, quantitative relationships derived in our paper, together with the wind speed obtained from another source for a given fire case, we can constrain the amount of emission required in the model to reproduce the observations. The results of this paper are useful to the developers of BB emission inventories, as they show the strengths and weaknesses of individual emission inventories in different regions of the globe, and also for modelers who use these inventories and wish to improve their model results.

Petrenko, Mariya↗