Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “comparable input ranges”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulations of short-term exposure to NO 2 and PM 2.5 to inform capture efficiency standards

The California building code requires all new or renovated residential dwelling units to have kitchen exhaust ventilation to manage air pollutants and moisture generated during cooking. Current performance requirements are specified as a minimum airflow and maximum sound level. This report presents an analysis to support consideration of adding a capture efficiency requirement to the code. The analysis uses a physics-based simulation model to calculate air pollutant concentrations in homes that result from cooking and the inflow of outdoor air, mediated by loss and removal process including deposition, dwelling unit ventilation, and the use of a range hood with varying levels of capture efficiency. Calculated pollutant concentrations are compared to relevant, health-based guidelines. The analysis considers the highest 1-hour concentration of nitrogen dioxide (NO 2 ), which is a product of natural gas combustion and the highest 24-hour concentrations of fine particulate matter (PM 2.5 ) which is emitted in substantial quantities from frying, broiling and grilling among other cooking activities. For NO 2 , the analysis considers cooking of a dinner for 3-4 persons consisting of pasta, meat sauce, a par-boiled vegetable, and baked garlic bread. For PM 2.5 , the analysis considers a day in which breakfast emitted particles at the 80th percentile and lunch and dinner entailed particle emissions at the 50th percentile of cooking emission events reported in the literature. Model simulations were performed to determine the level of range hood capture efficiency that will allow these cooking scenarios to occur in the vast majority (>99%) of new homes being built in California while maintaining pollutant concentrations below the health-based guidelines, if the range hood is used throughout cooking. All homes were assumed to have dwelling unit ventilation at the rate required in the building code. Simulation model input parameters were specified using a Monte Carlo approach to represent a range of housing characteristics, outdoor conditions, and indoor pollutant dynamics. Simulation results suggest that requiring a minimum capture efficiency of at least 70% is needed to avoid unacceptably high NO 2 (1-h average concentration of 100 ppb or higher) and at least 60% to avoid unacceptably high PM 2.5 (24-h average of 25 μg/m³ or higher). These results were driven by multi-family homes, which have a smaller volume of air to dilute any pollutants not captured and removed at the cooking area, resulting in higher concentrations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Thermal emission from particulate surfaces: A comparison of scattering models with measured spectra

Emissivity spectra of particulate mineral samples are highly dependent on particle size when that size is comparable to the wavelength of light emitted (5-50 micrometers for the midinfrared). Proper geologic interpretation of data from planetary infrared spectrometers will require that these particle size effects be well understood. To address this issue, samples of quartz powders were produced with narrow, well-characterized particle size distributions. Mean particle diameters in these samples ranged from 15 to 227 micrometers. Emission spectra of these powders allow the first detailed comparison of the complex spectral variations with particle size observed in laboratory data with the predictions of radiative transfer models. Four such models are considered here. Hapke's relectance theory (converted to emissivity via Kirchoff's law) is the first model tested. Hapke's more recently published emission theory is also employed. The third model, the 'Mie/Conel' model, uses Mie single scattering with a two-stream approximation for multiple scattering. This model, like the first, is a converted reflec- tance model. Mie scattering assumes particles are both spherical and well separated, which is not true for the quartz powders, but includes diffraction effects. The fourth model uses the Mie solution for single scattering by spheres and inputs those results into the multiple scattering formalism of Hapke's emission theory. The results of the four models are considered in relation to the values of the optical constants n and k. We have grouped these as class 1 (k large), class 2 (k moderate, n is approximately 2), class 3 (k small, n is approximately 2), and class 4 (k small, n is approximately 1). In general, the Mie/Hapke hybrid model does best at predicting variations with grain size. In particular, it predicts changes of the correct pattern, although incorrect magnitude, for class 1 bands, where large increases in emissivity with decreasing grain size are observed. This model also does an excellent job on moderate (class 2) and very weak and intraband (class 3) regions, and correctly predicts the emission maximum and its invariance with grain size near the Christiansen frequency (class 4). The Mie/Hapke hybrid model also has the fewest free parameters of the four models examined, while maintaining the most physical treatment of the radiative transfer. The Mie/Conel model performs as well as the Mie/Hapke hybrid model in strong bands (class 1) but does not accurately model the behavior of moderate (class 2) and very weak (class 3) bands.

Moersch, J. E.↗
DOE OSTI2023

Dimensionally reduced machine learning model for predicting single component octanol–water partition coefficients

Abstract MF-LOGP, a new method for determining a single component octanol–water partition coefficients ( $$LogP$$ LogP ) is presented which uses molecular formula as the only input. Octanol–water partition coefficients are useful in many applications, ranging from environmental fate and drug delivery. Currently, partition coefficients are either experimentally measured or predicted as a function of structural fragments, topological descriptors, or thermodynamic properties known or calculated from precise molecular structures. The MF-LOGP method presented here differs from classical methods as it does not require any structural information and uses molecular formula as the sole model input. MF-LOGP is therefore useful for situations in which the structure is unknown or where the use of a low dimensional, easily automatable, and computationally inexpensive calculations is required. MF-LOGP is a random forest algorithm that is trained and tested on 15,377 data points, using 10 features derived from the molecular formula to make $$LogP$$ LogP predictions. Using an independent validation set of 2713 data points, MF-LOGP was found to have an average $$RMSE$$ RMSE = 0.77 ± 0.007, $$MAE$$ MAE = 0.52 ± 0.003, and $${R}^{2}$$ R 2 = 0.83 ± 0.003. This performance fell within the spectrum of performances reported in the published literature for conventional higher dimensional models ( $$RMSE$$ RMSE = 0.42–1.54, $$MAE$$ MAE = 0.09–1.07, and $${R}^{2}$$ R 2 = 0.32–0.95). Compared with existing models, MF-LOGP requires a maximum of ten features and no structural information, thereby providing a practical and yet predictive tool. The development of MF-LOGP provides the groundwork for development of more physical prediction models leveraging big data analytical methods or complex multicomponent mixtures. Graphical Abstract

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
NASA NTRS2006

Binocular Coordination of the Human Vestibulo-Ocular Reflex during Off-axis Pitch Rotation

Head movements in the sagittal pitch plane typically involve off-axis rotation requiring both vertical and horizontal vergence ocular reflexes to compensate for angular and translational motion relative to visual targets of interest. The purpose of this study was to compare passive pitch VOR responses during rotation about an Earth-vertical axis (canal only cues) with off-axis rotation (canal and otolith cues). Methods. Eleven human subjects were oscillated sinusoidally at 0.13, 0.3 and 0.56 Hz while lying left-side down with the interaural axis either aligned with the axis of rotation or offset by 50 cm. In a second set of measurements, twelve subjects were also tested during sinusoidally varying centrifugation over the same frequency range. The modulation of vertical and horizontal vergence ocular responses was measured with a binocular videography system. Results. Off-axis pitch rotation enhanced the vertical VOR at lower frequencies and enhanced the vergence VOR at higher frequencies. During sinusoidally varying centrifugation, the opposite trend was observed for vergence, with both vertical and vergence vestibulo-ocular reflexes being suppressed at the highest frequency. Discussion. These differential effects of off-axis rotation over the 0.13 to 0.56 Hz range are consistent with the hypothesis that otolith-ocular reflexes are segregated in part on the basis of stimulus frequency. At the lower frequencies, tilt otolith-ocular responses compensate for declining canal input. At higher frequencies, translational otolith-ocular reflexes compensate for declining visual contributions to the kinematic demands required for fixating near targets.

Wood, S. J.↗
DOE OSTI2020

Dynamic Undervolting to Improve Energy Efficiency on Multicore X86 CPUs

Chip manufacturers introduce redundancy at various levels of CPU design to guarantee correct operation, even for worst-case combinations of non-idealities in process variation and system operating conditions. This redundancy is implemented partly in the form of voltage margins. However, for a wide range of real-world execution scenarios these margins are excessive and merely translate to increased power consumption, hindering the effort towards higher-energy efficiency in both HPC and general purpose computing. Our study on the x86-64 Haswell and Skylake multicore microarchitectures reveals-wide voltage margins, which vary across different microarchitectures, different chip parts of the same microarchitecture, and across different workloads. We find that it is necessary to quantify-voltage margins using multi-threaded and multi-instance workloads, as characterization with single-threaded and single-instance workloads that do not stress the CPU to its full capacity typically identifies overly optimistic margins that lead to errors when applied in realistic program execution scenarios. In addition, we introduce, deploy and evaluate a run-time governor that dynamically reduces the supply voltage of modern multicore x86-64 CPUs. Our governor employs a model that takes as input a set of performance metrics which are directly measurable via performance monitoring counters and have high predictive value for the minimum tolerable supply voltage (V min ), to predict and apply the appropriate reduction for the workload at hand. Compared with the conventional DVFS governor, our approach in this study achieves up to 42 percent energy savings for the Skylake family and 34 percent for the Haswell family for complex, real-world applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗
NASA NTRS2021

Sky-Scanning Sun-Tracking Airborne Radiometer (3STAR): Instrument Design, Flight Testing, and Tracking Performance

The Sky-Scanning, Sun-Tracking Airborne Radiometer (3STAR) adapts commercial radiometer technology developed for the ocean color research community to airborne measurement of spectrally resolved solar irradiance and sky radiance. These atmospheric observations are used to derive aerosol optical depth (AOD), supporting localized AOD inputs for atmospheric correction of satellite and airborne data over terrestrial and aquatic (including optically dark) targets. The ability to regionally “tune” atmospheric correction schemes with relevant spatial AOD supports constraining atmospheric correction of remote sensing reflectance. Very wide dynamic range has been achieved for multi-channel band-pass-filter-radiometers originally designed for deployment into the water column. By actively tracking and directly pointing to the Sun, the light attenuation by aerosol particles in the atmospheric column can be quantified. These measurements improve knowledge of atmospheric constituents and the atmospheric corrections required to improve remote sensing capabilities for interpreting reflectance measurements from the Earth surface. 3STAR incorporates a custom Sun-tracking/sky-scanning pointing head, a Sun-tracking camera, and a commercially available, cylindrical, 19-channel radiometer tube assembly. An accurate and responsive mount and tracking system has been developed and certified to NASA and Naval Air Systems Command (NAVAIR) airworthiness standards for deployment into the aircraft slipstream. Ground and flight testing indicate typical tracking errors of less than 0.1 degrees, well within the field of view of the radiometer as required to minimize measurement uncertainty from alignment error. Preliminary AOD measurements compare to within 0.013 with 15 measurements from an Aerosol Robotic Network (AERONET) Cimel instrument at 500 nm wavelength and low solar angle.

Atmosphere↗
DOE OSTI2024

Shorter function summaries for finite state machine-based high consequence systems using logic synthesis and tautologies (Final Report LDRD 24-1302)

Computer programs are often viewed as collections of functions – each function has parameters (inputs) and computes a return value, and each has potential side effects that modify program state (outputs). In this research, a Sandia symbolic execution tool designed to support “human-in-the-loop” analysis was modified to automatically create “function summaries,” and a new tool, “diaboolical,” was created to support enhancing readability of the summary using a novel approach to bit-vector simplification that leverages logic synthesis and tautologies. For this effort, students at Auburn University created several finite state machines (FSMs) to serve as exemplars for high-consequence systems. Function summaries for each of the machines were obtained, and then portions of the summaries were simplified using both diaboolical and the simplification procedure of a popular SMT solver. A comparison of the results shows that diaboolical can often produce smaller function summaries, with expression length improvements over the unsimplified function summaries ranging from 0% to 90% for diaboolical and 0% to 65% for the SMT solver, though diaboolical had a significantly greater cost in time. Diaboolical was evaluated against a collection of “arbitrary” C-code as well as FSM exemplars, and for both datasets it achieved an approximately 10% improvement in expression length compared to simplifications that could be obtained using existing techniques. Function summaries can assist assurance efforts that evaluate existing systems and their executable code. A smaller function summary is likely easier for humans to understand and could thus increase the ability and efficacy of assurance practices centered around the analysis of executable artifacts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Testing mechanisms of how mycorrhizal associations affect forest soil carbon and nitrogen cycling (Final Technical Report)

Trees are in a symbiotic partnership with mycorrhizal fungi in which they provide the fungi with carbon from photosynthesis and the fungi provide the trees with nutrients and water. In temperate forests, the vast majority of trees form symbioses with one of two types of mycorrhizal fungi—arbuscular mycorrhizal (AM) fungi or ectomycorrhizal (EcM) fungi. These fungi differ in their morphology, hyphal length, and nutrient acquisition strategies. Many studies have found systematic differences in soil organic matter and nitrogen availability between forest stands dominated by AM-associating trees versus EcM-associating trees. For instance, there is a larger proportion of organic matter that is mineral-associated, more available nitrogen, and lower soil carbon to nitrogen ratios in AM forest stands relative to EcM forests stands. However, the mechanisms driving these patterns are not known, which complicates our ability to model soil organic matter dynamics in forested ecosystems. The main objective of this research was to understand the degree to which the observed differences in soil C and N dynamics between AM and EcM dominated forests are driven by tree traits like litter decomposability and root exudation versus mycorrhizal fungal nutrient acquisition strategies. We investigated these mechanisms using observations and targeted experiments and incorporated this knowledge into a process-based soil organic matter model. The observational studies compared the importance of leaf litter decomposability versus fungal identity on soil organic matter processes. We found that often fungal identity and traits were more important drivers of soil organic matter patterns than leaf litter decomposability. We ran two novel experiments: 1) a growth chamber experiment across four EcM and four AM tree species using a 13 C-labeled atmosphere to trace seedling-derived C into hyphae, the rhizosphere, and soil; and 2) an in situ decomposition experiment of six different 13 C and 15 N labeled litters that ranged in decomposability incubated across a gradient of EcM dominance at three sites that capture important variation in climate, soils, and forest species composition. In the first experiment, we found no significant differences of seedling mycorrhizal association on soil carbon sequestration over a growing season, but we did find that mycorrhizal association affected rhizodeposition with EcM-associating seedlings depositing more carbon in response to increased nitrogen availability. The decomposition experiment is still ongoing, but thus far, we have found slower litter decomposition in only one of three EcM-dominated forests which suggests that differences between AM- and EcM-dominated forests depend on the environmental context and identity of the EcM fungi. Lastly, we explicitly incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model creating Myco-CORPSE. By including the different nutrient acquisition strategies of AM and EcM fungi, we explored the conditions under which EcM fungi can slow decomposition rates and lead to greater soil organic carbon accumulation compared to AM fungi. We found that the effect of EcM fungi was highly context dependent and that EcM fungi decreased decomposition in colder forests with recalcitrant litter inputs and when they produced oxidases and necromass-degrading enzymes. Our research highlights the importance of fungal nutrient acquisition in driving soil organic matter patterns and the need to move beyond the AM-EcM dichotomy to consider the identity and traits of the specific fungi participating in the symbiosis. Overall, this research has resulted in six, peer-reviewed published papers in journals such as Global Change Biology, Ecology (2), Soil Biology and Biochemistry, and Ecosystems (2). There are at least two more papers in progress on this research including one that was recently submitted to Global Change Biology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Energy Budget of Brittle Fracturing in Granite Under Stress Relaxation and Creep

Creep and relaxation are the two major time-dependent fracturing processes in rocks. While a considerable amount of research has been done in understanding these two mechanisms, critical gaps remain regarding how different energy components evolve during time-dependent fracturing processes in rocks. In this study, a series of relaxation and creep experiments were conducted on prismatic Barre granite specimens in the laboratory to estimate the energy budget of brittle fracturing in granite. For the input energy, the work done by the machine (W) is calculated and for the output energy the radiated seismic energy (E R ), released in the form of acoustic emissions (AEs), is calculated as the only measurable output energy component in the conducted experiments. The low-frequency plateau (Ω 0 ) and corner frequency (f 0 ) for each AE waveform was estimated by fitting the observed AE spectra with the theoretical spectra using the Omega model. These parameters were used to estimate the seismic moments (M 0 ) based on the radiation pattern for the double couple (shear) and non-double-couple (non-shear) events. The range of f 0 and M 0 varied from 150 to 750 kHz and 10 –4 to 10 –1 N m, respectively. Moment magnitude (M w ) varied in a wider range from –9 to –6 in creep and –8.5 to –7 in relaxation. Stress drops (Δ σ ) and source radius (r) were estimated for the AEs using Brune’s model. Here, the results report on three primary observations: (1) the effects of different source mechanisms on the estimated source parameters showed that M 0 and Δ σ were higher for DC events as compared to NDC in both relaxation and creep. (2) The radiation efficiency in the case of creep is 70% higher as compared to relaxation and, (3) the stress drop estimated in relaxation and creep demonstrated a breakdown in scaling with the seismic moment.

58 GEOSCIENCES↗
DOE OSTI2020

Vertex Reordering for Real-world Graphs and Applications: An Empirical Evaluation

Vertex reordering is a way to improve locality in graph computations. Given an input (or ``natural'') order, reordering aims to compute an alternate permutation of the vertices that is aimed at maximizing a locality-based objective. Given decades of research on this topic, there are tens of graph reordering schemes, and there are also several linear arrangement ``gap'' measures for treatment as objectives. However, a comprehensive empirical analysis of the efficacy of the ordering schemes against the different gap measures, and against real-world applications is currently lacking. In this study, we present an extensive empirical evaluation of up to 11 ordering schemes, taken from different classes of approaches, on a set of 34 real-world graphs emerging from different application domains. Our study is presented in two parts: a) a thorough comparative evaluation of the different ordering schemes on their effectiveness to optimize different linear arrangement gap measures, relevant to preserving locality; and b) extensive evaluation of the impact of the ordering schemes on two real-world, parallel graph applications, namely, community detection and influence maximization. Our studies show a significant divergence among the ordering schemes (up to $40\times$ between the best and the poor) in their effectiveness to reduce the gap measures; and a wide ranging impact of the ordering schemes on various aspects including application runtime (up to $4\times$), memory and cache use, load balancing, and parallel work and efficiency. The comparative study also help reveal the nuances of a parallel environment (compared to serial) on the ordering schemes and their role in optimizing applications.

Barik, Reet↗
DOE OSTI2023

A Deep Learning Pipeline for Optimizing Large-scale Phase Field Simulations

Phase field (PF) simulations are computationally expensive but remain a key analysis tool to understand the complex mechanisms of additive manufacturing (AM) processes. Each PF simulation-aided analysis requires thousands of node hours on leadership-class supercomputers. One of the main goals of these analyses is the study of microstructure evolution during the build process which begins with the onset of nucleation. Nucleation occurs under certain thermomechanical conditions which are not known a priori and many PF simulations are required to identify ranges of input thermo-mechanical parameters that can result in the onset of nucleation. Since many of the simulations do not result in nucleation, an analysis campaign often ends up wasting tremendous amounts of precious computing resources executing nucleation-absent simulations. The goal of this work is to design and train deep learning models to inform a PF simulation about the likelihood of the occurrence of nucleation in a future simulation time-step based on the state summary over a finite number of past time-steps of a running simulation. If the prediction determines that the running simulation is unlikely to reach nucleation in the allotted time, then its execution is stopped immediately ultimately resulting in vast reduction in wasted computations when accrued over all the PF simulations typically performed in a single or multiple analysis campaign(s). The paper presents the performance of a machine learning pipeline that uses a convolutional neural network (CNN) model to learn an embedding which is then used with a self-attention network to build a multi-task deep learning model to predict the likelihood of nucleation. The model also predicts the input parameters used in a simulation. Performance is compared with a baseline pipeline that uses an off-the-shelf LeNet-5 model to learn the initial embedding. Despite their smaller size, performance results indicate significant improvement in accuracy of the proposed models compared to the larger baseline models.

Kannan, Ramakrishnan {ramki}↗
NASA NTRS2021

Modeling Snow Ablation over the Western United States Mountains: Patterns and Controlling Factors

When compared with differences in snow accumulation predicted by widely used hydrological models, there is a much greater divergence among otherwise “good” models in their simulation of the snow ablation process. Here, we explore differences in the performance of the Variable Infiltration Capacity model (VIC), Noah land surface model with multiparameterization options (Noah-MP), the Catchment model, and the third-generation Simplified Simple Biosphere model (SiB3) in their ability to reproduce observed snow water equivalent (SWE) during the ablation season at 10 Snowpack Telemetry (SNOTEL) stations over 1992–2012. During the ablation period, net radiation generally has stronger correlations with observed melt rates than does air temperature. Average ablation rates tend to be higher (in both model predictions and observations) at stations with a large accumulation of SWE. The differences in the dates of last snow between models and observations range from several days to approximately a month (on average 5.1 days earlier than in observations). If the surface cover in the models is changed from observed vegetation to bare soil in all of the models, only the melt rate of the VIC model increases. The differences in responses of models to canopy removal are directly related to snowpack energy inputs, which are further affected by different algorithms for surface albedo and energy allocation across the models. We also find that the melt rates become higher in VIC and lower in Noah-MP if the shrub/grass present at the observation sites is switched to trees.

Snow↗
NASA NTRS2011

TPSAS-NF1676L-12387-DND

Fast-time wake prediction models are semi-empirical, and are designed to predict aircraft wake vortex information behind a generating aircraft. They are developed using theoretical concepts and often calibrated with data obtained from field measurements, laboratory studies, or numerical simulations. They rely upon inputs representing the generating aircraft (such as weight, airspeed, wingspan, lateral position, and altitude) and atmospheric environment (such as ambient winds, turbulence, and temperature lapse rate). Other input may include the type of aircraft (e.g., B-747- 400). Output from these models may include time histories of the wake vortex strength, and position, as well as a prediction of the uncertainties in strength and position. Fast-time models should provide solutions in a fraction of a second, and hence is an enabling technology able to perform in a real-time operational system, such as dynamic spacing systems envisioned by the Aircraft Vortex Spacing System or other Super Density Operations concepts. Other applications of fast-time models include 1) guidance for wake separation standards for new aircraft, 2) optimization of separation standards for existing aircraft, 3) guidance for setting vertical separation standards during cruise, 4) accident reconstructions, and 5) testing of new concepts for safe increases in airport capacity for both single and multiple runway configurations. Fast-time models are based on a number of assumptions which limit their accuracy and range of application. For example, they do not predict the meandering vortex paths that can result from the wake interaction with atmospheric turbulence. However, these models should account for the mean vortex position and may predict the level of uncertainty. This presentation describes a study in which comparative evaluations of fast-time wake transport and decay models were conducted using data sets from previous field programs taken at MEM (1995), DFW (1997, 1999, 2000), and Denver (2005, 2006) airports. These data sets contain, meteorological data, lidar measurements of wake vortices, and the aircraft parameters associated with these vortices.

Fred H Proctor↗
DOE OSTI2020

Coupled Gas-Exchange Model for C4 Leaves Comparing Stomatal Conductance Models

Plant simulation models are abstractions of plant physiological processes that are useful for investigating the responses of plants to changes in the environment. Because photosynthesis and transpiration are fundamental processes that drive plant growth and water relations, a leaf gas-exchange model that couples their interdependent relationship through stomatal control is a prerequisite for explanatory plant simulation models. Here, we present a coupled gas-exchange model for C4 leaves incorporating two widely used stomatal conductance submodels: Ball–Berry and Medlyn models. The output variables of the model includes steady-state values of CO2 assimilation rate, transpiration rate, stomatal conductance, leaf temperature, internal CO2 concentrations, and other leaf gas-exchange attributes in response to light, temperature, CO2, humidity, leaf nitrogen, and leaf water status. We test the model behavior and sensitivity, and discuss its applications and limitations. The model was implemented in Julia programming language using a novel modeling framework. Our testing and analyses indicate that the model behavior is reasonably sensitive and reliable in a wide range of environmental conditions. The behavior of the two model variants differing in stomatal conductance submodels deviated substantially from each other in low humidity conditions. The model was capable of replicating the behavior of transgenic C4 leaves under moderate temperatures as found in the literature. The coupled model, however, underestimated stomatal conductance in very high temperatures. This is likely an inherent limitation of the coupling approaches using Ball–Berry type models in which photosynthesis and stomatal conductance are recursively linked as an input of the other.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

Use of Design of Experiments and Rule-Based Inference in Determining Neural Network Architectures for Loss of Control Detection

In this work, we describe methods for selecting the neural network architectures and input spaces to implement belief state inference on generic commercial transport aircraft. First, we highlight a case study on the planning, execution, and analysis of a set of experiments to determine the configurations of a conditional variational autoencoder (CVAE). We present a structured method that can be used in a number of aerospace applications, to optimize the structure and training parameters of the CVAE for belief state inference, using Design of Experiments (DOE) statistical methodologies. The motivation for this specific DOE was to identify the appropriate hyperparameters for measuring the CVAE reconstruction probability and latent space, such that the measurements can be used to infer qualitative state changes for the aircraft. We demonstrate that this process yields information about a trained neural network’s utility for this specific application, along with a quantifiable range of certainty. We execute 84 experiments using loss-of-control flight maneuver data from a NASA T-2 aircraft, demonstrating that this empirical process allows us to construct cheap and simple models with specific attributes amenable to belief state inference in aerospace applications. While theoretically, we could create a single CVAE with an input space the size of all measurable flight variables and environmental dynamics, it becomes intractable to use such a neural network in an in-situ intelligent multi-agent system. Using the recommendations from our case study, we introduce a technical approach for feasibly describing the belief space by (1) identifying significant statistical relationships among flight variables using rule induction, (2) using a set of rules that cover all features to define the input space of multiple CVAEs, and (3) forming a belief space based on the joint probability density of their collective latent spaces. This results in a series of relatively small matrix multiplications that can be performed in real time, as opposed to large matrix computations in a single CVAE. We demonstrate the application of this approach on the T-2 flight loss-of control experiments, using the architecture and hyperparameter recommendations from the case study. We compare the utilities of an individual CVAE trained on all flight variables and multiple CVAEs defined on subsets of flight variables for detecting qualitative changes in flight. We demonstrate that the use of multiple CVAEs with smaller input spaces permits the CVAE to capture more granular relationships in the latent space, permitting better state space characterization and loss-of-control detection.

Design of experiments↗