Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Topology of the $O(3)$ non-linear sigma model under the gradient flow

The $O(3)$ non-linear sigma model (NLSM) is a prototypical field theory for QCD and ferromagnetism, featuring topological qualities. Though the topological susceptibility should vanish in physical theories, lattice simulations of the NLSM find that the topological susceptibility diverges in the continuum limit. We introduce a $\theta$-term and analyze the topological charge as a function of $\theta$ under the gradient flow. Our results show that divergence persists in the presence of the flow, even at non-zero $\theta$.

Thomas, Stuart↗
DOE OSTI2020

Partitioned exponential methods for coupled multiphysics systems

Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dynamically reconfigurable topological routing in nonlinear photonic systems

The propagation path of topologically protected states is bound to the interface between regions with different topology, and as such, the functionality of linear photonic devices leveraging these states is fixed during fabrication. Here, we propose a mechanism for dynamic control over a driven dissipative system’s local topology, yielding reconfigurable topological interfaces and thus tunable paths for protected routing. We illustrate our approach in non-resonantly pumped polariton lattices, where the nonlinear interaction between the polaritons and the exciton reservoir due to non-resonant pumping can yield picosecond-scale changes in the propagation paths of the chiral edge states. To analytically confirm the numerically observed topological dynamics, we generalize the spectral localizer framework to non-linear non-Hermitian Chern materials and apply this framework to a continuous model of the polariton system based on a driven-dissipative Gross-Pitaevskii equation. In doing so, we show that the local changes in the polariton lattice’s topology are captured by a local Chern marker. Looking forward, we anticipate such reconfigurable topological routing will enable the realization of novel classes of topological photonic devices.

Wong, Stephan [Sandia National Laboratories (SNL-N↗
DOE OSTI2021

Jet wake from linearized hydrodynamics

We explore how to improve the hybrid model description of the particles originating from the wake that a jet produced in a heavy ion collision leaves in the droplet of quark-gluon plasma (QGP) through which it propagates, using linearized hydrodynamics on a background Bjorken flow. Jet energy and momentum loss described by the hybrid model become currents sourcing linearized hydrodynamics. By solving the linearized hydrodynamic equations numerically, we investigate the development of the wake in the dynamically evolving droplet of QGP, study the effect of viscosity, scrutinize energy-momentum conservation, and check the validity of the linear approximation. We find that linearized hydrodynamics works better in the viscous case because diffusive modes damp the energy-momentum perturbation produced by the jet. We calculate the distribution of particles produced from the jet wake by using the Cooper-Frye prescription and find that both the transverse momentum spectrum and the distribution of particles in azimuthal angle are similar in shape in linearized hydrodynamics and in the hybrid model. Their normalizations are different because the momentum-rapidity distribution in the linearized hydrodynamics analysis is more spread out, due to sound modes. Since the Bjorken flow has no transverse expansion, we explore the effect of transverse flow by using local boosts to add it into the Cooper-Frye formula. After including the effects of transverse flow in this way, the transverse momentum spectrum becomes harder: more particles with transverse momenta bigger than 2 GeV are produced than in the hybrid model. Although we defer implementing this analysis in a jet Monte Carlo, as would be needed to make quantitative comparisons to data, we gain a qualitative sense of how the jet wake may modify jet observables by computing proxies for two example observables: the lost energy recovered in a cone of varying open angle, and the fragmentation function. We find that linearized hydrodynamics with transverse flow effects added improves the description of the jet wake in the hybrid model in just the way that comparison to data indicates is needed. Our study illuminates a path to improving the description of the wake in the hybrid model, highlighting the need to take into account the effects of both transverse flow and the broadening of the energy-momentum perturbation in spacetime rapidity on particle production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Deterministic and Monte Carlo Nuclear Data Adjustment Methods [Slides]

For the Bayesian Monte Carlo methodology, a need to understand convergence of the posterior moments as a function of the number of parameter realizations is required. In high-dimensional systems, it can be very costly to sample entire parameter space and perform functional evaluation for every realization. Bayesian Monte Carlo allows one to relax the GLLS approximations of model linearity and prior/posterior PDF shape. The Bayesian Stochastic Collocation Method is a deterministic approach to “sample” the parameter space. It allows one to relax the GLLS approximations of model linearity and posterior PDF shape. Higher-order posterior moments (i.e., skewness, kurtosis, etc.) can be studied through polynomial expansion. Tensor product quadrature scales poorly and can use sparse grid quadrature methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hidden zeros for particle/string amplitudes and the unity of colored scalars, pions and gluons

Recent years have seen the emergence of a new understanding of scattering amplitudes in the simplest theory of colored scalar particles — the Tr(Φ 3 ) theory — based on combinatorial and geometric ideas in the kinematic space of scattering data. In this paper we report a surprise: far from the toy model it appears to be, the “stringy” Tr(Φ 3 ) amplitudes secretly contains the scattering amplitudes for pions, as well as non-supersymmetric gluons, in any number of dimensions. The amplitudes for the different theories are given by one and the same function, related by a simple shift of the kinematics. This discovery was spurred by another fundamental observation: the tree-level Tr(Φ 3 ) field theory amplitudes have a hidden pattern of zeros when a special set of non-planar Mandelstam invariants is set to zero. These zeros are not manifest in Feynman diagrams but are made obvious by the connection of these amplitudes to the new understanding of associahedra arising from “causal diamonds” in kinematic space. Furthermore, near these zeros, the amplitudes simplify, by factoring into a non-trivial product of smaller amplitudes. Remarkably the amplitudes for pions and gluons are observed to also vanish in the same kinematical locus. These properties for Tr(Φ 3 ) amplitudes hold and further generalize to the “stringy” Tr(Φ 3 ) amplitudes. The “kinematic causal diamond” picture suggests a unique shift of the kinematic data that preserves the zeros, and this shift is precisely the one that unifies colored scalars, pions, and gluons into a single object. We will focus in this paper on explaining the hidden zeros and factorization properties and the connection between all the colored theories, working for simplicity at tree level. Subsequent works will describe this new formulation for the Non-linear Sigma Model and non-supersymmetric Yang-Mills theory, at all loop orders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Comfort Model Simplification for Tight Integration with Grid Service Optimizations

Localized generation, storage, and intelligent power electronics are increasingly being integrated into building. The excess storage, thermal, and generation capacity can be used to develop new markets for ancillary services such as peak reduction, voltage and frequency support, or load flattening. Widespread participation in these markets will make the power grid more resilient. However, widespread participation in these markets will require a simple and intuitive control over the excess capacity by the building occupants. Most importantly for adopting these technologies, participation in these markets should not require sacrificing occupant comfort. Many of the proposed control and optimization technologies for supplying the ancillary services do not explicitly incorporate occupant comfort into their models and when they do, generally utilize a rigid constraint on the indoor air temperature set a priori. This reduces flexibility in using excess thermal capacity. It also cannot take into account uncertainties in local factors that affect thermal comfort such as clothing insulation or relative humidity. In this paper we present and justify some assumptions that simplify the standard predictors of thermal comfort. We also develop a regression model of the thermal comfort that is linear with respect to the indoor air temperature. Next we develop the concept of thermal comfort variation that eliminates the dependence of the thermal comfort on uncertain thermal comfort factors. Then, we present a quadratic program utilizing the thermal comfort variation in both the objective function and as a constraint. Finally, we show the results of a time-of-use cost optimization that utilizes the simplified thermal comfort model.

Melin, Alexander↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2022

3D spherical functional expansion tallies in Serpent 2 Monte Carlo code

This work extends the application of functional expansion tallies to 3D spherical geometries. The 3D Zernike polynomials are set as an orthonormal polynomials basis for the functional reconstruction. The study describes the construction of the complete set of polynomials, a natural expansion of the spherical harmonics polynomials where 3D Zernike moments can be evaluated as a linear combination of the geometrical moments. The 3D Zernike polynomials formulation and the computational approach implemented in Serpent 2 are presented and tested through the Godiva model from the ICSBEP criticality benchmark test cases. The implementation results are in agreement with a reference solution described in a fine-resolution mesh, enhancing also the performance and memory demand. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Quantum time dynamics employing the Yang-Baxter equation for circuit compression

Quantum time dynamics (QTD) is considered a promising problem for quantum supremacy on near-term quantum computers. However, QTD quantum circuits grow with increasing time simulations. This study focuses on simulating the time dynamics of one-dimensional (1D) integrable spin chains with nearest-neighbor interactions. We have proved the existence of a reflection symmetry in the quantum circuit employed for simulating the time evolution of certain classes of 1D Heisenberg model Hamiltonians by virtue of the quantum Yang-Baxter equation, and how this symmetry can be exploited to compress and produce a shallow quantum circuit. With this compression scheme, the depth of the quantum circuit becomes independent of step size and only depends on the number of spins. We show that the depth of the compressed circuit is rigorously a linear function of the system size for the studied Heisenberg model Hamiltonians in the present work. As a consequence, the number of CNOT gates in the compressed circuit only scales quadratically with the system size, which allows for the simulations of time dynamics of very large 1D spin chains. We derive the compressed circuit representations for different special cases of the Heisenberg Hamiltonian. We compare and demonstrate the effectiveness of this approach by performing simulations on quantum computers.

1-dimensional spin chains↗
DOE OSTI2021

A new and alternative look at nonlinear Alfvénic states

The formulation for studying nonlinear Alfvenic states, sustainable in Hall Magnetohydrodynamics (HMHD), becomes considerably simpler and more tractable when circularly polarized Beltrami vectors (the eigenstates of linear HMHD) are used as the basis functions. Nonlinear HMHD is, then, reduced to a rather simple looking set of scalar equations from which a model problem of three interacting Beltrami modes is formulated and analytically solved. The triplet interactions span a variety of familiar nonlinear processes leading to a redistribution as well as periodic exchange of energy. The energy exchange processes (whose strength is measured by an energy exchange/depletion time) will, perhaps, play a dominant role in determining the spectral content of an eventual Alfvenic state. All nonlinearities (sensitive functions of the interacting wave vectors) operate at par, and none is dominant over any substantial region of k-space; their intricate interplay prevents a “universal” picture from emerging; few generalizations on the processes that may, for instance, lead to a turbulent state, are possible. However, the theory can definitely claim: (1) the energy tends to flow from lower to higher k and (2) the higher k z (in the direction of the ambient magnetic field) components of a mode with a given k are depleted/oscillate faster—in some cases much faster. It is noteworthy that the mode coupling is the strongest (with the shortest depletion time) when the participating wave vectors are nearly perpendicular; perhaps, an expected consequence of the curl (cross product) nonlinearities. Numerical simulations will be necessary to help create a fully reliable picture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Bouger's Law Shell Ionospheric Transfer Function

Ionospheric transfer function (ITF) algorithms determine the effects of the ionosphere on an electromagnetic (EM) radio-frequency (RF) signal as it propagates through. In this report, the Bouger’s law shell model is outlined. This ITF is very similar to the Snell’s Law Shell ITF; however, in this formulation, plasma parameters in the ionosphere are allowed to change radially. This algorithm is expressed in the frequency domain. In this way, it is applied as linear time invariant (LTI) filter function. Signals in this report are assumed to have only a single component (i.e. x, y or z in a rectangular coordinate system). Multi-component signals can be treated simply by applying the specific ITF to each component separately.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.

Angular momentum↗
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

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗