Numerical Simulation for Droplet Combustion Using Lagrangian Hydrodynamics

A predictive model of spray combustion must incorporate models for the wide variety of physical environments in a practical combustor. In regions where droplets are closely spaced, combustion resembles a diffusion flame; where they are well separated, an envelope or wake flame results. The relative velocity field between the fuel droplets and oxidizer in influences boundary layer development about the droplet, recirculating flow patterns, and droplet shape and stability. A model must encompass these interacting temporal and spatial effects as well as complicated combustor boundaries. The objective of the current work is to develop the triangular gridding method for describing the individual and collective properties of vaporizing and burning fuel droplets

Numerical simulations of fuel droplet flows using a Lagrangian triangular mesh

The incompressible, Lagrangian, triangular grid code, SPLISH, was converted for the study of flows in and around fuel droplets. This involved developing, testing and incorporating algorithms for surface tension and viscosity. The major features of the Lagrangian method and the algorithms are described. Benchmarks of the algorithms are given. Several calculations are presented for kerosene droplets in air. Finally, extensions which make the code compressible and three dimensional are discussed

Neutron-Diffraction Study of γ-Fe at High Pressure

Face-centered-cubic iron (γ-Fe) coherently precipitated in a single crystal of Cu has been studied by neutron diffraction at high pressures up to 0.9 GPa and in the temperature range from about 5 to 110 K. It undergoes antiferromagnetic ordering at 67±2 K at ambient pressure. The ordering temperature TN monotonically decreases with increasing pressure following a relation TN(K)=67-12.8p - 11.1p2 - 6.17p3, where p is in units of GPa. The change of TN in γ-Fe is discussed in terms of the volume change caused by pressurization or by alloying

Dynamic Screening by Plasticity in Amorphous Solids

In recent work, it was shown that elasticity theory can break down in amorphous solids subjected to nonuniform {\em static} loads. The elastic fields are screened by geometric dipoles; these stem from gradients of the quadrupole field associated with plastic responses. Here we study the dynamical responses induced by {\em oscillatory} loads. The required modification to classical elasticity is described. Exact solutions for the displacement field in circular geometry are presented, demonstrating that dipole screening results in essential departures from the expected predictions of classical elasticity theory. Numerical simulations are conducted to validate the theoretical predictions and to delineate their range of validity.Comment: 9 pages, 7 figure

Auroral Plasma Lines: A First Comparison of Theory and Experiment

In this preliminary report on low-energy (0.3 to 3 eV) secondary electrons in the auroral E layer (90 to 150 km), we compare intensities of plasma lines observed with the Chatanika radar to theoretical predictions obtained from a detailed numerical model. The model calculations are initiated with a flux of energetic auroral primary electrons which enter the atmosphere and lose energy to electrons, ions, and neutrals through a combination of elastic and inelastic collisions. This flux is chosen in order that the total calculated ionization rate matches one that is deduced from the radar measurements. From these same calculations the steady state secondary electron flux is deduced as a function of altitude, energy, and pitch angle. This flux is used to calculate plasma line intensities which are then compared with observed intensities. Initial comparisons suggest that the plasma line theory, when applied to low altitudes, must include the effect of electron-neutral collisions. When this is done, the good agreement obtained between theory and experiment indicates the promise of this approach for the study of low-energy auroral electrons

Jet-Induced Explosions of Core Collapse Supernovae

We numerically studied the explosion of a supernova caused by supersonic jets present in its center. The jets are assumed to be generated by a magneto-rotational mechanism when a stellar core collapses into a neutron star. We simulated the process of the jet propagation through the star, jet breakthrough, and the ejection of the supernova envelope by the lateral shocks generated during jet propagation. The end result of the interaction is a highly nonspherical supernova explosion with two high-velocity jets of material moving in polar directions, and a slower moving, oblate, highly distorted ejecta containing most of the supernova material. The jet-induced explosion is entirely due to the action of the jets on the surrounding star and does not depend on neutrino transport or re-acceleration of a stalled shock. The jet mechanism can explain the observed high polarization of Type Ib,c and Type II supernovae, pulsar kicks, very high velocity material observed in supernova remnants, indications that radioactive material was carried to the hydrogen-rich layers in SN1987A, and some others observations that are very difficult or impossible to explain by the neutrino energy deposition mechanism. The breakout of the jet from a compact, hydrogen- deficient core may account for the gamma-ray bursts and radio outburst associated with SN1998bw/GRB980425.Comment: 14 pages, LaTeX, aaspp4.sty, epsf.sty, submitted to ApJ Let