Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-idealities”

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

Passive Microwave Radiometry and Active Radar Sounding as Complementary Tools for Geophysical Investigations of Icy Ocean Worlds

Juno Microwave Radiometer (MWR) observations of Europa and Ganymede offer critical insights into the icy shells of these moons ahead of NASA's Europa Clipper and ESA's JUpiter ICy moons Explorer (JUICE) missions. Both missions are equipped with active radar sounders designed to address key unknowns such as ice shell thickness, thermal state, and composition. In this study, we explore how passive microwave radiometry and active radar sounding can constrain ice shell properties, focusing on Europa. Using modeled microwave brightness temperature observations at 0.6 and 1.2 GHz alongside simulated radar attenuation rate observations, we show that each instrument can independently produce robust ice shell thickness constraints under idealized conditions. We then relax these assumptions, quantifying how uncertainties from non-ideal properties—including convective layers, freezing-point depression, and chloride-doped ice—affect thickness estimates. Finally, we demonstrate how combining observations from these complementary techniques breaks degeneracies between ice shell properties, enabling more robust constraints than either method alone. This approach will maximize the science return of Europa Clipper and JUICE, advancing our understanding of the thermophysical structure and habitability of icy ocean worlds.

58 GEOSCIENCES↗

Leveraging prior mean models for faster Bayesian optimization of particle accelerators

Tuning particle accelerators is a challenging and time-consuming task that can be automated and carried out efficiently using suitable optimization algorithms, such as model-based Bayesian optimization techniques. One of the major advantages of Bayesian algorithms is the ability to incorporate prior information about beam physics and historical behavior into the model used to make control decisions. In this work, we examine incorporating prior accelerator physics information into Bayesian optimization algorithms by utilizing fast executing, neural network models trained on simulated or historical datasets as prior mean functions in Gaussian process models. We show that in ideal cases, this technique substantially increases convergence speed to optimal solutions in high-dimensional tuning parameter spaces. Additionally, we demonstrate that even in non-ideal cases, where prior models of beam dynamics do not exactly match experimental conditions, the use of this technique can still enhance convergence speed. Finally, we demonstrate how these methods can be used to improve optimization in practical applications, such as transferring information gained from beam dynamics simulations to online control of the LCLS injector, and transferring knowledge gained from experimental measurements across different operating modes, such as accelerating different ion species at the ATLAS heavy ion accelerator.

43 PARTICLE ACCELERATORS↗

Modeling ethanol/water adsorption in all-silica zeolites using the real adsorbed solution theory

A comprehensive set of single-component and binary isotherms were collected for ethanol/water adsorption into the siliceous forms of 185 known zeolites using grand-canonical Monte Carlo simulations. Using these data, a systematic analysis of ideal/real adsorbed-solution theory (IAST/RAST) was conducted and activity coefficients were derived for ethanol/water mixtures adsorbed in different zeolites based on RAST. It was found that activity coefficients of ethanol are close to unity while activity coefficients of water are larger in most zeolites, indicating a positive excess free energy of the mixture. This observation can be attributed to water/ethanol interactions being less favorable than water/water interactions in the single-component adsorption of water at comparable loadings. The deviation from ideal behavior can be highly structure-dependent but no clear correlation with pore diameters was identified. Furthermore, our analysis also demonstrates the following: (1) accurate unary isotherms in the low-loading regime are critical for obtaining physically sensible activity coefficients; (2) the global regression scheme to solve for activity model parameters performs better than fitting activity models to activity coefficients calculated locally at each binary state point; and (3) including the dependence on adsorption potential offers only a minor benefit for describing binary adsorption at the lowest fugacities. Finally, the Margules activity model was found incapable of capturing the non-ideal adsorption behavior over the entire range of fugacities and compositions in all zeolites, but for conditions typical of solution-phase adsorption, RAST predictions using zeolite-specific or even bulk Margules parameters provide an improved description compared to IAST.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sparsified time-dependent Fourier neural operators for fusion simulations

This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Testing strong-field QED with the avalanche precursor

A two-beam high-power laser facility is essential for the study of one of the most captivating phenomena predicted by strong-field quantum electrodynamics (QED) and yet unobserved experimentally: the avalanche-type cascade. In such a cascade, the energy of intense laser light can be efficiently transformed into high-energy radiation and electron–positron pairs. The future 50-petawatt-scale laser facility NSF OPAL will provide unique opportunities for studying such strong-field QED effects, as it is designed to deliver two ultra-intense, tightly focused laser pulses onto the interaction point. In this work, we investigate the potential of such a facility for studying elementary particle and plasma dynamics deeply in the quantum radiation-dominated (RD) regime, and the generation of QED avalanches. With three-dimensional particle-in-cell simulations, we demonstrate that QED avalanche precursors can be reliably triggered under realistic laser parameters and layout (namely, focusing f/2, tilted optical axes, and non-ideal co-pointing) with the anticipated capabilities of NSF OPAL. We demonstrate that seed electrons can be efficiently injected into the laser focus by using targets of three types: a gas of heavy atoms, an overcritical plasma, and a thin foil. A strong positron and high-energy photon signal is generated in all cases. The cascade properties can be identified from the final particle distributions, which have a clear directional pattern. At increasing laser field intensity, such distributions provide signatures of the transition, first, to the RD interaction regime, and then to a QED avalanche. Our findings can also be used for designing related future experiments.

Lasers↗

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative x-ray scattering of free molecules

Advances in x-ray free electron lasers have made ultrafast scattering a powerful method for investigating molecular reaction kinetics and dynamics. Accurate measurement of the ground-state, static scattering signals of the reacting molecules is pivotal for these pump-probe x-ray scattering experiments as they are the cornerstone for interpreting the observed structural dynamics. Here, this article presents a data calibration procedure, designed for gas-phase x-ray scattering experiments conducted at the Linac Coherent Light Source x-ray Free-Electron Laser at SLAC National Accelerator Laboratory, that makes it possible to derive a quantitative dependence of the scattering signal on the scattering vector. A self-calibration algorithm that optimizes the detector position without reference to a computed pattern is introduced. Angle-of-scattering corrections that account for several small experimental non-idealities are reported. Their implementation leads to near quantitative agreement with theoretical scattering patterns calculated with ab-initio methods as illustrated for two x-ray photon energies and several molecular test systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Resilience–runtime tradeoff relations for quantum algorithms

Abstract A leading approach to algorithm design aims to minimize the number of operations in an algorithm’s compilation. One intuitively expects that reducing the number of operations may decrease the chance of errors. This paradigm is particularly prevalent in quantum computing, where gates are hard to implement and noise rapidly decreases a quantum computer’s potential to outperform classical computers. Here, we find that minimizing the number of operations in a quantum algorithm can be counterproductive, leading to a noise sensitivity that induces errors when running the algorithm in non-ideal conditions. To show this, we develop a framework to characterize the resilience of an algorithm to perturbative noises (including coherent errors, dephasing, and depolarizing noise). Some compilations of an algorithm can be resilient against certain noise sources while being unstable against other noises. We condense these results into a tradeoff relation between an algorithm’s number of operations and its noise resilience. We also show how this framework can be leveraged to identify compilations of an algorithm that are better suited to withstand certain noises.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling optical systematics for the Taurus CMB experiment

We simulate a variety of optical systematics for Taurus, a balloon-borne cosmic microwave background (CMB) polarisation experiment, to assess their impact on large-scale E-mode polarisation measurements and constraints of the optical depth to reionisation τ. We model a one-month flight of Taurus from Wanaka, New Zealand aboard a super-pressure balloon (SPB). We simulate night-time scans of both the CMB and dust foregrounds in the 150 GHz band, one of Taurus's four observing bands. We consider a variety of possible systematics that may affect Taurus's observations, including non-gaussian beams, pointing reconstruction error, and half-wave plate (HWP) non-idealities. For each of these, we evaluate the residual power in the difference between maps simulated with and without the systematic, and compare this to the expected signal level corresponding to Taurus's science goals. Our results indicate that most of the HWP-related systematics can be mitigated to be smaller than sample variance by calibrating with Planck's TT spectrum and using an achromatic HWP model, with a preference for five layers of sapphire to ensure good systematic control. However, additional beam characterization will be required to mitigate far-sidelobe pickup from dust on larger scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spatial Mapping of Chiral-Induced Spin Selectivity in Chiral Perovskite via Spin-Schottky Junction

Chiral halide perovskite (c-HP) semiconductors exhibit on average a large chiral-induced spin selectivity (CISS) effect. Nevertheless, the microscopic details of CISS and its integration in opto-spintronic constructs remain nascent. Reliable reporting of CISS performance characteristics represents a significant challenge in providing the necessary design rules. We show a Kelvin probe force microscopy (KPFM) method that can quantitatively evaluate and spatially map the chirality-dependent surface contact potential difference resulting from the formation of a spin-Schottky junction. We revealed inhomogeneity in the CISS response, where low-CISS regions in the c-HP films reduce the overall macroscopic average, likely serving as a key factor in optimizing macroscopic performance. We also observed that although c-HP films made from higher precursor concentrations lead to thicker films and higher carrier concentrations with subsequent larger barrier heights in the Schottky junction, stronger spin relaxation due to non-ideal film quality reduces spin polarization.

14 SOLAR ENERGY↗

C60-Based Ionic Salt Electron Shuttle for High-Performance Inverted Perovskite Solar Modules

Although C60 is usually the electron transport layer (ETL) in inverted perovskite solar cells, its molecular nature of C60 leads to weak interfaces that lead to non-ideal interfacial electronic and mechanical degradation. Here, we synthesized an ionic salt from C60, 4-(1',5'-dihydro-1'-methyl-2'H-[5,6] fullereno-C60-Ih-[1,9-c]pyrrol-2'-yl) phenylmethanaminium chloride (CPMAC), and used it as the electron shuttle in inverted PSCs. The CH2-NH3+ head group in the CPMA cation improved the ETL interface and the ionic nature enhanced the packing, leading to ~3-fold increase in the interfacial toughness compared to C60. Using CPMAC, we obtained ~26% power conversion efficiencies (PCEs) with ~2% degradation after 2,100 hours of 1-sun operation at 65degrees C. For minimodules (four subcells, 6 centimeters square), we achieved the PCE of ~23% with <9% degradation after 2,200 hours of operation at 55degrees C.

14 SOLAR ENERGY↗

Quantitative Voltammetry Measurements of High-Concentration Actinides in Molten Chloride Fuel Salts

Molten chloride mixtures containing high concentrations of UCl 3 (i.e., >70 wt%) are candidates for next-generation molten salt reactor (MSR) fuel salts that enable new plant optimization and fuel cycle strategies. MSR licensing, however, is contingent on the development of effective material control and accounting (MC&A) practices for these salts. To close this gap, we have developed an electroanalytical approach that enables accurate in situ measurements of uranium concentrations in these highly loaded fuel salts. This methodology uses a multielectrode array voltammetry sensor combined with digital simulations that account for non-idealities such as uncompensated resistance effects. The theory-based numerical simulations of voltammetry responses yielded correction factors that are agnostic of electrochemical cell geometry and produced consistent voltammetry peak current results when correcting for varying amounts of ohmic resistance. In doing so, we have demonstrated accurate concentration measurements for UCl 3 in LiCl-KCl-UCl 3 across a range of mass fractions between 55 and 72.3 wt% at 550 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Impact of Catalyst Layer Properties and Transport Layer Interactions on Rutile Iridium Oxide Anode Durability in Proton Exchange Membrane Electrolysis

The performance and durability of proton exchange membrane electrolyzers are governed by the interplay between catalyst inks and integration, the resulting catalyst layer (CL) structure, and the interface between the CL and porous transport layer (PTL). In this study, we examine the impact of these factors on cell performance and degradation using rutile IrO2 as the anode catalyst. The two commercial rutile IrO2 catalysts differ from each other and from the amorphous benchmark in agglomerate size and ink stability. This has implications for CL morphology and electronic connectivity. In the context of long-term pressurized operation for up to 1000 h, the selection of the PTL has a marginal effect on the durability of CLs with good structural and compositional properties, whereas poorly structured CLs demonstrate evident degradation when used with certain PTLs. Poor CL properties amplify the negative impact of non-ideal CL/PTL interfaces, demonstrating the importance of CL/PTL interface engineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analog Systems for Edge Optimization

Over the past decade, analog computing has the subject of substantial research interest providing a path toward improved computational efficiency in the post-Dennard era. Analog matrix vector multiplication (MVM) accelerators provide a popular approach given the ubiquity of MVM operations in numerous applications. However, historically analog computing systems can struggle with applications requiring high precision due to the inherent susceptibility of these systems to analog non-idealities. Therefore, prior work on analog systems has focused either on applications known to be tolerant of limited precision (e.g., neural network inference), or using expensive techniques to emulate high-precision using many analog MVM operations. In this work, we propose an alternative approach. Motivated by recent advances in inexact nonlinear solvers and optimizers, we explore the potential of co-designing optimization algorithms which can take full advantage of the fundamentally inexact analog MVM operations. To enable these co-designed algorithms we also develop a general mathematical theory of the precision and energy efficiency of analog operations, and a new system architecture for tightly-coupled analog and digital computation. Finally, we examine the applicability of analog computing to a wider class of symmetric positive definite systems and find potential in using analog operations as a sparse approximate inverse preconditioner. With these core innovations, this project provides a path toward effectively implementing optimization algorithms on power-constrained autonomous and semi-autonomous systems.

97 MATHEMATICS AND COMPUTING↗

Improved Muon Energy Estimation Using a Detailed Model of Multiple Coulomb Scattering in the MicroBooNE LArTPC

We present an improved technique for estimating a muon's energy by measuring the deflections along its path inside the MicroBooNE detector from multiple Coulomb scattering (MCS). This approach implements several innovations that better capture detector non-idealizations compared to previous MCS-based muon energy estimators. As a result, it achieves improved resolution, reduced bias, and better data-model agreement. Using model simulation, for fully contained events the estimated bias is within 1\% and the estimated resolution narrows from 10\% to 4.3\% as muon energy increases from 0.1\,GeV to 2\,GeV. For events with particles exiting the detector volume, at least a meter of reconstructed muon track, and a muon energy below 2\,GeV, the estimated bias is less than 2\% and the estimated resolution varies from 7\% to 17\% over muon energy. These demonstrate significant improvements over the performance of previous work using an MCS-based energy estimator at MicroBooNE~\cite{mcs_2017}, which exhibited approximately twice worse resolution and a bias of 20\% over the same energy region. Data-model goodness-of-fit studies are used to validate the estimator's performance on data, showing good agreement within model uncertainties.

Cooper-Troendle, London [U. Pittsburgh (main); Fer↗