Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 1,693 records · Page 94

CFD Simulation on the J-2X Engine Exhaust in the Center-Body Diffuser and Spray Chamber at the B-2 Facility

A computational fluid dynamic (CFD) code is used to simulate the J-2X engine exhaust in the center-body diffuser and spray chamber at the Spacecraft Propulsion Facility (B-2). The CFD code is named as the space-time conservation element and solution element (CESE) Euler solver and is very robust at shock capturing. The CESE results are compared with independent analysis results obtained by using the National Combustion Code (NCC) and show excellent agreement.

Wang, Xiao-Yen↗

Realistic Modeling of Interaction of Quiet-Sun Magnetic Fields with the Chromosphere

High­‐resolution observations and 3D MHD simulations reveal intense interaction between the convection zone dynamics and the solar atmosphere on subarcsecond scales. To investigate processes of the dynamical coupling and energy exchange between the subsurface layers and the chromosphere we perform 3D radiative MHD modeling for a computational domain that includes the upper convection zone and the chromosphere, and investigate the structure and dynamics for different intensity of the photospheric magnetic flux. For comparison with observations, the simulation models have been used to calculate synthetic Stokes profiles of various spectral lines. The results show intense energy exchange through small­‐scale magnetized vortex tubes rooted below the photosphere, which provide extra heating of the chromosphere, initiate shock waves, and small­‐scale eruptions.

quiet-sun↗

Dynamics and X-ray emission of a galactic superwind interacting with disk and halo gas

There is a general agreement that the conspicuous extranuclear X-ray, optical-line, and radio-contiuum emission of starbursts is associated with powerful galactic superwinds blowing from their centers. However, despite the significant advances in observational studies of superwinds, there is no consensus on the nature of the emitting material and even on the emission mechanisms themselves. This is to a great extent a consequence of a poor understanding of dynamical processes in the starburst superwind regions. To address this issue, we have conducted two-dimensional hydrodynamical simulations of galactic superwinds. While previous similar studies have used a single (disk) component to represent the ISM of the starburst galaxy, we analyze the interaction of the wind with a two-component disk-halo ambient interstellar medium and argue that this two-component representation is crucial for adequate modeling of starbursts. The emphasis of this study is on the geometry and structure of the wind region and the X-ray emission arising in the wind material and the shocked gas in the disk and the halo of the galaxy. The simulation results have shown that a clear-cut bipolar wind can easily develop under a range of very different conditions. On the other hand, a complex 'filamentary' structure associated with the entrained dense disk material is found to arise within the hot bubble blown out by the wind. The flow pattern within the bubble is dominated equally by the central biconic outflow and a system of whirling motions r elated to the origin and development of the 'filaments'. The filament parameters make them a good candidate for optical-emission-line filamentary gas observed in starburst halos. We find that the history of mass and energy deposition in the starburst region of the galaxy is crucial for wind dynamics. A 'mild' early wind, which arises as a result of the cumulative effect of stellar winds from massive stars, produces a bipolar vertical cavity in the disk and halo gas without strongly affecting the gaseous disk, thus creating conditions for virtually free vertical escape of the hot gas at the later, much more violent supernova-dominated phases of the starburst. We calculate the luminosity, mass, and effective temperature of the X-ray emitting gas in the 'soft' (0.1 to 0.7 keV, 0.7 to 2.2 keV, and 0.1 to 2.2 keV) and 'hard' (1.6 to 8.3 keV) energy bands and estimate the contribution of different gaseous components to the X-ray flux in these bands. Analysis of these parameters enables us to make conclusions regarding the nature of the X-ray-emitting material. We have inferred that the bulk of the soft thermal X-ray emission from starbursts arises in the wind-shocked material of the disk and halo gas rather than in the wind material itself. This enables us to predict that the integrated soft X-ray spectra of starbursts need not show an overabundance of heavy elements which are believed to be produced copiously in the centers of starbursts. Unlike soft X-ray emission, the hard component of thermal X-ray emission is found to originate in the wind material ejected from the starburst region. However, the derived ratio of hard-to-soft X-ray luminosities is too small compared to that observed in starbursts. We conclude therefore that the observed hard X-ray emission of starbursts is probably not associated with the thermal emission of hot wind or ambient shocked gas. Typical temperatures of the bulk of the soft X-ray-emitting material in our very different models have been found to agree well with the ones estimated on the basis of the ROSAT data for the soft component of X-ray emission of nearby starbursts. We predict that temperatures of the extranuclear soft X-ray-emitting gas in starburst galaxies with heavy element abundances near solar should be close to T(sub Xs = 2 to 5 x 10(exp 6)K.

Suchkov, Anatoly A.↗

He(2+) and H(+) dynamics in the subsolar magnetosheath and plasma depletion layer

A magnetosheath interval monitored by a pair of close upstream spacecraft is analyzed. Composition measurements are used to determine separately many of the properties of H(+) and He(2+) in the magnetosheath, plasma depletion layer, and accelerated flow layer. Magnetosheath ion densities and density ratios are well correlated with those in the upstream and on average the He(2+)/H(+) temperature ratio remains constant across the shock. Temperatures decrease and anisotropies increase as the magnetosheath plasma is lost parallel to the magnetic field in the depletion layer. Compared to the magnetosheath, there is a 20-30-percent enhancement in the depletion layer He(2+)/H(+) density ratio and an enhancement in the depletion layer temperature ratio. By monitoring upstream conditions, a methodology is developed for studying the plasma depletion layer and magnetosheath separately.

Fuselier, S. A.↗

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

Shock compression of liquid helium to 360 GPa

Data for the shock equation of state of helium are obtained up to 360⁢G⁡Pa, nearly doubling the pressure of previous experimental measurements. The helium samples are first precompressed to 2.7 g⁡c⁢m −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.

Physics - Plasma physics↗

Atomic Ejecta Source Optical Probe (AESOP) (L3 Milestone Report FY2021 Work Package)

As part of its mission to ensure the reliability of the nation's nuclear weapons stockpile, the NNSA has funded a broad range of projects to characterize the mass of material ejected from shocked metal surfaces with special emphasis placed on determining the size distribution of macroscopic "ejecta" particles. Substantial work has been done for particles with diameters down to roughly 1 μm, but little data is available for smaller particles, and no data exist for the amount of atomic-scale ejecta produced in such dynamic environments. Such data are important due to the implications for weapons systems behavior but also for distinguishing between different potential ejecta production mechanisms, such as Rayleight-Taylor instability or the shallow-bubble-collapse phenomenon recently proposed by G. Maskally.

36 MATERIALS SCIENCE↗

Solar wind flow about the terrestrial planets. I - Modeling bow shock position and shape

A three-parameter method for modeling the position and shape of planetary bow waves was chosen to model the near portion of the Venus, earth and Mars bow shocks, and its results were compared with those of models using one to six free variables. It was found that the relative effective shapes of the near Martian, Cytherean, and terrestrial bow shocks are ellipsoidal, paraboloidal, and hyperboloidal, respectively, in response to the increasing bluntness of the obstacles that the planets present to the solar wind. No significant deviations from axial symmetry were found when the near bow waves of the earth and Venus were mapped into the aberrated terminator plane, in agreement with gas dynamic theory predictions neglecting the effects of the IMF because of their minuteness.

Slavin, J. A.↗

Recent advances in engineering science; Proceedings of the A. Cemal Eringen Symposium, University of California, Berkeley, June 20-22, 1988

Various papers on recent advances in engineering science are presented. Some individual topics addressed include: advances in adaptive methods in computational fluid mechanics, mixtures of two medicomorphic materials, computer tests of rubber elasticity, shear bands in isotropic micropolar elastic materials, nonlinear surface wave and resonator effects in magnetostrictive crystals, simulation of electrically enhanced fibrous filtration, plasticity theory of granular materials, dynamics of viscoelastic media with internal oscillators, postcritical behavior of a cantilever bar, boundary value problems in nonlocal elasticity, stability of flexible structures with random parameters, electromagnetic tornadoes in earth's ionosphere and magnetosphere, helicity fluctuations and the energy cascade in turbulence, mechanics of interfacial zones in bonded materials, propagation of a normal shock in a varying area duct, analytical mechanics of fracture and fatigue.

Koh, Severino L.↗

Thermal instabilities in protogalactic clouds

The means by which a protogalaxy can fragment to form the first generation of stars and globular clusters remains an important problem in astrophysics. Gravitational instabilities grow on timescales too long to drive fragmentation before the background density grows by many orders of magnitude (see Murray and Lin 1989a, and references therein). Thermal instability provides a much more likely mechanism. After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature approx. 10(exp 6) K. Cooling by H and He+ below 10(exp 6) K has a negative slope, so that the cloud is subject to strong thermal instabilities. Density enhancements may then grow rapidly, fragmenting the protogalaxy as it cools to lower temperatures. The role of dynamical effects upon the growth of perturbations is considered here. The method used is similar to that used in Murray and Lin (1989a; see also the Erratum to appear September 15), which examined the growth of thermal instabilities with a one-dimensional Lagrangian hydrodynamics code, written for spherical symmetry. Perturbed regions therefore take the form of shells. The dynamical variables are integrated explicitly, while the temperature, ionization fraction, and molecular fraction are integrated implicitly, and account is taken for non-equilibrium values of these quantities.

Murray, S. D.↗

Vibrational relaxation of anharmonic oscillators in expanding flows

Although the Landau-Teller vibrational model accurately predicts vibrational excitation process in post-shock and compressing flows, it under-predicts the rate of de-excitation in cooling and expanding flows. In the present paper, detailed calculations of the vibrational relaxation process of N2 and CO in cooling flows are conducted. A coupled set of vibrational transition rate equations and quasi-one-dimensional fluid dynamic equations is solved. Multiple quantum level transition rates are computed using SSH theory. The SSH transition rate results are compared with available experimental data and other theoretical models. Vibration-vibration exchange collisions are responsible for some vibrational relaxation acceleration in situations of high vibrational temperature and low translational temperature. The present results support the relaxation mechanisms proposed by Bray and by Treanor Rich and Rehm.

Ruffin, Stephen M.↗

Navier-Stokes, dynamics and aeroelastic computations for vortical flows, buffet and flutter applications

Research on Navier-Stokes, dynamics, and aeroelastic computations for vortical flows, buffet, and flutter applications was performed. Progress during the period from 1 Oct. 1992 to 30 Sep. 1993 is included. Papers on the following topics are included: vertical tail buffet in vortex breakdown flows; simulation of tail buffet using delta wing-vertical tail configuration; shock-vortex interaction over a 65-degree delta wing in transonic flow; supersonic vortex breakdown over a delta wing in transonic flow; and prediction and control of slender wing rock.

Kandil, Osama A.↗

Hypervelocity impact on shielded plates

A ballistic limit equation for hypervelocity impact on thin plates is derived analytically. This equation applies to cases of impulsive impact on a plate that is protected by a multi-shock shield, and it is valid in the range of velocity above 6 km/s. Experimental tests were conducted at the NASA Johnson Space Center on square aluminum plates. Comparing the center deflections of these plates with the theoretical deflections of a rigid-plastic plate subjected to a blast load, one determines the dynamic yield strength of the plate material. The analysis is based on a theory for the expansion of the fragmented projectile and on a simple failure criterion. Curves are presented for the critical projectile radius versus the projectile velocity, and for the critical plate thickness versus the velocity. These curves are in good agreement with curves that have been generated empirically.

Smith, James P.↗

Formulation for Simultaneous Aerodynamic Analysis and Design Optimization

An efficient approach for simultaneous aerodynamic analysis and design optimization is presented. This approach does not require the performance of many flow analyses at each design optimization step, which can be an expensive procedure. Thus, this approach brings us one step closer to meeting the challenge of incorporating computational fluid dynamic codes into gradient-based optimization techniques for aerodynamic design. An adjoint-variable method is introduced to nullify the effect of the increased number of design variables in the problem formulation. The method has been successfully tested on one-dimensional nozzle flow problems, including a sample problem with a normal shock. Implementations of the above algorithm are also presented that incorporate Newton iterations to secure a high-quality flow solution at the end of the design process. Implementations with iterative flow solvers are possible and will be required for large, multidimensional flow problems.

Hou, G. W.↗

MAVRIC Flutter Model Transonic Limit Cycle Oscillation Test

The Models for Aeroelastic Validation Research Involving Computation semi-span wind-tunnel model (MAVRIC-I), a business jet wing-fuselage flutter model, was tested in NASA Langley's Transonic Dynamics Tunnel with the goal of obtaining experimental data suitable for Computational Aeroelasticity code validation at transonic separation onset conditions. This research model is notable for its inexpensive construction and instrumentation installation procedures. Unsteady pressures and wing responses were obtained for three wingtip configurations clean, tipstore, and winglet. Traditional flutter boundaries were measured over the range of M = 0.6 to 0.9 and maps of Limit Cycle Oscillation (LCO) behavior were made in the range of M = 0.85 to 0.95. Effects of dynamic pressure and angle-of-attack were measured. Testing in both R134a heavy gas and air provided unique data on Reynolds number, transition effects, and the effect of speed of sound on LCO behavior. The data set provides excellent code validation test cases for the important class of flow conditions involving shock-induced transonic flow separation onset at low wing angles, including Limit Cycle Oscillation behavior.

Edwards, John W.↗

MAVRIC Flutter Model Transonic Limit Cycle Oscillation Test

The Models for Aeroelastic Validation Research Involving Computation semi-span wind-tunnel model (MAVRIC-I), a business jet wing-fuselage flutter model, was tested in NASA Langley's Transonic Dynamics Tunnel with the goal of obtaining experimental data suitable for Computational Aeroelasticity code validation at transonic separation onset conditions. This research model is notable for its inexpensive construction and instrumentation installation procedures. Unsteady pressures and wing responses were obtained for three wingtip configurations of clean, tipstore, and winglet. Traditional flutter boundaries were measured over the range of M = 0.6 to 0.9 and maps of Limit Cycle Oscillation (LCO) behavior were made in the range of M = 0.85 to 0.95. Effects of dynamic pressure and angle-of-attack were measured. Testing in both R134a heavy gas and air provided unique data on Reynolds number, transition effects, and the effect of speed of sound on LCO behavior. The data set provides excellent code validation test cases for the important class of flow conditions involving shock-induced transonic flow separation onset at low wing angles, including LCO behavior.

Edwards, John W.↗

Impact and Collisional Processes in the Solar System

In the past year, we have successfully developed the techniques necessary to conduct impact experiments on ice at very low temperatures. We employ the method of embedding gauges within a target to measure the shock wave and material properties. This means that our data are not model dependent; we directly measure the essential parameters needed for numerical simulations of impact cratering. Since then we have developed a new method for temperature control of icy targets that ensures temperature equilibrium throughout a porous target. Graduate student, Sarah Stewart-Mukhopadhyay, is leading the work on ices and porous materials as the main thrust of her thesis research. Our previous work has focused on icy materials with no porosity, and we propose to extend our research to include porous ice and porous ice-silicate mixtures. There is little shockwave data for porous ice, and none of the data was acquired under conditions applicable to the outer solar system. The solid ice Hugoniot is only defined for initial temperatures above -20 C. Our program uniquely measures the properties of ice at temperatures directly applicable to the solar system. Previous experiments were conducted at ambient temperatures soon after removing the target from a cold environment, usually just below freezing, or in a room just below freezing. Since ice has an extremely complicated phase diagram, it is important to conduct experiments at lower temperatures to determine the true outcome of impacts in the outer solar system. This research is complementary to other programs on icy materials. Our work focuses on the inherent material properties by measuring the shock wave directly; this complements the macroscopic observations and immediately provides the parameters necessary to extend this research to the gravity regime. Our numerical simulations of impacts in porous ice under very low gravity conditions, such as found on comets, show that the final crater size and shape is very dependent on the dynamic strength of the material.

Ahrens, Thomas J.↗

New York City Panel on Climate Change 2019 Report Chapter 9: Perspectives on a City in a Changing Climate 2008-2018

Cities experience multiple environmental shifts, stresses, and shocks—such as air and water pollution—and a variety of extreme events simultaneously and continuously. Current urban programs have focused on limiting the impacts of these conditions through a portfolio of multifaceted strategies, such as regulations and codes, management and restoration projects, and citizen engagement. Global climate change represents a new environmental dynamic to which cities now have to respond. While global climate change by definition has impacts world wide, residents and managers of cities, like New York, typically perceive changes in their own local environments. In most cities, temperature is warming with increasingly hotter and longer heatwaves, and heavier downpours are leading to more frequent inland flooding. In coastal cities, sea levels are rising, exacerbating coastal flooding. Analyzing and understanding the impacts of climate change on cities is important because of the dramatic growth in urban populations throughout the world. An estimated nearly 4.0 billion people reside in urban areas, accounting for 52% of the world’s population (UN, 2017). That percentage will increase dramatically in the coming decades as almost all of the growth to take place up to 2050 will be in urban areas (UN, 2017). The New York City metropolitan region (NYMR)—the five boroughs (equivalent to counties) of New York City and the adjacent 26 counties in the states of New York, New Jersey, and Connecticut—is an ideal model of an urban agglomeration. Approximately 8.6 million people live in the five boroughs and more than 15 million people live in the neighboring smaller cities, towns, and villages (City of New York, 2018a; US Census, 2017). The population of the five boroughs is projected to add 1 million people by 2030, while the total region is projected to reach 26.1 million (NYTC, 2015).

Solecki, William↗