Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “history variables”

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 631 records · Page 35

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗

Shield volcanism and lithospheric structure beneath the Tharsis plateau, Mars

The heights of four great shield volcanoes, when interpreted as reflecting the local hydrostatic head on a common source of upwelling magma, provide significant constraints on models of lithospheric structure beneath the Tharsis plateau. If Bouguer gravity anomalies are modeled in terms of a variable thickness crust, and a two-component (crust/mantle) earth-like structure is assumed for the Martian lithosphere, the derived model lithosphere beneath the Tharsis plateau has the following properties: (1) the upper low-density 'crustal' component is thickened beneath the Tharsis plateau; (2) the lower high-density 'mantle' component is thinned beneath the Tharsis plateau; and (3) there is a net gradient on the base of the Martian lithosphere directed downward away from beneath the summit of the Tharsis plateau. A long history of magmatic intrusion is hypothesized to have been the cause of the updoming of the Tharsis plateau and the maintenance of the plateau in a state of only partial compensation.

Blasius, K. R.↗

A quantitative measure of the structure of gamma-ray burst time profiles

A cursory examination of cosmic gamma-ray burst time profiles indicates an inhomogeneous distribution of structure. In the first approximation, there seem to be two types of profiles; smooth ones with little structure and highly variable ones with lots of structure. To put this observation to the test, we have examined the statistical nature of the profile derivative to choose which parameter might best be called the burst 'spikiness'. We have found that a good estimator is given by a count of the number of 'spikes' (defined by a specific numerical recipe) and not by the rms deviations from either a pre-burst background or any type of moving average background. The application of this parameter to 30 burst time histories shows it to be consistent over a wide range of profile types. The analysis also reveals a preferred average time between spikes of approximately 1.5 seconds.

Lestrade, John P.↗

A New Hybrid Spatio-temporal Model for Estimating Daily Multi-year PM2.5 Concentrations Across Northeastern USA Using High Resolution Aerosol Optical Depth Data

The use of satellite-based aerosol optical depth (AOD) to estimate fine particulate matter PM(sub 2.5) for epidemiology studies has increased substantially over the past few years. These recent studies often report moderate predictive power, which can generate downward bias in effect estimates. In addition, AOD measurements have only moderate spatial resolution, and have substantial missing data. We make use of recent advances in MODIS satellite data processing algorithms (Multi-Angle Implementation of Atmospheric Correction (MAIAC), which allow us to use 1 km (versus currently available 10 km) resolution AOD data.We developed and cross validated models to predict daily PM(sub 2.5) at a 1X 1 km resolution across the northeastern USA (New England, New York and New Jersey) for the years 2003-2011, allowing us to better differentiate daily and long term exposure between urban, suburban, and rural areas. Additionally, we developed an approach that allows us to generate daily high-resolution 200 m localized predictions representing deviations from the area 1 X 1 km grid predictions. We used mixed models regressing PM(sub 2.5) measurements against day-specific random intercepts, and fixed and random AOD and temperature slopes. We then use generalized additive mixed models with spatial smoothing to generate grid cell predictions when AOD was missing. Finally, to get 200 m localized predictions, we regressed the residuals from the final model for each monitor against the local spatial and temporal variables at each monitoring site. Our model performance was excellent (mean out-of-sample R(sup 2) = 0.88). The spatial and temporal components of the out-of-sample results also presented very good fits to the withheld data (R(sup 2) = 0.87, R(sup)2 = 0.87). In addition, our results revealed very little bias in the predicted concentrations (Slope of predictions versus withheld observations = 0.99). Our daily model results show high predictive accuracy at high spatial resolutions and will be useful in reconstructing exposure histories for epidemiological studies across this region.

Air pollution↗

Ar-40/Ar-39 age of the Shergotty achondrite and implications for its post-shock thermal history

Ar-40/Ar-39 measurements are used to determine the age of the Shergotty achondrite and the chronology of the shock event responsible for the complete conversion of its plagioclase to maskelynite is discussed. Apparent ages are found to vary between 240 and 640 million years for the whole rock sample, with a plateau age of 254 million years for a maskelynite separate. The Rb-Sr age of 165 million years determined by Nyquist at al (1978) suggests that the maskelynite as well as the whole rock was incompletely degassed. Argon diffusion characteristics indicate a post-shock cooling time greater than 1000 years and a burial depth greater than 300 m for a thermal model of a cooling ejecta blanket of variable thickness. It is concluded that the shock event which degassed the argon and reset the Rb-Sr systematics occurred between 165 and 250 million years ago when the parent body experienced a collision in the asteroid belt.

Bogard, D. D.↗

Node Distortions in UiO-66 Inform Negative Thermal Expansion Mechanisms: Kinetic Effects, Frustration, and Lattice Hysteresis

In metal–organic frameworks (MOFs) the interplay between the dynamics of individual components and how these are constrained by the extended lattice can yield unusual emergent phenomena. For the archetypal Zr-MOF, UiO-66, we explore the cooperative dynamics of a Zr-node transformation that gives rise to negative thermal expansion (NTE). Here, using in situ synchrotron X-ray scattering, with powder diffraction and pair distribution function (PDF) analyses, we identify lattice hysteresis and a thermal ramp-rate-dependence of the thermal expansion. Specifically, kinetic trapping of distorted node states formed at high temperature, leads to broad variability in the apparent thermal expansion which ranges from large positive to large negative thermal expansion with coefficients of thermal expansion (CTE) from +45 to –80 × 10 –6 K –1 . Time-resolved relaxation studies at selected temperatures suggest that when equilibrated UiO-66 is intrinsically NTE, with a CTE of –35 × 10 –6 K –1 . Kinetic trapping of the node-distorted state following high temperature activation has broad implications for characterization and applications of these Zr-MOFs; the nonequilibrium node state depends on the thermal history of the sample with quench vs slow cooling likely to impact gas binding, pore volume, and accessible catalytic sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimentation in machine discovery

KEKADA, a system that is capable of carrying out a complex series of experiments on problems from the history of science, is described. The system incorporates a set of experimentation strategies that were extracted from the traces of the scientists' behavior. It focuses on surprises to constrain its search, and uses its strategies to generate hypotheses and to carry out experiments. Some strategies are domain independent, whereas others incorporate knowledge of a specific domain. The domain independent strategies include magnification, determining scope, divide and conquer, factor analysis, and relating different anomalous phenomena. KEKADA represents an experiment as a set of independent and dependent entities, with apparatus variables and a goal. It represents a theory either as a sequence of processes or as abstract hypotheses. KEKADA's response is described to a particular problem in biochemistry. On this and other problems, the system is capable of carrying out a complex series of experiments to refine domain theories. Analysis of the system and its behavior on a number of different problems has established its generality, but it has also revealed the reasons why the system would not be a good experimental scientist.

Kulkarni, Deepak↗

Computational Study of Flow Establishment in Hypersonic Pulse Facilities

This paper presents a study of the temporal evolution of the combustion flowfield established by the interaction of ram-accelerator-type projectiles with an explosive gas mixture accelerated to hypersonic speeds in an expansion tube. The Navier-Stokes equations for a chemically reacting gas are solved in a fully coupled manner using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order, total variation diminishing (TVD) scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of expansion tube flowfields are then presented for two different configurations. In particular, the development of the shock-induced combustion process is followed. In one case, designed to ensure ignition only in the boundary layer, the lateral extent of the combustion front during the initial transient phase was surprisingly large. The time histories of the calculated thrust and drag forces on the ram accelerator projectile are also presented.

Yungster, S.↗

Temporal Study 2021-2022: Sample-Based Surface Water Chemistry and Organic Matter Characterization across Watersheds in the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of temporal variability in sediment respiration rates in the Yakima River Basin. The dataset provides geochemistry and organic matter characterization data generated from samples collected at weekly or bi-weekly intervals at six sites across the Yakima River Basin in Washington, USA. Related sensor data will be published separately and can be used to link sediment respiration rates to biogeochemical processing rates. The data package was originally published in November 2022. It was updated in November 2023 (v2; modified files) and April 2025 (new and modified files). See the change history section in readme for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) dissolved inorganic carbon (DIC), dissolved organic carbon (DOC; reported as non-purgeable organic carbon; NPOC), total nitrogen (TN), total suspended solids (TSS), and ions; (5) averaged values from water chemistry data; (6) surface water sampling protocol; (7) sensor protocol; (8) readme; (9) methods codes; (10) international generic sample number (IGSN) mapping file; and (11) folder of high resolution characterization of organic via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). This folder contains two subfolders, one containing the .xml data files and the other containing instructions for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .R, .ref, or .xml. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected our data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

Development of a New Data Tool for Computing Launch and Landing Availability with Respect to Surface Weather

The Marshall Space Flight Center Natural Environments Branch has a long history of expertise in the modeling and computation of statistical launch availabilities with respect to weather conditions. Their existing data analysis product, the Atmospheric Parametric Risk Assessment (APRA) tool, computes launch availability given an input set of vehicle hardware and/or operational weather constraints by calculating the climatological probability of exceeding the specified constraint limits, APRA has been used extensively to provide the Space Shuttle program the ability to estimate impacts that various proposed design modifications would have to overall launch availability. The model accounts for both seasonal and diurnal variability at a single geographic location and provides output probabilities for a single arbitrary launch attempt. Recently, the Shuttle program has shown interest in having additional capabilities added to the APRA model, including analysis of humidity parameters, inclusion of landing site weather to produce landing availability, and concurrent analysis of multiple sites, to assist in operational landing site selection. In addition, the Constellation program has also expressed interest in the APRA tool, and has requested several additional capabilities to address some Constellation-specific issues, both in the specification and verification of design requirements and in the development of operations concepts. The combined scope of the requested capability enhancements suggests an evolution of the model beyond a simple revision process. Development has begun for a new data analysis tool that will satisfy the requests of both programs. This new tool, Probabilities of Atmospheric Conditions and Environmental Risk (PACER), will provide greater flexibility and significantly enhanced functionality compared to the currently existing tool.

Burns, K. Lee↗

Unconventional Wells Interference: Supervised Machine Learning for Detecting Fracture Hits

The primary objective of the study was development of a machine learning (ML)-based workflow for fracture hit (“frac hit”) detection and monitoring using shale oil-field data such as drilling surveys, production history (oil and produced water), pressure, and fracking start time and duration records. The ML method takes advantage of long short-term memory (LSTM) and multilayer perceptron (MLP) neural networks to identify the frac hits due to hydraulic communication between the fracking child well(s) and the producing parent well(s) within the same pad (intra-pad interaction) and/or on different pads (inter-pad interaction). It utilizes time series of pressure and production data from within a pad and from adjacent pads. The workflow can capture time variable features of frac hits when the model architecture is deep and wide enough, with enough trainable parameters for deep learning and feature extraction, as demonstrated in this paper by using training and testing subsets of the field data from selected neighboring pads with over a couple of hundred wells. The study was focused on frac-hit interaction among paired wells and demonstrated that the ML model, once trained, can predict the frac-hit probability.

58 GEOSCIENCES↗

Sampling problems: The small scale structure of precipitation

The quantitative measurement of precipitation characteristics for any area on the surface of the Earth is not an easy task. Precipitation is rather variable in both space and time, and the distribution of surface rainfall data given location typically is substantially skewed. There are a number of precipitation process at work in the atmosphere, and few of them are well understood. The formal theory on sampling and estimating precipitation appears considerably deficient. Little systematic attention is given to nonsampling errors that always arise in utilizing any measurement system. Although the precipitation measurement problem is an old one, it continues to be one that is in need of systematic and careful attention. A brief history of the presently competing measurement technologies should aid us in understanding the problem inherent in this measurement task.

Crane, R. K.↗

Progressive Damage Analysis of Laminated Composite (PDALC)-A Computational Model Implemented in the NASA COMET Finite Element Code

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete list of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occur during the load history. Residual strength predictions made with this information compared favorably with experimental measurements.

Lo, David C.↗

Progressive Damage Analysis of Laminated Composite (PDALC) (A Computational Model Implemented in the NASA COMET Finite Element Code)

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged damage variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete listing of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occurs during the loading history. Residual strength predictions made with this information compared favorably with experimental measurements.

Coats, Timothy W.↗

Kinetics Modeling for Design of Continuous Enzymatic Hydrolysis

Enzymatic hydrolysis of cellulose to monomeric sugars continues to be a limiting step in cost-effectively producing sugar and fermentation-based biofuels from biomass. In particular, the high cost of enzymes coupled with the long time-scale of reaction pose challenges to economic viability. Performing enzymatic hydrolysis in a continuous mode with enzyme recycle may provide a path towards substantially reduced costs for sugar production, but design and analysis are complicated by a lack of suitable kinetics models. Computationally attainable models, such as fractal-based models, require knowledge of the reaction-history of the biomass, and are thus only suitable for describing the batch reactions from which they are derived. Fundamental models, while potentially more generalizable, are often too computationally intensive to use in process or reactor modeling. In this work, a phenomenological rate model is proposed based on a two-phase substrate representation. Good agreement is seen between batch and continuous enzymatic hydrolysis (CEH) experiment data, which validates the model and enables us to solve for reactor design parameters, such as CEH reactor size and stream flow rates, based on process variables like yield. This model is integrated with techno-economic analysis software to explore economic sensitivities. Important design optimizations and tradeoffs are identified and quantified, including the relative cost imposed by rate slowdown from sugar inhibition versus the cost to remove and concentrate sugars at a lower concentration. It also identifies, high-leverage avenues for further exploration, such as increasing the maximum feasible solids concentration, and sustaining high membrane flux and reliability.

09 BIOMASS FUELS↗

Using AVIRIS images to measure temporal trends in abundance of photosynthetic and nonphotosynthetic canopy components

The Jasper Ridge Biological Preserve, Stanford, California is a good example of hardwood rangeland ecosystems in California. Structurally, it is composed of a mosaic of serpentine grasslands, oak savannah, coastal chaparral, and mixed evergreen woodland, representing a broad cross-section of physiognomic classes. The Mediterranean climate produces an extended seasonal drought lasting throughout most of the growing season and has significant impact on the expression of divergent phenological patterns related to contrasting ecological strategies of these taxa. The region is well understood biologically due to the rich history of ecological research at the site. Thus, community characteristics, physiological characteristics, phenology, and temporal dynamics are reasonably well understood for many of the dominant species. Because of its proximity to NASA Ames Research Center, it has been subject to a large number of aircraft data acquisitions over many years. A more complete examination of this database would provide an opportunity to test current remote sensing hypotheses for measurement and detection of ecological attributes, particularly those involving canopy chemistry and physiology. Better definition of ecological rules might permit development of remotely sensed surrogate variables for biological properties that cannot be directly measured or measured with sufficient accuracy.

Ustin, Susan L.↗

S-process krypton of variable isotopic composition in the Murchison meteorite

Data are reported which, for the first time, permit the derivation of the full isotopic spectrum of s-process krypton with reasonable precision. It is shown that this s-Kr in a residue from the Murchison meteorite did not originate in one single s-process but rather is a mixture of contributions from stellar environments where the density of free neutrons was not the same. The astrophysical conditions under which this krypton was produced were distinct from those that have been invoked to explain the solar system s-process abundance. Similar to the C-13-rich carbon component in an aliquot of the same residue, the s-process Kr from different astrophysical sites has retained its identity during the accumulation and subsequent history of the meteorite.

Ott, Urlich↗

A Generic Nonlinear Aerodynamic Model for Aircraft

A generic model of the aerodynamic coefficients was developed using wind tunnel databases for eight different aircraft and multivariate orthogonal functions. For each database and each coefficient, models were determined using polynomials expanded about the state and control variables, and an othgonalization procedure. A predicted squared-error criterion was used to automatically select the model terms. Modeling terms picked in at least half of the analyses, which totalled 45 terms, were retained to form the generic nonlinear aerodynamic (GNA) model. Least squares was then used to estimate the model parameters and associated uncertainty that best fit the GNA model to each database. Nonlinear flight simulations were used to demonstrate that the GNA model produces accurate trim solutions, local behavior (modal frequencies and damping ratios), and global dynamic behavior (91% accurate state histories and 80% accurate aerodynamic coefficient histories) under large-amplitude excitation. This compact aerodynamics model can be used to decrease on-board memory storage requirements, quickly change conceptual aircraft models, provide smooth analytical functions for control and optimization applications, and facilitate real-time parametric system identification.

Grauer, Jared A.↗