Search NASA⌕ Search

SEARCH · Search NASA

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Validation of Delamination Growth Predictions under Fatigue Loading using the Single Leg Bending Specimen

The Single Leg Bending Specimen (SLB) was chosen for a study to compare experimentally determined delamination growth with analysis predictions. First, fatigue tests were performed on SLB specimens with several different starter delamination lengths. Specimens were cycled in load control and delamination length increase was continuously recorded with increasing number of cycles. Test results from different load levels were combined to capture the entire delamination growth law (Paris Law) over a wide range of energy release rates. Second, a recently developed interpolation scheme was used to predict the Paris Law associated with the SLB specimen. This interpolation scheme allows the calculation of the growth rates for a given mixed-mode ratio from known single-mode growth laws. Third, finite element analyses were performed using two-dimensional models of the SLB specimen. New Virtual Crack Closure Technique based analysis tools for automated delamination growth analysis under cyclic loading in Abaqus were used for growth predictions. Fourth, the predictions from analyses were compared to experimental fatigue data. Remaining analysis shortcomings were highlighted.

Ronald Krueger↗
NASA NTRS2020

User's Guide for GAA_JET_FV (v1): A Jet Noise Prediction Code Based on the Generalized Acoustic Analogy

This document is a user’s guide for the jet noise prediction code GAA_JET_FV, which can be used to make predictions of turbulent mixing noise in high-speed free jets (ie. in the absence of any solid surfaces) of arbitrary cross section. The code requires a Reynolds-averaged Navier-Stokes (RANS) solution for the mean flow and turbulence as input. A script is provided in the code package which can be used to interpolate structured or unstructured RANS solutions onto a structured grid suitable for the noise calculations. Output file formats for two commonly used RANS solvers are currently supported by this script. The document describes how the code can be obtained and installed on a user’s system. A simple test case is provided that can be run with minimal user knowledge of the code details. General instructions for running the interpolation script and the main code are given along with descriptions of the input and output data files and individual code modules. Several additional test cases are provided which allow the user to exercise additional features of the code. This document is Version 1, Revision 0 of the User’s Guide, which contains examples of round and non-axisymmetric unheated jet test cases. Future versions are planned which will include additional functionality for the code and more complex test cases.

Jet Noise↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRS2021

Waypoint Following Dynamics of a Quaternion Error Feedback Attitude Control System

Closed-loop attitude steering can be used to implement alternative attitude maneuvers by using a conventional attitude control system to track a non-standard attitude profile sampled as a series of sequential commands. The idea has been employed in practice to perform zero-propellant maneuvers on the International Space Station and minimum time maneuvers on NASA’s TRACE space telescope. A challenge for operational implementation of the idea is the finite capacity of a space vehicle’s command storage buffer. One approach to mitigate the problem is to downsample-and-hold the attitude commands so that the attitude control system to follows a set of waypoints. This paper explores the waypoint following dynamics of a quaternion error feedback control law for downsample-and-hold. It is shown that downsample-and-hold induces a ripple between downsamples that causes the satellite angular rate to significantly overshoot the desired limit. Analysis in the z-domain is carried out in order to understand the phenomenon. An interpolating Chebyshev-type filter is proposed that allows attitude commands to be encoded in terms of a small set of filter coefficients. Using the interpolating filter, commands can be issued at the ACS rate but with significantly reduced memory requirements. The attitude control system of NASA’s Lunar Reconnaissance Orbiter is used as an example to illustrate the behavior of a practical attitude control system.

Mark Karpenko↗
NASA NTRS2020

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗
NASA NTRSDate not supplied

Angle-of-Attack Sweep with Mesh Adaptation for High-Lift Configurations

This paper reports an angle-of-attack (alpha) sweep study with mesh adaptation for the high-lift configuration - JAXA Standard Model (JSM) from the 3rd AIAA High-Lift Prediction Workshop. The method of angle-of-attack (alpha) continuation, where the flow-field at a certain angle of attack is initialized from the solution obtained at the previous angle of attack, is applied to adapted unstructured meshes. Since the meshes are continuously adapted, flow-field initialization on a particular mesh for a certain angle is done through linear interpolation of the solution from the previous mesh. Interpolated initial conditions usually prevent massively separated solutions that may otherwise occur with impulsive initial conditions. Alpha continuations are done in both forward (increasing alpha) and reverse (decreasing alpha) directions, resulting in the formation of hysteresis loops. Skin frictions are compared with the wind-tunnel oil flow images and pressure coefficients are compared with the wind-tunnel values, at selected angles of attack. Final paper will include results also for the high-lift version of the Common Research Model (CRM-HL) from the 4th AIAA High-Lift Prediction Workshop. A thorough comparison of the results will also be made with the available wind-tunnel data.

Mesh Adaptation↗
NASA NTRS2021

Probabilistic Machine Learning Estimation of Ocean Mixed Layer Depth from Dense Satellite and Sparse In-Situ Observations

The ocean mixed layer plays an important role in the coupling between the upper ocean and atmosphere across a wide range of time scales. Estimation of the variability of the ocean mixed layer is therefore important for atmosphere-ocean prediction and analysis. The increasing coverage of in situ Argo profile data allows for an increasingly accurate analysis of the mixed layer depth (MLD) variability associated with deviations from the seasonal climatology. However, sampling rates are not sufficient to fully resolve subseasonal (<90 day) MLD variability. Yet, many multivariate observations-based analyses include implicit modeled subseasonal MLD variability. One analysis method is optimal interpolation of in situ data, but the interior analysis can be improved by leveraging surface data with regression or variational approaches. Here, we demonstrate how machine learning methods and satellite sea surface temperature, salinity, and height facilitate MLD estimation in a pilot study of two regions: the mid-latitude southern Indian and the eastern equatorial Pacific Oceans. We construct multiple machine learning architectures to produce weekly 1/2° gridded MLD anomaly fields (relative to a monthly climatology) with uncertainty estimates. We test multiple traditional and probabilistic machine learning techniques to compare both accuracy and probabilistic calibration. We validate our methodology by applying it to ocean model simulations. We find that incorporating sea surface data through a machine learning model improves the performance of spatiotemporal MLD variability estimation compared to optimal interpolation of Argo observations alone. These preliminary results are a promising first step for the application of machine learning to MLD prediction.

Machine Learning↗
NASA NTRSDate not supplied

Cross-Validation of Computational and Experimental Distributed Surface Pressures on the Space Launch System

This paper presents a new workflow for comparing experimental pressure-sensitive paint (PSP) data to computational fluid dynamic (CFD) simulations by way of mapping data from corresponding grids utilizing interpolation methods. In addition to generating quantitative and qualitative point-to-point comparisons between PSP and CFD data, this workflow extracts sectional loading data from both grids and generates lineload comparison charts for corresponding PSP and CFD runs. Experimental PSP data presented in this paper were taken from a 2016 NASA Ames Research Center Unitary Plan Wind Tunnel 11- by 11-Foot Transonic WindTunnel Facility test of the NASA Space Launch System. CFD simulation data for comparison purposes were generated using the FUN3D code. Overall, interpolation onto PSP grids versus CFD grids yields comparable surface pressure fields. However, lineload comparisons are easier to make on the CFD grid-mapped data due to the grid topology and the current capabilities of the lineload analysis tools at NASA Langley Research Center. This workflow is written using contemporary software (Python, Tecplot, PyTecplot), is compatible with existing tools at NASA Langley, and is developed to be adaptable depending on the situation.

SLS↗
NASA NTRSDate not supplied

Secondary Sonic Boom Ray Path Prediction using sBOOM

Sonic boom ray path equations are solved numerically to show the extent of secondary booms compared to primary sonic boom carpet. Secondary booms are shown to reach altitudes of ~125 km before turning toward the ground and impinging at lateral distances of ~50 km and beyond. Exact closed-form solutions of the ray paths and their partial derivatives with respect to certain relevant parameters are derived for the simplified case of windless atmospheres and compared against numerical solutions. Ray tube areas are computed using the Jacobian of a single ray as opposed to the traditional technique of using 4 rays that are temporally and azimuthally perturbed. The ray tube areas corresponding to secondary booms can get 1000 times larger than those from primary booms. The numerically computed ray tube areas are verified using exact expressions for the special case of windless atmospheres. Finally, pressure scaling due to ray tube area changes and atmospheric stratification is presented and the shortcomings associated with the current atmospheric interpolation techniques are discussed. It is determined that piece-wise linear interpolation of the atmospheric temperature may be inadequate and could lead to numerical noise in the prediction of secondary booms.

Sriram K Rallabhandi↗
NASA NTRS2021

Weekly Mapping of Sea Ice Freeboard in the Ross Sea from ICESat-2

NASA’s ICESat-2 has been providing sea ice freeboard measurements across the polar regions since October 2018. In spite of the outstanding spatial resolution and precision of ICESat-2, the spatial sparsity of the data can be a critical issue for sea ice monitoring. This study employs a geostatistical approach (i.e., ordinary kriging) to characterize the spatial autocorrelation of the ICESat-2 freeboard measurements (ATL10) to estimate weekly freeboard variations in 2019 for the entire Ross Sea area, including where ICESat-2 tracks are not directly available. Three variogram models (exponential, Gaussian, and spherical) are compared in this study. According to the cross-validation results, the kriging-estimated freeboards show correlation coefficients of 0.56–0.57, root mean square error (RMSE) of ~0.12 m, and mean absolute error (MAE) of ~0.07 m with the actual ATL10 freeboard measurements. In addition, the estimated errors of the kriging interpolation are low in autumn and high in winter to spring, and low in southern regions and high in northern regions of the Ross Sea. The effective ranges of the variograms are 5–10 km and the results from the three variogram models do not show significant differences with each other. The southwest (SW) sector of the Ross Sea shows low and consistent freeboard over the entire year because of the frequent opening of wide polynya areas generating new ice in this sector. However, the southeast (SE) sector shows large variations in freeboard, which demonstrates the advection of thick multiyear ice from the Amundsen Sea into the Ross Sea. Thus, this kriging-based interpolation of ICESat-2 freeboard can be used in the future to estimate accurate sea ice production over the Ross Sea by incorporating other remote sensing data.

Satellite altimeter↗
NASA NTRSDate not supplied

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗
NASA NTRS2023

Kamodo: Simplifying Model Data Access and Utilization

To address the lack of user-friendly software needed to simplify the utilization of model data across Heliophysics, the Community Coordinated Modeling Center (CCMC) at NASA’s Goddard Space Flight Center has developed a model-agnostic method via Kamodo for users to easily access and utilize model data in their workflows. By abstracting away the broad range of file formats and the intricacies of interpolation on specialized grids, this approach significantly lowers the barrier to model data access and utilization for the community while adding exciting new capabilities to their tool boxes. This paper describes the direct interfaces to the model data, called model readers, and a basic introduction on how to use them. Additionally, we detail the planned approach for including custom interpolation codes, and include current progress on specialized visualization developments. The CCMC is maintaining Kamodo as an official NASA open-sourced software to enable and encourage community collaboration.

Heliophysics↗
NASA NTRSDate not supplied

Assessing the Uncertainty Impact of CERES Fast Longwave and SHortwave Radiative Flux (FLASHFlux) Level 3 Product With and Without Terra Observations

The Clouds and Earth’s Radiant Energy System (CERES) project provides satellite-based observations of how Earth’s energy flows are varying in time and space and how clouds and aerosols are affecting the Earth’s radiation budget. Nominally, CERES data products require months of validation and calibration before releasing a climate quality data. The Fast Longwave And SHortwave radiative Flux (FLASHFlux) data product was developed to provide data for applied science research involving the renewable energy and agricultural sectors within a week of observation. FLASHFlux achieves this by using simplified calibration, an operational meteorological product from Global Monitoring and Assimilation Office (GMAO), and a surface parameterizations model. The CERES FLASHFlux provides two data products: 1) an hourly Level 2 Single Scanner Footprint (SSF) data separately for Terra and NOAA-20 observations. 2) a daily Level 3 Time Interpolated and Spatially Averaged (TISA) gridded data that combines NOAA-20 and Terra observations on a one-degree equal angle grid. FLASHFlux TISA data product interpolate on a diurnal model that assumes a satellite equilateral crossing time of 10:30 AM and 1:30 PM from Terra and Aqua, respectively. Aqua was replaced by NOAA-20 starting on September 2022. Terra is planned to be replace by the Satellite ClOud and Radiation Property retrieval System (SatCORPS) soon. We are currently using the Terra observations as it continues to drift. We assess the impact of FLASHFlux TISA data when Terra is removed. An uncertainty estimate of the Top-Of-Atmosphere (TOA) fluxes are given of FLASHFlux Version4A (before Terra drift) and Version4B (current), and Version4B (no Terra) in comparison to the CERES EBAF and SYN1deg. In addition, we compare FLASHFlux Version4B and Version4B (no Terra) surface radiative fluxes to ground base measurements to determine the impact of running without Terra data.

PC Sawaengphokhai↗
NASA NTRS2024

Validation and Application of Satellite-Derived Sea Surface Temperature Gradients in the Bering Strait and Bering Sea

The Arctic is one of the most important regions in the world’s oceans for understanding the impacts of a changing climate. Yet, it is also difficult to measure because of extreme weather and ice conditions. In this work, we directly compare four datasets from the Group for High-Resolution Sea Surface Temperature (GHRSST) with a NASA Saildrone deployment along the Alaskan Coast and the Bering Sea and Bering Strait. The four datasets used are the Remote Sensing Systems Microwave Infrared Optimally Interpolated (MWIR) product, the Canadian Meteorological Center (CMC) product, the Daily Optimally Interpolated Product (DOISST), and the Operational Sea Surface Temperature and Ice Analysis (OSTIA) product. Spatial sea surface temperature (SST) gradients were derived for both the Saildrone deployment and GHRSST products, with the GHRSST products collocated with the Saildrone deployment. Overall, statistics indicate that the OSTIA product had a correlation of 0.79 and a root mean square difference of 0.11 °C/km when compared with Saildrone. CMC had the highest correlation of 0.81. Scatter plots indicate that OSTIA had the slope closest to one, thus best reproducing the magnitudes of the Saildrone gradients. Differences increased at latitudes > 65°N where sea ice would have a greater impact. A trend analysis was then performed on the gradient fields. Overall, positive trends in gradients occurred in areas along the coastal regions. A negative trend occurred at approximately 60°N. A major finding of this study is that future work needs to revolve around the impact of changing ice conditions on SST gradients. Another major finding is that a northward shift in the southern ice edge occurred after 2010 with a maxima at approximately 2019. This indicates that the shift of the southern ice edge is not gradual but has dramatically increased over the last decade. Future work needs to revolve around examining the possible causes for this northward shift.

Arctic↗
NASA NTRSDate not supplied

Derivation of Integrated Load Distributions from Resampled Computational Data

Principal component analysis (PCA) has been the center of many surrogate models used to characterize fluid flows in recent years. However, little work has been done to character- ize the uncertainty in the PCA transform itself and its effect on derived surrogate models. To explore the uncertainty, a typical interpolated surrogate model is constructed for a represen- tative aerodynamic body from computational data. The computational data is then resampled to explore the robustness of the PCA transformation. The variations of the model predictions during this resampling are analyzed to get a measure of confidence in the PCA transformation, which is then applied to the interpolated surrogate model to get uncertainty on integrated force predictions. An initial test case has been explored with promising results.

Principal Component Analysis↗
NASA NTRSDate not supplied

A Queuing Theory Approach to Pilot-Controller Coordination for m:N Operations

In recent years, attention and interest by industry and researchers has grown in a control paradigm for remotely piloted aircraft termed “m:N operations.” In an m:N operation, a team of m remote pilots in command (RIPCs) collaboratively manage the flights of N aircraft. A consequence of an m:N concept of operations is that the RPICs will have to switch attention from one aircraft to another and from one task to another. Previous research in m:N operations has focused on the workload experienced by an RPIC and their level of situation awareness on their flights. Researchers have found that RPIC workload and situation awareness are generally sensitive to increasing N, although NASA’s Multi-Vehicle (m:N) Working Group has suggested that the driver of workload/situation awareness is the number of exceptions requiring human intervention as opposed to the value of N itself. In any case, a natural antecedent of workload is task load. In this paper, queueing theory is applied to a 1:N Urban Air Mobility (UAM) air taxi operation in order to estimate pilot task load for managing radio communications with air traffic controllers (ATCs) under increasing N. An M/M/1 queueing system is used to model the RIPC’s servicing of calls and clearance requests (e.g., departure, arrival, or airspace transition) to ATC for the N aircraft. Important parameters for the queueing model are the task arrival rate and the average service time for task completion. Radio communication times from past human-in-the-loop simulation studies are used to measure service times for a 1:4 and 1:12 UAM operation and to interpolate service times for 4 < N < 12. A Monte Carlo method is then employed, using the measured and interpolated service times, to estimate arrival rate and related queueing statistics. The paper concludes by considering the estimated queuing statistics, particularly the RPIC’s utilization (i.e., proportion of time actively servicing tasks), the length of the task queue over time, and the implications for task-balanced system design.

task load↗
NASA NTRSDate not supplied

Derivation of Integrated Load Distributions from Resampled Computational Data

Principal component analysis (PCA) has been the center of many surrogate models used to characterize fluid flows in recent years. However, little work has been done to character- ize the uncertainty in the PCA transform itself and its effect on derived surrogate models. To explore the uncertainty, a typical interpolated surrogate model is constructed for a represen- tative aerodynamic body from computational data. The computational data is then resampled to explore the robustness of the PCA transformation. The variations of the model predictions during this resampling are analyzed to get a measure of confidence in the PCA transformation, which is then applied to the interpolated surrogate model to get uncertainty on integrated force predictions. An initial test case has been explored with promising results.

Computational Fluid Dynamics↗