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 235 records · Page 13

Theoretical Foundation for Weld Modeling

Differential equations describe physics of tungsten/inert-gas and plasma-arc welding in aluminum. Report collects and describes necessary theoretical foundation upon which numerical welding model is constructed for tungsten/inert gas or plasma-arc welding in aluminum without keyhole. Governing partial differential equations for flow of heat, metal, and current given, together with boundary conditions relevant to welding process. Numerical estimates for relative importance of various phenomena and required properties of 2219 aluminum included

Traugott, S.↗

Theory of the current-driven ion cyclotron instability in the bottomside ionosphere

A comprehensive treatment of the theory of electrostatic ion cyclotron instability in the collisional bottomside ionosphere is given. The linear dispersion relations for the ion cyclotron instability (electron and ion collisions) are derived. Analytical solutions to the dispersion relation are obtained in three limits: collisionless; weakly collisional; and collisional. In addition, the linear dispersion relation is solved numerically, using parameters typical of the high-latitude bottomside ionosphere. It is pointed out that the linear growth rates in the bottomside collisional regime can be of the same order as the corresponding topside collisionless growth rates. A nonlinear saturation mechanism for the electrostatic ion cyclotron instability is proposed on the basis of Dupree's (1966) resonance broadening theory.

Satyanarayana, P.↗

Conforming versus non-conforming boundary elements in three-dimensional elastostatics

A critical comparison of two basic formulations in three-dimensional elastostatics, using conforming and nonconforming boundary elements, is presented. The basic structure of the boundary element method is developed. The peculiarities that both types of boundary elements present in relation to the numerical implementation are discussed. Through selected examples, key issues such as the computational advantages and disadvantages of both formulations, mesh discretization and accuracy questions, and optimal location of the collocation nodes in the case of nonconforming elements are addressed. It is shown that conforming elements are able to produce more accurate results than nonconforming ones, with substantial economy in the final size of the system equations.

Manolis, G. D.↗

Aspects of supercritical wing theory

Concepts and techniques that have proved successful in the computational theory of shockless airfoils are examined. Shockless airfoils have been used in the design of supercritical wings for fuel-efficient aircraft. In transonic aerodynamics, it is has been found that test data agree well with the best results obtained from CFD. This is particularly true for two-dimensional flow, where accurate experiments are difficult to perform and numerical calculations are relatively easy. However, it is important to choose reasonable models in setting up the computations. Particular attention is given to the formulation of the relevant differential equations in ways that have been found to be satisfactory in practice.

Garabedian, P. R.↗

USSR Space Life Sciences Digest, issue 20

Abstracts of research in the areas of biological rhythms, body fluids, botany, endrocrinology, enzymology, exobiology, genetics, human performance, immunology, life support systems, mathematical modeling, and numerous other topics related to space and life sciences are given.

Hooke, Lydia Razran↗

Discrete events and solar wind energization

Based on a multiple-magnetic-reconnection picture, an estimation of the energy flux suggests that small-scale EUV exploding events may contribute a significant amount of energy (of order of 100,000 erg/sq cm sec) to solar atmospheric heating and solar-wind acceleration. Most of the dissipated magnetic energy is converted into thermal energy and plasma turbulence. On a related aspect, a numerical study based on the nonlinear one-fluid hydrodynamic equations shows a self-smoothing effect, whereby a multistream structure of the solar wind formed near the sun can be gradually smoothed during its propagation through interplanetary space. This calculation gives support for the possible contribution of discrete energetic events to high-speed solar wind streams.

Yang, W.-H.↗

A current disruption mechanism in the neutral sheet for triggering substorm expansions

Two main areas were addressed in support of an effort to understand mechanism responsible for the broadband electrostatic noise (BEN) observed in the magnetotail. The first area concerns the generation of BEN in the boundary layer region of the magnetotail whereas the second area concerns the occassional presence of BEN in the neutral sheet region. For the generation of BEN in the boundary layer region, a hybrid simulation code was developed to perform reliable longtime, quiet, highly resolved simulations of field aligned electron and ion beam flow. The result of the simulation shows that broadband emissions cannot be generated by beam-plasma instability if realistic values of the ion beam parameters are used. The waves generated from beam-plasma instability are highly discrete and are of high frequencies. For the plasma sheet boundary layer condition, the wave frequencies are in the kHz range, which is incompatible with the observation that the peak power in BEN occur in the 10's of Hz range. It was found that the BEN characteristics are more consistent with lower hybrid drift instability. For the occasional presence of BEN in the neutral sheet region, a linear analysis of the kinetic cross-field streaming instability appropriate to the neutral sheet condition just prior to onset of substorm expansion was performed. By solving numerically the dispersion relation, it was found that the instability has a growth time comparable to the onset time scale of substorm onset. The excited waves have a mixed polarization in the lower hybrid frequency range. The imposed drift driving the instability corresponds to unmagnetized ions undergoing current sheet acceleration in the presence of a cross-tail electric field. The required electric field strength is in the 10 mV/m range which is well within the observed electric field values detected in the neutral sheet during substorms. This finding can potentially account for the disruption of cross-tail current and its diversion to the ionosphere to form the substorm current wedge. Furthermore, a number of features associated with substorm expansion onset can be understood based on this substorm onset scenario.

Lui, A. T. Y.↗

The nondeterministic divide

The nondeterministic divide partitions a vector into two non-empty slices by allowing the point of division to be chosen nondeterministically. Support for high-level divide-and-conquer programming provided by the nondeterministic divide is investigated. A diva algorithm is a recursive divide-and-conquer sequential algorithm on one or more vectors of the same range, whose division point for a new pair of recursive calls is chosen nondeterministically before any computation is performed and whose recursive calls are made immediately after the choice of division point; also, access to vector components is only permitted during activations in which the vector parameters have unit length. The notion of diva algorithm is formulated precisely as a diva call, a restricted call on a sequential procedure. Diva calls are proven to be intimately related to associativity. Numerous applications of diva calls are given and strategies are described for translating a diva call into code for a variety of parallel computers. Thus diva algorithms separate logical correctness concerns from implementation concerns.

Charlesworth, Arthur↗

AMPR/SSMI data comparisons

The AMPR (Advanced Microwave Precipitation Radiometer) was flown for the first time with successful data collection over precipitation targets in Florida and off the Oregon coast. The AMPR met its expected performance levels with very low noise and relatively troublefree operation. Numerous rain cloud systems over land and ocean were overflown and the measurements at 10.7, 19.35, 37.1, and 85.5 GHz reveal a wide variety of microphysical conditions which exist within rain cloud systems. Although predicted by radiative transfer model calculations from cloud model simulations, this diversity was not observed before due to the poor spatial resolution of spaceborne microwave radiometers. Saturation of the 19.35 GHz rain emission signal was frequently observed in the oceanic rain systems, supporting the desirability of a 10 GHz channel on the TRMM (Tropical Rain Measuring Mission) microwave radiometer for sensitivity to the higher rain rates.

Spencer, Roy W.↗

A fractional step solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

The time-dependent, three-dimensional incompressible Navier-Stokes equations are presently solved in generalized coordinate systems by means of a fractional-step method whose primitive variable formulation uses as dependent variables, in place of the Cartesian components of the velocity: (1) pressure (defined at the center of the computational cell), and (2) volume fluxes across the faces of the cells. The momentum equations are solved by means of an approximate factorization method. A novel 'ZEBRA' scheme incorporating four-color ordering efficiently solves the Poisson equation. Illustrative two- and three-dimensional laminar flow test cases are computed and evaluated relative to extant numerical and experimental results, and good agreement is obtained.

Rosenfeld, Moshe↗

Autonomous vehicle motion control, approximate maps, and fuzzy logic

Progress on research on the control of actions of autonomous mobile agents using fuzzy logic is presented. The innovations described encompass theoretical and applied developments. At the theoretical level, results of research leading to the combined utilization of conventional artificial planning techniques with fuzzy logic approaches for the control of local motion and perception actions are presented. Also formulations of dynamic programming approaches to optimal control in the context of the analysis of approximate models of the real world are examined. Also a new approach to goal conflict resolution that does not require specification of numerical values representing relative goal importance is reviewed. Applied developments include the introduction of the notion of approximate map. A fuzzy relational database structure for the representation of vague and imprecise information about the robot's environment is proposed. Also the central notions of control point and control structure are discussed.

Ruspini, Enrique H.↗

Integrated Power and Attitude Control System (IPACS) technology developments

Integrated Power and Attitude Control System (IPACS) studies performed over a decade ago established the feasibility of storing electrical energy in flywheels and utilizing the resulting angular momentum for spacecraft attitude control. Such a system has been shown to have numerous attractive features relative to more contemporary technology, and is appropriate to many applications (including high-performance slewing actuators). Technology advances over the last two decades in composite rotors, motor/generator/electronics, and magnetic bearings are found to support the use of IPACS for increasingly sophisticated applications. It is concluded that the concept offers potential performance advantages as well as savings in mass and life-cycle cost. Viewgraphs and discussion on IPACS are included.

Eisenhaure, David B.↗

Discrete ordinates solutions of nongray radiative transfer with diffusely reflecting walls

Nongray gas radiation in a plane parallel slab bounded by gray, diffusely reflecting walls is studied using the discrete ordinates method. The spectral equation of transfer is averaged over a narrow wavenumber interval preserving the spectral correlation effect. The governing equations are derived by considering the history of multiple reflections between two reflecting wails. A closure approximation is applied so that only a finite number of reflections have to be explicitly included. The closure solutions express the physics of the problem to a very high degree and show relatively little error. Numerical solutions are obtained by applying a statistical narrow-band model for gas properties and a discrete ordinates code. The net radiative wail heat fluxes and the radiative source distributions are obtained for different temperature profiles. A zeroth-degree formulation, where no wall reflection is handled explicitly, is sufficient to predict the radiative transfer accurately for most cases considered, when compared with increasingly accurate solutions based on explicitly tracing a larger number of wail reflections without any closure approximation applied.

Menart, J. A.↗

Parallel execution of a three-dimensional, chemically reacting, Navier-Stokes code on distributed-memory machines

This paper describes the parallel version of the three-dimensional, chemically reacting, computational fluid dynamics (CFD) code, SPARK. This work was performed on the Intel iPSC/860-based parallel computers. The SPARK code utilizes relatively simple explicit numerical algorithms, but models complex chemical reactions. The code solves the equations over a regular structured mesh so a simple dam decomposition is used to assign work to the individual processors. The explicit nature of the algorithm, combined with the computational intensity of the chemistry calculations, results in a very low communication-to-computation ratio when compared to typical CFD codes. The efficiency of the parallel code is examined and shown to be about 65 percent when the problem size is scaled with the number of processors. Two low-angle wall-jet injection cases are solved to demonstrate the capability of the parallel code for solving large problems efficiently.

Otto, John C.↗

Application of incremental unknowns to the Burgers equation

In this article, we make a few remarks on the role that attractors and inertial manifolds play in fluid mechanics problems. We then describe the role of incremental unknowns for approximating attractors and inertial manifolds when finite difference multigrid discretizations are used. The relation with direct numerical simulation and large eddy simulation is also mentioned.

Choi, Haecheon↗

Computational chemistry research

Task 41 is composed of two parts: (1) analysis and design studies related to the Numerical Aerodynamic Simulation (NAS) Extended Operating Configuration (EOC) and (2) computational chemistry. During the first half of 1987, Dr. Levin served as a member of an advanced system planning team to establish the requirements, goals, and principal technical characteristics of the NAS EOC. A paper entitled 'Scaling of Data Communications for an Advanced Supercomputer Network' is included. The high temperature transport properties (such as viscosity, thermal conductivity, etc.) of the major constituents of air (oxygen and nitrogen) were correctly determined. The results of prior ab initio computer solutions of the Schroedinger equation were combined with the best available experimental data to obtain complete interaction potentials for both neutral and ion-atom collision partners. These potentials were then used in a computer program to evaluate the collision cross-sections from which the transport properties could be determined. A paper entitled 'High Temperature Transport Properties of Air' is included.

Levin, Eugene↗

An Evaluation of the Iosipescu Specimen for Composite Materials Shear Property Measurement

A detailed evaluation of the suitability of the Iosipescu specimen tested in the modified Wyoming fixture is presented. An experimental investigation using conventional strain gage instrumentation and moire interferometry is performed. A finite element analysis of the Iosipescu shear test for unidirectional and cross-ply composites is used to assess the uniformity of the shear stress field in the vicinity of the notch, and demonstrate the effect of the nonuniform stress field upon the strain gage measurements used for the determination of composite shear moduli. From the test results for graphite-epoxy laminates, it is shown that the proximity of the load introduction point to the test section greatly influences the individual gage readings for certain fiber orientations but the effect upon shear modulus measurement is relatively unimportant. A numerical study of the load contact effect shows the sensitivity of some fiber configurations to the specimen/fixture contact mechanism and may account for the variations in the measured shear moduli. A comparison of the strain gage readings from one surface of a specimen with corresponding data from moire interferometry on the opposite face documented an extreme sensitivity of some fiber orientations to eccentric loading which induced twisting and yielded spurious shear stress-strain curves. In the numerical analysis, it is shown that the Iosipescu specimens for different fiber orientations have to be modeled differently in order to closely approximate the true loading conditions. Correction factors are needed to allow for the nonuniformity of the strain field and the use of the average shear stress in the shear modulus evaluation. The correction factors, which are determined for the region occupied by the strain gage rosette, are found to be dependent upon the material orthotropic ratio and the finite element models. Based upon the experimental and numerical results, recommendations for improving the reliability and accuracy of the shear modulus values are made, and the implications for shear strength measurement discussed. Further application of the Iosipescu shear test to woven fabric composites is presented. The limitations of the traditional strain gage instrumentation on the satin weave and high tow plain weave fabrics is discussed. Test results of a epoxy based aluminum particulate composite is also presented. A modification of the Iosipescu specimen is proposed and investigated experimentally and numerically. It is shown that the proposed new specimen design provides a more uniform shear stress field in the test section and greatly reduces the normal and shear stress concentrations in the vicinity of the notches. While the fabrication and the material cost of the proposed specimen is tremendously reduced, it is shown the accuracy of the shear modulus measurement is not sacrificed.

Ho, Henjen↗

Seal Leakages for Honeycomb or Smooth Configurations

Three dimensional-CFD simulations were attempted to analyze the flow field in a honeycomb flat plate tester. This discussion reviews some of the numerical difficulties and relations, including those relating the honeycombs to labyrinth throttles and the consistency of selected seal and tester data sets.

Athavale, Mahesh M.↗