A Vectorized Algorithm for Molecular Dynamics of Short Range Interacting Particles

We report on a lattice based algorithm, completely vectorized for molecular dynamics simulations. Its algorithmic complexity is of the order O(N), where $N$ is the number of particles. The algorithm works very effectively when the particles have short range interaction, but it is applicable to each kind of interaction. The code was tested on a Cray ymp el in a simulation of flowing granular material.Comment: 9 pages, 6 figures, Late

Advanced electrostatic ion thruster for space propulsion

The suitability of the baseline 30 cm thruster for future space missions was examined. Preliminary design concepts for several advanced thrusters were developed to assess the potential practical difficulties of a new design. Useful methodologies were produced for assessing both planetary and earth orbit missions. Payload performance as a function of propulsion system technology level and cost sensitivity to propulsion system technology level are among the topics assessed. A 50 cm diameter thruster designed to operate with a beam voltage of about 2400 V is suggested to satisfy most of the requirements of future space missions

Towards the Gravity Dual of Quarkonium in the Strongly Coupled QCD Plasma

We build a "bottom-up" holographic model of charmonium by matching the essential spectral data. We argue that this data must include not only the masses but also the decay constants of the J/psi and psi' mesons. Relative to the "soft-wall" models for light mesons, such a matching requires two new features in the holographic potential: an overall upward shift as well as a narrow "dip" near the holographic boundary. We calculate the spectral function as well as the position of the complex singularities (quasinormal frequencies) of the retarded correlator of the charm current at finite temperatures. We further extend this analysis by showing that the residues associated with these singularities are given by the boundary derivative of the appropriately normalized quasinormal mode. We find that the "melting" of the J/psi spectral peak occurs at a temperature of about 540 MeV, or 2.8 T_c, in good agreement with lattice results.Comment: 13 pages, 9 figure

Extended performance solar electric propulsion thrust system study. Volume 4: Thruster technology evaluation

Several thrust system design concepts were evaluated and compared using the specifications of the most advanced 30 cm engineering model thruster as the technology base. Emphasis was placed on relatively high power missions (60 to 100 kW) such as a Halley's comet rendezvous. The extensions in thruster performance required for the Halley's comet mission were defined and alternative thrust system concepts were designed in sufficient detail for comparing mass, efficiency, reliability, structure, and thermal characteristics. Confirmation testing and analysis of thruster and power processing components were performed, and the feasibility of satisfying extended performance requirements was verified. A baseline design was selected from the alternatives considered, and the design analysis and documentation were refined. The baseline thrust system design features modular construction, conventional power processing, and a concentrator solar array concept and is designed to interface with the Space Shuttle

A numerical study of the Navier-Stokes transport coefficients for two-dimensional granular hydrodynamics

A numerical study that aims to analyze the thermal mechanisms of unsteady, supersonic granular flow by means of hydrodynamic simulations of the Navier–Stokes granular equation is reported in this paper. For this purpose, a paradigmatic problem in granular dynamics such as the Faraday instability is selected. Two different approaches for the Navier–Stokes transport coefficients for granular materials are considered, namely the traditional Jenkins–Richman theory for moderately dense quasi-elastic grains and the improved Garzó–Dufty–Lutsko theory for arbitrary inelasticity, which we also present here. Both the solutions are compared with event-driven simulations of the same system under the same conditions, by analyzing the density, temperature and velocity field. Important differences are found between the two approaches, leading to interesting implications. In particular, the heat transfer mechanism coupled to the density gradient, which is a distinctive feature of inelastic granular gases, is responsible for a major discrepancy in the temperature field and hence in the diffusion mechanisms

Coefficient of restitution for elastic disks

We calculate the coefficient of restitution, $\epsilon$, starting from a microscopic model of elastic disks. The theory is shown to agree with the approach of Hertz in the quasistatic limit, but predicts inelastic collisions for finite relative velocities of two approaching disks. The velocity dependence of $\epsilon$ is calculated numerically for a wide range of velocities. The coefficient of restitution furthermore depends on the elastic constants of the material via Poisson's number. The elastic vibrations absorb kinetic energy more effectively for materials with low values of the shear modulus.Comment: 25 pages, 12 Postscript figures, LaTex2

Improving I/O Performance for Exascale Applications through Online Data Layout Reorganization

The applications being developed within the U.S. Exascale Computing Project (ECP) to run on imminent Exascale computers will generate scientific results with unprecedented fidelity and record turn-around time. Many of these codes are based on particle-mesh methods and use advanced algorithms, especially dynamic load-balancing and mesh-refinement, to achieve high performance on Exascale machines. Yet, as such algorithms improve parallel application efficiency, they raise new challenges for I/O logic due to their irregular and dynamic data distributions. Thus, while the enormous data rates of Exascale simulations already challenge existing file system write strategies, the need for efficient read and processing of generated data introduces additional constraints on the data layout strategies that can be used when writing data to secondary storage. We review these I/O challenges and introduce two online data layout reorganization approaches for achieving good tradeoffs between read and write performance. We demonstrate the benefits of using these two approaches for the ECP particle-in-cell simulation WarpX, which serves as a motif for a large class of important Exascale applications. We show that by understanding application I/O patterns and carefully designing data layouts we can increase read performance by more than 80 percent

Measurement of 30-Centimeter Ion Thruster Discharge Cathode Erosion

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76209/1/AIAA-22982-649.pd

Dormant Cathode Erosion in a Multiple-Cathode Gridded Ion Thruster

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76860/1/AIAA-37031-330.pd