Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Design, Analysis, and Experimental Testing of Hydrogen Lean Direct Injection Nozzles at Elevated Pressure

Abstract There are many challenges of commissioning a hydrogen combustor into future gas turbine engines; especially regarding achieving emissions goals. Previously, Escudero et al. and Tran et al. conducted a study to adapt the liquid fuel Lean Direct Injection (LDI) concept from Jet-A to gaseous natural gas-hydrogen blends and pure hydrogen [1], [2]. Experimental data was collected at atmospheric conditions using a Box Behnken design of experiments. The design of experiments suggested that biasing the air split in favor of the inner air circuit and increasing the swirl strength of this inner air passage resulted in improved NOx emissions, while the inverse was true for stability, which was quantified by studying the lean blowoff point (LBO) [1], [2]. The trends revealed by the original experiment [1], [2] provided a design direction for further iterations of the experimental hardware. The study presented herein describes the further investigation of such LDI injectors through experimental methods and computational fluid dynamic (CFD) simulations at atmospheric conditions, which were used to identify potential flow behaviors driving enhanced emissions performance. Further evaluation of select injectors from both studies was then conducted at elevated pressures up to 6 atmospheres. The results from both experiments are presented in this study, which include flame observations, emissions measurements, and operational challenges. NOx emissions results are reported on a volume basis in ppmvd corrected to 15% O2 and corrected for fuel. A predictive model for relating NOx emissions to test conditions at atmospheric conditions show high significance to adiabatic flame temperature while little to no significance to fuel composition for the best performing configurations. The results illustrate the connection between atmospheric testing and testing elevated pressures. The design direction indicated by the initial tests and CFD results in promising configurations for implementation into a Multi-point LDI array.

08 HYDROGEN↗
DOE OSTI2022

Atomically engineered cobaltite layers for robust ferromagnetism

Emergent phenomena at heterointerfaces are directly associated with the bonding geometry of adjacent layers. Effective control of accessible parameters, such as the bond length and bonding angles, offers an elegant method to tailor competing energies of the electronic and magnetic ground states. In this study, we construct unit-thick syntactic layers of cobaltites within a strongly tilted octahedral matrix via atomically precise synthesis. The octahedral tilt patterns of adjacent layers propagate into cobaltites, leading to a continuation of octahedral tilting while maintaining substantial misfit tensile strain. These effects induce severe rumpling within an atomic plane of neighboring layers, further triggering the electronic reconstruction between the splitting orbitals. First-principles calculations reveal that the cobalt ions transit to a higher spin state level upon octahedral tilting, resulting in robust ferromagnetism in ultrathin cobaltites. This work demonstrates a design methodology for fine-tuning the lattice and spin degrees of freedom in correlated quantum heterostructures by exploiting epitaxial geometric engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

NFE-24-10417: Durability and Life-time Energy Performance Enhancement of Vapor Compression Systems

This project investigated the use of a Centrifugal Particle Separation (CPS) system to improve the durability and operational reliability of outdoor heat exchangers in vapor-compression based heating, ventilation, air conditioning, and refrigeration (HVACR) systems. The developed method introduced a novel dust and fouling mitigation approach by integrating an ambient air pre-cleaner that employs centrifugal forces to remove particulate matter from incoming air before it reaches the heat exchanger. By reducing the accumulation of dust and debris on coil surfaces, the CPS system aims to lower maintenance requirements and extend equipment lifespan. The project specifically involves coupling a centrifugal air pre-cleaner with a single-zone ductless mini-split system and evaluating its effectiveness in maintaining cleaner coil surfaces. Additionally, the study examined the impact of CPS integration on system energy performance. Through controlled experimental testing and performance monitoring, this work seeks to validate CPS-based air pre-cleaning as a practical, efficient solution for mitigating outdoor coil dust accumulation and fouling in both residential and commercial HVACR applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗
NASA NTRS1989

Advances in fluid dynamics

Various papers on advances in fluid dynamics are presented. Individual topics addressed include: nonunique solutions of the Euler equations; numerical solution of the biharmonic equation using a spectral multigrid method; finite-volume shape and placement for efficiency, accuracy, and physics; NND schemes and their application to numerical simulation of complex plume flow problems; symbolic analysis and computational algorithm construction; and velocity and mass flux distribution measurements of spherical glass beads in air flow in a 90-deg vertical-to-horizontal bend. Also discussed are: hyperbolic solution adaptive mesh generation; progress and future directions in CFD; numerical simulation of shock focusing; microbubble formation and splitting in a turbulent boundary layer for turbulence reduction; numerical study of spherical Taylor-Couette flow, trajectories of vortex lines beneath a free surface or above a plane; and a mathematical model for response of the coronary circulation to high sustained gravitational force fields.

Ballhaus, W. F., Jr.↗
NASA NTRS2024

Spatiotemporal Evolution of Marine Heatwaves Globally

The spatiotemporal evolution of marine heatwaves (MHWs) is explored using a tracking algorithm termed Ocetrac that provides objective characterization of MHW spatiotemporal evolution. Candidate MHW grid points are defined in detrended gridded sea temperature data using a seasonally varying temperature threshold. Identified MHW points are collected into spatially distinct objects using edge detection with weak sensitivity to edge detection and size threshold criteria. These MHW objects are followed in space and time while allowing objects to split and merge. Ocetrac is applied to monthly satellite sea surface temperature data from September 1981 through January 2021. The resulting MHWs are characterized by their intensity, duration, and total area covered. The global analysis shows that MHWs in the Gulf of Maine and Mediterranean Sea evolve within a relatively small region, while major MHWs in the Pacific and Indian Oceans are linked in space and time. The largest and most long lasting MHW using this method lasts for 60 months from November 2013 to October 2018, encompassing previously identified MHW events including those in the Northeast Pacific (2014-2015), the Tasman Sea (2015-2016, 2017-2018), and the Great Barrier Reef (2016).

spatiotemporal↗
DOE OSTI2017

Pressure Distortion of the H{sub 2}–He Collision-induced Absorption at the Photosphere of Cool White Dwarf Stars

Collision-induced absorption (CIA) from molecular hydrogen is a dominant opacity source in the atmosphere of cool white dwarfs. It results in a significant flux depletion in the near-IR and IR parts of their spectra. Because of the extreme conditions of helium-rich atmospheres (where the density can be as high as a few g cm{sup −3}), this opacity source is expected to undergo strong pressure distortion and the currently used opacities have not been validated at such extreme conditions. To check the distortion of the CIA opacity, we applied state-of-the-art ab initio methods of computational quantum chemistry to simulate the CIA opacity at high densities. The results show that the CIA profiles are significantly distorted above densities of 0.1 g cm{sup −3} in a way that is not captured by the existing models. The roto-translational band is enhanced and shifted to higher frequencies as an effect of the decrease of the interatomic separation of the H{sub 2} molecule. The vibrational band is blueward shifted and split into Q {sub R} and Q {sub P} branches, separated by a pronounced interference dip. Its intensity is also substantially reduced. The distortions result in a shift of the maximum of the absorption from 2.3 μm to 3–7 μm, which could potentially explain the spectra of some very cool, helium-rich white dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Salt Vessel-Sample Generator Interface

The goal of the salt sampling program at Argonne is to develop and deploy automated molten salt sampling approaches for interfacing relevant unit operations with salt analysis to improve the timeliness and accuracy of sampling-based accountancy measurements. One technology under development is a vacuum sampler loop module. In this system, molten salt is drawn from a process vessel through a heated sampling line and into a charge vessel. Next, one or more samples are split from the volume in the charge vessel and the remaining salt is returned to the process. The vacuum sampler loop module is being developed to interface with other sample processing and analysis modules to enable rapid at-line sample characterization. The main purpose of the vacuum sampling loop module is to enable the collection of samples that are more representative of the bulk salt by replacing traditional point samplers (i.e., dip probes) with a sampling approach that captures a larger cross section of salt. Additionally, the vacuum sampling approach eliminates the risk of dross contamination of samples and avoids the use of moving parts in the salt. Two methods of interfacing the vacuum sampling loop with a precision pneumatic sample generator were investigated in FY21. This report covers the testing of fluidic coupling between the two modules. Two iterations of the fluidically coupled modules were tested. The first iteration system coupled the two modules using a freeze valve to seal the vacuum sampler during filling and to control flow into the pneumatic sample generator. While this integrated system functioned as intended, some changes were implemented to make the system more robust and better suited to remote deployment. Specifically, the system was made to be more modular and active control of the vacuum filling operation was replaced with a passive control mechanism. For passive filling, the salt charge vessel was vented to a small gas tank that was at negative pressure, causing salt to be drawn into the vessel until the force of the fluid head was in equilibrium with the gas pressure. The passive control system performed well and will be used in future iterations. Another change in the second system was a newly configured pneumatic sample generator in which sample ejection occurred through a hole in the reservoir’s stainless-steel side wall instead of through a non-wetted sapphire orifice on the bottom. This alternate configuration may be better suited for near-process deployment because it enables on-line orifice maintenance and an orifice bypass drain back to the process. A third change in the second iteration system was the transition to a two-chamber charge vessel which split off a fraction of the sampled salt as a liquid aliquot. The goal was to create a buffer mechanism that would allow reproducible aliquoting of samples, independently of variability in the charge vessel fill height. While the two-chamber design was functional, the sample size reproducibility was below target values. To improve reproducibility and overcome many of the impediments to remote deployment of the vacuum sampling loop module, separate work was conducted to replace the two-chamber liquid aliquoting mechanism with aliquoting into single-use sample tubes. Solid salt transfer in the sample tubes will replace fluidic coupling for integrating the vacuum sampler with downstream modules. Because the proposed operations can all be executed with simple overhead actuation mechanisms or other existing hot-cell technology, there will be no need for large investments in novel hot cell sample handling technologies using this alternate approach. As such, near-term deployment of the vacuum sampling technology will be achievable. This new approach for automated coupling of sample tubes with down-stream modules is covered in a separate FY21 report, and a remotely operated version of the vacuum sampler loop module with the tube aliquoting feature is planned for FY22.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

On the Energy of a “One-Dimensional” Two-Electron Atom

Based on the perturbation theory and variational method long known for a “three-dimensional” atom, the ground and first excited state energies are calculated for a “one-dimensional” two-electron atom in the “one-dimensional ortho-helium” configuration, which can be obtained experimentally in principle, as has been already done for a Na Bose condensate, or produced in a super strong magnetic field B ≫ (2α){sup 2}B{sub 0} (B{sub 0} = m{sup 2}c{sup 3}/eħ ≈ 4.41 × 10{sup 13} G). The “screening constant” σ for this atom in the ground and excited states was about 0.20 and 0.17, 0.18, respectively, depending on the relative parity PP' of the electronic states, which is somewhat smaller than in “two-dimensional” and “three-dimensional” variants (in these cases, this constant in the ground state is almost the same and about 0.3). The frequencies of the main spectral lines of a “onedimensional” He atom representing a doublet split over the relative parity PP' are found. The presence of the close lines of this doublet in the emission spectrum of magnetars at frequencies ω{sub 1,2} ≈ {1.15; 1.17}α{sup 2}(c/λC) (α = e{sup 2}/ħc, λ{sub C} =ħ/mc) corresponding to the “one-dimensional ortho-helium” would suggest the existence of a superstrong magnetic field in such astrophysical objects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Four‐Coordinate Pr 4+ Imidophosphorane Complex

Abstract The imidophosphorane ligand, [NP t Bu 3 ] − ( t Bu= tert ‐butyl), enables isolation of a pseudo‐tetrahedral, tetravalent praseodymium complex, [Pr 4+ (NP t Bu 3 ) 4 ] ( 1‐Pr ), which is characterized by a suite of physical characterization methods including single‐crystal X‐ray diffraction, electron paramagnetic resonance, and L 3 ‐edge X‐ray near‐edge spectroscopies. Variable‐temperature direct‐current magnetic susceptibility data, supported by multiconfigurational quantum chemical calculations, demonstrate that the electronic structure diverges from the isoelectronic Ce 3+ analogue, driven by increased crystal field. The four‐coordinate environment around Pr 4+ in 1‐Pr , which is unparalleled in reported extended solid systems, provides a unique opportunity to study the interplay between crystal field splitting and spin‐orbit coupling in a molecular tetravalent lanthanide within a pseudo‐tetrahedral coordination geometry.

Boggiano, Andrew C.↗
DOE OSTI2024

A Four‐Coordinate Pr 4+ Imidophosphorane Complex

Abstract The imidophosphorane ligand, [NP t Bu 3 ] − ( t Bu= tert ‐butyl), enables isolation of a pseudo‐tetrahedral, tetravalent praseodymium complex, [Pr 4+ (NP t Bu 3 ) 4 ] ( 1‐Pr ), which is characterized by a suite of physical characterization methods including single‐crystal X‐ray diffraction, electron paramagnetic resonance, and L 3 ‐edge X‐ray near‐edge spectroscopies. Variable‐temperature direct‐current magnetic susceptibility data, supported by multiconfigurational quantum chemical calculations, demonstrate that the electronic structure diverges from the isoelectronic Ce 3+ analogue, driven by increased crystal field. The four‐coordinate environment around Pr 4+ in 1‐Pr , which is unparalleled in reported extended solid systems, provides a unique opportunity to study the interplay between crystal field splitting and spin‐orbit coupling in a molecular tetravalent lanthanide within a pseudo‐tetrahedral coordination geometry.

Boggiano, Andrew C.↗
DOE OSTI2018

Mössbauer Analysis and Cation Distribution of Zn Substituted BaFe{sub 12}O{sub 19} Hexaferrites

Barium hexaferrite is a well-known hard magnetic material. Doping using nonmagnetic cation such as Zn{sup 2+} were found to enhance magnetization owing to preferential tetrahedral site (4 f{sub 1}) occupancy of the zinc. However, the distribution of cations in hexaferrites depends on many factors such as the method of preparation, nature of the cation, and chemical composition. Here, Zn-doped barium hexaferrites (Ba{sub 1−x}Zn{sub x}Fe{sub 12}O{sub 19}) were synthesized by sol-gel method. In this study, we summarized the magnetic properties of Ba{sub 1−x}Zn{sub x}Fe{sub 12}O{sub 19} (x = 0, 0.1, 0.2, 0.3) BaM, investigated by Mössbauer spectroscopy. Moreover, cation distribution was also calculated for all the products. Mössbauer parameters were determined from 57Fe Mössbauer spectroscopy and according to it, the replacement of Ba-Zn affects all parameters such as isomer shift, the variation in line width, hyperfine magnetic field, and quadrupole splitting. Cation distribution revealed the relative area of undoped BaM, 12k, 2a, and 4 f{sub 2} positions which are close to theoretical values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Valley splitting of single-electron Si MOS quantum dots

Here, silicon-based metal-oxide-semiconductor quantum dots are prominent candidates for high-fidelity, manufacturable qubits. Due to silicon's band structure, additional low-energy states persist in these devices, presenting both challenges and opportunities. Although the physics governing these valley states has been the subject of intense study, quantitative agreement between experiment and theory remains elusive. Here, we present data from an experiment probing the valley states of quantum dot devices and develop a theory that is in quantitative agreement with both this and a recently reported experiment. Through sampling millions of realistic cases of interface roughness, our method provides evidence that the valley physics between the two samples is essentially the same.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Density functional methods for the magnetism of transition metals: SCAN in relation to other functionals

We report tests of various density functionals for ferromagnetic Fe, Co, and Ni with a focus on characterizing the behavior of the so-called strongly constrained and appropriately normed (SCAN) functional. It is found that SCAN is closer in behavior to functionals that yield localized behavior, such as hybrid functionals, than other semilocal functionals that are tested. The results are understood in terms of a tendency to differentiate orbitals, favoring integer occupation, which is necessary for a correct description of atomic systems, but inappropriate for the open shell metallic ferromagnetic metals studied here. Not only is the exchange splitting for open shells enhanced with SCAN, but also as seen in Ni, there is much more band dependence with significantly more overestimation for bands corresponding to the partially filled orbitals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dimension-free path-integral molecular dynamics without preconditioning

Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics—including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods—is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann–Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of “dimension-free” numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi–Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

Korol, Roman (ORCID:0000000193076351)↗