Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized linear model”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Emergence of local scaling relations in adsorption energies on high-entropy alloys

Alloying has been proposed to circumvent scaling relations between the adsorption energies thus allowing for the complete optimization of multistep reactions. Herein the fidelity of scaling rules on high-entropy alloy (HEA) surfaces is assessed focusing on hydrogen-containing molecules, *AH x for A = C and N (x = 0, 1, 2, 3), A = S (x = 0, 1, 2) and A = O (x = 0, 1). Using an adsorbate- and site-specific deep learning model to rapidly compute the adsorption energies on CoMoFeNiCu HEA surfaces, the energies of *AH x and *A are shown to be linearly correlated if *A and *AH x have identical adsorption site symmetry. However, a local linear dependence emerges between the configuration-averaged adsorption energies irrespective of the site symmetry. Although these correlations represent a weaker form of the scaling relationships, they are sufficient to prohibit the optimization of multistep reactions. The underpinning of this behavior is twofold (1) the nearsightedness principle and (2) the narrow distribution of the adsorption energies around the mean-field value. While the nearsightedness is general for all electronic systems, the second criterion applies in HEAs with relatively strong reactive elements. The present findings strongly suggest that alloys may not generally enable the breaking of scaling relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Novel Probes Project: Tests of gravity on astrophysical scales

The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists, is introduced. This review focuses on tests of gravity and is intended to be of use primarily to observers, as well as theorists with an interest in the development of experimental tests. It is twinned with a separate upcoming review on dark matter self-interactions. The review focuses on astrophysical tests of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. This regime is complementary to both strong-field tests of gravity and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified-gravity models lead to characteristic signatures, specifically on astrophysical scales. The potential of these probes is not limited by cosmic variance but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, and large-scale structure. The groundwork is laid for a thorough exploration of the weak-field, nonlinear regime, with an eye to using the current and next generation of observations for tests of gravity. The scene is set by showing how gravitational theories beyond general relativity are expected to behave, focusing primarily on screening mechanisms. Analytic and numerical techniques for exploring the relevant astrophysical regime are described, as are the pertinent observational signals. With these in hand a range of astrophysical tests of gravity are presented, and prospects for future measurements and theoretical developments are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Collaborative Research: Particle Simulation of Collisionless Magnetic Reconnection Under Finite Guide Field

Magnetic reconnection is believed to be one of the most fundamental physical processes in laboratory and space plasmas. Fundamental physics of collisionless reconnection under general conditions of finite guide magnetic field, realistic electron-to-ion mass ratio, and in 3D geometries is still not adequately understood. The major goal of this project is to study the physics of collisionless magnetic reconnection in a current sheet under a finite guide field and with a realistic ion-to-electron mass ratio mi/me. Under the support of this DoE grant, we have made significant progress in (1) the investigation of linear and nonlinear current sheet instabilities under a wide range of guide field using the 3-D gyro kinetic electron and fully kinetic ion (GeFi) particle simulation model and (2) the investigation of fully nonlinear reconnection process with a strong guide field with the GeFi simulations. Furthermore, (3) we have extended the study to the realistic geometry of the magnetopause using 3D global hybrid simulations. This study, from local to global scales and for 3-D physics at plasma boundaries, has provides advanced theoretical understanding of the fundamental physics of magnetic reconnection in laboratory and space plasma regimes on multi-scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Large-scale simulations of Floquet physics on near-term quantum computers

Abstract Periodically driven quantum systems exhibit a diverse set of phenomena but are more challenging to simulate than their equilibrium counterparts. Here, we introduce the Quantum High-Frequency Floquet Simulation (QHiFFS) algorithm as a method to simulate fast-driven quantum systems on quantum hardware. Central to QHiFFS is the concept of a kick operator which transforms the system into a basis where the dynamics is governed by a time-independent effective Hamiltonian. This allows prior methods for time-independent simulation to be lifted to simulate Floquet systems. We use the periodically driven biaxial next-nearest neighbor Ising (BNNNI) model, a natural test bed for quantum frustrated magnetism and criticality, as a case study to illustrate our algorithm. We implemented a 20-qubit simulation of the driven two-dimensional BNNNI model on Quantinuum’s trapped ion quantum computer. Our error analysis shows that QHiFFS exhibits not only a cubic advantage in driving frequency ω but also a linear advantage in simulation time t compared to Trotterization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Deformed Fredkin model for the ν = 5 / 2 Moore-Read state on thin cylinders

We propose a frustration-free model for the Moore-Read quantum Hall state on sufficiently thin cylinders with circumferences ≲ 7 magnetic lengths. While the Moore-Read Hamiltonian involves complicated long-range interactions between triplets of electrons in a Landau level, our effective model is a simpler one-dimensional chain of qubits with deformed Fredkin gates. We show that the ground state of the Fredkin model has high overlap with the Moore-Read wave function and accurately reproduces the latter's entanglement properties. Moreover, we demonstrate that the model captures the dynamical response of the Moore-Read state to a geometric quench, induced by suddenly changing the anisotropy of the system. We elucidate the underlying mechanism of the quench dynamics and show that it coincides with the linearized bimetric field theory. The minimal model introduced here can be directly implemented as a first step towards quantum simulation of the Moore-Read state, as we demonstrate by deriving an efficient circuit approximation to the ground state and implementing it on an IBM quantum processor. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Value Proposition of UV-Absorbers in PV Module Encapsulation

Various common crystalline silicon cell technologies were exposed to UVA radiation (1.24 Wm-2 nm-1 at 340 nm peak) on the front and back faces at 45 degrees Celsius for periods up to 3000 h, representing about 3 y of solar exposure in Phoenix, Arizona, USA. The resulting degradation of the open-circuit voltage and short-circuit current is presented. Of the various cell types examined, significant levels of degradation were seen in all cases. Less degradation was generally found after UV irradiation of cell fronts and older cell types, whereas a bifacial PERC type exposed on the rear showed about 25% degradation in short circuit current, attributable to lack of a diffused surface field. Selected cells were exposed to UV irradiation with the addition of long pass UV filters to replicate UV-absorbers in encapsulants. Modern cell designs are sensitive to UV-ID because of reduced or eliminated front and back surface field and increased dependence on high quality surface passivation. Single transformation of the independent variable (t, kW h/m^2) could be used to achieve a linear model of the data to extrapolate to 50 y. Solar Advisor Model (SAM) shows appropriate filtering of UV-irradiation can improve LCOE and net present value of plant. Some advanced cell types are seen to be UV-resistant (cell level solutions also exist). Solutions therefore exist on the cell, glass, and encapsulant level. Changes over time in each of these would also need to be considered (solarization, encapsulant browning...).

ENGINEERING,SOLAR ENERGY↗
DOE OSTI2021

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗
DOE OSTI2020

SympNets: Intrinsic structure-preserving symplectic networks for identifying Hamiltonian systems

Here, we propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantum dynamics of non-Hermitian many-body Landau-Zener systems

Here, we develop a framework to solve a large class of linearly driven non-Hermitian quantum systems. Such a class of models in the Hermitian scenario is commonly known as multistate Landau-Zener models. The non-Hermiticity is due to the anti-Hermitian couplings between the diabatic levels. We find that there exists a conservation law, unique to this class of models, that describes the simultaneous growth of the unnormalized wave functions. These models have practical applications in Bose-Einstein condensates, and they can describe the dynamics of multispecies bosonic systems. The conservation law relates to a pair-production mechanism that explains the dissociation of diatomic molecules into atoms. We provide a general framework for both solvable and semiclassically solvable non-Hermitian Landau-Zener models. Our findings will open alternative avenues for a number of diverse emergent phenomena in explicitly time-dependent non-Hermitian quantum systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Signatures of cysteine oxidation on muscle structural and contractile proteins are associated with physical performance and muscle function in older adults: Study of Muscle, Mobility and Aging (SOMMA)

Oxidative stress is considered a contributor to declining muscle function and mobility during aging; however, the underlying molecular mechanisms remain poorly described. We hypothesized that greater levels of cysteine (Cys) oxidation on muscle proteins are associated with decreased measures of mobility. Herein, we applied a novel redox proteomics approach to measure reversible protein Cys oxidation in vastus lateralis muscle biopsies collected from 56 subjects in the Study of Muscle, Mobility and Aging (SOMMA), a community-based cohort study of individuals aged 70 years and older. We tested whether levels of Cys oxidation on key muscle proteins involved in muscle structure and contraction were associated with muscle function (leg power and strength), walking speed, and fitness (VO 2 peak on cardiopulmonary exercise testing) using linear regression models adjusted for age, sex, and body weight. Higher oxidation levels of select nebulin Cys sites were associated with lower VO 2 peak, while greater oxidation of myomesin-1, myomesin-2, and nebulin Cys sites was associated with slower walking speed. Higher oxidation of Cys sites in key proteins such as myomesin-2, alpha-actinin-2, and skeletal muscle alpha-actin were associated with lower leg power and strength. We also observed an unexpected correlation (R = 0.48) between a higher oxidation level of eight Cys sites in alpha-actinin-3 and stronger leg power. Despite this observation, the results generally support the hypothesis that Cys oxidation of muscle proteins impairs muscle power and strength, walking speed, and cardiopulmonary fitness with aging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High- temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes ad- vantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Resonance three-wave interactions and strange attractor

We report it is shown that the incorporation of linear sink/source terms in the three-wave resonance interaction model results in the time dependence of the wave amplitudes, which could exhibit the properties of a strange attractor. This finding demonstrates that the transition to turbulent dynamics of the waves could be related not only to the coupling of wave triads but also to the establishing of the strange attractor-like dynamics within individual wave triads.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The Generalized Green’s function Cluster Expansion: A Python package for simulating polarons

We present an efficient implementation of the Generalized Green’s function Cluster Expansion (GGCE), which is a new method for computing the ground-state properties and dynamics of polarons (single electrons coupled to lattice vibrations) in model electron-phonon systems. The GGCE works at arbitrary temperature and is well suited for a variety of electron-phonon couplings, including, but not limited to, site and bond Holstein and Peierls (Su-SchriefferHeeger) couplings, and couplings to multiple phonon modes with different energy scales and coupling strengths. Quick calculations can be performed efficiently on a laptop using solvers from NumPy and SciPy, or in parallel at scale using the PETSc sparse linear solver engine.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Joint Topology Identification and State Estimation in Unobservable Distribution Grids

Many distribution system operations (e.g., state estimation, control, fault detection/localization) rely on the assumption that the underlying topology is accurately defined. In general, topology identification is a challenging problem in distribution systems as these systems are unobservable with a very limited number of available measurements. In this paper, we tackle this problem by designing a compressive sensing framework that jointly estimates the systems states and network topology via an integrated mixed integer nonlinear program (MINLP) formulation. Here, two reformulations of the original MINLP problems are investigated. Firstly, in order to remove the nonlinearity in the MINLP formulation, a mixed integer linear programming (MILP) problem that employs auxiliary variables is derived. Furthermore, to achieve a faster solution, convex relaxation of the original formulation is derived. Finally, using a Markovian model for topology changes, prior information about system topology is used to improve topology identification particularly when a limited amount of measurements is available. Simulation results on IEEE 37-bus test feeder and IEEE 123-bus test feeder illustrate the efficiency and scalability of the proposed approaches from both state estimation and topology identification point of view (even with 30% of available data).

42 ENGINEERING↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Implementation of Creep-Damage Model for Concrete Fracture in MOOSE

Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

3D-printed mirror-less helicity preserving metasurface “mirror” for THz applications

Stimulated by seminal works on generalized reflection and refraction laws, metasurfaces have evolved to a highly promising research direction, as they allow a multitude of different functionalities by optically thin wave-control elements/structures. Among them, structures functioning in the low THz regime are of extensive research interest, due to their high potential in communication and sensing applications. In this paper we propose a simple THz metasurface design that exhibits ideally perfect cross-polarized reflection for linear polarization and perfect helicity-preserving reflection with a geometric (Pancharatnam–Berry) phase for circular polarization, without the presence of any back-reflector. Numerical calculations demonstrating the structure response are justified by analytical models, which provide physical insights on this response. The designed metasurfaces are fabricated using the direct laser writing 3D-printing technology, metallized with electroless silver plating, and are characterized by THz time domain spectroscopy, with the experimental results validating the corresponding theoretical ones. Applications including beam steering and focusing, exploring the Pancharatnam–Berry phase, are also demonstrated numerically. Besides those applications, the helicity preserving mirror response of our metasurfaces can be valuable, among others, also in molecular chirality sensing applications, an issue that is also highlighted here.

3D-printed metasurfaces↗
DOE OSTI2024

NoRA: A Tensor Network Ansatz for Volume-Law Entangled Equilibrium States of Highly Connected Hamiltonians

Motivated by the ground state structure of quantum models with all-to-all interactions such as mean-field quantum spin glass models and the Sachdev-Ye-Kitaev (SYK) model, we propose a tensor network architecture which can accomodate volume law entanglement and a large ground state degeneracy. We call this architecture the non-local renormalization ansatz (NoRA) because it can be viewed as a generalization of MERA, DMERA, and branching MERA networks with the constraints of spatial locality removed. We argue that the architecture is potentially expressive enough to capture the entanglement and complexity of the ground space of the SYK model, thus making it a suitable variational ansatz, but we leave a detailed study of SYK to future work. We further explore the architecture in the special case in which the tensors are random Clifford gates. Here the architecture can be viewed as the encoding map of a random stabilizer code. We introduce a family of codes inspired by the SYK model which can be chosen to have constant rate and linear distance at the cost of some high weight stabilizers. We also comment on potential similarities between this code family and the approximate code formed from the SYK ground space.

Physics↗