Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis functions”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Further Development of Rotating Rake Mode Measurement Data Analysis

The Rotating Rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. After analysis of the measured data, the mode amplitudes and phases were quantified. For low-speed fans within axisymmetric ducts, mode power levels computed from rotating rake measured data would agree with the far-field power levels on a tone by tone basis. However, this agreement required that the sound from the noise sources within the duct propagated outward from the duct exit without reflection at the exit and previous studies suggested conditions could exist where significant reflections could occur. To directly measure the modes propagating in both directions within a duct, a second rake was mounted to the rotating system with an offset in both the axial and the azimuthal directions. The rotating rake data analysis technique was extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode levels at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode amplitudes for the modes propagating in both directions within the duct. The fit equations were also modified to allow evanescent mode amplitudes to be computed. This extension of the rotating rake data analysis technique was tested using simulated data, numerical code produced data, and preliminary in-duct measured data.

Dahl, Milo D.↗
NASA NTRS2013

Advanced Optical Diagnostics for Ice Crystal Cloud Measurements in the NASA Glenn Propulsion Systems Laboratory

A light extinction tomography technique has been developed to monitor ice water clouds upstream of a direct connected engine in the Propulsion Systems Laboratory (PSL) at NASA Glenn Research Center (GRC). The system consists of 60 laser diodes with sheet generating optics and 120 detectors mounted around a 36-inch diameter ring. The sources are pulsed sequentially while the detectors acquire line-of-sight extinction data for each laser pulse. Using computed tomography algorithms, the extinction data are analyzed to produce a plot of the relative water content in the measurement plane. To target the low-spatial-frequency nature of ice water clouds, unique tomography algorithms were developed using filtered back-projection methods and direct inversion methods that use Gaussian basis functions. With the availability of a priori knowledge of the mean droplet size and the total water content at some point in the measurement plane, the tomography system can provide near real-time in-situ quantitative full-field total water content data at a measurement plane approximately 5 feet upstream of the engine inlet. Results from ice crystal clouds in the PSL are presented. In addition to the optical tomography technique, laser sheet imaging has also been applied in the PSL to provide planar ice cloud uniformity and relative water content data during facility calibration before the tomography system was available and also as validation data for the tomography system. A comparison between the laser sheet system and light extinction tomography resulting data are also presented. Very good agreement of imaged intensity and water content is demonstrated for both techniques. Also, comparative studies between the two techniques show excellent agreement in calculation of bulk total water content averaged over the center of the pipe.

Bencic, Timothy J.↗
NASA NTRS2013

Implementation of a Wavefront-Sensing Algorithm

A computer program has been written as a unique implementation of an image-based wavefront-sensing algorithm reported in "Iterative-Transform Phase Retrieval Using Adaptive Diversity" (GSC-14879-1), NASA Tech Briefs, Vol. 31, No. 4 (April 2007), page 32. This software was originally intended for application to the James Webb Space Telescope, but is also applicable to other segmented-mirror telescopes. The software is capable of determining optical-wavefront information using, as input, a variable number of irradiance measurements collected in defocus planes about the best focal position. The software also uses input of the geometrical definition of the telescope exit pupil (otherwise denoted the pupil mask) to identify the locations of the segments of the primary telescope mirror. From the irradiance data and mask information, the software calculates an estimate of the optical wavefront (a measure of performance) of the telescope generally and across each primary mirror segment specifically. The software is capable of generating irradiance data, wavefront estimates, and basis functions for the full telescope and for each primary-mirror segment. Optionally, each of these pieces of information can be measured or computed outside of the software and incorporated during execution of the software.

Smith, Jeffrey S.↗
NASA NTRS2012

Magnetization of the Lunar Crust

Magnetic fields measured by the satellite Lunar Prospector show large scale features resulting from remanently magnetized crust. Vector data synthesized at satellite altitude from a spherical harmonic model of the lunar crustal field, and the radial component of the magnetometer data, have been used to produce spatially continuous global magnetization models for the lunar crust. The magnetization is expressed in terms of localized basis functions, with a magnetization solution selected having the smallest root-mean square magnetization for a given fit to the data, controlled by a damping parameter. Suites of magnetization models for layers with thicknesses between 10 and 50 km are able to reproduce much of the input data, with global misfits of less than 0.5 nT (within the uncertainties of the data), and some surface field estimates. The magnetization distributions show robust magnitudes for a range of model thicknesses and damping parameters, however the magnetization direction is unconstrained. These global models suggest that magnetized sources of the lunar crust can be represented by a 30 km thick magnetized layer. Average magnetization values in magnetized regions are 30-40 mA/m, similar to the measured magnetizations of the Apollo samples and significantly weaker than crustal magnetizations for Mars and the Earth. These are the first global magnetization models for the Moon, providing lower bounds on the magnitude of lunar crustal magnetization in the absence of multiple sample returns, and can be used to predict the crustal contribution to the lunar magnetic field at a particular location.

Lunar↗
NASA NTRS2014

Combined Uncertainty and A-Posteriori Error Bound Estimates for CFD Calculations: Theory and Implementation

Simulation codes often utilize finite-dimensional approximation resulting in numerical error. Some examples include, numerical methods utilizing grids and finite-dimensional basis functions, particle methods using a finite number of particles. These same simulation codes also often contain sources of uncertainty, for example, uncertain parameters and fields associated with the imposition of initial and boundary data,uncertain physical model parameters such as chemical reaction rates, mixture model parameters, material property parameters, etc.

CFD↗
NASA NTRS2014

Algorithm for Detection of Ground and Canopy Cover in Micropulse Photon-Counting Lidar Altimeter Data in Preparation for the ICESat-2 Mission

NASA's Ice, Cloud and Land Elevation Satellite-II (ICESat-2) mission is a decadal survey mission (2016 launch). The mission objectives are to measure land ice elevation, sea ice freeboard, and changes in these variables, as well as to collect measurements over vegetation to facilitate canopy height determination. Two innovative components will characterize the ICESat-2 lidar: 1) collection of elevation data by a multibeam system and 2) application of micropulse lidar (photon-counting) technology. A photon-counting altimeter yields clouds of discrete points, resulting from returns of individual photons, and hence new data analysis techniques are required for elevation determination and association of the returned points to reflectors of interest. The objective of this paper is to derive an algorithm that allows detection of ground under dense canopy and identification of ground and canopy levels in simulated ICESat-2 data, based on airborne observations with a Sigma Space micropulse lidar. The mathematical algorithm uses spatial statistical and discrete mathematical concepts, including radial basis functions, density measures, geometrical anisotropy, eigenvectors, and geostatistical classification parameters and hyperparameters. Validation shows that ground and canopy elevation, and hence canopy height, can be expected to be observable with high accuracy by ICESat-2 for all expected beam energies considered for instrument design (93.01%-99.57% correctly selected points for a beam with expected return of 0.93 mean signals per shot (msp), and 72.85%-98.68% for 0.48 msp). The algorithm derived here is generally applicable for elevation determination from photoncounting lidar altimeter data collected over forested areas, land ice, sea ice, and land surfaces, as well as for cloud detection.

Algorithms↗
NASA NTRS2014

Analysis of Dual Rotating Rake Data from the NASA Glenn Advanced Noise Control Fan Duct with Artificial Sources

The Rotating Rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. Initially, the mode amplitudes and phases were quantified from a single rake measurement at one axial location. To directly measure the modes propagating in both directions within a duct, a second rake was mounted to the rotating system with an offset in both the axial and the azimuthal directions. The rotating rake data analysis technique was then extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode levels at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode amplitudes for the modes propagating in both directions within the duct. Validation experiments have been conducted using artificial acoustic sources. Results are shown for the measurement of the standing waves in the duct from sound generated by one and two acoustic sources that are separated into the component modes propagating in both directions within the duct. Measured reflection coefficients from the open end of the duct are compared to analytical predictions.

Noise↗
NASA NTRS2013

Nonlinear Modeling by Assembling Piecewise Linear Models

To preserve nonlinearity of a full order system over a parameters range of interest, we propose a simple modeling approach by assembling a set of piecewise local solutions, including the first-order Taylor series terms expanded about some sampling states. The work by Rewienski and White inspired our use of piecewise linear local solutions. The assembly of these local approximations is accomplished by assigning nonlinear weights, through radial basis functions in this study. The efficacy of the proposed procedure is validated for a two-dimensional airfoil moving at different Mach numbers and pitching motions, under which the flow exhibits prominent nonlinear behaviors. All results confirm that our nonlinear model is accurate and stable for predicting not only aerodynamic forces but also detailed flowfields. Moreover, the model is robustness-accurate for inputs considerably different from the base trajectory in form and magnitude. This modeling preserves nonlinearity of the problems considered in a rather simple and accurate manner.

Yao, Weigang↗
NASA NTRS2016

High Resolution Image Reconstruction from Projection of Low Resolution Images DIffering in Subpixel Shifts

In this paper, we demonstrate a simple algorithm that projects low resolution (LR) images differing in subpixel shifts on a high resolution (HR) also called super resolution (SR) grid. The algorithm is very effective in accuracy as well as time efficiency. A number of spatial interpolation techniques using nearest neighbor, inverse-distance weighted averages, Radial Basis Functions (RBF) etc. used in projection yield comparable results. For best accuracy of reconstructing SR image by a factor of two requires four LR images differing in four independent subpixel shifts. The algorithm has two steps: i) registration of low resolution images and (ii) shifting the low resolution images to align with reference image and projecting them on high resolution grid based on the shifts of each low resolution image using different interpolation techniques. Experiments are conducted by simulating low resolution images by subpixel shifts and subsampling of original high resolution image and the reconstructing the high resolution images from the simulated low resolution images. The results of accuracy of reconstruction are compared by using mean squared error measure between original high resolution image and reconstructed image. The algorithm was tested on remote sensing images and found to outperform previously proposed techniques such as Iterative Back Projection algorithm (IBP), Maximum Likelihood (ML), and Maximum a posterior (MAP) algorithms. The algorithm is robust and is not overly sensitive to the registration inaccuracies.

0↗
NASA NTRS2016

On the Assessment of Acoustic Scattering and Shielding by Time Domain Boundary Integral Equation Solutions

Based on the time domain boundary integral equation formulation of the linear convective wave equation, a computational tool dubbed Time Domain Fast Acoustic Scattering Toolkit (TD-FAST) has recently been under development. The time domain approach has a distinct advantage that the solutions at all frequencies are obtained in a single computation. In this paper, the formulation of the integral equation, as well as its stabilization by the Burton-Miller type reformulation, is extended to cases of a constant mean flow in an arbitrary direction. In addition, a "Source Surface" is also introduced in the formulation that can be employed to encapsulate regions of noise sources and to facilitate coupling with CFD simulations. This is particularly useful for applications where the noise sources are not easily described by analytical source terms. Numerical examples are presented to assess the accuracy of the formulation, including a computation of noise shielding by a thin barrier motivated by recent Historical Baseline F31A31 open rotor noise shielding experiments. Furthermore, spatial resolution requirements of the time domain boundary element method are also assessed using point per wavelength metrics. It is found that, using only constant basis functions and high-order quadrature for surface integration, relative errors of less than 2% may be obtained when the surface spatial resolution is 5 points-per-wavelength (PPW) or 25 points-per-wavelength squared (PPW2).

Hu, Fang Q.↗
NASA NTRS2017

Surface-Height Determination of Crevassed Glaciers-Mathematical Principles of an Autoadaptive Density-Dimension Algorithm and Validation Using ICESat-2 Simulator (SIMPL) Data

Glacial acceleration is a main source of uncertainty in sea-level-change assessment. Measurement of ice-surface heights with a spatial and temporal resolution that not only allows elevation-change calculation, but also captures ice-surface morphology and its changes is required to aid in investigations of the geophysical processes associated with glacial acceleration.The Advanced Topographic Laser Altimeter System aboard NASAs future ICESat-2 Mission (launch 2017) will implement multibeam micropulse photon-counting lidar altimetry aimed at measuring ice-surface heights at 0.7-m along-track spacing. The instrument is designed to resolve spatial and temporal variability of rapidly changing glaciers and ice sheets and the Arctic sea ice. The new technology requires the development of a new mathematical algorithm for the retrieval of height information.We introduce the density-dimension algorithm (DDA) that utilizes the radial basis function to calculate a weighted density as a form of data aggregation in the photon cloud and considers density an additional dimension as an aid in auto-adaptive threshold determination. The auto-adaptive capability of the algorithm is necessary to separate returns from noise and signal photons under changing environmental conditions. The algorithm is evaluated using data collected with an ICESat-2 simulator instrument, the Slope Imaging Multi-polarization Photon-counting Lidar, over the heavily crevassed Giesecke Braer in Northwestern Greenland in summer 2015. Results demonstrate that ICESat-2 may be expected to provide ice-surface height measurements over crevassed glaciers and other complex ice surfaces. The DDA is generally applicable for the analysis of airborne and spaceborne micropulse photon-counting lidar data over complex and simple surfaces.

Herzfeld, Ute C.↗
NASA NTRS2016

Tensor-Product Preconditioners for Higher-Order Space-Time Discontinuous Galerkin Methods

space-time discontinuous-Galerkin spectral-element discretization is presented for direct numerical simulation of the compressible Navier-Stokes equat ions. An efficient solution technique based on a matrix-free Newton-Krylov method is developed in order to overcome the stiffness associated with high solution order. The use of tensor-product basis functions is key to maintaining efficiency at high order. Efficient preconditioning methods are presented which can take advantage of the tensor-product formulation. A diagonalized Alternating-Direction-Implicit (ADI) scheme is extended to the space-time discontinuous Galerkin discretization. A new preconditioner for the compressible Euler/Navier-Stokes equations based on the fast-diagonalization method is also presented. Numerical results demonstrate the effectiveness of these preconditioners for the direct numerical simulation of subsonic turbulent flows.

Diosady, Laslo T.↗
NASA NTRS2018

Dual Rotating Rake Measurements of Higher-Order Duct Modes: Validation Using Experimental and Numerical Data

A rotating rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. After analysis of the measured data, the mode coefficient amplitudes and phases were quantified. Early studies using this system found that mode power levels computed from rotating rake measured data would agree with the far-field power levels. However, this agreement required that the sound from the noise sources within the duct propagated outward from the duct exit without reflection and previous studies suggested conditions could exist where significant reflections could occur. This paper shows that mounting a second rake to the rotating system, with an offset in both the axial and the azimuthal directions, measures the data necessary to determine the modes propagating in both directions within a duct. The rotating rake data analysis technique was extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode coefficients at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode coefficients for the modes propagating in both directions within the duct while accounting for the presence of evanescent modes. The validation of the dual-rotating-rake measurements was conducted using data from a combination of experiments and numerical calculations to compute reflection coefficients and other mode coefficient ratios. Compared to results from analytical and numerical computations, the results from dual-rotating-rake measured data followed the expected trends when frequency, mode number, and duct termination geometry were changed.

Validation↗
NASA NTRS2018

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗
NASA NTRS2018

Adjustments and Uncertainty Quantification for SLS Aerodynamic Sectional Loads

This paper presents a method for adjusting sectional loads to match target values for integrated force and moment coefficients. In a typical application, the sectional load profile for one flight condition is calculated from Computational Fluid Dynamics (CFD) while the integrated forces and moments are measured in a wind tunnel experiment. These two methods do not generally result in identical predictions, and this leads to an inherent inconsistency between different data products. This paper aims to provide a procedure to remove that inconsistency. A sectional load profile for a launch vehicle splits the rocket into slices along its length and calculates the aerodynamic loading on each slice, which leads to a one-dimensional aerodynamic load profile that is used for structural analysis. Adjusting sectional loads, also known as line loads, is a nontrivial matter due to several consistency constraints. For example, the adjusted sectional normal force profile must be consistent with both the integrated normal force and pitching moment. To avoid such inconsistency issues, this paper presents a method using a Proper Orthogonal Decomposition (POD) to generate basis functions to adjust the sectional load profiles. As a corollary, this correction method enables the creation of an uncertainty quantification for sectional loads that is consistent with the dispersed integrated force and moment database and its uncertainty quantification. Several extensions to this technique, such as applying the method to the surface pressures, are considered.

Quantification↗
NASA NTRS2018

Combined Error and Uncertainty Estimates for CFD Problems

Given input sources of uncertainty, non-intrusive uncertainty propagation methods quantify the uncertainty in output quantities of interest (QoI) by performing a nite number of CFD (Computational Fluid Dynamics) instance realizations needed in the calculation of output statistics. It is well known that this introduces multiple sources of error. CFD codes often utilize finite-dimensional approximation (grids, basis functions, etc.) thus incurring CFD numerical errors often approximately reinterpreted as a statistical bias. Uncertainty propagation methods calculate uncertainty statistics for output quantities of interest using a numerical method (e.g. deterministic quadrature, sampling, etc.) thus incurring UQ (Uncertainty Quantification) numerical errors. Importance of quantifying these errors in large scale scientific computing: How accurate is an output statistic?; How should additional computational resources be invested to further reduce the error in a statistic?

Posteriori↗
NASA NTRS2017

The Failure of Correlation to Describe Out-of-Plane Carbon=Carbon Bending

Carbon-carbon multiply bonded systems are improperly described with standard correlation methods and basis sets. For computations of vibrational modes, the out-of- plane bends can be reported as imaginary at worst or simply too low at best. Utilizing the simplest of aromatic structures (cyclopropenylidene) and various levels of theory, this work diagnoses this known behavior for the first time. A combined 1-particle and n-particle basis set effect conspire to produce these non-physical results. When moving from sp2 to sp3 hybridization in the carbon atoms, the larger number of basis functions overcorrects the energy. This is exacerbated by correlation methods. These allow for occupation of the π and π∗ orbitals in the expanded wave function that combine with the hydrogen s orbitals. As a result, the improperly described space can be further and non-physically stabilized by post-Hartree-Fock correlation. This represents a fundamental problem with at least Hartree-Fock based methods of all flavors in describing carbon. Beyond being a flaw in quantum chemical theory, other repercussions will be present in computations regarding spectroscopy as well as energy and environmental studies where highly-accurate hydrocabon vibrational transitions or thermochemical data are needed.

Fortenberry, Ryan C.↗
NASA NTRS2018

Simulation of Sound Absorption by Scattering Bodies Treated with Acoustic Liners Using a Time-Domain Boundary Element Method

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction, and are characterized by a complex valued frequency-dependent impedance, Z(w). Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering of geometric bodies treated with acoustic liners. This work uses an admittance boundary condition where the admittance, Y(w), is defined to be the inverse of impedance, i.e., Y(w) = 1/Z(w). An admittance boundary condition will be derived and coupled with a time domain boundary integral equation. The solution will be obtained iteratively using spatial and temporal basis functions and will allow for acoustic scattering problems to be modeled with geometries consisting of both unlined and soft surfaces. Stability will be demonstrated through eigenvalue analysis.

Pizzo, Michelle E.↗