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 487 records · Page 27

Error analysis of finite element solutions for postbuckled cylinders

A general method of error analysis and correction is investigated for the discrete finite-element results for cylindrical shell structures. The method for error analysis is an adaptation of the method of successive approximation. When applied to the equilibrium equations of shell theory, successive approximations derive an approximate continuous solution from the discrete finite-element results. The advantage of this continuous solution is that it contains continuous partial derivatives of an order higher than the basis functions of the finite-element solution. Preliminary numerical results are presented in this paper for the error analysis of finite-element results for a postbuckled stiffened cylindrical panel modeled by a general purpose shell code. Numerical results from the method have previously been reported for postbuckled stiffened plates. A procedure for correcting the continuous approximate solution by Newton's method is outlined.

Sistla, Rajaram↗

Error analysis for reducing noisy wide-gap concentric cylinder rheometric data for nonlinear fluids - Theory and applications

This work discusses the propagation of errors for the recovery of the shear rate from wide-gap concentric cylinder viscometric measurements of non-Newtonian fluids. A least-square regression of stress on angular velocity data to a system of arbitrary functions is used to propagate the errors for the series solution to the viscometric flow developed by Krieger and Elrod (1953) and Pawlowski (1953) ('power-law' approximation) and for the first term of the series developed by Krieger (1968). A numerical experiment shows that, for measurements affected by significant errors, the first term of the Krieger-Elrod-Pawlowski series ('infinite radius' approximation) and the power-law approximation may recover the shear rate with equal accuracy as the full Krieger-Elrod-Pawlowski solution. An experiment on a clay slurry indicates that the clay has a larger yield stress at rest than during shearing, and that, for the range of shear rates investigated, a four-parameter constitutive equation approximates reasonably well its rheology. The error analysis presented is useful for studying the rheology of fluids such as particle suspensions, slurries, foams, and magma.

Borgia, Andrea↗

Error Estimate of the Ares I Vehicle Longitudinal Aerodynamic Characteristics Based on Turbulent Navier-Stokes Analysis

Numerical predictions of the longitudinal aerodynamic characteristics for the Ares I class of vehicles, along with the associated error estimate derived from an iterative convergence grid refinement, are presented. Computational results are based on the unstructured grid, Reynolds-averaged Navier-Stokes flow solver USM3D, with an assumption that the flow is fully turbulent over the entire vehicle. This effort was designed to complement the prior computational activities conducted over the past five years in support of the Ares I Project with the emphasis on the vehicle s last design cycle designated as the A106 configuration. Due to a lack of flight data for this particular design s outer mold line, the initial vehicle s aerodynamic predictions and the associated error estimates were first assessed and validated against the available experimental data at representative wind tunnel flow conditions pertinent to the ascent phase of the trajectory without including any propulsion effects. Subsequently, the established procedures were then applied to obtain the longitudinal aerodynamic predictions at the selected flight flow conditions. Sample computed results and the correlations with the experimental measurements are presented. In addition, the present analysis includes the relevant data to highlight the balance between the prediction accuracy against the grid size and, thus, the corresponding computer resource requirements for the computations at both wind tunnel and flight flow conditions. NOTE: Some details have been removed from selected plots and figures in compliance with the sensitive but unclassified (SBU) restrictions. However, the content still conveys the merits of the technical approach and the relevant results.

Aerodynamic characteristics↗

Reducing Errors in Velocity–Azimuth Display (VAD) Wind and Deformation Retrievals from Airborne Doppler Radars in Convective Environments

The present study describes methods to reduce the uncertainty of velocity–azimuth display (VAD) wind and deformation retrievals from downward-pointing, conically scanning, airborne Doppler radars. These retrievals have important applications in data assimilation and real-time data processing. Several error sources for VAD retrievals are considered here, including violations to the underlying wind field assumptions, Doppler velocity noise, data gaps, temporal variability, and the spatial weighting function of the VAD retrieval. Specific to airborne VAD retrievals, we also consider errors produced due to the radar scans occurring while the instrument platform is in motion. While VAD retrievals are typically performed using data from a single antenna revolution, other strategies for selecting data can be used to reduce retrieval errors. Four such data selection strategies for airborne VAD retrievals are evaluated here with respect to their effects on the errors. These methods are evaluated using the second hurricane nature run numerical simulation, analytic wind fields, and observed Doppler radar radial velocities. The proposed methods are shown to reduce the median absolute error of the VAD wind retrievals, especially in the vicinity of deep convection embedded in stratiform precipitation. The median absolute error due to wind field assumption violations for the along-track and for the across-track wind is reduced from 0.36 to 0.08 m s−1 and from 0.35 to 0.24 m s−1, respectively. Although the study focuses on Doppler radars, the results are equally applicable to conically scanning Doppler lidars as well.

Charles N Helms↗

Numerical analysis of discrete fractional integrodifferential structural dampers

This paper develops solution algorithms enabling the handling of the dynamic response of nonlinear structures contained discretely attached dampers modeled by fractional integrodifferential operators of the Grunwald-Liouville-Riemann type. The development consists of two levels of formulation, namely: (1) numerical approximations of fractional operators and, (2) the establishment of global level implicit schemes enabling the solution to nonlinear structural formulations. To generalize the overall results, error estimates are derived for the fractional operator approximation algorithm. These enable an ongoing optimization of solution efficiency for a given error tolerance. To benchmark the scheme, the results of several numerical experiments are presented. These illustrate the numerical characteristics of the overall formulation.

Padovan, Joseph↗

Optimum employment of satellite indirect soundings as numerical model input

The characteristics of satellite-derived temperature soundings that would significantly affect their use as input for numerical weather prediction models were examined. Independent evaluations of satellite soundings were emphasized to better define error characteristics. Results of a Nimbus-6 sounding study reveal an underestimation of the strength of synoptic scale troughs and ridges, and associated gradients in isobaric height and temperature fields. The most significant errors occurred near the Earth's surface and the tropopause. Soundings from the TIROS-N and NOAA-6 satellites were also evaluated. Results again showed an underestimation of upper level trough amplitudes leading to weaker thermal gradient depictions in satellite-only fields. These errors show a definite correlation to the synoptic flow patterns. In a satellite-only analysis used to initialize a numerical model forecast, it was found that these synoptically correlated errors were retained in the forecast sequence.

Horn, L. H.↗

The analysis of star catalogs. 1: An intercomparison among the FK3, the FK4, and the FK5

Progress in the construction of the fundamental catalogs allows one to ask the following questions: 'Now that we have the FK5 with which to measure by, how good was the FK4? The FK3? How realistic were the error estimates they gave for their equatorial coordinates and their proper motions? Have their residual systematic differences continued to decrease in amplitude?' and so forth. Similar questions could be asked of the GC and the N30, and a companion paper will similarly analyze them. Finally, the identical series of questions could be asked of the AGK2, AGK3, AGK3U sequence. In this paper we develop a general method to objectively investigate, post facto, the true quality of star catalog error estimates. We present a complete mathematical formulation of the questions above and delineate the calculations needed to address them. Our initial numerical applications are restricted to a straightforward examination of the mean errors of the coordinates and proper motions in the FK3 (using the FK4 and the FK5 for comparison) and in the FK4 (using the FK5 for comparison). However, we can simply and consistently interpret the results of these computations -- for the FK catalogs fail to follow the predictions of mathematical statistics -- only if the FK3 and the FK4 contain higher angular frequency systematic errors. These errors, at the mean epochs of place of the older catalogs, are larger than the random errors in the FK5 system. We expect that our new ability to discern these features means that we can also prevent them during compilation. Finally, note that the concepts behind our procedures are completely general and can be applied to any set of compiled data -- astrometric, photometric, geodetic, and so forth.

Bucciarelli, B.↗

An hp-adaptivity and error estimation for hyperbolic conservation laws

This paper presents an hp-adaptive discontinuous Galerkin method for linear hyperbolic conservation laws. A priori and a posteriori error estimates are derived in mesh-dependent norms which reflect the dependence of the approximate solution on the element size (h) and the degree (p) of the local polynomial approximation. The a posteriori error estimate, based on the element residual method, provides bounds on the actual global error in the approximate solution. The adaptive strategy is designed to deliver an approximate solution with the specified level of error in three steps. The a posteriori estimate is used to assess the accuracy of a given approximate solution and the a priori estimate is used to predict the mesh refinements and polynomial enrichment needed to deliver the desired solution. Numerical examples demonstrate the reliability of the a posteriori error estimates and the effectiveness of the hp-adaptive strategy.

Bey, Kim S.↗

Distance corrections for single- and dual-color lasers by ray tracing

The physical arc length difference between the red and blue beam of a dual-color laser generates an error term in determining Delta-n as well as n (the refractive index) and D (the line length). Numerical ray trace examples and theoretical approximations show that the resulting relative error Delta-D/D increases as D-squared. Error reduction by a factor of nearly 10 is possible by using the pressure, temperature, and humidity of both endpoints to calculate n. The present study results from a desire to improve the precision of monitoring the local and regional position of the Lunar Laser Ranging Observatory atop Haleakala, Maui, Hawaii; the monitoring of the earth's crustal deformation and strain is a primary concern of the study.

Berg, E.↗

Guidance and trajectory considerations in lunar mass transportation

Flight-mechanics problems associated with large-scale transport of lunar mass to a space colony or manufacturing facility are discussed. The proposed transport method involves launch of payloads from a mass-driver on the lunar surface, onto ballistic trajectories to a passive mass-catcher located near the L2 libration point, with the caught mass subsequently being transported to the colony. Arrival velocities at L2, sensitivities in arrival dispersion due to launch errors, and effects of launch site location are treated, via numerically integrated orbits in the restricted three-body problem. From any launch site it is possible to define a target point reached with zero dispersion due to errors in a selected component of launch velocity. Effects of lunar geometrical librations and of obliquity, as well as the conditions for biasing a trajectory away from L2 so as to reduce stationkeeping costs, are dealt with along with transfer orbits from L2 to the colony. The theory of capture and the theory of resonance lead to a colony orbit, with period approximately two weeks, reached from L2 with velocity increment as low as 9.02 m/sec.

Heppenheimer, T. A.↗

ON THE ACCURACY OF MEASUREMENTS MADE UPON FILMS PHOTOGRAPHED BY BAKER-NUNN SATELLITE TRACKING CAMERAS

The photographs taken by the Baker-Nunn cameras of the Smithsonian Astrophysical Observatory were formerly measured with Van Biesbroeck goniometers, but are now measured exclusively with Mann two-screw comparators. Instead of measuring point images, we usually have to measure oblong images, either trails or breaks (the latter produced by interruption of the exposure). We can therefore expect that the accuracy of the determination of their positions is more influenced by the personal errors than it is in conventional astrometry. Approximately 800 measurements were made on 34 images of different length, and, from these, the frequency distribution of settings was determined. Then the relationship between this distribution and the length of the image (trail or break) was examined, and the relationship between the "magnitude error" and the length of image was determined. Measurements were made at different angles formed by the trails and the reticle lines, and it was found that the accuracy was not affected by the angle of the setting. Position determinations have been made using reference stars at different distances. On the basis of these results, we can conclude that the accuracy is neither influenced by the distortion of the emulsion nor by the optical distortion within an area of a diameter of S cm (5°.8). We believe that we are justified in using the linear plate constant method, at least within this area. From numerous double and multiple measurements, the following standard errors were determined for a single position.

ASTROMETRY↗

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

Control of the errors of discretization and idealization in finite element analysis

Understanding of the basic principles which control errors of discretization in finite element analysis has increased very substantially since 1980. The main milestones were: (1) development of the theoretical basis of p-extensions (1981); (2) understanding of the proper interplay between mesh design and assignment of polynomial degree to elements. Practical realization of exponential convergence rates, independently of the smoothness of the exact solution (1984); and (3) industrial experience with the new finite element technology known as the p- or hp-version of the finite element method: General Dynamics reported thirty- to forty-fold savings in terms of human time and large savings in computer time (1986). Lockheed reported favorably on their evaluation of error estimation and quality control capabilities of the p-version in industrial settings (1987). The gains in our understanding of how to control the errors of discretization represent only half of the control necessary to ensure that a numerical model is in fact an accurate representation of the corresponding physical system. Control of the errors of idealization is equally important. A brief overview of the main ideas of how to ensure the quality and reliability of mathematical models of structural systems is presented.

Szabo, Barna A.↗

Error analysis of an aircraft-mounted backscatter ultraviolet experiment

An experiment has been proposed to measure the zenith sky intensity and directly transmitted solar flux density at several wavelengths in the near ultraviolet from an aircraft platform located at about 15 km. From these measurements, the vertical distribution of ozone above the platform is to be inferred, by inversion of the equation of radiative transfer. The present study is an analysis of the errors to be expected in a real experiment, due to atmospheric fluctuations, numerical procedures, and ignorance of various parameters necessary for the calculation of the intensities and fluxes to be expected. These errors are found to vary from 17% at the shorter wavelengths to about 3% at the longer wavelengths used. It is shown elsewhere (Twomey et al., 1977) that unacceptable levels of instability and non-uniqueness of results will be found if the error levels in the measurements upon which we wish to invert exceed 3%. On the basis of this, it is concluded that with the present instrumentation and the present knowledge of necessary ancillary parameters, this experiment will not be feasible.

Rabinoff, R.↗

The error associated with use of various forms of the divergence equation to diagnose geopotential and temperature

An experimental study was made of errors arising from diagnosing the geopotential and temperature for numerical weather prediction models by means of an inviscid divergence equation (DE). The experiment consisted of a simulation of a snow storm which occurred in 1979, with observational data being input for every point of the model grid and for the geopotential boundary condition, thereby eliminating errors in the horizontal and vertical analysis. The DE terms were recorded at each grid point, and were found to include synoptic-scale geopotential height errors up to 23 m and rms temperature errors exceeding 26 C in the PBL. The errors were minimal above 800 mb. The various unsuccessful attempts to reduce or eliminate the errors are described.

Modica, G. M.↗

On the measurement of lateral velocity derivatives in turbulent flows

Direct numerical simulation data for the lateral velocity derivative delta(u)/delta(y) at the centerline of a fully developed turbulent channel flow provide reasonable support for Wyngaard's analysis of the error involved in measuring this quantity using parallel hot wires. Numerical data in the wall region of the channel flow also provide a useful indication of how to select the separation between the wires. Justification for this choice is obtained by comparing several measured statistics of delta(u)/delta(y) with the corresponding numerical data.

Antonia, R. A.↗

The role of global cloud climatologies in validating numerical models

The purpose of this work is to estimate sampling errors of area-time averaged rain rate due to temporal samplings by satellites. In particular, the sampling errors of the proposed low inclination orbit satellite of the Tropical Rainfall Measuring Mission (TRMM) (35 deg inclination and 350 km altitude), one of the sun synchronous polar orbiting satellites of NOAA series (98.89 deg inclination and 833 km altitude), and two simultaneous sun synchronous polar orbiting satellites--assumed to carry a perfect passive microwave sensor for direct rainfall measurements--will be estimated. This estimate is done by performing a study of the satellite orbits and the autocovariance function of the area-averaged rain rate time series. A model based on an exponential fit of the autocovariance function is used for actual calculations. Varying visiting intervals and total coverage of averaging area on each visit by the satellites are taken into account in the model. The data are generated by a General Circulation Model (GCM). The model has a diurnal cycle and parameterized convective processes. A special run of the GCM was made at NASA/GSFC in which the rainfall and precipitable water fields were retained globally for every hour of the run for the whole year.

HARSHVARDHAN↗

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗