Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive capture”

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.

53 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Formaldehyde Production from Isoprene Oxidation Across NOx Regimes

The chemical link between isoprene and formaldehyde (HCHO) is a strong, non-linear function of NOx (= 27 NO + NO2). This relationship is a linchpin for top-down isoprene emission inventory verification from orbital HCHO column observations. It is also a benchmark for overall mechanism performance with regard to VOC oxidation. Using a comprehensive suite of airborne in situ observations over the Southeast U.S., we quantify HCHO production across the urban-rural spectrum. Analysis of isoprene and its major first-generation oxidation products allows us to define both a "prompt" yield of HCHO (molecules of HCHO produced per molecule of freshly-emitted isoprene) and the background HCHO mixing ratio (from oxidation of longer-lived hydrocarbons). Over the range of observed NOx values (roughly 0.1 - 2 ppbv), the prompt yield increases by a factor of 3 (from 0.3 to 0.9), while background HCHO increases by more than a factor of 2 (from 1.5 to 3.3 ppbv). We apply the same method to evaluate the performance of both a global chemical transport model (AM3) and a measurement-constrained 0-D chemical box model. Both models reproduce the NOx dependence of the prompt HCHO yield, illustrating that models with updated isoprene oxidation mechanisms can adequately capture the link between HCHO and recent isoprene emissions. On the other hand, both models under-estimate background HCHO mixing ratios, suggesting missing HCHO precursors, inadequate representation of later-generation isoprene degradation and/or under-estimated hydroxyl radical concentrations. Moreover, we find that the total organic peroxy radical production rate is essentially independent of NOx, as the increase in oxidizing capacity with NOx is largely balanced by a decrease in VOC reactivity. Thus, the observed NOx dependence of HCHO mainly reflects the changing fate of organic peroxy radicals.

VOC↗
NASA NTRS2010

LES/RANS Simulation of a Supersonic Reacting Wall Jet

This work presents results from large-eddy / Reynolds-averaged Navier-Stokes (LES/RANS) simulations of the well-known Burrows-Kurkov supersonic reacting wall-jet experiment. Generally good agreement with experimental mole fraction, stagnation temperature, and Pitot pressure profiles is obtained for non-reactive mixing of the hydrogen jet with a non-vitiated air stream. A lifted flame, stabilized between 10 and 22 cm downstream of the hydrogen jet, is formed for hydrogen injected into a vitiated air stream. Flame stabilization occurs closer to the hydrogen injection location when a three-dimensional combustor geometry (with boundary layer development resolved on all walls) is considered. Volumetric expansion of the reactive shear layer is accompanied by the formation of large eddies which interact strongly with the reaction zone. Time averaged predictions of the reaction zone structure show an under-prediction of the peak water concentration and stagnation temperature, relative to experimental data and to results from a Reynolds-averaged Navier-Stokes calculation. If the experimental data can be considered as being accurate, this result indicates that the present LES/RANS method does not correctly capture the cascade of turbulence scales that should be resolvable on the present mesh. Instead, energy is concentrated in the very largest scales, which provide an over-mixing effect that excessively cools and strains the flame. Predictions improve with the use of a low-dissipation version of the baseline piecewise parabolic advection scheme, which captures the formation of smaller-scale structures superimposed on larger structures of the order of the shear-layer width.

Edwards, Jack R.↗
NASA NTRSDate not supplied

Predicting Lightning Initiation using Deep Learning

Lightning occurrence presents safety challenges to people and property. The main challenge with lightning safety is that the majority of guidance is reactive. In other words, lightning has to have already occurred nearby before a person will respond and take shelter. Further, most injuries or fatalities occur as the storm approaches, or as it's moving away, when rainfall may not be present at the time of the flash. Thus, this project develops a physically-based deep learning model to produce lightning probabilities out to 15 minutes. The deep learning model combines a Convolutional Neural Network (CNN) with a Long Short-Term Memory (LSTM) network to capture both the spatial and temporal evolution of storms to predict the probability that lightning initiation will occur in the next 15 minutes. The model combines radar reflectivity, correlation coefficient and differential reflectivity to inferred storm hydrometer type and precipitation phase, which aids in the identification of electrification processes. The model is trained with data from the Geostationary Lightning Mapper (GLM), which is a near infrared sensor onboard the GOES-R series of satellites that measures optical brightness from lightning. This presentation will provide an overview of the project.

Andrew T White↗
NASA NTRSDate not supplied

Advancing Autonomy in Distributed Space Systems: Insights From on-Orbit Testing with the Starling 1.0 Mission

Autonomous decision-making is crucial for enhancing mission effectiveness in Distributed Space Systems (DSS), particularly in multi-spacecraft operations where communication constraints and mission complexity pose challenges. The Distributed Spacecraft Autonomy (DSA) team at NASA’s Ames Research Center is advancing autonomy in DSS through five key technical areas: distributed resource and task management, reactive operations, system modeling and simulation, human-swarm interaction, and ad hoc network communications. The DSA experiment onboard the Starling 1.0 Mission showcases collaborative resource allocation for multi-point science data collection with four small spacecraft. Autonomy in decision-making is highlighted as a crucial factor for multi-spacecraft missions, enabling spacecraft to operate independently, reducing reliance on ground control. This capability is particularly significant for future deep-space missions, where communication delays and limited data transmission capacity make traditional command and control approaches impractical. This demonstration focuses on a GPS Channel Selection Experiment, leveraging emergent capabilities like "shared sampling" and "simultaneous sampling" to optimize channel selection across the spacecraft swarm. The experiment aims to capture ionospheric phenomena such as the Equatorial Ionization Anomaly and Polar Patches. The DSA system's autonomous reconfiguration ability is showcased, emphasizing its adaptability to natural phenomena without significant integration efforts. The GPS Channel Selection Experiment utilizes a dual-band GPS receiver to estimate plasma density in the ionosphere. Explorative and exploitative channel selections are employed based on the nature of observed phenomena. The performance of DSA algorithms is evaluated in terms of optimal channel allocations and responsiveness to changes in observed features. The DSA Flight Software utilizes the Core Flight System (cFS) framework, ensuring compatibility with the Starling 1.0 flight mission software. DSA showcases results from RTI’s Connext DDS Micro communication middleware, enabling message routing over the Ad-Hoc Network of Starling 1.0. This paper provides a comprehensive overview of the DSA experiment's initial results, emphasizing the advancements in autonomy for Distributed Space Systems and the successful collaboration with the Starling 1.0 mission.

Caleb Ashmore Adams↗
NASA NTRS2010

Reactive Sequencing for Autonomous Navigation Evolving from Phoenix Entry, Descent, and Landing

Virtual Machine Language (VML) is an award-winning advanced procedural sequencing language in use on NASA deep-space missions since 1997, and was used for the successful entry, descent, and landing (EDL) of the Phoenix spacecraft onto the surface of Mars. Phoenix EDL utilized a state-oriented operations architecture which executed within the constraints of the existing VML 2.0 flight capability, compatible with the linear "land or die" nature of the mission. The intricacies of Phoenix EDL included the planned discarding of portions of the vehicle, the complex communications management for relay through on-orbit assets, the presence of temporally indeterminate physical events, and the need to rapidly catch up four days of sequencing should a reboot of the spacecraft flight computer occur shortly before atmospheric entry. These formidable operational challenges led to new techniques for packaging and coordinating reusable sequences called blocks using one-way synchronization via VML sequencing global variable events. The coordinated blocks acted as an ensemble to land the spacecraft, while individually managing various elements in as simple a fashion as possible. This paper outlines prototype VML 2.1 flight capabilities that have evolved from the one-way synchronization techniques in order to implement even more ambitious autonomous mission capabilities. Target missions for these new capabilities include autonomous touch-and-go sampling of cometary and asteroidal bodies, lunar landing of robotic missions, and ultimately landing of crewed lunar vehicles. Close proximity guidance, navigation, and control operations, on-orbit rendezvous, and descent and landing events featured in these missions require elaborate abort capability, manifesting highly non-linear scenarios that are so complex as to overtax traditional sequencing, or even the sort of one-way coordinated sequencing used during EDL. Foreseeing advanced command and control needs for small body and lunar landing guidance, navigation and control scenarios, work began three years ago on substantial upgrades to VML that are now being exercised in scenarios for lunar landing and comet/asteroid rendezvous. The advanced state-based approach includes coordinated state transition machines with distributed decision-making logic. These state machines are not merely sequences - they are reactive logic constructs capable of autonomous decision making within a well-defined domain. Combined with the JPL's AutoNav software used on Deep Space 1 and Deep Impact, the system allows spacecraft to autonomously navigate to an unmapped surface, soft-contact, and either land or ascend. The state machine architecture enabled by VML 2.1 has successfully performed sampling missions and lunar descent missions in a simulated environment, and is progressing toward flight capability. The authors are also investigating using the VML 2.1 flight director architecture to perform autonomous activities like rendezvous with a passive hypothetical Mars sample return capsule. The approach being pursued is similar to the touch-and-go sampling state machines, with the added complications associated with the search for, physical capture of, and securing of a separate spacecraft. Complications include optically finding and tracking the Orbiting Sample Capsule (OSC), keeping the OSC illuminated, making orbital adjustments, and physically capturing the OSC. Other applications could include autonomous science collection and fault compensation.

asteroid↗
NASA NTRS2018

Decadal Changes in Summertime Reactive Oxidized Nitrogen and Surface Ozone over the Southeast United States

Widespread efforts to abate ozone (O3) smog have significantly reduced emissions of nitrogen oxides (NOx) over the past 2 decades in the Southeast US, a place heavily influenced by both anthropogenic and biogenic emissions. How reactive nitrogen speciation responds to the reduction in NOx emissions in this region remains to be elucidated. Here we exploit aircraft measurements from ICARTT (International Consortium for Atmospheric Research on Transport and Transformation - July-August 2004), SENEX (Southeast Nexus - June-July 2013), and SEAC4RS (Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys - August-September 2013) and long-term ground measurement networks alongside a global chemistry-climate model to examine decadal changes in summertime reactive oxidized nitrogen (RON) and ozone over the Southeast US. We show that our model can reproduce the mean vertical profiles of major RON species and the total (NO (sub y)) in both 2004 and 2013. Among the major RON species, nitric acid (HNO3) is dominant (approximately 42-45 percent), followed by NOx (31 percent), total peroxy nitrates (Sigma PNs; 14 percent), and total alkyl nitrates (Sigma ANs; 9-12 percent) on a regional scale. We find that most RON species, including NOx, Sigma PNs, and HNO3, decline proportionally with decreasing NOx emissions in this region, leading to a similar decline in NO (sub y). This linear response might be in part due to the nearly constant summertime supply of biogenic VOC (Volatile Organic Compounds) emissions in this region. Our model captures the observed relative change in RON and surface ozone from 2004 to 2013. Model sensitivity tests indicate that further reductions of NOx emissions will lead to a continued decline in surface ozone and less frequent high-ozone events.

Li, Jingyi↗
NASA NTRS2015

Fracture of Carbon Nanotube - Amorphous Carbon Composites: Molecular Modeling

Carbon nanotubes (CNTs) are promising candidates for use as reinforcements in next generation structural composite materials because of their extremely high specific stiffness and strength. They cannot, however, be viewed as simple replacements for carbon fibers because there are key differences between these materials in areas such as handling, processing, and matrix design. It is impossible to know for certain that CNT composites will represent a significant advance over carbon fiber composites before these various factors have been optimized, which is an extremely costly and time intensive process. This work attempts to place an upper bound on CNT composite mechanical properties by performing molecular dynamics simulations on idealized model systems with a reactive forcefield that permits modeling of both elastic deformations and fracture. Amorphous carbon (AC) was chosen for the matrix material in this work because of its structural simplicity and physical compatibility with the CNT fillers. It is also much stiffer and stronger than typical engineering polymer matrices. Three different arrangements of CNTs in the simulation cell have been investigated: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. The SWNT and MWNT array systems are clearly idealizations, but the SWNT bundle system is a step closer to real systems in which individual tubes aggregate into large assemblies. The effect of chemical crosslinking on composite properties is modeled by adding bonds between the CNTs and AC. The balance between weakening the CNTs and improving fiber-matrix load transfer is explored by systematically varying the extent of crosslinking. It is, of course, impossible to capture the full range of deformation and fracture processes that occur in real materials with even the largest atomistic molecular dynamics simulations. With this limitation in mind, the simulation results reported here provide a plausible upper limit on achievable CNT composite properties and yield some insight on the influence of processing conditions on the mechanical properties of CNT composites.

Jensen, Benjamin D.↗
NASA NTRS2009

A Summary of NASA and USAF Hypergolic Propellant Related Spills and Fires

Several unintentional hypergolic fluid related spills, fires, and explosions from the Apollo Program, the Space Shuttle Program, the Titan Program, and a few others have occurred over the past several decades. Spill sites include the following government facilities: Kennedy Space Center (KSC), Johnson Space Center (JSC), White Sands Test Facility (WSTF), Vandenberg Air Force Base (VAFB), Cape Canaveral Air Force Station (CCAFS), Edwards Air Force Base (EAFB), Little Rock AFB, and McConnell AFB. Until now, the only method of capturing the lessons learned from these incidents has been "word of mouth" or by studying each individual incident report. The root causes and consequences of the incidents vary drastically; however, certain "themes" can be deduced and utilized for future hypergolic propellant handling. Some of those common "themes" are summarized below: (1) Improper configuration control and internal or external human performance shaping factors can lead to being falsely comfortable with a system (2) Communication breakdown can escalate an incident to a level where injuries occur and/or hardware is damaged (3) Improper propulsion system and ground support system designs can destine a system for failure (4) Improper training of technicians, engineers, and safety personnel can put lives in danger (5) Improper PPE, spill protection, and staging of fire extinguishing equipment can result in unnecessary injuries or hardware damage if an incident occurs (6) Improper procedural oversight, development, and adherence to the procedure can be detrimental and quickly lead to an undesirable incident (7) Improper materials cleanliness or compatibility and chemical reactivity can result in fires or explosions (8) Improper established "back-out" and/or emergency safing procedures can escalate an event The items listed above are only a short list of the issues that should be recognized prior to handling hypergolic fluids or processing vehicles containing hypergolic propellants. The summary of incidents in this report is intended to cover many more issues than those listed above.

Nufer, Brian M.↗
NASA NTRS2007

Summertime Influence of Asian Pollution in the Free Troposphere over North America

We analyze aircraft observations obtained during INTEX-A (1 July 14 - August 2004) to examine the summertime influence of Asian pollution in the free troposphere over North America. By applying correlation analysis and Principal Component Analysis (PCA) to the observations between 6-12 km, we find dominant influences from recent convection and lightning (13 percent of observations), Asia (7 percent), the lower stratosphere (7 percent), and boreal forest fires (2 percent), with the remaining 71 percent assigned to background. Asian airmasses are marked by high levels of CO, O3, HCN, PAN, acetylene, benzene, methanol, and SO4(2-). The partitioning of reactive nitrogen species in the Asian plumes is dominated by peroxyacetyl nitrate (PAN) (approximately 600 pptv), with varying NO(x)/HNO3 ratios in individual plumes consistent with different plumes ages ranging from 3 to 9 days. Export of Asian pollution in warm conveyor belts of mid-latitude cyclones, deep convection, and lifting in typhoons all contributed to the five major Asian pollution plumes. Compared to past measurement campaigns of Asian outflow during spring, INTEX-A observations display unique characteristics: lower levels of anthropogenic pollutants (CO, propane, ethane, benzene) due to their shorter summer lifetimes; higher levels of biogenic tracers (methanol and acetone) because of a more active biosphere; as well as higher levels of PAN, NO(x), HNO3, and O3 (more active photochemistry possibly enhanced by injection of lightning NO(x)). The high delta O3/delta CO ratio (0.76 mol mol(exp -1)) of Asian plumes during INTEX-A is due to a combination of strong photochemical production and mixing with stratospheric air along isentropic surfaces. The GEOS-Chem global chemical transport model captures the timing and location of the Asian plumes remarkably well. However, it significantly underestimates the magnitude of the enhancements.

Liang, Q.↗
NASA NTRS2005

Capture of Hydrogen Using ZrNi

Water, as ice, is thought to reside in craters at the lunar poles along with CH4 and H2 . A proposed robotic mission for 2012 will utilize metal/metal hydrides for H2 recovery. Specifications are 99% capture of H2 initially at 5 bar and 100C (or greater), and degassing completely at 300C. Of 47-systems examined using the van't Hoff equation, 4 systems, Mg/MgH2, Mg2Ni/Mg2NiH4, ZrNi/ZrNiH2.8, and Pd/PdH0.77, were considered likely candidates for further examination. It is essential, when selecting a system, to also examine questions regarding activation, kinetics, cyclic stability, and gas impurity effects. After considering those issues, ZrN1 was selected as the most promising candidate, as it is easily activated and rapidly forms ZrNiH 2.8 . In addition, it resists oxide poisoning by CO2, and H2O, while some oxidation by O2 is recommended for improved activation . The presence of hydrogen in the as received Zr-Ni alloy from Alfa Aesar posed additional technical problems. X-ray diffraction of the Zr-Ni powder (-325 mesh), with a Zr:Ni wt% ratio of 70:30, was found to consist of ZrH2, ZrNiH2.8, and ZrNi. ZrH2 in the alloy presented the risk that after degassing that both Zr and ZrNi would be present, and thus lead to erroneous results regarding the reactivity of ZrNi with H2 . Fortunately, ZrH2 is a highly stable hydride that does not degas H2 to any significant extent at temperatures below 300C. Based on equilibrium calculations for the decomposition of ZrH2, only 1 millionth of the hydride decomposed at 300C under a N2 atmosphere flowing at 25 ccm for 64 hours, the longest time for pretreatment employed in the investigation. It was possible, from the X-ray results and knowledge of the Zr:Ni ratio, to compute the composition of a pretreated specimen as being 76 wt% ZrNi and the balance ZrH2.

Patton, Lisa↗
NASA NTRS2006

Opposed Jet Burner Approach for Characterizing Flameholding Potentials of Hydrocarbon Scramjet Fuels

Opposed Jet Burner (OJB) tools have been used extensively by the authors to measure Flame Strength (FS) extinction limits of laminar H2/N2 air and (recently) hydrocarbon (HC) air Counterflow Diffusion Flames (CFDFs) at one atm. This paper details normalization of FSs of N2- diluted H2 and HC systems to account for effects of fuel composition, temperature, pressure, jet diameter, inflow Reynolds number, and inflow velocity profile (plug, contoured nozzle; and parabolic, straight tube). Normalized results exemplify a sensitive accurate means of validating, globally, reduced chemical kinetic models at approx. 1 atm and the relatively low temperatures approximating the loss of non-premixed idealized flameholding, e.g., in scramjet combustors. Laminar FS is defined locally as maximum air input velocity, U(sub air), that sustains combustion of a counter-jet of g-fuel at extinction. It uniquely characterizes a fuel. And global axial strain rate at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or (sub t)) can be compared directly with computed extinction limits, determined using either a 1-D Navier Stokes stream-function solution, using detailed transport and finite rate chemistry, or (better yet) a detailed 2-D Navier Stokes numerical simulation. The experimental results define an idealized flameholding reactivity scale that shows wide ranging (50 x) normalized FS s for various vaporized-liquid and gaseous HCs, including, in ascending order: JP-10, methane, JP-7, n-heptane, n-butane, propane, ethane, and ethylene. Results from H2 air produce a unique and exceptionally strong flame that agree within approx. 1% of a recent 2-D numerically simulated FS for a 3 mm tube-OJB. Thus we suggest that experimental FS s and/or FS ratios, for various neat and blended HCs w/ and w/o additives, offer accurate global tests of chemical kinetic models at the Ts and Ps of extinction. In conclusion, we argue the FS approach is more direct and fundamental, for assessing, e.g., idealized scramjet flameholding potentials, than measurements of laminar burning velocity or blowout in a Perfectly Stirred Reactor, because the latter characterize premixed combustion in the absence of aerodynamic strain. And FS directly measures a chemical kinetic characteristic of non-premixed combustion at typical flameholding temperatures. It mimics conditions where gfuels are typically injected into a subsonic flameholding recirculation zone that captures air, where the effects of aerodynamic strain and associated multi-component diffusion become important.

Pellett, Gerald L.↗
NASA NTRS2017

Using Advanced Mass Spectrometry Techniques to Fully Characterize Atmospheric Organic Carbon: Current Capabilities and Remaining Gaps

Organic compounds in the atmosphere vary widely in their molecular composition and chemical properties, so no single instrument can reasonably measure the entire range of ambient compounds. Over the past decade, a new generation of in-situ, field-deployable mass spectrometers has dramatically improved our ability to detect, identify, and quantify these organic compounds, but no systematic approach has been developed to assess the extent to which currently available tools capture the entire space of chemical identity and properties that is expected in the atmosphere. Reduced-parameter frameworks that have been developed to describe atmospheric mixtures are exploited here to characterize the range of chemical properties accessed by a suite of instruments. Multiple chemical spaces (e.g. oxidation state of carbon vs. volatility, and oxygen number vs. carbon number) were populated with ions measured by several mass spectrometers, with gas- and particle-phase α-pinene oxidation products serving as the test mixture of organic compounds. Few gaps are observed in the coverage of the parameter spaces by the instruments employed in this work, though the full extent to which comprehensive measurement was achieved is difficult to assess due to uncertainty in the composition of the mixture. Overlaps between individual ions and regions in parameter space were identified, both between gas- and particle-phase measurements, and within each phase. These overlaps were conservatively found to account for little (<10%) of the measured mass. However, challenges in identifying overlaps and in accurately converting molecular formulas into chemical properties (such as volatility or reactivity) highlight a continued need to incorporate structural information into atmospheric measurements.

G. Isaacman-VanWertz↗
NASA NTRSDate not supplied

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗
NASA NTRSDate not supplied

Improved Assessment of Recent Trends in NOx and VOC Emissions and Ozone Production Sensitivity Regimes Using Satellite Data

This presentation highlights results from a NASA Aura Science Team and Atmospheric Composition Modeling and Analysis Program (ACMAP) project which study the capability to observe and model trends in ozone (O3) production regimes using spaceborne sensors. Ultraviolet– visible (UV–Vis) tropospheric column satellite retrievals of formaldehyde (HCHO) (a proxy for volatile organic compound [VOC] reactivity) and nitrogen dioxide (NO2) (a proxy for nitrogen oxides [NOx]) are frequently used to investigate the sensitivity of O3 production to emissions of NOx and VOCs. There are challenges that come from using satellite-derived ratios of HCHO and NO2 (FNR) to study O3 production sensitivity with the largest uncertainties associated with specific spaceborne sensor’s retrieval biases and errors. This study quantifies the differences and improvements in satellite retrievals of O3 production sensitivity regimes using FNRs when moving from legacy polar orbiting satellites such as the Ozone Monitoring Instrument (OMI) onboard NASA’s Aura satellite and Ozone Mapping and Profiler Suite Nadir Mapper (OMPS-NM) onboard the NASA/NOAA Suomi-NPP platform to newer, higher spatiotemporal resolution satellite sensors TROPOspheric Monitoring Instrument (TROPOMI) and eventually the recently launched NASA geostationary sensor Tropospheric Emissions: Monitoring of Pollution (TEMPO). Furthermore, we investigate how using retrievals of NO2 and HCHO from these different satellites to constrain model predictions impacts the ability to accurately simulate O3 chemistry including chemical production regimes. To this end, we have conducted inverse model simulations, using the WRF-CMAQ-DDM data assimilation system at 12 km × 12 km, to constrain emissions of NOx and VOCs over the contiguous United States (CONUS) when assimilating OMI and TROPOMI retrievals of NO2 and HCHO. Two advantages of this are that we a) account for each satellite’s errors/biases in the emission estimation and b) update the prior profile to ensure that only radiance information is used for optimizing the emissions. This presentation will demonstrate: a) the varying accuracy of different satellite retrieved FNRs and ability to capture known sub-annual emission trends (e.g., seasonal, weekend/weekday) and emission anomalies during the COVID-19 lockdown of 2020, b) the differences and improvements in top-down emission estimates of NOx and VOCs when constrained by newer satellite sensors compared to legacy systems, and c) multi-sensor optimized emission estimates of summer-time NOx and VOCs between 2019-2021.

Data↗
NASA NTRS2023

Characterization of Errors in Satellite-Based HCHO/NO 2 Tropospheric Column Ratios With Respect to Chemistry, Column-to-PBL Translation, Spatial Representation, and Retrieval Uncertainties

The availability of formaldehyde (HCHO) (a proxy for volatile organic compound reactivity) and nitrogen dioxide (NO 2 ) (a proxy for nitrogen oxides) tropospheric columns from ultraviolet–visible (UV–Vis) satellites has motivated many to use their ratios to gain some insights into the near-surface ozone sensitivity. Strong emphasis has been placed on the challenges that come with transforming what is being observed in the tropospheric column to what is actually in the planetary boundary layer (PBL) and near the surface; however, little attention has been paid to other sources of error such as chemistry, spatial representation, and retrieval uncertainties. Here we leverage a wide spectrum of tools and data to quantify those errors carefully. Concerning the chemistry error, a well-characterized box model constrained by more than 500 h of aircraft data from NASA's air quality campaigns is used to simulate the ratio of the chemical loss of HO 2 + RO 2 (LRO x ) to the chemical loss of NO x (LNO x ). Subsequently, we challenge the predictive power of HCHO/NO 2 ratios (FNRs), which are commonly applied in current research, in detecting the underlying ozone regimes by comparing them to LRO x /LNO x . FNRs show a strongly linear (R 2 =0.94) relationship with LRO x /LNO x , but only on the logarithmic scale. Following the baseline (i.e., ln(LRO x /LNO x ) = −1.0 ± 0.2) with the model and mechanism (CB06, r2) used for segregating NO x -sensitive from VOC-sensitive regimes, we observe a broad range of FNR thresholds ranging from 1 to 4. The transitioning ratios strictly follow a Gaussian distribution with a mean and standard deviation of 1.8 and 0.4, respectively. This implies that the FNR has an inherent 20 % standard error (1σ) resulting from not accurately describing the RO x –HO x cycle. We calculate high ozone production rates (PO 3 ) dominated by large HCHO × NO 2 concentration levels, a new proxy for the abundance of ozone precursors. The relationship between PO 3 and HCHO × NO 2 becomes more pronounced when moving towards NO x -sensitive regions due to nonlinear chemistry; our results indicate that there is fruitful information in the HCHO × NO 2 metric that has not been utilized in ozone studies. The vast amount of vertical information on HCHO and NO 2 concentrations from the air quality campaigns enables us to parameterize the vertical shapes of FNRs using a second-order rational function permitting an analytical solution for an altitude adjustment factor to partition the tropospheric columns into the PBL region. We propose a mathematical solution to the spatial representation error based on modeling isotropic semivariograms. Based on summertime-averaged data, the Ozone Monitoring Instrument (OMI) loses 12 % of its spatial information at its native resolution with respect to a high-resolution sensor like the TROPOspheric Monitoring Instrument (TROPOMI) (> 5.5 × 3.5 km 2 ). A pixel with a grid size of 216 km 2 fails at capturing ∼ 65 % of the spatial information in FNRs at a 50 km length scale comparable to the size of a large urban center (e.g., Los Angeles). We ultimately leverage a large suite of in situ and ground-based remote sensing measurements to draw the error distributions of daily TROPOMI and OMI tropospheric NO 2 and HCHO columns. At a 68 % confidence interval (1σ), errors pertaining to daily TROPOMI observations, either HCHO or tropospheric NO 2 columns, should be above 1.2–1.5 × 10 16 molec. cm −2 to attain a 20 %–30 % standard error in the ratio. This level of error is almost non-achievable with the OMI given its large error in HCHO. The satellite column retrieval error is the largest contributor to the total error (40 %–90 %) in the FNRs. Due to a stronger signal in cities, the total relative error (< 50 %) tends to be mild, whereas areas with low vegetation and anthropogenic sources (e.g., the Rocky Mountains) are markedly uncertain (> 100 %). Our study suggests that continuing development in the retrieval algorithm and sensor design and calibration is essential to be able to advance the application of FNRs beyond a qualitative metric.

Amir H. Souri↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗