Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical performance”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗
DOE OSTI2023

Compressible Rayleigh–Taylor instability subject to isochoric initial background state

The effect of compressibility on the single-mode Rayleigh–Taylor instability is examined using two (2D) and three-dimensional (3D) direct numerical simulations. To isolate compressibility from background stratification effects, this work employs a constant density profile on each side of the interface. The numerical simulations are performed at various Reynolds numbers using the gas kinetic method for static Mach numbers up to M = 0.4. The most important finding is that compressibility acting in isolation enhances the instability and perturbations grows faster with increasing Mach number, unlike previous results with background isothermal state, which show suppression of the instability at higher static Mach numbers. In addition, compressibility is also shown to increase the bubble-spike asymmetry. While the instability grows faster for the 3D case, the findings are qualitatively similar in 2D and 3D. Furthermore, the dynamical reasons underlying the effect of compressibility are elucidated by examining the evolution of vorticity and turbulent kinetic energy transport equations.

42 ENGINEERING↗
DOE OSTI2024

Nonplanar effects in simulations of laser-driven ejecta microjet experiments

Recent experiments of laser-driven ejecta microjets performed at OMEGA 60 reveal tortuous jets whereby the jets appear to deviate from their initial trajectory as they travel in vacuum. To understand these data, we perform two dimensional numerical simulations, considering different target thicknesses, pressures, and models of the drive conditions. In particular, modeling the finite laser spot size appears essential in reproducing qualitatively the non-planar shock observed in the experiment. Simulations capture jet deflection by accounting for a slight misalignment of the laser pointing with respect to the groove axis along with spatial variation of the laser pulse intensity. The principal physical mechanism appears to be that lateral momentum is imparted by release waves arising from the non-planar drive. The induced off-axis velocity is small in comparison to the jet axial velocity but integrates into a pronounced deflection over the course of the experiment. The analysis of jet axial and lateral mass distributions is found to be reproduced reasonably by the simulations. Simulated radiographs are in qualitative agreement with the experiments, though their differences point to potential shortcomings in modeling strictly three-dimensional experiments using two-dimensional hydrodynamic simulations. The simple analysis is able to explain part of the observed discrepancy in simulated vs experimental jet masses.

42 ENGINEERING↗
DOE OSTI2026

Locally purified maximally mixed states at scale: Entanglement pruning and symmetries

Locally Purified Density Operators (LPDOs) are state-of-the-art tensor network ansatze candidates that efficiently represent mixed quantum states at scale. However, given their non-uniqueness, their representational complexity is generally sub-optimal in practical computations. Here, in this work we perform a comprehensive numerical and analytical analysis and resolve this issue in the experimentally relevant limit where noise depolarizes the density operator into a maximally mixed state. To resolve the sub-optimality issue, we analyze two numerical tools, one analytic method, and detail the relations between them. The numerical tools used are fidelity-preserving truncations and isometric gauge transformations leveraging Riemannian optimizations over entropic objective functions. In addition, by invoking the injectivity and symmetry constraints of the maximally mixed LPDO, we also present analytical closed-form expressions for the disentangler and discuss their relation to numerical optimizers. Further, away from the maximally mixed state, our simulations highlight how the truncation threshold smoothly interpolate, as a function of depolarization, between established matrix product results and our new results. Our work shows how, by minimizing the resources required to represent key states of practical interest in experiment, the efficiency of tensor network algorithms can be substantially increased. This paves the path for uncovering tensor network’s fundamental scalability limits and latent potential in representing the wide locus of mixed quantum states that are accessible on near-term quantum devices.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗
DOE OSTI2023

Numerical integration of stochastic contact Hamiltonian systems via stochastic Herglotz variational principle

Within this work, we establish a stochastic contact variational integrator and its discrete version via stochastic Herglotz variational principle for stochastic contact Hamiltonian systems. A general structure-preserving stochastic contact method is provided to seek the stochastic contact variational integrators. Numerical experiments are performed to verify the validity of this approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Energy and particle balance during plasma detachment in a long-leg divertor configuration

Abstract Comprehensive studies of energy and particle balances in the transition to plasma detachment in an alternative divertor configuration with long outer legs are shown. Numerical simulations are performed with the 2D code suite SOLPS 4.3, using a disconnected double null grid with narrow, tightly baffled long poloidal leg divertors at the outer lower target and outer upper target. A particle count scan is performed using the ‘closed gas box’ model, where the tunable parameter in the simulations is the total number of deuterium particles in the simulation space and all other parameters are held fixed, including a constant input power and trace neon impurity radiation, to assess the physics of the transition to detachment in the system as the particle count increases. Three main aspects of the physics of divertor detachment are addressed: the criteria for the local onset of divertor detachment in each of the divertors, the distribution of heat flux and other plasma parameters between the four divertors as each divertor transitions to detachment, and the role of perpendicular transport in the transition to the detached regime. A synergistic mechanism by which the cross-field transport is reduced by factors associated with the onset of plasma recombination effects is identified. These results are compared to the existing understanding of the physics of the transition to plasma detachment in standard divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Characterization of W production during ICRF operations: experiments and modeling

For successfully heating plasma with waves in the ion cyclotron range of frequencies (ICRFs), mitigating impurity production is just as crucial as maximizing power coupling, especially in high-Z environments. ICRF can effectively deposit energy on ions, modify turbulence-driven transport, and enhance fusion reaction efficiency, but only when its power coupling has minimal impact on impurity production. To do so, one must rely on a toroidal array of at least three active elements excited with appropriate phasing and power ratio to reduce the currents induced on the antenna frame below levels critical for physical sputtering. In contrast to classic two-strap antennas, which are optimized for dipole phasing with equal power on both straps, three-strap antennas in ASDEX Upgrade (AUG)—but also four-strap antennas in JET, Alcator C-Mod, SPARC and ITER—offer the possibility to act also on the power ratio between the central and outer straps. With optimal settings, impurity production can be reduced substantially, making the ICRF compatible with the high-Z wall. This paper explores the characteristics of the AUG three-strap antennas in terms of impurity production, as well as the key role of plasma composition in this process. Numerical simulations were performed using SSWICH and Petra-M (finite element codes) to quantify impurity production and compare with experimental results. Energies of ions falling on antenna limiters (measured with probes) are well predicted by both codes. These tools are then used to further describe the source of the impurity, namely the gross erosion of tungsten from an ICRF antenna, for different plasma mixtures. Results are also compared to spectroscopy data. Ultimately, we show that deleterious effects of the ICRF on plasma surface interactions will be weaker in plasmas containing larger fractions of highly ionized heavier low-Z impurity, which is typically relevant for experiments relying on impurity seeding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The hybrid topological longitudinal transmon qubit

We introduce a new hybrid qubit consisting of a Majorana qubit interacting with a transmon longitudinally coupled to a resonator. To do so, we equip the longitudinal transmon qubit with topological quasiparticles, supported by an array of heterostructure nanowires, and derive charge- and phase-based interactions between the Majorana qubit and the resonator and transmon degrees of freedom. Inspecting the charge coupling, we demonstrate that the Majorana self-charging can be eliminated by a judicious choice of charge offset, thereby maintaining the Majorana degeneracy regardless of the quasiparticles spatial arrangement and parity configuration. We perform analytic and numerical calculations to derive the effective qubit-qubit interaction elements and discuss their potential utility for state readout and quantum error correction. Further, we find that select interactions depend strongly on the overall superconducting parity, which may provide a direct mechanism to characterize deleterious quasiparticle poisoning processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Charge-density-wave melting in the one-dimensional Holstein model

We study the Holstein model of spinless fermions, which at half filling exhibits a quantum phase transition from a metallic Tomonaga-Luttinger liquid phase to an insulating charge-density-wave (CDW) phase at a critical electron-phonon coupling strength. In our work, we focus on the real-time evolution starting from two different types of initial states that are CDW ordered: (i) ideal CDW states with and without additional phonons in the system and (ii) correlated ground states in the CDW phase. We identify the mechanism for CDW melting in the ensuing real-time dynamics and show that it strongly depends on the type of initial state. We focus on the far-from-equilibrium regime and emphasize the role of electron-phonon coupling rather than dominant electronic correlations, thus complementing a previous study of photoinduced CDW melting [H. Hashimoto and S. Ishihara, Phys. Rev. B 96, 035154 (2017)]. The numerical simulations are performed by means of matrix-product-state based methods with a local basis optimization (LBO). Within these techniques, one rotates the local (bosonic) Hilbert spaces adaptively into an optimized basis that can then be truncated while still maintaining a high precision. In this work, we extend the time-evolving block decimation (TEBD) algorithm with LBO, previously applied to single-polaron dynamics, to a half-filled system. We demonstrate that in some parameter regimes, a conventional TEBD method without LBO would fail. Furthermore, we introduce and use a ground-state density-matrix renormalization group method for electron-phonon systems using local basis optimization. In our examples, we account for up to M ph =40 bare phonons per site by working with O(10) optimal phonon modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Higher-order Van Hove singularities in kagome topological bands

Motivated by the growing interest in band structures featuring higher-order Van Hove singularities (HOVHS), we investigate a spinless fermion kagome system characterized by nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping amplitudes. While NN hopping preserves time-reversal symmetry, NNN hopping, akin to chiral hopping on the Haldane lattice, breaks time-reversal symmetry and leads to the formation of topological bands with Chern numbers ranging from 𝐶 = ±1 to ±4. We perform analytical and numerical analysis of the energy bands near the high-symmetry points Γ, ±𝐊, and 𝐌 𝑖 (𝑖 = 1, 2, and 3), which uncover a rich and complex landscape of HOVHS, controlled by the magnitude and phase of the NNN hopping. We observe power-law divergences in the density of states (DOS), 𝜌⁡(𝜀)∼|𝜀| −𝜈 , with exponents 𝜈 = 1/2, 1/3, 1/4, which can significantly affect the anomalous Hall response at low temperatures when the Fermi level crosses the HOVHS. Additionally, the NNN hopping induces the formation of higher Chern number bands 𝐶 = ±2, ±4 in the middle of the spectrum obeying a sublattice interference whereupon electronic states are maximally localized in each of the sublattices when the momentum approaches the three high-symmetry points 𝐌 𝑖 (𝑖 = 1, 2, and 3) on the Brillouin zone boundary. Finally, this classification of HOVHS in kagome systems provides a platform to explore unconventional electronic orders induced by electronic correlations.

Chern insulators↗
DOE OSTI2022

Optimal Sizing and Operation of a Hybrid Clean Energy Center

In this paper, we consider the deployment of battery energy storage systems (BESSs) paired with photovoltaics (PV) to form a hybrid clean energy center (HCEC), where the BESS is equipped for two main purposes: i) reducing the volatility of PV output and ii) shifting energy to meet an export limit. We develop a linear dynamical model to characterize the output of an HCEC and express the volatility and export power as functions of the battery's dispatch operations. We also propose an optimal sizing method via second-order cone programming (SOCP) to maximize the net benefits of the HCEC. Numerical studies are performed to illustrate and validate the proposed method.

Battery energy storage systems, photovoltaics, sec↗
DOE OSTI2020

Real-Time Modal Analysis of Electric Power Grids - The Need for Dynamic State Estimation

We articulate the reason why dynamic state estimation is needed to push the boundaries of the modal analysis of electric power grids in real-time operation. Then, we demonstrate how to unravel linear and nonlinear modes by using the extended dynamic mode decomposition along with estimates of the synchronous generators’ rotor angles and rotor speed deviations from nominal speed. The estimated modes are associated with electromechanical oscillations that take place continuously in electric power grids because of imbalances between power generation and demand. The numerical simulations are performed on a synthetic, albeit realistic, 2,000-bus network that was designed to resemble the electric power grid of Texas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Design, Optimization, and In-Depth Understanding of Front and Rear Junction Screen-Printed Double-Side Passivated Contacts Solar Cells

In this work, detailed numerical modeling is performed for front junction (FJ) and rear junction (RJ) n-type Si solar cells with screen-printed double-side poly-Si based tunnel oxide passivated contacts (TOPCon). Fundamental understanding of both structures is crucial to explore factors that limit the efficiency potential of the two designs. By investigating several key parameters such as front poly sheet resistance and thickness, bulk material properties, and carrier transport in our simulation model, we determine and explain why RJ cells outperform FJ cells. Our findings reveal that FJ suffers from present technological limitations of p-poly based passivated contacts, namely, i) large recombination observed in textured p-TOPCon emitter and ii) low solid solubility and hole mobility in boron-doped poly-Si resulting in very high sheet resistance of front p-poly layers which must be made thin to reduce absorption losses. In contrast, RJ desensitized the cell efficiency to front sheet resistance allowing application of ultra-thin front n-poly layers and is therefore ideally suited for double-side TOPCon cells. Lastly, we have developed a roadmap for efficiency projections of commercial-type screen-printed RJ and FJ topologies reaching ~24% and >22.5% efficiencies respectively, quantifying and explaining the impact of technological innovations on the performance of each design.

14 SOLAR ENERGY↗
DOE OSTI2020

A Fast Flexibility-Driven Generation Portfolio Planning Method for Sustainable Power Systems

We are witnessing an acceleration in the uptake of renewable energy in power systems. Because of the associated variability and uncertainty of renewables, power systems need to have an adequate supply of flexibility to allow for suitable management of short-term operations. So far most of the work in this area has neglected how flexibility needs associated with renewables are fulfilled as part of dispatchable generation capital investments decisions. To address this challenge, we propose an approach to plan the dispatchable generation mix of a power system as needed to counteract variability and uncertainty associated with significant shares of variable renewable generation. The approach exploits the linear time-invariant feature of variable generation variability using historical phase planes of capacity (in MW) and ramp (in MW/min) to bridge the gap between long-term capacity planning and short-term intra-hour flexibility needs. This approach is much more computationally tractable than other proposals, while also being able to capture adequately short-term operational features like ramping and net load variability. Numerical tests are performed on realistic datasets to substantiate the effectiveness of the approach.

bulk power system planning↗