Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximations and Expansions”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

NO PLIF flow visualization and time-resolved temperature distributions in laser induced breakdown plumes

NO planar laser induced fluorescence (PLIF) is used to obtain images of laser-induced breakdown plasma plumes in NO-seeded nitrogen and dry air at near atmospheric pressure. Single-shot PLIF-images show that the plume development 5–50 μs after the breakdown pulse is fairly reproducible shot-to-shot, although the plume becomes increasingly stochastic on longer timescales, 100–500 μs. The stochastic behavior of the plume is quantified using probability distributions of the loci of the plume boundary. Analysis of the single-shot images indicates that the mixing of the plume with ambient gas on sub-ms time scale is insignificant. The induced flow velocity in the plume is fairly low, up to 30 m s –1 , suggesting that laser breakdowns are ineffective for mixing enhancement in high speed flows. The ensemble-averaged PLIF images indicate the evolution of the plume from an initially elongated shape to near-spherical to toroidal shape, with a subsequent radial expansion and formation of an axial jet in the center. Temperature distributions in the plume in air are obtained from the NO PLIF images, using two rotational transitions in the NO(X, v' = 0 → A, v'' = 0) band, J'' = 6.5 and 12.5 of the Q R 12 + Q 2 branch. The results indicate that the temperature in the plume remains high, above 1000 K, for approximately 100 μs, after which it decays gradually, to below 500 K at 500 μs. The residual NO fraction in the plume is ~0.1%, indicating that repetitive laser-assisted ignition may result in significant NO-generation. Furthermore, these measured temperature and velocity distributions can be used for detailed validation of kinetic models of laser-induced breakdown and assessment of their predictive capability.

42 ENGINEERING↗
DOE OSTI2024

Generalized models for inflationary preheating: Oscillations and symmetries

The paradigm of the inflationary universe provides a possible explanation for several observed cosmological properties. In order for such solutions to be successful, the universe must convert the energy stored in the inflaton potential into standard model particles through a process known as reheating. In this paper, we reconsider the reheating process for the case where the inflaton potential respects an approximate (but spontaneously broken) conformal symmetry during the reheating epoch. After reviewing the Effective Field Theory of Reheating, we present solutions for the nonlinear oscillations of the inflaton field, derive the corresponding Hill’s equation for the coupled reheating field, and determine the stability diagram for parametric resonance. For this class of models —the simplest realization being a scalar field with a quartic term—the expansion of the universe drives the coupled field toward a more unstable part of parameter space, in contrast to the standard case. We also generalize this class of models to include quadratic breaking terms in the potential during the reheating epoch and address the process of stability in that universality class of models.

Barrowes, Leia [University of Michigan, Ann Arbor]↗
DOE OSTI2021

Ground Vibrational State SiO Emission in the VLA BAaDE Survey

Using a subsample of the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey of stellar SiO masers, we explore the prevalence and characteristics of {sup 28}SiO J = 1 − 0, v = 0 emission. We identify 90 detections of maser, thermal, or composite {sup 28}SiO J = 1 − 0, v = 0 emission out of approximately 13,000 candidate spectra from the National Science Foundation’s Karl G. Jansky Very Large Array (VLA). We find that the detected sources are likely asymptotic giant branch stars belonging to a bright, foreground Milky Way stellar disk population. For the 32 sources showing thermal components, we extract values for outflow velocity by fitting thermal line profiles. We find a range of circumstellar envelope expansion velocities, and compare to previously recorded OH and CO expansion velocities. This preliminary survey is already the largest study of stellar ground vibrational state SiO masers to date, and will be expanded to include the entire VLA BAaDE data set when data reduction for the 18,988 target sources is completed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

What do physics-informed DeepONets learn? Understanding and improving training for scientific computing applications

Physics-informed deep operator networks (DeepONets) have emerged as a promising approach toward numerically approximating the solution of partial differential equations (PDEs). In this work, we aim to develop further understanding of what is being learned by physics-informed DeepONets by assessing the universality of the extracted basis functions and demonstrating their potential toward model reduction with spectral methods. Results provide clarity about measuring the performance of a physics-informed DeepONet through the decays of singular values and expansion coefficients. In addition, we propose a transfer learning approach for improving training for physics-informed DeepONets between parameters of the same PDE as well as across different, but related, PDEs where these models struggle to train well. This approach results in significant error reduction and learned basis functions that are more effective in representing the solution of a PDE.

Deep operator networks↗
DOE OSTI2021

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. In this study, we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗
DOE OSTI2023

Mitigating the noise of DESI mocks using analytic control variates

In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological N-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Light Water Reactor Sustainability Program: Technical and Economic Considerations for Uprate of Existing Nuclear Reactors with Cogeneration

The United States nuclear reactor fleet consists of 63 pressurized water reactors and 31 boiling water reactors and is a pivotal component in the nation's energy infrastructure, supplying approximately 97 GW e of clean power. With the country's commitment to decarbonization by 2050, these reactors are not only instrumental in decarbonizing the electricity grid but also play a critical role in decarbonizing industrial processes, producing clean fuels, and scaling up CO 2 removal. This report delves into the potential for power uprates in the existing fleet to contribute to these decarbonization efforts, focusing on the expansion of capacity for applications such as hydrogen production and carbon capture and sequestration. Building on previous research, the report explores regional market demands for hydrogen, oxygen, and carbon dioxide, financial implications of oxygen and CO 2 sales from high-temperature steam electrolysis systems, and the potential for direct air capture systems paired with uprated nuclear plants.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Electronic scattering off a magnetic hopfion

We study scattering of itinerant electrons off a magnetic hopfion in a three-dimensional metallic magnet described by a magnetization vector S(r). A hopfion is a confined topological soliton of S(r) characterized by an emergent magnetic field Bγ(r)≡ε αβγ S∙(∇ α S×∇ β S)/4≠0 with vanishing average value [B(r)]=0. We evaluate the scattering amplitude in the opposite limits of large and small hopfion radius R using the eikonal and Born approximations, respectively. In both limits, we find that the scattering cross section contains a skew-scattering component giving rise to the Hall effect within a hopfion plane. That conclusion contests the popular notion that the Hall effect in noncollinear magnetic structures necessarily implies [B(r)]≠0. In the limit of small hopfion radius pR<<1, we expand the Born series in powers of momentum p and identify different expansion terms corresponding to the hopfion anisotropy, toroidal moment, and skew scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An Infinite Domain 3D Poisson Solver Based on the Barnes-Hut Algorithm

We present a domain decomposition method for the solution of the 3D Poisson equation with infinite domain boundary conditions. The method is based on an application of the Barnes-Hut tree particle scheme adapted to gridded data. Long range interactions are computed using the first two terms in the Cartesian multipole expansion of Green’s function convoluted with the charge while short range computations are performed using Hockney’s domain doubling algorithm. A standard domain decomposition strategy requires O(N 2 ) applications of Hockney’s algorithm, where N is the number of subdomains that intersect that charge support, while in the present approach only O(Nlog 2 N) such computations suffice. The discretization scheme employed is a sixth order Mehrstellen approximation of the 3D Laplace opera tor. The method exhibits satisfactory accuracy at a substantially reduced computational cost compared to the full domain decomposition Hockney’s algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Liquid Structure of Iron and Iron–Nitrogen–Carbon Alloys Within the Cores of Small Terrestrial Bodies

Abstract Nitrogen has been proposed to be stored within planetary cores, but its effects on the structure and density of molten Fe–alloys have not been explored experimentally. Using energy‐dispersive X‐ray diffraction, we determined the structure of Fe–N(–C) liquids at core conditions (1–7 GPa and 1700–1900°C) within a Paris‐Edinburgh press. Variation of N up to 7 wt.% and C up to 1.5 wt.% results in near‐linear changes in Fe–Fe atom distances and structure factor with increasing light element content. We did not observe a significant pressure‐driven structural transition in Fe–N(–C) liquids. We model the expansion of the Fe–Fe bonds using a modified Birch‐Murnaghan equation of state. With this model, we demonstrate that N or C contamination could lead to an overestimation of the Fe–Fe distances of pure Fe. We observe that the incorporation of 1 wt.% N or C into Fe results in a change in Fe–Fe distances that is twice as significant as the effect of 1 GPa. By approximating the change in volume, we infer that N and C incorporated in liquid iron could contribute to the density deficit observed in the cores of terrestrial bodies.

Pease, Allison [Department of Earth and Environmen↗
DOE OSTI2023

An Assessment of the Laminar Hypersonic Double-Cone Experiments in the LENS-XX Tunnel

This is an investigation on two experimental datasets of laminar hypersonic flows, over a double-cone geometry, acquired in Calspan—University at Buffalo Research Center’s Large Energy National Shock (LENS)-XX expansion tunnel. These datasets have yet to be modeled accurately. A previous paper suggested that this could partly be due to mis-specified inlet conditions. The authors of this paper solved a Bayesian inverse problem to infer the inlet conditions of the LENS-XX test section and found that in one case they lay outside the uncertainty bounds specified in the experimental dataset. However, the inference was performed using approximate surrogate models. Here in this paper, the experimental datasets are revisited and inversions for the tunnel test-section inlet conditions are performed with a Navier–Stokes simulator. The inversion is deterministic and can provide uncertainty bounds on the inlet conditions under a Gaussian assumption. It was found that deterministic inversion yields inlet conditions that do not agree with what was stated in the experiments. An a posteriori method is also presented to check the validity of the Gaussian assumption for the posterior distribution. This paper contributes to ongoing work on the assessment of datasets from challenging experiments conducted in extreme environments, where the experimental apparatus is pushed to the margins of its design and performance envelopes.

42 ENGINEERING↗
DOE OSTI2026

Multi-Objective design of interlocking metasurfaces using conditional diffusion models

Unit cell design remains a major challenge for interlocking metasurfaces, a promising joining technology for dissimilar materials, due to the complex, competing, multivariate design space and the need for rapid adaptation to varying performance requirements. This study explores Conditional Diffusion Models as a design optimization tool for interlocking metasurfaces. Given the complex, competing, multivariate design space for interlocking metasurfaces, unit cell design remains a major challenge for this joining technology. We trained a conditional diffusion model on 25,000 finite element analysis-simulated interlocking metasurface unit cells to generate designs with tailored thermo-mechanical properties (tensile strength, shear strength, and thermal conductivity) based on specified performance criteria. The model demonstrated a success rate of approximately 72 % in producing designs that met specified property bounds. The conditional diffusion model generated both thermally resistive and conductive designs, revealing clear trends in design characteristics: taller, dendritic structures were advantageous for tensile loads, while shorter, robust designs excelled in shear applications. Our findings indicate that the model's performance is more influenced by the breadth of the design space than by the quantity of training data, highlighting the importance of expansive design domains for generating innovative solutions. This work establishes conditional diffusion models as a highly efficient and adaptable tool for rapid interlocking metasurface unit cell design, paving the way for advancements in multi-material joining technologies, as well as highlighting the justification to leverage conditional diffusion models as design tools across complex design domains.

Conditional diffusion models↗
DOE OSTI2025

Design of an alumina forming coating for Nb-base refractory alloys

Refractory multi-principal element alloys (RMPEAs) promise to significantly enhance gas turbine engine efficiency, but their poor oxidation performance inhibits their implementation. Alumina-forming bond coat alloys can provide oxidation protection, but discovering suitable chemistries remains a challenge. We employed a design methodology that screens for alumina-formation capability using Al activity and phase constitution predictions from CalPhaD (Thermo-Calc). Alloy down-selection from approximately 7,800 alloys in the Nb-Si-Ti-Al-Hf system was conducted via analysis of calculated thermodynamic properties with number-density topology style maps. This approach is validated by creating and testing the composition Nb 12 Si 23 Ti 24 Al 36 Hf 5 , which forms protective alumina scales up to 1400 °C and resists pesting at 800 °C. Further, the alloy has an average coefficient of thermal expansion of ~10.1 ppm/K, making it well matched to Nb-based refractory alloys. The methodology will be useful for the design of coatings for RMPEAs, enabling their implementation and significant efficiency benefits sooner.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Moments of the Boltzmann collision operator for Coulomb interactions

Exact moments of the Boltzmann collision operator are calculated in the irreducible Hermitian moment expansion written in terms of the random-velocity variable of each species. The formulas are presented in closed, algebraic form and can be straightforwardly implemented in computer algebra systems. They are valid for two arbitrary masses, temperatures, and flow velocities, and hence include all other existing results derived for distribution functions expanded with respect to reference states of one temperature and flow velocity. In comparison, the Landau collisional moments are good approximations for large Coulomb logarithm and small relative flow velocity, but they fail to predict the correct behavior of most collisional moments for large relative flow even for weakly coupled plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quarkonium Polarization Kinetic Equation from Open Quantum Systems and Effective Field Theories

Recent measurements of polarization phenomena in relativistic heavy ion collisions have aroused a great interest in understanding dynamical spin evolution of the QCD matter. In particular, the spin alignment signature of J/ψ has been recently observed in Pb-Pb collisions at LHC, which may infer nontrivial spin transport of quarkonia in quark gluon plasmas. Motivated by this, we study the spin-dependent in-medium dynamics of quarkonia by using the potential nonrelativistic QCD (pNRQCD) and the open quantum system framework. By applying the Markovian approximation and Wigner transformation, we systematically derive the Boltzmann transport equation for vector quarkonia with polarization dependence in the quantum optical limit. As opposed to the previous study for the spin-independent case where the collision terms depend on chromoelectric correlators, the new kinetic equation incorporates gauge invariant correlators of chromomagnetic fields that determine the recombination and dissociation terms with polarization dependence at the order we are working in the multipole expansion. In the quantum Brownian motion limit, the Lindblad equation with new transport coefficients defined in terms of the chromomagnetic field correlators have also been derived. Our formalism is generic and valid for both weakly-coupled and strongly-coupled quark gluon plasmas. It may be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗