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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fast Approximate Methods for Modified Gravity Cosmological Simulations

The accelerated expansion of the Universe poses a major theoretical puzzle. Although the assumption of a non-zero cosmological constant provides a minimal extension of general relativity that is consistent with observational data, many theories of modified gravity have been suggested as possible alternatives. Predictions of structure formation for these models in the fully nonlinear regime are very expensive and it is difficult, if not impossible, to explore such a huge space of models and parameters using high-resolution N-body simulations. Even in the mildly nonlinear regime, perturbative methods can become extremely complex. We explore whether simplified dynamical approximations, applicable for a certain set of cosmological probes, can be used to investigate models of modified gravity with acceptable accuracy in the latter instance. For the case of chameleon gravity, we find that these methods can indeed be used to explore the region around the baryon acoustic oscillation scale, k ~ 0:1 hMpc-1 but not much further.

cosmology↗
DOE OSTI2024

Sodium-Cooled Fast Reactor Reference Plant Model

This report details the progress of Idaho National Laboratory (INL) in creating a reference plant multiphysics model for the Advanced Burner Test Reactor (ABTR). This model was developed under Task 13 of the U.S. Nuclear Regulatory Commission project “Development and Modeling Support for Advanced Non-Light Water Reactors,” and is an extension of the reference plant model developed in Task 4b, which was improved upon in the following ways. (1) The discrete ordinates method was used in lieu of the super-homogenization (SPH)-corrected diffusion approximation in order to better capture the anisotropic scattering contribution and the neutron leakage change due to thermal expansion. (2) The novel neutronic spatial discretization approach, termed the ring-heterogeneous (RH) approximation, was conceptualized and introduced to capture the differential expansion of the materials in the core. This new technique proved capable of preserving fission rates and maintaining the eigenvalue within 2.5% and 266 pcm with 9 neutron energy groups, respectively. Separating the different materials in the core enables the differential expansion of materials to be explicitly accounted for, eliminating the need for problem-specific cross-section functionalization techniques. (3) The SAM model for the core and system thermal-hydraulics analysis was updated to include 61 channels instead of just four representative ones. This enables users to obtain improved spatial resolution for sodium temperature and density scalar fields. (4) All the mesh files were created via the Multiphysics Object-Oriented Simulation Environment (MOOSE) Reactor module, eliminating all reliance on external tools for mesh creation. (5) Finally, the fuel axial expansion now leverages the HT9 and UPuZr material properties that have been validated against experimental data. The reference plant model was used to perform a full-core unprotected loss of flow (ULOF) transient calculation, including neutronics, thermal and mechanical feedback mechanisms. Future work will be devoted to further enhancements of the model. Potential improvements to the model include the addition of the control rod driveline expansion feedback and the upgrading of the support plate model so as to explicitly include 3D effects. Additionally, a Nuclear Energy Advanced Modeling and Simulation funded parallel effort has completely automated the creation of the ring-heterogeneous (RH) mesh from the fully heterogeneous (FH) geometry, thus maximizing user friendliness for the sodium fast reactor sodium-cooled fast reactor (SFR) workflow and will be incorporated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

First-principles thermodynamics of Al 10 ⁢V: An analytical treatment of localized anharmonic modes

Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 ⁢V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10⁢ V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 ⁢V phase at elevated temperatures.

anharmonic lattice dynamics↗
DOE OSTI2022

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Field moment expansion method for interacting bosonic systems

In this work, we introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗
DOE OSTI2020

Off-shell effective energy theory: A unified treatment of the Hubbard model from $\textit{d} = 1$ to $\textit{d = ∞}$

Here we propose an exact formalism, off-shell effective energy theory (OET), which provides a thermodynamic description of a generic quantum Hamiltonian. The OET is based on a partitioning of the Hamiltonian and a corresponding density matrix ansatz constructed from an off-shell extension of the equilibrium density matrix; and there are dual realizations based on a given partitioning. To approximate OET, we introduce the central point expansion (CPE), which is an expansion of the density matrix ansatz, and we renormalize the CPE using a standard expansion of the ground-state energy. In this work we showcase the OET for the one-band Hubbard model in $\textit{d}$ = 1 , 2, and $\infty$, using a partitioning between kinetic and potential energy, yielding two realizations denoted as $\mathcal{K}$ and $\mathcal{X}$. OET shows favorable agreement with exact or state-of-the-art results over all parameter space, and has a negligible computational cost. Physically, $\mathcal{K}$ describes the Fermi liquid, while $\mathcal{X}$ gives an analogous description of both the Luttinger liquid and the Mott insulator. Our approach should find broad applicability in lattice model Hamiltonians, in addition to real materials systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Systematic Extraction of Pion Parton Distributions Using Threshold Resummation and Applications to Pion TMDs

Following our recent Monte Carlo determination of the pion?s PDFs from Drell-Yan (DY) and leading neutron electroproduction data from the Jefferson Lab Angular Momentum (JAM) collaboration, we extend the analysis by including effects from threshold resummation. At higher orders in the strong coupling, ?S, soft gluon emissions cause large logarithmic corrections, which become important in the qq¯ channel of the DY partonic cross section near threshold. These corrections can be summed over all orders of ?S. However, different prescriptions exist for how the threshold resummation is implemented, for instance, using varying levels of approximation in the Minimal Prescription with cosine, expansion, and double Mellin methods. We present the Monte Carlo results of the first simultaneous fit of the valence, sea, and gluon distributions in the pion taking into account the ambiguities in the resummation calculations. The wide ranges of valence distributions at large x and the effective behavior of the valence distribution as x approaches 1 is discussed. While the PDFs are extracted through collinear factorization, we additionally present a dedicated study to the low and large lepton pair transverse momentum distributions (TMDs) in the same DY experimental data. By adjusting the nonperturbative TMD components, we analyze the degree of compatibility between data and theory across all regions of transverse momentum dependence. We attempt to match the description of the low transverse momentum region using Collins-Soper-Sterman (CSS) formalism and the large transverse momentum region, which is best described using collinear factorization.

Sato, Nobuo↗
DOE OSTI2024

The Thermophysical Properties of TcO2

Technetium-99 is a highly radioactive isotope with a long half-life that is common in nuclear waste. It volatizes at a low temperature, which poses a significant challenge to the clean-up and containment processes. Due to difficulties in purifying technetium compounds, their thermophysical properties have not been measured or calculated. Here, first principle methods are used along with the quasi quasi-harmonic harmonic approximation to compute the Debye temperature, volumetric thermal expansion coefficient, bulk modulus, and heat capacity of rutile TcO2 for temperatures ranging from 0 to 1500 K and applied pressures ranging from 0 to 255 GPa. The computed atomic structures agree well with the results from diffraction measurements. The computed thermophysical properties are in the neighborhood of other rutile metal oxides and, in particular, are within approximately 10–13% of rutile ReO2, which is frequently used as a substitute for TcO2 in experimental studies.

Zhong, Hong↗
DOE OSTI2027

From zonal to nodal capacity expansion planning: Spatial aggregation impacts on a realistic test-case

Solving power system capacity expansion planning (CEP) problems at realistic spatial resolutions is computationally challenging. Thus, a common practice is to solve CEP over zonal models with low spatial resolution rather than over full-scale nodal power networks. Due to improvements in solving large-scale stochastic mixed integer programs, these computational limitations are becoming less relevant, and the assumption that zonal models are realistic and useful approximations of nodal CEP is worth revisiting. Here, this work is the first to conduct a systematic computational study on the assumption that spatial aggregation can reasonably be used for ISO-scale CEP. By considering a realistic, large-scale test network based on the state of California with over 8000 buses, we find that well-designed small spatial aggregations can yield good approximations but that coarser zonal models may result in large distortions of investment decisions, e.g., capacity under-investment of up to 41% for the lowest resolution model considered.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Classical eikonal from Magnus expansion

In a classical scattering problem, the classical eikonal is defined as the generator of the canonical transformation that maps in-states to out-states. It can be regarded as the classical limit of the log of the quantum S-matrix. In a classical analog of the Born approximation in quantum mechanics, the classical eikonal admits an expansion in oriented tree graphs, where oriented edges denote retarded/advanced worldline propagators. The Magnus expansion, which takes the log of a time-ordered exponential integral, offers an efficient method to compute the coefficients of the tree graphs to all orders. We exploit a Hopf algebra structure behind the Magnus expansion to develop a fast algorithm which can compute the tree coefficients up to the 12th order (over half a million trees) in less than an hour. In a relativistic setting, our methods can be applied to the post-Minkowskian (PM) expansion for gravitational binaries in the worldline formalism. We demonstrate the methods by computing the 3PM eikonal and find agreement with previous results based on amplitude methods. Importantly, the Magnus expansion yields a finite eikonal, while the naïve eikonal based on the time-symmetric propagator is infrared-divergent from 3PM on.

Black Holes↗
DOE OSTI2021

Climate Change Mitigation Potential of Wind Energy

Global wind resources greatly exceed current electricity demand and the levelized cost of energy from wind turbines has shown precipitous declines. Accordingly, the installed capacity of wind turbines grew at an annualized rate of about 14% during the last two decades and wind turbines now provide ~6–7% of the global electricity supply. This renewable electricity generation source is thus already playing a role in reducing greenhouse gas emissions from the energy sector. Here we document trends within the industry, examine projections of future installed capacity increases and compute the associated climate change mitigation potential at the global and regional levels. Key countries (the USA, UK and China) and regions (e.g., EU27) have developed ambitious plans to expand wind energy penetration as core aspects of their net-zero emissions strategies. The projected climate change mitigation from wind energy by 2100 ranges from 0.3–0.8 °C depending on the precise socio-economic pathway and wind energy expansion scenario followed. The rapid expansion of annual increments to wind energy installed capacity by approximately two times current rates can greatly delay the passing of the 2 °C warming threshold relative to pre-industrial levels. To achieve the required expansion of this cost-effective, low-carbon energy source, there is a need for electrification of the energy system and for expansion of manufacturing and installation capacity.

17 WIND ENERGY↗
DOE OSTI2022

A fast approximation to supershell partition functions

A formula for supershell partition functions, which play a major role in the Super Transition Array approach to radiative-opacity calculations, is derived as a functional of the distribution of energies within the supershell. It consists in an alternative expansion for an arbitrary number of electrons or holes which also allows for quick approximate evaluations with truncated number of terms in the expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Jets from shocked metal surfaces with grooves: Missing experiments

Many studies have investigated the mass outflows generated when a planar shock transits an imperfect (“defected”) metal surface, where the defects are symmetric triangular or sinusoidal grooves. Yet a fundamental question remains unanswered: how does the quantity of outflow mass and its maximum velocity vary as a function of the groove cross-sectional aspect ratio? We identify two sets of missing experiments that must be addressed to answer the question. The aspect ratio (groove depth over width) is equivalently represented by θ, the cross-sectional half angle, or by η 0 k, the amplitude multiplied by an effective wavenumber. Low θ (high η 0 k) grooves comprise the first set of missing experiments, which are necessary to determine the validity of theoretical predictions of the nonlinear regime (η 0 k ≥ 1, θ < 57.5°). The second set of missing experiments are those in which the volume of the groove (or equivalently, the axial cross-sectional area) has been held constant as θ or η 0 k are varied. Such experiments are necessary to independently measure the effects of variations in groove volume and groove aspect ratio on the resulting jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗
DOE OSTI2024

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Cosmic Growth Rate as Probe of Dark Matter and Gravity. Final Report

Some 95% of the energy content of the universe lies outside the standard model of particle physics, with approximately 70% arising from physics that accelerates cosmic expansion in opposition to normal gravity and about 25% dark matter that clusters under gravity but does not interact directly with light. This project made significant advances in developing methods to extract information on the former – dark energy or modified gravity – and latter – dark matter through their effects on the expansion history of the universe and the cosmic growth of large scale structures under the influence of gravity and dark energy. These results directly benefit experimental/observational cosmic surveys of large scale structure imaging hundreds of millions of galaxies and mapping spectroscopically the 3D positions of tens of millions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Multi-axis Accelerometry and Rotation Sensing using a Point Source Atom Interferometer

A point source atom interferometer (PSI) is a device where atoms are split and recombined by applying a temporal sequence of Raman pulses during the expansion of a cloud of cold atoms behaving approximately as a point source. Unlike a conventional light pulse atom interferometer, the PSI can produce a signal that corresponds to multi-axis rotation only, independent of acceleration. In addition, it can be used to measure acceleration along one direction, independent of rotation. Here, we describe a modified PSI that can be used to measure multi-axis rotation and multi-axis acceleration. Specifically, this type of PSI can be used to measure two-axes rotation around the directions perpendicular to the light pulses, as well as the acceleration in all three directions, with only one pair of Raman beams. Using two pairs of Raman beams in orthogonal directions sequentially, such a scheme would enable the realization of a complete atom interferometric inertial measurement unit.

Li, Jinyang↗