Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functions and mappings”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Modern Era Retrospective Restrospective-Analysis for Research and Applications (MERRA) Data and Services at the GES DISC

The Modern Era Retrospective-analysis for Research and Applications (MERRA) dataset is a NASA satellite era, 30 year (1979 - present), reanalysis using the Goddard Earth Observing System Data Assimilation System, Version 5 (GEOS-5). The project, run out of NASA's Global Modeling and Assimilation Office at Goddard Space Flight Center, provides the science and application communities with a state-of-the-art global analysis with emphasis on improved estimates of the hydrological cycle over a broad range of weather and climate time scales. MERRA products are generated as a long-term synthesis that places the NASA EOS suite of observations in a climate context. The MERRA analysis is performed at a horizontal resolution of 2/3 longitude x 1/2 latitude (540x361 global gridpoints) with observational analyses every 6 hours. The MERRA output data will include 3 dimensional state fields for every 6 hourly analysis cycle on 42 pressure levels (or 72 terrain following model coordinate levels) from the surface through the stratosphere. Several data products are specifically designed to support chemistry and stratosphere transport modeling. The 2 dimensional surface and atmospheric diagnostics (numbering 259) are being stored on the native grid at 1 hourly intervals. These include radiation and vertical integrals of the atmosphere for water and energy budget studies and also surface diagnostics where the diurnal cycle is important. The one hourly surface and near surface data product will also facilitate research on the integrated analysis of Earth system observations in the land, ocean and cryosphere. The MERRA products are archived and distributed by the Goddard Earth Sciences Data and Information Services Center (GES DISC) through its Modeling DISC Web (MDISC) portal. Multiple data access methods and services are available for MERRA data through MDISC: (1) Mirador offers a quick, comprehensive search of MERRA and all GES DISC archived data holdings, allowing searches on keywords, location names or latitude/longitude box, and date/time, with responses within a few seconds. (2) Giovanni is a GES DISC developed Web application that provides data visualization and analysis online. Giovanni features popular visualizations such as latitude-longitude maps, animations, cross sections, profiles, time series, etc. and some basic statistical analysis functions such as scatter plots and correlation coefficient maps. Users are able to download results in several different formats, including Google Earth. (3) On-the-fly parameter subsetting of data within a spatial/temporal window is provided through a simple select and click Web page. (4) MERRA data are also available via OPeNDAP, GrADS Data Server (GDS) and can be converted to netCDF on the fly.

Berrick, Stephen W.↗
DOE OSTI2023

On the influence of over-parameterization in manifold based surrogates and deep neural operators

Constructing accurate and generalizable approximators (surrogate models) for complex physico-chemical processes exhibiting highly non-smooth dynamics is challenging. The main question is what type of surrogate models we should construct and should these models be under-parameterized or over-parameterized. In this work, we propose new developments and perform comparisons for two promising approaches: manifold-based polynomial chaos expansion (m-PCE) and the deep neural operator (DeepONet), and we examine the effect of over-parameterization on generalization. While m-PCE enables the construction of a mapping by first identifying low-dimensional embeddings of the input functions, parameters, and quantities of interest (QoIs), a neural operator learns the nonlinear mapping via the use of deep neural networks. Here, we demonstrate the performance of these methods in terms of generalization accuracy by solving the 2D time-dependent Brusselator reaction-diffusion system with uncertainty sources, modeling an autocatalytic chemical reaction between two species. We first propose an extension of the m-PCE by constructing a mapping between latent spaces formed by two separate embeddings of the input functions and the output QoIs. To further enhance the accuracy of the DeepONet, we introduce weight self-adaptivity in the loss function. We demonstrate that the performance of m-PCE and DeepONet is comparable for cases of relatively smooth input-output mappings. However, when highly non-smooth dynamics is considered, DeepONet shows higher approximation accuracy. We also find that for m-PCE, modest over-parameterization leads to better generalization, both within and outside of distribution, whereas aggressive over-parameterization leads to over-fitting. In contrast, an even highly over-parameterized DeepONet leads to better generalization for both smooth and non-smooth dynamics. Furthermore, we compare the performance of the above models with another recently proposed operator learning model, the Fourier Neural Operator, and show that its over-parameterization also leads to better generalization. Taken together, our studies show that m-PCE can provide very good accuracy at very low training cost, whereas a highly over-parameterized DeepONet can provide better accuracy and robustness to noise but at higher training cost. In both methods, the inference cost is negligible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Kernel Manifolds: Nonlinear‐Augmentation Dimensionality Reduction Using Reproducing Kernel Hilbert Spaces

This paper generalizes recent advances on quadratic manifold (QM) dimensionality reduction by developing kernel methods-based nonlinear-augmentation dimensionality reduction. QMs, and more generally feature map-based nonlinear corrections, augment linear dimensionality reduction with a nonlinear correction term in the reconstruction map to overcome approximation accuracy limitations of purely linear approaches. While feature map-based approaches typically learn a least squares optimal polynomial correction term, we generalize this approach by learning an optimal nonlinear correction from a user-defined reproducing kernel Hilbert space. Our approach allows one to impose arbitrary nonlinear structure on the correction term, including polynomial structure, and includes feature map and radial basis function-based corrections as special cases. Furthermore, our method has relatively low training cost and has monotonically decreasing error as the latent space dimension increases. In conclusion, we compare our approach to proper orthogonal decomposition and several recent QM approaches on data from several example problems.

kernel methods↗
DOE OSTI2023

Grain boundary softening from stress assisted helium cavity coalescence in ultrafine-grained tungsten

The formation of helium cavities in coarse-grained materials produces hardening proportional to the number density and size of the cavities and due to the interaction of dislocations with intragranular helium defects. In nanostructured metals containing a high density of interfacial sinks, preferential cavity formation in the grain boundaries instead produces softening that is often attributed to enhanced interfacial plasticity. Here, employing two grades of ultrafine-grained tungsten, we explore this effect using targeted implantation studies to map cavity evolution as a function of the irradiation conditions and quantify its impact on the mechanical response through nanoindentation. Softening is reported at implantation temperatures above the threshold for preferential grain boundary cavity formation but at a sufficiently low fluence prior to the growth of intragranular cavities. Collective changes in the mean cavity size, density, and morphology beneath a residual impression on an implanted surface indicate that cavity coalescence accompanied the reduction in hardness. Complementary atomistic simulations demonstrate that, in tungsten grain structures exhibiting softening, grain boundary bubble coalescence is driven by stress concentrations that further act to localize strain in the grain boundaries through cooperative deformation processes involving local atomic shuffling and sliding, dislocation emission, and even the nucleation of unstable twinning events.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reliable extrapolation of deep neural operators informed by physics or sparse observations

Deep neural operators can learn nonlinear mappings between infinite-dimensional function spaces via deep neural networks. As promising surrogate solvers of partial differential equations (PDEs) for real-time prediction, deep neural operators such as deep operator networks (DeepONets) provide a new simulation paradigm in science and engineering. Pure data-driven neural operators and deep learning models, in general, are usually limited to interpolation scenarios, where new predictions utilize inputs within the support of the training set. However, in the inference stage of real-world applications, the input may lie outside the support, i.e., extrapolation is required, which may result to large errors and unavoidable failure of deep learning models. Here, we address this challenge of extrapolation for deep neural operators. First, we systematically investigate the extrapolation behavior of DeepONets by quantifying the extrapolation complexity, via the 2-Wasserstein distance between two function spaces and propose a new strategy of bias–variance trade-off for extrapolation with respect to model capacity. Subsequently, we develop a complete workflow, including extrapolation determination, and we propose five reliable learning methods that guarantee a safe prediction under extrapolation by requiring additional information—the governing PDEs of the system or sparse new observations. The proposed methods are based on either fine-tuning a pre-trained DeepONet or multifidelity learning. We demonstrate the effectiveness of the proposed framework for various types of parametric PDEs. Furthermore, our systematic comparisons provide practical guidelines for selecting a proper extrapolation method depending on the available information, desired accuracy, and required inference speed.

42 ENGINEERING↗
DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗
DOE OSTI2025

Dual-season thermal energy storage-integrated heat pump system for load shifting: experiments, modeling and nationwide analysis

Here, this paper proposes and validates a dual-season thermal energy storage-integrated heat pump (TES-HP) system that shifts building electric load in both cooling and heating modes using a single phase-change TES. The design uses the TES as a heat sink for the condenser in summer and a heat source for the evaporator in winter, thereby reducing on-peak electric demand while ensuring compatible with existing air-distribution systems. A 14-kW TES-HP prototype with six operating modes and novel refrigerant charge management strategies was experimentally tested under varied conditions. Subsequently, data-driven polynomial performance curves were developed and validated against experiments, then coupled with a phase change material (PCM) model and a calibrated building model within a rule-based controller. Single-day case studies demonstrated effective on-peak demand reduction while maintaining thermal comfort. During cooling, hourly power was lowered by up to 1.5 kW; in heating, average hourly power decreased from over 3.5 kW to 1.8 kW, compared to an HP-only system, achieving 64.5 % load shifting, and the need for resistance heating was eliminated. Seasonal simulations showed typical on-peak electricity savings of 25–35 % in cooling and 40–65 % in heating, with the largest benefits on the hottest and coldest days. Extended response-surface analysis and nationwide mapping quantified load shifting as functions of ambient and TES temperatures, revealing a seasonal trade-off in TES phase-change temperature. These results demonstrate the TES-HP system as a practical and scalable solution for grid-interactive buildings that reduces on-peak demand and electricity use while maintaining thermal comfort and enhancing grid flexibility and reliability.

Grid-interactive buildings↗
DOE OSTI2025

Gapless dispersive continuum in a modulated quantum kagome antiferromagnet

The pursuit of quantum spin liquid (QSL) states in condensed matter physics has drawn attention to kagome antiferromagnets (AFM) where a two-dimensional corner-sharing network of triangles frustrates conventional magnetic orders. While quantum kagome AFMs based on Cu 2+ (3d 9 , s = 1/2) ions have been extensively studied, there is so far little work beyond copper-based systems. Here we present our bulk magnetization, specific heat and neutron scattering studies on single crystals of a new titanium fluorides Cs 8 RbK 3 Ti 12 F 48 where Ti 3+ (3d 1 , s = 1/2) ions form a modulated quantum kagome antiferromagnet that does not order magnetically down to 1.5 K. Our comprehensive map of the dynamic response function $S$($Q$,$\hslash$$\omega$) acquired at 1.5 K where the heat capacity is T-linear reveals a dispersive continuum emanating from soft lines that extend along (100). The data indicate fractionalized spinon-like excitations with quasi-one-dimensional dispersion within a quasi-two-dimensional spin system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Transient pattern formation in an active matter contact poisoning model

Abstract One of the most notable features in repulsive particle based active matter systems is motility-induced-phase separation (MIPS) where a dense, often crystalline phase and low density fluid coexist. Most active matter studies involve time-dependent activity; however, there are many active systems where individual particles transition from living or moving to dead or nonmotile due to lack of fuel, infection, or poisoning. Here we consider an active matter particle system at densities where MIPS does not occur. When we add a small number of infected particles that can poison other particles, rendering them nonmotile, we find a rich variety of time dependent pattern formation, including MIPS, a wetting phase, and a fragmented state formed when mobile particles plow through a nonmotile packing. We map the patterns as a function of time scaled by epidemic duration, and show that the pattern formation is robust for a wide range of poisoning rates and activity levels. We also show that pattern formation does not occur in a random death model, but requires the promotion of nucleation by contact poisoning. Our results should be relevant to biological and active matter systems where there is some form of poisoning, death, or transition to nonmotility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Fluctuating charge-density-wave correlations in the three-band Hubbard model

The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model’s fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.

42 ENGINEERING↗
DOE OSTI2022

Directed motion of liquid crystal skyrmions with oscillating fields

Abstract Using continuum simulations, we show that under a sinusoidal electric field, liquid crystal skyrmions undergo periodic shape oscillations which produce controlled directed motion. The speed of the skyrmion is non-monotonic in the frequency of the applied field, and exhibits multiple reversals of the motion as a function of changing frequency. We map out the dynamical regime diagram of the forward and reverse motion for two superimposed ac driving frequencies, and show that the reversals and directed motion can occur even when only a single ac driving frequency is present. Using pulsed ac driving, we demonstrate that the motion arises due to an asymmetry in the relaxation times of the skyrmion shape. We discuss the connection between our results and ratchet effects observed in systems without asymmetric substrates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Mapping variations of redshift distributions with probability integral transforms

We present a method for mapping variations between probability distribution functions and apply this method within the context of measuring galaxy redshift distributions from imaging survey data. This method, which we name PITPZ for the probability integral transformations it relies on, uses a difference in curves between distribution functions in an ensemble as a transformation to apply to another distribution function, thus transferring the variation in the ensemble to the latter distribution function. This procedure is broadly applicable to the problem of uncertainty propagation. In the context of redshift distributions, for example, the uncertainty contribution due to certain effects can be studied effectively only in simulations, thus necessitating a transfer of variation measured in simulations to the redshift distributions measured from data. We illustrate the use of PITPZ by using the method to propagate photometric calibration uncertainty to redshift distributions of the Dark Energy Survey Year 3 weak lensing source galaxies. For this test case, we find that PITPZ yields a lensing amplitude uncertainty estimate due to photometric calibration error within 1 per cent of the truth, compared to as much as a 30 per cent underestimate when using traditional methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Photometry and astrometry with JWST – III. A NIRCam-Gaia DR3 analysis of the open cluster NGC 2506

ABSTRACT In the third paper of this series aimed at developing the tools for analysing resolved stellar populations using the cameras on board of the James Webb Space Telescope (JWST), we present a detailed multiband study of the 2 Gyr Galactic open cluster NGC 2506. We employ public calibration data sets collected in multiple filters to: (i) derive improved effective Point Spread Functions (ePSFs) for 10 NIRCam filters; (ii) extract high-precision photometry and astrometry for stars in the cluster, approaching the main sequence (MS) lower mass of ∼0.1 M⊙; and (iii) take advantage of the synergy between JWST and Gaia DR3 to perform a comprehensive analysis of the cluster’s global and local properties. We derived a MS binary fraction of ∼57.5 per cent, extending the Gaia limit (∼0.8 M⊙) to lower masses (∼0.4 M⊙) with JWST. We conducted a study on the mass functions (MFs) of NGC 2506, mapping the mass segregation with Gaia data, and extending MFs to lower masses with the JWST field. We also combined information on the derived MFs to infer an estimate of the cluster present-day total mass. Lastly, we investigated the presence of white dwarfs (WDs) and identified a strong candidate. However, to firmly establish its cluster membership, as well as that of four other WD candidates and of the majority of faint low-mass MS stars, further JWST equally deep observations will be required. We make publicly available catalogues, atlases, and the improved ePSFs.

Nardiello, D. (ORCID:0000000311493659)↗
DOE OSTI2025

Seedless nonresonant gas-flow velocimetry with single-shot coherent Rayleigh-Brillouin scattering

We demonstrate neutral gas-flow velocimetry by exploiting the optical dipole forces exerted by precisely tailored high-intensity laser fields on polarizable particles, in the form of single-shot coherent Rayleigh-Brillouin scattering (CRBS). This enables mapping of the velocity distribution function of the particles in a single laser shot of approximately 200 ns duration, allowing us to perform spectroscopic measurements in the time, rather than the frequency, domain. Importantly, CRBS results in a coherent signal beam, which renders the technique ideal for measurements in optically noisy environments. By correlating the time shift of the resulting spectrum to the Doppler shift caused by the bulk motion of particles, we perform neutral gas-flow velocimetry without the need to seed the flow, nondestructively, for a variety of gases, atomic or molecular. We demonstrate gas-flow velocimetry at flow velocities as low as 1 m ⁢s −1 , with similar velocity resolution, while temperature and gas density can also be estimated from the same single-shot spectra. To prove the effectiveness of the measurement technique, we benchmark it against Pitot tube measurements, which show excellent agreement.

Gerakis, Alexandros [Luxembourg Institute of Scien↗
DOE OSTI2023

Cosmological constraints from the Subaru Hyper Suprime-Cam year 1 shear catalogue lensing convergence probability distribution function

Here we utilize the probability distribution function (PDF) of normalized convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) year 1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude S 8 = σ 8 $\sqrt {Ω_{m}/0.3}$. The large-scale structure’s statistical properties are incompletely described by the traditional two-point statistics, motivating our investigation of the PDF—a complementary higher-order statistic. By defining the PDF over the standard-deviation-normalized convergence map, we are able to isolate the non-Gaussian information. We use tailored simulations to compress the data vector and construct a likelihood approximation. We mitigate the impact of survey and astrophysical systematics with cuts on smoothing scales, redshift bins, and data vectors. We find S 8 = $0.860^{+0.066}_{–0.109}$ from the PDF alone and S 8 = 0.798$^{+0.029}_{–0.042}$ from the combination of the PDF and power spectrum (68% confidential level (CL)). The PDF improves the power-spectrum-only constraint by about 10%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Development of a machine-learning-based ionic-force correction model for quantum molecular dynamic simulations of warm dense matter

In this work Δ learning is used to map orbital-free density functional theory (OF-DFT) ionic forces to the corresponding Kohn-Sham (KS) DFT ionic forces. The development of the approximate force difference in terms of the ion positions is constructed and serves as a stand in for the ground truth force difference. Descriptor vectors for ion configurations are constructed using all distance between ions in conjunction with an indexing based on a nearest neighbor ranking. It is demonstrated that such a scheme of descriptors can uniquely describe an ionic configuration up to a rotation and reflection when no ambiguity in the nearest neighbor ranking exists. How to handle the case when an ambiguity exists in the nearest neighbor ranking is discussed. As a proof of principle, the model is trained and tested on warm dense hydrogen at temperatures between 1 and 15 eV. Once tested, the model was used to perform molecular dynamic simulations of warm dense hydrogen. Furthermore, the resulting energies and pressures are within 1% and 2% of their respective target KS values.

36 MATERIALS SCIENCE↗
DOE OSTI2022

MIONet: Learning Multiple-Input Operators via Tensor Product

As an emerging paradigm in scientific machine learning, neural operators aim to learn operators, via neural networks, that map between infinite-dimensional function spaces. Several neural operators have been recently developed. However, all the existing neural operators are only designed to learn operators defined on a single Banach space; i.e., the input of the operator is a single function. Here, for the first time, we study the operator regression via neural networks for multiple-input operators defined on the product of Banach spaces. We first prove a universal approximation theorem of continuous multiple-input operators. We also provide a detailed theoretical analysis including the approximation error, which provides guidance for the design of the network architecture. Based on our theory and a low-rank approximation, we propose a novel neural operator, MIONet, to learn multiple-input operators. MIONet consists of several branch nets for encoding the input functions and a trunk net for encoding the domain of the output function. Here, we demonstrate that MIONet can learn solution operators involving systems governed by ordinary and partial differential equations. In our computational examples, we also show that we can endow MIONet with prior knowledge of the underlying system, such as linearity and periodicity, to further improve accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Compressor is A Sensor: A Universal Refrigerant Charge Fault Detection and Diagnostics Method Based on Pump Down Operation

The primary goal of this report is to develop a universal charge fault detection method that requires only a few experimental data with high prediction accuracy. Currently, pump-down operations is typical practices by HVAC technicians when they need to open the refrigerant circuit to make a repair. In addition, compressors have a built-in low-pressure cut-off protection function, and the compressor performance maps are commonly available from manufacturers. The new charge fault detection and diagnostics method innovatively utilizes the typical pump down operation, the compressor’s low-pressure cut-off protection, and the compressor performance map. It does not require any geometry information of heat exchangers, refrigerant lines, or charge buffers.

42 ENGINEERING↗