Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical simulation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗
DOE OSTI2026

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

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

Pore-scale visualization of natural hydrate-bearing sediments

Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.

Liu, Mengwei [Oak Ridge Institute for Science and ↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Project Title: Demonstration High Temperature Superconducting NonPlanar Stellarator Magnet with Advanced Manufactured Assemblies

This is the final report for the project “Demonstration High Temperature Superconducting Non- Planar Stellarator Magnet with Advanced Manufactured Assemblies”, funded by DOE, and performed by Type One Energy from September, 2020, to March 2024 involving the Fusion Technology Institute at the University of Wisconsin–Madison, the Plasma Science and Fusion Center (PSFC) at the Massachusetts Institute of Technology (MIT) and Commonwealth Fusion Systems (CFS) to design and fabricate the first non-planar HTS (REBCO) coil for a high-field stellarator based on the SPARC tokamak’s VIPER cable concept. Stellarators at high fields make high-temperature superconducting magnets necessary for a compact fusion device. But the asymmetric and non-planar nature of its components, especially the magnets make it difficult for scalable producibility. To address these challenges, two promising technologies have emerged: advanced manufacturing (AM) for the supporting plates for forming the magnets, and high-temperature superconducting (HTS) cables inside the plates. AM has advanced enough to produce stellarator components with the necessary geometric complexity, size, and the precision, leading to potentially significant reductions in production time, cost, and waste. The cost of HTS tape has decreased dramatically, and progress in HTS planar magnet development has reached a point where it can be proposed for application to complex 3D non-planar magnets. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, non-planar HTS coil shape that remains superconducting to achieve production scalable reductions in time and cost and performance. The proposed technology is based on the novel concept of a precision sub-scale HTS nonplanar coil assembly. This project focuses on the design, fabrication, material optimization of cable design, and validation and demonstration of the high current carrying capability of superconducting magnets and their support in a complex 3D shape needed for application to stellarator magnetic plasma confinement. The specific objectives of this research program include: (1) The successful application of metal AM to build a precision sub-scale HTS nonplanar coil, (2) An HTS cable and cross-section design that can conform to the required nonplanar coil shape (bend radii as tight as 10-cm) and remains superconducting at an engineering current density of 1.35 kA/cm 2 at 77 K and 1 tesla at the conductor (5 kA in the cable). To achieve the above challenging goals, we have formed a multidisciplinary research team consisting of members from Type One Energy and UW-Madison, MIT PSFC and CFS with complementary skills and strong facilities. The team worked collaboratively on fundamental and applied research on the following three major technical areas: (1) Design, fabrication, and optimization of non-planar HTS Cable The ultimate goal of the project is to determine if commercial REBCO tapes and additive manufacturing can be used to fabricate high field (≥ 10T) non-planar coils with tight bending radii (≃ 100mm) and with a degradation of the critical current (Ic) smaller than 20% with respect to the expected performance. We started with shorter length cable to evaluate the scalability of the production process and eventually reached multiple turns for higher magnetic fields. Our findings suggest that a stellarator coil system of a relevant size, characterized by its asymmetric and non-planar components, can be fabricated using a formed cable in plate method. This system can be simulated using a large-scale modeling approach. The use of hybrid modeling 3 techniques will be pivotal in reducing the complexity of the model and in assessing expected performance in designs. (2) Modeling and simulation of the non-planar HTS Cable Multiphysics simulations are performed using the commercial software and are carried out in self-field conditions, involving 2D and 3D models and twisted around one slot of twist-pitched VIPER cable. Multiphysics simulations are mainly focused on the coil for the critical current evaluation, the magnetic field map, self-Lorentz forces and mechanical, and magnetothermal behavior and the quench dynamics. The detailed model and prediction of the superconducting performance of a stellarator-relevant demonstration cable from numerical simulations supports the results from the actual testing backing the results. A detailed description and results are provided in the later sections. (3) Design, fabrication, and optimization of support for the non-planar HTS Cable The team developed an additive manufactured (AM) coil positioning plate that formed into the required non-planar geometry (with bend radii as tight as 10-cm) and to acceptable tolerances required for a stellarator magnet: (+0.25-mm from ideal on dimensions of coil positioning plates and up to +1-mm from ideal for position of wound coil). The plate materials is also included in this selection process from fabrication and 3D printing perspective and commensurate with eventual application to a fusion reactor. From the cost effectiveness point of view, the HTS coil and plate has the potential to cost less than that made in conventional methods with less waste (<75% waste) reducing time (<50%) and cost (<50%), especially as the AM field matures. The application of advanced manufacturing in the construction of the support plates will also lead to cost reduction, as the cables can be easily replaced, thereby making the assembly modular. With such high primary cost and time savings, high current densities and magnetic field, the funded R&D work has validated the designs, proven the feasibility, and characterized the performance of the HTS coil and plate assembly, paving the way for a relevant-size stellarator coil system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation

The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.

entropy closure↗
DOE OSTI2025

Development of an Advanced Multiphysics Simulation Capability for Radiant's Microreactor Design

Argonne National Laboratory and Idaho National Laboratory, through a Department of Energy Gateway for Accelerated Innovation in Nuclear Voucher, supported key analysis needs of Radiant related to (i) air jacket thermal fluid performance, (ii) evaluation of decay heat source terms defining air jacket technical requirements, and (iii) assessment of modeling methodologies employed for core analysis. All of these activities center on numerical simulation of various aspects of Kaleidos using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, the Cardinal multiphysics application, the OpenMC Monte Carlo code, and the Nek5000 computational fluid dynamics (CFD) code. This project builds upon an earlier Nuclear Energy Advanced Modeling and Simulation (NEAMS) Thermal-Hydraulic (T/H) Center of Excellence (CoE) project focused on initial demonstration of Cardinal multiphysics simulation of High Temperature Gas Reactors (HTGRs) and now focuses on Radiant’s Kaleidos concept.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Idealized simulations of wind farm interactions with intermittent turbulence in stable boundary layer conditions

Stable atmospheric boundary layer conditions typically correspond to weak turbulence levels, but intermittent periods of elevated turbulence can occur during otherwise quiescent conditions. The interaction between intermittent turbulence and wind turbines is not well understood because of sparse observations, as well as the difficulty in realistically resolving small-scale turbulence during strongly stable conditions with numerical simulations. In this study, an explicit filtering and reconstruction approach for large-eddy simulation (LES) is used to simulate weakly and strongly stable conditions, with surface cooling rates of −0.2 and −2.0 K h −1 , respectively. This approach can sustain resolved background turbulence at relatively coarse grid spacing and stronger stratification compared to conventional closures, permitting more realistic intermittent stable boundary layer (SBL) turbulence. The idealized LES capability of the Weather Research and Forecasting model is employed with turbine rotors parameterized using generalized actuator disks to examine (1) how the presence of turbine wakes affects SBL evolution and (2) the effect of intermittent turbulence on power production and wake recovery. Wakes increase mixing and deepen the SBL, with a stronger effect under strongly stable conditions, primarily because the SBL is shallower and closer to the top of the wind turbine rotor layer. Intermittent turbulence does not have a significant impact on mean power generation and wake recovery because the relevant intermittent turbulent structures in this study only affect the bottom half of the rotor disk. Power production is, however, more variable during periods of elevated turbulence, demonstrating the impact of SBL intermittency. This study uses an idealized configuration, focusing on LES model performance and physical understanding, with the goal of informing future simulations of the conditions observed during the American Wake Experiment.

Energy - Wind↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗
DOE OSTI2024

False vacuum decay and nucleation dynamics in neutral atom systems

Metastable states of quantum many-body systems with confinement offer a means to simulate false vacuum phenomenology, including nonequilibrium dynamical processes like decay by nucleation, in truncated limits. Recent work has examined the decay process in one-dimensional (1D) ferromagnetic Ising spins and superfluids. Here, in this paper, we study nucleation dynamics in 1D antiferromagnetic neutral atom chains with Rydberg interactions, using both numerical simulations and analytic modeling. We apply a staggered local detuning field to generate the metastable and ground states. Our efforts focus on two dynamical regimes: decay and annealing. In the first, we corroborate the phenomenological decay rate scaling and determine the associated parameter range for the decay process; in the second, we uncover and elucidate a procedure to anneal the metastable state from the initial to the final system, with intermediate nucleation events. We further propose experimental protocols to prepare the required states and perform quenches on near-term neutral atom quantum simulators, examining the experimental feasibility of our proposed setup and parameter regime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The influence of kinematics of deformation on polycrystalline halite dynamic recrystallization: Full-field simulation of simple shear versus pure shear

Rock salt, composed mainly of halite, flows viscoplastically over a wide range of geological conditions, strongly impacting the dynamic evolution of sedimentary basins and orogens. Understanding how dislocation creep, which involves dislocation glide, intracrystalline recovery and dynamic recrystallization, influences the microstructure and rheology of halite under various deformation kinematics and temperatures is crucial for enhancing knowledge of salt flow dynamics. Here, this study employs a full-field numerical simulation method to compare the viscoplastic deformation of polycrystalline halite under simple shear and pure shear conditions up to a natural strain of ε = 1.5 at temperatures ranging from 100 °C to 300 °C. The results are presented in terms of crystallographic preferred orientation (CPO), grain shape preferred orientation (SPO), subgrain boundary direction, grain size and strain rate distribution. The results indicate that the crystallographic anisotropy of individual halite crystals is transferred to the polycrystalline scale, resulting in strain localization, particularly in simple shear simulations. The kinematics of deformation affect the evolution and distribution of high strain-rate bands, determining the direction of intragranular substructures and the morphology of strain-induced grain boundaries, with minimal impact on grain size. The intensity of grain boundary migration increases with temperature, significantly influencing grain morphology and size, thereby obscuring strain localization, while having little effect on CPOs. At low strain (ε < 1.0), CPOs relative to both the maximum shortening direction and the grain SPO are similar regardless of the deformation kinematics. At high strain (ε > 1.0), simple shear CPOs exhibit three stronger {100} maxima with a monoclinic symmetry relative to the grain SPO compared to the six {100} maxima with an orthotropic symmetry relative to the grain SPO generated under pure shear. Therefore, microstructures and CPOs can serve as indicators of the strain path in polycrystalline halite under various conditions, aiding in determining the shear sense and elucidating the deformation kinematics of salt structures.

58 GEOSCIENCES↗
DOE OSTI2024

Two-dimensional coherent spectrum of high-spin models via a quantum computing approach

Here in this work we present and benchmark a quantum computing approach to calculate the two-dimensional coherent spectrum (2DCS) of high-spin models. Our approach is based on simulating their real-time dynamics in the presence of several magnetic field pulses, which are spaced in time. We utilize the adaptive variational quantum dynamics simulation algorithm for the study due to its compact circuits, which enables simulations over sufficiently long times to achieve the required resolution in frequency space. Specifically, we consider an antiferromagnetic quantum spin model that incorporates Dzyaloshinskii-Moriya interactions and single-ion anisotropy. The obtained 2DCS spectra exhibit distinct peaks at multiples of the magnon frequency, arising from transitions between different eigenstates of the unperturbed Hamiltonian. By comparing the one-dimensional coherent spectrum with 2DCS, we demonstrate that 2DCS provides a higher resolution of the energy spectrum. We further investigate how the quantum resources scale with the magnitude of the spin using two different binary encodings of the high-spin operators: the standard binary encoding and the Gray code. At low magnetic fields both encodings require comparable quantum resources, but at larger field strengths the Gray code is advantageous. Numerical simulations for spin models with increasing number of sites indicate a polynomial system-size scaling for quantum resources. Lastly, we compare the numerical 2DCS with experimental results on a rare-earth orthoferrite system. The observed strength of the magnonic high-harmonic generation signals in the 2DCS of the quantum high-spin model aligns well with the experimental data, showing significant improvement over the corresponding mean-field results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Efficient Simulation of Logical Magic State Preparation Protocols

Developing space- and time-efficient logical magic state preparation (MSP) protocols will likely be an essential step toward building a large-scale fault-tolerant quantum computer. Motivated by this need, we introduce a scalable method for simulating logical MSP protocols under the standard circuit-level noise model. When applied to protocols based on code-switching, magic state cultivation, and magic state distillation, our method yields a complexity polynomial in (i) the number of qubits and (ii) the nonstabilizerness, e.g., stabilizer rank or Pauli rank, of the target encoded magic state. The efficiency of our simulation method is rooted in a curious fact: every circuit-level Pauli error in these protocols propagates to a Clifford error at the end. This property is satisfied by a large family of protocols, including those that repeatedly measure a transversal Clifford that squares to a Pauli. We provide a proof-of-principle numerical simulation that prepares a magic state using such logical Clifford measurements. Our work enables practical simulation of logical MSP protocols without resorting to approximations or resource-intensive state-vector simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Hybrid quantum simulations with qubits and qumodes on trapped-ion platforms

We explore the feasibility of gate-based hybrid quantum computing using both discrete (qubit) and continuous (qumode) variables on trapped-ion platforms. Trapped-ion systems have demonstrated record one- and two-qubit gate fidelities and long qubit coherence times, while qumodes, which can be represented by the collective vibrational modes of the ion chain, have remained relatively unex- plored for their use in computing. Using numerical simulations, we show that high-fidelity hybrid gates and measurement operations can be achieved for existing trapped-ion quantum platforms. As an exemplary application, we consider quantum simulations of the Jaynes-Cummings-Hubbard model, which is given by a one-dimensional chain of interacting spin and boson degrees of free- dom. Using classical simulations, we study its real-time evolution and develop a suitable variational quantum algorithm for ground state preparation. Furthermore, our results motivate further studies of hybrid quantum computing in this context, which may lead to direct applications in condensed matter and fundamental particle and nuclear physics.

Lower-dimensional field theories↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗