SEARCH · Engineering Papers

Results for “linear embedding”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Local-Field Effects in Linear Response Properties within a Polarizable Frozen Density Embedding Method

In this work, we present a polarizable frozen density embedding (FDE) method for calculating polarizabilities of coupled subsystems. The method (FDE-pol) combines a FDE method with an explicit polarization model such that the expensive freeze/thaw cycles can be bypassed, and approximate nonadditive kinetic potentials are avoided by enforcing external orthogonality between the subsystems. To describe the polarization of the frozen environment, we introduce a Hirshfeld partition-based density-dependent method for calculating the atomic polarizabilities of atoms in molecules, which alleviates the need to fit the atomic parameters to a specific system of interest or to a larger general set of molecules. Further, we show that the Hirshfeld partition-based method predicts molecular polarizabilities close to the basis set limit, and thus, a single basis set-dependent scaling parameter can be introduced to improve the agreement against the reference polarizability data. To test the model, we characterized the uncoupled and coupled response of small interacting molecular complexes. Here, the coupled response properties include the perturbation of the frozen system due to the external perturbation which is ignored in the uncoupled response. We show that FDE-pol can accurately reproduce both the exact uncoupled polarizability and the coupled polarizabilities of the supermolecular systems. Using damped response theory, we also demonstrate that the coupled frequency-dependent polarizability can be described by including local field effects. The results emphasize the necessity of including local-field effects for describing the response properties of coupled subsystems, as well as the importance of accurate atomic polarizability models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Quantum computing for fusion energy science applications

This is a review of recent research exploring and extending present-day quantum computing capabilities for fusion energy science applications. We begin with a brief tutorial on both ideal and open quantum dynamics, universal quantum computation, and quantum algorithms. Then, we explore the topic of using quantum computers to simulate both linear and nonlinear dynamics in greater detail. Because quantum computers can only efficiently perform linear operations on the quantum state, it is challenging to perform nonlinear operations that are generically required to describe the nonlinear differential equations of interest. In this work, we extend previous results on embedding nonlinear systems within linear systems by explicitly deriving the connection between the Koopman evolution operator, the Perron–Frobenius evolution operator, and the Koopman–von Neumann evolution (KvN) operator. We also explicitly derive the connection between the Koopman and Carleman approaches to embedding. Extension of the KvN framework to the complex-analytic setting relevant to Carleman embedding, and the proof that different choices of complex analytic reproducing kernel Hilbert spaces depend on the choice of Hilbert space metric are covered in the appendixes. Finally, we conclude with a review of recent quantum hardware implementations of algorithms on present-day quantum hardware platforms that may one day be accelerated through Hamiltonian simulation. We discuss the simulation of toy models of wave–particle interactions through the simulation of quantum maps and of wave–wave interactions important in nonlinear plasma dynamics.

Joseph, I. (ORCID:0000000255400840)↗

DOE OSTI2022

Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures

Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2024

Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent 3-dimensional simulation of plasma-cascade instability using code SPACE.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Topological network analysis of patient similarity for precision management of acute blood pressure in spinal cord injury

Background: Predicting neurological recovery after spinal cord injury (SCI) is challenging. Using topological data analysis, we have previously shown that mean arterial pressure (MAP) during SCI surgery predicts long-term functional recovery in rodent models, motivating the present multicenter study in patients. Methods: Intra-operative monitoring records and neurological outcome data were extracted (n = 118 patients). We built a similarity network of patients from a low-dimensional space embedded using a non-linear algorithm, Isomap, and ensured topological extraction using persistent homology metrics. Confirmatory analysis was conducted through regression methods. Results: Network analysis suggested that time outside of an optimum MAP range (hypotension or hypertension) during surgery was associated with lower likelihood of neurological recovery at hospital discharge. Logistic and LASSO (least absolute shrinkage and selection operator) regression confirmed these findings, revealing an optimal MAP range of 76–[104-117] mmHg associated with neurological recovery. Conclusions: We show that deviation from this optimal MAP range during SCI surgery predicts lower probability of neurological recovery and suggest new targets for therapeutic intervention. Funding: NIH/NINDS: R01NS088475 (ARF); R01NS122888 (ARF); UH3NS106899 (ARF); Department of Veterans Affairs: 1I01RX002245 (ARF), I01RX002787 (ARF); Wings for Life Foundation (ATE, ARF); Craig H. Neilsen Foundation (ARF); and DOD: SC150198 (MSB); SC190233 (MSB); DOE: DE-AC02-05CH11231 (DM).

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

AMReX v2024

The software framework, AMReX, supports the development of block-structured adaptive mesh refinement (AMR) algorithms for solving systems of partial differential equations. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving the essential local descriptions of different physical processes in complex multiphysics algorithms. AMR uses a hierarchical representation of the solution at multiple levels of resolution where the solution on each level is defined on the union of data containers at that resolution. These data containers, which represent the solution over a logically rectangular subregion of the domain, can contain field data defined on a mesh, Lagrangian particles or combinations of both. In addition to these basic data types, AMReX supports a multilevel embedded boundary representation of complex geometry; linear solvers for cell-centered and nodal data; asynchronous I/O in a native format readable by ParaView, VisIt and yt; and interfaces to hypre and PETSc solvers. AMReX enables applications to run on distributed memory architectures with multicore CPUs and with GPU accelerators. AMReX uses a lightweight abstraction layer that effectively hides the details of the architecture from the application. The framework currently supports CUDA, HIP and SYCL for GPU acceleration and OpenMP for multi-core CPU architectures.

Almgren, Ann↗

DOE OSTI2021

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗

DOE OSTI2024

Novel plasma actuator for mitigation of dynamic stall

A novel plasma actuator, the Linear Counter-flow using a Point Embedded Electrode (LCPEE), is developed for the prevention of dynamic stall for a sinusoidal pitching movement between α = 4° and α = 18°. The LCPEE is implemented on a NACA0012 airfoil and tested at a Reynolds Number of Re c = 2 × 105 and reduced frequency of k = π/16. Prior investigations using a standard linear actuator showed that the exposed electrode introduced perturbations passively which delayed dynamic stall when the actuator was off. For the LCPEE actuator, when turned off, there is least passive delay. When the LCPEE is turned on at St f = 50, the dynamic stall is prevented for the sinusoidal pitching motion of the airfoil. The LCPEE actuator is also tested for the same sinusoidal motion between α = 6° and α = 20°. Four cases are considered for the higher α range of motion: actuator off, actuator on at St f = 50 with a sinusoidal input waveform, actuator on at St f = 50 with a triangular input waveform, and actuator on at Stf = 100 with a sinusoidal input waveform. Specifically for the last case, experiment shows no flow reversal demonstrating the efficacy of LCPEE in controlling the dynamic stall. The effects of LCPEE on the flow energy distribution have also been studied by using proper orthogonal decomposition (POD) method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

Here, we investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet–triplet excitation energy is 2.31 eV and that for the optically inactive singlet–singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV – defect using NEVPT2 for periodic supercell models.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A Mixed integer linear programming‐based distributed energy management for networked microgrids considering network operational objectives and constraints

Abstract Mixed integer linear programming (MILP)–based distributed energy management for networked microgrids embedded modern distribution systems is proposed. Considering the diverse ownership of microgrids, distributed energy resources (DERs) that interface directly with utilities and responsive loads, an alternating direction method of multipliers–based distributed framework was formulated for the scheduling of networked microgrids embedded modern distribution systems by adjusting nodal price signals iteratively. In addition, to make the formulated optimization problems resolvable through more accessible and popular MILP solvers, different linearisation techniques were employed to transform the nonlinear terms into linear or mixed integer linear formats. The proposed MILP‐based distributed method preserves all participants' autonomy (e.g., microgrids, DERs that interface directly with utilities and responsive loads), while incentivising them to actively participate in the distribution system operation with price signals. The proposed method is validated with results of numerical simulation using a modern distribution system consisting of multiple networked microgrids, DERs that interface directly with utilities, as well as responsive loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Why it is Unfortunate that Linear Machine Learning “Works” so well in Electromechanical Switching of Ferroelectric Thin Films

Machine learning (ML) is relied on for materials spectroscopy. It is challenging to make ML models fail because statistical correlations can mimic the physics without causality. Here, using a benchmark band-excitation piezoresponse force microscopy polarization spectroscopy (BEPS) dataset the pitfalls of the so-called “better”, “faster”, and “less-biased” ML of electromechanical switching are demonstrated and overcome. Using a toy and real experimental dataset, it is demonstrated how linear nontemporal ML methods result in physically reasonable embedding (eigenvalues) while producing nonsensical eigenvectors and generated spectra, promoting misleading interpretations. A new method of unsupervised multimodal hyperspectral analysis of BEPS is demonstrated using long-short-term memory (LSTM) β-variational autoencoders (β-VAEs) . By including LSTM neurons, the ordinal nature of ferroelectric switching is considered. Further, to improve the interpretability of the latent space, a variational Kullback–Leibler-divergency regularization is imposed . Finally, regularization scheduling of β as a disentanglement metric is leveraged to reduce user bias. Combining these experiment-inspired modifications enables the automated detection of ferroelectric switching mechanisms, including a complex two-step, three-state one. Ultimately, this work provides a robust ML method for the rapid discovery of electromechanical switching mechanisms in ferroelectrics and is applicable to other multimodal hyperspectral materials spectroscopies.

36 MATERIALS SCIENCE↗

DOE OSTI2023

SpeckleNN: a unified embedding for real-time speckle pattern classification in X-ray single-particle imaging with limited labeled examples

With X-ray free-electron lasers (XFELs), it is possible to determine the three-dimensional structure of noncrystalline nanoscale particles using X-ray single-particle imaging (SPI) techniques at room temperature. Classifying SPI scattering patterns, or `speckles', to extract single-hits that are needed for real-time vetoing and three-dimensional reconstruction poses a challenge for high-data-rate facilities like the European XFEL and LCLS-II-HE. Here, we introduce SpeckleNN, a unified embedding model for real-time speckle pattern classification with limited labeled examples that can scale linearly with dataset size. Trained with twin neural networks, SpeckleNN maps speckle patterns to a unified embedding vector space, where similarity is measured by Euclidean distance. We highlight its few-shot classification capability on new never-seen samples and its robust performance despite having only tens of labels per classification category even in the presence of substantial missing detector areas. Without the need for excessive manual labeling or even a full detector image, our classification method offers a great solution for real-time high-throughput SPI experiments.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Contrastive Machine Learning with Gamma Spectroscopy Data Augmentations for Detecting Shielded Radiological Material Transfers

Data analysis techniques can be powerful tools for rapidly analyzing data and extracting information that can be used in a latent space for categorizing observations between classes of data. Machine learning models that exploit learned data relationships can address a variety of nuclear nonproliferation challenges like the detection and tracking of shielded radiological material transfers. The high resource cost of manually labeling radiation spectra is a hindrance to the rapid analysis of data collected from persistent monitoring and to the adoption of supervised machine learning methods that require large volumes of curated training data. Instead, contrastive self-supervised learning on unlabeled spectra can enhance models that are built on limited labeled radiation datasets. This work demonstrates that contrastive machine learning is an effective technique for leveraging unlabeled data in detecting and characterizing nuclear material transfers demonstrated on radiation measurements collected at an Oak Ridge National Laboratory testbed, where sodium iodide detectors measure gamma radiation emitted by material transfers between the High Flux Isotope Reactor and the Radiochemical Engineering Development Center. Label-invariant data augmentations tailored for gamma radiation detection physics are used on unlabeled spectra to contrastively train an encoder, learning a complex, embedded state space with self-supervision. A linear classifier is then trained on a limited set of labeled data to distinguish transfer spectra between byproducts and tracked nuclear material using representations from the contrastively trained encoder. The optimized hyperparameter model achieves a balanced accuracy score of 80.30%. Any given model—that is, a trained encoder and classifier—shows preferential treatment for specific subclasses of transfer types. Regardless of the classifier complexity, a supervised classifier using contrastively trained representations achieves higher accuracy than using spectra when trained and tested on limited labeled data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

FTL: Transfer Learning Nonlinear Plasma Dynamic Transitions in Low Dimensional Embeddings (FTL) v1.0

Fusion Transfer Learning (FTL) model provides a new paradigm to study high-dimensional dynamical behaviors, such as those in fusion plasma systems. The knowledge transfer process leverages a pre-trained neural encoder-decoder network, initially trained on linear simulations, to effectively capture nonlinear dynamics. The low-dimensional embeddings extract the coherent structures of interest, while preserving the inherent dynamics of the complex system. Experimental results highlight FTL's capacity to capture transitional behaviors and dynamical features in plasma dynamics -- a task often challenging for conventional methods. The model developed in this study is generalizable and can be extended broadly through transfer learning to address various magnetohydrodynamics (MHD) modes.

Bai, Zhe↗

DOE OSTI2024

Sequential Kalman tuning of the t -preconditioned Crank-Nicolson algorithm: efficient, adaptive and gradient-free inference for Bayesian inverse problems

Ensemble Kalman Inversion (EKI) has been proposed as an efficient method for the approximate solution of Bayesian inverse problems with expensive forward models. However, when applied to the Bayesian inverse problem EKI is only exact in the regime of Gaussian target measures and linear forward models. Here, in this work we propose embedding EKI and Flow Annealed Kalman Inversion, its normalizing flow (NF) preconditioned variant, within a Bayesian annealing scheme as part of an adaptive implementation of the t-preconditioned Crank-Nicolson (tpCN) sampler. The tpCN sampler differs from standard pCN in that its proposal is reversible with respect to the multivariate t-distribution. The more flexible tail behaviour allows for better adaptation to sampling from non-Gaussian targets. Within our Sequential Kalman Tuning (SKT) adaptation scheme, EKI is used to initialize and precondition the tpCN sampler for each annealed target. The subsequent tpCN iterations ensure particles are correctly distributed according to each annealed target, avoiding the accumulation of errors that would otherwise impact EKI. We demonstrate the performance of SKT for tpCN on three challenging numerical benchmarks, showing significant improvements in the rate of convergence compared to adaptation within standard SMC with importance weighted resampling at each temperature level, and compared to similar adaptive implementations of standard pCN. The SKT scheme applied to tpCN offers an efficient, practical solution for solving the Bayesian inverse problem when gradients of the forward model are not available. Code implementing the SKT schemes for tpCN is available at https://github.com/RichardGrumitt/KalmanMC.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Ring Formation in Protoplanetary Disks Driven by an Eccentric Instability

We find that, under certain conditions, protoplanetary disks may spontaneously generate multiple, concentric gas rings without an embedded planet through an eccentric cooling instability. Using both linear theory and nonlinear hydrodynamics simulations, we show that a variety of background states may trap a slowly processing, one-armed spiral mode that becomes unstable when a gravitationally stable disk rapidly cools. The angular momentum required to excite this spiral comes at the expense of nonuniform mass transport that generically results in multiple rings. For example, one long-term hydrodynamics simulation exhibits four long-lived, axisymmetric gas rings. We verify the instability evolution and ring-formation mechanism from first principles with our linear theory, which shows remarkable agreement with the simulation results. Dust trapped in these rings may produce observable features consistent with observed disks. Additionally, direct detection of the eccentric gas motions may be possible when the instability saturates, and any residual eccentricity left over in the rings at later times may also provide direct observational evidence of this mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

TRUST Nonlinear Dynamics Testbed Assessment

The following assessment evaluates the efficacy of the control script for carrying out a linear signal to generate a linear mechanical response of the system: the TRUST nonlinear dynamics (TRUSTND) testbed. The hardware has three main components: the controller (NI PXIe 8861) embedded in the National Instruments chassis (NI PXIe 1092), the signal amplifier (The Modal Shop Linear Power Amplifier 2050E09) to amplify the output of a custom LabVIEW script, and a shaker (The Modal Shop Electrodynamic Exciter 2075E) where wave spring specimens live and are tested within the attached aluminum testbed (center mass). Different kinds of signals (single-tone, swept frequency and white noise) were sent through this hardware in the E-1 lab space at TA-53. First, the linearity of the setup was tested by specifying a sine wave in the control script and checking the quality of the oscillations in a wave spring (McMaster-Carr 9714K19 [1]). When passing a sine wave through the amplifier was succeeded, the frequency response of the wave spring was tested with white noise for comparison to finite-element model predictions. The testbed outputs a repeatable, linear response to the input excitation when used with the wave spring. These measurements reduce the uncertainty associated with the testbed and test procedure to better characterize uncertainty associated with material behavior (SX358 foam samples) in future tests

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗