Search NASA⌕ Search

SEARCH · Search NASA

Results for “first principles simulations”

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

First-principles theory of the rate of magnetic reconnection in magnetospheric and solar plasmas

Abstract The rate of magnetic reconnection is of the utmost importance in a variety of processes because it controls, for example, the rate energy is released in solar flares, the speed of the Dungey convection cycle in Earth’s magnetosphere, and the energy release rate in harmful geomagnetic substorms. It is known from numerical simulations and satellite observations that the rate is approximately 0.1 in normalized units, but despite years of effort, a full theoretical prediction has not been obtained. Here, we present a first-principles theory for the reconnection rate in non-relativistic electron-ion collisionless plasmas, and show that the same prediction explains why Sweet-Parker reconnection is considerably slower. The key consideration of this analysis is the pressure at the reconnection site (i.e., the x-line). We show that the Hall electromagnetic fields in antiparallel reconnection cause an energy void, equivalently a pressure depletion, at the x-line, so the reconnection exhaust opens out, enabling the fast rate of 0.1. If the energy can reach the x-line to replenish the pressure, the exhaust does not open out. In addition to heliospheric applications, these results are expected to impact reconnection studies in planetary magnetospheres, magnetically confined fusion devices, and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dynamic separation of gases using microsieves

Separation of light weight molecules, such as nitrogen, argon, and oxygen, from heavier compounds can have significant impacts on energy capture, environmental monitoring, or isotopic applications. Large-scale gas separation techniques, like gas centrifugation and membrane mitigation, can be problematic as they impart tremendous energy and induce high mechanical stress onto the instrumentation. Microsieves, also known as micronozzles or microfunnels, are developed to create physical barriers to separate specific isotopes and gases. Separation is achieved using a converging and diverging micronozzle to impose supersonic gas flow around a curved wall, and it has been used for the separation of heavy actinide isotopes in low weight gas as well as separation of low weight gas compositions of nitrogen and argon back in 1900s. However, systematic reviews of this unique technology are lacking. The application of the Laval style nozzle, which has a converging/diverging entrance fundamental to the micronozzle, is included in this review due to its importance in industrial applications in uranium (U) isotope refinement. Using advanced computational fluid dynamic (CFD) simulations, the extent of gas separation can be modelled. Herein, we first examine the literature and survey recent advances on fabrication techniques for creating curved micronozzles, methods and separation principles used to design devices. Furthermore, we then follow with highlights of CFD simulations applied to evaluate the separation effects using microsieves. Finally, identification of the gap and recommendation for future development and applications are suggested for using intrinsic molecular features and fluidic dynamics in formulating separation strategies.

30 Microfluidics↗
DOE OSTI2024

Sparsified time-dependent Fourier neural operators for fusion simulations

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗
DOE OSTI2020

Linking electronic structure calculations to generalized stacking fault energies in multicomponent alloys

Abstract The generalized stacking fault energy is a key ingredient to mesoscale models of dislocations. Here we develop an approach to quantify the dependence of generalized stacking fault energies on the degree of chemical disorder in multicomponent alloys. We introduce the notion of a “configurationally-resolved planar fault” (CRPF) energy and extend the cluster expansion method from alloy theory to express the CRPF as a function of chemical occupation variables of sites surrounding the fault. We apply the approach to explore the composition and temperature dependence of the unstable stacking fault energy (USF) in binary Mo–Nb alloys. First-principles calculations are used to parameterize a formation energy and CRPF cluster expansion. Monte Carlo simulations show that the distribution of USF energies is significantly affected by chemical composition and temperature. The formalism is broadly applicable to arbitrary crystal structures and alloy chemistries and will enable the development of rigorous models for deformation mechanisms in high-entropy alloys.

Natarajan, Anirudh Raju (ORCID:0000000202529480)↗
DOE OSTI2022

Negative capacitance regime in antiferroelectric PbZrO 3

The potential of antiferroelectrics to exhibit a negative capacitance regime has been largely overlooked as all the attention focused on their polar counterparts, ferroelectrics. We use nonequilibrium first-principles-based molecular dynamics to probe a negative capacitance regime in prototypical antiferroelectric PbZrO 3 . Simulations predict that this antiferroelectric can exhibit a negative capacitance/susceptibility regime in response of polarization to an internal electric field, which is a superposition of applied and residual depolarizing fields. Consequently, the regime emerges when the polarization surface charge in the polar phase of antiferroelectric is not fully screened, as is often the case in thin films and nanostructures. The negative capacitance regime occurs below the Curie temperature and disappears in the paraelectric phase. We find that the time the material spends in the negative capacitance regime is proportional to the time needed to complete antipolar–polar (or its reverse) transition and shortens as the frequency of the applied field increases. Furthermore, a negative susceptibility value exhibits strong dependence on the quality of surface charge screening with the largest in magnitude values occurring in the vicinity of the transition into a negative capacitance regime.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of nuclear quantum effects on the electronic properties of amorphous carbon

We carry out quantum simulations to study the physical properties of diamond-like amorphous carbon by coupling first-principles molecular dynamics with a quantum thermostat, and we analyze multiple samples representative of different defective sites present in the disordered network. We show that quantum vibronic coupling is critical in determining the electronic properties of the system, in particular its electronic and mobility gaps, while it has a moderate influence on the structural properties. We find that despite localized electronic states near the Fermi level, the quantum nature of the nuclear motion leads to a renormalization of the electronic gap surprisingly similar to that found in crystalline diamond. We also discuss the notable influence of nuclear quantum effects on band-like and variable-hopping mechanisms contributing to electrical conduction. Our calculations indicate that methods often used to evaluate electron–phonon coupling in ordered solids are inaccurate to study the electronic and transport properties of amorphous semiconductors composed of light atoms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Enhancing charge ratio sensitivity to hadronization effects via jet selections on resolved SoftDrop splitting

The study of quantum chromodynamics (QCD) at ultrarelativistic energies can be performed in a controlled environment through lepton-hadron deep inelastic scatterings. In such collisions, the high-energy partonic emissions that follow from the ejected hard partons are accurately described by perturbative QCD. However, the lower energy scales at which quarks and gluons experience color confinement, i.e., hadronization mechanism, fall outside the validity regions for perturbative calculations, requiring phenomenological models tuned to data to describe it. As such, hadronization physics cannot be currently derived from first principles alone. Monte Carlo event generators are useful tools to describe these processes as they simulate both the perturbative and the nonperturbative interactions, with model-dependent energy scales that control parton dynamics. This work employs jets—experimental reconstructions of final-state particles likely to have a common partonic origin—to inspect this transition further. Although originally proposed to circumvent hadronization effects, we show that jets can be utilized as probes of nonperturbative phenomena via their substructure. The charge correlation ratio was recently shown to be sensitive to hadronization effects. Our work further improves this sensitivity to nonperturbative scales by introducing a new selection based on the relative placement of the within the clustering tree, defined as the unclustering that resolves the jet’s leading charged particles. Published by the American Physical Society 2025

Apolinário, Liliana (ORCID:0000000335009681)↗
DOE OSTI2024

Bid-in Price-Sensitive Load: A Comparison with Fixed Load in a Realistic Test System

Bulk power systems operations, in both simulation and practice, are often built on the principle of minimizing production costs, which assumes to first order that demand is fixed and that the only degrees of freedom for balancing generation and load are supply-side unit commitment and dispatch decisions. The purpose of this report is to interrogate the value of welfare-maximizing markets that not only physically balance demand and supply but also ensure that market outcomes provide consumer value by directly maximizing welfare, defined as the value of consumption minus production costs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Learning Grain-Boundary Segregation: From First Principles to Polycrystals

The segregation of solute atoms at grain boundaries (GBs) can strongly impact the structural and functional properties of polycrystals. Yet, due to the limited availability of simulation tools to study polycrystals at the atomistic scale (i.e., interatomic potentials), there is a minimal understanding of the variation of solute segregation tendencies across the very complex space of GB microenvironments and the large range of alloys in which it can occur. Here, we develop an algorithmic framework that can directly learn the full spectrum of segregation energies for a metal solute atom in a metal polycrystal from ab initio methods, bypassing the need for alloy interatomic potentials. This framework offers a pathway to a comprehensive catalog of GB solute segregation with quantum accuracy, for the entire alloy space. As an initial demonstration in this pursuit, we build an extensive GB segregation database for aluminum-based alloys across the periodic table, including dozens of alloys for which there are substantially no prior data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Derivation of A Representative Elementary Volume (REV) for Upscaled Two-Phase Flow in Porous Media

Relative permeability plays an important role in the upscaling of multiphase flow in porous media from the pore scale to the Darcy scale. The entire concept of relative permeability is contingent on the existence of a representative elementary volume (REV). As we move to smaller samples to measure relative permeability, such as with digital core analysis, the concept of a classical REV has become increasingly unlikely when using the conventional approach to defining a representative volume. The “‘conventional”’ understanding of an REV is that a large enough volume must be considered such that spatial variability averages out. In digital rock methods, such as pore-scale simulations based on micro-computed tomography (CT) images, the domain size is typically 2 to 4 mm. This is approximately the length scale of a single-phase flow REV using the classic REV approach. However, the single-phase perspective does not consider the complex dynamics and fluctuations often observed in multiphase flow systems, even at centimeter-scale experiments and/or simulations. A fundamental question is, therefore, whether the domain size commonly used in digital rock simulations can provide a consistent energy budget such that the concept of relative permeability exists. Based on first principles, relative permeability accounts for the rate of energy dissipated in a stationary process. If the dynamics are fluctuating, the energy dissipated can vary but will average out over a long enough timescale. The key to determining the validity of the relative permeability is the timescale of the measurement, not the spatial scale. The conventional REV theory assumes that spatial, temporal, and ensemble averages are equivalent in an ergodic system, but it does not provide a way to test this assumption. We provide a formal way to identify the timescale where the relative permeability accurately captures energy dissipation as a way to validate relative permeability measurements and quantitatively assess their accuracy. This result will be tested for a practical SCAL test, determining how long a flow experiment needs to be run to accurately characterize the rate of energy dissipation by the flow. The outcome will be a best practice guide for the determination of relative permeability from core-scale experiments and/or digital core simulations that ensure the energy budget is fully accounted for in the relative permeability coefficient.

Mcclure, James [Virginia Tech, Blacksburg]↗
DOE OSTI2023

Spin dynamics of [1,2- 13 C 2 ]pyruvate hyperpolarization by parahydrogen in reversible exchange at micro Tesla fields

Hyperpolarization of 13 C-pyruvate via Signal Amplification By Reversible Exchange (SABRE) is an important recent discovery because of both the relative simplicity of hyperpolarization and the central biological relevance of pyruvate as a biomolecular probe for in vitro or in vivo studies. Here, we analyze the [1,2- 13 C 2 ]pyruvate-SABRE spin system and its field dependence theoretically and experimentally. We provide first-principles analysis of the governing 4-spin dihydride- 13 C 2 Hamiltonian and numerical spin dynamics simulations of the 7-spin dihydride- 13 C 2 –CH 3 system. Furthermore, the analytical and the numerical results are compared to matching systematic experiments. With these methods we unravel the observed spin state mixing of singlet states and triplet states at microTesla fields and we also analyze the dynamics during transfer from micro-Tesla field to high field for detection to understand the resulting spectra from the [1,2- 13 C 2 ]pyruvate-SABRE system.

13C pyruvate↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗
DOE OSTI2020

A First-principle Model for Polarization Swings during Reconnection-powered Flares

We show that magnetic reconnection in a magnetically dominated fast-cooling plasma can naturally produce bright flares accompanied by rotations in the synchrotron polarization vector. With particle-in-cell simulations of reconnection, we find that flares are powered by efficient particle acceleration at the interface of merging magnetic flux ropes, or “plasmoids.” The accelerated particles stream through the post-merger plasmoid toward the observer, thus progressively illuminating regions with varying plane-of-sky field direction, and so leading to a rotation in the observed polarization vector. Our results provide evidence for magnetic reconnection as the physical cause of highenergy flares from the relativistic jets of blazars (which recent observations have shown to be frequently associated with polarization rotations), and provide a first-principle physical mechanism for such flares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gating of Hydrophobic Nanopores with Large Anions

In this report, understanding ion transport in nanoporous materials is critical to a wide variety of energy and environmental technologies, ranging from ion-selective membranes, drug delivery, and biosensing, to ion batteries and supercapacitors. While nanoscale transport is often described by continuum models that rely on a point charge description for ions and a homogeneous dielectric medium for the solvent, here, we show that transport of aqueous solutions at a hydrophobic interface can be highly dependent on the size and hydration strength of the solvated ions. Specifically, measurements of ion current through single silicon nitride nanopores that contain a hydrophobic-hydrophilic junction show that transport properties are dependent not only on applied voltage but also the type of anion. We find that in Cl - -contained solutions the nanopores only conducted ionic current above a negative voltage threshold. On the other hand, introduction of large polarizable anions, such as Br - and I - , facilitated the pore wetting, making the pore conductive at all examined voltages. Molecular dynamics simulations revealed that the large anions, Br - and I - , have a weaker solvation shell compared to that of Cl - , and consequently were prone to migrate from the aqueous solution to the hydrophobic surface, leading to the anion accumulation responsible for pore wetting. The results are essential for designing nanoporous systems that are selective to ions of the same charge, for realization of ion-induced wetting in hydrophobic pores, as well as for fundamental understanding on the role of ion hydration shell on the properties of solid/liquid interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Thickness-Dependent Thermal Conductivity and Phonon Mean Free Path Distribution in Single-Crystalline Barium Titanate

Nanosized perovskite ferroelectrics are widely employed in several electromechanical, photonics, and thermoelectric applications. Scaling of ferroelectric materials entails a severe reduction in the lattice (phonon) thermal conductivity, particularly at sub-100 nm length scales. Such thermal conductivity reduction can be accurately predicted using the information of phonon mean free path (MFP) distribution. The current understanding of phonon MFP distribution in perovskite ferroelectrics is still inconclusive despite the critical thermal management implications. Here, high-quality single-crystalline barium titanate (BTO) thin films, a representative perovskite ferroelectric material, are grown at several thicknesses. Using experimental thermal conductivity measurements and first-principles based modeling (including four-phonon scattering), the phonon MFP distribution is determined in BTO. The simulation results agree with the measured thickness-dependent thermal conductivity. The results show that the phonons with sub-100 nm MFP dominate the thermal transport in BTO, and phonons with MFP exceeding 10 nm contribute ≈35% to the total thermal conductivity, in significant contrast to previously published experimental results. The experimentally validated phonon MFP distribution is consistent with the theoretical predictions of other complex crystals with strong anharmonicity. This work paves the way for thermal management in nanostructured and ferroelectric-domain-engineered systems for oxide perovskite-based functional materials.

36 MATERIALS SCIENCE↗