Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetization Process of Atacamite: A Case of Weakly Coupled S = 1/2 Sawtooth Chains

We present a combined experimental and theoretical study of the mineral atacamite Cu 2 Cl(OH) 3 . Density-functional theory yields a Hamiltonian describing anisotropic sawtooth chains with weak 3D connections. Experimentally, we fully characterize the antiferromagnetically ordered state. Here, magnetic order shows a complex evolution with the magnetic field, while, starting at 31.5 T, we observe a plateaulike magnetization at about M sat /2. Based on complementary theoretical approaches, we show that the latter is unrelated to the known magnetization plateau of a sawtooth chain. Instead, we provide evidence that the magnetization process in atacamite is a field-driven canting of a 3D network of weakly coupled sawtooth chains that form giant moments.

1-dimensional spin chains↗
DOE OSTI2021

Spectral Denoising for Accelerated Analysis of Correlated Ionic Transport

Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Dual Mechanism for Transient Capacitance Anomaly in Improper Ferroelectrics

Negative capacitance (NC) effects in ferroelectrics can potentially break fundamental limits of power dissipation known as “Boltzmann tyranny.” However, the origin of transient NC of ferroelectrics, which is attributed to two different mechanisms involving free-energy landscape and nucleation, is under intense debate. Here, we report the coexistence of transient NC and an S-shaped anomaly during the switching of ferroelectric hexagonal ferrites capacitor in an RC circuit. The early-stage NC arises from the nucleation process, while the late-stage S-shaped anomaly corresponds to a nascent NC associated with the free-energy landscape. The entire waveform can be reproduced using a hybrid model that simultaneously incorporates these two mechanisms. These results highlight the multivariable free-energy landscape of hexagonal ferrites that enables an abrupt change of the internal field and demonstrate that the two mechanisms are not mutually exclusive, resolving the long-standing debate. In conclusion, the behavior of the S-shaped anomaly also provides a pathway to extract parameters of free-energy landscape and switching dynamics.

Capacitance↗
DOE OSTI2022

Transitions in the morphology and critical stresses of gliding dislocations in multiprincipal element alloys

Refractory multiprincipal element alloys (MPEAs) are promising material candidates for high-temperature, high-strength applications. However, their deformation mechanisms, particularly at the dislocation scale, are unusual compared to those of conventional alloys, making it challenging to understand the origin of their high strength. Here, using an atomistically informed phase-field dislocation dynamics model, we study transitions in the morphology and critical stresses of long gliding screw dislocations over extended distances in a refractory MPEA. The model MPEA crystal accounts for the atomic-scale fluctuations in chemical composition across the glide planes via spatially correlated lattice energies and for the differences in glide resistances between screw and edge dislocations of unit length. Here, we show that the dislocation moves in a stop-start motion, alternating between wavy morphology in free flight and nearly recovered straight screw orientation in full arrest. The periods of wavy glide are due to variable kink-pair formation and migration rates along the length, where portions with higher rates glide more quickly. The critical stress to initiate motion corresponds to the stress required to form and migrate a kink pair at the weakest region along the length of the dislocation. Heterogeneity in lattice energy leads to variability in the local stress-strain response and to a strain hardening-like response, in which the critical stress to reactivate glide increases with glide distance. Statistical assessment of hundreds of realizations of dislocations indicates that the amount of hardening directly scales with the dispersion in underlying lattice energy.

3-dimensional systems↗
DOE OSTI2025

Quantifying the impact of precision errors on quantum approximate optimization algorithms

The quantum approximate optimization algorithm (QAOA) is a hybrid quantum-classical algorithm that seeks to achieve approximate solutions to optimization problems by iteratively alternating between intervals of controlled quantum evolution. Here, we examine the effect of analog precision errors on QAOA performance from the perspective of both algorithmic training and performance guarantees. Leveraging cumulant expansions, we recast the faulty QAOA as a control problem in which precision errors are expressed as multiplicative control noise and derive bounds on the performance of QAOA. We show using both analytical techniques and numerical simulations that fixed precision implementations of QAOA circuits are subject to an exponential degradation in performance dependent upon the number of optimal QAOA layers and magnitude of the precision error. Despite this significant reduction, we show that it is possible to mitigate precision errors in QAOA via digitization of the variational parameters at the cost of increasing circuit depth.

quantum algorithms↗
DOE OSTI2021

Many-Body Physics in the NISQ Era: Quantum Programming a Discrete Time Crystal

Recent progress in the realm of noisy intermediate-scale quantum (NISQ) devices represents an exciting opportunity for many-body physics by introducing new laboratory platforms with unprecedented control and measurement capabilities. We explore the implications of NISQ platforms for many-body physics in a practical sense: we ask which physical phenomena, in the domain of quantum statistical mechanics, they may realize more readily than traditional experimental platforms. While a universal quantum computer can simulate any system, the eponymous noise inherent to NISQ devices practically favors certain simulation tasks over others in the near term. As a particularly well-suited target, we identify discrete time crystals (DTCs), novel nonequilibrium states of matter that break time translation symmetry. These can only be realized in the intrinsically out-of-equilibrium setting of periodically driven quantum systems stabilized by disorder-induced many-body localization. While promising precursors of the DTC have been observed across a variety of experimental platforms—ranging from trapped ions to nitrogen-vacancy centers to NMR crystals—none have all the necessary ingredients for realizing a fully fledged incarnation of this phase, and for detecting its signature long-range spatiotemporal order. We show that a new generation of quantum simulators can be programmed to realize the DTC phase and to experimentally detect its dynamical properties, a task requiring extensive capabilities for programmability, initialization, and readout. Specifically, the architecture of Google’s Sycamore processor is a remarkably close match for the task at hand. We also discuss the effects of environmental decoherence, and how they can be distinguished from ‘internal’ decoherence coming from closed-system thermalization dynamics. Already with existing technology and noise levels, we find that DTC spatiotemporal order would be observable over hundreds of periods, with parametric improvements to come as the hardware advances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Fokker-Planck Equation Governing the Distribution of Walkers in Auxiliary-Field Quantum Monte Carlo

Auxiliary-field quantum Monte Carlo (AFQMC) is typically formulated as an open-ended random walk in an overcomplete space of Slater determinants, implemented through a Langevin equation. However, the explicit form of the underlying Fokker-Planck equation governing the walker population distribution has remained unknown. Here, in this Letter, we derive the Fokker-Planck equation for AFQMC and propose a novel numerical scheme to solve it. The solution of the Fokker-Planck equation reveals the wave function actually sampled by the AFQMC algorithm. Interestingly, we find that even when the exact ground state is used as a guiding wave function in constrained path AFQMC, contrary to the common assumption, the wave function sampled by AFQMC is not exact. Beyond clarifying several fundamental aspects of AFQMC, the availability of a Fokker-Planck equation formulation opens new avenues for systematically improving its accuracy, which we outline in this Letter.

Monte Carlo methods↗
DOE OSTI2023

Factorization Machine Learning for Disaggregation of Transmission Load Profiles with High Penetration of Behind-the-Meter Solar

The ever-growing high penetration of ubiquitously distributed energy resources, especially behind-the-meter solar (BTM) generations, has significant impacts on nodal load (i.e., net injection) profiles and consequently caused imperative operational challenges to system operators such as regional transmission organizations (RTOs). Illustrated by real-world nodal data and examples at PJM Interconnection, this paper first discusses the application and necessity of effectively extracting daily nodal load profiles in a non-intrusive manner. More importantly, a novel bi-level architecture, including Factorization Machines (FM) learning procedure has been proposed to effectively disaggregate not only one node but every node in an RTO service territory. Specifically, FM leaning is adopted to capture the interconnections between related features to better utilize the correlation between buses in the same region and between a single bus and the zonal load. The proposed bi-level technique is numerically validated using real-world, minute-level, normalized, and anonymized nodal data at PJM service territory.

behind the meter solar, load disaggregation, load ↗
DOE OSTI2020

Computationally efficient Bayesian quantum state tomography

We describe a method for Bayesian quantum state estimation combining efficient parameterization, a pseudo-likelihood, and advanced numerical sampling techniques. Examples reveal significant computational speedup, indicating the approach’s promise in practical quantum state tomography.

Lukens, Joseph↗
DOE OSTI2021

PMU-Based Decoupled State Estimation for Unsymmetrical Power Systems

Modal decomposition of measurement equations has already been shown to simplify the formulation and resulting computational complexity of three-phase state estimation of systems where all the transmission lines are three-phase and fully transposed. When there are non-transposed and/or mixed-phase lines, modal decomposition can no longer fully decouple the threephase measurement equations. Here, this paper addresses the above shortcoming by proposing a simple yet practical solution based on the commonly used numerical compensation techniques. Thus, it enables application of the powerful decoupling approach to any type of three-phase networks which may contain non-transposed or mixed-phase lines and are fully observable by PMUs. The proposed procedure modifies the measurement set by deriving additive terms that compensate for the neglected unsymmetrical effects. It will be shown that unbalanced systems including nontransposed and mixed-phase elements, can still be transformed into three decoupled subsystems and solved in parallel by the proposed approach. Performance of the proposed algorithm is validated against several IEEE test cases.

42 ENGINEERING↗
DOE OSTI2023

Insights into Plastic Localization by Crystallographic Slip from Emerging Experimental and Numerical Approaches

Advanced experimental and numerical approaches are being developed to capture the localization of plasticity at the nanometer scale as a function of the multiscale and heterogeneous microstructure present in metallic materials. These innovative approaches promise new avenues to understand microstructural effects on mechanical properties, accelerate alloy design, and enable more accurate mechanical property prediction. This article provides an overview of emerging approaches with a focus on the localization of plasticity by crystallographic slip. New insights into the mechanisms and mechanics of strain localization are addressed. The consequences of the localization of plasticity by deformation slip for mechanical properties of metallic materials are also detailed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Vibration-Based Sensor Design: A Grey-Box Approach

Knowledge of the internal structure of an object or device under investigation proceeds from the basic idea of constructing its dynamic behavioral relations governed by a set of differential/algebraic equations that characterize its response. These equations can be partial differential equations leading to finite element or finite difference relations requiring a complex numerical solution on a super computer or ordinary differential equations requiring sophisticated numerical integration techniques to obtain the desired solution. Discrete dynamic systems evolving from digitized data acquisition are typically captured by sampled-data (continuous-to-discrete) representations characterized by a set of difference equations specifying the underlying system dynamics. In any case, with a mathematical description in hand, Grey-Box modeling techniques have evolved, concerned with the estimation of model parameters embedded in a prescribed set of equations (the system) governing its behavior, while capturing the underlying physical phenomenology of the problem at hand.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fabrication of Biofunctionalized PVA/PEI Nanofibers for Rare Earth Element Extraction from Coal Mine Tailings

The need for rare earth elements (REEs) for new technology is becoming increasingly stiff with demand at a record high. Biological ligands like lanmodulin and lanthanide binding tags (LBT) have gained noticeable attention for their selective extraction of valuable lanthanide metals in aqueous solution. Building on previous work using microbeads as the extraction media, we propose that nanofibers can improve extraction capabilities. Nanofibers could provide more efficient flow and higher extraction yields due to their high surface area to volume ratio. Relying on thiolmaleimide click chemistry as the biological ligand immobilization route, electrospinning as the nanofiber manufacturing method, and numerous characterization techniques, we produced electrospinning polymer solutions and nanofibers with a maximum neodymium adsorbent capacity of 12.1 μmol/g of extraction media immobilized with LBT. Although the prior work produced an extraction media with an equivalent capacity of ~40 μmol/g, the electrospinning process is much more economical and scalable. This work paves the way for an industrially relevant, environmentally-friendly REE extraction material. Further work is required to produce higher quality nanofibers, to investigate common ion interference with REE adsorbance, and to scale production up to a relevant scale.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Ensemble Simulation Techniques and Fast Randomized Algorithms

The major goals of the project were to develop and analyze new ensemble simulation techniques, including trajectory stratification and preconditioned MCMC techniques, as well as develop fast numerical linear algebra techniques closely related to ensemble simulation ideas. The trajectory stratification techniques involve simulating in parallel short trajectory fragments of a Markov process confined to a specific region of space‐time and then patching together the statistics gathered to assemble estimates of very general dynamical properties. We have also developed this approach for rare event simulation and extended the techniques to applications requiring a more general framework (such as electronic structure calculations). The preconditioned MCMC techniques involve simulating multiple Markov chains in parallel and then using information from the ensemble to speed the mixing of each individual chain. The fast randomized linear algebra methods are motivated by the diffusion Monte Carlo technique, but are applicable to finding the dominant eigenvalue of (almost) general matrices. For most non‐negative matrices, the schemes result in an error (compared to the power method) that is constant in the dimension of the problem. For more general matrices, we see a very clear sublinear cost trend in computational tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Sample-efficient verification of continuously-parameterized quantum gates for small quantum processors

Most near-term quantum information processing devices will not be capable of implementing quantum error correction and the associated logical quantum gate set. Instead, quantum circuits will be implemented directly using the physical native gate set of the device. These native gates often have a parameterization (e.g., rotation angles) which provide the ability to perform a continuous range of operations. Verification of the correct operation of these gates across the allowable range of parameters is important for gaining confidence in the reliability of these devices. In this work, we demonstrate a procedure for sample-efficient verification of continuously-parameterized quantum gates for small quantum processors of up to approximately 10 qubits. This procedure involves generating random sequences of randomly-parameterized layers of gates chosen from the native gate set of the device, and then stochastically compiling an approximate inverse to this sequence such that executing the full sequence on the device should leave the system near its initial state. We show that fidelity estimates made via this technique have a lower variance than fidelity estimates made via cross-entropy benchmarking. This provides an experimentally-relevant advantage in sample efficiency when estimating the fidelity loss to some desired precision. We describe the experimental realization of this technique using continuously-parameterized quantum gate sets on a trapped-ion quantum processor from Sandia QSCOUT and a superconducting quantum processor from IBM Q, and we demonstrate the sample efficiency advantage of this technique both numerically and experimentally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Pressure Exchanger for Energy Recovery in a Trans-Critical CO2 Refrigeration System

Trans-critical CO2 vapor compression (VC) refrigeration cycles require a high compression ratio, which is associated with high expansion losses. To recover these expansion losses, a pressure exchange process between the low- and high-pressure sides of the VC cycle is proposed and examined in this study. The proposed pressure exchange system is an open type constant volume process where the high- and low-pressure flows mix inside the system. This prototype is inspired by the pressure exchangers used in reverse-osmosis (RO) desalination systems. In this system, a 2D model was generated and modeled using the computational fluid dynamics (CFD) technique. The numerical model ignored any losses due to leakage or hydraulic friction and the process is considered adiabatic. For the modeling, it was assumed that the inlet conditions for the two pressure exchanger flows are similar to the flow conditions at the evaporator and gas cooler outlets in a VC cycle. Two parameters are examined to test the validity of the system and understand their effect on the performance, including the inlet flow rate represented by the inlet velocity and the process time represented by the speed of rotation. A total of nine cases were simulated and analyzed in this study.

30 DIRECT ENERGY CONVERSION↗