Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

A Priori Subgrid Analysis of Temporal Mixing Layers with Evaporating Droplets

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5. In the DNS, the gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. The Large Eddy Simulation (LES) equations corresponding to the DNS are first derived, and key assumptions in deriving them are first confirmed by computing the terms using the DNS database. Since LES of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be the sum of the filtered variables and a correction based on the filtered standard deviation; this correction is then computed from the Subgrid Scale (SGS) standard deviation. This model predicts the unfiltered variables at droplet locations considerably better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: the Smagorinsky approach, the Gradient model and the Scale-Similarity formulation. When the proportionality constant inherent in the SGS models is properly calculated, the Gradient and Scale-Similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗
NASA NTRS2019

Real-Time Estimation of Bare-Airframe Frequency Responses from Closed-Loop Data and Multisine Inputs

A method is presented for computing frequency responses of multiple-input multiple-output bare-airframe dynamics from flight test data containing feedback control and/or mixing of control effectors. Orthogonal phase-optimized multisines are used to simultaneously excite each input with unique harmonic frequencies, at which frequency responses are computed as ratios of output-to-input Fourier transform data. The confounding effects of feedback and mixing for frequency response estimation are resolved by interpolating the frequency responses among all the harmonic frequencies. The method can be run in batch for post-flight analysis, or in real time as the aircraft is flying. The effectiveness of the method was verified using closed-loop simulations of the subscale NASA T-2 generic transport airplane. The method was also demonstrated using flight test data from the X-56A MUTT aeroelastic airplane, which was flown with feedback control and mixing.

Jared A Grauer↗
NASA NTRS1991

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗
NASA NTRS2017

Development of Multi-Sensor Global Cloud and Radiance Composites for Earth Radiation Budget Monitoring from DSCOVR

The Deep Space Climate Observatory (DSCOVR) enables analysis of the daytime Earth radiation budget via the onboard Earth Polychromatic Imaging Camera (EPIC) and National Institute of Standards and Technology Advanced Radiometer (NISTAR). Radiance observations and cloud property retrievals from low earth orbit and geostationary satellite imagers have to be co-located with EPIC pixels to provide scene identification in order to select anisotropic directional models needed to calculate shortwave and longwave fluxes. A new algorithm is proposed for optimal merging of selected radiances and cloud properties derived from multiple satellite imagers to obtain seamless global hourly composites at 5-kilometer resolution. An aggregated rating is employed to incorporate several factors and to select the best observation at the time nearest to the EPIC measurement. Spatial accuracy is improved using inverse mapping with gradient search during reprojection and bicubic interpolation for pixel resampling. The composite data are subsequently remapped into EPIC-view domain by convolving composite pixels with the EPIC point spread function (PSF) defined with a half-pixel accuracy. PSF-weighted average radiances and cloud properties are computed separately for each cloud phase. The algorithm has demonstrated contiguous global coverage for any requested time of day with a temporal lag of under 2 hours in over 95 percent of the globe.

Khlopenkov, Konstantin↗
NASA NTRS2025

Robust Trajectory Optimization for NRHO Rendezvous Using SPICE Kernel Relative Motion

In this paper, robust optimization is performed on trajectory correction maneuvers during the lunar lander return phase of an Artemis mission, treating the trajectory from one hour after low lunar orbit departure to arrival in the vicinity of the lunar Gateway as a relative motion problem. To enable rapid stochastic optimization techniques requiring many candidate trajectories, SPICE kernel relative motion as implemented by the Quadratic Interpolated State Transition (QIST) system is used as the underlying dynamics propagation. The optimization is performed with a genetic optimizer using linear covariance (LinCov) software in a simplified operational context, taking into account the availability of navigation sensors with varying measurement models, ranges, and accuracies. No numerical integration is used, since the relative motion around Gateway is fully characterized with the a priori computation of the QIST coefficients. Maneuver placements are computed to optimize the minimum 3σ delta-v of the trajectory, the position dispersion at a target point, and a convex combination of these two metrics. An order of magnitude runtime improvement is provided over legacy methods with less than 10% error introduced. All QIST results are shown to be in-family with legacy methods. The tradespace for optimal delta-v design is found to range from 77.0 to 93.9 m/s, while the range of optimal dispersion is between 1.4 and 11.7 km.

Relative Motion↗
NASA NTRS2011

In Situ Trace Element Measurements on Roda and the Origin of Diogenites

The origin of diogenites remains poorly understood. A recent model interprets many diogenites to have been formed from melts that were derived by remelting initial magma ocean cumulates, and these penultimate parent melts were then contaminated by melts derived from remelting of the basaltic (eucritic) crust to form the ultimate diogenite parent melts [1] (hereafter the remelting model). This is a very complicated petrogenesis that has profound implications for the geological evolution of 4 Vesta if correct. This model was developed based on trace element analyses of bulk rock samples that had been leached in acids to remove phosphates; the compositions of the residues were interpreted to be close to those of cumulus orthopyroxenes plagioclase, chromite and olivine [1]. In situ measurements of phases in diogenites can be used to test this model. We have begun a campaign of laser ablation ICP-MS of orthopyroxene grains in diogenites for this purpose. Here we report our first results on one diogenite, Roda. We have determined a suite of trace lithophile elements on nine, mm-sized pyroxene grains separated from Roda that have previously been studied [2, 3]. A key observation supporting the remelting model is the very low Eu/Eu* of leached residues; values too low to represent orthopyroxene that crystallized from melts with chondritic Sm/Eu and Gd/Eu [1]. (Eu* = Eu interpolated from REE diagrams.) Crustal remelts have low Sm/Eu and Gd/Eu, and orthopyroxenes that crystallized from parent melts contaminated by them would have very low Eu/Eu* [1]. Roda grains have Eu/Eu* of 0.243 to 0.026; the latter a value lower than any measured on bulk diogenite leached residues (0.041) [1]. There is a general negative correlation between Eu/Eu* and some incompatible elements (Zr, Nb, Hf), but not others (LREE). This appears inconsistent with the remelting model as it would suggest an evolving parent melt with La de-creasing as Zr increased and Eu/Eu* decreased. Grain R-15 includes trace-element-rich trapped melt phases [2, 3]. This grain has the highest Eu/Eu* and LREE contents, indicating that the trapped melt had a high Eu/Eu*. Thus, our first data on one diogenite do not provide support for the remelting model [1]. Roda is unusual in that its orthopyroxene grains show wide ranges in trace element contents [4]. Previous in situ REE analyses of grain R-15 did not reveal evidence for subsolidus equilibration with trace-element-rich trapped melt phases, and led to the suggestion that Roda may be polymict, with different grains representing different lithologies of diverse compositions [3]. Thus, based on our results on Roda, it is perhaps premature to abandon the remelting model. In situ measurements on a suite of diogenites is planned to further address this issue.

Mittlefehldt, David W.↗
NASA NTRS1992

Tomographic separation of composite spectra. I - The components of Plaskett's Star

The UV photospheric lines of Plaskett's Star (HD 47129), a 14.4 day period, double lined O-type spectroscopic binary were analyzed. Archival data from IUE (17 spectra well distributed in orbital phase) were analyzed with several techniques. A cross correlation analysis, which showed that the secondary produces significant lines in the UV, indicates that the mass ratio is q = 1.18 +/- 0.12 (secondary slightly more massive). A tomography algorithm was used to produce the separate spectra of the two stars in six spectral regions. The interpolated spectral classifications of the primary and secondary, 07.3 I and 06.2 I, respectively, were estimated through a comparison of UV line ratios with those in spectral standard stars. The intensity ratio of the stars in the UV is 0.53 +/- 0.05 (primary brighter). The secondary lines appear rotationally broadened, and the projected rotational velocity V sin i for this star is estimated to be 310 +/- 20 km/s. The possible evolutionary history of this system is discussed through a comparison of the positions of the components and evolutionary tracks in the H-R diagram.

Bagnuolo, William G., Jr.↗
NASA NTRS1991

Tomographic separation of composite spectra. The components of Plaskett's Star

The UV photospheric lines of Plaskett's Star (HD 47129), a 14.4 day period, double lined O-type spectroscopic binary were analyzed. Archival data from IUE (17 spectra well distributed in orbital phase) were analyzed with several techniques. A cross correlation analysis, which showed that the secondary produces significant lines in the UV, indicates that the mass ratio is q = 1.18 + or - 0.12 (secondary slightly more massive). A tomography algorithm was used to produce the separate spectra of the two stars in six spectral regions. The interpolated spectral classifications of the primary and secondary, 07.3 I and 06.2 I, respectively, were estimated through a comparison of UV line ratios with those in spectral standard stars. The intensity ratio of the stars in the UV is 0.53 + or - 0.05 (primary brighter). The secondary lines appear rotationally broadened, and the projected rotational velocity V sin i for this star is estimated to be 310 + or - 20 km/s. The possible evolutionary history of this system is discussed through a comparison of the positions of the components and evolutionary tracks in the H-R diagram.

Bagnuolo, William G., Jr.↗
NASA NTRS2002

Diffusion NMR methods applied to xenon gas for materials study

We report initial NMR studies of (i) xenon gas diffusion in model heterogeneous porous media and (ii) continuous flow laser-polarized xenon gas. Both areas utilize the pulsed gradient spin-echo (PGSE) techniques in the gas phase, with the aim of obtaining more sophisticated information than just translational self-diffusion coefficients--a brief overview of this area is provided in the Introduction. The heterogeneous or multiple-length scale model porous media consisted of random packs of mixed glass beads of two different sizes. We focus on observing the approach of the time-dependent gas diffusion coefficient, D(t) (an indicator of mean squared displacement), to the long-time asymptote, with the aim of understanding the long-length scale structural information that may be derived from a heterogeneous porous system. We find that D(t) of imbibed xenon gas at short diffusion times is similar for the mixed bead pack and a pack of the smaller sized beads alone, hence reflecting the pore surface area to volume ratio of the smaller bead sample. The approach of D(t) to the long-time limit follows that of a pack of the larger sized beads alone, although the limiting D(t) for the mixed bead pack is lower, reflecting the lower porosity of the sample compared to that of a pack of mono-sized glass beads. The Pade approximation is used to interpolate D(t) data between the short- and long-time limits. Initial studies of continuous flow laser-polarized xenon gas demonstrate velocity-sensitive imaging of much higher flows than can generally be obtained with liquids (20-200 mm s-1). Gas velocity imaging is, however, found to be limited to a resolution of about 1 mm s-1 owing to the high diffusivity of gases compared with liquids. We also present the first gas-phase NMR scattering, or diffusive-diffraction, data, namely flow-enhanced structural features in the echo attenuation data from laser-polarized xenon flowing through a 2 mm glass bead pack. c2002 John Wiley & Sons, Ltd.

NASA Discipline Life Sciences Technologies↗
NASA NTRS1988

Least-Squares Frequency-Acquisition Algorithm

Algorithm finds frequency and phase of sinudoidal signal in presence of noise. Algorithm is special case of more-general, adaptive-paramenter-estimation techniques. Computational requirements of algorithm comparable to corresponding fast-Fourier-transform (FFT) algorithm. Algorithm works directly in time domain, whereas FFT algorithm transforms data into frequency domain for estimation and detection and requires secondary algorithm to interpolate between frequencies.

Kumar, Rajendra↗
NASA NTRS1994

Modeling the bidirectional reflectance distribution function of mixed finite plant canopies and soil

An analytical model of the bidirectional reflectance for optically semi-infinite plant canopies has been extended to describe the reflectance of finite depth canopies contributions from the underlying soil. The model depends on 10 independent parameters describing vegetation and soil optical and structural properties. The model is inverted with a nonlinear minimization routine using directional reflectance data for lawn (leaf area index (LAI) is equal to 9.9), soybeans (LAI, 2.9) and simulated reflectance data (LAI, 1.0) from a numerical bidirectional reflectance distribution function (BRDF) model (Myneni et al., 1988). While the ten-parameter model results in relatively low rms differences for the BRDF, most of the retrieved parameters exhibit poor stability. The most stable parameter was the single-scattering albedo of the vegetation. Canopy albedo could be derived with an accuracy of less than 5% relative error in the visible and less than 1% in the near-infrared. Sensitivity were performed to determine which of the 10 parameters were most important and to assess the effects of Gaussian noise on the parameter retrievals. Out of the 10 parameters, three were identified which described most of the BRDF variability. At low LAI values the most influential parameters were the single-scattering albedos (both soil and vegetation) and LAI, while at higher LAI values (greater than 2.5) these shifted to the two scattering phase function parameters for vegetation and the single-scattering albedo of the vegetation. The three-parameter model, formed by fixing the seven least significant parameters, gave higher rms values but was less sensitive to noise in the BRDF than the full ten-parameter model. A full hemispherical reflectance data set for lawn was then interpolated to yield BRDF values corresponding to advanced very high resolution radiometer (AVHRR) scan geometries collected over a period of nine days. The resulting parameters and BRDFs are similar to those for the full sampling geometry, suggesting that the limited geometry of AVHRR measurements might be used to reliably retrieve BRDF and canopy albedo with this model.

Schluessel, G.↗
NASA NTRS1981

A program to evaluate a control system based on feedback of aerodynamic pressure differentials. Part 2: Data report for phase 1 wind tunnel test

Wind tunnel test data are presented. The purpose of the test was to determine the sensor location for sensitivity to flap deflection and to determine the range required of the sensor. The calculations followed to correct the raw pressure measurements, the corrected measurements, which represent the static pressures along the airfoil surface, and the static pressures tabulated and in coefficient form are described. However, the chordwise locations of the upper and lower surface taps do not coincide. An interpolation routine is performed to find the lower surface pressure coefficients at the upper surface tap locations. The difference, C sub P(lower) - C sub P(upper), is calculated and listed and the data are reorganized with respect to each tap location. A numerical regression performed to quantify the slopes of the graphs allows more accurate comparison of sensitivity and linearity.

Hrabak, R. R.↗
NASA NTRS2011

Land Data Assimilation of Satellite-Based Soil Moisture Products Using the Land Information System Over the NLDAS Domain

This presentation will include results from data assimilation simulations using the NASA-developed Land Information System (LIS). Using the ensemble Kalman filter in LIS, two satellite-based soil moisture products from the AMSR-E instrument were assimilated, one a NASA-based product and the other from the Land Parameter Retrieval Model (LPRM). The domain and land-surface forcing data from these simulations were from the North American Land Data Assimilation System Phase-2, over the period 2002-2008. The Noah land-surface model, version 3.2, was used during the simulations. Changes to estimates of land surface states, such as soil moisture, as well as changes to simulated runoff/streamflow will be presented. Comparisons over the NLDAS domain will also be made to two global reference evapotranspiration (ET) products, one an interpolated product based on FLUXNET tower data and the other a satellite- based algorithm from the MODIS instrument. Results of an improvement metric show that assimilating the LPRM product improved simulated ET estimates while the NASA-based soil moisture product did not.

Mocko, David M.↗
NASA NTRS2006

Rapid Aeroelastic Analysis of Blade Flutter in Turbomachines

The LINFLUX-AE computer code predicts flutter and forced responses of blades and vanes in turbomachines under subsonic, transonic, and supersonic flow conditions. The code solves the Euler equations of unsteady flow in a blade passage under the assumption that the blades vibrate harmonically at small amplitudes. The steady-state nonlinear Euler equations are solved by a separate program, then equations for unsteady flow components are obtained through linearization around the steady-state solution. A structural-dynamics analysis (see figure) is performed to determine the frequencies and mode shapes of blade vibrations, a preprocessor interpolates mode shapes from the structural-dynamics mesh onto the LINFLUX computational-fluid-dynamics mesh, and an interface code is used to convert the steady-state flow solution to a form required by LINFLUX. Then LINFLUX solves the linearized equations in the frequency domain to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. A post-processor uses the unsteady pressures to calculate generalized aerodynamic forces, response amplitudes, and eigenvalues (which determine the flutter frequency and damping). In comparison with the TURBO-AE aeroelastic-analysis code, which solves the equations in the time domain, LINFLUX-AE is 6 to 7 times faster.

Trudell, J. J.↗
NASA NTRS1986

Objective analysis of tidal fields in the Atlantic and Indian Oceans

An objective analysis technique has been developed to extrapolate tidal amplitudes and phases over entire ocean basins using existing gauge data and the altimetric measurements which are now beginning to be provided by satellite oceanography. The technique was previously tested in the Lake Superior basin. The method has now been developed and applied in the Atlantic-Indian ocean basins using a 6 deg x 6 deg grid to test its essential features. The functions used in the interpolation are the eigenfunctions of the velocity potential (Proudman functions) which are computed numerically from a knowledge of the basin's bottom topography, the horizontal plan form and the necessary boundary conditions. These functions are characteristic of the particular basin. The gravitational normal modes of the basin are computed as part of the investigation, they are used to obtain the theoretical forced solutions for the tidal constituents, the latter provide the simulated data for the testing of the method and serve as a guide in choosing the most energetic modes for the objective analysis. The results of the objective analysis of the M2 and K1 tidal constituents indicate the possibility of recovering the tidal signal with a degree of accuracy well within the error bounds of present day satellite techniques.

Sanchez, B. V.↗
NASA NTRS2011

Multiangle Implementation of Atmospheric Correction (MAIAC):: Radiative Transfer Basis and Look-up Tables - 1

This paper describes a radiative transfer basis of the algorithm MAIAC which performs simultaneous retrievals of atmospheric aerosol and bidirectional surface reflectance from the Moderate Resolution Imaging Spectroradiometer (MODIS). The retrievals are based on an accurate semianalytical solution for the top-of-atmosphere reflectance expressed as an explicit function of three parameters of the Ross-Thick Li-Sparse model of surface bidirectional reflectance. This solution depends on certain functions of atmospheric properties and geometry which are precomputed in the look-up table (LUT). This paper further considers correction of the LUT functions for variations of surface pressure/height and of atmospheric water vapor, which is a common task in the operational remote sensing. It introduces a new analytical method for the water vapor correction of the multiple ]scattering path radiance. It also summarizes the few basic principles that provide a high efficiency and accuracy of the LUT ]based radiative transfer for the aerosol/surface retrievals and optimize the size of LUT. For example, the single-scattering path radiance is calculated analytically for a given surface pressure and atmospheric water vapor. The same is true for the direct surface-reflected radiance, which along with the single-scattering path radiance largely defines the angular dependence of measurements. For these calculations, the aerosol phase functions and kernels of the surface bidirectional reflectance model are precalculated at a high angular resolution. The other radiative transfer functions depend rather smoothly on angles because of multiple scattering and can be calculated at coarser angular resolution to reduce the LUT size. At the same time, this resolution should be high enough to use the nearest neighbor geometry angles to avoid costly three ]dimensional interpolation. The pressure correction is implemented via linear interpolation between two LUTs computed for the standard and reduced pressure levels. A linear mixture and a modified linear mixture methods are used to represent different aerosol types in the aerosol/surface retrievals from several base models of the fine and coarse aerosol fractions. In summary, the developed LUT algorithm allows fast high-accuracy simulations of the outgoing radiance with full variability of the atmospheric and surface bidirectional reflectance properties for the aerosol/surface remote sensing.

Lyapustin, Alexei↗
NASA NTRS2002

Parametric Identification of Nonlinear Dynamical Systems

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Feeny, Brian↗
NASA NTRS2000

The Linear Bicharacteristic Scheme for Computational Electromagnetics

The upwind leapfrog or Linear Bicharacteristic Scheme (LBS) has previously been implemented and demonstrated on electromagnetic wave propagation problems. This paper extends the Linear Bicharacteristic Scheme for computational electromagnetics to treat lossy dielectric and magnetic materials and perfect electrical conductors. This is accomplished by proper implementation of the LBS for homogeneous lossy dielectric and magnetic media, and treatment of perfect electrical conductors (PECs) are shown to follow directly in the limit of high conductivity. Heterogeneous media are treated through implementation of surface boundary conditions and no special extrapolations or interpolations at dielectric material boundaries are required. Results are presented for one-dimensional model problems on both uniform and nonuniform grids, and the FDTD algorithm is chosen as a convenient reference algorithm for comparison. The results demonstrate that the explicit LBS is a dissipation-free, second-order accurate algorithm which uses a smaller stencil than the FDTD algorithm, yet it has approximately one-third the phase velocity error. The LBS is also more accurate on nonuniform grids.

Beggs, John H.↗