Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical method”

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 1,531 records · Page 85

A simple numerical orthogonal coordinate generator for fluid dynamic applications

An application of a simple numerical technique which allows for the rapid construction of orthogonal coordinate systems about two dimensional and axisymmetric bodies is presented. The technique is based on a predictor corrector numerical method. It can be used to generate orthogonal meshes which have unequally spaced points in two directions. These orthogonal meshes in their transformed computational plane are, however, equally spaced so that the differencing for the metric coefficients and the fluid dynamic equation terms can be easily determined using equally spaced central finite differences. Solutions to the Navier-Stokes equations for flow over blunt bodies with reverse curvature are presented. The coupling of the time dependent fluid dynamic equations and the coordinate generator worked well with no undesirable effects noted.

Graves, R. A., Jr.↗

Two-dimensional zonal mean flow model

The interactions of dynamics, photochemistry and radiation in the stratosphere can be described by the continuity, momentum, and energy equations. A more comprehensive model for stratospheric transport theory was developed for the purpose of aiding predictions of changes in the stratospheric ozone content as a consequence of natural and anthropogenic processes. This model is time dependent and the dependent variables are zonal means of the relevant meteorological quantities which are functions of latitude and height. The detailed formulation of a numerical model both in physics and mathematics is given. A set of fundamental dynamical equations, and the numerical method used in the integration are described.

Chen, H. C.↗

Numerical solution of compressible viscous flows at high Reynolds numbers

A new numerical method which was used to reduce the computation time required in fluid dynamics to solve the Navier-Stokes equations at flight Reynolds numbers is described. The method is the implicit analogue of the explicit finite different method. It uses this as its first stage, while the second stage removes the restrictive stability condition by recasting the difference equations in an implicit form. The resulting matrix equations to be solved are either upper or lower block bidiagonal equations. The new method makes it possible and practical to calculate many important three dimensional, high Reynolds number flow fields on computers.

Maccormack, R. W.↗

An experimental and analytical investigation of the aerodynamics of a hypersonic vehicle at free stream Mach 6

The goal of this research is the assessment of the validity of existing three dimensional numerical programs in the prediction of the flow fields about general three dimensional hypersonic bodies. A detailed experimental research program was performed in which surface and flow field pressures were mapped. The results of the experimental work were compared with existing inviscid programs. Improvements were made on the existing numerical methods to include angle of attack. A summary of this work is presented.

Zakkay, V.↗

Analysis of starvation effects on hydrodynamic lubrication in nonconforming contacts

Numerical methods were used to determine the effects of lubricant starvation on the minimum film thickness under conditions of a hydrodynamic point contact. Starvation was effected by varying the fluid inlet level. The Reynolds boundary conditions were applied at the cavitation boundary and zero pressure was stipulated at the meniscus or inlet boundary. A minimum-fill-thickness equation as a function of both the ratio of dimensionless load to dimensionless speed and inlet supply level was determined. By comparing the film generated under the starved inlet condition with the film generated from the fully flooded inlet, an expression for the film reduction factor was obtained. Based on this factor a starvation threshold was defined as well as a critically starved inlet. The changes in the inlet pressure buildup due to changing the available lubricant supply are presented in the form of three dimensional isometric plots and also in the form of contour plots.

Brewe, D. E.↗

Hybrid-finite-element analysis of some nonlinear and 3-dimensional problems of engineering fracture mechanics

In this paper, efficient numerical methods for the analysis of crack-closure effects on fatigue-crack-growth-rates, in plane stress situations, and for the solution of stress-intensity factors for arbitrary shaped surface flaws in pressure vessels, are presented. For the former problem, an elastic-plastic finite element procedure valid for the case of finite deformation gradients is developed and crack growth is simulated by the translation of near-crack-tip elements with embedded plastic singularities. For the latter problem, an embedded-elastic-singularity hybrid finite element method, which leads to a direct evaluation of K-factors, is employed.

Atluri, S. N.↗

Large Eddy Simulation of turbulent shear flows

The conceptual foundation underlying Large Eddy Simulation (LES) is summarized, and the numerical methods developed for simulation of the time-developing turbulent mixing layer and turbulent plane Poiseuille flow are discussed. Computational results show that the average Reynolds stress profile nearly attains the equilibrium shape which balances the downstream mean pressure gradient in the regions away from the walls. In the vicinity of the walls, viscous stresses are shown to be significant; together with the Reynolds stresses, these stresses balance the mean pressure gradient. It is stressed that the subgrid scale contribution to the total Reynolds stress is significant only in the vicinity of the walls. The continued development of LES is urged.

Moin, P.↗

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗

Quantitative simulation of a magnetospheric substorm. I - Model logic and overview

Results of a comprehensive computer simulation of the behavior of the earth's inner magnetosphere during a substorm-type event are reported. It is pointed out that the computer model self-consistently computes electric fields, currents, and plasma distributions and velocities in the inner-magnetosphere/ionosphere system; parallel electric fields and ionospheric neutral winds, however, are not included. The basic equations of the model are derived, and the inputs are described. An overview of the results is also given. The first appendix contains derivations of general, useful laws of bounce-averaged gradient, curvature, and E x B drifts in a plasma with isotropic pitch angle distribution. The second appendix gives a description of the numerical method used in the simulation.

Harel, M.↗

Asymptotic solution for the reflectivity of a very rough surface

Asymptotic solutions for the reflectivity of a very rough surface are obtained in closed form up to second order in variance of slope. The zeroth-order solution gives the reflectivity of a specular surface and the sum of the higher-order terms yields the net effect of the rough surface. Within the common regions of validity of small rms slope and incident angle not close to grazing, the result agrees very well with those obtained with numerical methods.

Tsang, L.↗

Hot oxygen atoms in the upper atmosphere of Venus

Using data obtained by an ultraviolet spectrometer on the Pioneer-Venus Orbiter, energy and altitude distribution of nonthermal oxygen atoms in the Venus atmosphere for typical daytime and nighttime conditions are calculated using two different numerical methods. Agreement is found between results obtained by the two-stream transport method and the diffusion equation technique. Dissociative recombination of molecular oxygen ions and charge exchange of atomic oxygen ions with the neutral hydrogen and oxygen gas are the source terms, and calculations show that the dissociative recombination term is the dominant one. It is shown that measured hot atom concentrations exceed estimates of the thermal hydrogen atom concentrations derived from Brinton et al. (1980) over most of the daytime exosphere, and the calculated concentrations do so to an even greater extent. It is concluded that hot oxygen atoms are an important component of the dayside exosphere of Venus, their most important effect being the absorption of solar wind.

Nagy, A. F.↗

Computation of simple three-dimensional wing-vortex interaction in transonic flow

An implicit delta form finite-difference algorithm for Euler equations in conservation law form has been used in preliminary calculations of three-dimensional wing-vortex interactions. Both steady and unsteady transonic flow wing-vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two- and three-dimensions.

Srinivasan, G. R.↗

Techniques for prediction of large nonanalytic reflector antenna performance

Projected spacecraft missions for communications, radio astronomy and earth observations indicate a need for very large (on the order of 100 meters) space deployable or erectable antennas. This paper discusses a numerical method of analyzing reflector characteristics which is applicable to large reflector antennas whose surfaces are subject to mechanical and thermal distortions as well as fabrication and deployment errors. A technique of describing nonanalytic surface configurations by a finite number of known (measured) points was successfully applied in a computer model to predict the RF radiation pattern of a reflector using aperture plane integration techniques. The accuracy of the results is shown to be both a function of the number of measured points used and of their distribution.

Moskowitz, S. M.↗

The design and development of a mixer compound exhaust system for a quiet clean general aviation turbofan

Lower engine fuel consumption has become a dominant factor in turbofan engine design due to rapidly increasing fuel costs. One engine component with a large impact on engine performance is the exhaust system. Previous exhaust system studies have demonstrated the significant exhaust system efficiency gains available through mixing of the core and bypass flows. Typically, a large, costly rig and engine program are required to develop and optimize these gains. The purpose of this paper is to present the results of the low-cost design system used for the quiet, clean general aviation turbofan mixer nozzle design and development. The scale model and full-scale engine test results confirm the predicted 3 to 5% reduction in cruise fuel consumption. This unique design system, which is based on integrating advanced three-dimensional viscous numerical methods with empirical optimization techniques, is summarized and detailed comparisons with test data are presented. The ability to accurately predict relative performance of mixer systems with substantially reduced development time and cost savings is demonstrated.

Blackmore, W. L.↗

Mn-53 and Al-26 evidence for solar cosmic ray constancy - An improved model for interpretation

A general, numerical method for calculating activities of solar cosmic ray produced radionuclides at any point within an irregularly shaped lunar rock of known surface contour and lunar surface orientation is described. This method is then used to predict the activities of Mn-53 andAl-26 as a function of postion within lunar rock 68815 for various assumed values of solar cosmic ray flux (J), rigidity (R sub 0), and rock erosion rate (ER). The predicted activities agree with the measured activities of Kohl et al. (1978) when values of R sub 0 = 100 MV, J = 70 p/sq cm-sec (4 pi, E greater than 10 MeV), ER not greater than 1 mm/m.y. and a total exposure time of 2 m.y. are assumed. These values are in agreement with those found for rocks exposed for not less than 10 m.y. and provide no evidence for variation of the average solar cosmic ray parameters between the last 2 and 10 m.y. intervals. When interpreted with the improved model the activity vs. depth profiles for three faces of 68815 show no evidence of SCR anisotropy or differential erosion.

Russ, G. P., III↗

Higher-level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and the approaches to it are compared. Subgrid scale models and the development of models for the Reynolds-averaged equations are discussed as well as the use of full simulation in testing these models. Numerical methods used in simulating large eddies, the simulation of homogeneous flows, and results from full and large scale eddy simulations of such flows are examined. Free shear flows are considered with emphasis on the mixing layer and wake simulation. Wall-bounded flow (channel flow) and recent work on the boundary layer are also discussed. Applications of large eddy simulation and full simulation in meteorological and environmental contexts are included along with a look at the direction in which work is proceeding and what can be expected from higher-level simulation in the future.

Ferziger, J. H.↗

Spherical models: Basic states and stability

Steady axisymmetric basic states in a spherical cap and their stability to azimuthally varying perturbations of different wave numbers were numerically investigated. The main features of the planned computer code and numerical methods are outlined.

Roberts, G. O.↗

Spectral albedos of midlatitude snowpacks

Spectral albedos of impure-nonhomogeneous snowpacks, typical of midlatitudes, from 400 to 2200 nm were modeled through a numerical solution of the radiative transfer equation in the two-stream approximation. Discrete depth-dependent values of density, grain size and impurity concentration were used to characterize the snowpacks. The model is for diffuse incident radiation, and the numerical method is based on doubling and invariant imbedding. The effect of soot impurities on snowpack albedos is illustrated when a snowpack is several centimeters deep and soot reduces the albedos at visible wavelengths, however, when a snowpack is only a few centimeters deep, soot may increase the albedos at visible wavelengths. By adjusting soot content and snow grain size, good quantitative agreement with some observations at the Cascade Mountains (Washington) and at Point Barrow (Alaska) are obtained; however, the model grain sizes are found to be fifty to four hundred percent larger than the measured values. For satellite snowcover observations, a model for effective albedo of partially snow-covered areas was developed and compared with some NOAA-2 observations of the southeastern United States.

Choudhury, B.↗