SEARCH · Search NASA

Results for “core 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.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Strong Effect of Nonpolar Solvent Molecular Structure on CdSe Nanoplatelet Stacking

We report a drastic difference in stacking behavior of oleic acid-stabilized 4-monolayer (4 ML) CdSe nanoplatelets (NPLs) in toluene and methylcyclohexane (MCH), two nonpolar solvents that differ in the conformational flexibility of their molecules. Using liquid cell transmission electron microscopy (TEM) and small angle scattering (SAXS) techniques, we show that NPLs form microns-long ribbons consisting of 4 ML CdSe NPLs in toluene, the solvent widely used to form stable colloidal solutions of a broad range of quasi-spherical nanoparticles. In contrast, 4 ML CdSe NPLs are well dispersed in MCH. The difference in stacking behavior of NPLs in toluene and MCH suggests that the conformational flexibility of the solvent molecules, such as the ability to adopt multiple chair conformations, modulates nanoplatelet interactions. Molecular dynamics (MD) simulations reveal that solvent molecules subtly alter the structure of the organic ligand shell. These solvent-dependent changes propagate to the inorganic core, modulating the degree of CdSe nanoplatelet (NPL) twisting and, consequently, the properties of the nanoparticles. We show that toluene better solvates oleate ligands while MCH induces a bimodal oleate span distribution, which can lead to increased solubility of CdSe NPLs. In addition, the solvent can also influence the inorganic core, which, in turn, can modify the nanoparticle properties. We demonstrate that destabilization of toluene solution containing ribbons of 4 ML CdSe NPLs without CdS shells results in the formation of NPL assemblies with amplified spontaneous emission (ASE) with a low threshold of 14 µJ cm−2 that is comparable with that of CdSe/CdS core/shell NPLs. Our results emphasize that the solvent plays a major role in mediating interactions between NPLs and hence their processability for fabrication of functional structures.

CdSe↗

DOE OSTI2025

Significantly enhanced near-field coupling via tip engineering

The ability to significantly enhance near-field coupling between light and matter at the nanoscale is crucial for advancing the fields of nanophotonics and nanopolariotonics. However, conventional probes face challenges in achieving optimal light–matter interaction. In this study, we propose a novel, to the best of our knowledge, simulation-based strategy that leverages tip engineering to dramatically amplify the scattering field through tailored double-layer geometries. By employing a core-shell structure with a thin shell layer optimized for specific dielectric permittivity and effective polarizability, we demonstrate a near-field enhancement of up to 10 times compared to conventional probes. Our findings highlight exciting new possibilities for optimizing near-field interactions through probe designs with customized resonances, paving the way for substantially improved nano-optical sensing, imaging, and detection.

Shiravi, H.↗

DOE OSTI2025

Probing autoionization decay lifetimes of the 4d−16ℓ core-excited states in xenon using attosecond noncollinear four-wave-mixing spectroscopy

The decay of core-excited states is a sensitive probe of autoionization dynamics and correlation effects in many-electron systems, occurring on ultrafast timescales. Xenon, with its dense manifold of autoionizing resonances that can be coupled with near-infrared light, provides a platform to investigate these processes. In this work, the autoionization decay lifetimes of 4d-16ℓ (ℓ = s, p, d, …) core-excited states in xenon atoms are probed with extreme ultraviolet (XUV) attosecond noncollinear four-wave-mixing (FWM) spectroscopy. The 4d{5/2,3/2}-16p XUV-bright states (optically dipole-allowed) exhibit decay lifetimes of ∼6 fs, which is consistent with spectator-type decay. In contrast, the 4d{5/2,3/2}-16s and 4d{5/2,3/2}-16d XUV-dark states (optically dipole forbidden) show longer decay lifetimes of ∼20 fs. Photoionization calculations confirm that all core-hole states with 4d character should decay via spectator channels in ≤6 fs, suggesting that the apparently longer dark-state decay times arise from an alternative mechanism. A few-level simulation of the FWM process shows that the inclusion of a nearby, longer-lived dark state can mimic the experimental FWM signal, suggesting population cycling with a second electronic state with non-4d character. Ab initio calculations support the presence of such multielectron excited states in the 60-70 eV range. These results demonstrate that FWM signals can encode coupled-state dynamics when probing complex systems, highlighting the importance of combining theoretical and experimental approaches to disentangle accurate core-level decay pathways and lifetimes.

Puskar, Nicolette G↗

DOE OSTI2025

High temporal frequency data from a four turbine, blade-resolved wind farm simulation with ExaWind

The data was generated with ExaWind (https://github.com/Exawind) which couples AMR-Wind (https://github.com/Exawind/amr-wind/), Nalu-Wind (https://github.com/Exawind/nalu-wind), TIOGA (https://github.com/Exawind/tioga), and OpenFAST (https://github.com/OpenFAST/openfast). This is a large-scale simulation of a blade-resolved wind farm using the ExaWind software stack. ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. Another application, OpenFAST, handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. This particular simulation includes four blade-resolved wind turbines operating in a turbulent atmospheric boundary layer. The AMR-Wind solver uses 500 million cells and is being solved on 256 AMD GPUs of the Oakridge Leadership Computing Facility Frontier supercomputer. Each turbine is assigned its own Nalu-Wind solver with over 13 million elements per turbine and solved using 448 CPU cores, for a total of 1792 CPU cores. For each node, 56 cores contain Nalu-Wind, while 8 cores correspond to AMR-Wind operations on the GPUs. Consequently, ExaWind is entirely utilizing the CPUs and the GPUs of the nodes concurrently. The data used in the visualization is full flow field data output from the simulation. It is lossy-compressed to a specific accuracy using ZFP and written to disk every 16 time-steps to enable real-time flow visualization. The flow fields are sampled at a high temporal frequency to enable real-time, 24fps visualization. The flow fields are sampled every 12 simulation time steps (every 0.04132s).

17 WIND ENERGY↗

DOE OSTI2025

Scaling Ultrahigh-Resolution E3SM Land Model for Leadership-Class Supercomputers

This paper presents advancements in scaling the ultrahigh-resolution E3SM Land Model (uELM) for deployment on leadership-class supercomputers, addressing the increased demand for km-scale Earth system modeling. By focusing on km-scale ELM simulations, we enhance predictive capabilities for climate interactions, facilitating improved responses to climate change impacts on energy systems, agriculture, and water resources. Our approach leverages innovative software architecture optimizations, sophisticated data handling techniques, and advanced parallel processing, achieving strong scalability on two leadership supercomputers (2400 nodes (105,600 cores) on Summit, and 1200 nodes (76,800 cores) on Frontier). Results from extensive scalability assessments on the Summit and Frontier also demonstrate outstanding I/O performance (close to 400 GB/s write throughput) and the model's ability to efficiently handle increasing computational demands. This study not only establishes uELM's capability for high-resolution simulations over vast geographical domains, but also sets a foundation for future Earth system modeling breakthroughs.

Wang, Dali [ORNL] (ORCID:0000000168065108)↗

DOE OSTI2024

Simulation of metal nanoparticles growth in methane atmosphere of arc discharge: comparison to experiment

A direct current arc discharge in a methane atmosphere is a scalable and sustainable method to produce metal-carbon core–shell nanoparticles and single-walled carbon nanotubes, where a metal catalyst can be continuously supplied through evaporation of an anode made from the catalyst material. The size of catalyst particles is of critical importance as it can affect the synthesis yield and properties of nanotubes and core–shell nanoparticles. This study presents a numerical model describing the formation and growth of metal particles for the conditions representative of the arc discharge with an evaporating iron anode at near-atmospheric pressure of a methane-rich atmosphere. The model incorporates carbon adsorption to the metal surface and explains the limiting effect of carbon coverage on the size of metal nanoparticles. The predicted particle sizes are compared with experimental observations. The model also predicts higher concentrations of metal particles with the increasing partial pressure of methane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2024

The Core-collapse Supernova Sensitivity Machine (Final Technical Report)

Core-collapse supernovae are vital phenomena for understanding the extreme physics of stellar evolution, nuclear interactions, and the origins of elements in the universe. This project developed high-resolution 3D simulations and innovative 1D turbulence models to study the explosion mechanisms, yielding critical insights into the role of turbulence, magnetohydrodynamics, and neutrino physics. The outcomes include enhanced sensitivity analyses connecting nuclear physics inputs to observable supernovae characteristics, tools for the broader research community, and substantial contributions to the training of future scientists.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Significant Efficiency Enhancements in Non‐Y Series Acceptors by the Addition of Outer Side Chains

Abstract Most current highly efficient organic solar cells utilize small molecules like Y6 and its derivatives as electron acceptors in the photoactive layer. In this work, a small molecule acceptor, SC8‐IT4F, is developed through outer side chain engineering on the terminal thiophene of a conjugated 6,12‐dihydro‐dithienoindeno[2,3‐d:2′,3′‐d′]‐s‐indaceno[1,2‐b:5,6‐b′]dithiophene (IDTT) central core. Compared to the reference molecule C8‐IT4F, which lacks outer side chains, SC8‐IT4F displays notable differences in molecule geometry (as shown by simulations), thermal behavior, single‐crystal packing, and film morphology. Blend films of SC8‐IT4F and the polymer donor PM6 exhibit larger carrier mobilities, longer carrier lifetimes, and reduced recombination compared to C8‐IT4F, resulting in improved device performance. Binary photovoltaic devices based on the PM6:SC8‐IT4F films reveal an optimal efficiency over 15%, which is one of the best values for non‐Y type small molecule acceptors (SMAs). The resultant devices also show better thermal and operational stability than the control PM6:L8‐BO devices. SC8‐IT4F and its blend exhibit a higher relative degree of crystallinity and π coherence length, compared to C8‐IT4F samples, beneficial for charge transport and device performance. The results indicate that outer side chain engineering on existing small electron acceptors can be a promising molecular design strategy for further pursuing high‐performance organic solar cells.

He, Qiao [Department of Chemistry and Centre for P↗

DOE OSTI2024

An Open-Source Coupling for Depletion During Fuel Cycle Modeling

Fuel depletion is an important aspect of fuel cycle modeling, allowing a user to account for how loaded fuel compositions affect in-core and spent fuel compositions and their related fuel cycle metrics. Therefore, multiple methods have been developed to account for depletion within fuel cycle simulations. This work adds to that list of methods by introducing an open-source coupling between C$\scriptsize{YCLUS}$ and OpenMC to perform fuel depletion during a fuel cycle simulation, called OpenMCyclus. This work explains the methodology of OpenMCyclus and presents a benchmark comparison between the performance of OpenMCyclus and another C$\scriptsize{YCLUS}$ archetype that uses recipes to define spent fuel compositions. In conclusion, the development of this coupling expands the functionalities possible through C$\scriptsize{YCLUS}$ by providing real-time fuel depletion that is reactor agnostic and open source.

C$\scriptsize{YCLUS}$↗

DOE OSTI2025

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic↗

DOE OSTI2026

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

First tungsten radiation studies in DIII-D’s ITER baseline demonstration discharges

Abstract ITER Baseline Scenario plasmas were studied in DIII-D using krypton and xenon gases as a proxy for the tungsten that will be present in ITER. These impurities were chosen for having the same radiative loss rate Lz as tungsten would exhibit in the hotter ITER core. Results show that the scenario with these core radiators spans the range of impurity concentration and W radiated fraction expected for ITER, and up to 50% higher values, explored at zero injected torque, as well as 1 Nm and full co-torque injection with T ∼ 3 Nm. Stationary discharges with duration >2–4 τ R are achieved with f rad ⩾ 30% leading to a reduction in confinement of ∼10%, and a comparison with real metal radiators in the same range of f rad shows that the higher Lz at the lower temperatures in these plasmas yields too pessimistic results on the survivability and performance of this scenario in ITER. Simulations of ITER power balance including W radiation show that with concentration up to three times higher than in the DIII-D plasmas the scenario can be stationary, remaining at acceptable core radiated fraction values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Cross slip of extended dislocations in face-centered cubic metals through phase-field modeling

Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.

36 MATERIALS SCIENCE↗

DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Optical diagnostic design for measuring the radiation front with a mid-leg pumped divertor on DIII-D

The design of an optical diagnostic system to localize the radiation front in a mid-leg pumped divertor configuration on DIII-D is presented. Divertor detachment is a key mechanism for handling power exhaust in tokamaks, and mid-leg pumping offers a promising approach to radiatively dissipate power while maintaining acceptable core performance. To predict the location of the radiation front and design a spectroscopic diagnostic to measure its position, a database of SOLPS-ITER simulations across a range of input powers and gas puffing rates representative of DIII-D operation was employed. These simulations provide self-consistent plasma backgrounds for a Cherab-Raysect synthetic diagnostic framework, which incorporates detailed tokamak geometry and physically accurate ray-tracing. Synthetic simulations of optical sightlines and viewing cone geometry were used to calculate line-integrated emission from the plasma, particularly of the C III 465 nm line, to serve as a proxy for the T e ≈ 7 − 10 eV temperature region associated with the onset and evolution of divertor detachment. Neutral deuterium emission is also evaluated for comparison. The synthetic diagnostic is used to assess the performance of existing DIII-D optical systems, including filterscopes and the Multichord Divertor Spectrometer, and to optimize line-of-sight placement within mechanical and installation constraints. The results provide quantitative guidance for diagnostic implementation on DIII-D and demonstrate the advantage of integrated synthetic diagnostics for divertor design studies and future advanced divertor concepts.

Cherab↗

DOE OSTI2026

Probing the Mechanism of Cadmium Sulfide Cluster Nucleation and Growth via Sequential Infiltration Synthesis

Few-atom metal chalcogenide clusters may be realized through the sequential infiltration of metal–organic precursors in polymer films. However, the underlying nucleation and growth mechanisms that allow for cluster synthesis with near atomic-scale precision are not fully resolved. The kinetics of the sequential infiltration synthesis (SIS) method that control the nucleation and growth mechanisms of primarily Cd 4 S 4 -core clusters within a poly(4-vinylpyridine) (P4VP) matrix are probed with in situ UV–visible absorbance spectroscopy. Density functional theory (DFT) calculations that allow simulation of the optical properties of cluster fragments further reveal the thermodynamics that guide cluster growth within the P4VP matrix. Here, we conclude that a reactive capture mechanism for cluster nucleation and growth is dominant, although transient dimethyl cadmium adduction to the polymer backbone may contribute to cluster nucleation under shorter metal–organic purge process conditions. Grazing incidence X-ray diffraction (GI-XRD) and X-ray absorption spectroscopy (XAS) analyses further corroborate the cluster size and atom connectivity throughout the stepwise synthesis.

Jayaweera, Nuwanthaka P. [Argonne National Laborat↗

DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗