SEARCH · Search NASA

Results for “Problem solving”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

Summary Computational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

2,5‐Dimercapto‐1,3,4‐Thiadiazole (DMCT)‐Based Polymers for Rechargeable Metal–Sulfur Batteries

Organosulfur materials are a sustainable alternative to the present-day layered oxide cathodes in lithium-based batteries. One such organosulfur material that was intensely explored from the 1990s to early 2010s is 2,5-dimercapto-1,3,4-thiadiazole (DMCT). However, research interest declined as the electrode reactions with DMCT were assumed to be too sluggish to be practical. Armed with the advances in metal–sulfur batteries, we revisit DMCT-based materials in the form of poly[tetrathio-2,5-(1,3,4-thiadiazole)], referred to as pDMCT-S. With an appropriate choice of electrode design and electrolyte, pDMCT-S cathode paired with a Li-metal anode shows a capacity of 715 mA h g −1 and a Coulombic efficiency of 97.7% at a C/10 rate, thus quelling the concerns of sluggish reactions. Surprisingly, pDMCT-S shows significantly improved long-term cyclability compared to a sulfur cathode. Investigations into the origin of the stability reveals that the discharge product Li-DMCT in its mesomeric form can strongly bind to polysulfides, preventing their dissolution into the electrolyte and shuttling. This unique mechanism solves a critical problem faced by sulfur cathodes. Encouragingly, this mechanism results in a stable performance of pDMCT-S with Na-metal cells as well. In conclusion, this study opens the potential for exploring other organic materials that have inherent polysulfide sequestering capabilities, enabling long-life metal–sulfur batteries.

2,5-dimercapto-1,3,4-thiadiazole↗

DOE OSTI2026

Extremized nonlinear and linearized responses in soft metamaterials enabled by gradient-based design and grayscale digital light processing

In this study, we develop a gradient-based design approach that exploits grayscale digital light processing (DLP) 3D printing for extremizing the nonlinear and linearized response of soft metamaterials — materials that harness engineered geometric instabilities to undergo large and programmable changes in configuration. Grayscale DLP approaches modulate local mechanical properties at the pixel scale by tuning the light intensity within a single grayscale image, unlocking an exceptionally large design space. To effectively navigate this space, we develop smooth mappings between local light intensity values and global quantities of interest that characterize the behavior of soft metamaterials. Enabling these smooth mappings are robust and differentiable nonlinear finite element simulations powered by a trust region solver. A PDE-constrained optimization problem is then solved to invert these mappings and produce light intensity distributions that endow the printed part with varying stiffness and flexibility in distinctive regions. It is shown that optimizing the distribution of soft and stiff phases throughout a metamaterial structure results in markedly different buckling and self-contact configurations to drive extremized nonlinear compression and linearized vibration responses. Optimized light intensity distributions are translated to grayscale images and directly used to print soft metamaterial samples, showing remarkable agreement between the buckling and self-contact response in simulated and measured deformed configurations.

Additive manufacturing↗

DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Multidimensional modeling of fuel-cladding friction in an LWR fuel rod

Solving finite element problems with friction often increases the level of difficulty to obtain properly converged solutions. This type of challenge becomes more salient when advanced, possibly multiscale, material models are employed to capture the thermomechanical behavior of fuel and cladding materials. Here the present work details our recent developments in a nuclear fuel performance finite element code for the systematic consideration of friction in nuclear reactor finite element simulations. We show the application of friction and its effects on the mechanics of light-water reactor rods accounting for various fuel constitutive modeling techniques, model dimensionalities, pellet assumed geometries, and power conditions. In particular, we focus on the fuel rod mechanical behavior as it relates to fuel constitutive models, sensitivity to the coefficient of friction, pellet states of stress, and rod elongation. We discuss the trade-offs between the various multidimensional modeling options and highlight the relevance of frictional effects in the prediction of the fuel rod deformation and interfacial stresses. To relate our modeling results with actual reactor operation, simulations including frictional effects are compared with fuel rod elongation experimental data and some challenges for carrying out a full validation of the axial mechanics are discussed.

42 ENGINEERING↗

DOE OSTI2025

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗

DOE OSTI2024

Simulating nonlinear optical processes on a superconducting quantum device

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Characterization of admissible quasisymmetries

We solve ‘half’ the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy magnetohydrostatic force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined nonlinear partial differential equations involving second derivatives of the symmetry's infinitesimal generator.

Burby, J. W. (ORCID:0000000317727647)↗

DOE OSTI2024

Ising meson spectroscopy on a noisy digital quantum simulator

Abstract Quantum simulation has the potential to be an indispensable technique for the investigation of non-perturbative phenomena in strongly-interacting quantum field theories (QFTs). In the modern quantum era, with Noisy Intermediate Scale Quantum (NISQ) simulators widely available and larger-scale quantum machines on the horizon, it is natural to ask: what non-perturbative QFT problems can be solved with the existing quantum hardware? We show that existing noisy quantum machines can be used to analyze the energy spectrum of several strongly-interacting 1+1D QFTs, which exhibit non-perturbative effects like ‘quark confinement’ and ‘false vacuum decay’. We perform quench experiments on IBM’s quantum simulators to compute the energy spectrum of 1+1D quantum Ising model with a longitudinal field. Our results demonstrate that digital quantum simulation in the NISQ era has the potential to be a viable alternative to numerical techniques such as density matrix renormalization group or the truncated conformal space methods for analyzing QFTs.

Science & Technology - Other Topics↗

DOE OSTI2025

Measuring error rates of mid-circuit measurements

High-fidelity mid-circuit measurements, which read out the state of specific qubits in a multiqubit processor without destroying them or disrupting their neighbors, are a critical component for useful quantum computing. They enable fault-tolerant quantum error correction, dynamic circuits, and other paths to solving classically intractable problems. But there are few methods to assess their performance comprehensively. In this work, we address this gap by introducing the first randomized benchmarking protocol that measures the rate at which mid-circuit measurements induce errors in many-qubit circuits. Using this protocol, we detect and eliminate previously undetected measurement-induced crosstalk in a 20-qubit trapped-ion quantum computer. Then, we use the same protocol to measure the rate of measurement-induced crosstalk error on a 27-qubit IBM Q processor, and quantify how much of that error is eliminated by dynamical decoupling.

Hothem, Daniel [Sandia National Laboratories (SNL-↗

DOE OSTI2026

A mathematical framework for thermodynamic computing with applications to chemical reaction networks

The widespread adoption of energy-intensive computing applications has led to a growing need for energy-efficient computing approaches. Thermodynamic computing offers a promising approach for low-energy computation by leveraging the intrinsic computational capabilities of physical, chemical, or biological systems. However, the mathematical foundations of thermodynamic computing require further development to fully realize the potential energy efficiencies, as well as to assess factors like noise and operational speed. In this paper, we establish a mathematical framework for utilizing thermodynamic processes to perform fundamental operations, including addition, subtraction, multiplication, and division. We highlight the use of chemical reactions as potential computational units and explore synthetic chemical and biochemical systems as practical implementations. Additionally, we demonstrate how these principles can be applied to solving complex mathematical problems, such as ordinary differential equations (ODEs) and suggest the necessary components to implement the thermodynamic computing framework using chemical reactions based in a microfluidic device. This work enhances our understanding of thermodynamic processes for natural computing as a basis for scalable, energy-efficient computation in paradigm disruptive next-generation systems.

Cannon, William R. [Pacific Northwest National Lab↗

DOE OSTI2025

Experimental progress and future plans on spherical tokamak, QUEST

QUEST (Q-shu university experiment with steady state spherical tokamak) aims at effective plasma current start-up and stable maintenance of plasma discharge. To solve the inherent problems in a spherical tokamak (ST) arising from insufficient space for placing the inductive center solenoid, electron cyclotron current drive (ECCD) and transient coaxial helicity injection (T-CHI) are implemented as a non-inductive plasma start-up method in QUEST. Efficient ECCD assisted by energetic electrons could be achieved. By combining control of the wave injection angle and application of a negative toroidal electric field, the bulk electron temperature could be raised up to 1 keV due to selective wave power absorption in the bulk electrons. The plasma current of over 50 kA contained within the closed flux surface could be obtained with a floating single biased electrode placed on lower divertor plates for T-CHI. Long-pulse operations on QUEST are impeded frequently due to wall saturation and subsequent density runaway caused by fuel particle imbalance. Since 2014, a unique tool called the ‘hot wall’ has been implemented to overcome the imbalance. The hot wall has a capability to regulate its surface temperature using a heater and two water cooling channels. With the help of the hot wall, 6 h discharges were obtained in 2020. Cooling down of the surface of the hot wall was significantly effective in recovering the wall pumping and was useful to extend the pulse duration. Augmentation of the toroidal magnetic field, B T up to 0.5 T from 0.25 T and a continuous wave (CW) gyrotron of 28 GHz are planned for QUEST in the near future. As raising B T provides a fundamental resonance of electron cyclotron waves (ECWs) with 28 GHz, more effective plasma current start-up and heating will be performed. Long-pulse operations with higher plasma parameters are expected.

QUEST↗

DOE OSTI2024

Control Co-Design Studies for a 22 MW Semisubmersible Floating Wind Turbine Platform

We present a control co-design software framework that can be used to optimize floating wind turbines and their controllers. Because this framework has many options for design variables, constraints, and merit figures, along with modeling fidelity levels, we seek to demonstrate best practices for using the tool while designing a floating platform for the new 22 MW offshore reference wind turbine developed within the International Energy Agency Wind Technology Commercialization Programme 55 on Reference Wind Turbines and Farms. During these studies, we evaluate the use of different simulation fidelity levels, the effect of using different load cases for controller tuning, and the difference between sequential and simultaneous control co-design solutions. Based on these efforts, we suggest using an algorithm that performs an initial search of the design space before optimization. We find that solving smaller optimization problems, in a sequential manner, leads to more reliable outcomes in fewer iterations than larger, simultaneous control co-design solutions. However a simultaneous CCD solution produces a platform with a 2% lower mass than the sequential CCD outcome.

17 WIND ENERGY↗

DOE OSTI2024

Uncovering the truth about M101, NGC 3938, and their significant others through radiative transfer

ABSTRACT Solving the inverse problem in spiral galaxies, that allows the derivation of the spatial distribution of dust, gas, and stars, together with their associated physical properties, directly from panchromatic imaging observations, is one of the main goals of this work. To this end, we used radiative transfer models to decode the spatial and spectral distributions of the nearby face-on galaxies M101 and NGC 3938. In both cases, we provide excellent fits to the surface-brightness distributions derived from GALEX, SDSS, 2MASS, Spitzer, and Herschel imaging observations. Together with previous results from M33, NGC 628, M51, and the Milky Way, we obtain a small statistical sample of modelled nearby galaxies that we analyse in this work. We find that in all cases Milky Way-type dust with Draine-like optical properties provide consistent and successful solutions. We do not find any ‘submm excess’, and no need for modified dust-grain properties. Intrinsic fundamental quantities like star-formation rates (SFR), specific SFR (sSFR), dust opacities, and attenuations are derived as a function of position in the galaxy and overall trends are discussed. In the SFR surface density versus stellar mass surface density space, we find a structurally resolved relation (SRR) for the morphological components of our galaxies, that is steeper than the main sequence (MS). Exception to this is for NGC 628, where the SRR is parallel to the MS.

Pricopi, D.↗

DOE OSTI2026

High-quality axion from exact SUSY chiral dynamics

We use supersymmetric chiral dynamics perturbed by anomaly-mediated supersymmetry breaking to obtain a high-quality, composite axion that solves the strong problem. The strong dynamics arises from a supersymmetric SU(10) chiral gauge theory with massless matter chiral superfields. This leads to a stable, nonsupersymmetric vacuum, calculated exactly, where a spontaneously broken global symmetry, identified with the Peccei-Quinn (PQ) symmetry, gives rise to a composite QCD axion. The chiral gauge theory also admits a discrete (or ) gauge symmetry that forbids PQ-violating operators up to dimension eight. An extension to an chiral gauge theory incorporates grand unification where the unification scale is identified with the PQ-breaking scale and PQ-violating operators are forbidden up to dimension 20. The supersymmetry breaking scale, near 100 TeV, ameliorates the Higgs hierarchy problem, while colored Nambu–Goldstone bosons may be detected at future colliders via decays to gluons or form heavy isotopes.

Gherghetta, Tony↗

DOE OSTI2024

Topological Kondo effect with spinful Majorana fermions

Motivated by the importance of studying topological superconductors beyond the mean-field approximation, we here investigate mesoscopic islands of time-reversal-invariant topological superconductors. We characterize the spectrum in the presence of strong order-parameter fluctuations in the presence of an arbitrary number of Kramers pairs of Majorana edge states and study the effect of coupling the Coulomb blockaded island to external leads. In the case of an odd fermionic parity on the island, we derive an unconventional Kondo Hamiltonian in which metallic leads couple to both topological Majorana degrees of freedom (which keep track of the parity in different leads) and the overall spin 1 2 in the island. For the simplest case of a single wire (two pairs of Majorana edge states), we demonstrate that anisotropies are irrelevant in the weak coupling renormalization group flow. This permits us to solve the Kondo problem in the vicinity of a Toulouse-type point using Abelian bosonization. We demonstrate a residual ground-state entropy of ln ( 2 ) , which is protected by spin-rotation symmetry, but reduced to ln ( 2 ) (as in the spinless topological Kondo effect) by symmetry-breaking perturbations. In the symmetric case, we further demonstrate the simultaneous presence of both Fermi-liquid and non-Fermi-liquid-like thermodynamics (depending on the observable) and derive charge and spin transport signatures of the Coulomb blockaded island. Published by the American Physical Society 2024

Materials Science↗

DOE OSTI2024

Spontaneous symmetry breaking, gauge hierarchy, and electroweak vacuum metastability

The so-called metastability bound asserts that an unnaturally small Higgs mass is a necessary condition for electroweak vacuum metastability, offering a new approach toward solving the hierarchy problem. So far, this result relies on the assumption of a negative Higgs mass parameter, or equivalently, on electroweak spontaneous symmetry breaking. We derive a new, corresponding bound for the case of a positive mass parameter. When the new bound is significantly more restrictive than or comparable to its established counterpart, it may offer an explanation for the sign of the Higgs mass parameter, and thus, spontaneous symmetry breaking itself. New physics at scales O ( 1 – 10 ) TeV can lower these bounds as far as the TeV scale. As an illustration, we consider vacuum stability in the presence of additional TeV-scale fermions with Yukawa couplings to the Higgs, as well as a dimension-six term parametrizing new physics in the UV. This scenario requires new physics that couples strongly to the Higgs, and can potentially be probed at future colliders. Finally, to allow for comparison with concrete mechanisms predicting metastability, we provide the mass-dependent lifetime of the electroweak vacuum for this model. Published by the American Physical Society 2024

Benevedes, Sean (ORCID:000000016390198X)↗

DOE OSTI2025

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗