SEARCH · Engineering Papers
Results for “Structure preservation”
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.
Code Demonstration: Approximation of nearly-periodic symplectic maps via structure-preserving neural networks
In this notebook, we use the symplectic gyroceptron architecture from [Duruisseaux et al., 2022] to learn a surrogate map for the nearly-periodic symplectic flow map associated to a nearly-periodic Hamiltonian system composed of two nonlinearly coupled oscillators, where one of them oscillates significantly faster than the other: $\Big\lbrace$ $ {\dot{q}_1 = p_1 \hspace{5mm} \dot{p}_1 = -q_1 - \epsilon \partial_{q1}U(q_1 , q_2 )} \atop {\dot{q}_2 = \epsilon p_2 \hspace{4mm} \dot{p}_2 = -\epsilon q_2 - \epsilon \partial_{q2}U(q_1 , q_2 )}$. These equations of motion are the Hamilton's equations associated to the Hamiltonian $H_{\epsilon}(q_1,q_2,p_1,p_2)=\frac{1}{2}(q^2_1+p^2_1)+\frac{1}{2}\epsilon(q^2_2+p^2_2)+\epsilon$$U(q_1,q_2)$.
Structure-preserving, high-order and oscillation-limited transport operators.
Abstract not provided.
Data-Driven Whitney Forms for Structure-Preserving Control Volume Analysis.
Abstract not provided.
Discovery of Whitney Formsfor construction of structure preserving data-driven models.
Abstract not provided.
Machine-Learned Whitney Forms for Structure Preservation
Explore the source record for details and available documents.
Data-Driven Compatible Schemes for Structure Preserving Surrogate Models
Explore the source record for details and available documents.
W24-09 Geometric Mechanics Formulations and Structure Preserving Discretizations: An Introductory Course
Explore the source record for details and available documents.
Structure-preserving model reduction for marginally stable LTI systems.
Abstract not provided.
State Space Reconstruction from Embeddings of Partial Observables in Structural Dynamic Systems for Structure-Preserving Data-Driven Methods.
Abstract not provided.
An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks
The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.
Data-driven reduced-order models for port-Hamiltonian systems with operator inference
Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.
Phase-Field Modeling of Materials Interfaces and Nanostructures
Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.
Differential geometric approaches to momentum-based formulations for fluids [Slides]
This SAND report documents CIS Late Start LDRD Project 22-0311, "Differential geometric approaches to momentum-based formulations for fluids". The project primarily developed geometric mechanics formulations for momentum-based descriptions of nonrelativistic fluids, utilizing a differential geometry/exterior calculus treatment of momentum and a space+time splitting. Specifically, the full suite of geometric mechanics formulations (variational/Lagrangian, Lie-Poisson Hamiltonian and Curl-Form Hamiltonian) were developed in terms of exterior calculus using vector-bundle valued differential forms. This was done for a fairly general version of semi-direct product theory sufficient to cover a wide range of both neutral and charged fluid models, including compressible Euler, magnetohydrodynamics and Euler-Maxwell. As a secondary goal, this project also explored the connection between geometric mechanics formulations and the more traditional Godunov form (a hyperbolic system of conservation laws). Unfortunately, this stage did not produce anything particularly interesting, due to unforeseen technical difficulties. There are two publications related to this work currently in preparation, and this work will be presented at SIAM CSE 23, at which the PI is organizing a mini-symposium on geometric mechanics formulations and structure-preserving discretizations for fluids. The logical next step is to utilize the exterior calculus based understanding of momentum coupled with geometric mechanics formulations to develop (novel) structure-preserving discretizations of momentum. This is the main subject of a successful FY23 CIS LDRD "Structure-preserving discretizations for momentum-based formulations of fluids".
Microbial Preservation in Sulfates in the Haughton Impact Structure Suggests Target in Search for Life on Mars
Microbes in Haughton Crater Sulfates: Impact craters are of high interest in planetary exploration because they are viewed as possible sites for evidence of life [1]. Hydrothermal systems in craters are particularly regarded as sites where primitive life could evolve. Evidence from the Miocene Haughton impact structure shows that crater hydrothermal deposits may also be a preferred site for subsequent colonization and hence possible extant life: Hydrothermal sulfates at Haughton are colonized by viable cyanobacteria [2]. The Haughton impact structure, Devon Island, Canadian High Arctic, is a 24 km-diameter crater of mid-Tertiary age. The structure preserves an exceptional record of impact-induced hydrothermal activity, including sulfide, and sulfate mineralization [3]. The target rocks excavated at the site included massive gypsum-bearing carbonate rocks of Ordovician age. Impact-remobilized sulfates occur as metre-scale masses of intergrown crystals of the clear form of gypsum selenite in veins and cavity fillings within the crater s impact melt breccia deposits [4]. The selenite is part of the hydrothermal assemblage as it was precipitated by cooling hot waters that were circulating as a result of the impact. Remobilization of the sulfate continues to the present day, such that it occurs in soil crusts (Fig. 1) including sandy beds with a gypsum cement. The sulfate-cemented beds make an interesting comparison with the sulfate-bearing sandy beds encountered by the Opportunity MER [5]. The selenite crystals are up to 0.3 m in width, of high purity, and transparent. They locally exhibit frayed margins where cleavage surfaces have separated. This exfoliation may be a response to freeze-thaw weathering. The selenite contains traces of rock detritus, newly precipitated gypsum, and microbial colonies. The rock detritus consists of sediment particles which penetrated the opened cleavages by up to 2cm from the crystal margins. Some of the detritus is cemented into place by gypsum, which must have been dissolved and reprecipitated from the host selenite.
Early Life on Earth and the Search for Extraterrestrial Biosignatures
In the last 2 years, scientists within the ARES Directorate at JSC have applied the technology of Secondary Ion Mass Spectrometry (SIMS) to individual organic structures preserved in Archean (approximately 3 billion years old) sediments on Earth. These organic structures are among the oldest on Earth that may be microfossils - structurally preserved remnants of ancient microbes. The SIMS work was done to determine the microfossils' stable carbon isotopic composition (delta C-13 values). This is the first time that such ancient, potential microfossils have been successfully analyzed for their individual delta C-13 values. The results support the interpretation that these structures are remnants of early life on Earth and that they may represent planktonic organisms that were widely distributed in the Earth's earliest oceans. This study has been accepted for publication in the journal Geology.
Quadrature Based Neural Network Learning of Stochastic Hamiltonian Systems
Hamiltonian Neural Networks (HNNs) provide structure-preserving learning of Hamiltonian systems. In this paper, we extend HNNs to structure-preserving inversion of stochastic Hamiltonian systems (SHSs) from observational data. We propose the quadrature-based models according to the integral form of the SHSs’ solutions, where we denoise the loss-by-moment calculations of the solutions. The integral pattern of the models transforms the source of the essential learning error from the discrepancy between the modified Hamiltonian and the true Hamiltonian in the classical HNN models into that between the integrals and their quadrature approximations. This transforms the challenging task of deriving the relation between the modified and the true Hamiltonians from the (stochastic) Hamilton–Jacobi PDEs, into the one that only requires invoking results from the numerical quadrature theory. Meanwhile, denoising via moments calculations gives a simpler data fitting method than, e.g., via probability density fitting, which may imply better generalization ability in certain circumstances. Numerical experiments validate the proposed learning strategy on several concrete Hamiltonian systems. The experimental results show that both the learned Hamiltonian function and the predicted solution of our quadrature-based model are more accurate than that of the corrected symplectic HNN method on a harmonic oscillator, and the three-point Gaussian quadrature-based model produces higher accuracy in long-time prediction than the Kramers–Moyal method and the numerics-informed likelihood method on the stochastic Kubo oscillator as well as other two stochastic systems with non-polynomial Hamiltonian functions. Moreover, the Hamiltonian learning error εH arising from the Gaussian quadrature-based model is lower than that from Simpson’s quadrature-based model. These demonstrate the superiority of our approach in learning accuracy and long-time prediction ability compared to certain existing methods and exhibit its potential to improve learning accuracy via applying precise quadrature formulae.
Cognitive aging outcomes are related to both tau pathology and maintenance of cingulate cortex structure
Abstract INTRODUCTION Successful cognitive aging is related to both maintaining brain structure and avoiding Alzheimer's disease (AD) pathology, but how these factors interplay is unclear. METHODS A total of 109 cognitively normal older adults (70+ years old) underwent amyloid beta (Aβ) and tau positron emission tomography (PET) imaging, structural magnetic resonance imaging (MRI), and cognitive testing. Cognitive aging was quantified using the cognitive age gap (CAG), subtracting chronological age from predicted cognitive age. RESULTS Lower CAG (younger cognitive age) was related to slower decline in episodic memory, multi‐domain cognition, and atrophy of the midcingulate cortex (MCC). Lower entorhinal cortical tau was linked to slower decline in episodic memory, multi‐domain cognition, and hippocampal atrophy. DISCUSSION These results suggest that aging outcomes may be influenced by two independent pathways: one associated with tau accumulation, affecting primarily memory and hippocampal atrophy, and another involving tau‐independent structural preservation of the MCC, benefiting multi‐domain cognition over time. Highlights Younger cognitive age (lower cognitive age gap [CAG]) is related to slower cognitive decline. Lower CAG is linked to slower midcingulate cortex (MCC) atrophy. Reduced tau in the entorhinal cortex is related to less hippocampal atrophy and cognitive decline. Structural preservation of the MCC benefits multi‐domain cognition over time. Two independent pathways influence cognitive aging: tau accumulation and MCC preservation.