SEARCH · Search NASA
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.
Drift waves parametric excitation in resistive plasma derived from fluid equations, assuming plane geometry with periodic boundary conditions in magnetic field and azimuthal directions
One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.
A technique is reviewed for calculating thermodynamic quantities for mixtures of light elements at high pressure, in the metallic state. Ensemble averages are calculated with Monte Carlo techniques and periodic boundary conditions. Interparticle potentials are assumed to be coulombic, screened by the electrons in dielectric function theory. This method is quantitatively accurate for alloys at pressures above about 10 Mbar. An alloy of equal parts hydrogen and helium by mass appears to remain liquid and mixed for temperatures above about 3000 K, at pressures of about 15 Mbar. The additive volume law is satisfied to within about 10%, but the Gruneisen equation of state gives poor results. A calculation at 1300 K shows evidence of a hydrogen-helium phase separation.
Package of three computer programs performs two-dimensional electrostatic analysis to determine efficiency of charge transfer. One program can be used to analyze three-electrode charge-coupled device input/output gates, and other two programs can be used to analyze two-phase structures containing two or four electrodes with periodic boundary conditions.
A hybrid plasma simulation model is described and applied to the study of electrostatic wave propagation in a one-dimensional Maxwellian plasma with periodic boundary conditions. The model employs a cloud-in-cell scheme which can drastically reduce the fluctuations in particle simulation models and greatly ease the computational difficulties of the Vlasov equation approach. A grid in velocity space is introduced and the particles are represented by points in the x-v phase space. The model is tested first in the absence of an applied signal and then in the presence of a small-amplitude perturbation. The method is also used to study propagation of an essentially monochromatic plane wave. Results on amplitude oscillation and nonlinear frequency shift are compared with available theories.
The equilibrium properties of a two-dimensional plasma are examined theoretically within the framework of the random-phase approximation. The guiding-center case is of most interest because of the formal equivalence to two-dimensional discrete-line vortex fluids and because of the formal admissibility of negative temperatures. Using a judicious treatment of the periodic boundary conditions, extensive thermodynamics for nonnegative temperatures and a new value for the energy at which the temperature becomes negative are found. An explicit form for the structure function is derived and compared with the result of a Monte Carlo calculation. Good agreement is found for the negative-temperature-threshold prediction and for the general shape of the structure function. An equilibrium energy spectrum is calculated, which for positive temperatures, is the Debye-Hueckel result; for negative temperatures, it predicts an enhanced excitation of the lowest mode, corresponding to macroscopic charge separation.
The three dimensional, time dependent (incompressible) vorticity equations were used to simulate numerically the decay of isotropic box turbulence and time developing mixing layers. The vorticity equations were spatially filtered to define the large scale turbulence field, and the subgrid scale turbulence was modeled. A general method was developed to show numerical conservation of momentum, vorticity, and energy. The terms that arise from filtering the equations were treated (for both periodic boundary conditions and no stress boundary conditions) in a fast and accurate way by using fast Fourier transforms. Use of vorticity as the principal variable is shown to produce results equivalent to those obtained by use of the primitive variable equations.
The unaveraged Navier-Stokes equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three dimensional cosine velocity fluctuations and periodic boundary conditions are used. No mean gradients are present. The three components of the mean square velocity fluctuations are equal for the initial conditions chosen. The resulting solution shows characteristics of turbulence, such as the nonlinear excitation of small scale fluctuations. For the higher Reynolds numbers the initially nonrandom flow develops into an apparently random turbulence.
A procedure was developed using the differential equation approach for generating boundary fitted coordinates for regions with highly curved boundaries as well as reentrant boundaries, such as those encountered in breaking surface waves. The resulting coordinates are nearly orthogonal and can provide adequate resolution even in the reentrant region. Consistent treatment of end boundaries and the use of a systematic initialization scheme and advanced implicit numerical solution techniques make the procedure highly efficient. The method developed for implicit enforcement of the periodicity boundary condition should be beneficial in the analysis of turbomachinery flow applications.
To study the nonlinear physics of incompressible turbulent flow, the unaveraged Navier-Stokes equations are solved numerically. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used. No mean gradients are present. The three components of the mean-square velocity fluctuations are equal for the initial conditions chosen. The resulting solution shows characteristics of turbulence, such as the nonlinear excitation of small-scale fluctuations. For the higher Reynolds numbers the initially nonrandom flow develops into an apparently random turbulence.
A numerical technique for solving the Thomas-Fermi-Dirac (TED) equation in three dimensions, for an array of ions obeying periodic boundary conditions, is presented. The technique is then used to calculate deviations from ideal mixing for an alloy of hydrogen and helium at zero temperature and high presures. Results are compared with alternative models which apply perturbation theory to calculation of the electron distribution, based upon the assumption of weak response of the electron gas to the ions. The TFD theory, which permits strong electron response, always predicts smaller deviations from ideal mixing than would be predicted by perturbation theory. The results indicate that predicted phase separation curves for hydrogen-helium alloys under conditions prevailing in the metallic zones of Jupiter and Saturn are very model dependent.
Rotating baroclinic flow for Richardson number lower than about 1 is studied by means of a finite difference Navier-Stokes model assuming no variations except in the vertical plane that completely contains the density gradient vector. The horizontally infinite channel to which attention is given further assumes periodic boundary conditions at the vertical computational boundaries and no-slip conducting horizontal boundaries. Two configurations are considered. Symmetric baroclinic waves developed in the flows in a manner consistent with linear theory, and it is noted that the structures and energetics of the fully developed waves were dependent on the Prandtl number Pr. For Pr greater than 1, the conversion from potential energy to wave kinetic energy was direct, via temperature and vertical motion correlation, while for Pr of less than 1, conversion proceeded from potential energy to average kinetic energy by means of an induced meridional flow, and then to wave kinetic energy.
Analytic and numerical experiments involving the use of the space-marching method for the prediction of the viscous flow field in turbomachinery cascades indicate that the grid lines in the physical domain must be nearly orthogonal in order to arrive at accurate predictions. For staggered cascades, a nonperiodic grid system was used and appropriate approximations were applied in place of the periodicity boundary conditions upstream and downstream of the cascade. Predicted drag coefficient, turning angle, boundary layer momentum thicknesses, and velocity profiles are compared with experimental data. Good agreement is obtained in most cases.
The initial growth of thin Ge fims on the (100) surface of an Fe substrate is investigated theoretically by means of Monte Carlo simulations based on a potential-energy-scaling technique. The substrate is modeled as a 20 x 20-square array with periodic boundary conditions, as described by Heinbockel et al. (1983), and the movement of surface atoms under the influence of the substrate interaction potential and the lateral interaction of neighboring atoms is explored via continuous updating (on the time scale of single events) of the potential energy at each site in the array. Results for the clustering of nine dispersed atoms over 1.0 s at 600 K and for deposition at 5 x 10 to the -14th/sq cm s over 2.0 s at 500 K are presented graphically.
The three-dimensional, incompressible Navier-Stokes and energy equations with the Bousinesq assumption have been directly simulated at a Rayleigh number of 3.8 x 10 to the 5th power and a Prandtl number of 0.76. In the vertical direction, wall boundaries were used and in the horizontal, periodic boundary conditions were used. A spectral/finite difference numerical method was used to simulate the flow. The flow at these conditions is turbulent and a sufficiently fine mesh was used to capture all relevant flow scales. The results of the simulation are compared to experimental data to justify the conclusion that the small scale motion is adequately resolved.
The incompressible laminar flow around single and double rows of an infinite number of posts between two flat plates has been calculated numerically. A C-grid with periodic boundary conditions has been used. The angle of attack, measured from the line normal to the row of posts, was varied between zero and less than 90 deg. The pressure loading was computed for all of the posts in each of these cases. Most of these calculations have been carried out on the Numerical Aerodynamic Simulation Cray 2 at NASA Ames Research Center.
Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.
The formulation of the Dirac equation on a discrete lattice with half-integer spacing and periodic boundary conditions is investigated analytically. The importance of lattice formulations for problems in field theory and quantum mechanics is explained; the concept of half-integer Fourier representation is introduced; the discrete Dirac equation for the two-dimensional case is derived; dispersion relations for the four-dimensional case are developed; and the spinor formulation for the Dirac fields on the half-integer lattice and the discrete time variable for the four-dimensional time-dependent Dirac equation are obtained. It is argued that the half-integer lattice, because it takes the Dirac Lagrangian into account, is more than a mere relabeling of the integer lattice and may have fundamental physical meaning (e.g., for the statistics of fermions). It is noted that the present formulation does not lead to species doubling, except in the continuum limit.