Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physical systems”

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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Characterizing Structure-Dependent TiS 2 /Water Interfaces Using Deep-Neural-Network-Assisted Molecular Dynamics

As a promising layered electrode material, TiS 2 -based capacitive deionization (CDI) devices for water desalination have attracted significant attention. However, TiS 2 /H 2 O interfacial features, potentially important for device optimization, remain unidentified. Using Deep Potential Molecular Dynamics (DPMD), we characterized distinct aqueous interfaces introduced by four TiS 2 terminations expected to be present as water intercalates into TiS 2 , namely, Armchair, Zigzag, Zigzag-L, and Zigzag-R. First, we assessed important representative physical properties of the system to validate the deep potentials (DPs). DPMD simulations agree well with experiments and first-principles simulations, suggesting the DPs are accurate and reliable. Subsequent simulations of these TiS 2 /water interfaces revealed how TiS 2 surface termination influences the structure of interfacial water. This effect is most evident in the first and second water layers close to the TiS 2 surface, and more pronounced when spontaneous dissociative adsorption of water occurs. The extent of water dissociation on each surface was evaluated using enhanced sampling. Zigzag-L is the only interface where proton transfer from adsorbed water to TiS 2 surface S atoms is thermodynamically and kinetically favored. The coexistence of surface four-fold-coordinated Ti (Ti 4c ) and one-fold-coordinated S (S 1c ) is found to be essential to making proton transfer feasible on the Zigzag-L surface. Furthermore, remaining unprotonated S 1c atoms can act as good proton acceptors after water dissociation. Thus, TiS 2 with Zigzag-L termination may be a surface to avoid in CDI device construction, given that pH fluctuations adversely affect performance. Furthermore, this work provides new understanding of TiS 2 /H 2 O interfacial features that could aid future design and optimization of TiS 2 -based CDI devices for water desalination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Machine-Learning-Assisted Hybrid Earth System Modelling

We propose a novel hybrid modeling approach that combines a physics-based, numerical model of the Earth system with a machine learning model. We also envision to integrate the hybrid model training with a data assimilation cycle, so that new observations update both the state estimate used as initial conditions and the machine learning model parameters.

58 GEOSCIENCES↗
DOE OSTI2024

PRISTINE: An Emulation Platform for PCB-Level Hardware Trojans

Printed circuit Boards (PCBs) are becoming increasingly vulnerable to malicious design alteration, also known as Trojan attacks, due to a distributed business model that often involves various untrusted parties. Such attacks can be mounted at various stages in the PCB life cycle. The relative ease of alteration of PCB hardware even after fabrication (due to physical access to surface-mounted critical components and traces) makes them attractive for an adversary to manipulate their functional/physical behavior for malicious intent. There is a growing need to explore viable Trojan attacks in a PCB, analyze their functional and physical characteristics (e.g., impact on power or delay), and study the effectiveness of countermeasures against these attacks. While simulation-based approaches for PCB Trojan insertion are effective at creating a large population of possible Trojans, they fail to provide functional feasibility analysis with a realistic workload for a trigger circuit. Also, they cannot estimate a Trojan’s side-channel footprint due to the unavailability of physical models of diverse PCB components. To address these deficiencies, in this paper, we present PRISTINE, a PCB-level emulation system for any integrity or physical tampering issues, specifically, hardware Trojan insertion. The need for building such an emulation platform to resolve PCB trust issues in the supply chain is also surveyed and discussed. Both custom Hardware Hacking (HaHa) boards and multiple commercial PCBs are then used to test the ability of the proposed system to emulate various hardware Trojans specially designed to exploit board-specific hardware characteristics. Experimental results on emulated board-level Trojans show that a wide range of Trojans can be successfully activated, thus enabling the expected payload effects on both types of boards to be studied and quantified. The resulting data are further analyzed to create PCB-level Trojan benchmarks. In particular, a comparative evaluation of the experimental results is used to propose a risk level metric that quantifies the probability of detection and degree of payload impact of each Trojan on a given commercial PCB.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Local SU(2) × U(1) quark flavor symmetry in the RS bulk

We propose a model of quark flavor based on an additional SU(2) × U(1) local symmetry in a warped extra dimensional bulk. In contrast to other works, we break the additional gauge symmetry in the bulk via two complex scalars which acquire bulk vevs, rather than relying on brane-localized symmetry breaking. A gauge-covariant Kaluza-Klein decomposition of a theory with a bulk spontaneously broken gauge symmetry is performed, and exact expressions for the bulk profiles of all physical particles in such systems are given. The SM quark masses and mixings are then recreated using gauge-covariant bulk quark mass terms and Yukawa-like couplings to the new bulk scalars. A numerical sampling of points in the model parameter space that recreate the quark masses and mixings is performed at a KK scale of M KK = 5 TeV. We then compute the ΔF = 2 4-quark operators arising from our new flavor gauge bosons and scalars, and those arising from Kaluza-Klein modes of SM gauge bosons. By decoupling one of our bulk scalar fields to all quark fields except the right-handed up-like sector, we find that it is possible to greatly suppress tree-level contributions to the highly constrained Kaon mixing parameters. Instead, the dominant constraints on the model emerge from neutral Bd and D meson mixing. These constraints are explored with our numerical sampling of the model parameter space, and the specific contribution of the new flavor gauge bosons and scalars is discussed. We find that for a significant range of realistic flavor gauge couplings, the new gauge bosons compete with the normally dominant gluon flavor-changing currents, but flavor-changing operators emerging from the bulk scalar fields are highly suppressed. Finally, we briefly comment on flavor constraints that are independent of the flavor gauge sector arising from the Z$\overline{b}$ L b L coupling and rare top decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Structural and magnetic properties of Pt/Co/Mn-based multilayers

Magnetic multilayers are a rich class of material systems with numerous highly tunable physical parameters that determine both their magnetic and electronic properties. Here, we present a comprehensive experimental study of the system Pt/Co/Mn, which extends the group of Pt / Co / X ( X = metal) multilayers that have been investigated thus far. We demonstrate that an increasing Co layer thickness changes the magnetic anisotropy from out of plane to in plane, whereas the deposition of thicker Mn layers leads to a decrease in the saturation magnetization. Temperature-dependent magnetometry measurements reinforce the hypothesis of antiferromagnetic coupling at the Co/Mn interfaces being responsible for the observed Mn thickness dependence of the magnetization reversal. Moreover, magneto-optical imaging experiments indicate systematic changes in magnetic domain patterns as a function of the Co and Mn layer thickness, suggesting the existence of bubblelike domains—potentially even magnetic skyrmions—in the case of sufficiently thick Mn layers, which are expected to contribute to a sizable Dzyaloshinskii-Moriya interaction in the multilayer stacks. As a result, we identify Pt/Co/Mn as a highly complex multilayer system with strong potential for further fundamental studies and possible applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Polar System Analysis Package (Ver. 1 Specifications)

The Polar System Analysis Package (POLeSTAr) is a modern data fusion tool being designed to benchmark,evaluate and compare the polar physics and biogeochemistry of Earth System Models (ESMs). It is being written in Python to make use of libraries including Xarray, Matplotlib, SciPy, NumPy and PyTorch within algorithms developed specifically for polar system modeling. Initial POLeSTAr algorithms will be adapted from the MATLAB Ridgepack package that was designed to analyze the Model for Prediction Across Scales (MPAS) sea ice component of the Energy Exascale Earth System Model (E3SM) and sea ice model CICE within the Community Earth System Model (CESM).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The X-Ray Emission Effectiveness of Plasma Mirrors: Reexamining Power-Law Scaling for Relativistic High-Order Harmonic Generation

Ultrashort pulsed lasers provide uniquely detailed access to the ultrafast dynamics of physical, chemical, and biological systems, but only a handful of wavelengths are directly produced by solid-state lasers, necessitating efficient high-power frequency conversion. Relativistic plasma mirrors generate broadband power-law spectra, that may span the gap between petawatt-class infrared laser facilities and x-ray free-electron lasers; despite substantial theoretical work the ultimate efficiency of this relativistic high-order-harmonic generation remains unclear. We show that the coherent radiation emitted by plasma mirrors follows a power-law distribution of energy over frequency with an exponent that, even in the ultrarelativistic limit, strongly depends on the ratio of laser intensity to plasma density and exceeds the frequently quoted value of -8/3 over a wide range of parameters. The coherent synchrotron emission model, when adequately corrected for the finite width of emitting electron bunches, is not just valid for p-polarized light and thin foil targets, but generally describes relativistic harmonic generation, including at normal incidence and with finite-gradient plasmas. Our numerical results support the ω-4/3 scaling of the synchrotron emission model as a limiting efficiency of the process under most conditions. The highest frequencies that can be generated with this scaling are usually restricted by the width of the emitting electron bunch rather than the Lorentz factor of the fastest electrons. The theoretical scaling relations developed here suggest, for example, that with a 20-PW 800-nm driving laser, 1 TW/harmonic can be produced for 1-keV photons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coulomb Interaction-Driven Entanglement of Electrons on Helium

The generation and evolution of entanglement in many-body systems is an active area of research that spans multiple fields, from quantum information science to the simulation of quantum many-body systems encountered in condensed matter, subatomic physics, and quantum chemistry. Motivated by recent experiments exploring quantum information processing systems with electrons trapped above the surface of cryogenic noble gas substrates, we theoretically investigate the generation of entanglement between two electrons via their unscreened Coulomb interaction. The model system consists of two electrons confined in separate electrostatic traps that establish microwave-frequency quantized states of their motion. We compute the motional energy spectra of the electrons, as well as their entanglement, by diagonalizing the model Hamiltonian with respect to a single-particle Hartree product basis. We also compare our results with the predictions of an effective Hamiltonian. The computational procedure outlined here can be employed for device design and guidance of experimental implementations. In particular, the theoretical tools developed here can be used for fine-tuning and optimization of control parameters in future experiments with electrons trapped above the surface of superfluid helium or solid neon. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

Direct Bayesian inference for fault severity assessment in Digital-Twin-Based fault diagnosis

For applications in condition-based maintenance of nuclear systems, the assessment of fault severity is crucial. In this work, we developed a framework that allows for direct inference of the probability distributions of possible faults in a system. Here, we employed a model-based approach with model residuals generated from analytical redundancy relations provided by physics-based models of the system components. From real-time sensor readings, the values of the model residuals can be calculated, and the posterior probability distributions of the faults can be computed directly using the methods of Bayesian networks. From the posterior distribution of each fault, one can estimate the fault probability based on a chosen threshold and assess the severity of the fault. By eliminating the discretization and simplifications in middle steps, this approach allows us to leverage the available computational resources to provide more accurate fault probability estimates and severity assessments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Dynamic Parameter Estimation with Physics-based Neural Ordinary Differential Equations

Accurate estimation of dynamic parameters of gen-erators is crucial to building a reliable model for dynamical studies and reliable operation of the power system. This paper develops a physics-based neural ordinary differential equations (ODE) approach to learn the parameters of generator dynamic model using phasor measurement units (PMU) data. We design a physics-based neural network to represent the swing equations of the power system dynamics. A loss function is defined as the difference between dynamic simulation results from the physics-based neural networks and pseudo PMU measurements. The parameters of generator dynamic model are iteratively updated using the neural ODEs and the adjoint method. By exploiting the mini-batch scheme in neural ODE training, the parameter estimation performance is significantly improved. Numerical study results on a 3-machine 9-bus system show that the proposed algorithm outperforms state-of-the-art baseline method in both computation time and dynamic parameter estimation accuracy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Development of data‐driven exponential integrators with application to modeling of delay photocurrents

Abstract Radiation‐induced photocurrent effects represent a threat to microelectronic components operating in space, manmade and terrestrial radiation environments. Analysis of these threats by circuit simulations requires accurate and computationally efficient compact models. Most existing compact models are based on closed form analytic solutions of the governing equations and require empirical assumptions and idealizations that can limit their validity. In this paper we formulate an alternative numerical, data‐driven approach that learns a compact model from data representative of the type of measurements one can obtain in an experimental facility. To develop the model we start from a generic discrete‐time dynamical system and then use physics knowledge to refine its structure. Numerical studies demonstrate the potential of the model and establish some empirical guidelines for its training.

Bochev, Pavel↗
DOE OSTI2021

Reflectometry with Polarized Neutrons on In Situ Grown Thin Films

Originating from the demand for obtaining depth-resolved magnetization profiles from thin films and heterostructures, polarized neutron reflectometry (PNR) has developed into a unique research tool, which also finds application in the analysis of superconducting or soft matter thin films. While certain in situ sample environments such as gas-loading or humidity cells were quickly realized after PNR first emerged, preparing and growing thin magnetic films directly in the neutron beam could only be realized in recent years. Herein, a dedicated insight is given on the history and development of in situ thin film growth capabilities for PNR, from early pioneering experiments to the present day. The scientific and technological challenges as well as the advances of neutron sources, neutronics, and data treatment that have led to its realization are highlighted together with the unique research opportunities that it provides and recently obtained experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Generalized symmetries of the graviton

We find the set of generalized symmetries associated with the free graviton theory in four dimensions. These are generated by gauge invariant topological operators that violate Haag duality in ring-like regions. As expected from general QFT grounds, we find a set of “electric” and a dual set of “magnetic” topological operators and compute their algebra. To do so, we describe the theory using phase space gauge-invariant electric and magnetic dual variables constructed out of the curvature tensor. Electric and magnetic fields satisfy a set of constraints equivalent to the ones of a stress tensor of a 3d CFT. The constraints give place to a group $\mathbb{R}^{20}$ of topological operators that are charged under space-time symmetries. Finally, we discuss similarities and differences between linearized gravity and tensor gauge theories that have been introduced recently in the context of fractonic systems in condensed matter physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Spectral fluctuations in the Sachdev-Ye-Kitaev model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N ) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q= 4 or higher even q-body interactions, but also for the q= 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Building bulk geometry from the tensor Radon transform

Using the tensor Radon transform and related numerical methods, we study how bulk geometries can be explicitly reconstructed from boundary entanglement entropies in the specific case of AdS 3 /CFT 2 . We find that, given the boundary entanglement entropies of a 2d CFT, this framework provides a quantitative measure that detects whether the bulk dual is geometric in the perturbative (near AdS) limit. In the case where a well-defined bulk geometry exists, we explicitly reconstruct the unique bulk metric tensor once a gauge choice is made. We then examine the emergent bulk geometries for static and dynamical scenarios in holography and in many-body systems. Apart from the physics results, our work demonstrates that numerical methods are feasible and effective in the study of bulk reconstruction in AdS/CFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dynamic modelling of slip in a wind turbine spherical roller main bearing

This paper considers the problem of the dynamic modelling of macro slip in spherical roller bearings. By revisiting the fundamental physics which drive these systems, potential issues in existing models have been identified. Furthermore, in pure rolling conditions it was found that governing differential equations become "stiff", requiring the use of implicit methods of time integration. The problem of individual roller macro slip in a wind turbine main bearing is then investigated using a simplified representation of system dynamics. Model results indicate clear links between slip/friction and the operational strategy of the wind turbine, as well as significantly higher frictional effects in the downwind main bearing row. Due to modelling simplifications, these results should not yet be considered conclusive, with further work required.

17 WIND ENERGY↗