SEARCH · Engineering Papers

Results for “interpolation”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Dimensional Interpolation for Random Walk

In this work, we employ a simple and accurate dimensional interpolation formula for the shapes of random walks at D = 3 and D = 2 based on the analytically known solutions at both limits D = ∞ and D = 1. The results obtained for the radius of gyration of an arbitrary shaped object have about 2% error compared with accurate numerical results at D = 3 and D = 2. We also calculated the asphericity for a three-dimensional random walk using the dimensional interpolation formula. The results agree very well with the numerically simulated results. The method is general and can be used to estimate other properties of random walks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Potential energy surface fitting by a statistically localized, permutationally invariant, local interpolating moving least squares method for the many-body potential: Method and application to N{sub 4}

Fitting potential energy surfaces to analytic forms is an important first step for efficient molecular dynamics simulations. Here, we present an improved version of the local interpolating moving least squares method (L-IMLS) for such fitting. Our method has three key improvements. First, pairwise interactions are modeled separately from many-body interactions. Second, permutational invariance is incorporated in the basis functions, using permutationally invariant polynomials in Morse variables, and in the weight functions. Third, computational cost is reduced by statistical localization, in which we statistically correlate the cutoff radius with data point density. We motivate our discussion in this paper with a review of global and local least-squares-based fitting methods in one dimension. Then, we develop our method in six dimensions, and we note that it allows the analytic evaluation of gradients, a feature that is important for molecular dynamics. The approach, which we call statistically localized, permutationally invariant, local interpolating moving least squares fitting of the many-body potential (SL-PI-L-IMLS-MP, or, more simply, L-IMLS-G2), is used to fit a potential energy surface to an electronic structure dataset for N{sub 4}. We discuss its performance on the dataset and give directions for further research, including applications to trajectory calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

The Grassmann interpolation method for spin-unrestricted open-shell systems

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis.

Chemistry↗

DOE OSTI2025

Kinematically enhanced interpolating operators for boosted hadrons

We propose to use interpolating operators for lattice quantum chromodynamics calculations of highly boosted pions and nucleons with kinematically enhanced ground-state overlap factors at large momentum. Because this kinematic enhancement applies to the signal but not the variance of the correlation function, these interpolating operators can achieve better signal-to-noise ratios at large momentum. We perform proof-of-principle calculations with boosted pions and nucleons using close-to-physical and larger quark masses to explore the utility of our proposal. Results for effective energies and matrix elements, as well as Lanczos ground-state energy estimators, are consistent with theoretical expectations for signal-to-noise improvement at large momenta.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Chiral properties of the nucleon interpolating current and θ-dependent observables

We revisit the chiral properties of nucleon interpolating currents, and show that of the two leading order currents j 1 and j 2 , only two linear combinations j 1 ± j 2 transform covariantly under the anomalous U⁢(1) A symmetry. As a result, calculations of quantities which vanish by symmetry in the chiral limit may produce unphysical results if carried out with different linear combinations of the currents. This includes observables such as electric dipole moments, induced by the quantum chromodynamics (QCD) parameter θ, and the θ-dependence of the nucleon mass. For completeness, we also exhibit the leading order results for nucleon electric dipole moments (d n,p ) induced by θ, and the nucleon magnetic moments (μ n,p ), when calculated using QCD sum rules for both the covariant choices of the nucleon interpolating current. The results in each channel, conveniently expressed as the ratios, d n,p /μ n,p , are numerically consistent, and reflect the required physical dependence on θ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Tuning the Interpolation Basis in a Multigrid Decomposition for Local Error Control

In the compression of scientific data, error-controlled compressors enable to considerably decrease the size of the dataset while maintaining adequate levels of accuracy. In this paper, we note that multi-level refactoring scheme such as MGARD i) rely on an approximation of the data based on the interpolation of coefficients, ii) estimate the resulting error with global metrics on the dataset. To improve on these two aspects, we propose a method that aims to divide the original dataset into blocks based on their smoothness and refactors each block separately with the most relevant interpolation order. We show the relevance of such a method on tailored datasets and the benefits and challenges when applying it to large scientific data.

Vidal, Nicolas [ORNL]↗

DOE OSTI2016

Plug-and-Play Priors for Bright Field Electron Tomography and Sparse Interpolation

Many material and biological samples in scientific imaging are characterized by nonlocal repeating structures. These are studied using scanning electron microscopy and electron tomography. Sparse sampling of individual pixels in a two-dimensional image acquisition geometry, or sparse sampling of projection images with large tilt increments in a tomography experiment, can enable high speed data acquisition and minimize sample damage caused by the electron beam. In this paper, we present an algorithm for electron tomographic reconstruction and sparse image interpolation that exploits the nonlocal redundancy in images. We adapt a framework, termed plug-and-play priors, to solve these imaging problems in a regularized inversion setting. The power of the plug-and-play approach is that it allows a wide array of modern denoising algorithms to be used as a “prior model” for tomography and image interpolation. We also present sufficient mathematical conditions that ensure convergence of the plug-and-play approach, and we use these insights to design a new nonlocal means denoising algorithm. Finally, we demonstrate that the algorithm produces higher quality reconstructions on both simulated and real electron microscope data, along with improved convergence properties compared to other methods.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Combined Dimensional and Topology Optimization of Synchronous Machine Rotors Using a Material Density Interpolation Method

This article presents a magneto-structural combined dimensional and topology optimization technique for interior permanent magnet synchronous machine (IPMSM) and wound field synchronous machine (WFSM) rotors. Dimensional changes to the permanent magnet (PM) and rotor winding location or size are accomplished by interpolating or projecting a smoothed Heaviside rectangular function representing the presence of PM material onto the IPMSM rotor design domain mesh, whereas the presence of copper material for WFSM rotor is interpolated in the same manner. A density based Solid Isotropic Material with Penalization (SIMP) topology optimization approach is then used to vary the presence of electrical steel in mesh elements to form flux barriers around the PM. The proposed method enforces a defined shape for the PM without requiring the mesh in the design domain to be deformed. Four examples are presented to demonstrate the technique: Two flat bar IPMSM, one V-shaped IPMSM, and a WFSM. Here, a comparative study is performed on one of the flat bar IPMSMs using Metamodel of Optimal Prognosis (MOP) based method.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Fast and Invertible Simplicial Approximation of Magnetic‐Following Interpolation for Visualizing Fusion Plasma Simulation Data

We introduce a fast and invertible approximation for fusion plasma simulation data represented as 2D planar meshes with connectivities approximating magnetic field lines along the toroidal dimension in deformed 3D toroidal spaces. Scientific variables (e.g., density and temperature) in these fusion data are interpolated following a complex magnetic-field-line-following scheme in the toroidal space represented by a cylindrical coordinate system. This deformation in the 3D space poses challenges for root-finding and interpolation. To this end, we propose a novel paradigm for visualizing and analyzing such data based on a newly developed algorithm for constructing a 3D simplicial mesh within the deformed 3D space. Our algorithm generates a tetrahedral mesh that connects the 2D meshes using tetrahedra while adhering to the constraints on node connectivities imposed by the magnetic field-line scheme. Specifically, we first divide the space into smaller partitions to reduce complexity based on the input geometries and constraints on connectivities. Then, we independently search for a feasible tetrahedralization of each partition, considering nonconvexity. We demonstrate our method with two X-Point Gyrokinetic Code (XGC) simulation datasets on the International Thermonuclear Experimental Reactor (ITER) and Wendelstein 7-X (W7-X), and use an ocean simulation dataset to substantiate broader applicability of our method. An open source implementation of our algorithm is available at https://github.com/rcrcarissa/DeformedSpaceTet.

Ren, Congrong [The Ohio State Univ., Columbus, OH ↗

DOE OSTI2014

SU-D-12A-01: An Inter-Projection Interpolation (IPI) Approach for the Synchronized Moving Grid (SMOG) to Reduce Dose in Cone Beam CT

Purpose: Synchronized moving grid is a promising technique to reduce scatter and ghost artifacts in cone beam computed tomography (CBCT). However, it requires 2 projections in the same gantry angle to obtain full information due to signal blockage by the grid. We proposed an inter-projection interpolation (IPI) method to estimate blocked signals, which may reduce the scan time and the dose. This study aims to provide a framework to achieve a balance between speed, dose and image quality. Methods: The IPI method is based on the hypothesis that an abrupt signal in a projection can be well predicted by the information in the two immediate neighboring projections if the gantry angle step is small. The study was performed on a Catphan and a head phantom. The SMOG was simulated by erasing the information (filling with “0”) of the areas in each projection corresponding to the grid. An IPI algorithm was applied on each projection to recover the erased information. FDK algorithm was used to reconstruct CBCT images for the IPI-processed projections, and compared with the original image in term of signal to noise ratio (SNR) measured in the whole reconstruction image range. The effect of gantry angle step was investigated by comparing the CBCT images from projection sets of various gantry intervals, with IPI-predicted projections to fill the missing projection in the interval. Results: The IPI procession time was 1.79s±0.53s for each projection. SNR after IPI was 29.0db and 28.1db for the Catphan and head phantom, respectively, comparing to 15.3db and 22.7db for an inpainting based interpolation technique. SNR was 28.3, 28.3, 21.8, 19.3 and 17.3 db for gantry angle intervals of 1, 1.5, 2, 2.5 and 3 degrees, respectively. Conclusion: IPI is feasible to estimate the missing information, and achieve an reasonable CBCT image quality with reduced dose and scan time. This study is supported by NIH/NCI grant 1R01CA166948-01.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Interpolating Low Time-Resolution Forecast Data

Methodology that interpolates low time-resolution data (e.g., hourly) to high time-resolution (e.g., minutely) with variability patterns extracted from historical records. Magnitude of the variability inserted into the low time resolution data can be adjusted according to the installed capacity represented by the low time-resolution data compared to that by historical records. This approach enables detailed analysis of the impacts from wind and solar on power system intra-hour operations and balancing reserve requirements even with only hourly data. It also allows convenient creation of high resolution wind or solar generation data with various degree of variability to investigate their operational impacts. The methodology comprises of the following steps: 1. Smooth the historical data (set A) with an appropriate window length l to get its trend (set B); l can be a fraction of an hour (e.g., 15 minutes) or longer than an hour, of which the length of the variability patterns will be; 2. Extract the variable component (set C) of historical data by subtracting the smooth trend from it, i.e. set C = set A - set B 3. For each window length l of the variable component data set, find the average value x (will call it base component) of the corresponding window of the historical data set; 4. Define a series of segments (set D) that the values of data will be grouped into, e.g. (0, 0.1), (0.1, 0.2), ..., (0.9, 1.0) after normalization; Link each variability pattern to a data segment based on its corresponding base component x; after this step, each data segment should be linked to multiple variability patterns after this step; 5. Use spline function to interpolate the low time-resolution forecast data (set E) to become a high time-resolution smooth curve (set F); 6. Based on the window length l , calculate the average value y in each window length of set F; find the data segment that y belongs to; then randomly select one of the variability patterns linked to this data segment; 7. Add the selected variability pattern in each window to set F to get the high time-resolution data (set G) with variability observed from historical data set A; a coefficient k can be multiplied to the selected variability pattern to reflect the potential difference in variability between forecast set E and historical data set A caused by installation capacity and other factors. 8. After set G is constructed, variability statistics can be checked for set A and set G for comparison. Adjust the coefficient k if necessary to achieve desired variability in the reconstructed data.

Lu, Shuai↗

DOE OSTI2022

Comparing machine learning and interpolation methods for loop-level calculations

The need to approximate functions is ubiquitous in science, either due to empirical constraints or high computational cost of accessing the function. In high-energy physics, the precise computation of the scattering cross-section of a process requires the evaluation of computationally intensive integrals. A wide variety of methods in machine learning have been used to tackle this problem, but often the motivation of using one method over another is lacking. Comparing these methods is typically highly dependent on the problem at hand, so we specify to the case where we can evaluate the function a large number of times, after which quick and accurate evaluation can take place. We consider four interpolation and three machine learning techniques and compare their performance on three toy functions, the four-point scalar Passarino-Veltman D_0 D 0 function, and the two-loop self-energy master integral M. We find that in low dimensions (d = 3), traditional interpolation techniques like the Radial Basis Function perform very well, but in higher dimensions (d=5, 6, 9) we find that multi-layer perceptrons (a.k.a neural networks) do not suffer as much from the curse of dimensionality and provide the fastest and most accurate predictions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Interpolated Sounding and Gridded Sounding Value-Added Products

Standard Atmospheric Radiation Measurement (ARM) Climate Research Facility sounding files provide atmospheric state data in one dimension of increasing time and height per sonde launch. Many applications require a quick estimate of the atmospheric state at higher time resolution. The INTERPOLATEDSONDE (i.e., Interpolated Sounding) Value-Added Product (VAP) transforms sounding data into continuous daily files on a fixed time-height grid, at 1-minute time resolution, on 332 levels, from the surface up to a limit of approximately 40 km. The grid extends that high so the full height of soundings can be captured; however, most soundings terminate at an altitude between 25 and 30 km, above which no data is provided. Between soundings, the VAP linearly interpolates atmospheric state variables in time for each height level. In addition, INTERPOLATEDSONDE provides relative humidity scaled to microwave radiometer (MWR) observations.The INTERPOLATEDSONDE VAP, a continuous time-height grid of relative humidity-corrected sounding data, is intended to provide input to higher-order products, such as the Merged Soundings (MERGESONDE; Troyan 2012) VAP, which extends INTERPOLATEDSONDE by incorporating model data. The INTERPOLATEDSONDE VAP also is used to correct gaseous attenuation of radar reflectivity in products such as the KAZRCOR VAP.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Interpolated Sounding and Gridded Sounding Value-Added Products

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Polynomial Similarity Transformation Theory: A smooth interpolation between coupled cluster doubles and projected BCS applied to the reduced BCS Hamiltonian

We present a similarity transformation theory based on a polynomial form of a particle-hole pair excitation operator. In the weakly correlated limit, this polynomial becomes an exponential, leading to coupled cluster doubles. In the opposite strongly correlated limit, the polynomial becomes an extended Bessel expansion and yields the projected BCS wavefunction. In between, we interpolate using a single parameter. The e ective Hamiltonian is non-hermitian and this Polynomial Similarity Transformation Theory follows the philosophy of traditional coupled cluster, left projecting the transformed Hamiltonian onto subspaces of the Hilbert space in which the wave function variance is forced to be zero. Similarly, the interpolation parameter is obtained through minimizing the next residual in the projective hierarchy. We rationalize and demonstrate how and why coupled cluster doubles is ill suited to the strongly correlated limit whereas the Bessel expansion remains well behaved. The model provides accurate wave functions with energy errors that in its best variant are smaller than 1% across all interaction stengths. The numerical cost is polynomial in system size and the theory can be straightforwardly applied to any realistic Hamiltonian.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Grain Boundary Energy Interpolation for Uranium Dioxide

Grain boundary (GB) energies have an effect on the frequency with which certain GBs are present in a material. Since GBs have a direct impact on a material's properties, knowledge of GB energies is useful in determining those properties. This knowledge is difficult to establish due to an incomplete understanding of the atomic structure of a GB. Previous research has successfully created an interpolation function for the GB energies of face-centered cubic (fcc) metals. An extension of this work was recently done by creating an interpolation function for uranium dioxide (UO 2 ), which has a fluorite crystal structure. This work improves the accuracy of that function. Molecular dynamics (MD) simulations were used to calculate more accurate GB energies for UO 2 , and new fitting parameters were calculated from those results. Comparison with previous work shows an overall improvement.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Block smoothers and generalized ideal interpolation in AMG (Final Report)

The Pennsylvania State University (“Subcontractor”) worked on developing new parallel algebraic multilevel methods suitable for solving PDEs. Specifically, work on the design of multigrid solvers for coupled systems of partial differential equations arising in numerical modeling of various applications was completed. A main emphasis was on the design of new ideal algebraic multigrid interpolation for problems such as Maxwell’s equations where block smoothers are needed and the standard form of ideal interpolation is not an effective choice.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Interpolation as a means of shift selection for multilevel Monte Carlo with lattice displacements

The calculation of disconnected diagram contributions to physical signals is a computationally expensive task in Lattice QCD. To extract the physical signal, the trace of the inverse Lattice Dirac operator, a large sparse matrix, must be stochastically estimated. Because the variance of the stochastic estimator is typically large, variance reduction techniques must be employed. Multilevel Monte Carlo (MLMC) methods reduce the variance of the trace estimator by utilizing a telescoping sequence of estimators. Frequency Splitting is one such method that uses a sequence of inverses of shifted operators to estimate the trace of the inverse of the lattice Dirac operator, however there is no a priori way to select the shifts that minimize the cost of the multilevel trace estimation. We present a sampling and interpolation scheme that is able to predict the variances associated with Frequency Splitting under displacements of the underlying space time lattice. The interpolation scheme is able to predict the variances to high accuracy and therefore choose shifts that correspond to an approximate minimum of the cost for the trace estimation. We show that Frequency Splitting with the chosen shifts displays significant speedups over multigrid deflation

Whyte, Travis↗