Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “learning dynamics”

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 217 records · Page 12

Unstructured mesh generation and adaptivity

An overview of current unstructured mesh generation and adaptivity techniques is given. Basic building blocks taken from the field of computational geometry are first described. Various practical mesh generation techniques based on these algorithms are then constructed and illustrated with examples. Issues of adaptive meshing and stretched mesh generation for anisotropic problems are treated in subsequent sections. The presentation is organized in an education manner, for readers familiar with computational fluid dynamics, wishing to learn more about current unstructured mesh techniques.

Mavriplis, D. J.↗

Multigrid techniques for unstructured meshes

An overview of current multigrid techniques for unstructured meshes is given. The basic principles of the multigrid approach are first outlined. Application of these principles to unstructured mesh problems is then described, illustrating various different approaches, and giving examples of practical applications. Advanced multigrid topics, such as the use of algebraic multigrid methods, and the combination of multigrid techniques with adaptive meshing strategies are dealt with in subsequent sections. These represent current areas of research, and the unresolved issues are discussed. The presentation is organized in an educational manner, for readers familiar with computational fluid dynamics, wishing to learn more about current unstructured mesh techniques.

Mavriplis, D. J.↗

Measurement of Meteor Impact Experiments Using Three-Component Particle Image Velocimetry

The study of hypervelocity impacts has been aggressively pursued for more than 30 years at Ames as a way to simulate meteoritic impacts. Development of experimental methods coupled with new perspectives over this time has greatly improved the understanding of the basic physics and phenomenology of the impact process. These fundamental discoveries have led to novel methods for identifying impact craters and features in craters on both Earth and other planetary bodies. Work done at the Ames Vertical Gun Range led to the description of the mechanics of the Chicxualub crater (a.k.a. K-T crater) on the Yucatan Peninsula, widely considered to be the smoking gun impact that brought an end to the dinosaur era. This is the first attempt in the world to apply three-component particle image velocimetry (3-D PIV) to measure the trajectory of the entire ejecta curtain simultaneously with the fluid structure resulting from impact dynamics. The science learned in these experiments will build understanding in the entire impact process by simultaneously measuring both ejecta and atmospheric mechanics.

Heineck, James T.↗

Southern California Megacity CO2, CH4, and CO Flux Estimates Using Ground- and Space-Based Remote Sensing and a Lagrangian Model

We estimate the overall CO2, CH4, and CO flux from the South Coast Air Basin using an inversion that couples Total Carbon Column Observing Network (TCCON) and Orbiting Carbon Observatory-2 (OCO-2) observations, with the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model and the Open-source Data Inventory for Anthropogenic CO2 (ODIAC). Using TCCON data we estimate the direct net CO2 flux from the So-CAB to be 104±26 Tg CO2 yr(exp -1) for the study period of July 2013–August 2016. We obtain a slightly higher estimate of 120±30 Tg CO2 yr(exp -1) using OCO-2 data. These CO2 emission estimates are on the low end of previous work. Our net CH4 (360±90 Gg CH4 y(exp -1)) flux estimate is in agreement with central values from previous top-down studies going back to 2010 (342–440 Gg CH4 yr(exp -1)). CO emissions are estimated at 487±122 Gg CO yr(exp -1), much lower than previous top-down estimates (1440 Gg CO yr(exp -1)). Given the decreasing emissions of CO, this finding is not unexpected. We perform sensitivity tests to estimate how much errors in the prior, errors in the covariance, different inversion schemes, or a coarser dynamical model influence the emission estimates. Overall, the uncertainty is estimated to be 25%, with the largest contribution from the dynamical model. Lessons learned here may help in future inversions of satellite data over urban areas.

Total Carbon Column Observing Network (TCCON)↗

Noble Gases in Samples Returned From Asteroids Ryugu and Bennu

Introduction: The asteroid sample return missions, JAXA’s Hayabusa2 [1] and NASA’s OSIRIS-REx [2], fromC-rich asteroids Ryugu and Bennu offer unique opportunities to sample pristine early solar system material, unaffected by terrestrial weathering and contamination, and with a known geological context. The comparison with most primi-tive members of carbonaceous –and inner solar system ordinary and enstatite– chondrite classes allows us to (i) assess the variation (or the lack of thereof) of the distribution of the most volatile elements throughout the entire early solar system, (ii) determine the effects of parent body thermal and aqueous alteration, and (iii) learn about the dynamic evolution and history of the various parent bodies, as sampled by both return missions and meteorite delivery. Here we will present noble gas data obtained from newly received particles from Ryugu and Bennu. Samples & Analysis: We received nine particles allocated by JAXA, with masses of 0.5-1.7 mg mass: four col-lected during Hayabusa2’s first touchdown (chamber A) and five from a small artificial crater (chamber C). We received three particles from Bennu of 0.07, 0.89 and 0.94 mg mass (OREX-800032-102/3/4). All particles were processed, allocated, weighed, and transferred into our UHV sample chamber within pure N2, to prevent any contact with air. We examined all He-Xe isotopes with our custom-built mass spectrometer “Albatros” ([3,4]). Gases were extracted individually from each particle by melting induced by a 1064 μm Nd:YAG IR laser in 2–4 heating steps. Large amounts of water and other reactive volatiles were released, consistent with the presence of abundant volatiles from both Ryugu and Bennu. Residues of Ryugu particles remaining after lasering, subsequently exposed to air, were additionally extracted in a crucible to verify complete gas extraction. Only the particle that still showed a “crystalline” structure after lasering contained detectable Ryugu noble gases, whereas all other particles turned immediately into glassy spheres by laser heating and lacked any gas. Results & Discussion: Ryugu and Bennu matter contains noble gases similar to the most pristine, unheated aqueously altered meteorites, e.g., CM or CI chondrites [5,6]. Remarkably, all four chamber A samples that originated from Ryugu’s surface show solar wind (SW) Ne (trapped (20Ne/22Ne)tr ~12–13), whereas all five chamber C samples do not ((20Ne/22Ne)tr ~9.0–10.7 consistent with the presence of non-solar trapped Ne components such as Q and HL). This demonstrates that JAXA’s sampling strategy succeeded: Chamber C samples analyzed in this work originated from areas shielded from SW, slightly below the surface, that was sampled after the impactor excavated material. The three Bennu particles populate the same range as the Ryugu particles in Ne isotope space with one particle containing SW suggesting that it was part of the uppermost surface layer of Bennu during sampling, while Ne in the smallest sample can be explained best by mixing Ne from presolar diamonds and SiC or graphite, and that in the third one by mixing of, e.g., Q and HL. Using cosmic ray production rates given by [7] yields preliminary exposure ages for Ryugu particles in the same range as determined in [7]. There is, as expected, also no systematic difference between samples from chambers A and C because both sample sets were likely collected within the uppermost layer of Ryugu, which were not shielded from cosmic rays. The concentration of cosmogenic (cos) Ne in the Bennu particle, where Necos was resolvable, is in the same range as found for Ryugu, suggesting roughly similar exposure time to cosmic rays for 5–8 Ma for all particles examined here. The Xe isotopic compositions of all particles are consistent with those of phase Q, with small additions of Xe-HL from presolar diamonds and excess 129Xe* from short-lived 129I decay. The comparably high 129Xe*/129XeQ ratio is similar in all samples, suggesting similar initial I concentrations and closure to Xe loss –or the incorporation of a well-mixed Xe reservoir. The Bennu material was in brief contact with air in the capsule, after atmospheric entry and before storage in pure N2 in a clean room. Both, the –most diagnostic– Xe isotopes and the Ar/Xe vs. Kr/Xe systematics prove that the noble gases are not affected by terrestrial contamination, in contrast to many CI, CY, and CM chondrites found on Earth [5,6], illustrating the importance of sample return. He-Xe concentrations in all particles unaffected bySW are comparable to but at the upper end of ranges reported [5-7]. In accordance with the strong aqueous alteration experienced on the parent bodies of Ryugu and Bennu [1,2], all samples lack the Ar-rich and water-susceptible noble gas components found in less aqueously altered CMs or CRs [e.g., 5,8] and all COs.

Henner Busemann↗

Frameworks, Algorithms, and Scalable Technologies for Mathematics (FASTMath) SciDAC Institute

As computational models scale to larger computers, the rate at which they produce data has far outstripped the same computers ability to write that data and further the file systems ability to store that data. Almost all of the SciDAC applications, but especially those related to fusion solve very large scale PDEs whose scientific output his impacted by this problem. To gain access to dynamics in an exascale simulation that are not identifiable a priori and to make that dynamical data available to machine learning requires fundamental research in the area of in situ data data analytics. Here data analytics includes compression, visualization, uncertainty quantification, and machine learning. This in situ data analytics will enable on-the-fly spatial and temporal compression of solution dynamics, expose that space-time compressed field to machine learning algorithms that have been specialized to work with dynamically evolving data (existing machine learning algorithms treat data sets as static), greatly improving the opportunity for machine learning to provide feedback to the compression, all within an ongoing simulation, without the need to write data to files. The same concepts are also being applied to uncertainty quantification and multi-fidelity modeling which have similar needs for spatial and temporal compression of the ongoing exascale simulation to perform either without the typical, unacceptable writing of data to files.

97 MATHEMATICS AND COMPUTING↗

Mars Atmospheric Dynamics

The Martian atmosphere is dynamically similar to the Earth's. Its spin-axis rotation rate is only minutes longer than Earth's so the Coriolois force is nearly identical to Earth's. The inclination of its spin axis is also similar to Earth's giving it similarity in seasonal change. And the Martian atmosphere is nearly transparent to solar radiation (except during dust periods) such that it is heated primarily by upwelling infrared radiation from the surface. These characteristics make Mars an ideal laboratory for studying the dynamics of rapidly rotating differentially heated atmospheres. This talk reviews what we have learned about Mars atmospheric dynamics and how if compares with Earth. The source of information to make such a comparison comes from observations and models. The former are sparse and that the latter have played a major role in shaping our thinking about the general circulation on Mars. However, the models need validation. Fortunately, the first two orbiters in NASA's Mars Surveyor Program have instrumentation to address many of the issues related to the general circulation and climate of Mars. The first, Mars Global Surveyor, is already at Mars gathering data. The second, the Mars 98 Orbiter to be launched later this year, carries a dedicated atmospheric sounder. Thus, much will be learned about Mars' atmosphere in the next few years.

Haberle, Robert↗

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗

Machine Learning Enabled Advanced Electron Tomography for Resolving Chemical Inhomogeneity and Materials Dynamics in Lithium-Ion Battery Electrodes

The objective of this project is to develop machine learning-assisted electron microscopy, together with three-dimensional, cryogenic, and in-situ imaging techniques, for resolving chemical inhomogeneity and materials dynamics in lithium battery electrodes and interfaces. The project aims to expand the spatial, temporal, and dimensional resolution of transmission electron microscopy and to enable quantitative analysis of beam-sensitive battery materials.

25 ENERGY STORAGE↗

Molecular simulation using transfer-learned potentials for the disordered nanoscale structure of nitrogen-doped nanoporous carbons

Machine learning (ML)-based molecular dynamics (MD) simulations of the formation of a class of N-doped nanoporous carbons are performed to assess their disordered partially graphitized nanoscale structure. The study is motivated by the effectiveness of so-called nitrogen assembly carbons (NACs) for catalysis applications. Benchmark simulations for pure-C disordered graphitic systems reveal the importance of reliably capturing the vdW component of the potentials in order to accurately describe the tendency for layering of disordered graphene-like sheets. In our modeling, this is achieved by a transfer learning strategy incorporating features of the energetics from the optB88-vdW DFT functional into potentials initially trained with a less expensive functional, thereby providing a superior description of the pure-C systems. Generation from MD simulations of realistic partially graphitized structures is significantly more challenging for N-doped versus for pure C systems. However, such structures are achieved by a tailored MD simulation protocol mimicking the experimental synthesis process and in particular incorporating an annealing and subsequent quenching stages. Simulated PXRD patterns effectively reproduce the features of experimental observations for NACs, including the appearance of a prominent but broad (002) peak at around 25, and the development of another weaker feature associated with in-layer ordering of mixed C-N graphene-like sheets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of structural dynamic data from Skylab. Volume 1: Technical discussion

A compendium of Skylab structural dynamics analytical and test programs is presented. These programs are assessed to identify lessons learned from the structural dynamic prediction effort and to provide guidelines for future analysts and program managers of complex spacecraft systems. It is a synopsis of the structural dynamic effort performed under the Skylab Integration contract and specifically covers the development, utilization, and correlation of Skylab Dynamic Orbital Models.

Demchak, L.↗

Dynamic Range and Sensitivity Requirements of Satellite Ocean Color Sensors: Learning from the Past

Sensor design and mission planning for satellite ocean color measurements requires careful consideration of the signal dynamic range and sensitivity (specifically here signal-to-noise ratio or SNR) so that small changes of ocean properties (e.g., surface chlorophyll-a concentrations or Chl) can be quantified while most measurements are not saturated. Past and current sensors used different signal levels, formats, and conventions to specify these critical parameters, making it difficult to make cross-sensor comparisons or to establish standards for future sensor design. The goal of this study is to quantify these parameters under uniform conditions for widely used past and current sensors in order to provide a reference for the design of future ocean color radiometers. Using measurements from the Moderate Resolution Imaging Spectroradiometer onboard the Aqua satellite (MODISA) under various solar zenith angles (SZAs), typical (L(sub typical)) and maximum (L(sub max)) at-sensor radiances from the visible to the shortwave IR were determined. The Ltypical values at an SZA of 45 deg were used as constraints to calculate SNRs of 10 multiband sensors at the same L(sub typical) radiance input and 2 hyperspectral sensors at a similar radiance input. The calculations were based on clear-water scenes with an objective method of selecting pixels with minimal cross-pixel variations to assure target homogeneity. Among the widely used ocean color sensors that have routine global coverage, MODISA ocean bands (1 km) showed 2-4 times higher SNRs than the Sea-viewing Wide Field-of-view Sensor (Sea-WiFS) (1 km) and comparable SNRs to the Medium Resolution Imaging Spectrometer (MERIS)-RR (reduced resolution, 1.2 km), leading to different levels of precision in the retrieved Chl data product. MERIS-FR (full resolution, 300 m) showed SNRs lower than MODISA and MERIS-RR with the gain in spatial resolution. SNRs of all MODISA ocean bands and SeaWiFS bands (except the SeaWiFS near-IR bands) exceeded those from prelaunch sensor specifications after adjusting the input radiance to L(sub typical). The tabulated L(sub typical), L(sub max), and SNRs of the various multiband and hyperspectral sensors under the same or similar radiance input provide references to compare sensor performance in product precision and to help design future missions such as the Geostationary Coastal and Air Pollution Events (GEO-CAPE) mission and the Pre-Aerosol-Clouds-Ecosystems (PACE) mission currently being planned by the U.S. National Aeronautics and Space Administration (NASA).

chlorophyll fluorescence↗

Active learning of a crystal plasticity flow rule from discrete dislocation dynamics simulations

Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

36 MATERIALS SCIENCE↗

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electron-phonon interaction on optical spectra of nanoelectronic devices

Information obtained on the solid-state lattice dynamics by electron-phonon interaction between lattice phonons and electrons could open up to learn more about lattice dynamics and to apply it in nanoelectronic devices including software reliability, nano-size capacitors, master clock sources, as well as non-contact temperature probes on nano-electronic and photonicdevices.

electron-phonon interactions lattice dynamics SrTi↗

Learning of networked spreading models from noisy and incomplete data

Recent years have seen a lot of progress in algorithms for learning parameters of spreading dynamics from both full and partial data. Some of the remaining challenges include model selection under the scenarios of unknown network structure, noisy data, missing observations in time, as well as an efficient incorporation of prior information to minimize the number of samples required for an accurate learning. Here, in this work, we introduce a universal learning method based on a scalable dynamic message-passing technique that addresses these challenges often encountered in real data. The algorithm leverages available prior knowledge on the model and on the data, and reconstructs both network structure and parameters of a spreading model. We show that a linear computational complexity of the method with the key model parameters makes the algorithm scalable to large network instances.

97 MATHEMATICS AND COMPUTING↗

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗