Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Numerical and laboratory simulation of fault motion and earthquake occurrence

This paper reviews the simulation of earthquake occurrence by numerical and laboratory mechanical block models. Simple linear rheological elements are used with elastic forces driving the main events and viscoelastic forces being important for aftershock and creep occurrence. Friction and its dependence on velocity, stress, and displacement also play a key role in determining how, when, and where fault motion occurs. The discussion of the qualitative behavior of the simulators focuses on the manner in which energy is stored in the system and released by the unstable and stable sliding processes. The numerical results emphasize the statistics of earthquake occurrence and the correlations among source parameters.

Cohen, S. C.↗
NASA NTRS1979

The response of an equatorial ocean to simple wind stress patterns. I - Model formulation and analytic results. II - Numerical results

A time-dependent, primitive equation, beta plane model that is two-dimensional in the horizontal has been developed to model wind-driven equatorial ocean circulation. A simple vertical structure consisting of two layers above the thermocline with the same constant density permits a steady-state undercurrent in the model. An analytical study of the linear dynamics of the model suggests that the addition of inertial effects is needed to simulate the undercurrent properly. Also, both linear and nonlinear dynamics of the model are investigated numerically. Such nonlinear response to wind stress as a strong eastward equatorial undercurrent and an intense eastward 'countercurrent' at three deg N are noted in the numerical results.

Cane, M. A.↗
NASA NTRS1979

Numerical study of an explosion in a non-homogeneous medium with and without magnetic fields

A version of the two step Lax-Wendroff difference method with second order accuracy is used to seek solutions of the unsteady magnetohydrodynamic (MHD) equations for the study of an explosion in a non-homogeneous medium with and without magnetic fields. The explosion is initiated by introducing a finite amount of energy in a small volume of gas which leads to an instantaneous increase of gas temperature. Two types of magnetic field configurations (i.e. open and closed) are considered to illustrate the dependence and differences of magnetic and a non-magnetic flow motion. Numerical results show the strong dependence of induced MHD flow field on the magnetic field configuration and strength. The development and decay of the fast and slow MHD shock waves, as well as the ordinary gasdynamic shock waves, are also represented well by the present computer simulations. In conclusion, we have demonstrated that the numerical scheme we have used is a reliable one for studying 2-dimensional MHD problems in non-homogeneous medium.

Han, S. M.↗
NASA NTRS1978

Numerical solutions of a control problem governed by functional differential equations

A numerical procedure is proposed for solving optimal control problems governed by linear retarded functional differential equations. The procedure is based on the idea of 'averaging approximations', due to Banks and Burns (1975). For illustration, numerical results generated on an IBM 370/158 computer, which demonstrate the rapid convergence of the method are presented.

Banks, H. T.↗
NASA NTRS1979

Numerical solution for the flow field of a body with jet

This paper presents a method of computing the laminar flow field of a body with stern-mounted jet. Problems involved are the numerical solutions of boundary layer and Navier-Stokes equations as well as the viscous-inviscid interaction. The adequacy of a simple viscous-inviscid interaction approach is investigated and a technique for connecting hybrid mesh systems without overlapping of meshes is presented. The numerical results indicate that while the external wake flow affects the development of the jet flow, the presence of a jet also influences the flow in the stern region of the body.

Liu, N.-S.↗
NASA NTRS1978

Application of a trigonometric finite difference procedure to numerical analysis of compressive and shear buckling of orthotropic panels

A numerical analysis developed for the buckling of rectangular orthotropic layered panels under combined shear and compression is described. This analysis uses a central finite difference procedure based on trigonometric functions instead of using the conventional finite differences which are based on polynomial functions. Inasmuch as the buckle mode shape is usually trigonometric in nature, the analysis using trigonometric finite differences can be made to exhibit a much faster convergence rate than that using conventional differences. Also, the trigonometric finite difference procedure leads to difference equations having the same form as conventional finite differences; thereby allowing available conventional finite difference formulations to be converted readily to trigonometric form. For two-dimensional problems, the procedure introduces two numerical parameters into the analysis. Engineering approaches for the selection of these parameters are presented and the analysis procedure is demonstrated by application to several isotropic and orthotropic panel buckling problems. Among these problems is the shear buckling of stiffened isotropic and filamentary composite panels in which the stiffener is broken. Results indicate that a break may degrade the effect of the stiffener to the extent that the panel will not carry much more load than if the stiffener were absent.

Stein, M.↗
NASA NTRS1979

Numerical experiments on the theta pinch

Numerical calculation of theta pinch problems are presented. Physical processes in theta pinch systems are considered in a one dimensional, two temperature magnetohydrodynamic, approximation with allowance for end losses by longitudinal heat conductivity. The numerical calculations are compared with results of earlier experiments.

Volosevich, P. P.↗
NASA NTRS1979

Experimental and numerical results of sound scattering by a body

The interaction of aerodynamic noise with a fuselage shaped body is discussed. A numerical technique is presented which permits the computation of the scattering of an acoustic source by a body at rest for frequencies of aeroacoustic interest. A parallel experiment is described which confirms the results of the computations. A numerical study of varying the geometry of the scattering is presented. In addition, the effect of forward motion on the mean velocity and static pressure profiles in the wake of such a body with a jet exiting from it is simulated. Experimental results are presented and a similarity law is given.

Maestrello, L.↗
NASA NTRS1980

Analysis and testing of numerical formulas for the initial value problem

Three computer programs for evaluating and testing numerical integration formulas used with fixed stepsize programs to solve initial value systems of ordinary differential equations are described. A program written in PASCAL SERIES, takes as input the differential equations and produces a FORTRAN subroutine for the derivatives of the system and for computing the actual solution through recursive power series techniques. Both of these are used by STAN, a FORTRAN program that interactively displays a discrete analog of the Liapunov stability region of any two dimensional subspace of the system. The derivatives may be used by CLMP, a FORTRAN program, to test the fixed stepsize formula against a good numerical result and interactively display the solutions.

Brown, R. L.↗
NASA NTRS1979

Applying integrals of motion to the numerical solution of differential equations

A method is developed for using the integrals of systems of nonlinear, ordinary differential equations in a numerical integration process to control the local errors in these integrals and reduce the global errors of the solution. The method is general and can be applied to either scaler or vector integrals. A number of example problems, with accompanying numerical results, are used to verify the analysis and support the conjecture of global error reduction.

Jezewski, D. J.↗
NASA NTRS1980

Application of a numerically generated orthogonal coordinate system to the solution of inviscid axisymmetric supersonic flow over blunt bodies

A numerically generated orthogonal coordinate system (with the body surface and shock wave as opposite boundaries) was applied with a time asymptotic method to obtain steady flow solutions for axisymmetric inviscid flow over several blunt bodies including spheres, paraboloids, ellipsoids, hyperboloids, hemisphere cylinders, spherically blunted cones, and a body with a concavity in the stagnation region. Comparisons with experimental data and with the results of other computational methods are discussed. The numerically generated orthogonal coordinate system is described and applications of the method to complex body shapes, particularly those with concave regions, are discussed.

Hamilton, H. H., II↗
NASA NTRS1980

Applying integrals of motion to the numerical solution of differential equations

A method is developed for using the integrals of systems of nonlinear, ordinary, differential equations in a numerical integration process to control the local errors in these integrals and reduce the global errors of the solution. The method is general and can be applied to either scalar or vector integrals. A number of example problems, with accompanying numerical results, are used to verify the analysis and support the conjecture of global error reduction.

Vezewski, D. J.↗
NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility demonstrated. Application of the surrogate equation technique allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations or their various approximations. Computational results are presented for the full Euler equations and for the transonic disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗
NASA NTRS1980

Numerical calculation of transonic axial turbomachinery flows

A numerical method and the results of a computer program are presented for solving an exact, three dimensional, full potential equation that models rotating and nonrotating inviscid, absolutely irrotational, homentropic flows. Besides calculating the flows through an arbitrarily shaped rotor or stator blade row mounted on an axisymmetric hub and confined in an axisymmetric duct, the computer program is also capable of analyzing flow fields about arbitrarily shaped wing body combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The governing equation is solved numerically in a fully conservative form by using an artificial time concept, a finite volume technique, rotated type dependent differencing, successive line overrelaxation, and sequential boundary conforming grid refinement. An artificial viscosity is added in fully conservative form, and an initial guess for the potential field is applied, as determined by a two dimensional cascade analysis.

Dulikravich, D. S.↗
NASA NTRS1980

Numerical techniques in linear duct acoustics

Both finite difference and finite element analyses of small amplitude (linear) sound propagation in straight and variable area ducts with flow, as might be found in a typical turboject engine duct, muffler, or industrial ventilation system, are reviewed. Both steady state and transient theories are discussed. Emphasis is placed on the advantages and limitations associated with the various numerical techniques. Examples of practical problems are given for which the numerical techniques have been applied.

Baumeister, K. J.↗
NASA NTRS1980

Numerical analysis of natural convection in two-dimensional square and circular containers in low gravity

A numerical study of natural convection in circular cylinder and square enclosures shows that the analytic low Rayleigh number theory of previous investigators is valid for Rayleigh numbers up to 1000. For a Rayleigh number of 5000, steady state values of maximum fluid velocity differ by 20 percent. This deviation between analytic theory and numerical results increases for higher Rayleigh numbers. In addition, the low Rayleigh number theory is shown to be valid for higher Rayleigh numbers for a portion of the transient phase before significant deviation becomes apparent. It is also shown that square shaped experimental configurations may be analytically approximated with good accuracy by circular cylinders of equal cross sectional area for the prediction of convection velocities and flow patterns at low Rayleigh number.

Robertson, S. J.↗
NASA NTRS1980

Numerical simulation of advection fog formation on multi-disperse aerosols due to combustion-related pollutants

The effects of multi-disperse distribution of the aerosol population are presented. Single component and multi-component aerosol species on the condensation/nucleation processes which affect the reduction in visibility are described. The aerosol population with a high particle concentration provided more favorable conditions for the formation of a denser fog than the aerosol population with a greater particle size distribution when the value of the mass concentration of the aerosols was kept constant. The results were used as numerical predictions of fog formation. Two dimensional observations in horizontal and vertical coordinates, together with time-dependent measurements were needed as initial values for the following physical parameters: (1)wind profiles; (2) temperature profiles; (3) humidity profiles; (4) mass concentration of aerosol particles; (5) particle size distribution of aerosols; and (6) chemical composition of aerosols. Formation and dissipation of advection fog, thus, can be forecasted numerically by introducing initial values obtained from the observations.

Hung, R. J.↗