Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase errors, capture effect”

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

Quantum criticality using a superconducting quantum processor

Quantum criticality emerges from the collective behavior of many interacting quantum particles, often at the transition between different phases of matter. It is one of the cornerstones of condensed matter physics, which we access on noisy intermediate-scale (NISQ) quantum devices by leveraging a dynamically driven phenomenon. Here we probe the critical properties of the one-dimensional quantum Ising model on a programmable superconducting quantum chip via a Kibble-Zurek process, obtain scaling laws, and estimate critical exponents despite inherent sources of errors on the hardware. In addition, we investigate how the improvement of NISQ computers (more qubits, less noise) will consolidate the computation of those universal physical properties. A one-parameter noise model captures the effect of imperfections and reproduces the experimental data. Its systematic study reveals that the noise, analogously to temperature, induces a new length scale in the system. We introduce and successfully verify modified scaling laws, directly accounting for the noise without any prior knowledge. It makes data analyses for extracting physical properties transparent to noise. By understanding how imperfect quantum hardware modifies the genuine properties of quantum states of matter, we enhance the power of NISQ processors considerably for addressing quantum criticality and potentially other phenomena and algorithms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Wave Prediction using X-Band Radar (Final Summary Report)

Model Predictive Control (MPC) is the best controls framework available to wave energy today. It allows for the maximization of energy yield from a WEC system while respecting various constraints and losses in the system. This type of constrained optimization is a critical capability for the techno-economic optimization of WEC devices and is difficult to do effectively with causal controls frameworks. MPC requires a future prediction of wave excitation forces, requiring wave prediction on a future time horizon of up to 30 seconds. The future horizon required is highly dependent on the WEC and PTO type available. To our knowledge no-one has been able to implement MPC on a WEC device at sea, due to the fact that phase-resolved wave-prediction is not a capability that has been sufficiently developed to date. The key objectives of this project are focused on developing “industry-ready” wave prediction technology building blocks that can be applied to a wide range of different wave energy conversion topologies. Our collaborative work with device developers has demonstrated that the level of improvement attainable for any particular device is heavily dependent on the device and PTO topology chosen. The range of annual average power capture improvements is on the order of 10% to over 100%. However, the important aspect is that we can get within about 10% of theoretical upper limits even when considering errors in the wave prediction and more importantly, we can optimize performance while respecting various constraints such as motion amplitude or peak structural loads. The combination of improved performance while limiting structural loads has a net effect of reducing the levelized cost of energy from these emerging technologies and allows us to design control laws that meet an economic optimum. Efforts under a previous project on controls has focused on developing a set of wave-prediction algorithms that work by leveraging a network of wave sensing measurement buoys that are equipped with real-time telemetry. The present work is focusing primarily on using X-band radar to measure the wave-field around the WEC to predict the future excitation forces on the structure and combining the measurements of these different sources to create improvements in the wave-prediction accuracy.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

A machine learning approach to quantify degradation of nuclear fuels and the effects of fission products

Nuclear fuel performance is critically dependent on understanding the evolution of fuel properties under operational conditions, a complex challenge driven by chemical changes and substantial radiation damage during fission. Traditionally, property evolution has been determined via empirical data collected following irradiation. However, these empirical correlations are limited in their applicability beyond the specific conditions in which they were obtained. This study explores a novel approach to address this challenge by applying materials informatics to develop a machine learning random forest (ML-RF) model that captures the effects of fission products on fuel compounds. The model predicts formation enthalpy (ΔH f ) by leveraging extensive quantum materials property data and correlating it with material descriptors such as composition, atomic and site features, and crystal lattice properties. This ML-RF model enables rapid interpolation across the compositional and structural spaces covered by the training data, thus supporting high-throughput screening and energetic ranking of candidate phases. The model demonstrates the ability to predict ΔH f with a mean absolute error (MAE) of approximately 0.1 to 0.2 eV/atom across a wide range of compounds, including key nuclear fuel systems (U-O, U-N, U-C, U-Si, and U-Mo). For example, it was used to assess shifts in stoichiometry for UO 2 (O/M) and UN (N/M) fuels, revealing their distinct tendencies in chemical potential variation and enabling preliminary convex hull analyses. Furthermore, the model provides insights into how individual fission products affect fuel properties. Results indicate that larger fission products (e.g., Nd, Pu, Ce) have a more pronounced impact on UO 2 , while lighter ones (e.g., Zr) strongly influence UN. Here, the model developed in this work can be used to support the Accelerated Fuel Qualification approach by facilitating preliminary evaluations prior to extensive materials modeling and experimentation. To this end, the trained model has been made available to the fuel community to support ongoing fuel development efforts.

Accelerated fuel qualification↗
DOE OSTI2019

Charge transfer excitation energies from ground state density functional theory calculations

Calculating charge transfer (CT) excitation energies with high accuracy and low computational cost is a challenging task. Kohn-Sham density functional theory (KS-DFT), due to its efficiency and accuracy, has achieved great success in describing ground state problems. To extend to excited state problems, our group recently demonstrated an approach with good numerical results to calculate low-lying and Rydberg excitation energies of an N-electron system from a ground state KS or generalized KS calculations of an (N-1)-electron system via its orbital energies. In the present work, we explore further the same methodology to describe CT excitations. Numerical results from this work show that performance of conventional density functional approximations (DFAs) is not as good for CT excitations as for other excitations due to the delocalization error. Applying localized orbital scaling correction (LOSC) to conventional DFAs, a recently developed method in our group to effectively reduce the delocalization error, can improve the results. Overall, the performance of this methodology is better than time dependent DFT (TDDFT) with conventional DFAs. In addition, it shows that results from LOSC-DFAs in this method reach similar accuracy to other methods, such as ΔSCF, G 0 W 0 with Bethe-Salpeter equations, particle-particle random phase approximation, and even high-level wavefunction methods like CC2. Our analysis shows that the correct 1/R trend for CT excitation can be captured from LOSC-DFA calculations, stressing that the application of DFAs with the minimal delocalization error is essential within this methodology. This work provides an efficient way to calculate CT excitation energies from ground state DFT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Fast Quasi-Static Time-Series (QSTS) for yearlong PV impact studies using vector quantization

The rapidly growing penetration levels of distributed photovoltaic (PV) systems requires more comprehensive studies to understand their impact on distribution feeders. IEEE P.1547 highlights the need for Quasi-Static Time Series (QSTS) simulation in conducting distribution impact studies for distributed resource interconnection. Unlike conventional scenario-based simulation, the time series simulation can realistically assess time-dependent impacts such as the operation of various controllable elements (e.g. voltage regulating tap changers) or impacts of power fluctuations. However, QSTS simulations are still not widely used in the industry because of the computational burden associated with running yearlong simulations at a 1-s granularity, which is needed to capture device controller effects responding to PV variability. This study presents a novel algorithm that reduces the number of times that the non-linear 3-phase unbalanced AC power flow must be solved by storing and reassigning power flow solutions as it progresses through the simulation. Each unique power flow solution is defined by a set of factors affecting the solution that can easily be queried. We demonstrate a computational time reduction of 98.9% for a yearlong simulation at 1-s resolution with minimal errors for metrics including: number of tap changes, capacitor actions, highest and lowest voltage on the feeder, line losses, and ANSI voltage violations. Finally, the key contribution of this work is the formulation of an algorithm capable of: (i) drastically reducing the computational time of QSTS simulations, (ii) accurately modeling distribution system voltage-control elements with hysteresis, and (iii) efficiently compressing result time series data for post-simulation analysis.

14 SOLAR ENERGY↗
DOE OSTI2016

SU-G-IeP4-09: Method of Human Eye Aberration Measurement Using Plenoptic Camera Over Large Field of View

Purpose: The measurement based on Shack-Hartmann wave-front sensor(WFS), obtaining both the high and low order wave-front aberrations simultaneously and accurately, has been applied in the detection of human eyes aberration in recent years. However, Its application is limited by the small field of view (FOV), slight eye movement leads the optical bacon image exceeds the lenslet array which result in uncertain detection error. To overcome difficulties of precise eye location, the capacity of detecting eye wave-front aberration over FOV much larger than simply a single conjugate Hartmann WFS accurately and simultaneously is demanded. Methods: Plenoptic camera’s lenslet array subdivides the aperture light-field in spatial frequency domain, capture the 4-D light-field information. Data recorded by plenoptic cameras can be used to extract the wave-front phases associated to the eyes aberration. The corresponding theoretical model and simulation system is built up in this article to discuss wave-front measurement performance when utilizing plenoptic camera as wave-front sensor. Results: The simulation results indicate that the plenoptic wave-front method can obtain both the high and low order eyes wave-front aberration with the same accuracy as conventional system in single visual angle detectionand over FOV much larger than simply a single conjugate Hartmann systems. Meanwhile, simulation results show that detection of eye aberrations wave-front in different visual angle can be achieved effectively and simultaneously by plenoptic method, by both point and extended optical beacon from the eye. Conclusion: Plenoptic wave-front method possesses the feasibility in eye aberrations wave-front detection. With larger FOV, the method can effectively reduce the detection error brought by imprecise eye location and simplify the eye aberrations wave-front detection system comparing with which based on Shack-Hartmann WFS. Unique advantage of the plenoptic method lies in obtaining wave-front in different visual angle simultaneously, which provides an approach in building up 3-D model of eye refractor tomographically. Funded by the key Laboratory of High Power Laser and Physics, CAS Research Project of National University of Defense Technology No. JC13-07-01; National Natural Science Foundation of China No. 61205144.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2013

Covariance Analysis Tool (G-CAT) for Computing Ascent, Descent, and Landing Errors

G-CAT is a covariance analysis tool that enables fast and accurate computation of error ellipses for descent, landing, ascent, and rendezvous scenarios, and quantifies knowledge error contributions needed for error budgeting purposes. Because GCAT supports hardware/system trade studies in spacecraft and mission design, it is useful in both early and late mission/ proposal phases where Monte Carlo simulation capability is not mature, Monte Carlo simulation takes too long to run, and/or there is a need to perform multiple parametric system design trades that would require an unwieldy number of Monte Carlo runs. G-CAT is formulated as a variable-order square-root linearized Kalman filter (LKF), typically using over 120 filter states. An important property of G-CAT is that it is based on a 6-DOF (degrees of freedom) formulation that completely captures the combined effects of both attitude and translation errors on the propagated trajectories. This ensures its accuracy for guidance, navigation, and control (GN&C) analysis. G-CAT provides the desired fast turnaround analysis needed for error budgeting in support of mission concept formulations, design trade studies, and proposal development efforts. The main usefulness of a covariance analysis tool such as G-CAT is its ability to calculate the performance envelope directly from a single run. This is in sharp contrast to running thousands of simulations to obtain similar information using Monte Carlo methods. It does this by propagating the "statistics" of the overall design, rather than simulating individual trajectories. G-CAT supports applications to lunar, planetary, and small body missions. It characterizes onboard knowledge propagation errors associated with inertial measurement unit (IMU) errors (gyro and accelerometer), gravity errors/dispersions (spherical harmonics, masscons), and radar errors (multiple altimeter beams, multiple Doppler velocimeter beams). G-CAT is a standalone MATLAB- based tool intended to run on any engineer's desktop computer.

Boussalis, Dhemetrios↗
NASA NTRS1994

Development of a Model Following Control Law for Inflight Simulation and Flight Controls Research

The U.S. Army and NASA are currently developing the Rotorcraft Aircrew Systems Concepts Airborne Laboratory (RASCAL) at the Ames Research Center. RASCAL, shown in Figure 1, is a UH-60, which is being modified in a phased development program to have a research fly-by-wire flight control system, and an advanced navigation research platform. An important part of the flight controls and handling qualities research on RASCAL will be an FCS design for the aircraft to achieve high bandwidth control responses and disturbance rejection characteristics. Initially, body states will be used as feedbacks, but research into the use of rotor states will also be considered in later stages to maximize agility and maneuverability. In addition to supporting flight controls research, this FCS design will serve as the inflight simulation control law to support basic handling qualities, guidance, and displays research. Research in high bandwidth controls laws is motivated by the desire to improve the handling qualities in aggressive maneuvering and in severely degraded weather conditions. Naturally, these advantages will also improve the quality of the model following, thereby improving the inflight simulation capabilities of the research vehicle. High bandwidth in the control laws provides tighter tracking allowing for higher response bandwidths which can meet handling qualities requirements for aggressive maneuvering. System sensitivity is also reduced preventing variations in the response from the vehicle due to changing flight conditions. In addition, improved gust rejection will result from this reduced sensitivity. The gust rejection coupled with a highly stable system will make more precise maneuvering and pointing possible in severely degraded weather conditions. The difficulty in achieving higher bandwidths from the control laws in the feedback and in the responses arises from the complexity of the models that are needed to produce a satisfactory design. In this case, high quality models that include rotor dynamics in a physically meaningful context must be available. A non-physical accounting of the rotor, such as lumping the effect as a time delay, is not likely to produce the desired results. High order simulation models based on first principals are satisfactory for the initial design phase in order to work out the control law design concept and get an initial set of gains. These models, however, have known deficiencies, which must be resolved in the final control law design. The error in the pitch-roll cross coupling is one notable deficiency that even sophisticated rotorcraft models including complex wake aerodynamics have yet to capture successfully. This error must be accounted for to achieve the desired decoupling. The approach to design the proposed inflight simulation control law is based on using a combination of simulation and identified models. The linear and nonlinear higher order models were used to develop an explicit model following control structure. This structure was developed to accommodate the design of control laws compliant to many of the quantitative requirements in ADS-33C. Furthermore, it also allows for control law research using rotor-state feedback and other design methodologies such as Quantitative Feedback and H-Infinity. Final gain selection will be based on higher order identified models which include rotor degrees of freedom.

Takahashi, Mark↗
NASA NTRS2002

Segmented mirror coarse phasing with a dispersed fringe sensor: experiments on NGST's wavefront control testbed

A piston sensing and control algorithm for the segmented mirror coarse phasing using a dispersed fringe sensor (DFS) has been developed for Next Generation Space Telescope (NGST) wavefront sensing and control. The DFS can detect residual piston errors as large as order of a depth-of-focus and can phase the segment mirrors with an accuracy less than 0.1 micron, which is within the capture range of the fine phasing for NGST. A series of experiments has been carried out on the NGST's Wavefront Control Testbed (WCT) to validate the modeling results, evaluate the DFS performance, and systematically explore the factors that affect the DFS performance. This paper reports the testbed results for several critical issues of DFS performance, including DFS dynamic range, accuracy, fringe visibility, and the effects of segment mirror aberrations.

Space↗
DOE OSTI2017

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Highlights: • Two numerical methods are presented to accelerate DMD simulations. • The methods greatly reduce the computation time without much loss of accuracy. • The accelerated DMD model is capable of simulating long-term mass transport at atomistic length scales. • Demonstrated for the propagation of an atomically sharp phase boundary over 30 seconds. - Abstract: Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Late-time near-infrared observations of SN 2005df

We present late-time near-infrared (NIR) spectral evolution, at 200–400 days, for the Type Ia supernova SN 2005df. The spectra show numerous strong emission features of [Co ii], [Co iii], and [Fe ii] throughout the 0.8–1.8 μm region. As the spectrum ages, the cobalt features fade as would be expected from the decay of {sup 56}Co to {sup 56}Fe. We show that the strong and isolated [Fe ii] emission line at 1.644 μm provides a unique tool to analyze NIR spectra of SNe Ia. Normalization of spectra to this line allows the separation of features produced by stable versus unstable isotopes of iron group elements. We develop a new method of determining the initial central density, ρ{sub c}, and the magnetic field, B, of the white dwarf (WD) using the width of the 1.644 μm line. The line width (LW) is sensitive because of electron capture in the early stages of burning, which increases as a function of density. The sensitivity of the LW to B increases with time, and the effects of the magnetic field shift toward later times with decreasing ρ{sub c}. Through comparison with spherical models, the initial central density for SN 2005df is measured as ρ{sub c}=0.9(±0.2)×10{sup 9} g cm{sup −3}, which corresponds to a WD close to the Chandrasekhar mass, with M{sub WD}=1.31(±0.03) M{sub ⊙} and systematic error less than 0.04 M{sub ⊙}. This error estimate is based on spherical models. We discuss the potential uncertainties due to multi-dimensional effects, mixing, and rotation. The latter two effects would increase the estimate of the WD mass. Within M{sub Ch} explosions, however, the central density found for SN 2005df is very low for a H-accretor, possibly suggesting a helium star companion or a tidally disrupted WD companion. As an alternative, we suggest mixing of the central region. We find some support for high initial magnetic fields of strength 10{sup 6} G for SN 2005df, however, 0 G cannot be ruled out because of noise in the spectra combined with low ρ{sub c}. We discuss our findings in the context of mixing by Rayleigh–Taylor instabilities during deflagration burning and a wide variety of explosion scenarios. Observations strongly support a very limited amount of mixing during a deflagration phase and high central densities characteristic of a M{sub Ch} WD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Modeling Analysis of Ball-Milling Process for Battery-Electrode Synthesis

The mechanical alloying process is a promising method for synthesizing electrode materials for batteries owing to its benefits such as the ability to produce nanostructured, high-performing electrode alloys, no adverse effects on the solid electrolyte for solid-state batteries, stable production of thick electrodes, simple processing steps, and low processing costs. It is gaining intensive attention in the battery industry as one of the best methods to replace the conventional wet-slurry-solvent method, and its application is rapidly increasing these days. However, the operation is currently conducted purely based on trial-and-error methods without fully utilizing the features of its functions. Here, this may be attributed to a lack of understanding of the effect of operating parameters on the alloying process and final products. Surprisingly, there is a scarcity of the literature conducting fundamental research to comprehend the underlying physics of the entire mechanical alloying process, resulting in a significant knowledge gap. To address this knowledge gap, extensive research was conducted. The existing literature on mechanical alloying was reviewed to comprehend the current state of understanding and to discuss the direction for future research. Mathematical expressions were developed to create physics-based models capable of capturing the entire mechanical alloying process, including milling kinetics and defect-enhanced phase evolution. These methods were then applied to investigate the impact of operating parameters such as milling frequency, initial mole ratio of the alloyed materials, density of grinding balls, and energy required for the powders to become amorphous (i.e., the amorphization energy threshold). This research aimed not only to comprehend the direct effects of these operating parameters but also to unveil the physics underlying the ball-milling process. The results of our study can serve as crucial information for the battery industry in designing or operating the ball-milling process.

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

Comparison of Distributive Consensus Algorithms for AnomalyDetection on the Power-Grid

As an increasing number of devices in the power-grid rely on high-precision clocks to measure electrical properties and time stamp data, GPS receivers have become a vulnerability on the power-grid. GPS spoofing attacks, in which a fake signal is transmitted, can cause sensor measurement errors that can misinform protective equipment and grid operators for both the transmission and distribution grid. The technology being developed to secure precision clocks in the transmission grid are not appropriate for the distribution grid. Thus, there is a need to developed anomaly detection algorithms that are suitable for the hardware available on the distribution grid. Traditionally, a distributed detection framework is comprised of spatially distributed nodes which acquire observations about a quantity of interest and sends the data to the fusion center, where a global decision is made. Instead the distribution grid has low-power devices, such as smart inverters and microPMUs, that could be used to perform highly parallelized calculations. Having on node calculations is beneficial since it is more network and computation efficient, as the more calculations can be performed where the data is being captured instead of communicating across a slow network to a centralized fusion center. Instead, we could employ an alternate peer-to-peer type of information exchange on the distribution grid that can still reach a global decision. A decentralized approach to anomaly detection makes efficient use of the computational and network bandwidth that exists, and it also lacks the single point of failure of a centralized approach. Thus, a decentralized local information exchange and/or inference is preferred. Lastly, the unreliability of the distribution grid network connections and potential noise or error in the measurements must be considered. In order to produce an effective detection tool, a Byzantine-fault tolerant and asynchronous algorithm must be developed. Thus, the first phase of our project is to develop a distributed, decentralized, asynchronous, Byzantine-fault tolerant, collaborative autonomy-based anomaly detection algorithm to collectively analyze clusters of data for signs of a malicious attack.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Improving the MJO Forecast of S2S Operation Models by Correcting Their Biases in Linear Dynamics

The operational dynamic subseasonal to seasonal (S2S) models for Madden-Julian oscillation (MJO) forecasting mostly still suffer from systematic errors in capturing the MJO's key dynamic features, such as its growth rate and propagation speed. By deriving the linear dynamic operators using the linear inverse modeling (LIM) approach, we propose a method to partly correct the errors in MJO linear dynamic operators to improve the MJO predictions of three operational dynamic S2S models. Correcting the deficiencies of the too-fast decay rates and the unrealistic propagating phase speeds lead to MJO prediction skills being extended by approximately 2–4 days. The improvements are more significant for the models with larger biases in MJO amplitude and propagation. This approach in principle may be extendable to predictions of other types of climate variability such as ENSO on one hand, and possible inclusions of nonlinear dynamics effects on the other hand.

58 GEOSCIENCES↗
DOE OSTI2018

Closed‐loop hierarchical control strategies for connected and autonomous hybrid electric vehicles with random errors

This study presents a closed‐loop hierarchical control strategy for a group of connected and autonomous hybrid electric vehicles (HEVs) with the purpose of optimising fuel economy while guaranteeing traffic mobility and vehicle safety. In the hierarchical control architecture, a decentralised stochastic model predictive control (SMPC)‐based higher level controller incorporating signal phase and timing information is formulated to optimise the velocity profile of each vehicle, and an adaptive equivalent consumption minimisation strategy based lower level controller is employed for the energy management control of the HEVs. Creatively, random errors of the control variable are considered in the SMPC framework. The errors are discretised and modelled as a Markov process and the state transition matrix of the errors is generated randomly to capture the error transition dynamics. To solve the SMPC problem more efficiently, a scenario‐based SMPC is employed. Moreover, the propulsion and recuperation efficiencies of the lower level controller are calculated at each time step with measurable variables, and fed back to the higher level controller for velocity optimisation in next time step. Simulation results validate the control effectiveness and advantages of the proposed control architecture.

Qiu, Shaolin↗
DOE OSTI2019

Auxiliary field quantum Monte Carlo for multiband Hubbard models: Controlling the sign and phase problems to capture Hund's physics

In the study of strongly correlated, many-electron systems, the Hubbard-Kanamori (HK) model has emerged as one of the prototypes for transition metal oxide physics. The model is multiband in nature and contains Hund's coupling terms, which have pronounced effects on metal-insulator transitions, high-temperature superconductivity, and other physical properties. In the following, we present a complete theoretical framework for treating the HK model using the ground-state auxiliary field quantum Monte Carlo (AFQMC) method and analyze its performance on few-band models whose parameters approximate those observed in ruthenates, rhodates, and other materials exhibiting Hund's physics. Unlike previous studies, the constrained path and phaseless approximations are used to respectively control the sign and phase problems, which enables high-accuracy modeling of the HK model's ground-state properties within parameter regimes of experimental interest. We demonstrate that, after careful consideration of the Hubbard-Stratonovich transformations and trial wave functions employed, relative errors in the energy of less than 1% can routinely be achieved for moderate to large values of the Hund's coupling constant. Crucially, our methodology also accurately predicts magnetic ordering and phase transitions. Furthermore, the results presented open the door to more predictive modeling of Hund's physics within a wide range of strongly correlated materials using AFQMC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Temperature dependence of magnetic anisotropy and magnetoelasticity from classical spin-lattice calculations

Here we present a classical molecular-spin dynamics (MSD) methodology that enables accurate computations of the temperature dependence of the magnetocrystalline anisotropy as well as magnetoelastic properties of magnetic materials. The nonmagnetic interactions are accounted for by a spectral neighbor analysis potential (SNAP) machine-learned interatomic potential, whereas the magnetoelastic contributions are accounted for using a combination of an extended Heisenberg Hamiltonian and a Néel pair interaction model, representing both the exchange interaction and spin-orbit-coupling effects, respectively. All magnetoelastic potential components are parameterized using a combination of first-principles and experimental data. Our framework is applied to the α phase of iron. Initial testing of our MSD model is done using a 0 K parametrization of the Néel interaction model. After this, we examine how individual Néel parameters impact the $B$ 1 and $B$ 2 magnetostrictive coefficients using a moment-independent δ sensitivity analysis. The results from this study are then used to initialize a genetic algorithm optimization which explores the Néel parameter phase space and tries to minimize the error in the B 1 and B 2 magnetostrictive coefficients in the range of 0–1200 K. Our results show that while both the 0 K and genetic algorithm optimized parametrization provide good experimental agreement for $B$ 1 and $B$ 2 , only the genetic algorithm optimized results can capture the second peak in the $B$ 1 magnetostrictive coefficient which occurs near approximately 800 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗