Search NASA⌕ Search

SEARCH · Search NASA

Results for “Applied Computing, Transportation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Thermodynamic and electron transport properties of Ca 3 Ru 2 O 7 from first-principles phonon calculations and Boltzmann transport theory

Here, this work demonstrates a first-principles-based approach to obtaining finite temperature thermal and electronic transport properties which can be employed to model and understand mesoscale structural evolution during electronic, magnetic, and structural phase transitions. A computationally tractable model was introduced to estimate electron relaxation time and its temperature dependence. The model is applied to Ca 3 Ru 2 O 7 with a focus on understanding its electrical resistivity across the electronic phase transition at 48 K. A quasiharmonic phonon approach to the lattice vibrations was employed to account for thermal expansion while the Boltzmann transport theory including spin-orbit coupling was used to calculate the electron-transport properties, including the temperature dependence of electrical conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Prediction of defect properties in concentrated solid solutions using a Langmuir-like model

The alleged existence of sluggish diffusion in high-entropy alloys has drawn controversy. In high-entropy alloys and, in general, in all solids, transport properties are controlled by point defect concentration, which must be known before performing atomistic simulations to compute transport coefficients. In this work, we present a general Langmuir-like model for defect concentration in an arbitrarily complex solid solution and apply this model to generate expressions for concentrations of vacancies and small interstitial atoms. We then calculate the vacancy concentration as a function of temperature in the equiatomic CoNiCrFeMn and FeAl alloys with modified embedded-atom-method potentials for various chemical orderings, showing there is no clear correlation between vacancy thermodynamics and chemical ordering in the CoNiCrFeMn alloy, but clear systematic patterns for FeAl. We believe this is due to the high stability of disordered, random, and ordered intermetallic phases, respectively, in the CoNiCrFeMn and FeAl systems. Finally, this work provides future avenues to the prediction of thermal interstitials and vacancies in solid solutions, which is necessary for models of nonequilibrium behavior of solid solutions.

composition↗
DOE OSTI2025

Radiation GRMHD Models of Accretion onto Stellar-mass Black Holes. I. Survey of Eddington Ratios

We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Precise 3D reactor core calculation using spherical harmonics and discontinuous Galerkin finite element methods

We study the use of P{sub N} method in angle and discontinuous Galerkin is space to solve 3D neutron transport problem. P{sub N} method consists in developing the angular flux on truncated spherical harmonics basic. In this paper, we couple this method with the discontinuous finite elements in space to obtain a complete discretization of the multigroup neutron transport equation. To investigate its precision, the method was applied to Takeda and C5G7 benchmark problems. These calculations point out that the proposed P{sub N}-DG method is capable of producing accurate solutions in small computational time, and that it is able to handle complex 3D geometries. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK Utility Programs for DIF3D

The DIF3D-VARIANT solver is part of the DIF3D code which was developed at ANL for engineering level predictive analysis of nuclear reactor cores. Over the course of the last 20 years, the DIF3D-VARIANT option of DIF3D has become increasingly more used because of improvements in computer memory capabilities and performance. This is beneficial to the users of DIF3D as DIF3D-VARIANT is able to apply both p- and h- spatial mesh refinement strategies in diffusion and transport theory. Given its greater use, utility programs were created to assist in post-processing the DIF3D output. This document discusses the basic execution process of DIF3D and the output files it creates but is primarily focused on being the manual for the utility programs EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK Utility Programs for DIF3D (Rev. 1)

The DIF3D-VARIANT solver is part of the DIF3D code which was developed at ANL for engineering level predictive analysis of nuclear reactor cores. Over the course of the last 20 years, the DIF3D-VARIANT option of DIF3D has become increasingly more used because of improvements in computer memory capabilities and performance. This is beneficial to the users of DIF3D as DIF3D-VARIANT is able to apply both p- and h-spatial mesh refinement strategies in diffusion and transport theory. Given its greater use, utility programs were created to assist in post-processing the DIF3D output. This document discusses the basic execution process of DIF3D and the output files it creates but is primarily focused on being the manual for the utility programs EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Entangling-gate error from coherently displaced motional modes of trapped ions

Entangling gates in trapped-ion quantum computers are most often applied to stationary ions with initial motional distributions that are thermal and close to the ground state, while those demonstrations that involve transport generally use sympathetic cooling to reinitialize the motional state prior to applying a gate. Future systems with more ions, however, will face greater nonthermal excitation due to increased amounts of ion transport and exacerbated by longer operational times and variations over the trap array. In addition, pregate sympathetic cooling may be limited due to time costs and laser access constraints. In this paper, we analyze the impact of such coherent motional excitation on entangling-gate error by performing simulations of Mølmer-Sørenson (MS) gates on a pair of trapped-ion qubits with both thermal and coherent excitation present in a shared motional mode at the start of the gate. Here, we quantify how a small amount of coherent displacement erodes gate performance in the presence of experimental noise, and we demonstrate that adjusting the relative phase between the initial coherent displacement and the displacement induced by the gate or using Walsh modulation can suppress this error. We then use experimental data from transported ions to analyze the impact of coherent displacement on MS-gate error under realistic conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Orbit-modulated transport and sources in time-dependent plasmas

The orbit-modulation framework for steady-state ion-orbit boundary fluxes is reformulated to apply to time-dependent plasmas. The new formulation can quantify the extent to which orbit fluxes are affected by the time dependence of the distribution function, and may be applied both to time-dependent physical scenarios and to computationally expensive simulations that are unable to reach a true steady state. As the system tends to a time-independent state, the orbit-loss fluxes approach their steady-state values, which are determined entirely by upstream transport (turbulent and collisional) and upstream sources and sinks. This result holds for the fluxes of any moments that are expressible in terms of the invariants that determine the orbits, including parallel toroidal angular momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Uncertainty Quantification for Neutron Shield Using Convolutional Neural Networks

Uncertainty quantification from radiation transport calculations was conducted using a Bayesian inference approach. A surrogate model, using a convolutional neural network, was employed to emulate the neutron fluence, which was simulated with a Monte Carlo radiation transport model. This allowed for a computationally cheap approach to evaluate input parameters and to sample their corresponding posterior probability distributions. Experimental data from the literature were employed to perform uncertainty quantification studies for concrete shields. As a result, the method is a nonintrusive approach that enables studies with multiple input parameters and can be applied to any radiation transport model.

Bayesian inference↗
DOE OSTI2023

Multifunctional Chiral Chemically‐Powered Micropropellers for Cargo Transport and Manipulation

Practical applications of synthetic self-propelled nano and microparticles for microrobotics, targeted drug delivery, and manipulation at the nanoscale are rapidly expanding. However, fabrication limitations often hinder progress, resulting in relatively simple shapes and limited functionality. Here, taking advantage of 3D nanoscale printing, chiral micropropellers powered by the hydrogen peroxide reduction reaction are fabricated. Due to their chirality, the propellers exhibit multifunctional behavior controlled by an applied magnetic field: spinning in place (loitering), directed migration in the prescribed direction, capture, and transport of polymer cargo particles. Design parameters of the propellers are optimized by computation modeling based on mesoscale molecular dynamics. It is predicted by computer simulations, and confirmed experimentally, that clockwise rotating propellers attract each other and counterclockwise repel. These results shed light on how chirality and shape optimization enhance the functionality of synthetic autonomous micromachines.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Electric-field-tunable electronic nematic order in twisted double-bilayer graphene

Graphene-based moiré systems have attracted considerable interest in recent years as they display a remarkable variety of correlated phenomena. Besides insulating and superconducting phases in the vicinity of integer fillings of the moiré unit cell, there is growing evidence for electronic nematic order both in twisted bilayer graphene and twisted double-bilayer graphene (tDBG), as signaled by the spontaneous breaking of the threefold rotational symmetry of the moiré superlattices. Here, we combine symmetry-based analysis with a microscopic continuum model to investigate the structure of the nematic phase of tDBG and its experimental manifestations. First, we perform a detailed comparison between the theoretically calculated local density of states and recent scanning tunneling microscopy data (arXiv:2009.11645) to resolve the internal structure of the nematic order parameter in terms of the layer, sublattice, spin, and valley degrees of freedom. We find strong evidence that the dominant contribution to the nematic order parameter comes from states at the moiré scale rather than at the microscopic scale of the individual graphene layers, which demonstrates the key role played by the moiré degrees of freedom and confirms the correlated nature of the nematic phase in tDBG. Secondly, our analysis reveals an unprecedented tunability of the orientation of the nematic director in tDBG by an externally applied electric field, allowing the director to rotate away from high-symmetry crystalline directions. Here, we compute the expected fingerprints of this rotation in both STM and transport experiments, providing feasible ways to probe it. Rooted in the strong sensitivity of the flat bands of tDBG to the displacement field, this effect opens an interesting route to the electrostatic control of electronic nematicity in moiré systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

SISGR: Chemomechanics of Far-From-Equilibrium Interfaces (COFFEI)

Portable, reliable, and deployable devices for energy storage and conversion require fundamental changes in design of solid-state composites comprising ceramics and metals. These materials comprise the electrodes and electrolytes of next-generation solid-oxide fuel cells and all solid-state batteries, forming solid-state functional composites. The advent of solid-state batteries – which replace liquid electrolytes with solid electrolytes capable of lithium ion transport for reliable energy storage in portable batteries – and the increased demand space for all solid-state fuel cells capable of oxygen reduction at intermediate temperatures remain important challenges for improved material stability and decreased system cost. However, little is understood about three fundamental facets of materials that enable such solid-state energy applications. First, how do such materials deform, fracture, or delaminate under operando conditions? Second, how does such mechanical deformation limit or facilitate electronic and ionic transport within and across such material interfaces? Third, how we can predictably design interface-rich composites to engineer both structural and electrochemical stability? This COFFEI Group comprised expertise from Materials Science & Engineering and Nuclear Science & Engineering to integrate unique in situ experiments, simulations, and fabricated interfaces that address these fundamental questions in solid-state interfaces of nanoscale composites that will guide solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric materials that enable such applications. In particular, we built on COFFEI’s understanding of chemomechanical coupling among defect concentrations, ionic transport, electron transport, and stored elastic energy that is particularly acute in the far-from-equilibrium conditions typical of energy device applications. By tailoring our focus to solid-state interfaces, we addressed these important issues by (a) developing and applying advanced in situ and ex situ characterization tools to characterize model materials and interfaces synthesized with molecular-level control, under both laboratory-controlled and extreme environments representative of energy device operation; and (b) employing computational modeling and simulation frameworks to predict transport mechanisms, reactivity and stability of these model materials and interfaces under significant chemical strains typical of energy device operation. Recent progress provided insights to additional materials systems and electrochemomechanical fatigue and fracture that were not fully envisioned when the program was initiated. Specifically, in the final three years of COFFEI we pursued two integrated thrusts, with complementary focus. Thrust I focused on failure-resistant electrochemomechanical composites, while Thrust II focused on strain-modulated conductivity and reactivity across interfaces. In contrast to our initial COFFEI focus, these thrusts concentrated wholly on solid-state material interfacial interactions and included greater integration of multiscale visualization including in situ electron microscopy of strained structures/interactions and mesoscale simulations. Successful development of functionally superior and long-lived battery and fuel cell systems and stress adaptable oxides requires a deeper, fundamental understanding of the coupling among the historically important subfields of solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric metal oxide electrodes. In this program, the understanding and the application of chemomechanical coupling of defect concentrations, ionic transport, electro-catalytic activity and stored elastic energy, particularly acute in the far-from-equilibrium conditions typical of energy device applications, are being refined and implications for device operation clarified, including for miniaturized solid-state batteries and fuel cells.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

risingMicrobubbleLattice

This code is a simulation case to be run with OpenFOAM, an open-source computational fluid dynamics software. A gmsh mesh file is also included. Specifically, this simulation demonstrates the transport of a single microbubble rising through an ordered lattice due to an applied flow. The bubble deforms as it squeezes through the pores of the lattice. Is

Guo, Jack [Lawrence Livermore National Laboratory ↗
DOE OSTI2024

Effect of Cathode Cooling in Three-Dimensional Simulations of an Atmospheric Pressure Glow Discharge

The Atmospheric Pressure Glow Discharge (APGD) is a relatively simple and versatile plasma source used in a wide range of applications. Active cooling of the cathode can effectively mitigate instabilities, leading to glow-to-arc transitions. This study investigates the effect of varying the degree of cathode cooling in APGD with a planar cathode in helium. The plasma flow model incorporates mass conservation, chemical species transport, momentum conservation, conservation of thermal energy of heavy species and of electrons, and electrostatics. The model is applied to time-dependent simulations through a three-dimensional computational domain describing the whole discharge, without geometric symmetry or steady-state assumptions. Simulations of an experimentally characterized APGD explore the effects of electric current and cathode cooling—ranging from thermally insulated to extreme convective cooling. Results show the formation of an annular region with high electric field over the cathode surface under conditions of high current and low cooling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗