SEARCH · Engineering Papers

Results for “numerical relativity”

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 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1992

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗

NASA NTRS1994

On the dynamics of a shock-bubble interaction

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Quirk, James J.↗

NASA NTRS1993

Tensoral: A system for post-processing turbulence simulation data

Many computer simulations in engineering and science -- and especially in computational fluid dynamics (CFD) -- produce huge quantities of numerical data. These data are often so large as to make even relatively simple post-processing of this data unwieldy. The data, once computed and quality-assured, is most likely analyzed by only a few people. As a result, much useful numerical data is under-utilized. Since future state-of-the-art simulations will produce even larger datasets, will use more complex flow geometries, and will be performed on more complex supercomputers, data management issues will become increasingly cumbersome. My goal is to provide software which will automate the present and future task of managing and post-processing large turbulence datasets. My research has focused on the development of these software tools -- specifically, through the development of a very high-level language called 'Tensoral'. The ultimate goal of Tensoral is to convert high-level mathematical expressions (tensor algebra, calculus, and statistics) into efficient low-level programs which numerically calculate these expressions given simulation datasets. This approach to the database and post-processing problem has several advantages. Using Tensoral the numerical and data management details of a simulation are shielded from the concerns of the end user. This shielding is carried out without sacrificing post-processor efficiency and robustness. Another advantage of Tensoral is that its very high-level nature lends itself to portability across a wide variety of computing (and supercomputing) platforms. This is especially important considering the rapidity of changes in supercomputing hardware.

Dresselhaus, Eliot↗

NASA NTRS1993

Improvement of method for experimental determination of flutter speed by parameter identification

The method investigated in the current research work is based on the identification of the equations of motion during test flights, followed by the solution of these identified equations to numerically compute the flutter dynamic pressure. The current research work is aimed at overcoming two main difficulties which arise when: (1) A continuous system is truncated into a system with a finite number of degrees of freedom. (2) The desired range of frequencies is wide; in the latter case, numerical difficulties occur that can only partly be overcome by giving relatively more weight to the high frequency modes. Results pertaining to the above first point will be presented under section entitled 'Truncation Effects', and results pertaining to the above second point will be presented under section entitled 'Frequency Range Effects'. The numerical example used herein consists of a continuous simply supported beam. It is further assumed that the object of the identification procedure is to identify the beam's first seven modes.

Nissim, E.↗

NASA NTRS1979

Thermonuclear flashes on accreting neutron stars

Observations of X-ray bursts from binary pulsars and globular clusters are reviewed. The previously proposed hypothesis is considered that such X-ray bursts result from thermonuclear flashes on accreting neutron stars. A general scenario for this mechanism is outlined, and numerical computations of the evolution of the surface layers of an accreting neutron star are discussed. The relation of these calculations to X-ray bursts and other phenomena is examined. Possible improvements in the numerical calculations are suggested.

Joss, P. C.↗

NASA NTRS1982

Numerical computation of exponential matrices using the Cayley-Hamilton theorem

A method for computing exponential matrices, which often arise naturally in the solution of systems of linear differential equations, is developed. An exponential matrix is generated as a linear combination of a finite number (equal to the matrix order) of matrices, the coefficients of which are scalar infinite sums. The method can be generalized to apply to any formal power series of matrices. Attention is focused upon the exponential function, and the matrix exponent is assumed tri-diagonal in form. In such cases, the terms in the coefficient infinite sums can be extracted, as recursion relations, from the characteristic polynomial of the matrix exponent. Two numerical examples are presented in some detail: (1) the three dimensional infinitesimal rotation rate matrix, which is skew symmetric, and (2) an N-dimensional tri-diagonal and symmetric finite difference matrix which arises in the numerical solution of the heat conduction partial differential equation. In the second example, the known eigenvalues and eigenvectors of the finite difference matrix permit an analytical solution for the exponential matrix, through the theory of diagonalization and similarity transformations, which is used for independent verification. The convergence properties of the scalar infinite summations are investigated for finite difference matrices of various orders up to ten, and it is found that the number of terms required for convergence increases slowly with the order of the matrix.

Walden, H.↗

NASA NTRS1988

Numerical methods for turbulent flow

It has generally become accepted that the Navier-Strokes equations predict the dynamic behavior of turbulent as well as laminar flows of a fluid at a point in space away form a discontinuity such as a shock wave. Turbulence is also closely related to the phenomena of non-uniqueness of solutions of the Navier-Strokes equations. These second order, nonlinear partial differential equations can be solved analytically for only a few simple flows. Turbulent flow fields are much to complex to lend themselves to these few analytical methods. Numerical methods, therefore, offer the only possibility of achieving a solution of turbulent flow equations. In spite of recent advances in computer technology, the direct solution, by discrete methods, of the Navier-Strokes equations for turbulent flow fields is today, and in the foreseeable future, impossible. Thus the only economically feasible way to solve practical turbulent flow problems numerically is to use statistically averaged equations governing mean-flow quantities. The objective is to study some recent developments relating to the use of numerical methods to study turbulent flow.

Turner, James C., Jr.↗

NASA NTRS1972

Governing equations for vibrating constrained-layer damping sandwich plates and beams.

A simple differential equation is derived to describe constrained-layer damping in nonsymmetric sandwich plates and beams composed of isotropic and homogeneous layers. The natural boundary conditions related to this equation are determined and some typical numerical results obtained by this equation are given. The equation is valid within the linear theories of elasticity and viscoelasticity in the absence of any constraints on thicknesses, positions, symmetries, and densities of the layers.

Yan, M.-J.↗

NASA NTRS1973

Some services of the Time and Frequency Division of the National Bureau of Standards

The Time and Frequency Division of the National Bureau of Standards (NBS) provides several services to the general public. The radio broadcasts of WWV, WWVH, and WWVB supply reliable, unambiguous time signals to many users. The NBS telephone time-of-day service attracts several hundreds of thousands of calls each year. Periodically, the NBS provides courses on specific topics relating to time and frequency technology. In addition to numerous technical papers published each year, the NBS has prepared the first volume of a comprehensive monograph on time and frequency. The results of research in the Time and Frequency Division of the NBS have had significant impact. An active TV time system capable of serving most of the U.S. currently awaits a ruling by the FCC on a petition filed last year on behalf of the NBS by the Department of Commerce. Three more recent developments are: (1) a TV frequency comparator (patent applied for); (2) a method to perform an independent (absolute) frequency evaluation of commercial cesium beam oscillators; and (3) a method of removing one source of frequency drift in commercial cesium beam oscillators.

Barnes, J. A.↗

NASA NTRS1975

Experimental observation of the ion-ion beam instability in a cylindrical geometry

The ion-ion beam instability is investigated in a cylindrical geometry. In the presence of a cylindrically symmetric ion beam a coherent ion-ion beam instability is observed for a well-defined range of beam velocity. This instability gives rise to cylindrical standing waves whose frequency depends upon beam velocity. The beam distribution function is measured to be a ring in velocity space. The dispersion relation appropriate for the cylindrical standing waves is solved numerically and is shown to be in good agreement with the experimentally observed instability.

Romesser, T.↗

NASA NTRS1976

Investigation of slot cooling at high subsonic speeds

The paper discusses the results of an experimental and numerical investigation of tangential slot injection film cooling with zero pressure gradients in subsonic boundary layers at freestream Mach numbers of 0.4, 0.6, and 0.8. The results are compared with the predictions obtained from a finite-difference boundary-layer program developed by NASA for slot injection into turbulent supersonic boundary layers. The two sets of results are found to compare favorably. The numerical results point to the existence of a unique relation between isothermal effectiveness and adiabatic effectiveness, thereby confirming the existence of similarity conditions for temperature and velocity profiles.

Ferri, A.↗

NASA NTRS1977

Effect of ambient temperature and humidity on emissions of an idling gas turbine

The effects of inlet pressure, temperature, and humidity on the oxides of nitrogen produced by an engine operating at takeoff power setting were investigated and numerous correction factors were formulated. The effect of ambient relative humidity on gas turbine idle emissions was ascertained. Experimentally, a nonvitiating combustor rig was employed to simulate changing combustor inlet conditions as generated by changing ambient conditions. Emissions measurements were made at the combustor exit. For carbon monoxide, a reaction kinetic scheme was applied within each zone of the combustor where initial species concentrations reflected not only local combustor characteristics but also changing ambient conditions.

Kauffman, C. W.↗

NASA NTRS1977

Meteorite impact ejecta - Dependence of mass and energy lost on planetary escape velocity

The amounts of material and energy which escape a planet in a meteorite impact event is calculated as functions of impact and escape velocities. Results are obtained from the computed flow induced by the impact of iron and gabbroic anorthosite spheres onto a half-space of anorthosite at impact velocities of 5 to 45 km/sec. The impact-induced flows were determined by a numerical method using the mass, momentum, and energy conservation relations in finite-difference approximation, within an Eulerian computational grid. The impact velocities at which ejecta losses equal meteorite mass gains are found to be approximately 20, 35, and 45 km/sec for anorthosite objects and approximately 25, 35, and 40 km/sec for iron objects striking anorthosite surfaces for the gravity fields of the moon, Mercury and Mars.

Okeefe, J. D.↗

NASA NTRS1979

Numerical methods for problems involving the Drazin inverse

The objective was to try to develop a useful numerical algorithm for the Drazin inverse and to analyze the numerical aspects of the applications of the Drazin inverse relating to the study of homogeneous Markov chains and systems of linear differential equations with singular coefficient matrices. It is felt that all objectives were accomplished with a measurable degree of success.

Meyer, C. D., Jr.↗

NASA NTRS1981

Vehicle concepts and technology requirements for buoyant heavy-lift systems

Several buoyant-vehicle (airship) concepts proposed for short hauls of heavy payloads are described. Numerous studies identified operating cost and payload capacity advantages relative to existing or proposed heavy-lift helicopters for such vehicles. Applications involving payloads of from 15 tons up to 800 tons were identified. The buoyant quad-rotor concept is discussed in detail, including the history of its development, current estimates of performance and economics, currently perceived technology requirements, and recent research and technology development. It is concluded that the buoyant quad-rotor, and possibly other buoyant vehicle concepts, has the potential of satisfying the market for very heavy vertical lift but that additional research and technology development are necessary. Because of uncertainties in analytical prediction methods and small-scale experimental measurements, there is a strong need for large or full-scale experiments in ground test facilities and, ultimately, with a flight research vehicle.

Ardema, M. D.↗

NASA NTRS1981

Transportation systems evaluation methodology development and applications, phase 3

Transportation systems or proposed changes in current systems are evaluated. Four principal evaluation criteria are incorporated in the process, operating performance characteristics as viewed by potential users, decisions based on the perceived impacts of the system, estimating what is required to reduce the system to practice; and predicting the ability of the concept to attract financial support. A series of matrix multiplications in which the various matrices represent evaluations in a logical sequence of the various discrete steps in a management decision process is used. One or more alternatives are compared with the current situation, and the result provides a numerical rating which determines the desirability of each alternative relative to the norm and to each other. The steps in the decision process are isolated so that contributions of each to the final result are readily analyzed. The ability to protect against bias on the part of the evaluators, and the fact that system parameters which are basically qualitative in nature can be easily included are advantageous.

Kuhlthau, A. R.↗

NASA NTRS1981

Study of Advanced Propulsion Systems for Small Transport Aircraft Technology (STAT) Program

Definitions of takeoff gross weight, performance, and direct operating cost for both a 30 and 50 passenger airplane were established. The results indicate that a potential direct operating cost benefit, resulting from advanced technologies, of approximately 20 percent would be achieved for the 1990 engines. Of the numerous design features that were evaluated, only maintenance-related items contributed to a significant decrease in direct operating cost. Recommendations are made to continue research and technology programs for advanced component and engine development.

Baerst, C. F.↗

NASA NTRS1985

Generation of capillary instabilities by external disturbances in a liquid jet

The receptivity problem in a circular liquid jet is considered. A time harmonic axial pressure gradient is imposed on the steady, parallel flow of a jet of liquid emerging from a circular duct. Using a technique developed in plasma physics a casual solution to the forced problem is obtained over certain ranges of Weber number for a number of mean velocity profiles. This solution contains a term which grows exponentially in the downstream direction and can be identified with a capillary instability wave. Hence, it is found that the externally imposed disturbances can indeed trigger instability waves in a liquid jet. The amplitude of the instability wave generated relative to the amplitude of the forcing is computed numerically for a number of cases.

Leib, S. J.↗