Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical errors”

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 181 records · Page 10

Chemical Source Inversion using Assimilated Constituent Observations in an Idealized Two-dimensional System

We present a source inversion technique for chemical constituents that uses assimilated constituent observations rather than directly using the observations. The method is tested with a simple model problem, which is a two-dimensional Fourier-Galerkin transport model combined with a Kalman filter for data assimilation. Inversion is carried out using a Green's function method and observations are simulated from a true state with added Gaussian noise. The forecast state uses the same spectral spectral model, but differs by an unbiased Gaussian model error, and emissions models with constant errors. The numerical experiments employ both simulated in situ and satellite observation networks. Source inversion was carried out by either direct use of synthetically generated observations with added noise, or by first assimilating the observations and using the analyses to extract observations. We have conducted 20 identical twin experiments for each set of source and observation configurations, and find that in the limiting cases of a very few localized observations, or an extremely large observation network there is little advantage to carrying out assimilation first. However, in intermediate observation densities, there decreases in source inversion error standard deviation using the Kalman filter algorithm followed by Green's function inversion by 50% to 95%.

Tangborn, Andrew↗

Calculation of confined swirling flows

Swirling jets are an important constituent flow of many types of combustion equipment. Flame size, shape, stability and combustion intensity are all favorably influenced by swirling some portion of the incoming air with the fuel. Fluid dynamics computer codes are being developed for the mathematical simulation of flows in practical combustors. It is important that these codes have the ability to accurately and reliably calculate swirling flows. Turbulent flow calculations with swirl can be inaccurate due to at least three major reasons: limitations of the turbulence modeling used, incorrect or inadequate specifications of inlet boundary conditions, and, error introduced through the numerics. These sources of error are described, and examples provided of each. The state-of-the-art in such calculations is reviewed. The controlling character of current numerics is demonstrated, and it is concluded that as a first priority, better numerics must be arrived at, and that an improved understanding of the discharge flow from swirl generators is essential as a second priority. When these priorities are satisfied an improved turbulence model is desirable. It is also concluded that an improved understanding of swirl-induced recirculation would be an asset.

Sturgess, G. J.↗

Numerical computation of transient coaxial entry tube flows

A numerical program was developed to compute transient laminar flows in two dimensions including multicomponent mixing and chemical reaction. The program can compute both incompressible flows and compressible flows at all speeds, and it is applied to describe transient and steady state solutions for low subsonic, coaxial entry, tue flows. Single component, nonreacting flows comprise most of the solutions, but one steady state solution is presented for trace concentration constituents engaging in a second order reaction. Numerical stability was obtained by adding at each calculation point a correction for numerical diffusion errors caused by truncation of the Taylor series used to finite difference the conservation equations. Transient computations were made for fluids initially at rest, then subjected to step velocity inputs that were uniform across each region of the entry plane and were held constant throughout the computation period. For center tube to annulus velocity ratios of 0.5 and 2.0, the bulk fluid in the tube initially moved in plug flow, but strong radial flows developed near the injection plane which moved the fluid into the high shear region between the jets and away from the tube wall.

Wieber, P. R.↗

Evaluation of SMART sensor displays for multidimensional precision control of Space Shuttle remote manipulator

An enhanced proximity sensor and display system was developed at the Jet Propulsion Laboratory (JPL) and tested on the full scale Space Shuttle Remote Manipulator at the Johnson Space Center (JSC) Manipulator Development Facility (MDF). The sensor system, integrated with a four-claw end effector, measures range error up to 6 inches, and pitch and yaw alignment errors within + or 15 deg., and displays error data on both graphic and numeric displays. The errors are referenced to the end effector control axes through appropriate data processing by a dedicated microcomputer acting on the sensor data in real time. Both display boxes contain a green lamp which indicates whether the combination of range, pitch and yaw errors will assure a successful grapple. More than 200 test runs were completed in early 1980 by three operators at JSC for grasping static and capturing slowly moving targets. The tests have indicated that the use of graphic/numeric displays of proximity sensor information improves precision control of grasp/capture range by more than a factor of two for both static and dynamic grapple conditions.

Bejczy, A. K.↗

Linear discriminant analysis with misallocation in training samples

Linear discriminant analysis for a two-class case is studied in the presence of misallocation in training samples. A general appraoch to modeling of mislocation is formulated, and the mean vectors and covariance matrices of the mixture distributions are derived. The asymptotic distribution of the discriminant boundary is obtained and the asymptotic first two moments of the two types of error rate given. Certain numerical results for the error rates are presented by considering the random and two non-random misallocation models. It is shown that when the allocation procedure for training samples is objectively formulated, the effect of misallocation on the error rates of the Bayes linear discriminant rule can almost be eliminated. If, however, this is not possible, the use of Fisher rule may be preferred over the Bayes rule.

Chhikara, R.↗

Spline-based distributed system identification with application to large space antennas

A parameter and state estimation technique for distributed models is demonstrated through the solution of a problem generic to large space antenna system identification. Assuming the position of the reflective surface of the maypole (hoop/column) antenna to be approximated by the static two-dimensional, stretched-membrane partial differential equation with variable-stiffness coefficient functions, a spline-based approximation procedure is described that estimates the shape and stiffness functions from data set observations. For given stiffness functions, the Galerkin projection with linear spline-based functions is applied to project the distributed problem onto a finite-dimensional subspace wherein algebraic equations exist for determining a static shape (state) prediction. The stiffness functions are then parameterized by cubic splines and the parameters estimated by an output error technique. Numerical results are presented for data descriptive of a 100-m-diameter maypole antenna.

Banks, H. T.↗

Radial orbit error reduction and sea surface topography determination using satellite altimetry

A method is presented in satellite altimetry that attempts to simultaneously determine the geoid and sea surface topography with minimum wavelengths of about 500 km and to reduce the radial orbit error caused by geopotential errors. The modeling of the radial orbit error is made using the linearized Lagrangian perturbation theory. Secular and second order effects are also included. After a rather extensive validation of the linearized equations, alternative expressions of the radial orbit error are derived. Numerical estimates for the radial orbit error and geoid undulation error are computed using the differences of two geopotential models as potential coefficient errors, for a SEASAT orbit. To provide statistical estimates of the radial distances and the geoid, a covariance propagation is made based on the full geopotential covariance. Accuracy estimates for the SEASAT orbits are given which agree quite well with already published results. Observation equations are develped using sea surface heights and crossover discrepancies as observables. A minimum variance solution with prior information provides estimates of parameters representing the sea surface topography and corrections to the gravity field that is used for the orbit generation. The simulation results show that the method can be used to effectively reduce the radial orbit error and recover the sea surface topography.

Engelis, Theodossios↗

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.↗

Interface procedures for overlapping grids

The solution of partial differential equations by means of either implicit or conservative explicit finite-difference methods, on the one hand, or multicomponent overlapping grid schemes on the other, is presently subjected to interpolation at grid boundaries in an attempt to eliminate the time lag. The accuracy of the transient solution of a hyperbolic or parabolic partial differential equation is dependent on the procedures used to transfer information between grids in a composite grid system. The present techniques reduce numerical solution errors.

Mastin, C. Wayne↗

Geometric bias and time coincidence in 3-dimensional laser Doppler velocimeter systems

The measurement by a three-dimensional laser Doppler velocimeter of a turbulent flow has been numerically simulated. Errors associated with the probe volume geometry and the coincidence time window concept are revealed. One type of error occurs for high system data rates when multiple particles lead to system realizations. Another error occurs associated with a geometric bias discovered in the present study. This three-dimensional ldv geometric bias exists even for single-particle realizations and regardless of the system data rate. A technique for the elimination of the geometric bias is presented.

Brown, J. L.↗

Calculation of turbulent three-dimensional jet-induced flow in a rectangular enclosure

The application of a flux-spline scheme to three-dimensional fluid flow is reported. A comparison is made of the performance of this scheme with that of the Power-law differencing scheme. The numerical results are compared with experimental data. For the problem considered in this study, the flux-spline scheme is more accurate than the Power-law. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential of providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid.

Nikjooy, M.↗

The application of Encke's method to long arc orbit determination solutions

The Laser Geodynamics Satellite (LAGEOS) was launched on May 4, 1976 to provide geophysical measurements by means of laser ranging techniques. To date, over twelve years of laser range measurements have been collected from various tracking stations located around the world. Laser range measurements to LAGEOS have contributed to studies of earth rotation, plate tectonics, global baseline, and the gravity field as well as many other areas. Some of these studies are based upon the determination of a single, continuous orbit for LAGEOS for time spans on the order of several years. Current studies at the University of Texas Center for Space Research include the precision orbit determination of LAGEOS for arc lengths of up to 12.8 years which represents over 31,000 orbital revolutions. These long arc studies have led to the implementation of Encke's method to improve the convergence of the batch filter while reducing numerical integration errors. While the technique has been successfully applied to arc lengths of up to 12.8 years, the results presented focus on the solution of a six-year arc.

Lundberg, J. B.↗

Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms

The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.

Lee, Jeh Won↗

Improving the numerical integration solution of satellite orbits in the presence of solar radiation pressure using modified back differences

The method of modified back differences, a technique that significantly reduces the numerical integration errors associated with crossing shadow boundaries using a fixed-mesh multistep integrator without a significant increase in computer run time, is presented. While Hubbard's integral approach can produce significant improvements to the trajectory solution, the interpolation method provides the best overall results. It is demonstrated that iterating on the point mass term correction is also important for achieving the best overall results. It is also shown that the method of modified back differences can be implemented with only a small increase in execution time.

Lundberg, J. B.↗

Ellipsometric study of metal-organic chemically vapor deposited III-V semiconductor structures

An ellipsometric study of MOCVD-grown layers of AlGaAs and InGaAs in thick films and strained layer complex structures is presented. It is concluded that the ternary composition of thick nonstrained layers can be accurately determined to within experimental errors using numerical algorithms. In the case of complex structures, thickness of all layers and the alloy composition of nonstrained layers can be determined simultaneously, provided that the correlations between parameters is no higher than 0.9.

Alterovitz, Samuel A.↗

An improved far field drag calculation method for nonlinear CFD codes

An improved method to calculate drag based on far field conditions is presented and demonstrated. The method is illustrated using two CFD codes, PARC and CFL3D, and examining their ability to preserve lift and drag coefficients along grid line defined integration paths. The flow fields were generated by solving the Euler equations for a NACA 0012 airfoil at a free stream Mach number of 0.8 and angle of attack of 1.25 degrees. In comparison to force coefficients obtained by surface pressure integration, neither code acceptably preserved both force coefficients throughout the near and far fields. An investigation into the relationship between numerical prediction error and calculated force coefficients revealed a direct connection between solution mass conservation error and force coefficients error. A method to correct the predicted force coefficients based on integrated mass conservation error is described and demonstrated. The corrected force coefficients for both codes are shown to be more accurate than the uncorrected values. The correction method is applicable to both two and three dimensions and is independent of the algorithm used to generate the flow field.

Sirbaugh, James R.↗

Dynamic analysis and control of lightweight manipulators with flexible parallel link mechanisms

The flexible parallel link mechanism is designed for increased rigidity to sustain the buckling when it carries a heavy payload. Compared to a one link flexible manipulator, a two link flexible manipulator, especially the flexible parallel mechanism, has more complicated characteristics in dynamics and control. The objective of this research is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model. The step response of the analytical model and the TREETOPS model match each other well. The nonlinear dynamics is studied using a sinusoidal excitation. The actuator dynamic effect on a flexible robot was investigated. The effects are explained by the root loci and the Bode plot theoretically and experimentally. For the base performance for the advanced control scheme, a simple decoupled feedback scheme is applied.

Lee, Jeh Won↗