Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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.

At least 487 records · Page 27

Molecular Dynamics Calculations

The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two winning videos in a competition held at the meeting of the American Physical Society's Division of Fluid Dynamics, held in Atlanta, Georgia, in November 1994. Of great interest was the result that in every shock there were a few high-speed precursor particles racing ahead of the shock, carrying information about its impending arrival. Most recently, Dr. Woo has been applying molecular dynamics techniques to the problem of determining the drag produced by the space station truss structure as it flies through the thin residual atmosphere of low-Earth orbit. This problem is made difficult by the complex structure of the truss and by the extreme supersonic nature of the flow. A fully filled section of the truss has already been examined, and drag predictions have been made. Molecular dynamics techniques promise to make realistic drag calculations possible even for very complex partially filled truss segments flying at arbitrary angles.

Source record↗

Theoretical study of the X-ray emission from astrophysical shock waves

Theoretical X-ray emission spectra are needed to interpret the X-ray emission observed by many low and moderate resolution X-ray instruments, and to provide diagnosis of physical conditions for high resolution spectra. Over the past decade, a set of model codes which compute the X-ray and XUV emission for a wide set of physical conditions, including high or low densities, photoionized gas, and time-dependent ionization balance was developed. In the past year, the atomic rate coefficients in the code was improved. Further capabilities were added, and applied to several astrophysical problems.

Raymond, J.↗

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗

The Effect of Shock on the Amorphous Component in Altered Basalt

Investigation of the geochemical and mineralogical composition of the Martian surface provides insight into the geologic history of the predominantly basaltic crust. The Chemistry and Mineralogy (CheMin) instrument onboard the Curiosity rover has returned the first X-Ray diffraction data from the Martian surface. However, large proportions (27 +/- 14 with some estimates as high as 50 weight percentage) of an amorphous component have been reported. As a remedy to this problem, mass balance equations using geochemistry, volatile chemistry, and mineralogy have been employed to constrain the geochemistry of the amorphous component. However, "the nature and number of amorphous phases that constitute the amorphous component is not unequivocally known". Multiple hypotheses have been proposed to explain the origin of this amorphous component: Allophane (Al2O); Basaltic glass (Volcanic and impact); Palagonite (Altered basaltic glass); Hisingerite (Fe (sup 3 plus)-bearing phyllosilicate); S/Cl-rich component (sulfates and/or akaganeite); Nanophase ferric oxide component (npOx). Establishing a multi-phase amorphous component from a basaltic precursor that has undergone physical and chemical weathering within geochemical constraints is of paramount importance to better understand the composition of a large portion of the Martian surface (up to 50 weight percentage). Shocked basalts from Lonar Crater in India are valuable analogs for the Martian surface because it is a well-preserved impact crater in a basaltic target. Having undergone pre- and post-shock aqueous alteration, these rocks provide crucial data regarding the effect of shock on the amorphous component in altered basalt. By conducting mass balance equations similar to what has been performed for Gale crater materials, we attempt to calculate the geochemistry of the amorphous component in altered basalts ranging from unshocked to Class 5 (Table 1). This has the potential to reveal the nature and origin (i.e. primary igneous, shock metamorphic, and/or aqueous alteration occurring before or after the impact event) of the amorphous component in shocked basalt with the goal of unravelling the history of the Martian surface.

Eckley, S. A.↗

Modeling Solar Energetic Particle Events Using ENLIL Heliosphere Simulations

Solar energetic particle (SEP) event modeling (SEPMOD) has gained renewed attention in part because of the availability of a decade of multipoint measurements from STEREO (Solar TErrestrial RElations Observatory) and L1 (Lagrangian point 1) spacecraft at 1 AU (Astronomical Unit). These observations are coupled with improving simulations of the geometry and strength of heliospheric shocks obtained by using coronagraph images to send erupted material into realistic solar wind backgrounds. The STEREO and ACE (Aerosol, Cloud systems, ocean Ecosystems) measurements in particular have highlighted the sometimes surprisingly widespread nature of SEP events. It is thus an opportune time for testing SEP models, which typically focus on protons approximately 1-100 megaelectronvolts, toward both physical insight to these observations and potentially useful space radiation environment forecasting tools. Some approaches emphasize the concept of particle acceleration and propagation from close to the Sun, while others emphasize the local field line connection to a traveling, evolving shock source. Among the latter is the previously introduced SEPMOD treatment, based on the widely accessible and well-exercised WSA-ENLIL (Wang-Sheeley-Arge-ENLIL)-cone model. SEPMOD produces SEP proton time profiles at any location within the ENLIL domain. Here we demonstrate a SEPMOD version that accommodates multiple, concurrent shock sources occurring over periods of several weeks. The results illustrate the importance of considering longer-duration time periods and multiple CME (Coronal Mass Ejection) contributions in analyzing, modeling, and forecasting SEP events.

Luhmann, J. G.↗

The violent interstellar medium

Observational evidence for high-velocity and high-temperature interstellar gas is reviewed. The physical processes that characterize this gas are described, including the ionization and emissivity of coronal gas, the behavior and appearance of high-velocity shocks, and interfaces between coronal gas and cooler interstellar gas. Hydrodynamical models for the action of supernova explosions and stellar winds on the interstellar medium are examined, and recent attempts to synthesize all the processes considered into a global model for the interstellar medium are discussed.

Mccray, R.↗

NGC 7213 - A key to the nature of liners?

New optical spectra of the nearby S0 galaxy NGC 7213 are analyzed. The observations are described and a brief overview of the spectrum is given. The subtraction of the stellar component is presented, and the measurement of various emission lines is discussed in detail. Important physical parameters are derived, including those pertaining to atomic data, line widths, physical conditions, geometry and kinematics, the Balmer decrement, and ultraviolet emission lines. The Seyfert characteristics of NGC 7213 are reviewed, and the galaxy's shock features are reconciled with photoionization models. The presented arguments may apply to classical liners as well. Theoretical models that can help explain the relationships discovered among the emission lines are examined.

Filippenko, A. V.↗

Supersonic propeller noise in a uniform flow

The sound field produced by a supersonic propeller operating in a uniform flow is investigated. The main interest is the effect of the finite forward flight speed on the directivity of the sound field as seen by an observer on the aircraft. It is found that there are cones of silence on the axis of the propeller. The semiapex angles on these cones are equal fore and aft of the propeller plane, and depend on the tip Mach number only. The Fourier coefficients of the acoustic pressure contain the Doppler amplification factor. The sound field weakens in the upstream direction and strengthen downstream. Kinematic considerations of the emitted Mach waves not only confirm these results, but also provide physical insight into the sound generation mechanism. The predicted zone of silence and the Doppler amplification factor are compared to the theoretical prediction of shock wave formation and the flight test of the SR3 propeller.

Jou, Wen-Huei↗

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗

Two- and three-dimensional blade vortex interactions

A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction (BVI) problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The two-dimensional study considered the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The 2-D BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge. The three-dimensional interaction study considered the case of a vortex at 20 deg incidence to the blade leading edge. Although the qualitative results were similar to that of the two-dimensional interaction, details clearly showed the three-dimensional nature of the interaction process.

Davoudzadeh, F.↗

Efficient Multidisciplinary Analysis Procedure Using Multi-Level Parallelization Approach

Multidisciplinary applications are suitable for parallel computing environment by adopting the domain decomposition method. Immediately, a multidisciplinary application can be parallelized by solving each discipline separately. In order to perform coupled multidisciplinary analysis, coupling of each discipline can be accomplished by exchanging boundary data at the interfaces. This is regarded as discipline-level parallelization. Next level could be a "coarse-grain" parallelization of each discipline, which mainly depends on the physical geometry and nature of each discipline. For example, it is almost impossible for structured-grid based computational fluid dynamics codes to do flow analysis of an aircraft by using a single grid because of the complexity of its configuration. Thus, multi-block grid is commonly used to describe the details of complex geometry. Similarly, in structural analysis, the structure is frequently subdivided into substructures. Thus, the computation of each subdomain can be easily parallelized since each subdomain is solved separately independent of other domains. The parallelization is accomplished by solving each subdomain separately on a separate processor and exchanging the boundary conditions at domain interfaces periodically. However, the physical decomposition of the domain introduces explicit boundary conditions at the domain interfaces. This is not desirable for critical areas such as those containing shock waves or flow separations. Thus, a "fine-grain" parallelization is introduced to overcome this problem. The "fine-grain" parallelization is one that solves exactly the same system of equations of a subdomain by using more than one processors without introducing any explicit boundary conditions. An efficient multidisciplinary analysis procedure can be accomplished by successfully combining the above multi-level parallelism. A multidisciplinary analysis code, ENSAERO developed at NASA Ames Research Center is used in this study to implement the proposed approach. The communication data structure required for the proposed approach will be studied in detail. This work will demonstrate the feasibility of using multi-level parallelization approach in multidisciplinary analysis applications.

Byun, Chansup↗

PIV Measurements of Chevrons on F400 Tactical Aircraft Nozzle Model

Previous talks at this meeting have covered our collaborative work on high-energy jets such as present in tactical aircraft (those with supersonic plumes). The emphasis of this work is improving our understanding of flow physics and our prediction tools. In this presentation we will discuss recent flow diagnostics acquired using Particle Image Velocimetry (PIV) made on an underexpanded shocked jet plume from a tactical aircraft nozzle. In this presentation we show cross-sectional and streamwise cuts of both mean and turbulent velocities of an F404 engine nozzle with various chevron designs applied. The impact of chevron penetration, length, and width are documented. The impact of the parameters is generally nonlinear in measures considered here, a surprising result given the relatively smooth behavior of the noise to variations in these chevron parameters.

Bridges, James↗

Expanded Very Large Array Nova Project Observations of the Classical NovaV1723 Aquilae

We present radio light curves and spectra of the classical nova VI723 Aql obtained with the Expanded Very Large Array (EVLA). This is the first paper to showcase results from the EVLA Nova Project, which comprises a team of observers and theorists utilizing the greatly enhanced sensitivity and frequency coverage of EVLA radio observations, along with observations at other wavelengths, to reach a deeper understanding of the energetics, morphology, and temporal characteristics of nova explosions. Our observations of VI723 Aql span 1-37 GHz in frequency, and we report on data from 14 to 175 days following the time of the nova explosion. The broad frequency coverage and frequent monitoring show that the radio behavior of VI723 Aql does not follow the classic Hubble-flow model of homologous spherically expanding thermal ejecta. The spectra are always at least partially optically thin, and the flux rises on faster timescales than can be reproduced with linear expansion. Therefore, any description of the underlying physical processes must go beyond this simple picture. The unusual spectral properties and light curve evolution might be explained by multiple emitting regions or shocked material. Indeed, X-ray observations from Swift reveal that shocks are likely present.

Krauss, Miriam I.↗

Luminosity-dependent Changes of the Cyclotron Line Energy and Spectral Hardness in Cepheus X-4

Context. X-ray spectra of accreting pulsars are generally observed to vary with their X-ray luminosity. In particular, the hardness of the X-ray continuum is found to depend on luminosity. In a few sources, the correlation between the energy of the cyclotron resonance scattering feature (CRSF) and the luminosity is clear. Different types (signs) of the correlation are believed to reflect different accretion modes. Aims. We analyse two NuSTAR observations of the transient accreting pulsar Cep X-4 during its 2014 outburst. Our analysis is focused on a detailed investigation of the dependence of the CRSF energy and of the spectral hardness on X-ray luminosity, especially on short timescales. Methods. To investigate the spectral changes as a function of luminosity within each of the two observations, we used the intrinsic variability of the source on the timescale of individual pulse cycles (tens of seconds), the so-called pulse-to-pulse variability. Results. We find that the NuSTAR spectrum of Cep X-4 contains two CRSFs: the fundamental line at ~30 keV and its harmonic at ~55 keV. We find for the first time that the energy of the fundamental CRSF increases and the continuum becomes harder with increasing X-ray luminosity not only between the two observations, that is, on the long timescale, but also within an individual observation, on the timescale of a few tens of seconds. We investigate these dependencies in detail including their non-linearity. We discuss a possible physical interpretation of the observed behaviour in the frame of a simple one-dimensional model of the polar emitting region with a collisionless shock formed in the infalling plasma near the neutron star surface. With this model, we are able to reproduce the observed variations of the continuum hardness ratio and of the CRSF energy with luminosity.

Vybornov, V.↗

Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations

Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.

CFD↗

Radiation Modeling with Direct Simulation Monte Carlo

Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available rate data. The new method is compared with the current modeling technique and both techniques are compared with available experimental data. The differences in the results are evaluated. The test case is based on experimental measurements from the AVCO-Everett Research Laboratory electric arc-driven shock tube of a normal shock wave in air at 10 km/s and .1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to di erent ow conditions.

Carlson, Ann B.↗

Noise Characteristics of a Four-Jet Impingement Device Inside a Broadband Engine Noise Simulator

The noise generation mechanisms for four directly impinging supersonic jets are investigated employing implicit large eddy simulations with a higher-order accurate weighted essentially non-oscillatory shock-capturing scheme. Impinging jet devices are often used as an experimental apparatus to emulate a broadband noise source. Although such devices have been used in many experiments, a detailed investigation of the noise generation mechanisms has not been conducted before. Thus, the underlying physical mechanisms that are responsible for the generation of sound waves are not well understood. The flow field is highly complex and contains a wide range of temporal and spatial scales relevant for noise generation. Proper orthogonal decomposition of the flow field is utilized to characterize the unsteady nature of the flow field involving unsteady shock oscillations, large coherent turbulent flow structures, and the sporadic appearance of vortex tubes in the center of the impingement region. The causality method based on Lighthill's acoustic analogy is applied to link fluctuations of flow quantities inside the source region to the acoustic pressure in the far field. It will be demonstrated that the entropy fluctuation term in the Lighthill's stress tensor plays a vital role in the noise generation process. Consequently, the understanding of the noise generation mechanisms is employed to develop a reduced-order linear acoustic model of the four-jet impingement device. Finally, three linear acoustic FJID models are used as broadband noise sources inside an engine nacelle and the acoustic scattering results are validated against far-field acoustic experimental data.

Broadband↗

Relating Jet Structure to Photometric Variability: the Herbig Ae Star HD 163296

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) spectra of the jet from the Herbig Ae star HD 163296, obtained with VLT/X-shooter. We determine the physical conditions in the knots and also their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ∼16 yr. The velocity dispersion versus jet velocity and the energy input are comparable between both lobes. However, the mass-loss rate, velocity, and shock conditions are asymmetric. We find ˙M(jet)/ ˙M(acc) is approximately 0.01−0.1, which is consistent with magneto-centrifugal jet launching models. No evidence of any dust is found in the high-velocity jet, suggesting a launch region within the sublimation radius (greater than 0.5 au). The jet inclination measured from proper motions and radial velocities confirms that it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the central source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are very likely caused by dust lifted high above the disk plane, and this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

kinematic launch epochs↗