Search NASA⌕ Search

SEARCH · Search NASA

Results for “analyzable function”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Finite-temperature kinetic ferromagnetism in the square-lattice Hubbard model

While the exact phase diagram of the Fermi-Hubbard model remains poorly understood despite decades of progress, nearly 60 years ago Nagaoka proved that a single dopant in an otherwise halffilled Hubbard system can bring about ferromagnetism through kinetic means. The phenomenon was recently observed with ultracold atoms in triangular optical lattices. Here, we explore the kinetic ferromagnetism within the square lattice Hubbard model and its strong-coupling counterpart, the t − J model, at finite temperatures in the thermodynamic limit via numerical linked-cluster expansions. We find evidence of ferromagnetic Nagaoka polarons at dopings up to ∼ 30% away from half filling for a variety of interaction strengths and at temperatures as low as 0.2 of the hopping energy. We map out the boundaries of this phase through analyzing various correlation functions.

Newby, Robin C. [San José State University] (ORCID↗
DOE OSTI2025

Diffusion power spectra as a window into dynamic materials architecture

Chemical recycling of commodity and specialty polymers presents a multifaceted challenge for industrial societies. On one hand, macromolecular architectures must be engineered to yield durable products that, on the other hand, rapidly deconstruct to recyclable monomers under pre-determined conditions. Polymer deconstruction is a chemical process that requires deep understanding of molecular reactivity in heterogeneous media, where porous material architectures evolve in both space and time. To build this understanding, we develop herein experimental and analytical methods describing sets of diffusive eigenmodes that exist within time-varying, non-Euclidean boundary conditions, a situation commonly encountered in the reactive deconstruction of polymers where chain fragments splay, alter their local dynamics, and evolve in their confinement of reacting media. Diffusion power spectra, discerned experimentally by NMR, yield polymer and solvent frequency-domain velocity autocorrelation functions that are analyzed in the context of physical models for chemical reactions parameterized with fractal mathematics. The results connect local motion in polymers to chemical reactivity during acidolysis of circular elastomers.

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

An Improved Double-Gaussian Closure for the Subgrid Vertical Velocity Probability Distribution Function

Abstract The vertical velocity probability distribution function (PDF) is analyzed throughout the depth of the lower atmosphere, including the subcloud and cloud layers, in four large-eddy simulation (LES) cases of shallow cumulus and stratocumulus. Double-Gaussian PDF closures are examined to test their ability to represent a wide range of turbulence statistics, from stratocumulus cloud layers characterized by Gaussian turbulence to shallow cumulus cloud layers displaying strongly non-Gaussian turbulence statistics. While the majority of the model closures are found to perform well in the former case, the latter presents a considerable challenge. A new model closure is suggested that accounts for high skewness and kurtosis seen in shallow cumulus cloud layers. The well-established parabolic relationship between skewness and kurtosis is examined, with results in agreement with previous studies for the subcloud layer. In cumulus cloud layers, however, a modified relationship is necessary to improve performance. The new closure significantly improves the estimation of the vertical velocity PDF for shallow cumulus cloud layers, in addition to performing well for stratocumulus. In particular, the long updraft tail representing the bulk of cloudy points is much better represented and higher-order moments diagnosed from the PDF are also greatly improved. However, some deficiencies remain owing to fundamental limitations of representing highly non-Gaussian turbulence statistics with a double-Gaussian PDF.

Fitch, A. C.↗
DOE OSTI2026

Topsoil bulk geochemical compositions - An updated harmonized global dataset

Mineral weathering is a key biogeochemical process because of the capacity of minerals to stabilize organic matter. However, predicting soil weathering status across large spatial areas still isn’t possible due to a lack of global data and theoretical frameworks. To address this knowledge gap, multiple global datasets of bulk topsoil geochemical compositions have been harmonized using R. These datasets document topsoil bulk geochemical compositions across five continents (n = ~16,000 observations). Source data for these observations include the EuroGEOSurveys Geochemical Baseline Database (FOREGS), the US Geological Survey National Geochemical Database (NASGLP), the Geochemical Atlas of Australia (GAA), the US Geological Survey Alaska Geochemical Database (AGD84), the National Cooperative Soil Survey (NCSS), the European Geochemical Mapping of Agricultural Soil (GEMAS), Ecorespira-Amazon (ERA), the New Zealand Geochemical Baseline Survey (NZ_GBS), and the African Soil Information Service (AFSIS). Major elements observed include Aluminum (Al), Calcium (Ca), Iron (Fe), Potassium (K), Magnesium (Mg), Sodium (Na), Titanium (Ti), Manganese (Mn), Phosphorus (P), Carbon (C), and Sulfur (S). This data package includes the harmonized dataset itself, and the R scripts necessary to harmonize these datasets, in addition to metadata that describes all columns, files, and databases used in this project. Methods & Sampling Step 1 – Databases of geochemical data identified This study aimed to leverage existing measurements of topsoil geochemical data. Databases were first identified and deemed appropriate for inclusion if they were measuring soils and performed these measurements on the <2mm soil fraction. Databases such as NCSS and AGD84 needed more post processing to include in the database and this was done using the NCSS_datamerge_031626 R file and Alaska_USGSmerge_031626 R file, respectively. Step 2 – Database harmonization Once appropriate databases were identified, they were harmonized for ease of analysis using the R script Database_Harmonization_031826. This included removing columns from original datasets that would not be used in analysis (removed columns are noted in the code). Then, data cleaning procedures specific to each dataset were undertaken. This includes standardizing columns to include units and adding metadata columns regarding procedures for analyzing specific elements. Functions for standardizing measurements and units are outline in R files: calculate element_mg_kg_031626, calculate_oxide_wt_perc_031626, change_oxide_caps_031626, and conv_2_numeric_031626. This also included adding a unique identifier for each sample to identify it with its respective database (see CD_ID in data dictionary). Geographic information: Data reflect a compilation of datasets collected globally. Geographic areas covered by each of the datasets include: - EuroGEOSurveys Geochemical Baseline Database (FOREGS) - European continent - North American Soil Geochemical Landscapes (NASGLP) - continental United States and limited parts of Canada (see database key for more details) - National Geochemical Survey of Australia (GAA) - Australia - Alaska geochemical database (AGDB4) - Alaska - National Cooperative Soil Survey (NCSS) - Global measurements, but concentrated in the continental United States - Geochemical data for arable land and land under permanent grass cover in continental Europe (GEMAS) - continental Europe - Ecorespira-Amazon (ERA) - Geochemical data from the Amazon basin - Geochemical baseline data for New Zealand (NZGBS) - New Zealand - Geochemical data collected across continental Africa (AfSIS) - Measurements across Africa

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2024

Uncertainty Quantification in GADRAS Inverse Modeling

The Gamma Detector Response and Analysis Software (GADRAS) package includes an inverse modeling tool that is helpful in identifying characteristics of unknown radioactive materials. Traditionally, uncertainties in this analysis were derived solely from measurement data quality and the fit of synthetic spectra. This paper aims to rigorously quantify additional sources of uncertainty, focusing on uncertainties arising from measurements being analyzed, Detector Response Function (DRF) characterization, and DRF extrapolation. Applying these findings to the BeRPBall benchmark data set, we demonstrated the impact of these uncertainties on plutonium and polyethylene estimates. The results underscore the importance of incorporating diverse uncertainty sources to enhance the accuracy and reliability of GADRAS’s inverse modeling capabilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Microstructure-Based Modeling of Inner Oxygen Pressure in Solid Oxide Electrolysis Cells

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here, a microstructure-based electrochemical model which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

Xue, Fei↗
DOE OSTI2019

Production of Hollow Alloy Blanks with Shape Memory Effect

Abstract—The manufacture of blanks from alloys based on titanium nickelide with shape memory is proposed by means of reverse extrusion with subsequent thermomechanical treatment and crimping. The alloy plasticity is analyzed as a function of the temperature and strain rate. The shaping of such blanks is investigated.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Search for CPT and Lorentz invariance violation in the muon $g -2$ experiment at Fermilab

The Muon g-2 experiment at Fermilab (E989) aims to measure the anomalous magnetic moment of the muon, $a_{\mu}= (g-2)/2$, to a groundbreaking precision of $140$ ppb, obtaining a near four-fold increase in precision over the previous experiment, E821, at the Brookhaven National Laboratory (BNL). The value of $a_{\mu}$ from BNL currently differs from the Standard Model prediction by $\sim 3.7$ standard deviations, suggesting the potential for new physics and therefore, motivating a new experiment.Because the theory predicts this number with high precision, testing the g-factor through experiment provides a stringent test of the SM and can suggest physics beyond the Standard Model. The goal of the Fermilab Muon $g-2$ experiment is to increase the statistical precision by more than a factor of 20 and reduce systematic errors by a factor of 3. By measuring muon precession rate ($\omega_a$) in an external magnetic field, the anomalous magnetic moment will be calculated. This is an incredibly challenging experiment with a unique opportunity to provide new insight into nature. \\ The $g-2$ data also provides a great opportunity for setting the most stringent limits on some of the Standard Model Extension CPT Lorentz violating (LV) parameters in the muon sector. One of the CPT and Lorentz violating signatures that we can look for using $g-2$ data is a sidereal variation of $\omega_a(t)$. Extensive simulation studies confirm that the sensitivity regarding the sidereal varation roughly scales with $\omega_a$ uncertainty. Hence, the $g-2$ experiment at FNAL should be able to reach limits of $\sim 5\times10^{-25}$ GeV. Because the CPT and LV analyses are essentially studies of variations in $\omega_a$ as a function of time and charge, performing an $\omega_a$ analysis sets the stage for the CPT and LV measurement. This dissertation focuses on the methodology of a fully functioning framework and analyzing the Fermilab Muon $g - 2$ Run 2 data containing $\sim 11$ billion events above an energy threshold of $1.7$~GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Exploring the exact limits of the real-time equation-of-motion coupled cluster cumulant Green’s functions

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green’s function approach [Rehr et al., J. Chem. Phys. 152, 174113 (2020)]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz, ensuring the exactness in the expansion limit. In addition, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green’s function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green’s function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green’s function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green’s function approach and other Green’s function approaches in the numerical exact limit.

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

Use of Cardiac Imaging to Evaluate Cardiac Function and Pulmonary Hemodynamics in Patients with Heart Failure

Purpose of Review: Noninvasive hemodynamic assessments in patients with heart failure (HF) are essential for appropriate diagnosis and establishment of the best treatment strategies. Recently, the impact of pulmonary circulation and right ventricular function on prognosis in HF patients has drawn increasing attention. In this article, we explore the usefulness of cardiac imaging for hemodynamic assessments, mainly focusing on echocardiographic evaluation. Recent Findings: The reliability of Doppler echocardiography as a noninvasive alternative to Swan-Ganz catheterization has been well investigated with higher than 80% accuracy for estimating pulmonary artery pressure. Strain measurement and three-dimensional echocardiography are useful for evaluating right ventricular function together with pulmonary circulation. The accuracy of analyzing left and right ventricular functions by cardiac computed tomography and cardiac magnetic resonate imaging has also been established. These modalities can provide myocardial tissue information and allow calculation of the extracellular volume fraction as well. Summary: According to the rapid improvement of technologies, cardiac imaging has become an essential tool for hemodynamic evaluation in HF management.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Monitoring of asteroids in cometary orbits and activated asteroids through archival images and new observations

Transitional objects are minor bodies that share some characteristics with asteroids and others with comets. These objects include asteroids in cometary orbits (ACOs), which behave dynamically like comets, but lack observed activity, while activated asteroids (AAs) follow typical asteroidal orbits, but have shown dust ejections. The monitoring of a set of these objects carried out in 2015 and 2016 is continued using archival images from various observatories and new data from the IMPACTON telescope in Brazil. Two techniques were applied to detect activity: (i) surface brightness profiles were compared with those of field stars to identify widening, and (ii) the magnitudes reported in the Minor Planet Center, combined with our observations, were reduced and analyzed to identify abrupt brightness increases as a function of heliocentric distance. We analyzed the surface brightness profiles of 133 ACOs and 7 AAs. To study the reduced magnitude, we obtained data from the 705 ACOs that were known at the time of the analysis. Together with the data from our previous work, our analysis covered 23% of the total known ACOs; 8 deviated slightly in the surface brightness profile, 6 brightened in the reduced magnitude, and one object is in common in both samples. A very low percentage of objects might show activity (4% of the sample with brightness profiles and <1% in the reduced magnitudes). These results would rule out a slow transition from active to inert. Regarding AAs, 4 showed activity, and 3 of them matched previously reported periods, while the data we analyzed for P/2015 X 6 were obtained 19 days before the first existing activity report. The activity episodes of these objects are very restricted in time and do not always occur in the same region of the orbit.

asteroids: general↗
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2021

Resource selection functions based on hierarchical generalized additive models provide new insights into individual animal variation and species distributions

Habitat selection studies are designed to generate predictions of species distributions or inference regarding general habitat associations and individual variation in habitat use. Such studies frequently involve either individually indexed locations gathered across limited spatial extents and analyzed using resource selection functions (RSFs) or spatially extensive locational data without individual resolution typically analyzed using species distribution models. Both analytical methodologies have certain desirable features, but analyses that combine individual- and population-level inference with flexible non-linear functions may provide improved predictions while accounting for individual variation. Here, we describe how RSFs can be fit using hierarchical generalized additive models (HGAMs) using widely available software, providing a means to explore individual variation in habitat associations and to generate species distribution maps. We used GPS tracking data from golden eagles Aquila chrysaetos from across eastern North America with four environmental predictors to generate monthly distribution models. We considered three model structures that assumed different amounts of individual variation in the functional relationship between predictors and habitat use and used k-fold cross-validation to compare model performance. Models accounting for individual variability in shape and smoothness of functional responses performed best. Eagles exhibited the least amount of individual variation in response to land cover variables during winter months, with most individuals more closely adhering to the population-level trend. During the summer months, eagles exhibited more substantial individual variation in shape and smoothness of the functional relationships, suggesting some need to account for individual variation in eagle habitat use for both inferential and predictive purposes, during this time of year. Because they allow users to blend flexible functions with random effects structures and are well-supported by a variety of software platforms, we believe that HGAMs provide a useful addition to the suite of analyses used for modeling habitat associations or predicting species distributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Sequential Quadratic Programming Algorithm for Nonsmooth Problems with Upper- \({\boldsymbol{\mathcal{C}^2}}\) Objective

An optimization algorithm for nonsmooth nonconvex constrained optimization problems with upper- \({\boldsymbol{\mathcal{C}^2}}\) objective functions is proposed and analyzed. Upper- \({\boldsymbol{\mathcal{C}^2}}\) is a weakly concave property that exists in difference of convex (DC) functions and arises naturally in many applications, particularly certain classes of solutions to parametric optimization problems e.g., recourse of stochastic programming and projection onto closed sets. The algorithm can be viewed as an extension of sequential quadratic programming (SQP) to nonsmooth problems with upper- \({\boldsymbol{\mathcal{C}^2}}\) objectives or a simplified bundle method. It is globally convergent with bounded algorithm parameters that are updated with a trust-region criterion. The algorithm handles general smooth constraints through linearization and uses a line search to ensure progress. The potential inconsistencies from the linearization of the constraints are addressed through a penalty method. In conclusion, the capabilities of the algorithm are demonstrated by solving both simple upper- \({\boldsymbol{\mathcal{C}^2}}\) problems and a real-world optimal power flow problem used in current power grid industry practices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Analysis of Electric Vehicle Charging Behavior Patterns with Function Principal Component Analysis Approach

This manuscript focused on analyzing electric vehicles’ (EV) charging behavior patterns with a functional data analysis (FDA) approach, with the goal of providing theoretical support to the EV infrastructure planning and regulation, as well as the power grid load management. 5-year real-world charging log data from a total of 455 charging stations in Kansas City, Missouri, was used. The focuses were placed on analyzing the daily usage occupancy variability, daily energy consumption variability, and station-level usage variability. Compared with the traditional discrete-based analysis models, the proposed FDA modeling approach had unique advantages in preserving the smooth function behavior of the data, bringing more flexibility in the modeling process with little required assumptions or background knowledge on independent variables, as well as the capability of handling time series data with different lengths or sizes. In addition to the patterns revealed in the EV charging station’s occupancy and energy consumption, the differences between EV driver’s charging time and parking time were analyzed and called for the needs for parking regulation and enforcement. The different usage patterns observed at charging stations located on different land-use types were also analyzed.

Engineering↗
DOE OSTI2025

Growth functions of periodic space tessellations

This work analyzes the rules governing the growth of the numbers of vertices, edges and faces in all possible periodic tessellations of the 2D Euclidean space, and encodes those rules in several types of polynomial growth functions. These encodings map the geometric, combinatorial and topological properties of the tessellations into sets of integer coefficients. Several general statements about these encodings are given with rigorous mathematical proof. The variation of the growth functions is represented graphically and analyzed in orphic diagrams, so named because of their similarity to orphic art. Several examples of 3D space groups are included, to emphasize the complexity of the growth functions in higher dimensions. A freely available Python library is presented to facilitate the discovery of the growth functions and the generation of orphic diagrams.

Chemistry↗
DOE OSTI2024

Simulating ‘Two Ribbon’ Type Solar Flares [Slides]

The underlying mechanism in solar flares in a process known as magnetic reconnection. This astrophysical process is when magnetic field lines, running anti-parallel, break and reconnect. The change in configuration of field lines results in an explosive release of energy as the leftover magnetic energy is converted to kinetic and thermal energies. In this project, we are using an astrophysical magnetohydrodynamic (MHD) simulation code known as Athena++ with a reconnection problem generator file created by Li et al. 2018. Our current research is to successfully implement a radiative cooling term into the MHD equations as it could have important physical effects on the plasma. We use Klimchuk et al. 2008 and SPEX_DM as our chosen cooling functions. The results from the SPEX_DM function show there is a condensation forming that we had predicted but had not seen before with our simulation code. We will continue to analyze the SPEX_DM function and potentially implement thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗