SEARCH · Engineering Papers
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.
When a fusion plasma suddenly intercepts a solid surface (wall or pellet), thermal collapse is distinctly kinetic & has novel physics. VPIC simulations & theory revealed the fascinating dynamics of four propagating fronts controlling parallel electron temperature cooling.
Abstract It was recently proposed that the electric field oscillation as a result of self-consistent e ± pair production may be the source of coherent radio emission from pulsars. Direct particle-in-cell simulations of this process have shown that the screening of the parallel electric field by this pair cascade manifests as a limit cycle, as the parallel electric field is recurrently induced when pairs produced in the cascade escape from the gap region. In this work, we develop a simplified time-dependent kinetic model of e ± pair cascades in pulsar magnetospheres that can reproduce the limit-cycle behavior of pair production and electric field screening. This model includes the effects of a magnetospheric current, the escape of e ± , as well as the dynamic dependence of pair production rate on the plasma density and energy. Using this simple theoretical model, we show that the power spectrum of electric field oscillations averaged over many limit cycles is compatible with the observed pulsar radio spectrum.
Homogeneous plane-parallel model atmospheres for solar flares have been constructed to approximately simulate observations of flares. The wings of the Ca II lines have been used to derive flare upper photosphere models, which indicate temperature increases of about 100 K over the temperature distribution in the pre-existing facula at a height of 300 km. In the case of flares covering sunspots the temperature rise seems to occur much higher in the atmosphere. We find that with increasing flare importance the heights of the upper chromosphere and transition region occur lower in the solar atmosphere, accounting for the factor of 60-600 increase in pressure in these regions relative to the quiet sun. The Ca II line profiles agree with observations only by assuming a macrovelocity distribution that increases with height. Also, the chromospheric parts of flares appear to be highly inhomogeneous. We show that shock and particle heated flare models do not agree with the observations and propose a thermal response model for flares. In particular, it appears that heating in the photosphere is an essential aspect of flares.
The high latitude, or downstream flank, and dayside low latitude magnetospheric boundaries are modeled in an MHD simulation of Kelvin-Helmholtz instablities in a compressible plasma for parallel and transverse configurations. Detailed attention is given to the nonlinear consequences of the instabilities for several different Alfven and sound Mach numbers in both configurations. Emphasis is given to the anomalous transport of momentum and energy by the Kelvin-Helmholtz instabilities across the magnetospheric boundary, which is crucial in evaluating the instabilities' contribution to magnetospheric convection. It is concluded that the anomalous tangential stress at the magnetospheric boundary caused by the instability may be of the order of 1 percent of the magnetosheath flow momentum adjacent to the boundary, and gives a contribution to the convection potential drop over the polar cap of 10-30 kV.
An attempt is made to relate the reflection from a broken cloud field to that from a plane parallel cloud of the same optical properties by using Monte Carlo simulations for a very simple model of an array of regularly spaced, infinitely long bar clouds. Reflection from the array is represented by an effective cloud fraction, and it is shown that for many cases this may be obtained by treating the cloud elements to be noninteracting azimuthally symmetric reflectors. The effective cloud fraction is a function of solar zenith angle, the aspect ratio of the individual cloud elements, and the normal or zenith cloud fraction of the array. Model computations for reflection from an array of randomly distributed identical cylinders are also presented as an example of a possible application of this technique.
WAVDRAG calculates supersonic zero-lift wave drag of complex aircraft configurations. Numerical model of aircraft used throughout design process from concept to manufacturing. Incorporates extended geometric input capabilities to permit use of more accurate mathematical model. Engineer defines aircraft components as fusiform or nonfusiform in terms of traditional parallel contours or nonintersecting contours in any direction. Laterally asymmetric configurations simulated. WAVDRAG calculates total drag and wave-drag coefficient of specified aircraft configuration.
A simulation of magnetohydrodynamic Kelvin-Helmholtz instabilities has been performed for parallel and transverse configurations, modeling high latitude (or downstream flanks) and dayside low latitude magnetospheric boundaries. In the parallel configuration, a super-Alfvenic and transsonic shear flow develops into small eddies, which strongly compresses, twists, and hence amplifies the magnetic field by the dynamo action with an amplification factor of M(A)/2. In the nonlinear stage, however large the initial Alfven mach number M(A) may be, the magnetic field amplified and twisted by the hydromagnetic flow vortices reacts back upon the flow evolution, and the flow vortices cascade into smaller structures. In the transverse configuration, the instability leads to the formation of a fast shock discontinuity from an initially subfast shear flow. Anomalous tangential stress by the instability in the transverse configuration reaches 1 percent of the magnetosheath momentum flux, but for the parallel configuration, the anomalous transport is 2-3 times larger than the anomalous transport in the transverse configuration. The anomalous transport for both configurations satisfies the requirement of the viscous-like interaction at the magnetospheric boundary.
Transversely localized double layers evolve randomly in turbulent regions of strongly magnetized plasma carrying current along the magnetic field. Results from numerical simulations and spacecraft observations in the auroral plasma indicate that the parallel electric field in such regions is microscopically intermittent or stochastic. The implications of stochastic double layer fields on electron acceleration will be discussed in terms of a statistical process involving ensemble averages over test particle motion. A Fokker-Planck equation can be derived for the electron phase space density, which depends on the mean and rms amplitudes of the double layers, the mean double layer density, and the initial electron velocity distribution. It is shown that the resulting electron acceleration is very sensitive to the ratio of the initial electron energy to the rms double layer amplitude. When this ratio is large, the acceleration process differs little from that expected in a dc electric field. When it is small, stochastic heating competes with directed acceleration. Evidence for both cases can be found in the auroral ionosphere in association with so-called inverted-V precipitation and collimated edge precipitation.
The use of cinematography or holographic interferometry for dynamic flow visualization in an internal combustion engine requires a control device that globally synchronizes camera and light source timing at a predefined shaft encoder angle. The device is capable of 0.35 deg resolution for rotational speeds of up to 73 240 rpm. This was achieved by implementing the shaft encoder signal addressed look-up table (LUT) and appropriate latches. The developed digital signal processing technique achieves 25 nsec of high speed triggering angle detection by using direct parallel bit comparison of the shaft encoder digital code with a simulated angle reference code, instead of using angle value comparison which involves more complicated computation steps. In order to establish synchronization to an AC reference signal whose magnitude is variant with the rotating speed, a dynamic peak followup synchronization technique has been devised. This method scrutinizes the reference signal and provides the right timing within 40 nsec. Two application examples are described.
The Massively Parallel Processor (MPP) is an ideal machine for computer experiments with simulated neural nets as well as more general cellular automata. Experiments using the MPP with a formal model neural network are described. The results on problem mapping and computational efficiency apply equally well to the neural nets of Hopfield, Hinton et al., and Geman and Geman.
In pursuit of real time animation of computer simulated space plasma phenomena, the code was rewritten for the Massively Parallel Processor (MPP). The program creates a dynamic representation of the global bowshock which is based on actual spacecraft data and designed for three dimensional graphic output. This output consists of time slice sequences which make up the frames of the animation. With the MPP, 16384, 512 or 4 frames can be calculated simultaneously depending upon which characteristic is being computed. The run time was greatly reduced which promotes the rapid sequence of images and makes real time animation a foreseeable goal. The addition of more complex phenomenology in the constructed computer images is now possible and work proceeds to generate these images.
Ion-acoustic solitary waves in a magnetized plasma containing an arbitrary mixture of H(+) and O(+) ions are studied. A nonlinear wave equation has been derived from the Poisson-Vlasov equations, including a uniform magnetic field and dissipation due to reflected electrons. When dissipation is ignored, the equation has soliton solutions associated with both oxygen and hydrogen acoustic modes, which can be either rarefactive or compressive depending on the ion concentrations and the electron/ion temperature ratio and, more weakly, on the bulk drifts of the species. If electron reflection is included, the solitary wave can be intensified. Under somewhat restrictive conditions the oxygen solitary wave is rarefactive and propagates with a velocity comparable to that observed by the Viking satellite. The three-dimensional solitons obey a relation of scales parallel to the magnetic field and in the transverse direction. Computer simulations of one-dimensional versions of the nonlinear wave equation are presented.
Hot flow anomalies (HFAs) are localized plasma structures observed in the solar wind and magnetosheath near the earth's quasi-parallel bow shock. This paper presents one-dimensional hybrid computer simulations illustrating a formation mechanism for HFAs in which the single hot ion population results from a spatial separation of two counterstreaming ion beams. The higher-density cooler regions are dominated by the background (solar wind) ions, and the lower-density hotter internal regions are dominated by the beam ions. The spatial separation of the beam and background is caused by the deflection of the ions in large-amplitude magnetic fields which are generated by ion/ion streaming instabilities.
A three dimensional, finite rate chemistry, Navier-Stokes code was extended to a multi-block code with mismatched interface for practical calculations of supersonic combustors. To ensure global conservation, a conservative algorithm was used for the treatment of mismatched interfaces. The extended code was checked against one test case, i.e., a generic supersonic combustor with transverse fuel injection, examining solution accuracy, convergence, and local mass flux error. After testing, the code was used to simulate the chemically reacting flow fields in a scramjet combustor with parallel fuel injectors (unswept and swept ramps). Computational results were compared with experimental shadowgraph and pressure measurements. Fuel-air mixing characteristics of the unswept and swept ramps were compared and investigated.
A three-dimensional, finite-rate chemistry, Navier-Stokes code has been extended to a multi-block with mismatched interface code for practical calculations of supersonic combustors. To ensure global conservation, a conservative algorithm was used for the treatment of mismatched interfaces. The extended code was checked against one test case, i.e., a generic supersonic combustor with transverse fuel injection, examining solution accuracy, convergence, and local mass flux error. After testing, the code was used to simulate the chemically reacting flow fields in a scramjet combustor with parallel fuel injectors (unswept and swept ramps). Computational results were compared with experimental shadowgraph and pressure measurements. Fuel-air mixing characteristics of the unswept and swept ramps were compared and investigated.
Acceleration control of robotic devices can provide improvements to many space-based operations using flexible manipulators and to ground-based operations requiring better precision and efficiency than current industrial robots can provide. This paper reports on a preliminary study of acceleration measurement on robotic motion during parabolic flights on the NASA KC-135 and a parallel study of accelerations with and without gravity arising from computer simulated motions using TREETOPS software.
The fault-tolerance of analog parallel distributed implementations of a multivariable aircraft neurocontroller is analyzed by simulating weight and neuron failures in a simplified scheme of analog processing based on the functional architecture of the ETANN chip (Electrically Trainable Artificial Neural Network). The neural information processing is found to be only partially distributed throughout the set of weights of the neurocontroller synthesized with the backpropagation algorithm. Although the degree of distribution of the neural processing, and consequently the fault-tolerance of the neurocontroller, could be enhanced using Locally Distributed Weight and Neuron Approaches, a satisfactory level of fault-tolerance could only be obtained by retraining the degrated VLSI neurocontroller. The possibility of maintaining neurocontrol performance and stability in the presence of single weight of neuron failures was demonstrated through an automated retraining procedure of the neurocontroller based on a pre-programmed choice and sequence of the training parameters.
The Space Station/Space Operations Mechanism Test Bed consists of the following: a hydraulically driven, computer controlled Six Degree-of-Freedom Motion System (6DOF); a six degree-of-freedom force and moment sensor; remote driving stations with computer generated or live TV graphics; and a parallel digital processor that performs calculations to support the real time simulation. The function of the Mechanism Test Bed is to test docking and berthing mechanisms for Space Station Freedom and other orbiting space vehicles in a real time, hardware-in-the-loop simulation environment. Typically, the docking and berthing mechanisms are composed of two mating components, one for each vehicle. In the facility, one component is attached to the motion system, while the other component is mounted to the force/moment sensor fixed in the support structure above the 6DOF. The six components of the contact forces/moments acting on the test article and its mating component are measured by the force/moment sensor.