Automatic selection of interpolation segments for a machine tool with linear interpolator
Automatic selection of interpolation segments for machine tool with linear interpolator
MACHINE TOOL↗
SEARCH · Search NASA
Results for “linear 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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Random sampling of random processes - Optimum linear interpolation.
Optimum mean square linear interpolation for random sampling, discussing periodic and Poisson sampling
RANDOM SAMPLE↗
On the use of linear interpolation in opacity tables
Errors resulting from linear interpolation use in opacity tables for stellar interior calculations
Rood, R. T.↗
Data smoothing and error detection based on linear interpolation
A method is presented, based on linear interpolation, for detecting and correcting bad data points in a set of data without contaminating the good data points. The method used is not concerned with the small random errors usually attributed to a noisy system. It assumes that the data points which are in error are relatively isolated from each other and that the number of such points is small compared to the total number of data points.
Guerra, V. M.↗
INTERP3: A computer routine for linear interpolation of trivariate functions defined by nondistinct unequally spaced variables
A report on the computer routine INTERP3 is presented. The routine is designed to linearly interpolate a variable which is a function of three independent variables. The variables within the parameter arrays do not have to be distinct, or equally spaced, and the array variables can be in increasing or decreasing order.
Hill, D. C.↗
Linear interpolation of four-dimensional tabulated data for computers with single subscripted variable capability
Using only a one-dimensional subscripted variable, a FORTRAN computer subprogram was developed to linearly interpolate tabulated data of functions of four or less variables. The primary motivation was for faster computation.
Farr, W. R.↗
Differences between measured and linearly interpolated synoptic variables over a 12-h period during AVE 4
Results of analyses revealed that nonlinear changes or differences formed centers or systems, that were mesosynoptic in nature. These systems correlated well in space with upper level short waves, frontal zones, and radar observed convection, and were very systematic in time and space. Many of the centers of differences were well established in the vertical, extending up to the tropopause. Statistical analysis showed that on the average nonlinear changes were larger in convective areas than nonconvective regions. Errors often exceeding 100 percent were made by assuming variables to change linearly through a 12-h period in areas of thunderstorms, indicating that these nonlinear changes are important in the development of severe weather. Linear changes, however, accounted for more and more of an observed change as the time interval (within the 12-h interpolation period) increased, implying that the accuracy of linear interpolation increased over larger time intervals.
Dupuis, L. R.↗
Deriving Three Dimensional Reservoir Bathymetry from Multi-Satellite Datasets
We evaluate different techniques that rebuild reservoir bathymetry by combining multi-satellite imagery of surface water elevation and extent. Digital elevation models (DEMs) are processed in two distinct ways in order to determine 3-D reservoir bathymetry. They are defined as (a) linear extrapolation and (b) linear interpolation. The first one linearly extrapolates the land slope, defining the bottom as the intersection of all extrapolated lines. The second linearly interpolates the uppermost and lowermost pixels of the reservoir's main river, repeating the process for all other tributaries. A visible bathymetry, resulting from the combination of radar altimetry and water extent masks, can be coupled with the DEM, improving the accuracy of techniques (a) and (b). Envisat-and Altika-based altimetric time series is combined to a Landsat-based water extent database over the2002-2016 period in order to generate the visible bathymetry, and topography is derived from the 3-arcsecHydroSHEDS DEM. Fourteen 3-D bathymetries derived from the combination of these techniques and datasets, plus the inclusion of upstream and downstream riverbed elevations, are evaluated over Lake Mead. Accuracy is measured using ground observations, and show that metrics improve as a function of added data requirement and processing. Best bathymetry estimates are obtained when the visible bathymetry, linear extrapolation technique and riverbed elevation are combined. Water storage variability is also evaluated and shows that best results are derived from the aforementioned combination. This study contributes to our understanding and representation of reservoir water impoundment impacts on the hydrological cycle.
HydroSHEDS↗
Noise and drift analysis of non-equally spaced timing data
Generally, it is possible to obtain equally spaced timing data from oscillators. The measurement of the drifts and noises affecting oscillators is then performed by using a variance (Allan variance, modified Allan variance, or time variance) or a system of several variances (multivariance method). However, in some cases, several samples, or even several sets of samples, are missing. In the case of millisecond pulsar timing data, for instance, observations are quite irregularly spaced in time. Nevertheless, since some observations are very close together (one minute) and since the timing data sequence is very long (more than ten years), information on both short-term and long-term stability is available. Unfortunately, a direct variance analysis is not possible without interpolating missing data. Different interpolation algorithms (linear interpolation, cubic spline) are used to calculate variances in order to verify that they neither lose information nor add erroneous information. A comparison of the results of the different algorithms is given. Finally, the multivariance method was adapted to the measurement sequence of the millisecond pulsar timing data: the responses of each variance of the system are calculated for each type of noise and drift, with the same missing samples as in the pulsar timing sequence. An estimation of precision, dynamics, and separability of this method is given.
Vernotte, F.↗
Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations
This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.
97 MATHEMATICS AND COMPUTING↗
Filling the Gaps: A Bayesian Mixture Model for Imputing Missing Soil Water Content Data
ABSTRACT Soil water content (SWC) data are central to evaluating how soil moisture varies over time and space and influences critical plant and ecosystem functions, especially in water‐limited drylands. However, sensors that record SWC at high frequencies often malfunction, leading to incomplete timeseries and limiting our understanding of dryland ecosystem dynamics. We developed an analytical approach to impute missing SWC data, which we tested at six eddy flux tower sites along an elevation gradient in the southwestern United States. We impute missing data as a mixture of linearly interpolated SWC between the observed endpoints of a missing data gap and SWC simulated by an ecosystem water balance model (SOILWAT2). Within a Bayesian framework, we allowed the relative utility (mixture weight) of each component (linearly interpolated vs. SOILWAT2) to vary by depth, site and gap characteristics. We explored “fixed” weights versus “dynamic” weights that vary as a function of cumulative precipitation, average temperature, and time since the start of the gap. Both models estimated missing SWC data well ( R 2 = 0.70–0.88 vs. 0.75–0.91 for fixed vs. dynamic weights, respectively), but the utility of linearly interpolated versus SOILWAT2 values depended on site and depth. SOILWAT2 was more useful for more arid sites, shallower depths, longer and warmer gaps and gaps that received greater precipitation. Overall, the mixture model reliably gap‐fills SWC, while lending insight into processes governing SWC dynamics. This approach to impute missing data could be adapted to accommodate more than two mixture components and other types of environmental timeseries.
Ogle, Kiona [School of Informatics, Computing, and↗
Three dimensional unstructured multigrid for the Euler equations
The three dimensional Euler equations are solved on unstructured tetrahedral meshes using a multigrid strategy. The driving algorithm consists of an explicit vertex-based finite element scheme, which employs an edge-based data structure to assemble the residuals. The multigrid approach employs a sequence of independently generated coarse and fine meshes to accelerate the convergence to steady-state of the fine grid solution. Variables, residuals and corrections are passed back and forth between the various grids of the sequence using linear interpolation. The addresses and weights for interpolation are determined in a preprocessing stage using linear interpolation. The addresses and weights for interpolation are determined in a preprocessing stage using an efficient graph traversal algorithm. The preprocessing operation is shown to require a negligible fraction of the CPU time required by the overall solution procedure, while gains in overall solution efficiencies greater than an order of magnitude are demonstrated on meshes containing up to 350,000 vertices. Solutions using globally regenerated fine meshes as well as adaptively refined meshes are given.
Mavriplis, D. J.↗
Geospatial modeling of near subsurface temperatures of the contiguous United States for assessment of materials degradation
Abstract Understanding subsurface temperature variations is crucial for assessing material degradation in underground structures. This study maps subsurface temperatures across the contiguous United States for depths from 50 to 3500 m, comparing linear interpolation, gradient boosting (LightGBM), neural networks, and a novel hybrid approach combining linear interpolation with LightGBM. Results reveal heterogeneous temperature patterns both horizontally and vertically. The hybrid model performed best achieving a root mean square error of 2.61 °C at shallow depths (50–350 m). Model performance generally decreased with depth, highlighting challenges in deep temperature prediction. State-level analyses emphasized the importance of considering local geological factors. This study provides valuable insights for designing efficient underground facilities and infrastructure, underscoring the need for depth-specific and region-specific modeling approaches in subsurface temperature assessment.
Science & Technology - Other Topics↗
An efficient algorithm for choosing scattering directions in Monte Carlo work with arbitrary phase functions
This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(exp -5).
Barkstrom, Bruce R.↗
Calculating the grain boundary inclination of voxelated grain structures using a smoothing algorithm
We have developed a flexible method for calculating the grain boundary (GB) inclinations of voxelated grain structure data using smoothing algorithms. We compared the performance of four algorithms: the linear interpolation, Allen–Cahn, level-set, and vertex algorithms. We assessed their accuracy using 2D and 3D cases with known inclinations. The vertex algorithm provided the best balance between accuracy and efficiency for 2D structures while the linear interpolation algorithm provided the best balance for 3D structures. We compared the GB inclinations calculated using our smoothing method on a 3D high energy X-ray diffraction microscopy (HEDM) dataset to those determined by meshing the GBs. The two approaches determined similar GB plane distributions, though they varied significantly at triple junctions. In conclusion, the smoothing method was demonstrated for two sources of 3D voxelated grain structures: HEDM data and results from Monte Carlo Potts grain growth simulations.
36 MATERIALS SCIENCE↗
A scalable framework for efficient coupling of thermal and microstructural simulations in additive manufacturing
Predicting microstructure evolution in metal additive manufacturing (AM) is important for process optimization, but spatiotemporal scale disparities between thermal transport and microstructure evolution create significant challenges for efficient data transfer between simulation codes. To address this, we present Stork, a scalable framework for coupling thermal and microstructural simulations. Stork uses a sparse data representation to identify and store active solidification sub-volumes, enabling highly parallel quad-linear interpolation from coarse thermal grids to fine microstructure grids without large intermediate storage. We demonstrate the framework by coupling the semi-analytic heat transfer code 3DThesis with the time-parallel cellular automata code Toucan. This approach achieves over two orders of magnitude reduction in data generation time and file size compared to prior workflows. Numerical studies show that quad-linear interpolation preserves grain morphology and crystallographic texture in laser powder bed fusion (LPBF) simulations for coarsening ratios up to 16. Overall, Stork provides a scalable pathway for high-throughput, component-scale AM simulations on modern high-performance computing systems.
36 MATERIALS SCIENCE↗
Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring
Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.
Salish Sea↗
Characteristics of Whipple Shield Performance in the Shatter Regime
Between the onset of projectile fragmentation and the assumption of rear wall failure due to an impulsive load, multi-wall ballistic limit equations are linearly interpolated to provide reasonable yet conservative predictions of perforation thresholds with conveniently simple mathematics. Although low velocity and hypervelocity regime predictions are based on analytical expressions, there is no such scientific foundation for predictions in the intermediate (or shatter) regime. As the debris flux in low earth orbit (LEO) becomes increasingly dominated by manmade pollution, the profile of micrometeoroid and orbital debris (MMOD) risk shifts continually towards lower velocities. For the International Space Station (ISS), encounter velocities below 7 km/s now constitute approximately 50% of the penetration risk. Considering that the transition velocity from shatter to hypervelocity impact regimes described by common ballistic limit equations (e.g. new non-optimum Whipple shield equation [1]) occurs at 7 km/s, 50% of station risk is now calculated based on failure limit equations with little analytical foundation. To investigate projectile and shield behavior for impact conditions leading to projectile fragmentation and melt, a series of hypervelocity impact tests have been performed on aluminum Whipple shields. In the experiments projectile diameter, bumper thickness, and shield spacing were kept constant, while rear wall thickness was adjusted to determine spallation and perforation limits at various impact velocities and angles. The results, shown in Figure 1 for normal and 45 impacts, demonstrated behavior that was not sufficiently described by the simplified linear interpolation of the NNO equation (also shown in Figure 1). Hopkins et al. [2] investigated the performance of a nominally-identical aluminum Whipple shield, identifying the effects of phase change in the shatter regime. The results (conceptually represented in Figure 2) were found to agree well with those obtained in this study at normal incidence, suggesting that shielding performance in the shatter regime could be well described by considering more complex phase conditions than currently implemented in most BLEs. Furthermore, evidence of these phase effects were found in the oblique test results, providing the basis for an empirical description of these effects that can be applied in MMOD risk assessment software. In this paper, results of the impact experiments are presented, and characteristics of target damage are evaluated. A comparison of intermediate velocity impact failure mechanisms in current BLEs are discussed and compared to the findings of the experimental study. Risk assessment calculations have been made on a simplified structure using currently implemented penetration equations and predicted limits from the experimental program, and the variation in perceived mission risk is discussed. It was found that ballistic limit curves that explicitly incorporated phase change effects within the intermediate regime lead to a decrease in predicted MMOD risk for ISS-representative orbits. When considered for all Whipple-based shielding configurations onboard the ISS, intermediate phase change effects could lead to significant variations in predicted mission risk.
Ryan, Shannon↗