Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical 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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Error indicators and accuracy improvements of finite element solutions

Practical and reliable estimators of the discretization errors in engineering problems are developed. Error indicators for identifying the regions or elements of the solution domain which are likely to have the largest discretization errors are presented, and a simple computational procedure for improving the accuracy of the finite element solutions for shell problems is given. The similarities between the proposed procedure and a preconditioned conjugate gradient (PCG) technique are identified and exploited to generate pointwise error indicators from the PCG technique. Numerical examples in the linear static analysis of shells are presented.

Noor, Ahmed K.↗
NASA NTRS1989

Transition regime droplet growth and evaporation - An integrodifferential variational approach

A variational principle based on the integrodifferential form of the Boltzmann equation, which allows very general forms of the collision integral and of the boundary operator, has been used to compute the mass flux to a spherical particle. For simplicity the problem has been restricted to the one-speed, constant cross-section approximation, the black sphere problem of neutron transport theory. Using the Hilbert expansion of the solution far from the sphere as a trial function leads to a rational expression in the inverse Knudsen number for the extrapolation distance, which is simply related to the number flux to the sphere. A five-term expansion gives results that are no more than 6 percent in error when compared to an accurate numerical solution. The technique can be applied to other forms of the collision model and other boundary conditions.

Cipolla, J. W., Jr.↗
NASA NTRS1990

Optimal model reduction and frequency-weighted extension

In this paper the quadratically optimal model reduction problem for single-input, single-output systems is considered. The reduced order model is determined by minimizing the integral of the magnitude-squared of the transfer function error. It is shown that the numerator coefficients of the optimal approximant satisfy a weighted least squares problem and, on this basis, a two-step iterative algorithm is developed combining a least squares solver with a gradient minimizer. The existence of globally optimal stable solutions to the optimization problem is established, and convergence of the algorithm to stationary values of the cost function is proved. The formulation is extended to handle the frequency-weighted optimal model reduction problem. Three examples demonstrate the optimization algorithm.

Spanos, J. T.↗
NASA NTRS1992

Lumped mass modelling for the dynamic analysis of aircraft structures

Aircraft structures may be modelled by lumping the masses at particular strategic points and the flexibility or stiffness of the structure is obtained with reference to these points. Equivalent moments of inertia for the section at these positions are determined. The lumped masses are calculated based on the assumption that each point will represent the mass spread on one half of the space on each side. Then these parameters are used in the differential equation of motion and the eigen characteristics are determined. A comparison is made with results obtained by other established methods. The lumped mass approach in the dynamic analysis of complicated structures provides an easier means of predicting the dynamic characteristics of these structures. It involves less computer time and avoids computational errors that are inherent in the numerical solution of complicated systems.

Abu-Saba, Elias G.↗
NASA NTRS1992

A new algorithm for L2 optimal model reduction

In this paper the quadratically optimal model reduction problem for single-input, single-output systems is considered. The reduced order model is determined by minimizing the integral of the magnitude-squared of the transfer function error. It is shown that the numerator coefficients of the optimal approximant satisfy a weighted least squares problem and, on this basis, a two-step iterative algorithm is developed combining a least squares solver with a gradient minimizer. Convergence of the proposed algorithm to stationary values of the quadratic cost function is proved. The formulation is extended to handle the frequency-weighted optimal model reduction problem. Three examples demonstrate the optimization algorithm.

Spanos, J. T.↗
NASA NTRS1993

Attitude control with realization of linear error dynamics

An attitude control law is derived to realize linear unforced error dynamics with the attitude error defined in terms of rotation group algebra (rather than vector algebra). Euler parameters are used in the rotational dynamics model because they are globally nonsingular, but only the minimal three Euler parameters are used in the error dynamics model because they have no nonlinear mathematical constraints to prevent the realization of linear error dynamics. The control law is singular only when the attitude error angle is exactly pi rad about any eigenaxis, and a simple intuitive modification at the singularity allows the control law to be used globally. The forced error dynamics are nonlinear but stable. Numerical simulation tests show that the control law performs robustly for both initial attitude acquisition and attitude control.

Paielli, Russell A.↗
NASA NTRS1947

Determination of elastic stresses in gas-turbine disks

A method is presented for the calculation of elastic stresses in symmetrical disks typical of those of a high-temperature gas turbine. The method is essentially a finite-difference solution of the equilibrium and compatibility equations for elastic stresses in a symmetrical disk. Account can be taken of point-to-point variations in disk thickness, in temperature, in elastic modulus, in coefficient of thermal expansion, in material density, and in Poisson's ratio. No numerical integration or trial-and-error procedures are involved and the computations can be performed in rapid and routine fashion by nontechnical computers with little engineering supervision. Checks on problems for which exact mathematical solutions are known indicate that the method yields results of high accuracy. Illustrative examples are presented to show the manner of treating solid disks, disks with central holes, and disks constructed either of a single material or two or more welded materials. The effect of shrink fitting is taken into account by a very simple device.

Manson, S S↗
NASA NTRS1996

Effects of modeling errors on trajectory predictions in air traffic control automation

Air traffic control automation synthesizes aircraft trajectories for the generation of advisories. Trajectory computation employs models of aircraft performances and weather conditions. In contrast, actual trajectories are flown in real aircraft under actual conditions. Since synthetic trajectories are used in landing scheduling and conflict probing, it is very important to understand the differences between computed trajectories and actual trajectories. This paper examines the effects of aircraft modeling errors on the accuracy of trajectory predictions in air traffic control automation. Three-dimensional point-mass aircraft equations of motion are assumed to be able to generate actual aircraft flight paths. Modeling errors are described as uncertain parameters or uncertain input functions. Pilot or autopilot feedback actions are expressed as equality constraints to satisfy control objectives. A typical trajectory is defined by a series of flight segments with different control objectives for each flight segment and conditions that define segment transitions. A constrained linearization approach is used to analyze trajectory differences caused by various modeling errors by developing a linear time varying system that describes the trajectory errors, with expressions to transfer the trajectory errors across moving segment transitions. A numerical example is presented for a complete commercial aircraft descent trajectory consisting of several flight segments.

Jackson, Michael R. C.↗
NASA NTRS2000

Computerized Generation and Simulation of Meshing and Contact of New Type of Novikov-Wildhaber Helical Gears

In this report, we propose a new geometry for low-noise, increased-strength helical gears of the Novikov-Wildhaber type. Contact stresses are reduced as a result of their convex-concave gear tooth surfaces. The gear tooth surfaces are crowned in the profile direction to localize bearing contact and in the longitudinal direction to obtain a parabolic function of transmission errors. Such a function results in the reduction of noise and vibrations. Methods for the generation of the proposed gear tooth surfaces by grinding and hobbing are considered, and a tooth contact analysis (TCA) computer program to simulate meshing and contact is applied. The report also investigates the influence of misalignment on transmission errors and shift of bearing contact. Numerical examples to illustrate the developed approaches are proposed. The proposed geometry was patented by Ford/UIC (Serial Number 09-340-824, pending) on June 28, 1999.

Litvin, Faydor L.↗
NASA NTRS2001

Improved Use of Satellite Imagery to Forecast Hurricanes

This project tested a novel method that uses satellite imagery to correct phase errors in the initial state for numerical weather prediction, applied to hurricane forecasts. The system was tested on hurricanes Guillermo (1997), Felicia (1997) and Iniki (1992). We compared the performance of the system with and without phase correction to a procedure that uses bogus data in the initial state, similar to current operational procedures. The phase correction keeps the hurricane on track in the analysis and is far superior to a system without phase correction. Compared to operational procedure, phase correction generates somewhat worse 3-day forecast of the hurricane track, but better forecast of intensity. It is believed that the phase correction module would work best in the context of 4-dimensional variational data assimilation. Very little modification to 4DVar would be required.

Louis, Jean-Francois↗
NASA NTRS1994

Impact of Uncertainties and Errors in Converting NWS Radiosonde Hygristor Resistances to Relative Humidity

A set of 2,600 6-second, National Weather Service soundings from NASA's FIRE-II Cirrus field experiment are used to illustrate previously known errors and new potential errors in the VIZ and SDD ) brand relative humidity (RH) sensors and the MicroART processing software. The entire spectrum of RH is potentially affected by at least one of these errors. (These errors occur before being converted to dew point temperature.) Corrections to the errors are discussed. Examples are given of the effect that these errors and biases may have on numerical weather prediction and radiative transfer. The figure shows the OLR calculated for the corrected and uncorrected soundings using an 18-band radiative transfer code. The OLR differences are sufficiently large to warrant consideration when validating line-by-line radiation calculations that use radiosonde data to specify the atmospheric state, or when validating satellite retrievals. in addition, a comparison of observations of RH during FIRE-II derived from GOES satellite, raman lidar, MAPS analyses, NCAR CLASS sondes, and the NWS sondes reveals disagreement in the RH distribution and underlines our lack of an understanding of the climatology of water vapor.

Westphal, Douglas L.↗
NASA NTRS1994

Impact of Uncertainties and Errors in Converting NWS Radiosonde Hygristor Resistances to Relative Humidity

A set of 2,600 6-second, National Weather Service soundings from NASA's FIRE-II Cirrus field experiment are used to illustrate previously known errors and new potential errors in the VIZ and SDD brand relative humidity (RH) sensors and the MicroART processing software. The entire spectrum of RH is potentially affected by at least one of these errors. (These errors occur before being converted to dew point temperature.) Corrections to the errors are discussed. Examples are given of the effect that these errors and biases may have on numerical weather prediction and radiative transfer. The figure shows the OLR calculated for the corrected and uncorrected soundings using an 18-band radiative transfer code. The OLR differences are sufficiently large to warrant consideration when validating line-by-line radiation calculations that use radiosonde data to specify the atmospheric state, or when validating satellite retrievals. In addition, a comparison of observations of RE during FIRE-II derived from GOES satellite, raman lidar, MAPS analyses, NCAR CLASS sondes, and the NWS sondes reveals disagreement in the RH distribution and underlines our lack of an understanding of the climatology of water vapor.

Westphal, Douglas L.↗
NASA NTRS2012

Dual-Pulse Pulse Position Modulation (DPPM) for Deep-Space Optical Communications: Performance and Practicality Analysis

Due to its simplicity and robustness against wavefront distortion, pulse position modulation (PPM) with photon counting detector has been seriously considered for long-haul optical wireless systems. This paper evaluates the dual-pulse case and compares it with the conventional single-pulse case. Analytical expressions for symbol error rate and bit error rate are first derived and numerically evaluated, for the strong, negative-exponential turbulent atmosphere; and bandwidth efficiency and throughput are subsequently assessed. It is shown that, under a set of practical constraints including pulse width and pulse repetition frequency (PRF), dual-pulse PPM enables a better channel utilization and hence a higher throughput than it single-pulse counterpart. This result is new and different from the previous idealistic studies that showed multi-pulse PPM provided no essential information-theoretic gains than single-pulse PPM.

Li, Jing↗
NASA NTRSDate not supplied

Simulations of the Modular Axisymmetric Scramjet Test Rig Under Reacting Flow Conditions

Simulations of the Air Force Research Laboratory (AFRL) Modular Axisymmetric Scramjet Test Rig (MASTeR) are presented. MASTeR is a parametric test article capable of investigating various scramjet cavity flameholder designs and fueling strategies with the goal to characterize and optimize flameholding capability. In the current work, three cavity aspect ratios, and two fueling strategies (upstream and in-cavity) with ethylene at a nominal facility pressure and temperature condition are evaluated. The reacting simulations are performed, and the resulting flow characteristics are discussed. For each configuration and fueling strategy, cavity residence time, entrainment rate, and fuel-air equivalence ratio are computed. The MASTeR geometry is defined in Engineering Sketch Pad (ESP) and the simulations use a sketch-to-solution (S2S) automated unstructured grid adaptation tool in VULCAN-CFD. This tool automatically generates a simulation grid from the ESP geometry and systematically adapts it to the numerical solution based on the Hessian error estimate of a specified flow field parameter. Reynolds averaged simulations (RAS) are used with a two-equation linear eddy viscosity and diffusivity model. The resulting database can be compared with the experimental data as those become available, and explored to develop models for cavity performance for scramjet propulsion design applications.

hypersonics↗
NASA NTRSDate not supplied

Simulations of the Modular Axisymmetric Scramjet Test Rig Under Reacting Flow Conditions

Simulations of the Air Force Research Laboratory (AFRL) Modular Axisymmetric Scramjet Test Rig (MASTeR) are presented. MASTeR is a parametric test article capable of investigating various scramjet cavity flameholder designs and fueling strategies with the goal to characterize and optimize flameholding capability. In the current work, three cavity aspect ratios, three depths, and two fueling strategies (upstream and in-cavity) with ethylene at a nominal facility pressure and temperature conditions are evaluated. The simulations are performed for mixing-only and reacting flows, and the resulting flow characteristics are compared. For each configuration and fueling strategy, cavity residence time, entrainment rate, and fuel-air equivalence ratio are computed. The MASTeR geometry is defined in the Engineering Sketch Pad (ESP) and the simulations use a sketch-to-solution (S2S) automated unstructured grid adaptation tool in VULCAN-CFD. This tool automatically generates a simulation grid from the ESP geometry and systematically adapts it to the numerical solution based on the Hessian error estimate of a specified flow field parameter. Reynolds averaged simulations (RAS) are used with typical two-equation linear eddy viscosity and diffusivity model. The resulting database can be compared with the experimental data as those becomes available and explored to develop models for cavity performance for scramjet propulsion design applications.

hypersonics↗
DOE OSTI2013

Recent developments in DYNSUB: New models, code optimization and parallelization

DYNSUB is a high-fidelity coupled code system consisting of the reactor simulator DYN3D and the sub-channel code SUBCHANFLOW. It describes nuclear reactor core behavior with pin-by-pin resolution for both steady-state and transient scenarios. In the course of the coupled code system's active development, super-homogenization (SPH) and generalized equivalence theory (GET) discontinuity factors may be computed with and employed in DYNSUB to compensate pin level homogenization errors. Because of the largely increased numerical problem size for pin-by-pin simulations, DYNSUB has bene fitted from HPC techniques to improve its numerical performance. DYNSUB's coupling scheme has been structurally revised. Computational bottlenecks have been identified and parallelized for shared memory systems using OpenMP. Comparing the elapsed time for simulating a PWR core with one-eighth symmetry under hot zero power conditions applying the original and the optimized DYNSUB using 8 cores, overall speed up factors greater than 10 have been observed. The corresponding reduction in execution time enables a routine application of DYNSUB to study pin level safety parameters for engineering sized cases in a scientific environment. (authors)

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗
NASA NTRS1978

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗