Convergent and structure preserving architectures for SciML.
Abstract not provided.
Trask, Nathaniel Albert↗
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.
Structure preserving machine learning for data-driven multiscale/multiphysics modeling.
Abstract not provided.
Trask, Nathaniel↗
Structure-preserving model reduction for marginally stable LTI systems.
Abstract not provided.
Peng, Liqian↗
State Space Reconstruction from Embeddings of Partial Observables in Structural Dynamic Systems for Structure-Preserving Data-Driven Methods.
Abstract not provided.
Najera-Flores, David↗
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.
97 MATHEMATICS AND COMPUTING↗
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.
77 NANOSCIENCE AND NANOTECHNOLOGY↗
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".
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Brilliant angle-independent structural colours preserved in weevil scales from the Swiss Pleistocene
Extant weevils exhibit a remarkable colour palette that ranges from muted monochromatic tones to rainbow-like iridescence, with the most vibrant colours produced by three-dimensional photonic nanostructures housed within cuticular scales. Although the optical properties of these nanostructures are well understood, their evolutionary history is not fully resolved, in part due to a poor knowledge of their fossil record. Here, we report three-dimensional photonic nanostructures preserved in brightly coloured scales of two weevils, belonging to the genus Phyllobius or Polydrusus , from the Pleistocene (16–10 ka) of Switzerland. The scales display vibrant blue, green and yellow hues that resemble those of extant Phyllobius/Polydrusus . Scanning electron microscopy and small-angle X-ray scattering analyses reveal that the subfossil scales possess a single-diamond photonic crystal nanostructure. In extant Phyllobius/Polydrusus , the near-angle-independent blue and green hues function primarily in crypsis. The preservation of far-field, angle-independent structural colours in the Swiss subfossil weevils and their likely function in substrate matching confirm the importance of investigating fossil and subfossil photonic nanostructures to understand the evolutionary origins and diversification of colours and associated behaviours (e.g. crypsis) in insects.
54 ENVIRONMENTAL SCIENCES↗
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.
Mathematics↗
Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm
Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm
Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.
Zhang, Xin↗
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.
Neurosciences & Neurology↗
Fast neural Poincaré maps for toroidal magnetic fields
Poincaré maps for toroidal magnetic fields are routinely employed to study gross confinement properties in devices built to contain hot plasmas. In most practical applications, evaluating a Poincaré map requires numerical integration of a magnetic field line, a process that can be slow and that cannot be easily accelerated using parallel computations. We propose a novel neural network architecture, the HénonNet, and show that it is capable of accurately learning realistic Poincaré maps from observations of a conventional field-line-following algorithm. After training, such learned Poincaré maps evaluate much faster than the field-line integration method. Moreover, the HénonNet architecture exactly reproduces the primary physics constraint imposed on field-line Poincaré maps: flux preservation. Furthermore, this structure-preserving property is the consequence of each layer in a HénonNet being a symplectic map. We demonstrate empirically that a HénonNet can learn to mock the confinement properties of a large magnetic island by using coiled hyperbolic invariant manifolds to produce a sticky chaotic region at the desired island location. This suggests a novel approach to designing magnetic fields with good confinement properties that may be more flexible than ensuring confinement using KAM tori.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Topology-Preserving Nonlinear Structured Light Generation from All-Dielectric Metasurfaces
We experimentally demonstrate the simultaneous generation of phase-locked structured light beams, including optical vortices, Hopf links, and trefoil knots at fundamental and tripled frequencies in all-dielectric nonlinear optical metasurfaces.
Gao, Jiannan↗
Discrete gravity with local Lorentz invariance
A novel structure-preserving algorithm for general relativity in vacuum is derived from a lattice gauge theoretic discretization of the tetradic Palatini action. Here, the resulting model of discrete gravity is demonstrated to preserve local Lorentz invariance and symplectic structure.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Canonical and noncanonical Hamiltonian operator inference
Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.
97 MATHEMATICS AND COMPUTING↗
Numerical integration of stochastic contact Hamiltonian systems via stochastic Herglotz variational principle
Within this work, we establish a stochastic contact variational integrator and its discrete version via stochastic Herglotz variational principle for stochastic contact Hamiltonian systems. A general structure-preserving stochastic contact method is provided to seek the stochastic contact variational integrators. Numerical experiments are performed to verify the validity of this approach.
97 MATHEMATICS AND COMPUTING↗