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.
Richtmyer-Meshkov Instability (RMI) is defined as the growth of irregularities (jets) on the interface of shocked metal with initial sinusoidal perturbations on the free surface. The presentation discusses the use of RMI to study solid media, including jet formation and a tamped environment.
Climate predictions using coupled models in different time scales, from intraseasonal to decadal, are usually affected by initial shocks, drifts, and biases, which reduce the prediction skill. These arise from inconsistencies between different components of the coupled models and from the tendency of the model state to evolve from the prescribed initial conditions toward its own climatology over the course of the prediction. Aiming to provide tools and further insight into the mechanisms responsible for initial shocks, drifts, and biases, this paper presents a novel data set developed within the Long Range Forecast Transient Intercomparison Project, LRFTIP. This data set has been constructed by averaging hindcasts over available prediction years and ensemble members to form a hindcast climatology, that is a function of spatial variables and lead time, and thus results in a useful tool for characterizing and assessing the evolution of errors as well as the physical mechanisms responsible for them. A discussion on such errors at the different time scales is provided along with plausible ways forward in the field of climate predictions
Data for the shock equation of state of helium are obtained up to 360GPa, nearly doubling the pressure of previous experimental measurements. The helium samples are first precompressed to 2.7 gcm −3 in a diamond anvil cell prior to laser-driven shock compression at the Omega Laser Facility. Time-resolved Doppler velocimetry and pyrometry reveal significant reflectivity and greater compressibility compared to the predictions of existing broad-range tabular equation of state models, which may be caused by the onset of ionization. These experimental observations, however, are largely captured with molecular-dynamics simulations based on density functional theory, affirming the ability of first-principles techniques to capture complex physics, while enabling critical insight for the behavior of warm dense helium in Jovian interiors and white dwarf atmospheres.
The lower-hybrid instabilities in high-beta (ratio of plasma to magnetic pressure) plasmas driven by ring-ion distributions in velocity space are investigated. A dispersion equation including electromagnetic effects is derived. In the low-beta limit, analytic expressions are obtained which illuminate the physical nature of the instabilities. The complete dispersion equation is solved numerically as a function of ring speed and plasma beta for several types of ring distribution. Electromagnetic effects are important for relatively energetic rings even in the low-beta regime, suppressing growth rates and shifting the angle of propagation to more oblique angles. Stabilization by thermal effects is also discussed. Application of these results to the earth's bow shock, AMPTE, comets and solar flares is suggested.
General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery, such as inlets, ramjets, and scramjets. The discussion is separated into four areas: (1) computational fluid dynamics models for the entire nonlinear system or high order nonlinear models; (2) high order linearized models derived from fundamental physics; (3) low order linear models obtained from the other high order models; and (4) low order nonlinear models (order here refers to the number of dynamic states). Included in the discussion are any special considerations based on the relevant control system designs. The methods discussed are for the quasi-one-dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, including moving normal shocks, hammershocks, simple subsonic combustion via heat addition, temperature dependent gases, detonations, and thermal choking. The report also contains a comprehensive list of papers and theses generated by this grant.
We have obtained low resolution IUE short wavelength spectra at several positions within a ridge of molecular gas within the Gunn Nebula. Optical absorption line measurements toward HD 62542 indicate moderately high densities within the ridge. The physical appearance of the ridge, along with H-alpha and IRAS images of the nebula, suggest that winds and radiation from Gamma Vel and Zeta Pup are heating and compressing the ridge. The IUE spectra revealed faint C II and C III emission at one position, while C IV was weakly detected at another. The former are interpreted as arising in photoionized gas while the latter probably arises in a shocked wind striking the ridge. No detections were made at other positions in the ridge.
A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.
In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.
Plastic-Bonded Explosive (PBX) prills are agglomerates a few millimeters in diameter composed of High Explosive crystals non-uniformly distributed in and surrounded by a polymer binder. When pressed at elevated temperatures at or above the binder melt temperature the polymer flows to uniformly coat the HE crystals to create a well-consolidated compaction. In reality it does so imperfectly, such that x-ray tomographic scans often look like a collection of prills mashed together. It is speculated that the larger the prill, the more non-uniform the binder distribution, but this has never been proven. Even if that were not so, the larger the prill the farther binder must flow toward the middle to homogenize during pressing. Thus, the degree of homogenization depends in part on prill size. The degree to which binder flows during pressing in turn affects the void distribution within pressed charges, which one suspects will affect shock sensitivity and material strength. In this paper we explore, for ~30 PBX 9502 formulation batches using the same TATB powder lot, how prill size effects the compressive mechanical properties.
Determination of elemental composition and thermophysical properties of materials at high temperatures, as visualized in the context of containerless materials processing in a microgravity environment, presents a variety of unusual requirements owing to the thermal hazards and interferences from electromagnetic control fields. In addition, such information is intended for process control applications and thus the measurements must be real time in nature. A new technique is described which was developed for real time, in-situ determination of the elemental composition of molten metallic alloys such as specialty steel. The technique is based on time-resolved spectroscopy of a laser produced plasma (LPP) plume resulting from the interaction of a giant laser pulse with a material target. The sensitivity and precision were demonstrated to be comparable to, or better than, the conventional methods of analysis which are applicable only to post-mortem specimens sampled from a molten metal pool. The LPP technique can be applied widely to other materials composition analysis applications. The LPP technique is extremely information rich and therefore provides opportunities for extracting other physical properties in addition to the materials composition. The case in point is that it is possible to determine thermophysical properties of the target materials at high temperatures by monitoring generation and transport of acoustic pulses as well as a number of other fluid-dynamic processes triggered by the LPP event. By manipulation of the scaling properties of the laser-matter interaction, many different kinds of flow events, ranging from shock waves to surface waves to flow induced instabilities, can be generated in a controllable manner. Time-resolved detection of these events can lead to such thermophysical quantities as volume and shear viscosities, thermal conductivity, specific heat, mass density, and others.
Who would have thought cosmic rays could be so hip? Although discovered 90 years ago on death-defying manned balloon flights hip even by twenty-first-century extremesport standards cosmic rays quickly lost popularity as way-cool telescopes were finding way-too-cool phenomena across the electromagnetic spectrum. Yet cosmic rays are back in vogue, boasting their own set of superlatives. Scientists are tracking them down with new resolve from the Arctic to Antarctica and even on the high western plains of Argentina. Theorists, too, now see cosmic rays as harbingers of funky physics. Cosmic rays are atomic and subatomic particles - the fastest moving bits of matter in the universe and the only sample of matter we have from outside the solar system (with the exception of interstellar dust grains). Lower-energy cosmic rays come from the Sun. Mid-energy particles come from stellar explosions - either spewed directly from the star like shrapnel, or perhaps accelerated to nearly the speed of light by shock waves. The highest-energy cosmic rays, whose unequivocal existence remains one of astronomy's greatest mysteries, clock in at a staggering 10(exp 19) to 10(exp 22) electron volts. This is the energy carried in a baseball pitch; seeing as how there are as many atomic particles in a baseball as there are baseballs in the Moon, that s one powerful toss. No simple stellar explosion could produce them. At a recent conference in Albuquerque, scientists presented the first observational evidence of a possible origin for the highest-energy variety. A team led by Elihu Boldt at NASA s Goddard Space Flight Center found that five of these very rare cosmic rays (there are only a few dozen confirmed events) come from the direction of four 'retired' quasar host galaxies just above the arm of the Big Dipper, all visible with backyard telescopes: NGC 3610, NGC 3613, NGC 4589, and NGC 5322. These galaxies are billions of years past their glory days as the brightest beacons in the universe. Yet they still harbor central, supermassive black holes, which could generate energetic particles if they are spinning.
The magnetosphere is often thought of as a cavity resonator, with eigenfrequencies dependent on the size of the magnetospheric cavity and the Alfven speed profile. A commonly invoked mechanism to "ring the bell" is a solar wind shock or dynamic pressure discontinuity. The sudden impulse feeds broadband compressional power into the magnetosphere, feeding energy into the cavity eigenmodes. Although the cavity mode model has strong theoretical and simulation support, there are few observations validating the scenario. We present several case studies indicating that damped oscillations observed in the magnetosphere following a shock are clearly directly driven by oscillations already present in the solar wind. We discuss these new observations in the context of the turbulent solar wind model, and discuss recent observations relating these periodicities to the solar surface. We further report on similar oscillations observed in magnetotail flow bursts, and relate them to the solar wind phenomena. We close with a discussion of the role of modulated magnetic reconnection, a potentially universal physical process.
Papers are presented in the areas of the radio characteristics of the quiet sun and active regions, the centimeter, meter and decameter wavelength characteristics of solar bursts, space observations of low-frequency bursts, theoretical interpretations of solar active regions and bursts, joint radio, visual and X-ray observations of active regions and bursts, and the similarities of stellar radio characteristics to solar radio phenomena. Specific topics include the centimeter and millimeter wave characteristics of the quiet sun, radio fluctuations arising upon the transit of shock waves through the transition region, microwave, EUV and X-ray observations of active region loops and filaments, interferometric observations of 35-GHz radio bursts, emission mechanisms for radio bursts, the spatial structure of microwave bursts, observations of type III bursts, the statistics of type I bursts, and the numerical simulation of type III bursts. Attention is also given to the theory of type IV decimeter bursts, Voyager observations of type II and III bursts at kilometric wavelengths, radio and whitelight observations of coronal transients, and the possibility of obtaining radio observations of current sheets on the sun.
Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.
High-Entropy Alloys and Electrospun Nanofiber materials are two classes of novel materials that can offer improved resistance to beam-induced radiation damage and thermal shock. Research to develop these new materials specifically for multi-megawatt accelerator target applications, such as beam windows and particle-production targets, are ongoing at Fermilab within the scope of a DOE Early Career Research Program. The research program combines in-beam experiments with complementary simulations to tailor the microstructures of these novel materials for use in next-generation accelerator target facilities. Iterative simulations to optimize the material composition, physics performance, and beam-induced thermomechanical response will guide the material design and fabrication processes based on established figures of merit. This will be followed by material irradiation experiments using low-energy ions and prototypic high-energy protons with extensive post-irradiation material characterization to assess and qualify the selected novel materials. This talk will describe the alloy design and synthesis, microstructural pre-characterization of the alloys, and plans for the eventual down selection following low-energy ion irradiation studies.
It is pointed out that a finite difference grid which moves with the solution of the partial differential equation being solved can improve the accuracy and efficiency of a numerical algorithm. This technique is particularly advantageous in the solution of problems involving boundary layers or shocks where a poorly chosen grid may give a numerical solution which is useless because of poor resolution or extreme oscillations. The present investigation is concerned with the development of a scheme which does not excessively distort the grid. The grid generation algorithm is based on the numerical solution of a system of elliptic differential equations. Holst and Brown (1981) have used a preliminary solution to move points on the boundary of the physical region and then resolved the problem on a new grid generated by an elliptic system. In the current investigation, the grid movement and the solution will develop simultaneously. The solution is used to modify the generating equations, thereby controlling the grid point distributions.
The Univ. of Maryland portion investigated the following areas. The Space Physics Group performed studies of data from the AMPTE/CCE spacecraft CHEM experiment and found that the ratio of solar wind to photospheric abundances decreased rather smoothly with the first ionization potential (FIP) of the ion with the low FIP ion being about a factor of two overabundant. Carbon and hydrogen fit this trend particularly well. Several occurrences were analyzed of field aligned beams observed when CCE was upstream of the Earth's bow shock. Also using CHEM data, ring current intensity and composition changes during the main and recovery phases of the great geomagnetic storm that occurred in February 1986 was examined in detail. Still using CHEM data, ring current characteristics were examined in a survey of 20 magnetic storms ranging in size from -50 nT to -312 nT. A study was done of energetic ion anisotropy characteristics in the Earth's magnetosheath region using data from the UMD/MPE experiment on ISEE-1. The properties were analyzed of approx. 30 to 130 keV/e protons and alpha particles upstream of six quasi-parallel interplanetary shocks that passed by the ISEE-3 spacecraft during 1978 to 1979. Work from NASA-Goddard include studies from the High Energy Cosmic Ray Group, Low Energy Cosmic Ray Group, Low Energy Gamma Ray Group, High Energy Astrophysics Theory Group, and the X ray Astronomy Group.
A family of high aspect ratio nozzles were designed to provide a parametric database of canonical embedded propulsion concepts. Nozzle throat geometries with aspect ratios of 2:1, 4:1, and 8:1 were chosen, all with convergent nozzle areas. The transition from the typical round duct to the rectangular nozzle was designed very carefully to produce a flow at the nozzle exit that was uniform and free from swirl. Once the basic rectangular nozzles were designed, external features common to embedded propulsion systems were added: extended lower lip (a.k.a. bevel, aft deck), differing sidewalls, and chevrons. For the latter detailed Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations were made to predict the thrust performance and to optimize parameters such as bevel length, and chevron penetration and azimuthal curvature. Seventeen of these nozzles were fabricated at a scale providing a 2.13 inch diameter equivalent area throat." ! The seventeen nozzles were tested for far-field noise and a few data were presented here on the effect of aspect ratio, bevel length, and chevron count and penetration. The sound field of the 2:1 aspect ratio rectangular jet was very nearly axisymmetric, but the 4:1 and 8:1 were not, the noise on their minor axes being louder than the major axes. Adding bevel length increased the noise of these nozzles, especially on their minor axes, both toward the long and short sides of the beveled nozzle. Chevrons were only added to the 2:1 rectangular jet. Adding 4 chevrons per wide side produced some decrease at aft angles, but increased the high frequency noise at right angles to the jet flow. This trend increased with increasing chevron penetration. Doubling the number of chevrons while maintaining their penetration decreased these effects. Empirical models of the parametric effect of these nozzles were constructed and quantify the trends stated above." Because it is the objective of the Supersonics Project that future design work be done more by physics-based computations and less by experiments, several codes under development were evaluated against these test cases. Preliminary results show that the RANS-based code JeNo predicts the spectral directivity of the low aspect ratio jets well, but has no capability to predict the non-axisymmetry. An effort to address this limitations, used in the RANS-based code of Leib and Goldstein, overpredicted the impact of aspect ratio. The broadband shock noise code RISN, also limited to axisymmetric assumptions, did a good job of predicting the spectral directivity of underexpanded 2:1 cold jet case but was not as successful on high aspect ratio jets, particularly when they are hot. All results are preliminary because the underlying CFD has not been validated yet. An effort using a Large Eddy Simulation code by Stanford University predicted noise that agreed with experiments to within a few dB.