Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Discretization error”

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 649 records · Page 36

Evolution of the Scope and Capabilities of Uplink Support Software for Mars Surface Operations

In January of 2004 both of the Mars Exploration Rover spacecraft landed safely, initiating daily surface operations at the Jet Propulsion Laboratory for what was anticipated to be approximately three months of mobile exploration. The longevity of this mission, still ongoing after ten years, has provided not only a tremendous return of scientific data but also the opportunity to refine and improve the methodology by which robotic Mars surface missions are commanded. Since the landing of the Mars Science Laboratory spacecraft in August of 2012, this methodology has been successfully applied to operate a Martian rover which is both similar to, and quite different from, its predecessors. For MER and MSL, daily uplink operations can be most broadly viewed as converting the combined interests of both the science and engineering teams into a spacecraft-safe set of transmittable command files. In order to accomplish these ends a discrete set of mission-critical software tools were developed which not only allowed for conformation to established JPL standards and practices but also enabled innovative technologies specific to each mission. Although these primary programs provided the requisite capabilities for meeting the high-level goals of each distinct phase of the uplink process, there was little in the way of secondary software to support the smooth flow of data from one phase to the next. In order to address this shortcoming a suite of small software tools was developed to aid in phase transitions, as well as to automate some of the more laborious and error-prone aspects of uplink operations. This paper describes the evolution of this software suite, from its initial attempts to merely shorten the duration of the operator's shift, to its current role as an indispensable tool enforcing workflow of the uplink operations process and agilely responding to the new and unexpected challenges of missions which can, and have, lasted many years longer than originally anticipated.

CoUGAR↗

Adaptive control based on retrospective cost optimization

A discrete-time adaptive control law for stabilization, command following, and disturbance rejection that is effective for systems that are unstable, MIMO, and/or nonminimum phase. The adaptive control algorithm includes guidelines concerning the modeling information needed for implementation. This information includes the relative degree, the first nonzero Markov parameter, and the nonminimum-phase zeros. Except when the plant has nonminimum-phase zeros whose absolute value is less than the plant's spectral radius, the required zero information can be approximated by a sufficient number of Markov parameters. No additional information about the poles or zeros need be known. Numerical examples are presented to illustrate the algorithm's effectiveness in handling systems with errors in the required modeling data, unknown latency, sensor noise, and saturation.

Santillo, Mario A.↗

An Algebraic Multigrid Solver for Navier-Stokes Problems in the Discrete Second-Order Approximation

An algebraic multigrid scheme is presented for solving the discrete Navier-Stokes equations to second-order accuracy using the defect-correction method. Solutions have been obtained for problems involving both structured and unstructured meshes, with the resolution and resolution grading controlled by global and local mesh refinements. The solver is efficient and robust to the extent that no underrelaxation of variables has been required to ensure convergence, but rates of convergence can be improved with small amounts of underrelaxation of the velocity-pressure coupling. Provided that the computational mesh can resolve the flow field, convergence characteristics are almost mesh independent. Rates of convergence actually improve with refinement, asymptotically approaching mesh independent values. For extremely coarse meshes where dispersive truncation errors would be expected to prevent convergence (or even induce divergence), solutions can still be obtained by using explicit underrelaxation in the iterative cycle.

Webster, R.↗

Retrieval of Water Vapor Column Abundance and Aerosol Properties from ChemCam Passive Sky Spectroscopy

We derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. This paper covers the methodology and initial results for water vapor and also provides preliminary results for aerosols. The data set presented here includes the results of 113 observations spanning from Mars Year 31 L(s) = 291 deg (March 30, 2013) to Mars Year 33 L(s) = 127 deg (March 24, 2016). Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption bands. The retrieval proceeds by first fitting the continuum of the ratio of the two elevation angles to solve for aerosol properties, and then fitting the continuum-removed ratio to solve for gas abundances. The final step of the retrieval makes use of the observed CO 2 absorptions and the known CO 2 abundance to correct the retrieved water vapor abundance for the effects of the vertical distribution of scattering aerosols and to derive an aerosol scale height parameter. Our water vapor results give water vapor column abundance with a precision of +/-0.6 precipitable microns and systematic errors no larger than +/-0.3 precipitable microns, assuming uniform vertical mixing. The ChemCam-retrieved water abundances show, with only a few exceptions, the same seasonal behavior and the same timing of seasonal minima and maxima as the TES, CRISM, and REMS-H data sets that we compare them to. However ChemCam-retrieved water abundances are generally lower than zonal and regional scale from-orbit water vapor data, while at the same time being significantly larger than pre-dawn REMS-H abundances. Pending further analysis of REMS-H volume mixing ratio uncertainties, the differences between ChemCam and REMS-H pre-dawn mixing ratios appear to be much too large to be explained by large scale circulations and thus they tend to support the hypothesis of substantial diurnal interactions of water vapor with the surface. Our preliminary aerosol results, meanwhile, show the expected seasonal pattern in dust particle size but also indicate a surprising interannual increase in water-ice cloud opacities.

Mcconnochie, Timothy H.↗

An Evaluation Technique for an F/A-18 Aircraft Loads Model Using F/A-18 Systems Research Aircraft Flight Data

A limited evaluation of the F/A-18 baseline loads model was performed on the Systems Research Aircraft at NASA Dryden Flight Research Center (Edwards, California). Boeing developed the F/A-18 loads model using a linear aeroelastic analysis in conjunction with a flight simulator to determine loads at discrete locations on the aircraft. This experiment was designed so that analysis of doublets could be used to establish aircraft aerodynamic and loads response at 20 flight conditions. Instrumentation on the right outboard leading edge flap, left aileron, and left stabilator measured the hinge moment so that comparisons could be made between in-flight-measured hinge moments and loads model-predicted values at these locations. Comparisons showed that the difference between the loads model-predicted and in-flight-measured hinge moments was up to 130 percent of the flight limit load. A stepwise regression technique was used to determine new loads derivatives. These derivatives were placed in the loads model, which reduced the error to within 10 percent of the flight limit load. This paper discusses the flight test methodology, a process for determining loads coefficients, and the direct comparisons of predicted and measured hinge moments and loads coefficients.

Olney, Candida D.↗

A performance comparison of integration algorithms in simulating flexible structures

Asymptotic formulas for the characteristic root errors as well as transfer function gain and phase errors are presented for a number of traditional and new integration methods. Normalized stability regions in the lambda h plane are compared for the various methods. In particular, it is shown that a modified form of Euler integration with root matching is an especially efficient method for simulating lightly-damped structural modes. The method has been used successfully for structural bending modes in the real-time simulation of missiles. Performance of this algorithm is compared with other special algorithms, including the state-transition method. A predictor-corrector version of the modified Euler algorithm permits it to be extended to the simulation of nonlinear models of the type likely to be obtained when using the discretized structure approach. Performance of the different integration methods is also compared for integration step sizes larger than those for which the asymptotic formulas are valid. It is concluded that many traditional integration methods, such as RD-4, are not competitive in the simulation of lightly damped structures.

Howe, R. M.↗

Mathematical Inversion of Lightning Data: Techniques and Applications

A survey of some interesting mathematical inversion studies dealing with radio, optical, and electrostatic measurements of lightning are presented. A discussion of why NASA is interested in lightning, what specific physical properties of lightning are retrieved, and what mathematical techniques are used to perform the retrievals are discussed. In particular, a relatively new multi-station VHF time-of-arrival (TOA) antenna network is now on-line in Northern Alabama and will be discussed. The network, called the Lightning Mapping Array (LMA), employs GPS timing and detects VHF radiation from discrete segments (effectively point emitters) that comprise the channel of lightning strokes within cloud and ground flashes. The LMA supports on-going ground-validation activities of the low Earth orbiting Lightning Imaging Sensor (LIS) satellite developed at NASA Marshall Space Flight Center (MSFC) in Huntsville, Alabama. The LMA also provides detailed studies of the distribution and evolution of thunderstorms and lightning in the Tennessee Valley, and offers interesting comparisons with other meteorological/geophysical datasets. In order to take full advantage of these benefits, it is essential that the LMA channel mapping accuracy (in both space and time) be fully characterized and optimized. A new channel mapping retrieval algorithm is introduced for this purpose. To characterize the spatial distribution of retrieval errors, the algorithm has been applied to analyze literally tens of millions of computer-simulated lightning VHF point sources that have been placed at various ranges, azimuths, and altitudes relative to the LMA network. Statistical results are conveniently summarized in high-resolution, color-coded, error maps.

Koshak, William↗

Utilization of O4 Slant Column Density to Derive Aerosol Layer Height from a Space-Borne UV-Visible Hyperspectral Sensor: Sensitivity and Case Study

The sensitivities of oxygen-dimer (O4) slant column densities (SCDs) to changes in aerosol layer height are investigated using the simulated radiances by a radiative transfer model, the linearized pseudo-spherical vector discrete ordinate radiative transfer (VLIDORT), and the differential optical absorption spectroscopy (DOAS) technique. The sensitivities of the O4 index (O4I), which is defined as dividing O4 SCD by 10(sup 40) molecules (sup 2) per centimeters(sup -5), to aerosol types and optical properties are also evaluated and compared. Among the O4 absorption bands at 340, 360, 380, and 477 nanometers, the O4 absorption band at 477 nanometers is found to be the most suitable to retrieve the aerosol effective height. However, the O4I at 477 nanometers is significantly influenced not only by the aerosol layer effective height but also by aerosol vertical profiles, optical properties including single scattering albedo (SSA), aerosol optical depth (AOD), particle size, and surface albedo. Overall, the error of the retrieved aerosol effective height is estimated to be 1276, 846, and 739 meters for dust, non-absorbing, and absorbing aerosol, respectively, assuming knowledge on the aerosol vertical distribution shape. Using radiance data from the Ozone Monitoring Instrument (OMI), a new algorithm is developed to derive the aerosol effective height over East Asia after the determination of the aerosol type and AOD from the MODerate resolution Imaging Spectroradiometer (MODIS). About 80 percent of retrieved aerosol effective heights are within the error range of 1 kilometer compared to those obtained from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) measurements on thick aerosol layer cases.

Air Pollution↗

MODTRAN3: An update and recent validations against airborne high resolution interferometer measurements

MODTRAN, the Moderate Resolution Atmospheric Radiance and Transmittance Model, encompasses all the capabilities of LOWTRAN 7, the widely used 20 cm(exp -1) resolution radiance code, but incorporates a much more sensitive molecular band model with 2 cm(exp -1) resolution. MODTRAN contains many important elements that other band model based radiative transfer codes do not incorporate. It shares with FASCODE: spherical geometry, single and multiple scattering default atmospheric profile descriptors (gases, aerosols, clouds, fogs, and rain), and molecular continua (H2O, CO2, O3, O2, N2). In addition, it can calculate the solar/lunar direct and scattered radiation. MODTRAN3 was released to the general public in November 1994. It has several important features that the previous version, MODTRAN2, does not have. Chloro-fluorocarbon (CFC) and related heavy molecules (whose spectroscopic properties first appear on the HITRAN92 data base as temperature-dependent cross sections) have been incorporated into pseudo-band models, with provision for using both default and user supplied profiles. The addition of SO2 and O2 in the UV, along with upgraded ozone Chappuis bands in the visible is also part of MODTRAN3. An improved multiple scattering algorithm, the DIScrete Ordinate Radiative Transfer (DISORT) has also been incorporated into MODTRAN3. MODTRAN is very fast: simple timing runs of MODTRAN3 vs. FASCOD3 show an improvement of more than a factor of 100 for a typical 500 cm(exp -1) spectral interval and comparable vertical layering. Speed is an important consideration in heating/cooling rates calculations, where a large number of radiative transfer calculations are needed. The MODTRAN3 used in this study is based on HITRAN92, but as mentioned, above, it will be upgraded to HITRAN94 upon its release at the end of 1994. MODTRAN has been adopted by: some researchers in the AVIRIS program as one radiative transfer code to derive surface reflectance from AVIRIS measurements. The accuracy of the code is very important because any errors in the radiative transfer calculation will directly translate into errors in the derived surface reflectance. In this paper, the new solar irradiance calculated by Kurucz, which is adopted in MODTRAN3, will be presented. Recent validations of MODTRAN3 with airborne high resolution interferometer measurements over ocean will be discussed. Good agreeement between model calculations and measurements was achieved.

Anderson, Gail P.↗

A Simple Tool for the Design and Analysis of Multiple-Reflector Antennas in a Multi-Disciplinary Environment

The process of designing and analyzing a multiple-reflector system has traditionally been time-intensive, requiring large amounts of both computational and human time. At many frequencies, a discrete approximation of the radiation integral may be used to model the system. The code which implements this physical optics (PO) algorithm was developed at the Jet Propulsion Laboratory. It analyzes systems of antennas in pairs, and for each pair, the analysis can be computationally time-consuming. Additionally, the antennas must be described using a local coordinate system for each antenna, which makes it difficult to integrate the design into a multi-disciplinary framework in which there is traditionally one global coordinate system, even before considering deforming the antenna as prescribed by external structural and/or thermal factors. Finally, setting up the code to correctly analyze all the antenna pairs in the system can take a fair amount of time, and introduces possible human error. The use of parallel computing to reduce the computational time required for the analysis of a given pair of antennas has been previously discussed. This paper focuses on the other problems mentioned above. It will present a methodology and examples of use of an automated tool that performs the analysis of a complete multiple-reflector system in an integrated multi-disciplinary environment (including CAD modeling, and structural and thermal analysis) at the click of a button. This tool, named MOD Tool (Millimeter-wave Optics Design Tool), has been designed and implemented as a distributed tool, with a client that runs almost identically on Unix, Mac, and Windows platforms, and a server that runs primarily on a Unix workstation and can interact with parallel supercomputers with simple instruction from the user interacting with the client.

Katz, Daniel S.↗

A Simple Tool for the Design and Analysis of Multiple-reflector Antennas in a Multi-disciplinary Environment

The process of designing and analyzing a multiple-reflector system has traditionally been time-intensive, requiring large amounts of both computational and human time. At many frequencies, a discrete approximation of the radiation integral may be used to model the system. The code which implements this physical optics (PO) algorithm was developed at the Jet Propulsion Laboratory. It analyzes systems of antennas in pairs, and for each pair, the analysis can be computationally time-consuming. Additionally, the antennas must be described using a local coordinate system for each antenna, which makes it difficult to integrate the design into a multi-disciplinary framework in which there is traditionally one global coordinate system, even before considering deforming the antenna as prescribed by external structural and/or thermal factors. Finally, setting up the code to correctly analyze all the antenna pairs in the system can take a fair amount of time, and introduces possible human error.

Katz, Daniel, S.↗

Pulse Response of the Short-Wave Infrared Detection System of the Ocean Color Instrument for the NASA PACE Mission

The Ocean Color Instrument (OCI) on NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission is a hyperspectral Earth imager with a spatial resolution of 1 km x 1 km and a spectral resolution of 5 nm in 2.5 nm steps over 342-887 nm. In addition, OCI provides 7 discrete bands in the 940-2260 nm Short-Wave InfraRed (SWIR) range. The front-end optical imager is a rotating mirror-based system that images the ground scene onto a slit with an instantaneous field of view of 16 km x 1 km. For the SWIR bands, the slit-image is re-imaged onto a 16x1 micro-lens array that effectively acts as the focal plane since each lens element is fiber coupled to wavelength filtered InGaAs and HgCdTe Photo Diodes (PDs). The pulse response of the detection system is critical to OCI SWIR performance. We find that PDs introduce an inherent slow tail in the pulse response due to slow diffusion moving carriers in their n and p regions. We show that this introduces response errors ranging from 1 down to 0.01 % for up to tens of science pixels after the pulse depending on the PD design and materials. It is shown that the response is distinctly different for the InGaAs and HgCdTe PDs. We explain how the front-end design can further increase this error. Finally, we detail the cause of the slow pulse response tail, how to model it, its impact on OCI performance and how it is characterized and corrected to meet OCI requirements.

ocean color↗

Pulse Response of the Short-Wave Infrared Detection System of the Ocean Color Instrument for the NASA Pace Mission

The Ocean Color Instrument (OCI) on NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission is a hyperspectral Earth imager with a spatial resolution of 1 km x 1 km and a spectral resolution of 5 nm in 2.5 nm steps over 342-887 nm. In addition, OCI provides 7 discrete bands in the 940-2260 nm Short-Wave InfraRed (SWIR) range. The front-end optical imager is a rotating mirror-based system that images the ground scene onto a slit with an instantaneous field of view of 16 km x 1 km. For the SWIR bands, the slit-image is re-imaged onto a 16x1 micro-lens array that effectively acts as the focal plane since each lens element is fiber coupled to wavelength filtered InGaAs and HgCdTe Photo Diodes (PDs). The pulse response of the detection system is critical to OCI SWIR performance. We find that PDs introduce an inherent slow tail in the pulse response due to slow diffusion moving carriers in their n and p regions. We show that this introduces response errors ranging from 1 down to 0.01 % for up to tens of science pixels after the pulse depending on the PD design and materials. It is shown that the response is distinctly different for the InGaAs and HgCdTe PDs. We explain how the front-end design can further increase this error. Finally, we detail the cause of the slow pulse response tail, how to model it, its impact on OCI performance and how it is characterized and corrected to meet OCI requirements.

ocean color↗

Efficient massively parallel simulation of dynamic channel assignment schemes for wireless cellular communications

Fast, efficient parallel algorithms are presented for discrete event simulations of dynamic channel assignment schemes for wireless cellular communication networks. The driving events are call arrivals and departures, in continuous time, to cells geographically distributed across the service area. A dynamic channel assignment scheme decides which call arrivals to accept, and which channels to allocate to the accepted calls, attempting to minimize call blocking while ensuring co-channel interference is tolerably low. Specifically, the scheme ensures that the same channel is used concurrently at different cells only if the pairwise distances between those cells are sufficiently large. Much of the complexity of the system comes from ensuring this separation. The network is modeled as a system of interacting continuous time automata, each corresponding to a cell. To simulate the model, conservative methods are used; i.e., methods in which no errors occur in the course of the simulation and so no rollback or relaxation is needed. Implemented on a 16K processor MasPar MP-1, an elegant and simple technique provides speedups of about 15 times over an optimized serial simulation running on a high speed workstation. A drawback of this technique, typical of conservative methods, is that processor utilization is rather low. To overcome this, new methods were developed that exploit slackness in event dependencies over short intervals of time, thereby raising the utilization to above 50 percent and the speedup over the optimized serial code to about 120 times.

Greenberg, Albert G.↗

Dust extinction and molecular gas in the dark cloud IC 5146

In this paper we describe a powerful method for mapping the distribution of dust through a molecular cloud using data obtained in large-scale, multiwavelength, infrared imaging surveys. This method combines direct measurements of near-infrared color excess and certain techniques of star counting to derive mean extinctions and map the dust column density distribution through a cloud at higher angular resolutions and greater optical depths than those achieved previously by optical star counting. We report the initial results of the application of this method to a dark cloud complex near the cluster IC 5146, where we have performed coordinated, near-infrared, JHK imaging and (13)CO, C(18)O, and CS millimeter-wave, molecular-line surveys of a large portion of the complex. More than 4000 stars were detected in our JHK survey of the cloud. Of these, all but about a dozen appear to be field stars not associated with the cloud. Star count maps at J band show a striking and detailed anticorrelation between the surface density of J-band sources and CO and CS molecular-line emission. We used the (H-K) colors and positions of nearly 1300 sources to directly measure and map the extinction and thus trace the dust column density through the cloud at an effective angular resolution of 1 min .5. We report an interesting correlation between the measured dispersion in our extinction determinations and the extinction. Modeling this relation indicates that effects of small-scale cloud structure dominate the uncertainties in our measurements. Moreover, we demonstrate that such observations can be used to place constraints on the nature of the spatial distribution of extinction on scales smaller than our resolution. In particular, we show that models in which the dust is distributed uniformly or in discrete high-extinction clumps on scales smaller than (1 min .5) are inconsistent with the observations. We have derived extinctions at the same positions and at the same angular resolution (1 min .7) as our molecular-line observations. This enabled a direct comparison of (13)CO, C(18)O, and CS integrated intensities and column densities with A(sub V) for more than 500 positions in the cloud, corresponding to a range in A(sub V) between 0 to 32 mag of extinction. We found the integrated intensities of (13)CO, C(18)O, and CS to be roughly linearly correlated with extinction over different ranges of extinction. However, for all three molecules we find the scatter in the observed relations to be larger than can be accounted for by instrumental error, suggesting that there are large intrinsic variations in the abundances or excitation of the molecules through the cloud. Mean abundances for all the molecules relative to hydrogen were directly derived from the data. The ratio of (13)CO to C(18)O abundances was found to be significantly higher than the terrestrial ratio in regions where extinction is less than 10 mag. In the same region, the dispersion in the abundance ratio is also found to be very large, suggesting that the abundances of one or both molecules are very unstable even at relatively large cloud optical depths. Beyond 10 mag of extinction the abundances of both species appear very stable with their ratio close to the terrestrial value.

Lada, Charles J.↗

Equivalent Plate Analysis of Aircraft Wing with Discrete Source Damage

An equivalent plate procedure is developed to provide a computationally efficient means of matching the stiffness and frequencies of flight vehicle wing structures for prescribed loading conditions. First, the equivalent plate is used to match the stiffness of a stiffened panel without damage and the stiffness of a stiffened panel with damage. For both stiffened panels, the equivalent plate models reproduce the deformation of a corresponding detailed model exactly for the given loading conditions. Once the stiffness was matched, the equivalent plate models were then used to predict the frequencies of the panels. Two analytical procedures using the lumped-mass matrix were used to match the first five frequencies of the corresponding detailed model. In both the procedures, the lumped-mass matrix for the equivalent plate is constructed by multiplying the diagonal terms of the consistent-mass matrix by a proportionality constant. In the first procedure, the proportionality constant is selected such that the total mass of the equivalent plate is equal to that of the detailed model. In the second method, the proportionality constant is selected to minimize the sum of the squares of the errors in a set of pre-selected frequencies between the equivalent plate model and the detailed model. The equivalent plate models reproduced the fundamental first frequency accurately in both the methods. It is observed that changing only the mass distribution in the equivalent plate model did not provide enough flexibility to match all of the frequencies.

Krishnamurthy, T.↗

Evaluation of discrete frequency sound in closed-test-section wind tunnels

The principal objective of this study is to assess the adequacy of linear acoustic theory with an impedance wall boundary condition for modeling the detailed sound field of an acoustic source in a duct. This study compares measurements and calculations of a simple acoustic source in a rectangular concrete duct lined with foam on the walls and anechoic end terminations. Measuring acoustic pressure for 12 wave numbers provides variation in frequency and absorption characteristics of the duct walls. The cases in this study contain low frequencies and low wall absorptions corresponding to measurements of low-frequency helicopter noise in a lined wind tunnel. This regime is particularly difficult to measure in wind tunnels due to high levels of the reverberant field relatively close to the source. Close to the source, where the interference of wall reflections is minimal, correlation is very good. Away from the source, correlation degrades, especially for the lower frequencies. Sensitivity studies show little effect on the predicted results for changes in impedance boundary condition values, source location, measurement location, temperature, and source model for variations spanning the expected measurement error.

Mosher, Marianne↗

An Exact Dual Adjoint Solution Method for Turbulent Flows on Unstructured Grids

An algorithm for solving the discrete adjoint system based on an unstructured-grid discretization of the Navier-Stokes equations is presented. The method is constructed such that an adjoint solution exactly dual to a direct differentiation approach is recovered at each time step, yielding a convergence rate which is asymptotically equivalent to that of the primal system. The new approach is implemented within a three-dimensional unstructured-grid framework and results are presented for inviscid, laminar, and turbulent flows. Improvements to the baseline solution algorithm, such as line-implicit relaxation and a tight coupling of the turbulence model, are also presented. By storing nearest-neighbor terms in the residual computation, the dual scheme is computationally efficient, while requiring twice the memory of the flow solution. The scheme is expected to have a broad impact on computational problems related to design optimization as well as error estimation and grid adaptation efforts.

Nielsen, Eric J.↗