SEARCH · Search NASA

Results for “Functional expansion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

iSPECTRON: a simulation interface for linear and nonlinear spectra with ab-initio quantum chemistry software

We introduce iSPECTRON, an open source (under the Educational Community License version 2.0) program that parses data from common quantum chemistry software (NWChem, OpenMolcas, Gaussian, Cobramm, etc.), produces the input files for the simulation of linear and nonlinear spectroscopy of molecules with the Spectron code, and analyzes the spectra with a broad range of tools. Vibronic spectra are expressed in term of the electronic eigenstates, obtained through quantum chemistry computations, and vibrational/bath effects are incorporated in the framework of the displaced harmonic oscillator model, where all required quantities are computed at the Franck-Condon point. The code capabilities are illustrated by simulating linear absorption, transient absorption and two dimensional electronic spectra of the pyrene molecule. Two levels of electronic structure theory, TDDFT (with NWChem) and RASSCF/RASPT2 (with OpenMolcas), are compared where possible. Acknowledgements: F.S., A.N., D.R.N., N.G., S.M, M.G. acknowledge support from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division under Award Nos. DE-SC0019484, KC-030103172684. The Spectron code was developed with support from the National Science Foundation (Grant CHE- 1953045). This research benefited from computational resources provided by EMSL, a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research and located at PNNL. PNNL is operated by Battelle Memorial Institute for the United States Department of Energy under DOE Contract No. DE-AC05-76RL1830.

Segatta, Francesco↗

DOE OSTI2023

Chemical heterogeneity modulated zero thermal expansion alloy over super-wide temperature range

Chemical heterogeneity is usually avoided in solution chemistry, but it may still occur with sometimes dramatic effects on target materials and their properties. Here, we propose chemical heterogeneity as a counterintuitive strategy to design high-performance zero thermal expansion (ZTE) alloys. We apply this approach in a Hf-Ti-Fe alloy with excess Fe in the Hf/Ti sublattice and produce Hf/Ti concentration alternations at the micro level. Such chemical heterogeneity regulates local magnetic interactions in alloy and triggers a dispersed magnetic phase transition that modulates the thermal expansions at the micro level and hence results in a remarkable ZTE behavior over a super-wide temperature window from 10 to 480 K. This mechanism is supported by comprehensive studies on morphological microstructures, crystal and magnetic structures, and theoretical calculations. The strategy of local chemical heterogeneity opens up an avenue to design ZTE and the related functional materials directly via microstructure engineering.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Bayesian-based response expansion and uncertainty quantification using sparse measurement sets

Systems subjected to dynamic loads often require monitoring of their vibrational response, but limitations on the total number and placement of the measurement sensors can hinder the data-collection process. Here, we present an indirect approach to estimate a system’s full-field dynamic response, including all uninstrumented locations, using response measurements from sensors sparsely located on the system. This approach relies on Bayesian inference that utilizes a system model to estimate the full-field response and quantify the uncertainty in these estimates. By casting the estimation problem in the frequency domain, this approach utilizes the modal frequency response functions as a natural, frequency-dependent weighting scheme for the system mode shapes to perform the expansion. This frequency-dependent weighting scheme enables an accurate expansion, even with highly correlated mode shapes that may arise from spatial aliasing due to the limited number of sensors, provided these correlated modes do not have natural frequencies that are closely spaced. Furthermore, the inherent regularization mechanism that arises in this Bayesian-based procedure enables the utilization of the full set of system mode shapes for the expansion, rather than any reduced subset. This approach can produce estimates when considering a single realization of the measured responses, and with some modification, it can also produce estimates for power spectral density matrices measured from many realizations of the responses from statistically stationary random processes. A simply supported beam provides an initial numerical validation, and a cylindrical test article excited by acoustic loads in a reverberation chamber provides experimental validation.

42 ENGINEERING↗

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 OSTI2020

In Situ High-Temperature Synchrotron Diffraction Studies of (Fe,Cr,Al) 3 O 4 Spinels

The modeling of a loss-of-coolant-accident scenario involving nuclear fuels with FeCrAl cladding materials in consideration to replace a Zircaloy requires knowledge of the thermodynamics of oxidized structures. At temperatures higher than 1500 degrees C, oxidation of FeCrAl alloys forms (Fe,Cr,Al) 3 O 4 spinels. In situ high-energy X-ray diffraction in a conical nozzle levitator installed at beamline 6-ID-D of the APS was used to study the structural evolution of the oxides as a function of the temperature. Single-phase (spinel) and multiphase (spinel-corundum-FeAlO 3 ) regions are mapped as a function of the temperature for three different compositions of FeCrAl oxidation products. The thermal expansion coefficients and cation distribution in the spinel structure have been refined. The temperature at which complete melting of the fuel cladding is expected has been determined by the liquidus temperatures of the oxidized products to be between 1657 and 1834 degrees C in a 20% O 2 /Ar atmosphere using the cooling trace method. The liquidus temperature increases with increasing Al and Cr content in the spinel phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Quantifying a partial polyamorphic transition in a cerium-based metallic glass during cooling

Cerium-based metallic glasses are prototype polyamorphous systems with pressure-induced polyamorphic transitions extensively reported. Cooling typically has a similar effect on materials as compression with regard to reducing volume. However, previous studies show dramatically different behavior of Ce-based metallic glasses between cooling and compression, whose origin remains unclear. Here, using in situ low-temperature synchrotron high-energy x-ray diffraction, the structural evolution of a Ce68Al10Cu20Co2 metallic glass is accurately determined and analyzed by a structure factor and a reduced pair distribution function (PDF) during cooling from 298 to 83 K. An unusually large linear thermal expansion coefficient is revealed, which is associated with both continuous but inconsistent structural changes between the two subpeaks of the first atomic shell in terms of average bond lengths and coordination numbers. These phenomena are suggested to be attributed to a gradual 4f electron delocalization of only a minimal amount (∼2.6% at 83 K) of Ce atoms by quantitative analysis of the PDF data. However, a previously expected global polymorphic transition from a low-density amorphous state to a high-density amorphous state with an abrupt volume collapse is not observed. Moreover, electrical resistivity also shows a continuous increase during cooling without any sharp change. It is clarified that cryogenic temperatures could facilitate but are not powerful enough alone to trigger a global polymorphic transition in the Ce68Al10Cu20Co2 metallic glass, suggesting a wide distribution of its local atomic environment.

Chen, Zhi↗

DOE OSTI2022

Permutation invariant polynomial neural network based diabatic ansatz for the ( E + A ) × ( e + a ) Jahn–Teller and Pseudo-Jahn–Teller systems

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of $H^3_+$ is studied, which is an (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. Further, the adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Dynamic asymmetric strain imprinted into substrates by an oxide thin film

In film-substrate systems, the substrate role is often considered to be limited to providing static mechanical constraints. Dynamic film-substrate interactions when a structural change in the film modifies the substrate are generally disregarded. Here, using combined x-ray and electron microscopies, we observed that an electrically induced filament in a vanadium dioxide film created strong asymmetric strain in an underlying sapphire substrate. This asymmetric substrate strain fed back into the film and defined the filament expansion direction, revealing the importance of film-substrate dynamic interactions in determining film functionality. Furthermore, the strain imprint propagated at least tens of micrometers deep into the substrate, exceeding the film thickness by more than 200-fold, potentially enabling substrate functionalization as an active mechanical coupling media in three-dimensional integrated microelectronic architectures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Astrometric Calibration and Performance of the Dark Energy Spectroscopic Instrument Focal Plane

Abstract The Dark Energy Spectroscopic Instrument, consisting of 5020 robotic fiber positioners and associated systems on the Mayall telescope at Kitt Peak, Arizona, is carrying out a survey to measure the spectra of 40 million galaxies and quasars and produce the largest 3D map of the universe to date. The primary science goal is to use baryon acoustic oscillations to measure the expansion history of the universe and the time evolution of dark energy. A key function of the online control system is to position each fiber on a particular target in the focal plane with an accuracy of 11 μ m rms 2D. This paper describes the set of software programs used to perform this function along with the methods used to validate their performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

New H(z) Measurement at Redshift = 0.12 with DESI Data Release 1

The Hubble parameter (H(z)) is a function of the redshift and a reliable measurement is very important to understand the expansion history of the Universe. In this work, we perform full-spectrum fitting using BAGPIPES on more than four thousand massive, passively evolving galaxies released by the DESI collaboration to estimate their cosmological-independent stellar ages and star formation histories, and derive a new measurement of H(z = 0.12) = 71.33 ± 4.20 km s −1 Mpc −1 , which is well consistent with those derived in other ways.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

The MOOSE fluid properties module

The fluid properties module in MOOSE~\cite{lindsay2022moose} serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, the MOOSE thermal hydraulics module, SAM~\cite{sam}, RELAP-7~\cite{relap7}, Sockeye~\cite{sockeye}, Pronghorn-subchannel~\cite{subchannel} and the MOOSE porous flow module~\cite{porous}. These applications are used to solve coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, heat pipe modeling~\cite{sockeye} and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled flow simulations~\cite{anl_sam_pgh,osti_1889653}. Each application has historically driven the implementation of several fluid properties, which were later extended to be compatible with other applications. The unique diversity of applications of the module, due to its presence in MOOSE, has driven its expansion to new fluids, such as advanced nuclear reactor coolants and, more recently, arbitrary functions or tables-based property definitions, as detailed in section~\ref{content}, as well as numerous thermophysical properties and variable sets, as detailed in subsection~\ref{sec:prop}. The need for different discretizations of flow equations based on the compressibility of the fluid has motivated support for both primitive (pressure- and temperature-based) and conservative (internal energy- and specific volume-based) flow variables in the module. Thermodynamic relations are used to convert between these two formulations, as needed. The module serves a dual purpose of providing fluid properties for direct use in flow simulations and facilitating the implementation of user-specific fluid properties. Contributions of new properties for existing fluids or new fluids are strongly encouraged.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Self-consistent optimization of the z expansion for B -meson decays

We discuss the self-consistency imposed by the analyticity of regular parts of form factors, appearing in the z expansion for semileptonic B -meson decays, when fitted in different kinematic regions. Relying on the uniqueness of functions defined by analytic continuation, we propose four metrics which measure the departure from the ideal analytic self-consistency. We illustrate the process using Belle data for B → D ℓ ν ℓ with the two kinematic regions chosen as the five low- z and the five high- z bins. For this specific example, the metrics provide consistent indications that some choices (order of truncation, Boyd-Grinstein-Lebed or Bourrely-Caprini-Lellouch) made in the form of the z expansion can be optimized. However, other choices ( z origin, location of isolated poles and threshold constraints) appear to have very little effect on these metrics. On the other hand, changing the kinetic regions affects the results and should also be considered in the optimization process. We briefly discuss the implication for optimization of the z expansion for nucleon form factors relevant for neutrino oscillation experiments. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Graphically Contracted Function Construction with the Recursive Pairwise Merge Algorithm

A new implementation of a recursive pairwise merge algorithm to construct a GCF from a list of CSF expansion coefficients is presented. The essential new feature is the preallocation of some work arrays used within the intermediate steps of the merge procedure. This results in roughly an order of magnitude improvement in overall efficiency and also approximately eliminates a factor of n, the molecular orbital dimension, from the original implementation. Initial application of this merge procedure to a series of H m molecules shows that the GCF wave functions can be represented well both with delocalized canonical Hartree-Fock orbitals and with localized molecular orbitals. Finally, for a given wave function complexity, as measured by the average facet count, $\bar{\text {f}}$, the delocalized Hartree-Fock orbitals show smaller errors for small $\bar{\text {f}}$ values, while the localized orbitals show smaller errors for larger $\bar{\text {f}}$ values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗