Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed 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 595 records · Page 33

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗

Thermophysical Properties of Liquid Tritium: A Path Integral Monte Carlo Study

Here, we present worm-algorithm, path integral Monte Carlo simulations of bulk liquid tritium. The simulations are benchmarked against empirically known thermophysical properties of liquid deuterium and liquid tritium. Results for the pair correlation function, chemical potential, isothermal compressibility, isochoric heat capacity, and single-particle momentum distributions are reported. Given the benchmark comparisons, our predictions of liquid tritium properties are expected to be accurate to within a few percent. Our simulations unambiguously demonstrate the significance of nuclear quantum effects to the properties of liquid tritium. In particular, under saturated vapor pressure, the average molecular kinetic energy of the liquid is found to be more than 60% higher than the value expected from the classical equipartition theorem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Xyce™ Parallel Electronic Simulator Users' Guide (V.7.9)

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: • Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. • A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. • Device models that are specifically tailored to meet Sandia’s needs, including some radiation-aware devices (for Sandia users only). • Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

42 ENGINEERING↗

Xyce™ Parallel Electronic Simulator Users’ Guide, Version 7.10

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: • Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. • A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. • Device models that are specifically tailored to meet Sandia’s needs, including some radiation-aware devices (for Sandia users only). • Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

97 MATHEMATICS AND COMPUTING↗

A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation

Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Q-Cluster: Quantum Error Mitigation Through Noise-Aware Unsupervised Learning

Quantum error mitigation (QEM) is critical in reducing the impact of noise in the pre-fault-tolerant era, and is expected to complement error correction in fault-tolerant quantum computing (FTQC). In this work, we propose a novel QEM approach, Q-Cluster, that uses unsupervised learning (clustering) to reshape the measured bit-string distribution. Our approach starts with a simplified bit-flip noise model. It first performs clustering on noisy measurement results, i.e., bit-strings, based on the Hamming distance. The centroid of each cluster is calculated using a qubit-wise majority vote. Next, the noisy distribution is adjusted with the clustering outcomes and the bitflip error rates using Bayesian inference. Our simulation results show that Q-Cluster can mitigate high noise rates (up to 40% per qubit) with the simple bit-flip noise model. However, real quantum computers do not fit such a simple noise model. To address the problem, we (a) apply Pauli twirling to tailor the complex noise channels to Pauli errors, and (b) employ a machine learning model, ExtraTrees regressor, to estimate an effective bit-flip error rate using a feature vector consisting of machine calibration data (gate & measurement error rates), circuit features (number of qubits, numbers of different types of gates, etc.) and the shape of the noisy distribution (entropy). Our experimental results show that our proposed Q-Cluster scheme improves the fidelity by a factor of 1.46x, on average, compared to the unmitigated output distribution, for a set of low-entropy benchmarks on five different IBM quantum machines. Our approach outperforms the state-of-art QEM approaches RZNE [28], M3 [24], Hammer [35], and QBEEP [33] by 1.26x,1.29x,1.47x, and 2.65 x, respectively.

42 ENGINEERING↗

Impact of super-Gaussian electron distributions on plasma K-shell emission

Electron distributions in laser-produced plasmas will be driven toward a super-Gaussian distribution due to inverse bremsstrahlung absorption [Langdon, Phys. Rev. Lett. 44, 575 (1980)]. Both theoretical and experimental evidence suggest that fundamental plasma properties are altered by the super-Gaussian distribution. Here, this paper examines how the super-Gaussian distribution affects the ionization balance and K-shell emission of atomic plasmas, utilizing approximate formulas and detailed collisional-radiative simulations. While the impact on plasma ionization is small, K-shell spectra can be significantly modified. Based on these findings, we demonstrate that K-shell spectroscopy can be used to infer super-Gaussian or other similar nonequilibrium electron distributions.

Atomic spectra↗

Analysis of Screening Current Effects in a Hybrid Nb$_{3}$Sn/REBCO Superconducting Accelerator Magnet Using a T - A Formulation

To explore the feasibility of using high-temperature superconducting (HTS) REBCO coated conductors in future accelerator magnets, two REBCO flat racetrack coils were fabricated using 4-mm wide EuBCO tapes at the High Energy Accelerator Research Organization (KEK). These coils were tested as an insert inside a Nb$_{3}$Sn common-coil dipole magnet, which provides a background field of up to $\sim$ 9.5T, at the Brookhaven National Laboratory (BNL). REBCO tapes offer exceptionally high critical current density under strong magnetic fields; however, they also exhibit significant magnetization due to screening currents, leading to magnetic field errors. Here, this study presents a 2D finite element model of screening current-induced fields (SCIF) in REBCO coils using the T-A formulation, along with the results obtained. Simulations were then performed for two KEK test cases: one where the REBCO conductors were oriented with the HTS tapes parallel to the background field, and another where the tapes were perpendicular to it. Since screening currents also influence the stress distribution and increase the peak stress in the coils, the mechanical effects of these currents were analyzed. The implications of these simulation and test results for the design of Nb$_{3}$Sn/REBCO superconducting accelerator magnets are discussed.

Accelerator magnets↗

Modeling the Effect of Material Properties on Liquid-Alkaline Water Electrolysis

Liquid-alkaline water electrolyzers (LAWEs) use electricity to drive the conversion of water to H 2 and O 2 gas. These devices benefit from the use of low-cost nickel electrodes and metal-oxide separators, but suffer from lower current densities and higher cell voltages than proton-exchange-membrane water electrolyzers. Identifying the inefficiencies that result in this poor performance is key to mitigating losses and optimizing LAWEs. Here, we report an experimentally-validated 1-D continuum model of a LAWE that elucidates the gradients within the cell, simulates H 2 crossover, and projects the energy improvements made possible by modulating the properties of the electrodes and separator. The model captures the Nernstian polarization losses and the distribution of gas- and liquid-phases within the electrodes, enabling quantification of energy losses associated with kinetic, ohmic, and bubble-induced (mass-transport) resistances. Simulations demonstrate that LAWE can achieve energy intensities of 50 kWh kg −1 of H 2 at 1 A cm −2 using improved electrode and separator properties.

08 HYDROGEN↗

Heterogeneous Multilayer Nanopores via Chemically Tuned Dielectric Breakdown for Single‐Molecule Sensing

Solid-state nanopores are powerful platforms for single-molecule sensing, yet their performance is often constrained by fabrication complexity, noise, and limited control over surface properties. Here we report a direct method to fabricate heterogeneous multilayer nanopores using chemically tuned controlled dielectric breakdown (CT-CDB). We integrate hBN, MoS 2 , or graphene atop a silicon nitride membrane to form five distinct bilayer and tri-layer architectures, with bare SiN x nanopore as a control. CT-CDB achieves pore formation reproducibly through material-stacks with high efficiency, good pore size control, and strong yield, validated by various characterizations. Transferrin protein translocation experiments, supported by simulations, reveal that multilayer configurations modulate protein conformations, ionic current blockade and dwell time distributions, reflecting combined effects of membrane type, interfacial chemistry, and local electric field gradients. A supervised machine learning framework is implemented to assist identifying multilayer structure effects embedded in signal signatures, with over 96% accuracy. This work presents a modular and scalable framework for functional nanopore engineering with complex structural integration, thereby expanding the potential of 2D materials in single-molecule sensing applications.

2D materials↗

Investigation of fuel film formation and soot emissions in a GDI engine during cold-start with Split-injection strategies

This study investigates the impact of fuel-film formation on engine-out soot emissions in a gasoline direct injection (GDI) engine under cold-start conditions. Split-injection strategies were applied by varying the number of injections, injection duration, and total fuel quantity to affect wall wetting and control the average in-cylinder equivalence ratio. A combined experimental and numerical approach was employed to analyze fuel-film deposition, combustion efficiency, and engine-out soot and unburnt hydrocarbons (UHC) emissions. In particular, fuel film distribution estimated by means of non-reacting, 3-D, computational fluid dynamics (CFD) simulations, together with experimentally measured soot data, were used to investigate fuel film formation and its role in soot generation. In the experiments, a skip-fired engine control strategy was applied to mimic the transient nature of engine cold-start operation. The results indicate that, under the same total number of injections, increasing the average in-cylinder equivalence ratio through longer injection durations improves combustion stability, as indicated by the decrease of the coefficient of variation of IMEP n (Net Indicated Mean Effective Pressure) from 6.0 % to 0.6 %. However, this strategy leads to higher soot emissions, which increased by nearly an order of magnitude, primarily due to enhanced wall-film formation. In contrast, increasing the number of injections while maintaining a constant equivalence ratio significantly impacts fuel-film deposition and, consequently, soot emissions, with a fivefold reduction of the measured engine-out soot, decreasing from 3.5 mg to 0.7 mg. A soot-film correlation was developed and achieved a high coefficient of determination (r 2 = 0.95) and was further extended to account for spark timing effects. These findings confirm the effectiveness of split-injection for avoiding wall film formation and soot emissions, and the critical role of fuel film in soot generation, supporting the hypothesis that pool fires play a crucial role in contributing to soot formation under these cold-start conditions. In conclusion, the study also indicates the value of a predictive soot-film correlation for developing cold-start emission control strategies.

Engine Cold-start↗

Coherent and Dynamic Small Polaron Delocalization in CuFeO 2

Small polarons remain a bottleneck in realizing efficient transition metal oxide devices. Routes to engineer small polaron coupling to electronic states and lattice modes to control carrier localization remain unclear. Here, we measure small polaron formation in CuFeO 2 using transient extreme ultraviolet reflection spectroscopy and compare to theoretical predictions in realistically parametrized Holstein models, demonstrating that polaron localization depends on coupling to high-frequency versus low-frequency phonon bath components. We measure small polaron formation on a comparable ∼100 fs timescale to other Fe(III) compounds. Dynamic delocalization of the polaron follows formation through a coherent lattice expansion between Fe–O layers and charge-sharing with surrounding Fe(IV) states. Simulations reveal two major factors dictate polaron formation timescales: phonon density and reorganization energy distributions between acoustic and optical modes, matching experimental findings. Our work shows how electronic-structural coupling in a polaron-host material can be leveraged to suppress polaronic effects for various applications.

Hematite↗

Investigation of Isobaric Mixing as a Mechanism for Boundary‐Layer Cloud Formation

This study investigates the potential role of isobaric mixing in the formation of marine boundary layer clouds. Cloud formation theory emphasizes uplift and adiabatic cooling, but recent observations support the existence of small clouds forming at various altitudes, even below the lifting condensation level. Isobaric mixing of air with different thermodynamic properties can generate localized supersaturation. A Gaussian mixing model is employed to simulate this process, considering the correlation between temperature and water vapor. Cloud droplet size distributions from aircraft measurements show a persistent and prominent mode of small droplets at 9 m, and the size of this mode compares favorably with predictions from the model. The results suggests that isobaric mixing plausibly contributes to the formation of clouds, particularly those observed at multiple altitudes with narrow droplet size distributions. This finding highlights the importance of considering isobaric mixing processes in understanding and modeling cloud formation.

54 ENVIRONMENTAL SCIENCES↗

Chlorine gas and anion radical reactivity in molten salts and the link to chlorobasicity

Next-generation nuclear power plants may include exciting novel designs in which molten salts are the coolant or a combination of the coolant and fuel. Whereas it is straightforward to see why having a low volatility coolant can be advantageous for safety, much is not understood about the production of volatile halogen gases as a result of radiation and even less is known about the distribution of these species at and away from interfaces. Using first principles molecular dynamics simulations, we investigate the product of the disproportionation reaction between chlorine anion radicals (nominally Cl 2 ˙ − ) in the bulk and slab configurations. We find that the product depends on the chlorobasicity of the medium. For example, in ZnCl 2 , Cl 2 forms, but in a eutectic mixture of LiCl and KCl, Cl 3 − is formed as a product. We also find that Cl 3 − prefers to form at the vapor interface and this may have implications for corrosion and reactivity. Furthermore, the mechanisms of the mobility of Cl 2 and Cl 3 − are radically different, the first one being vehicular and the second Grotthus-like. Chlorobasicity is linked to the electronic structure of the host melt; ZnCl 2 forms extended networks along which metal ions and anionic counterions have significant electronic orbital overlap forming long, linear, molecular-like constructs; the opposite is true for the alkali metal eutectic salt.

Nguyen, Hung H. [Univ. of Iowa, Iowa City, IA (Uni↗

Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+

Chemistry↗

Accurate field-level weak lensing inference for precision cosmology

We present miko, a catalog-to-cosmology pipeline for general flat-sky field-level inference, which provides access to cosmological information beyond the two-point statistics. In the context of weak lensing, we identify several new field-level analysis systematics (such as aliasing, Fourier mode-coupling, and density-induced shape noise), quantify their impact on cosmological constraints, and correct the biases to a percent level. Next, we find that model misspecification can lead to both absolute bias and incorrect uncertainty quantification for the inferred cosmological parameters in realistic simulations. The Gaussian map prior infers unbiased cosmological parameters, regardless of the true data distribution, but it yields overconfident uncertainties. The log-normal map prior quantifies the uncertainties accurately, although it requires careful calibration of the shift parameters for unbiased cosmological parameters. Here, we demonstrate systematics control down to the 2% level for both models, making them suitable for ongoing weak lensing surveys.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗

A Survey of Open-Source Tools for Transmission and Distribution Systems Research

This work presents a review of open-source electric power transmission and distribution systems analysis tools suitable for use by industry professionals and academic researchers. Due to the high complexity of the electric grid, there exist numerous tools and extensive research pertaining to nearly every aspect of the design, operation, and control of transmission and distribution networks. In addition to the commercial tools, a wide range of free, open-source tools, models, and data usable by the scientific community for related research have been developed by different organizations, including both international and US universities and national laboratories. However, due to the absence of a catalog of available tools, models and data, researchers often lack a knowledge of existing capabilities and may develop duplicative software and tools. Increasing awareness of these available resources seeks to accelerate their broader use, leading to more efficient and standardized grid analysis. This review paper (which is part of a larger survey effort that studied over 400 tools in the transmission, distribution, buildings, and electric vehicles space) outlines selected open-source resources that have been developed in power transmission and distribution systems research. It is anticipated that this work can serve as a guide for industry and academic researchers alike, ensuring that research efforts are well-channeled.

24 POWER TRANSMISSION AND DISTRIBUTION↗