Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

Global Acceleration of Coronal Mass Ejections

Using the observed relation between speeds of coronal mass ejections (CMEs) near the Sun and in the solar wind, we estimate a global acceleration acting on the CMEs. Our study quantifies the qualitative results of Gosling [1997] and numerical simulations that CMEs at 1 AU with speeds closer to the solar wind. We found a linear relation between the global acceleration and the initial speed of the CMEs and the absolute value of the acceleration is similar to the slow solar wind acceleration. Our study naturally divides CMEs into fast and slow ones, the dividing line being the solar wind speed. Our results have important implications to space weather prediction models which need to incorporate this effect in estimating the CME arrival time at 1 AU. We show that the arrival times of CMEs at 1 AU are drastically different from the zero acceleration case.

Gopalswamy, Nat↗

Numerical study of cathode emission constraints on cylindrical self-field MPD thruster performance

A stability equation relating thruster discharge current, argon propellant mass flow rate, and electrode geometry has been solved for a variety of cylindrical self-field MPD thruster configurations and discharge currents. Realistic values for cathode emission current densities were used to provide additional constraints on the thruster geometries. A two-dimensional MPD code was employed to provide better estimates of the maximum achievable specific impulse, thrust, and flow efficiency for cases of interest. The model results indicate that long life, cylindrical self-field MPD thrusters operated with argon propellant may not be able to provide specific impulse values in excess of 2100s. Alternate electrode geometries, applied magnetic fields, and/or low molecular weight propellants may be necessary to achieve higher values of specific impulse.

Lapointe, Michael R.↗

Identification of Gravity-Related Effects on Crystal Growth From Melts With an Immiscibility Gap

This work involves an experimental-numerical approach to study the effects of natural and Marangoni convections on solidification of single crystals from a silicate melt with a liquid-liquid immiscibility gap. Industrial use of crystals grown from silicate melts is becoming increasingly important in electronic, optical, and high temperature structural applications. Even the simplest silicate systems like Al203-SiO2 have had, and will continue to have, a significant role in the development of traditional and advanced ceramics. A unique feature of crystals grown from the silicate systems is their outstanding linear electro-optic properties. They also exhibit exceptionally high optical rotativity. As a result, these crystals are attractive materials for dielectric, optical, and microwave applications. Experimental work in our laboratory has indicated that directional solidification of a single crystal mullite appears to be preceded by liquid-liquid phase separation in the melt. Disruption of the immiscible state results in crystallization of a two phase structure. There is also evidence that mixing in the melt caused by density-driven convection can significantly affect the stability of the immiscible liquid layers and result in poly-crystalline growth. On earth, the immiscible state has only been observed for small diameter crystals grown in float zone systems where natural convection is almost negligible. Therefore, it is anticipated that growth of large single crystals from silicate melts would benefit from microgravity conditions because of the reduction of the natural convective mixing. The main objective of this research is to determine the effects of transport processes on the phase separation in the melt during growth of a single crystal while addressing the following issues: (1) When do the immiscible layers form and are they real?; (2) What are the main physical characteristics of the immiscible liquids?; and (3) How mixing by natural or Marangoni convection affects the stability of the phase separated melt.

Kassemi, M.↗

Role of optimization in interdisciplinary analyses of naval structures

The need for numerical design optimization of naval structures is discussed. The complexity of problems that arise due to the significant roles played by three major disciplines, i.e., structural mechanics, acoustics, and hydrodynamics are discussed. A major computer software effort that has recently begun at the David W. Taylor Naval Ship R&D Center to accommodate large multidisciplinary analyses is also described. In addition to primarily facilitating, via the use of data bases, interdisciplinary analyses for predicting the response of the Navy's ships and related structures, this software effort is expected to provide the analyst with a convenient numerical workbench for performing large numbers of analyses that may be necessary for optimizing the design performance. Finally, an example is included that investigates several aspects of optimizing a typical naval structure from the viewpoints of strength, hydrodynamic form, and acoustic characteristics.

Dhir, S. K.↗

Numerical method for the solution of large systems of differential equations of the boundary layer type

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Green, M. J.↗

Method for Approximating the Vacuum Motions of Spinning Symmetrical Bodies with Nonconstant Spin Rates

A method for approximating the vacuum motions of spinning rigid symmetrical bodies with varying spin rates and inertias has been completed. The analysis includes the effects of time varying thrust misalignments, mass unbalance, and jet damping. Results are given in the form of equations for space­ referenced Euler angles, flight-path angles, body­ referenced attitude rates, and earth-referenced vehicle-trajectory coordinates. The method consists of dividing the problem into intervals during which the time-dependent variables are assumed constant at their mean interval value. In order to check this procedure, solutions for various interval sizes are compared with solutions obtained with numerical methods. Although the method is somewhat lengthy for accurate hand computation in most cases, it is readily programed for machine solutions. Probably more important, the general solutions give insight into the separate effects of the variables and, in many cases, can be quickly used to determine the approxi­mate ranges of the variables required for the desired solution to a given problem. In this respect, equations for determining maximum wobble have been derived for certain input conditions. The method has been shown to compare closely with the numerical solutions of two sample problems. The sample problems also illustrated the relatively large effect of pitch and yaw jet damping on body motions.

VACUUM↗

The diurnal heat budget of the thermosphere

Detailed numerical calculations of thermospheric heat sources and sinks are presented and their relative importance is discussed in reference to the energy balance phenomena of the neutral atmosphere. It is shown that the thermal energy available from the absorption in the Schumann-Runge continuum leading to photodissociation of O2 is by far the largest energy source in the lower thermosphere. Other sources of varying importance in different altitude ranges are: (1) energy from photoelectrons, (2) energy exchange from thermal plasma, (3) chemical reactions (ion-electron dissociative recombination) energy gain, and (4) kinetic and dissipative energy associated with the neutral wind. The energy sinks of importance are thermal conduction at the lower boundary (120 km) and radiative cooling of atomic oxygen. It is shown that the combined energy from processes 2 to 4 constitute only a small fraction of the total energy available from photoelectrons and are in phase with the latter. It is suggested that a choice of a lower boundary much below 120 km, e.g. near the mesopause level (90 km), should be more appropriate for general thermospheric studies.

Chandra, S.↗

Comparison of characteristic schemes for three-dimensional, steady, isentropic flow.

The accuracy and order of three bicharacteristic schemes based on the numerical method of characteristics for three-dimensional, steady, isentropic, supersonic flow were compared based on numerical calculations. Comparisons were made for source flows and for flows in contoured axisymmetric nozzles. The first two methods were tetrahedral line networks, differing only in the manner in which the cross derivatives were evaluated. The third method was the pentahedral line network proposed by Butler. All three methods were numerically stable and produced efficient and accurate solutions. A discussion of the relative merits and disadvantages of each method is presented.

Cline, M. C.↗

The diurnal heat budget of the thermosphere.

Detailed numerical calculations of thermospheric heat sources and sinks are presented and their relative importance is discussed in reference to the energy balance phenomena of the neutral atmosphere. It is shown that the thermal energy available from the absorption in the Schumann-Runge continuum leading to photo-dissociation of O2 is by far the largest energy source in the lower thermosphere. Other sources of varying importance in different altitude ranges are: (1) energy from photoelectrons; (2) energy exchange from thermal plasma; (3) chemical reaction (ion-electron dissociative recombination) energy gain; (4) kinetic and dissipative energy associated with the neutral wind. The energy sinks of importance are (1) thermal conduction at the lower boundary (120 km); and (2) radiative cooling of atomic oxygen. It is shown that the combined energy from processes 2 to 4 constitutes only a small fraction of the total energy available from photoelectrons and is in phase with the latter. These secondary sources (processes 2 to 4), therefore, do not constitute a significant energy source and their contribution can be simply incorporated into photoelectron energy (process 1) by defining an effective photoionization heating efficiency.

Chandra, S.↗

Application of ERTS-1 data to the harvest model of the US menhaden fishery

The author has identified the following significant results. The project was conducted in Mississippi Sound in the north-central Gulf of Mexico. It utilized conventional surface data, obtained from fishing and other vessels, as well as aircraft and spacecraft remote data. A relationship was established between surface measured water transparency, temperature and salinity, and commercial fish-stock availability. Numerical models of the relationships were derived. A multiple regression was performed relating ERTS-1 MSS Band 5 image density to measured transparency and water depth. It is concluded that remotely acquired data can play a role in harvest decisions of commercial fisheries.

Maughan, P. M.↗

Finite element shell instability analysis

Formulation procedures and the associated computer program for finite element thin shell instability analysis are discussed. Data cover: (1) formulation of basic element relationships, (2) construction of solution algorithms on both the conceptual and algorithmic levels, and (3) conduction of numerical analyses to verify the accuracy and efficiency of the theory and related programs therein are described.

Source record↗

The intermediate anomaly

Time transformations of the equation dt = cr to the n ds, where s is a variable called the intermediate anomaly, are known to reduce global error in the solution of gravitational systems obtained by numerical integration. Attention is given to the Sundman time transformation, and its relation to equations of Keplerian elliptical motion.

Nacozy, P.↗

Onset of condensation effects as detected by total pressure probes in the Langley 0.3-meter transonic cryogenic tunnel

Total pressure probes mounted in the test section of a 0.3 meter transonic cryogenic tunnel were used to detect the onset of condensation effects for free stream Mach numbers of 0.50, 0.75, 0.85, and 0.95 and for total pressure between one and five atmospheres. The amount of supercooling was found to be about 3 K and suggests that condensation was occurring on pre-existing liquid nitrogen droplets resulting from incomplete evaporation of the liquid nitrogen injected to cool the tunnel. The liquid nitrogen injection process presently being used for the 0.3 m tunnel was found to result in a wide spectrum of droplet sizes being injected into the flow. Since the relatively larger droplets took much more time to evaporate than the more numerous smaller droplets, the larger ones reached the test section first as the tunnel operating temperature was reduced. However, condensation effects in the test section were not immediately measurable because there was not a sufficient number of the larger droplets to have an influence on the thermodynamics of the flow.

Hall, R. M.↗

The electric potential of particles in interstellar space released from a nuclear waste payload

Mechanisms for charging a grain in the interplanetary medium include: (1) capture of solar wind electrons; (2) capture of solar wind protons; (3) ejection of electrons through the photoelectric effect due to the solar radiation; (4) escape of beta particles from beta emitters in the grain; and (5) escape of alpha particles from alpha emitters in the grain. The potentials on both nonradioactive and radioactive grains are considered with relation to particle size and time, and the distance from the Sun. Numerical results are presented where the waste mix is assumed to be PW-4b.

Williams, A. C.↗

Project Galileo mission design

The mission design for Project Galileo is described. Details are presented on the design of the trajectory from launch through completion of the satellite tour, with emphasis on how this design is influenced both by science requirements on the mission and by the interaction of this design with the various elements of the Project. A number of features contribute to the uniqueness of the mission design for this Project. Among the more salient of these that are described in this paper are the use of numerous planetary or satellite encounters for gravity assisted trajectory shaping, and closely related to this, the threading of a trajectory through the Galilean satellite system which meets the mission objectives with a minimum of propulsive trajectory corrections.

Mitchell, R. T.↗

Numerical methods and calculations for droplet flow, heating and ignition

A numerical method was devised and employed to solve a variety of problems related to liquid droplet combustion. The basic transport equations of mass, momentum and energy were formulated in terms of generalized nonorthogonal coordinates, which allows for adaptive griding and arbitrary particle shape. Example problems are solved for internal droplet heating, droplet ignition and high Reynolds number flow over a droplet.

Dwyer, H. A.↗

Estimation of the thickness of the Tharsis lava flows and implications for the nature of the topography of the Tharsis plateau

An attempt to estimate the thickness of the volcanics which cover the large part of the Tharsis region is reported. Estimates for the thickness of the flows were made using the numerous occurrences of the partially buried craters throughout the region and a relation between the crater diameter and crater rim height. Estimates were made along the southern contact with the heavily cratered terrain in the Phoenicis Lacus and Memnonia quadrangles, and along the northeast contact with older plains units in the Tharsis and Lunae Palus quadrangles. In these distal regions the thickness estimates range from 0 to 600-700 m of lava. This would indicate that only 20-30% of the topography is the result of volcanic construction. If the 20-30% estimate is representative of Tharsis as a whole then lava thickness would reach only 2-3 km near the center of the region. Thus only a small fraction of the topography presently observed can be attributed to units readily identifiable as volcanic. This additional evidence indicates that the model for Tharsis of a pre-existing elevated region of heavily cratered crust, capped by a veneer of volcanics is correct.

Plescia, J. B.↗