Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function approximation”

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

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

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

Collisional Broadening of CO{sub 2} Transition Lines 10{sup 0}0–00{sup 0}1 by Inert Gas Atoms at 300–700 K

Using a tunable CO{sub 2} laser, unsaturated absorption coefficients in pure CO{sub 2} and binary CO{sub 2} gas mixtures with inert gases He, Ar, Kr, and Xe in a 1 : 2 ratio at 100 Torr were measured at the central frequencies of the R(8), R(22), P(8), P(22), and P(32) CO{sub 2} spectral lines of the 10{sup 0}0-00{sup 0}1 transition at 300–700 K. The coefficients of the collisional broadening of the CO{sub 2} spectral lines by He, Ar, Kr, and Xe are determined. The obtained temperature dependences of the collisional broadening coefficients are approximated by power functions with two different degrees.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Photoionization of the Bound Systems at High Energies

We consider photoionization of a system bound by the central potential V(r). We demonstrate that the high energy nonrelativistic asymptotics of the photoionization cross section can be obtained without solving the wave equation. The asymptotics can be expressed in terms of the Fourier transform of the potential by employing the Lippmann–Schwinger equation. We find the asymptotics for the screened Coulomb field. We demonstrate that the leading corrections to this asymptotics are described by the universal factor. The high energy nonrelativistic asymptotics is found to be determined by the analytic properties of the potential V(r). We show that the energy dependence of the asymptotics of photoionization cross sections of fullerenes is to large extent model-dependent. We demonstrate that if the fullerene field V(r) is approximated by the function with singularities in the complex plane, the power drop of the asymptotics is reached at the energies which are so high that the cross section becomes unobservably small. The preasymptotic behavior with a faster decrease of the cross sections becomes important in these cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Resolved photons in S HERPA

We present the first complete simulation framework, in the SHERPA event generator, for resolved photon interactions at next-to leading order accuracy. It includes photon spectra obtained through the equivalent-photon approximation, parton distribution functions to parametrize the hadronic structure of quasi-real photons, the matching of the parton shower to next-to leading order QCD calculations for resolved photon cross sections, and the modelling of multiple-parton interactions. We validate our framework against a wide range of photo-production data from LEP and HERPA experiments, observing good overall agreement. We identify important future steps relevant for high-quality simulations at the planned Electron–Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Manipulating Light with Transition Metal Clusters, Organic Dyes, and Metal Organic Frameworks

The primary goals of our research program is to develop and apply state-of-the-art first-principles methods to predict electronic and optical properties of three systems of significant scientific and technological interest: transition metal clusters, organic dyes, and metal-organic frameworks. These systems offer great opportunities to manipulate light for a wide ranging list of energy-related scientific problems and applications. During this grant period, we focused our investigations on the development, implementation, and benchmarking of many-body Green’s function methods (GW approximation and the Bethe-Salpeter equation) to examine excited-state properties of transition metal/transition-metal-oxide clusters and organic molecules that comprise the building blocks of dyes and metal-organic frameworks.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Distributed fiber sensing systems for 3D combustion temperature field monitoring in coal-fired boilers using optically generated acoustic waves (Final Report)

In this project, we have developed and tested three kinds of fiber optic sensing systems for real time monitoring of temperature variations within an industrial scale boiler furnace. The fiber optic sensing systems target spatial and temporal distributions of high temperature profiles in a boiler furnace in fossil power plants. The reconstructed temperature profile will provide critical input for the control mechanisms to optimize the combustion process. This temperature profile will address the essential problem for fossil power plants in achieving higher efficiency and fewer pollutant emissions. Acoustic pyrometer systems have been used to reconstruct temperature field of power plant boilers based on measuring TOF (times-of-flight) of sound waves along some straight paths in a 2D cross-section of the boiler. In this project, optically generated acoustic signals from a fiber optic sensing system have replaced the acoustic signals generated from an electrical transducer. A 3D reconstruction algorithm replaced the previous 2D model. In this project, three kinds of fiber optic sensing systems have been developed and tested. They are fiber optic sensing system I, fiber optic sensing system II (Distributed Sensing System I) and fiber optic sensing system III (Distributed Sensing System II). For fiber optic sensing system I, the fiber optic ultrasound generator acts as a signal generator. A microphone, hydrophone or other electronic devices serve as a signal receiver. In this system, there are one generator and one receiver. Distance test, water temperature test, air temperature test, air temperature reconstruction, and GE ISBF pilot test were performed by Fiber optic sensing system I. The fiber optic sensing system I successfully detected temperature in all these tests. We got 2D temperature reconstruction results by using the fiber optic sensing system I and it matched the reference data. The fiber optic sensing system I successfully survived in GE ISBF boiler environment (480 °F). For fiber optic sensing system II (Distributed Sensing System I), it is an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Fiber Bragg Grating (FBG) and Fabry-Perot (FP) sensor act as a signal receiver. In this system, there is one generator and one receiver. Aluminum plate temperature test, furnace high temperature test, and GE ISBF pilot test were performed by the fiber optic sensing system II. Fiber optic sensing system II successfully detected the temperature in all these tests. The fiber optic sensing system II successfully survived at up to 700 °C furnace environment and 320 °C GE ISBF boiler environment. For fiber optic sensing system III (Distributed Sensing System II), it is also an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Multiple FP fiber sensors act as signal receivers. In this system, there is one generator and three receivers. Three GE ISBF pilot tests were performed by fiber optic sensing system III. The fiber optic sensing system III survived in the cold flow tests in GE’s ISBF pilot test facility. However, we didn’t get high temperature data by using this system since the nanosecond laser issues. During the period of the project, test trials, data simulation and algorithm optimization was performed successfully. For real time temperature field construction, the sampling rate must be fast enough to capture the field variations. The technology of Code-division multiple access (CDMA) is well studied which could allow parallel multiplexing, even if signals overlap in time or frequencies. Moreover, it has been known that extending the length of signal significantly improves SNR. For acoustic signals, these multiplexing techniques have also been widely used, mainly for sonar and acoustic communications. The CDMA modulation technique has been proposed and studied to guarantee high network throughput, low channel access delay and low energy consumption. We have studied the temperature field reconstruction using Gaussian Radial Basis Functions (GRBF)-based approximation approach. Reconstruction of 3D temperature field using Neural Networks with measured TOF and known propagation paths is feasible. 2D and 3D temperature field reconstruction simulation results are achieved. The milestone status is shown in Table 1. We finished milestone 1-8 and milestone 10. For milestone 9, we did three pilot tests by using the fiber optic sensing system III (Distributed Sensing System II) at GE Power. However, due to the failure of the ns laser, we did not get the temperature results. We conducted some additional tasks that were not originally proposed: 1) We fabricated a fiber optic sensing system I and did a pilot test based on this system. 2) In the proposal, we proposed two pilot tests at GE Power. In reality, we finished at least seven pilot tests at GE Power. GE Power has made a lot of efforts for supporting the pilot tests. 3) We got a simulation results based on CDMA. In summary, most of the tasks have been accomplished. The outcome of this project removed a few barriers that hinder the achievement of the final product of the distributed sensing systems. With the successful accomplishment of this project, a prototype of the fiber optic sensing system can be fabricated to attract more interests from companies and other funding agencies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

HAIMOS Ensemble Forecasts for Intra-day and Day- Ahead GHI, DNI and Ramps

The objective of this research is to develop a hybrid physics-based/data-driven forecast model to improve direct normal and global horizontal irradiance (DNI and GHI) prediction for horizons ranging from 1 to 72 hours. Project objectives also address key gaps in state-of-the-art solar forecasting: accurate probabilistic solar forecasts and the forecasting of large irradiance ramps (ramp onset and magnitude). The proposed model ensembles Numerical Weather Prediction (NWP) forecasts, determinist physics-based algorithms, and new-generation cloud cover products (high-resolution rapid refresh satellite images and Large Eddy Simulations). The result is the Hybrid Adaptive Input Model Objective Selection (HAIMOS) ensemble model. HAIMOS blends state of the art machine learning methodologies with physics-based models for cloud cover and cloud optical depth forecasts. The technical activities followed a two-pronged strategy. First, the preprocessing of data, the selection of inputs to the nonlinear approximators, the type of approximator and objective functions, and post-processing ensembling techniques included in HAIMOS were all optimized adaptively to find the best model for a specific goal (reduce DNI/GHI forecast error, improve the prediction of ramp onset, etc.). Second, a large effort was put in improving cloud identification and the forecast of cloud cover and cloud optical depth. To this end, new-generation cloud parametrization products were developed in this work. These include improved algorithms to assist in cloud identification, cloud classification and cloud parametrization from satellite images – three key factors in the accuracy of 1 to 6-hours irradiance forecasts and prediction of ramp onset. Furthermore, we also included cloud information extracted from high resolution rapid refresh satellite images (GOES-16) and Large Eddy Simulations (LES). LES was used to model the atmosphere in detail over locations of interest and produce cloud optical depth forecasts. Once these data streams were validated, they were used as input data to the HAIMOS forecast. The model was developed using data from several climatologically distinct locations with potential for high solar penetration. In the last year of the project, we conducted a validation campaign according to the guidelines stipulated by the Topic Area 1 project as described in the FOA. This effort brings, for the first time, proven machine-learning methodologies for generating state-of-the-art solar forecasts interweaved with detailed physics-based models for cloud detection, and cloud optical depth forecasts. HAIMOS will generate accurate irradiance probabilistic forecast to assist in reducing solar generation prediction error. Globally optimized solar forecast models are more likely to impact solar energy stakeholders. The goal of this project was to increase the state-of-the-art forecast skill from their present values of 10 to 35%. At the end of the project, we achieved between 30% and 50% forecast skill across a wide range of horizons for both GHI and DNI.

14 SOLAR ENERGY↗
DOE OSTI2024

PhILMs: Collaboratory on Mathematics and Physics-Informed Learning Machines for Multiscale and Multiphysics Problems

The landscape of computational science and engineering is continually evolving, with the challenge of high-dimensional regression problems standing as a significant hurdle in numerous scientific endeavors. Addressing this challenge, our research, funded by this award, has led to the development of an innovative computational framework known as Probabilistic Partition of Unity Networks (PPOU-Nets). This initiative represents a collaborative effort to harness the potential of mathematics and physics-informed machine learning in tackling multiscale and multiphysics problems prevalent in high-dimensional spaces. Through this work, we have proposed a novel methodology that seamlessly integrates adaptive dimensionality reduction and a mixture of experts model, thereby facilitating a more efficient and accurate approximation of complex functions. This research effort has not only advanced the state of computational science but also opened new avenues for exploration in quantum computing and beyond. This report outlines the motivation, methodology, key findings, and implications of our work, underscoring our contributions to the broader scientific community and the potential pathways for future research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Solar and Storage Integration in the U.S. Southeast: Implications for Resource Adequacy [Slides]

In this study, we evaluate how an iterative portfolio approach compares to a more traditional capacity credit method, using the U.S. Southeast as a case study region. Using open-source planning and resource adequacy tools, we compare results using a capacity credit approximation method with those from iterating behind the two. We also explore how the iterative approach performs under a range of sensitivities, including higher load growth, regional coordination, and alternative weather years. We find that traditional capacity credit approximation methods can function well in today's system, but may face challenges for systems with higher levels of solar and storage. As such, integrating planning and resource adequacy models can address some of these gaps, helping planners deliver more reliable systems.

14 SOLAR ENERGY↗
DOE OSTI2017

Effect of Radiation on Chromospheric Magnetic Reconnection: Reactive and Collisional Multi-fluid Simulations

We study magnetic reconnection under chromospheric conditions in five different ionization levels from 0.5% to 50% using a self-consistent two-fluid (ions + neutrals) model that accounts for compressibility, collisional effects, chemical inequilibrium, and anisotropic heat conduction. Results with and without radiation are compared, using two models for the radiative losses: an optically thin radiation loss function, and an approximation of the radiative losses of a plasma with photospheric abundances. The results without radiation show that reconnection occurs faster for the weakly ionized cases as a result of the effect of ambipolar diffusion and fast recombination. The tearing mode instability appears earlier in the low ionized cases and grows rapidly. We find that radiative losses have a stronger effect than was found in previous results as the cooling changes the plasma pressure and the concentration of ions inside the current sheet. This affects the ambipolar diffusion and the chemical equilibrium, resulting in thin current sheets and enhanced reconnection. The results quantify this complex nonlinear interaction by showing that a strong cooling produces faster reconnections than have been found in models without radiation. The results accounting for radiation show timescales and outflows comparable to spicules and chromospheric jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Structure-Preserving Decorated Particle Method for the Vlasov-Poisson System

We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.

65M75, 70H05, 70G65↗
DOE OSTI2017

Modular Expression Language for Ordinary Differential Equation Editing

MELODEEis a system for describing systems of initial value problem ordinary differential equations, and a compiler for the language that produces optimized code to integrate the differential equations. Features include rational polynomial approximation for expensive functions and automatic differentiation for symbolic jacobians

Blake, Robert C.↗
DOE OSTI2019

Monte Carlo Control Loops for Cosmic Shear Cosmology with DES Year 1

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, Point Spread Function (PSF), an Approximate Bayesian Computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and the fast estimation of the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $\Omega_m$, $\sigma_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$, defined as $S_8D_{\rm{IA}} = \sigma_8(\Omega_m/0.3)^{0.5}D_{\rm{IA}}$, where $D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find $S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of roughly 60\% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Constrained curve fitting for semi-parametric models with radial basis function networks

Common to many analysis pipelines in lattice gauge theory and the broader scientific discipline is the need to fit a semi-parametric model to data. We propose a fit method that utilizes a radial basis function network to approximate the non-parametric component of such models. The approximate parametric model is fit to data using the basin hopping global optimization algorithm. Parameter constraints are enforced through Gaussian priors. The viability of our method is tested by examining its use in a finite-size scaling analysis of the $q$-state Potts model and $p$-state clock model with $q=2,3$ and $p=4,\infty$.

Peterson, Curtis T.↗
DOE OSTI2019

Multiplicity dependence of light-flavor hadron production in p p collisions at s = 7 TeV

Comprehensive results on the production of unidentified charged particles, $π^±, K^±, K^0_S, K^*(892)^0, p,$ $\bar{p}$, $\phi$ (1020), Λ, $\bar{Λ}$, $\Xi^-$, $\bar{\Xi}^+$, $Ω^-$, and $\bar{Ω}^+$ hadrons in proton-proton ($pp$) collisions at $\sqrt{s}$ = 7 TeV at midrapidity ($|y|$ < 0.5) as a function of charged-particle multiplicity density are presented. In order to avoid autocorrelation biases, the actual transverse momentum ($p$ T ) spectra of the particles under study and the event activity are measured in different rapidity windows. In the highest multiplicity class, the charged-particle density reaches about 3.5 times the value measured in inelastic collisions. While the yield of protons normalized to pions remains approximately constant as a function of multiplicity, the corresponding ratios of strange hadrons to pions show a significant enhancement that increases with increasing strangeness content. Furthermore, all identified particle-to-pion ratios are shown to depend solely on charged-particle multiplicity density, regardless of system type and collision energy. The evolution of the spectral shapes with multiplicity and hadron mass shows patterns that are similar to those observed in $p$-Pb and Pb-Pb collisions at Large Hadron Collider energies. The obtained $p$ T distributions and yields are compared to expectations from QCD-based $pp$ event generators as well as to predictions from thermal and hydrodynamic models. These comparisons indicate that traces of a collective, equilibrated system are already present in high-multiplicity $pp$ collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

The effect of Limber and flat-sky approximations on galaxy weak lensing

We review the effect of the commonly-used Limber and flat-sky approximations on the calculation of shear power spectra and correlation functions for galaxy weak lensing. These approximations are accurate at small scales, but it has been claimed recently that their impact on low multipoles could lead to an increase in the amplitude of the mass fluctuations inferred from surveys such as CFHTLenS, reducing the tension between galaxy weak lensing and the amplitude determined by Planck from observations of the cosmic microwave background. Here, we explore the impact of these approximations on cosmological parameters derived from weak lensing surveys, using the CFHTLenS data as a test case. We conclude that the use of small-angle approximations for cosmological parameter estimation is negligible for current data, and does not contribute to the tension between current weak lensing surveys and Planck.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

LTAU (Loss Trajectory Analysis for Uncertainty) is a technique for estimating the uncertainty in a model's predictions for any regression task by approximating the cumulative distribution function (CDF) of the model's errors over the course of training. The approximated CDF, combined with a similarity search algorithm in the model's descriptor space, can then be used to estimate the likelihood that the model's predictions on any given test point will fall below a chosen threshold. LTAU-FF (LTAU in atomistic Force Fields) is the application of LTAU specifically for use with atomistic force fields.

Vita, JoshuaA↗