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 397 records · Page 22

Mean excitation energies for stopping powers in various materials using local plasma oscillator strengths

The basic model of Lindhard and Scharff, known as the local plasma model, is used to study the effects on stopping power of the chemical and physical state of the medium. Unlike previous work with the local plasma model, in which individual electron shifts in the plasma frequency were estimated empirically, he Pines correction derived for a degenerate Fermi gas is shown herein to provide a reasonable estimate, even on the atomic scale. Thus, the model is moved to a complete theoretical base requiring no empirical adjustments, as characteristic of past applications. The principal remaining error is in the overestimation of the low-energy absorption properties that are characteristic of the plasma model in the region of the atomic discrete spectrum, although higher-energy phenomena are accurately represented, and even excitation-to-ionization ratios are given to fair accuracy. Mean excitation energies for covalent-bonded gases and solids, for ionic gases and crystals, and for metals are calculated using first-order models of the bonded states.

Wilson, J. W.↗

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗

Precision Oxygen Isotope Measurements of Two C-Rich Hydrated Interplanetary Dust Particles

Introduction: Chondritic-smooth IDPs (Interplanetary Dust Particles) are low porosity objects whose mineralogy is dominated by aqueous alteration products such as Mg-rich phyllosilicates (smectite and serpentine group) and Mg-Fe carbonate minerals. Their hydrated mineralogy combined with low atmospheric entry velocities have been used to infer an origin largely from asteroidal sources. Spectroscopic studies show that the types and abundance of organic matter in CS IDPs is similar to that in CP IDPs. Although CS IDPs show broad similarities to primitive carbonaceous chondrites, only a few particles have been directly linked to specific meteorite groups such as CM and CI chondrites based on the presence of diagnostic minerals. Many CS IDPs however, have carbon contents that greatly exceed that of known meteorite groups suggesting that they either may derive from comets or represent samples of more primitive parent bodies than do meteorites. It is now recognized that many large, dark primitive asteroids in the outer main belt, as well as some trans-Neptunian objects, show spectroscopic evidence for aqueous alteration products on their surfaces. Some CS IDPs exhibit large bulk D enrichments similar to those observed in the cometary CP IDPs. While hydrated minerals in comets have not been unambiguously identified to date, the presence of the smectite group mineral nontronite has been inferred from infrared spectra obtained from the ejecta from comet 9P/Tempel 1 during the Deep Impact mission. Recent observations of low temperature sulfide minerals in Stardust mission samples suggest that limited aqueous activity occurred on comet Wild-2. All of these observations, taken together, suggest that the high-carbon hydrated IDPs are abundant and important samples of primitive solar system objects not represented in meteorite collections. Oxygen isotopic compositions of chondrites reflect mixing between a 16O-rich reservoir and a 17O,18O-rich reservoir produced via mass-independent fractionation. The composition of the 16O-rich reservoir is well constrained but material representing the 17O,18O-rich end-member is rare. Self-shielding models predict that cometary water, presumed to represent this reservoir, should be enriched in 17O and 18O by greater than 200 per mille. The high-carbon hydrated IDPs may be among the best materials available to search for preserved "cometary" H2O signatures. In order to better understand the origin and evolution of these particles, we have obtained 10 hydrated interplanetary dust particles for coordinated mineralogical, isotopic and organic analyses. We have previously reported the results of mineralogical and O isotopic measurements of two hydrated IDPs; here we present results of O isotopic measurements of three additional IDPs. Samples and Methods: Three interplanetary dust particles (L2079C35, L2083D46 and L2083E46) were embedded in S and partially ultramicrotomed into approximately 70 nanometer sections for analysis via transmission electron microscopy (TEM). The remainders of the unsliced particles were removed from S and pressed into high purity Au foil that was cleaned with HF acid and annealed at 800 degrees Centigrade. The pressed IDPs were analyzed via electron microprobe analysis (EPMA) for quantitative bulk chemical analysis. After EPMA analysis, the IDPs were subjected to precision O isotope analysis with the UCLA Cameca IMS-1270 ion probe. A 20 kiloelectronvolt, 0.5 nanoangstrom Cs+ primary beam of approximately 15 micrometers diameter was used for each measurement. Small particles of San Carlos olivine and Burma spinel were pressed into the Au foil for use as standards to correct for instrumental mass fractionation. The detection system was configured for multicollection, with 16O measured on a Faraday cup, and 17O and 18O measured on electron multipliers (EMs). Individual analyses consisted of 15 cycles of 10 seconds per cycle. Additionally, two microtome thin sections were measured for H isotopic compositions with the JSC NanoSIMS 50L ion probe. An 8 picoangstrom, 16 kiloelectronvolt Cs plus primary beam was used. Measurements consisted of H, C, 12C, 16O, and 18O collected with EMs in multicollection. Terrestrial biotite and kerogen were used for isotopic standards. A significant challenge in O isotope measurement of hydrated minerals is the interference from 16OH at mass 17O. We ensured that the 17O and 16OH peaks were fully resolved by using a mass resolution of greater than 7000 and by careful analyses of San Carlos olivine, Burma spinel and chlorite hydrated mineral standards. The hydride was further suppressed with a cold finger attached to an LN2 dewar to trap volatiles in the sample chamber. All sputtered ions were counted (i.e. presputtering was not used); after applying background, yield and deadtime corrections, we performed a change-point analysis on our data via R in order to determine when the sample reached sputtering equilibrium; data points collected prior to the change point were excluded. Change-point analysis was also used to determine whether the IDP had completely sputtered. Results: Mineralogy. IDPs C35 and E46 exhibited hydrated mineralogies, Fe-Ni sulfide grains, nanoglobules and occasional enstatite grains distributed throughout a fine-grained Mg-Fe saponite matrix. C35 also contained breunnerite (Mg,Fe)CO3; solar flare tracks were observed in enstatite, indicating minimal atmospheric entry heating. The mineralogy of D46 is dominated by a large FeS grain with a minor component of adhering silicate material. D46 was strongly heated during atmorpheric entry as evidenced by a well-developed magnetite rim. EPMA analyses show that both C35 and E46 have high carbon contents of 20 weight percentage (approximately 6X CI). D46 contains approximately 6 wt.weight percentage C. A significant amount of the carbon is present as carbon nanoglobules. Results: Hydrogen isotopes: Although the bulk delta D values of both sections of L2079C35 were within error of SMOW (minus 33 plus or minus 19 per mille, 1 plus or minus 14 per mille 1 sigma), several delta D-rich hotspots were also identified, reaching 2000 per mille. As shown in Fig. 1, these hotspots are clearly associated with discrete carbonaceous inclusions that are akin to nanoglobules found in many meteorites and other IDPs. Oxygen Isotopes. Results of the oxygen isotope measurements are shown in Figure 1. The oxygen isotope composition for L2079C35 was delta 18O equals plus11.6 plus or minus 1.9per mille, delta 17O equals plus 7.9 plus or minus 1.9per mille (2 standard errors). The oxygen isotope composition for L2083D46 was delta 18O equals minus 8.1 plus or minus 1.9 per mille, delta 17O equals minus 6.4 plus or minus 3.1 per mille (2 standard errors). The oxygen isotope composition for L2083E46 was delta 18O equals plus 12.0 plus or minus 1.9 per mille, delta 17O equals plus 9.2 plus or minus 2.0 per mille (2 standard errors). Discussion: Despite mineralogical similarities to highly aqueously altered carbonaceous chondrites, the hydrated IDPs we analyzed have oxygen isotopic compositions that are distinct from matrix materials in the CI, CM, and CR chondrites. The IDPs plot along the Young-Russell line, with delta 17O values for C35 and E46 suggestive of interaction with a 16O-poor reservoir. However, we have thus far not observed evidence of extreme 16O-poor reservoirs expected from self-shielding models and observed in Acfer 094 simplectite. The high carbon contents of the IDPs also set them apart from known meteoritic samples. The lack of atmospheric entry heating effects are consistent with low encounter velocities and suggest either an asteroidal source, or a low inclination, low eccentricity cometary origin. Conclusions: The unusual oxygen isotopic compositions, high carbon contents, and the abundance of Drich nanoglobules, together, suggest that the high-carbon, hydrated IDPs are derived from a primitive source that is not yet represented in meteorite collections.

Snead, C. J.↗

Wind estimates from cloud motions: Phase 1 of an in situ aircraft verification experiment

An initial experiment was conducted to verify geostationary satellite derived cloud motion wind estimates with in situ aircraft wind velocity measurements. Case histories of one-half hour to two hours were obtained for 3-10km diameter cumulus cloud systems on 6 days. Also, one cirrus cloud case was obtained. In most cases the clouds were discrete enough that both the cloud motion and the ambient wind could be measured with the same aircraft Inertial Navigation System (INS). Since the INS drift error is the same for both the cloud motion and wind measurements, the drift error subtracts out of the relative motion determinations. The magnitude of the vector difference between the cloud motion and the ambient wind at the cloud base averaged 1.2 m/sec. The wind vector at higher levels in the cloud layer differed by about 3 m/sec to 5 m/sec from the cloud motion vector.

Hasler, A. F.↗

Wind estimates from cloud motions - Phase 1 of an in situ aircraft verification experiment

An initial experiment has been conducted to verify geostationary-satellite-derived cloud motion wind estimates with in situ aircraft wind velocity measurements. Case histories of 1/2 to 2 h were obtained for 3-10 km diameter cumulus cloud systems on 6 days. Also, one cirrus cloud case was obtained. In most cases the clouds were discrete enough that both the cloud motion and the ambient wind could be measured with the same aircraft Inertial Navigation System (INS). Since the INS drift error is the same for both the cloud motion and wind measurements, the drift error drops out of the relative motion determinations. The magnitude of the vector difference between the cloud motion and the ambient wind at the cloud base averaged 1.2 m/sec. The wind vector at higher levels in the cloud layer differed by about 3 to 5 m/sec from the cloud motion vector.

Hasler, A. F.↗

VLSI Reed-Solomon Encoder

Modular Reed-solomon encoder uses identical custom VLSI chips called "symbol slices." By cascading and properly interconnecting group of these chips, encoder is made for any desired error-correcting capability and interleaving level. VLSI encoder requires only one-tenth the number of chips required by conventional Reed-Solomon Circuit implemented with discrete IC's.

Liu, K. Y.↗

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Harten, A.↗

Prediction of helicopter rotor discrete frequency noise: A computer program incorporating realistic blade motions and advanced acoustic formulation

A computer program has been developed at the Langley Research Center to predict the discrete frequency noise of conventional and advanced helicopter rotors. The program, called WOPWOP, uses the most advanced subsonic formulation of Farassat that is less sensitive to errors and is valid for nearly all helicopter rotor geometries and flight conditions. A brief derivation of the acoustic formulation is presented along with a discussion of the numerical implementation of the formulation. The computer program uses realistic helicopter blade motion and aerodynamic loadings, input by the user, for noise calculation in the time domain. A detailed definition of all the input variables, default values, and output data is included. A comparison with experimental data shows good agreement between prediction and experiment; however, accurate aerodynamic loading is needed.

Brentner, K. S.↗

Improved numerical methods for turbulent viscous recirculating flows

The hybrid-upwind finite difference schemes employed in generally available combustor codes possess excessive numerical diffusion errors which preclude accurate quantative calculations. The present study has as its primary objective the identification and assessment of an improved solution algorithm as well as discretization schemes applicable to analysis of turbulent viscous recirculating flows. The assessment is carried out primarily in two dimensional/axisymetric geometries with a view to identifying an appropriate technique to be incorporated in a three-dimensional code.

Turan, A.↗

Digital visual communications using a Perceptual Components Architecture

The next era of space exploration will generate extraordinary volumes of image data, and management of this image data is beyond current technical capabilities. We propose a strategy for coding visual information that exploits the known properties of early human vision. This Perceptual Components Architecture codes images and image sequences in terms of discrete samples from limited bands of color, spatial frequency, orientation, and temporal frequency. This spatiotemporal pyramid offers efficiency (low bit rate), variable resolution, device independence, error-tolerance, and extensibility.

Watson, Andrew B.↗

Visually Lossless Data Compression for Real-Time Frame/Pushbroom Space Science Imagers

A visually lossless data compression technique is currently being developed for space science applications under the requirement of high-speed push-broom scanning. The technique is also applicable to frame based imaging and is error-resilient in that error propagation is contained within a few scan lines. The algorithm is based on a block transform of a hybrid of modulated lapped transform (MLT) and discrete cosine transform (DCT), or a 2-dimensional lapped transform, followed by bit-plane encoding; this combination results in an embedded bit string with exactly the desirable compression rate as desired by the user. The approach requires no unique table to maximize its performance. The compression scheme performs well on a suite of test images typical of images from spacecraft instruments. Flight qualified hardware implementations are in development; a functional chip set is expected by the end of 2001. The chip set is being designed to compress data in excess of 20 Msamples/sec and support quantizations from 2 to 16 bits.

Yeh, Pen-Shu↗

Image Correlation Pattern Optimization for Micro-Scale In-Situ Strain Measurements

The accuracy and precision of digital image correlation (DIC) is a function of three primary ingredients: image acquisition, image analysis, and the subject of the image. Development of the first two (i.e. image acquisition techniques and image correlation algorithms) has led to widespread use of DIC; however, fewer developments have been focused on the third ingredient. Typically, subjects of DIC images are mechanical specimens with either a natural surface pattern or a pattern applied to the surface. Research in the area of DIC patterns has primarily been aimed at identifying which surface patterns are best suited for DIC, by comparing patterns to each other. Because the easiest and most widespread methods of applying patterns have a high degree of randomness associated with them (e.g., airbrush, spray paint, particle decoration, etc.), less effort has been spent on exact construction of ideal patterns. With the development of patterning techniques such as microstamping and lithography, patterns can be applied to a specimen pixel by pixel from a patterned image. In these cases, especially because the patterns are reused many times, an optimal pattern is sought such that error introduced into DIC from the pattern is minimized. DIC consists of tracking the motion of an array of nodes from a reference image to a deformed image. Every pixel in the images has an associated intensity (grayscale) value, with discretization depending on the bit depth of the image. Because individual pixel matching by intensity value yields a non-unique scale-dependent problem, subsets around each node are used for identification. A correlation criteria is used to find the best match of a particular subset of a reference image within a deformed image. The reader is referred to references for enumerations of typical correlation criteria. As illustrated by Schreier and Sutton and Lu and Cary systematic errors can be introduced by representing the underlying deformation with under-matched shape functions. An important implication, as discussed by Sutton et al., is that in the presence of highly localized deformations (e.g., crack fronts), error can be reduced by minimizing the subset size. In other words, smaller subsets allow the more accurate resolution of localized deformations. Contrarily, the choice of optimal subset size has been widely studied and a general consensus is that larger subsets with more information content are less prone to random error. Thus, an optimal subset size balances the systematic error from under matched deformations with random error from measurement noise. The alternative approach pursued in the current work is to choose a small subset size and optimize the information content within (i.e., optimizing an applied DIC pattern), rather than finding an optimal subset size. In the literature, many pattern quality metrics have been proposed, e.g., sum of square intensity gradient (SSSIG), mean subset fluctuation, gray level co-occurrence, autocorrelation-based metrics, and speckle-based metrics. The majority of these metrics were developed to quantify the quality of common pseudo-random patterns after they have been applied, and were not created with the intent of pattern generation. As such, it is found that none of the metrics examined in this study are fit to be the objective function of a pattern generation optimization. In some cases, such as with speckle-based metrics, application to pixel by pixel patterns is ill-conditioned and requires somewhat arbitrary extensions. In other cases, such as with the SSSIG, it is shown that trivial solutions exist for the optimum of the metric which are ill-suited for DIC (such as a checkerboard pattern). In the current work, a multi-metric optimization method is proposed whereby quality is viewed as a combination of individual quality metrics. Specifically, SSSIG and two auto-correlation metrics are used which have generally competitive objectives. Thus, each metric could be viewed as a constraint imposed upon the others, thereby precluding the achievement of their trivial solutions. In this way, optimization produces a pattern which balances the benefits of multiple quality metrics. The resulting pattern, along with randomly generated patterns, is subjected to numerical deformations and analyzed with DIC software. The optimal pattern is shown to outperform randomly generated patterns.

Bomarito, G. F.↗

Digital simulation of continuous error models with application to an instrument landing system error

A digital simulation of the continuous error of the localized beam of a conventional instrument landing system is discussed. The digital simulation was developed during the analysis of space shuttle navigation capabilities. A discrete mathematical model for use on a digital computer is described. The model generates an output random sequence which is equivalent, for simulation purposes, to the desired random process. The model is a system of difference equations driven by a zero-mean Gaussian random sequence.

Merrick, R. B.↗

New method for propagating the square root covariance matrix in triangular form

A method proposed by Potter of applying the Kalman-Bucy filter to the problem of estimating the state of a dynamic system is described, in which the square root of the state error covariance matrix is used to process the observations. A new technique which propagates the covariance square root matrix in lower triangular form is given for the discrete observation case. The technique is faster than previously proposed algorithms and is well-adapted for use with the Carlson square root measurement algorithm.

Choe, C. Y.↗

A separation theorem for the stochastic sampled-data LQG problem

This paper considers the control of a continuous linear plant disturbed by white plant noise when the control is constrained to be a piecewise constant function of time; i.e. a stochastic sampled-data system. The cost function is the integral of quadratic error terms in the state and control, thus penalizing errors at every instant of time while the plant noise disturbs the system continuously. The problem is solved by reducing the constrained continuous problem to an unconstrained discrete one. It is shown that the separation principle for estimation and control still holds for this problem when the plant disturbance and measurement noise are Gaussian.

Halyo, N.↗

Dynamics Questions Associated with the AERCam Sprint Free-Flyer

The International Space Station will require the development of small robotic vehicles for such tasks as external inspection, monitoring of extravehicular activities (EVA's) and station build-up, and providing additional lighting at EVA worksites. The Autonomous EVA Robotic Camera (AERCam) family of free-flyers is currently being developed at NASA Johnson Space Center to perform these functions; the first member of the family is the AERCam Sprint vehicle. Many interesting dynamical questions are associated with the Sprint free-flyer. For instance, the reaction of a vehicle which is nearly spherically symmetric (such as Sprint) to a stuck-on thruster is significantly more complicated than that obtained for an idealized, perfectly spherical, spacecraft model. In particular, the real spacecraft will experience a form of forced nutation, with convergence towards either its major or minor principal axis, depending on both the applied torque and the mass properties of the vehicle. Furthermore, the body-fixed jet force vector may have a significant component along this principal axis, so giving rise to a considerable net linear acceleration of the spacecraft. The large velocity that can result is very important, as it may lead to collision with the nearby Orbiter, and is completely overlooked in the idealized analysis. This report will firstly briefly describe the stuck-on thruster dynamics of the real vehicle, and outline how the small products of inertia of the spacecraft determine the time constants of the motion. Secondly, the dynamical effects of a failed-off jet on the Sprint free-flyer will be described in more detail, and compared with the stuck-on thruster case. This will help to show whether the two malfunctions should be dealt with differently in flight. Finally, the stuck-on thruster detection software (known as the uncommanded motion algorithm) that is proposed to be flown on the Sprint vehicle will be analyzed, and all possible perturbation sources that may tend to give rise to false stuck-on thruster alarms quantified. It will be shown that false alarms can be triggered by discrete events such as a failed-off thruster, a glancing collision with Orbiter structure, or thruster saturation. They can also be triggered by a combination of errors introduced by constantly present sources such as thrust level errors, thruster misalignment, inertial cross-coupling, angular rate sensor noise, and structural vibrations induced by thruster firings. Based on this fact, the original plan to have the spacecraft enter safe mode whenever the uncommanded motion algorithm is triggered does not appear to be advisable. A better approach seems to be to provide the Sprint pilot with a warning whenever uncommanded motion is detected, and then allow the crew to determine whether safing the vehicle is appropriate.

Williams, Trevor↗

Flight Deck Surface Trajectory-Based Operations

Surface Trajectory-Based Operations (STBO) is a future concept for surface operations where time requirements are incorporated into taxi operations to support surface planning and coordination. Pilot-in-the-loop flight deck simulations have been conducted to study flight deck displays algorithms to aid pilots in complying with the time requirements of time-based taxi operations (i.e., at discrete locations in 3 12 D operations or at all points along the route in 4DT operations). The results of these studies (conformance, time-of-arrival error, eye-tracking data, and safety ratings) are presented. Flight deck simulation work done in collaboration with DLR is described. Flight deck research issues in future auto-taxi operations are also introduced.

flight deck↗

Performance Assessment of the Eastern Range Mosquito Lagoon 915-MHz Doppler Radar Wind Profiler

Please find enclosed the performance assessment for the Mosquito Lagoon 915-MHz Doppler Radar Wind Profiler (DRWP). This assessment was conducted following the installation of the new DRWP systems at the Eastern Range. DRWP profiles were compared to Automated Meteorological Profiling System balloons launched from Cape Canaveral Space Force Station. Launching balloons from the DRWP site was discussed, but was determined to be unfeasible. Therefore, these results include errors due to the spatial separation of winds. Mean wind component deltas generally ranged between -1.4-2.0 m/s between the balloons and DRWP systems. Data availability at discrete altitudes ranged from 7% to 98% with the least amount of available profiles being in the upper altitude ranges between 4000 m and 5000 m as well as below 300 m. The effective vertical resolution was not able to be estimated in this study, likely due to the relatively small dataset. The results from this analysis were requested by the United States Space Force for use of the acceptance of the DRWP system and for internal EV44 reference.

meteorology↗