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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Airfoil computation at high angles of attack, inviscid and viscous phenomena

An implicit central difference code is used to calculate two dimensional inviscid and thin-layer Navier-Stokes solutions for flow about an NACA0012 airfoil at high angles of attack. Among the issues addressed are whether separation can occur in an inviscid calculation and what the causes would be of such separation. Examples are shown of inviscid shocked flow with and without separation and shock-free flow with separation. An Euler solution with self-induced oscillation and separation driven by a strong shock is contrasted with a shock-free solution whose separation is caused by numerical error. Computed solutions to the Euler equations are compared to those of the potential equations. Comparisons are also made between experimental data from wind tunnel tests and viscous calculations at similar conditions.

Barton, J. T.↗
NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗
NASA NTRS1986

Computation of leading-edge vortex flows

The simulation of the leading edge vortex flow about a series of conical delta wings through solution of the Navier-Stokes and Euler equations is studied. The occurrence, the validity, and the usefulness of separated flow solutions to the Euler equations of particular interest. Central and upwind difference solutions to the governing equations are compared for a series of cross sectional shapes, including both rounded and sharp tip geometries. For the rounded leading edge and the flight condition considered, viscous solutions obtained with either central or upwind difference methods predict the classic structure of vortical flow over a highly swept delta wing. Predicted features include the primary vortex due to leading edge separation and the secondary vortex due to crossflow separation. Central difference solutions to the Euler equations show a marked sensitivity to grid refinement. On a coarse grid, the flow separates due to numerical error and a primary vortex which resembles that of the viscous solution is predicted. In contrast, the upwind difference solutions to the Euler equations predict attached flow even for first-order solutions on coarse grids. On a sufficiently fine grid, both methods agree closely and correctly predict a shock-curvature-induced inviscid separation near the leeward plane of symmetry. Upwind difference solutions to the Navier-Stokes and Euler equations are presented for two sharp leading edge geometries. The viscous solutions are quite similar to the rounded leading edge results with vortices of similar shape and size. The upwind Euler solutions predict attached flow with no separation for both geometries. However, with sufficient grid refinement near the tip or through the use of more accurate spatial differencing, leading edge separation results. Once the leading edge separation is established, the upwind solution agrees with recently published central difference solutions to the Euler equations.

Newsome, R. W.↗
NASA NTRS1986

A numerical model of acoustic choking. II - Shocked solutions

The one dimensional equations of gas dynamics are used to model subsonic acoustic choking. This model can accommodate non-linear distortion of waves and the eventual formation of shock waves. Several finite differencing schemes are adapted to obtain solutions. The results obtained with the various schemes are compared with the asymptotic results available. The results suggest that no one finite differencing scheme gives solutions significantly better than the others and that most of the difference solutions are close to the asymptotic results. If the acoustic shock wave is sufficiently strong it almost annihilates the acoustic wave; in this situation numerical errors may dominate the results. Such solutions involve very large acoustic attenuations.

Walkington, N. J.↗
NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗
NASA NTRS1986

Modified LTR robust control for flexible structures

A procedure is developed for dealing with performance and robustness issues in the design of multi-input multi-output compensators for lightly damped flexible structures. The procedure is based upon representing errors in the plant design model as structured uncertainties, and applying a modified version of the Loop Transfer Recovery (LTR) design method. Real parameter errors, such as frequency errors, damping errors or modal displacement errors can be treated. The approach involves adjusting the cost function in the regulator problem and the process noise model in the estimator problem in a particular manner which reflects the assumed structure of the modeling errors. Numerical examples dealing with the control of a large flexible space antenna with uncertain frequencies demonstrate a considerable improvement over standard LTR methods. Convenient design parameters can be varied until a satisfactory compromise is achieved between performance and robustness.

Blelloch, P. A.↗
NASA NTRS1988

Radiometric calibration analysis of SIR-B imagery

The second flight of the NASA Spaceborne Imaging Radar (SIR-B) collected nearly 8 hours of digital synthetic-aperture radar (SAR) data of the earth's surface. This paper analyzes the performance of the SAR to determine the accuracy of the relative radiometric calibration of these data. Procedures are presented for deriving accuracy estimates as well as limitations under which these estimates are valid. The impact on calibration of an RF breakdown in the primary antenna feed system cable is evaluated. An analysis of the platform stability based on the SAR echo data is also presented in conjunction with its resultant effect on the calibration accuracy. Finally, numerical error bounds are derived with guidelines on their utilization.

Wall, Stephen D.↗
NASA NTRS1988

Numerical simulation of the integrated space shuttle vehicle in ascent

A simulation of the flow about the integrated space shuttle vehicle in ascent mode has been undertaken for various flight conditions using the Chimera composite grid discretization approach. Overset body-conforming grids were used to represent each geometric component, and an implicit approximately factored finite-difference procedure was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computational results have been compared with both wind tunnel and flight test data. Although relatively good agreement is obtained with the experimental data, further refinement and evaluation of numerical error is under way.

Buning, P. G.↗
NASA NTRS1989

Toward a CFD nose-to-tail capability - Hypersonic unsteady Navier-Stokes code validation

Computational fluid dynamics (CFD) research for hypersonic flows presents new problems in code validation because of the added complexity of the physical models. This paper surveys code validation procedures applicable to hypersonic flow models that include real gas effects. The current status of hypersonic CFD flow analysis is assessed with the Compressible Navier-Stokes (CNS) code as a case study. The methods of code validation discussed to beyond comparison with experimental data to include comparisons with other codes and formulations, component analyses, and estimation of numerical errors. Current results indicate that predicting hypersonic flows of perfect gases and equilibrium air are well in hand. Pressure, shock location, and integrated quantities are relatively easy to predict accurately, while surface quantities such as heat transfer are more sensitive to the solution procedure. Modeling transition to turbulence needs refinement, though preliminary results are promising.

Edwards, Thomas A.↗
NASA NTRS1989

Small scale structure on cosmic strings

The current understanding of cosmic string evolution is discussed, and the focus placed on the question of small scale structure on strings, where most of the disagreements lie. A physical picture designed to put the role of the small scale structure into more intuitive terms is presented. In this picture it can be seen how the small scale structure can feed back in a major way on the overall scaling solution. It is also argued that it is easy for small scale numerical errors to feed back in just such a way. The intuitive discussion presented here may form the basis for an analytic treatment of the small scale structure, which argued in any case would be extremely valuable in filling the gaps in the present understanding of cosmic string evolution.

Albrecht, Andreas↗
NASA NTRS1989

Dynamic analysis of flexible multibody structures

A system composed of several interconnected elastic components that may experience large angular motion relative to each other during operation is referred to as a flexible multibody structure. Several formulations were proposed for the determination of the dynamic response of controlled flexible multibody structures. In general, these formulations consist of superposing elastic deformations of the component body onto the large rigid body motion of the component. It was shown that this particular methodology for combining linear structural deformations with nonlinear kinematics can lead to erroneous response predictions when either the beam member is very flexible or the rotational speed is high. In addition, previous formulations introduce constraint equations to define the interrelations among system components. This approach increases the number of equations that must be solved, and may result in contraint violation when numerical error accumulates during the integration process. In order to overcome the difficulties, a new approach was suggested. The approach is essentially a finite element formulation which takes advantage of the fact that many multibody structures are joint dominated. The Large Angle Transient Dynamic Analysis (LATDYN) program for clarity of documentation, ease of use, user friendliness, modeling generality, and accuracy of results was evaluated. This required gaining a working familiarity with the code and performing several case studies.

Hernried, Alan G.↗
NASA NTRS1990

Global simulation of the solar wind-magnetosphere interaction - The importance of its numerical validity

This paper examines the global interaction of the solar wind with the geomagnetic field, using a newly developed three-dimensional high-precision MHD simulation code which employs the fourth-order Runge-Kutta-Gill scheme in time and the direct finite-difference method in space, and which allows much less numerical error than the conventional code. The simulation reveals in detail the formation processes of the magnetosphere and the bow shock, as well as the plasma sheet formation process. It is shown that, contrary to the conventional understanding, no Kelvin-Helmholtz instability occurs along the magnetopause in an ideal MHD.

Watanabe, K.↗
NASA NTRS1989

Advanced rotor computations with a corrected potential method

An unsteady Full-Potential Rotor code (FPR) has been enhanced with modifications directed at improving its drag prediction capability. The potential code has been rewritten with modifications to increase the code accuracy. Also, the shock generated entropy has been included to provide solutions comparable to the Euler equations. Two different weakly interacted boundary layer models have also been coupled to FPR in order to estimate skin-friction drag. One is a two-dimensional integral method and the other is a three-dimensional finite-difference scheme. The new flow solver is able to find accurate inviscid drags without recourse to numerical error tares. This permits the resolution of drag distributions resulting from rotor geometric variations. Good comparisons have been obtained between computed and measured torque for a rectangular and a highly swept model rotor.

Bridgeman, John O.↗
NASA NTRS1990

Nonaxisymmetric instabilities in a slender torus - Two- and three-dimensional simulations

Nonaxisymmetric instabilities in accretion disks are investigated in terms of the slender, non-self-gravitating, ideal fluid torus model. Nonlinear simulations are carried out in two and tree dimensions, using Cartesian-grid finite differencing. The fastest growing instability is the principal mode. This mode saturates with the formation of ellipsoidal density distributions, designated as 'planets'. The simulations provide evidence for nonlinear mode-mode coupling. When several wavenumbers m are unstable, multiple 'planets' form that subsequently merge. Sources of numerical error are examined, and the effects are contrasted with those of the physical instability. The same qualitative evolution is seen in both the two- and the three-dimensional simulations, even though strict vertical hydrostatic equilibrium is no longer rigorously maintained after mode saturation.

Hawley, John F.↗
NASA NTRS1991

Optimal controllers for finite wordlength implementation

When a controller is implemented in a digital computer, with A/D and D/A conversion, the numerical errors of the computation can drastically affect the performance of the control system. There exists realizations of a given controller transfer function yielding arbitrarily large effects from computational errors. Since, in general, there is no upper bound, it is important to have a systematic way of reducing these effects. Optimum controller designs are developed which take account of the digital round-off errors in the controller implementation and in the A/D and D/A converters. These results provide a natural extension to the Linear Quadratic Gaussian (LQG) theory since they reduce to the standard LQG controller when infinite precision computation is used. But for finite precision the separation principle does not hold.

Liu, K.↗
NASA NTRS1992

A study of the control problem of the shoot side environment delivery system of a closed crop growth research chamber

The details of our initial study of the control problem of the crop shoot environment of a hypothetical closed crop growth research chamber (CGRC) are presented in this report. The configuration of the CGRC is hypothetical because neither a physical subject nor a design existed at the time the study began, a circumstance which is typical of large scale systems control studies. The basis of the control study is a mathematical model which was judged to adequately mimic the relevant dynamics of the system components considered necessary to provide acceptable realism in the representation. Control of pressure, temperature, and flow rate of the crop shoot environment, along with its oxygen, carbon dioxide, and water concentration is addressed. To account for mass exchange, the group of plants is represented in the model by a source of oxygen, a source of water vapor, and a sink for carbon dioxide. In terms of the thermal energy exchange, the group of plants is represented by a surface with an appropriate temperature. Most of the primitive equations about an experimental operating condition and a state variable representation which was extracted from the linearized equations are presented. Next, we present the results of a real Jordan decomposition and the repositioning of an undesirable eigenvalue via full state feedback. The state variable representation of the modeling system is of the nineteenth order and reflects the eleven control variables and eight system disturbances. Five real eigenvalues are very near zero, with one at zero, three having small magnitude positive values, and one having a small magnitude negative value. A Singular Value Decomposition analysis indicates that these non-zero eigenvalues are not results of numerical error.

Blackwell, C. C.↗
NASA NTRS1992

Computing high-speed flows past an oscillating cylinder near a vertical wall

A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.

Yen, Guan-Wei↗
NASA NTRS1992

Towards industrial-strength Navier-Stokes codes

In this paper we discuss our experiences with Navier-Stokes (NS) codes using central differencing (CD) and scalar artificial dissipation (SAD). The NS-CDSAD codes have been developed by several researchers. Our results confirm that for typical commercial transport wing and wing/body configurations flying at transonic conditions with all turbulent boundary layers, NS-CDSAD codes, when used with the Johnson-King turbulence model, are capable of computing pressure distributions in excellent agreement with experimental data. However, results are not as good when laminar boundary layers are present. Exhaustive 2-D grid refinement studies supported by detailed analysis suggest that the numerical errors associated with SAD severely contaminate the solution in the laminar portion of the boundary layer. It is left as a challenge to the CFD community to find and fix the problems with Navier-Stokes codes and to produce a NS code which converges reliably and properly captures the laminar portion of the boundary layer on a reasonable grid.

Jou, Wen-Huei↗