SEARCH · Engineering Papers
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.
Here, the Chebyshev Rational Approximation Method (CRAM) algorithm has become the method of choice for modern nuclear depletion calculations. A simple modification of the CRAM yields not only depleted nuclide amounts but also time-averaged nuclide amounts, with little additional calculation cost.
Guiding-centre (GC) and full-orbit (FO) simulations of the beam-injected fast ion distribution and the corresponding neutron emissivity have been carried out for magnetohydrodynamics-quiescent National Spherical Torus eXperiment (NSTX) plasmas, using a combination of ASCOT5 and DRESS, to assess the suitability of the GC approximation for fast ions in NSTX. It was found that GC and FO simulations predicted substantially different steady-state distributions in both position and velocity space and different neutron emissivity profiles, leading to a 15% reduction in the predicted global neutron rate for FO relative to GC. These changes accompany a higher magnetic moment in FO, and correspond to a change in particle orbits from co-passing to trapped and stagnation orbits. ASCOT5 was also benchmarked against TRANSP/NUBEAM with input loaded entirely from TRANSP/NUBEAM output files, with agreement found between the GC simulations when finite Larmor radius (FLR) corrections were omitted. ASCOT5 FO and TRANSP/NUBEAM with FLR produced fast ion distributions which differed in localised regions, but predicted global neutron rates which agree within 3%.
Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.
ABSTRACT Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade’s largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon, and K-mouflage. In particular, we compare the predictions from N-body simulations solving the full scalar field equation, two N-body codes with approximate time integration schemes, a parametrized modified N-body implementation, and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in 2 per cent agreement at $z\le 1$ and $k\le 1~h\,{\rm Mpc}^{-1}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.
Here, we present an approach for including relativistic corrections in lepton-nucleus scattering calculations within the short-time approximation (STA). Previous ab initio studies employed electromagnetic currents expanded in powers of 𝑞/𝑚, where 𝑞 is the momentum transfer and 𝑚 is the nucleon mass, restricting their validity to low-𝑞 kinematics. We adopt an expansion scheme that treats the initial nucleon momentum perturbatively while allowing for arbitrary momentum transfer, thereby extending the applicability of the STA to high-𝑞 regimes. Additionally, we incorporate a relativistic treatment of the two-nucleon final-state energies. Calculations for 3 He and 4 He inclusive electron-scattering cross sections show a substantial improvement over previous results, achieving good agreement with experimental data in the quasi-elastic region for both low- and high-momentum transfer.
The quantum approximate optimization algorithm (QAOA) is a quantum heuristic for combinatorial optimization that has been demonstrated to scale better than state-of-the-art classical solvers for some problems. For a given problem instance, QAOA performance depends crucially on the choice of the parameters. While average-case optimal parameters are available in many cases, meaningful performance gains can be obtained by fine-tuning these parameters for a given instance. This task is especially challenging, however, when the number of circuit executions (shots) is limited. In this work, we develop an end-to-end protocol that combines multiple parameter settings and fine-tuning techniques. We use large-scale numerical experiments to optimize the protocol for the shot-limited setting and observe that optimizers with the simplest internal model (linear) perform best. We implement the optimized pipeline on a trapped-ion processor using up to 32 qubits and 5 QAOA layers, and we demonstrate that the pipeline is robust to small amounts of hardware noise. To the best of our knowledge, these are the largest demonstrations of QAOA parameter fine-tuning on a trapped-ion processor in terms of two-qubit gate count.
We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.
We generalize a recent prescription for the relaxation time approximation for the relativistic Boltzmann equation for systems with multiple particle species at finite temperature. This is performed by adding counter-terms to the traditional Anderson-Witting ansatz for each particle species. Our approach allows for the use of momentum-dependent relaxation times and the obedience of local conservation laws regardless of the definition of the local equilibrium state. As an application, we derive the first order Chapman-Enskog corrections to the equilibrium distribution and display results for the hadron-resonance gas. We also demonstrate that our collision term ansatz obeys the second law of thermodynamics.
The approximate SACOT-𝜒 scheme for heavy quark production in deep-inelastic scattering was initially formulated for the neutral current structure functions 𝐹 2 and 𝐹 𝐿 . We extend this approach to the charged current case (also including 𝐹 3 ), and thereby complete the definitions for the most relevant inclusive structure functions. Furthermore, we implement these structure functions in the open-source code APFEL++ which provides fast numerical evaluations over a wide kinematic range; this addition to the APFEL++ code is publicly available, with details provided in the Appendix. This SACOT-𝜒 implementation enables detailed numerical insights on the mass dependence of the structure functions and cross sections in the (𝑥,𝑄 2 )-plane for both neutral and charged current processes. We consider kinematic regions relevant for the experimental measurements from fixed-target 𝜈 DIS experiments (NuTeV, CCFR, and Chorus) and HERA, and also projections for the upcoming EIC. In particular, the 𝜈 DIS experiments reveal a surprisingly strong dependence on the mass effects, offering valuable insights that may help resolve long-standing challenges in accurately describing these datasets.
We present a method to construct high-order polynomial approximate invariants (AI) for nonintegrable Hamiltonian dynamical systems and apply it to a modern ring-based particle accelerator. Taking advantage of a special property of one-turn transformation maps expressed as square matrices, AIs can be constructed order by order iteratively. Evaluating AI with simulation data, we observe that AI’s fluctuation is actually a measure of chaos. Through minimizing the fluctuations, the stable region of long-term motions, i.e., the dynamic aperture of the accelerator, could be enlarged.
The short-time approximation is a method introduced to evaluate electroweak nuclear response for systems with A ≥ 12, extending the reach of first-principle many-body quantum Monte Carlo calculations. Using realistic two- and three-body nuclear interactions and consistent one- and two-body electromagnetic currents, we calculate longitudinal and transverse response densities and response functions of 12 C. Here, we compare the resulting cross sections with experimental data for electron-nucleus scattering, finding good agreement.
The Quantum Approximate Optimization Algorithm (QAOA) has emerged as a promising tool for accelerating optimization processes in the Noisy Intermediate-Scale Quantum (NISQ) era. Compared to classical methods, QAOA efficiently solves optimization problems, often formulated as Quadratic Unconstrained Binary Optimization (QUBO) problems. Classical quantum simulators are crucial for evaluating quantum algorithms due to limited quantum resources. However, QAOA's performance can vary with different simulation methods. This study analyzes QAOA's performance using various quantum simulators (e.g., density _matrix, statevector, and matrix_product_state) and demonstrates the benefits of HPC-QC integrated systems in solving QUBO problems on an active learning workflow. By simulating QAOA on dense, large-matrix QUBO problems, we evaluate accuracy and problem-solving time. We also assess QAOA's performance on local computers and HPC-QC inte-grated systems, using Oak Ridge Leadership Computing Facility (OLCF)'s Frontier supercomputer with local Qiskit Aer and remote IBM Quantum simulators.
Solving optimization problems is critical across many research domains, but the high dimensionality of parameter spaces often poses significant challenges. The Quantum Approximate Optimization Algorithm (QAOA) has emerged as a promising approach for accelerating optimization in the Noisy Intermediate-Scale Quantum (NISQ) era by leveraging both classical and quantum computational resources. However, its performance can vary depending on the underlying quantum hardware architecture. In this work, we evaluate the performance of QAOA on different quantum hardware platforms, specifically, superconducting transmon qubits and trapped-ion qubits, targetting real-world optimization problems formulated as fully connected Quadratic Unconstrained Binary Optimization (QUBO) instances. We evaluate both the solution quality and time-to-solution using dense QUBO matrices. Furthermore, we show that large-scale problems, such as a 100-bit QUBO instance, can be effectively tackled by integrating quantum computing with high-performance computing (HPC) resources. This study provides practical insights into the strengths and limitations of different qubit technologies and advances the application of quantum computing in solving real-world optimization problems.
The quantum approximate optimization algorithm (QAOA) is a leading candidate algorithm for solving optimization problems on quantum computers. However, the potential of QAOA to tackle classically intractable problems remains unclear. Here, we perform an extensive numerical investigation of QAOA on the low autocorrelation binary sequences (LABS) problem, which is classically intractable even for moderately sized instances. We perform noiseless simulations with up to 40 qubits and observe that the runtime of QAOA with fixed parameters scales better than branch-and-bound solvers, which are the state-of-the-art exact solvers for LABS. The combination of QAOA with quantum minimum finding gives the best empirical scaling of any algorithm for the LABS problem. We demonstrate experimental progress in executing QAOA for the LABS problem using an algorithm-specific error detection scheme on Quantinuum trapped-ion processors. Our results provide evidence for the utility of QAOA as an algorithmic component that enables quantum speedups.
SF-25-115 COBRA (COmputed-tomography Based Random-field Approximation) is a Python application for generating statistically equivalent random fields from CT-scan imagery. It leverages Karhunen–Loève expansions to model microstructural variability, enabling users to: Preprocess CT scans (filtering and Gaussian transformation); Fit covariance kernels fromempirical data; Solve eigenproblems to obtain KL modes; Sample random fields onsistent with fitted statistics; Postprocess samples back into the physical domain.
The linear coefficient of thermal expansion has been measured for 11 polymers from approximately -100 to +100°C. The polymers tested were poly(tetrafluoroethylene) (PTFE), poly(chlorotrifluoroethylene) (PCTFE), three polyethylenes; high density (HDPE), ultra-high molecular weight (UHMWPE) and cross-linked (XPE). The commercial fluoropolymers Kel-F 800 (also known as FK-800®) and THV 500®were tested together with samples of poly(ether ether ketone) (PEEK), poly(methyl methyl methacrylate) (PMMA), polycarbonate (PC) and poly(vinylidene fluoride) (PVDF). Polynomial fits were made to the data over appropriate ranges allowing the expansivity (α) to be estimated.
We present a method to construct high-order polynomial approximate invariants (AI) for non integrable Hamiltonian dynamical systems, and apply it to a modern ring-based particle accelerator. Taking advantage of a special property of one-turn transformation maps in the form of a square matrix, AIs can be constructed order-by-order iteratively. Evaluating AI with simulation data, we observe that AI’s fluctuation is actually a measure of chaos. Through minimizing the fluctuations, the stable region of long-term motions, i.e., the dynamic aperture of the accelerator, could be enlarged.
Molecular dynamics simulations are performed to study the stability of an fcc metallic interlayer phase between an iron aluminide quasicrystal approximant (QCA) phase (84.2 Al, 5.7 Fe, 5.4 Ni, and 4.7 at.% Cr) and aluminum oxides (α-Al 2 O 3 and γ-Al 2 O 3 ). The presence of the interlayer phase was experimentally observed in some regions of the inner diameter surface of a TPBAR cladding tube. The simulations employed existing many-body (MEAM and EAM) interatomic potentials to describe atomic interactions in the interlayer and QCA phases, and Buckingham+Coulomb pair interactions for the alumina and between the alumina and the QCA (or fcc). Interface systems of QCA/alumina, fec/alumina, and QCA/fcc are modeled and the fracture energy of the interfaces is calculated. Based on the fracture energies, for the γ-Al 2 O 3 , both MEAM and EAM suggest that the fcc interlayer phase would not form. For the α-Al 2 O 3 , the MEAM predicts that the fcc phase would form, however, the EAM predicts that it forms only if the α-Al 2 O 3 is Al-terminated at the interface. For the O-terminated α-Al 2 O 3 , the EAM predicts the fcc phase would not form, consistent with the results for γ-Al 2 O 3 .