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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Double quadriphase modulation/demodulation technique for three-channel communication link

A modulation technique for a three-channel communication link is introduced. The structure of the modulator is such as to form an unbalanced quadriphase signal wherein the high rate data stream is bi-phase modulated on the in-phase carrier component, while the sum of the two lower rate signals is bi-phase modulated on the quadrature component of the same carrier. The sum of the two lower signals is, in turn, formed by modulating with the respective data streams the in-phase and the quadrature components of a square wave subcarrier. At the demodulator, the tracking of the carrier and the subcarrier is performed by two independent Costas loops. The demodulation of the high data rate signal is carried out after establishing the carrier reference signal, while the lower rate signals are demodulated after the subcarrier loop recovers the subcarrier. In this paper, the performance of the two loops is analyzed and the expressions for the tracking errors are derived. Finally, a numerical example pertaining to the Space Shuttle-to-TDRS Ku-band link is presented for illustration.

Alem, W. K.↗

NASA NTRS1985

Medium range forecasting: The experience of ECMWF

The results of five years of operational medium-range forecasts at ECMWF is discussed. It is shown that a considerable improvement in predictive skill has taken place, resulting in more accurate forecasts and an extension in time of useful predictive skill. For the Northern Hemisphere extratropics, the improvement in skill during the winter season is about three days for the 500-mb geopotential compared to the very first experiments in the early 1970s. Three are large geographical and temporal variations in the performance of the model, as well as large-scale errors of a systematic nature. Numerical experiments have indicated great sensitivity to the representation of orography and to tropical forcing.

Bengtsson, L.↗

NASA NTRS1986

A spline-based parameter estimation technique for static models of elastic structures

The problem of identifying the spatially varying coefficient of elasticity using an observed solution to the forward problem is considered. Under appropriate conditions this problem can be treated as a first order hyperbolic equation in the unknown coefficient. Some continuous dependence results are developed for this problem and a spline-based technique is proposed for approximating the unknown coefficient, based on these results. The convergence of the numerical scheme is established and error estimates obtained.

Dutt, P.↗

NASA NTRS1989

Structural damage assessment as an identification problem

Damage assessment of structural assemblies is treated as an identification problem. A brief review of identification methods is first presented with particular focus on the output error approach. The use of numerical optimization methods in identifying the location and extent of damage in structures is studied. The influence of damage on eigenmode shapes and static displacements is explored as a means of formulating a measure of damage in the structure. Preliminary results obtained in this study are presented and special attention is directed at the shortcomings associated with the nonlinear programming approach to solving the optimization problem.

Hajela, Prabhat↗

NASA NTRS1988

The Galerkin/least-squares method for advective-diffusive equations

Galerkin/least-squares finite-element methods are presented for advective-diffusive equations. Galerkin/least-squares represents a conceptual simplification of streamline-upwind Petrov-Galerkin methods, and is in fact applicable to a wide variety of other problem types. A convergence analysis and error estimates are presented. Some numerical results for compressible Navier-Stokes flows are presented.

Hughes, T. J. R.↗

NASA NTRS1989

A spline-based parameter estimation technique for static models of elastic structures

Dutt, P.↗

NASA NTRS1992

Analysis of stellar occultation data for planetary atmospheres. I - Model fitting, with application to Pluto

Consideration is given to an analytic model for a stellar-occultation light curve developed for a small, spherically symmetric planetary atmosphere that includes thermal and molecular weight gradients in a region that overlies an extinction layer. The model incorporates two equivalent sets of parameters. One set specifies the occultation light curve in terms of signal levels, times, and time intervals. The other set specifies physical parameters of the planetary atmosphere. Equations are given for the transforming between the sets of parameters, including their errors and correlation coefficients. Detailed numerical calculations are presented for a benchmark case. The results obtained are consistent with the isothermal prediction of the 'methane-thermostat' model of Pluto's atmosphere.

Elliot, J. L.↗

NASA NTRS1990

Simplified detection and correction of critical data for ill-conditioned systems

Ill-conditioned systems arising in analysis and optimization can display a high sensitivity to numerical precision for changes and errors in data input. Such data may be in the form of system parameter input or desired system response. The ill-conditioning referred to generally arises from the lack of sufficient independent data to define a complex system or the weak sensitivity of response to source input parameters. It is shown how small errors in data and assumed fixed and known parameters can lead to highly erroneous results in ill-conditioned linear algebraic equations. A simplified detection and correction of critical input data arising in the coefficient matrix and desired response (i.e., right hand side) is proposed.

Ojalvo, Irving U.↗

NASA NTRS1995

A compact solution to computational acoustics

This paper demonstrates that the linearized, dimensional Euler equations for acoustic computation can be accurately solved as a set of decoupled first-order wave equations, and that if ordered properly, this system of simple waves has unambiguous, easily implemented boundary conditions, allowing waves of same group speeds to pass through numerical boundaries or comply with wall conditions. Thus, the task of designing a complex multi-dimensional scheme with approximate far-field boundary conditions reduces to the design of higher order schemes for the one-dimensional simple wave equation. A compact finite-difference scheme and a characteristically exact but numerically n(th) order accurate boundary condition are introduced for solving the first order wave equation. Spanning a three-point two-level stencil, this low-dispersion implicit scheme has a third order spatial accuracy when used on nonuniform meshes, fourth order accurate on uniform meshes, and a temporal accuracy of second order due to the choice of trapezoidal integration for algorithmic simplicity. The robustness and accuracy of the scheme are demonstrated through a series of numerical experiments and comparisons with published results. When tested on the one-dimensional wave equation on a uniform grid, this scheme allows a Gaussian wave packet to pass through any finite domain with low numerical dispersion characteristic of a spatially fourth-order scheme and reflections at numerical boundaries maintained below truncation error. On highly stretched and irregular grids, only mild dispersions are found in the solution while solutions by other methods fail or are severely distorted. Yet, this scheme is no more sophisticated to solve or implement than the Crank-Nicolson scheme. This scheme has been tested on four categories of the ICASE/LaRC benchmark problems, which include propagation of acoustic and convective waves in Cartesian and cylindrical domains, reflection of acoustic wave at stationary/moving boundaries, and sound generation by gust-blade interaction.

Fung, K.-Y.↗

NASA NTRS1993

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗

NASA NTRS2000

General Purpose Fortran Program for Discrete-Ordinate-Method Radiative Transfer in Scattering and Emitting Layered Media: An Update of DISORT

This update reports a state-of-the-art discrete ordinate algorithm for monochromatic unpolarized radiative transfer in non-isothermal, vertically inhomogeneous, but horizontally homogeneous media. The physical processes included are Planckian thermal emission, scattering with arbitrary phase function, absorption, and surface bidirectional reflection. The system may be driven by parallel or isotropic diffuse radiation incident at the top boundary, as well as by internal thermal sources and thermal emission from the boundaries. Radiances, fluxes, and mean intensities are returned at user-specified angles and levels. DISORT has enjoyed considerable popularity in the atmospheric science and other communities since its introduction in 1988. Several new DISORT features are described in this update: intensity correction algorithms designed to compensate for the 8-M forward-peak scaling and obtain accurate intensities even in low orders of approximation; a more general surface bidirectional reflection option; and an exponential-linear approximation of the Planck function allowing more accurate solutions in the presence of large temperature gradients. DISORT has been designed to be an exemplar of good scientific software as well as a program of intrinsic utility. An extraordinary effort has been made to make it numerically well-conditioned, error-resistant, and user-friendly, and to take advantage of robust existing software tools. A thorough test suite is provided to verify the program both against published results, and for consistency where there are no published results. This careful attention to software design has been just as important in DISORT's popularity as its powerful algorithmic content.

Tsay, Si-Chee↗

NASA NTRS1994

Merging Intelligent Systems Technologies with CFD Analysis Strategies: Prototype Development at NASA Ames

A collaborative team of researchers from fields of Computational Fluid Dynamics (CFD), fluid physics, computer architectures, and computer science and knowledge engineering have begun work on a prototype system that addresses several of industry's concerns in using NASA-developed CFD codes as part of the design cycle. A major problem exists in the application of CFD technologies within the aeronautics design cycle due primarily to misunderstandings in the ranges of applicability of the various solver codes or turbulence models. Features that arise during the CFD solution process need to be discriminated and recognized as actual flow features with physical support in the geometry and flow conditions of the problem being solved, or as numerical or non-physical errors arising from mis-application of solver code and its parameters, gridding strategies, or discretization. interpolations. The fundamental concept is to develop an intelligent computational system that can accept the engineer's definition of the problem and construct an optimal CFD solution. To do this requires capturing both the knowledge of how to apply the various CFD tools and how to adapt the application of those tools to flow structures as they evolve during the flow simulation. Embedded within this adaptive system approach is the additional desire to automatically identify and quantify the quality of resolution of the pertinent flow structures, be they genuine or error-induced, and then to adjust the solution strategy accordingly. This paper discusses the status of that prototyping effort.

Thompson, David E.↗

NASA NTRS2003

Satellite Detection of Precipitation over the North Pacific

Intense extratropical winter cyclones often impact the West Coast of North America with strong winds and heavy precipitation. Several times during a winter season, short-term forecasts (24 - 48 hours) of these storms are seriously deficient with central pressure errors in the 10's of hPa and surface low position errors in the 100's of km. For example, 48-hr sea level pressure errors (forecast - observation) at buoy 46005 off the Oregon coast for the 2001 - 2002 winter season is plotted. In addition, two times the standard deviation (determined from pressure errors from the last four winter seasons) are also shown. It is evident from this figure that large forecast errors (i.e. greater than 10 hPa) occurred about 10 times this past winter at buoy 46005 with three events where the errors were 20 hPa. Beside large forecast errors of sea level pressure, numerical forecasts of precipitation for land falling cyclones can also be flawed. This is due in large part to the lack of accurate precipitation information over the ocean. Therefore, remote sensing techniques are the only viable option for obtaining accurate information on the distribution and intensity of precipitation over the North Pacific. Due to the radiative characteristics of precipitation sized hydrometeors at microwave frequencies, microwave sensors are able to detect precipitation over oceanic regions. Past studies have demonstrated the utility of passive microwave rainrate data for locating intense rainfall in rapidly deepening cyclones, in detecting developing polar mesocyclones and in determining frontal bands. There are currently many sources of microwave rainrate data: the Special Sensor Microwave Imager (SSM/I) (currently flying on three platforms), the Advanced Microwave Sounding Unit (AMSU-B) (currently flying on NOAA-15, NOAA-16, and NOAA-17), and the Tropical Rainfall Measuring Mission Microwave Imager (TMI). Data will soon be available from the Advanced Microwave Radiometer-EOS (AMSR-E) on the Aqua platform. In this paper, we present a new technique for mapping rainrate distributions over the North Pacific utilizing rainrate estimates from several microwave sensors and upper-tropospheric winds derived from geosynchronous satellite IR data. The goal of this work is to develop a way to obtain high temporal and spatial rainfall information over the North Pacific. This information will be used to support the verification of model derived precipitation distributions and to support the analysis of in situ measurements of rainfall during the Improvement of Microphysical Parameterization through Observational Verification Experiment (IMPROVE) field campaigns.

Smith, Jeremy↗

NASA NTRS2010

Reduced Order Modeling Incompressible Flows

The details: a) Need stable numerical methods; b) Round off error can be considerable; c) Not convinced modes are correct for incompressible flow. Nonetheless, can derive compact and accurate reduced-order models. Can be used to generate actuator models or full flow-field models

Helenbrook, B. T.↗

NASA NTRS2010

Multiple Velocity Profile Measurements in Hypersonic Flows Using Sequentially-Imaged Fluorescence Tagging

Nitric-oxide planar laser-induced fluorescence (NO PLIF) was used to perform velocity measurements in hypersonic flows by generating multiple tagged lines which fluoresce as they convect downstream. For each laser pulse, a single interline, progressive scan intensified CCD (charge-coupled device) camera was used to obtain two sequential images of the NO molecules that had been tagged by the laser. The CCD configuration allowed for sub-microsecond acquisition of both images, resulting in sub-microsecond temporal resolution as well as sub-mm spatial resolution (0.5-mm horizontal, 0.7-mm vertical). Determination of axial velocity was made by application of a cross-correlation analysis of the horizontal shift of individual tagged lines. A numerical study of measured velocity error due to a uniform and linearly-varying collisional rate distribution was performed. Quantification of systematic errors, the contribution of gating/exposure duration errors, and the influence of collision rate on temporal uncertainty were made. Quantification of the spatial uncertainty depended upon the signal-to-noise ratio of the acquired profiles. This velocity measurement technique has been demonstrated for two hypersonic flow experiments: (1) a reaction control system (RCS) jet on an Orion Crew Exploration Vehicle (CEV) wind tunnel model and (2) a 10-degree half-angle wedge containing a 2-mm tall, 4-mm wide cylindrical boundary layer trip. The experiments were performed at the NASA Langley Research Center's 31-Inch Mach 10 Air Tunnel.

Bathel, Brett F.↗

NASA NTRS2015

A Comparison of Three Algorithms for Orion Drogue Parachute Release

The Orion Multi-Purpose Crew Vehicle is susceptible to ipping apex forward between drogue parachute release and main parachute in ation. A smart drogue release algorithm is required to select a drogue release condition that will not result in an apex forward main parachute deployment. The baseline algorithm is simple and elegant, but does not perform as well as desired in drogue failure cases. A simple modi cation to the baseline algorithm can improve performance, but can also sometimes fail to identify a good release condition. A new algorithm employing simpli ed rotational dynamics and a numeric predictor to minimize a rotational energy metric is proposed. A Monte Carlo analysis of a drogue failure scenario is used to compare the performance of the algorithms. The numeric predictor prevents more of the cases from ipping apex forward, and also results in an improvement in the capsule attitude at main bag extraction. The sensitivity of the numeric predictor to aerodynamic dispersions, errors in the navigated state, and execution rate is investigated, showing little degradation in performance.

Matz, Daniel A.↗

NASA NTRSDate not supplied

Theodorsen's and Garrick's Computational Aeroelasticity, Revisited

This paper describes a multiyear effort to recompute all of the numerical flutter results contained in Theodore Theodorsen’s and I.E. Garrick’s (“T&G’s”) groundbreaking trilogy of National Advisory Committee for Aeronautics technical reports on aeroelastic flutter (NACA Reports 496, 685, and 741). The paper includes overviews of T&G’s aeroelastic equations and solution methods, comparisons between the original and recomputed numerical results, and instances of errors and tripping points (potentials for confusing the reader) in the original NACA reports.

Flutter↗

NASA NTRS2021

Evaluation of GEOS Precipitation Flagging for SMAP Soil Moisture Retrieval Accuracy

The precipitation flag in the Soil Moisture Active Passive (SMAP) Level 2 passive soil moisture (L2SMP) retrieval product indicates the presence or absence of heavy precipitation at the time of the SMAP overpass. The flag is based on precipitation estimates from the Goddard Earth Observing System (GEOS) Forward Processing numerical weather prediction system. An error in flagging during an active or recent precipitation event can produce either 1) an overestimation of soil moisture due to short-term surface wetting of vegetation and/or surface ponding (if soil moisture retrieval was attempted in the presence of rain) or 2) an unnecessary nonretrieval of soil moisture and loss of data (if retrieval is flagged due to an erroneous indication of rain). Satellite precipitation estimates from the Integrated Multisatellite Retrievals for GPM (IMERG), version 06, Early Run (latency of ~4 h) precipitationCal product are used here to evaluate the GEOS-based precipitation flag in the L2SMP product for both the 1800 local time (LT) ascending and 0600 LT descending SMAP overpasses over the first five years of the mission (2015–20). Consisting of blended precipitation measurements from the Global Precipitation Mission (GPM) satellite constellation, IMERG is treated as the “truth” when comparing to the GEOS model forecasts of precipitation used by SMAP. Key results include (i) IMERG measurements generally show higher spatial variability than the GEOS forecast precipitation, (ii) the IMERG product has a higher frequency of light precipitation amounts, and (iii) the effect of incorporating IMERG rainfall measurements in lieu of GEOS precipitation forecasts are minimal on the L2SMP retrieval accuracy (determined vs in situ soil moisture measurements at core validation sites). Our results indicate that L2SMP retrievals continue to meet the mission’s accuracy requirement [standard deviation of the unbiased RMSE (ubRMSE) less than 0.04 cu. m/cu. m].

IMERG-precipitationCal↗