Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

142 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2004

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

Spectral Data Fusion From Handheld Laser-Induced Breakdown Spectroscopy (LIBS) and X-ray Fluorescence (XRF) Analyzers for Improved Detection of Cerium in a Simulated Dispersal Accident

Here, this work implements a mid-level data fusion methodology on spectral data from handheld X-ray fluorescence and laser-induced breakdown spectroscopy analyzers to quantify plutonium surrogate (CeO 2 ) contamination in soil samples for the first time. Spectral data from each analyzer were used independently to train supervised machine learning regressions to predict Ce concentration. Fused features from both data sets were then used to train the same models, comparing prediction performance by evaluating model precision and sensitivity. Fusing principal component scores from the two sensors yielded an order of magnitude improvement in precision and sensitivity of predictions made with an artificial neural network, compared to predictions made by models trained on independent sensor data. As a result, a boosted ensemble trained on the fused spectral features yielded an ideal predictor with root-mean-squared error on the order of 10 –6 and calculated limit of detection order 10 –5 wt %.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1994

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2010

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

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

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

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↗
NASA NTRS1998

Application of a Planar Doppler Velocimetry System to a High Reynolds Number Compressible Jet

A Planar Doppler Velocimetry (PDV) system has been constructed and used to investigate the instantaneous turbulent velocity structure of a round high-speed compressible air jet with a low-speed co-flow. The exit condition was Mach=0.85 at ambient pressure, yielding a Reynolds number of about 650,000 on diameter. The PDV system was installed at NASA Langley Research Center in the Small Anechoic Jet Facility (SAJF), a chamber in which both the acoustic and aerodynamic properties of jets can be studied. For this test, the goal was to gather data which can be used to relate the turbulence structure of the jet to the levels and character of the acoustic noise produced by the jet. The current PDV system can acquire single-velocity-component, single-shot, planar images (15ns exposures) at 30 Hz. For this paper, the primary data set consists of 240 frames of velocity data acquired with both the jet and the low-speed co-flow seeded with light-scattering articles. Thus, velocities could be measured everywhere in the jet shear layer, both in the jet fluid and in the entrained co-flow. Some data were also taken with only the jet flow seeded. These provided mixing concentration images along with the reduced velocity fields. Other images were taken with only the co-flow seeded. These produced unique quantitative images of high speed entrainment. Optical "laser speckle" noise is the largest source of random noise in pulsed PDV systems. Components for the PDV imaging system were specifically selected to minimize speckle noise. To reduce systematic velocity errors due to laser drift, a frequency monitoring reference leg with a temperature-tuned reference iodine cell, was employed. In the course of this study, a novel flow seeder was developed. It enabled continuously variable seeding of the flow with particles of Sheared Pyrogenic Amorphous Hydrophobic Silica (SPAHS). The seeder comprised a dry fluidized bed hopper and a supersonic nozzle "pickup." Shearing action in the pickup dispersed the seed material in an exceptionally fine cloud (approximately 0.3 micron). These particles followed the flow well, did not clump or cake on screens or model surfaces, and were not susceptible to evaporation. Because of the refractory nature of the particles, SPAHS seeding should also be applicable to anticipated future testing at high temperatures.

Smith, Michael W.↗
NASA NTRS1999

Bubble Generation in a Flowing Liquid Medium and Resulting Two-Phase Flow in Microgravity

The present investigation reports a study of bubble generation under reduced gravity conditions, using both a co-flow and a cross-flow configuration. This study may be used in the conceptual design of a space-based thermal management system. Ensuing two-phase flow void fraction can be accurately monitored using a single nozzle gas injection system within a continuous liquid flow conduit, as utilized in the present investigation. Accurate monitoring of void fraction leads to precise control of heat and mass transfer coefficients related to a thermal management system; hence providing an efficient and highly effective means of removing heat aboard spacecraft or space stations. Our experiments are performed in parabolic flight aboard the modified DC-9 Reduced Gravity Research Aircraft at NASA Lewis Research Center, using an air-water system. For the purpose of bubble dispersion in a flowing liquid, we use both a co-flow and a cross-flow configuration. In the co-flow geometry, air is introduced through a nozzle in the same direction with the liquid flow. On the other hand, in the cross-flow configuration, air is injected perpendicular to the direction of water flow, via a nozzle protruding inside the two-phase flow conduit. Three different flow conduit (pipe) diameters are used, namely, 1.27 cm, 1.9 cm and 2.54 cm. Two different ratios of nozzle to pipe diameter (D(sub N))sup * are considered, namely (D(sub N))sup * = 0.1 and 0.2, while superficial liquid velocities are varied from 8 to 70 cm/s depending on flow conduit diameter. It is experimentally observed that by holding all other flow conditions and geometry constant, generated bubbles decrease in size with increase in superficial liquid velocity. Detached bubble diameter is shown to increase with air injection nozzle diameter. Likewise, generated bubbles grow in size with increasing pipe diameter. Along the same lines, it is shown that bubble frequency of formation increases and hence the time to detachment of a forming bubble decreases, as the superficial liquid velocity is in-creased. Furthermore, it is shown that the void fraction of the resulting two-phase flow increases with volumetric gas flow rate Q(sub d), pipe diameter and gas injection nozzle diameter, while they decrease with surrounding liquid flow. The important role played by flowing liquid in detaching bubbles in a reduced gravity environment is thus emphasized. We observe that the void fraction can be accurately controlled by using single nozzle gas injection, rather than by employing multiple port injection, since the later system gives rise to unpredictable coalescence of adjacent bubbles. It is of interest to note that empirical bubble size and corresponding void fraction are somewhat smaller for the co-flow geometry than the cross-flow configuration at similar flow conditions with similar pipe and nozzle diameters. In order to supplement the empirical data, a theoretical model is employed to study single bubble generation in the dynamic (Q(sub d) = 1 - 1000 cu cm/s) and bubbly flow regime within the framework of the co-flow configuration. This theoretical model is based on an overall force balance acting on the bubble during the two stages of generation, namely the expansion and the detachment stage. Two sets of forces, one aiding and the other inhibiting bubble detachment are identified. Under conditions of reduced gravity, gas momentum flux enhances, while the surface tension force at the air injection nozzle tip inhibits bubble detachment. In parallel, liquid drag and inertia can act as both attaching and detaching forces, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with our experimental results. However, at higher superficial liquid velocities, as the bubble loses its spherical form, empirical bubble size no longer matches the theoretical predictions. In summary, we have developed a combined experimental and theoretical work, which describes the complex process of bubble generation and resulting two-phase flow in a microgravity environment. Results of the present study can be used in a wide range of space-based applications, such as thermal energy and power generation, propulsion, cryogenic storage and long duration life support systems, necessary for programs such as NASA's Human Exploration for the Development of Space (HEDS).

Pais, S. C.↗
NASA NTRS2018

Towards an Enhanced Droplet Activation Scheme for Multi-Moment Bulk Microphysics Schemes

Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.

Loftus, Adrian M.↗
DOE OSTI2025

Prioritizing urban heat adaptation infrastructure based on multiple outcomes: Comfort, health, and energy

Globally, cities face increasing extreme heat, impacting comfort, health, and energy consumption. Infrastructure-based heat adaptation strategies can improve these outcomes, but each strategy has a unique mix of benefits and drawbacks. Here, we apply an urbanized meteorological model (WRF) with the newly integrated multilayer BEP-Tree street tree model to dynamically downscale Earth System Model projections and a 3-D microclimate model (TUF-Pedestrian) to simulate the street-scale radiation environment impacting pedestrians. We evaluate the performance of five heat adaptation strategies (street trees, cool roofs, green roofs, rooftop photovoltaics (PV), and reflective pavements) during extreme heat events in three cities with contrasting background climates (Toronto, Phoenix, and Miami), under contemporary and end-of-century projected climates, based on three metrics: outdoor heat stress, air conditioning (AC) energy use, and ventilation of vehicular air pollution. No single adaptation strategy improves all three outcomes. While street trees inhibit ventilation, they reduce outdoor heat stress four times more effectively than the next best strategy via shade provision, fully offsetting heat stress increases under a high-emissions end-of-century climate scenario in all cities studied. Cool roofs and green roofs moderately reduce heat stress and energy use. Alternatively, rooftop PV with energy storage can generate sufficient power for space cooling but have marginal effects on heat stress. Reflective pavements are the least effective across metrics. Where the ventilation of street-level emissions is of less concern, our results clearly support the combination of street trees and rooftop PV as a highly complementary and effective means of adaptive mitigation across different climates and neighborhood densities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2016

Passive Water Assurance Delivery System

This report provides an analysis and evaluation of proposed design improvements to the internal water delivery system of NASA’s passive vegetable growing system. The VEGGIE system is designed to provide fresh vegetables and psychological benefits for astronauts aboard the International Space Station. Preliminary testing conducted by NASA for the VEGGIE system revealed flaws in the water delivery system. In these tests, the interface between the plants and the water reservoir failed to provide passive water delivery through the system both in microgravity experiments and experiments conducted on Earth. Our team’s goal was to fix this interface such that water can be passively delivered from the water reservoir to the plant rooting pillows to grow vegetative crops from seed to harvest. The system requirements outlined by NASA include: minimal total mass, on-demand passive water delivery, same dimensions as existing fixtures, minimal pressure on system to prevent leaking, zero mold growth, minimal swelling or clogging for non-flammable capillary materials, pressure stabilization between the water reservoir and plant pillow bags, even dispersal of water during initial priming, minimal bubble obstruction of capillary interface, maximum gas availability to plant roots, growth in a mixed artificial media (50:50 Arcillite:Fafard #2), and to avoid overwatering and drought conditions for plants. The major proposed design improvements are to use a single-interface capillary cord design to directly connect the water reservoir to the plant-rooting pillows, and to alter the water reservoir to model a propellant management device (PMD) in order to ensure consistent and long term watering for the VEGGIE system. The plant rooting pillow required minimal changes outside of replacing the capillary mat on the bottom of the pillow with O-ring insertion points the single-interface capillary system. - Research was conducted to determine which materials are able to uptake water through capillary action to grow Outredgeous Lettuce plants from seed to harvest. The primary requirements to be met by the team's design and evaluation of the capillary interface were: continuous passive watering for 90 days using non-flammable capillary materials and a peak water delivery rate of 30 mL/hr./0.15 m 2 . The results of the team's experimentation showed that the capillary material Nomex displayed the highest capillary water delivery potential with a maximum flow rate of 3.6 mL/hr. The Nomex capillary systems were the only capillary material to consistently grow Outredgeous Lettuce plants from seed to harvest, and displayed the highest average flow rate for multiple experiment sets. The Nomex material was previously incorporated into the VEGGIE system using a matted version of the capillary material, though the team recommends using a cord configuration in the single-interface capillary system for greater system stability. Nomex has proven to be a promising material, as it has passed both health and fire standards for use aboard the International Space Station. PMDs are made of materials that utilize surface tension and adhesive forces to overcome adverse accelerations to improve stability and ensure fluid delivery. PMDs are typically used in fuel tanks to ensure fuel delivery. The team recommends using a sponge PMD in order to mitigate bubble obstruction, decrease system weight, and ensure reliable water delivery to the capillary interface. The PMD water reservoir requires a rigid water reservoir and a vent tube for pressure stabilization. Since PMD’s cannot be tested in 1-G (Earth conditions), further theoretical modeling and testing is required for the proposed water reservoir design. The intent of this proposed system is to passively water plants in microgravity, though the technology is not limited to microgravity applications. The testing at The Ohio State University has proven that the design is highly effective on Earth, demonstrating that it could serve as a simple water delivery system in home and office applications. This would make vegetative crops more accessible in all indoor applications, thereby improving indoor air quality and occupant comfort. Additionally, this technology has great potential to be utilized in greenhouse plant production, cutting back on more sophisticated watering system energy and time requirements.

Alexandria Jensen↗
NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗