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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

SWIFT Code Assessment for Two Similar Transonic Compressors

One goal of the NASA Fundamental Aeronautics Program is the assessment of computational fluid dynamic (CFD) codes used for the design and analysis of many aerospace systems. This paper describes the assessment of the SWIFT turbomachinery analysis code for two similar transonic compressors, NASA rotor 37 and stage 35. The two rotors have identical blade profiles on the front, transonic half of the blade but rotor 37 has more camber aft of the shock. Thus the two rotors have the same shock structure and choking flow but rotor 37 produces a higher pressure ratio. The two compressors and experimental data are described here briefly. Rotor 37 was also used for test cases organized by ASME, IGTI, and AGARD in 1994-1998. Most of the participating codes over predicted pressure and temperature ratios, and failed to predict certain features of the downstream flowfield. Since then the AUSM+ upwind scheme and the k- turbulence model have been added to SWIFT. In this work the new capabilities were assessed for the two compressors. Comparisons were made with overall performance maps and spanwise profiles of several aerodynamic parameters. The results for rotor 37 were in much better agreement with the experimental data than the original blind test case results although there were still some discrepancies. The results for stage 35 were in very good agreement with the data. The results for rotor 37 were very sensitive to turbulence model parameters but the results for stage 35 were not. Comparison of the rotor solutions showed that the main difference between the two rotors was not blade camber as expected, but shock/boundary layer interaction on the casing.

Chima, Rodrick V.↗
NASA NTRS2011

Passage of Bolides Through the Atmosphere

Different fragmentation models are applied to a number of events, including the entry of TC3 2008 asteroid in order to reproduce existing observational data. Keywords meteoroid entry fragmentation modeling 1 Introduction Fragmentation is a very important phenomenon which occurs during the meteoroid entry into the atmosphere and adds more drastic effects than mere deceleration and ablation. Modeling of bolide fragmentation (100 106 kg in mass) may be divided into several approaches. Detail fitting of observational data (deceleration and/or light curves) allows the determination of some meteoroid parameters (ablation and shape-density coefficients, fragmentation points, amount of mass loss) (Ceplecha et al. 1993; Ceplecha and ReVelle 2005). Observational data with high accuracy are needed for the gross-fragmentation model (Ceplecha et al. 1993), which is used for the analysis of European and Desert bolide networks data. Hydrodynamical models, which describe the entry of the meteoroid including evolution of its material, are applied mainly for large bodies (>106 kg) (Boslough et al. 1994; Svetsov et al. 1995; Shuvalov and Artemieva 2002, and others). Numerous papers were devoted to the application of standard equations for large meteoroid entry in the attempts to reproduce dynamics and/or radiation for different bolides and to predict meteorite falls. These modeling efforts are often supplemented by different fragmentation models (Baldwin and Sheaffer, 1971; Borovi.ka et al. 1998; Artemieva and Shuvalov, 2001; Bland and Artemieva, 2006, and others). The fragmentation may occur in different ways. For example, few large fragments are formed. These pieces initially interact through their shock waves and then continue their flight independently. The progressive fragmentation model suggests that meteoroids are disrupted into fragments, which continue their flight as independent bodies and may be disrupted further. Similar models were suggested in numerous papers, beginning with Levin (1956) and initial interaction of fragments started to be taken into account after the paper by Passey and Melosh (1980). The progressive fragmentation model with lateral spreading of formed fragments is widely used (Artemieva and Shuvalov, 1996; Nemtchinov and Popova, 1997; Borovi.ka et al. 1998; Bland and Artemieva, 2006).

Popova, O.↗
NASA NTRS2002

Photoevaporation of Clumps in Photodissociation Regions

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo external heating and photodissociation as they are exposed to the FUV field, resulting in a loss of cold, molecular lump mass as it is converted to warm atomic gas. The heating, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to three dimensionless parameters: Nc0, the ratio of the initial column density of the clump to the column N(0) approx. 10(exp 21) cm(exp -2) of a warm FUV-heated surface region; upsilon, the ratio of the sound speed in the heated surface to that in the cold clump material: and t(FUV)t(c), the ratio of the "turn-on time" t(FUV) of the heating flux on a clump to its initial sound crossing-time t(c). The evolution also depends on whether a confining interclump medium exists, or whether the interclump region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the dependence of clump photoevaporation on the physical parameters of the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photo evaporation timescales, approx. 10(exp 4-5) years for a 1 M(solar mass) clump in a typical star-forming region (Nc0 = 10, upsilon = 10). Clumps with insufficient magnetic pressure support, and in strong FUV fields may be driven to collapse by the compressional effect of converging shock waves. We also estimate the rocket effect on photoevaporating clumps and find that it is significant only for the smallest clumps, with sizes much less than the extent of the PDR itself. Clumps that are confined by all interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that, absorbs the incident FUV flux. We compare our results with observations of some well studied PDRs: the Orion Bar, M17SW NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin. turbulence and pressure confinement, with a slight indication for favouring the turbulent model for clumps over pressure-confined

Gorti, Uma↗
NASA NTRS1973

Dynamic response of Mach 2.5 axisymmetric inlet with 40 percent supersonic internal area contraction

Results of experimental tests conducted on a supersonic, mixed compression, axisymmetric inlet are presented. The inlet is designed for operation at Mach 2.5 with a turbofan engine (TF-30). The inlet was terminated with either a choked-orifice plate or a long pipe with variable area choked exit plug. Frequency responses were obtained for selected static pressures in the diffuser. These pressures were selected as potential control signals for terminal shock control. Frequency responses were obtained for the Mach 2 and 2.5 conditions for different terminations. Responses also were obtained with and without cowl bleed. Internal disturbances were produced by sinusoidally varying the inlet overboard bypass doors at frequencies out to 100 hertz.

Baumbick, R. J.↗
NASA NTRS1981

On the utility of finite element theory for computational fluid dynamics

An implicit finite element numerical solution algorithm is derived for the compressible Navier-Stokes equations expressed in generalized coordinates. The theoretical basis utilizes a Galerkin-Weighted Residuals formulation, and extremization of approximation error within the context of a multipole expansion. A von Neumann analysis for a simplified form indicates the algorithm fourth- to sixth-order phase accurate, with third-order dissipation for the elementary linear element construction. Performance is improved for the algorithm constructed using quadratic interpolation. Numerical experiments for shocked duct flows are employed to optimize the several algorithm parameters. Additional numerical solutions validate algorithm accuracy and utility for aerodynamics applications.

Baker, A. J.↗
NASA NTRS1992

The dynamics and control of fluctuating pressure loads in the reattachment region of a supersonic free shear layer

An experimental program to investigate the control of a turbulent, reattaching shear layer at Mach 2.9 is described. In preliminary experiments, schlieren photography and Rayleigh scattering were used to visualize the effects of air injection normal to the plane of the shear layer. Localized blowing was qualitatively found to increase the growth rate of the shear layer, the intensity of the turbulence, and the unsteadiness of the reattachment shock. However, uniformly distributed air injection did not appear to affect the flow strongly. These results suggest that three-dimensional disturbances are more effective for flow control than two-dimensional disturbances, and that there is a strong connection between the incoming turbulence and the shock motion in this flow.

Poggie, J.↗
DOE OSTI2026

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗
DOE OSTI2026

New Perspectives on Vibrational Energy Transfer in Energetic Materials: Insights from Pressure-Tuned Ultrafast Spectroscopy

Understanding the manner in which vibrational energy flows between molecular and lattice vibrations is of great interest in physical chemistry due to its central role in reactivity and energy dissipation in molecular materials. Here, in this feature article, we highlight our recent efforts employing ultrafast broadband infrared spectroscopy toward understanding the interplay between molecular and lattice vibrations in energetic materials, motivated by the open questions surrounding the role of vibrational energy transfer (VET) in reaction initiation in these materials. Our work addresses the ongoing debate on the participation of doorway modes in VET. We further present new results from high-pressure ultrafast experiments on RDX, a hydrogen-bonded material, and BNFF, a hydrogen-free material, to explore how intermolecular interaction strength governs VET pathways and time scales. Collectively, our findings reveal that vibrational dynamics in these systems occurs across three distinct time regimes, with VET being incomplete out to hundreds of picoseconds, suggesting the importance of considering nonstatistical reactions in the modeling of these materials. These time scales vary as intermolecular interaction strength is indirectly modified by application of static pressure, indicating dramatic changes to the vibrational structure of these materials under shock-relevant conditions. Our results thus shed light on how intermolecular interactions shape vibrational energy redistribution in molecular materials, and highlight the need for further theoretical and experimental investigation.

Carlson, Daniel Ryan [Sandia National Laboratories↗
NASA NTRS1994

Tearing mode instability in a multiple current sheet system

The tearing mode and magnetic reconnection are studied for multiple current sheet systems by two-dimensional magnetohydrodynamic (MHD) simulations. Both the linear and nonlinear evolution of this process are anaylsed for laminar perturbations. The results illustrate the existence of a linear regime with a symmetric and antisymmetric mode and agree with previous analytic results (Otto and Birk, 1992). The nonlinear evolution shows a number of interesting new features and may explain some properties in corresponding studies of turbulent reconnection. For wavelengths larger than twice the current sheet separation the evolution of antisymmetric modes leads to an entire reconfiguration of the magnetic field and converts a major portion of the magnetic energy into kinetic energy. Antisymmetric modes with smaller wavelengths and symmetric modes are found to saturate. The influence of the value of the resistivity on the reconnection rate decreases in the nonlinear evolution, and the ratio of current sheet separation to wavelength seems to be of major importance. A comparion of the dynamics of periodic current sheets with the evolution of only two current sheets indicates that some of the results for the periodic system also apply to the evolution of only two interacting current sheets. The results are discussed with respect to observations of large-scale plasma and magnetic field reconfigurations in the magnetosheath and near the Earth's bow shock.

Yan, M.↗
NASA NTRS1998

Multi-Dimensional High Order Essentially Non-Oscillatory Finite Difference Methods in Generalized Coordinates

This project is about the development of high order, non-oscillatory type schemes for computational fluid dynamics. Algorithm analysis, implementation, and applications are performed. Collaborations with NASA scientists have been carried out to ensure that the research is relevant to NASA objectives. The combination of ENO finite difference method with spectral method in two space dimension is considered, jointly with Cai [3]. The resulting scheme behaves nicely for the two dimensional test problems with or without shocks. Jointly with Cai and Gottlieb, we have also considered one-sided filters for spectral approximations to discontinuous functions [2]. We proved theoretically the existence of filters to recover spectral accuracy up to the discontinuity. We also constructed such filters for practical calculations.

Shu, Chi-Wang↗
NASA NTRS2001

Hypervelocity Capability of the HYPULSE Shock-Expansion Tunnel for Scramjet Testing

New hypervelocity capabilities for scramjet testing have recently been demonstrated in the HYPULSE Shock-Expansion Tunnel (SET). With NASA's continuing interests in scramjet testing at hypervelocity conditions (Mach 12 and above), a SET nozzle was designed and added to the HYPULSE facility. Results of tests conducted to establish SET operational conditions and facility nozzle calibration are presented and discussed for a Mach 15 (M15) flight enthalpy. The measurements and detailed computational fluid dynamics calculations (CFD) show the nozzle delivers a test gas with sufficiently wide core size to be suitable for free-jet testing of scramjet engine models of similar scale as, those tested in conventional low Mach number blow-down test facilities.

Foelsche, Robert O.↗
NASA NTRS2016

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗
NASA NTRS2013

Silent and Efficient Supersonic Bi-Directional Flying Wing

The supersonic bi-directional (SBiDir) flying wing (FW) concept has a great potential to achieve low sonic boom with high supersonic aerodynamic performance due to removal of performance conflict between high speed and low speed by rotating goo in flight. This NIAC Phase 1 research has achieved three objectives: 2) prove the concept based on simulation that it can achieve very low boom with smooth Sine wave ground over-pressure signature and excellent aerodynamic efficiency; 3) conduct trade study to correlate the geometric parameters with sonic boom and aerodynamic performance for further automated design optimization in Phase II. The design methodology developed in Phase I includes three parts: 1) an advanced geometry model, which can vary airfoil meanline angle distribution to control the expansion and shock waves on the airplane surface to mitigate sonic boom and improve aerodynamic efficiency. 2) a validated CFD procedure to resolve near field flow with accurate shock strength. The sonic boom propagation from near field to far field ground is simulated by NASA NF Boom code. The surface friction drag prediction is based on fiat plate correlation adopted by Seebass and supported by the experimental study of Winter and Smith, which is on the conservative side and is more reliable than CFD RANS simulation. 3) a mission analysis tool based on Corke's model that provides design requirements and constraints of supersonic airplanes for range, payload, volume, size, weight, etc. The design mission target is a supersonic transport with cruise Mach number 1.6, 100 passengers, and 4000nm range. The trade study has several very important findings: 1) The far field ground sonic boom signature is directly related to the smoothness of the flow on the airplane surface. The meanline angle distribution is a very effective control methodology to mitigate surface shock and expansion wave strength, and mitigating compression wave coalescing by achieving smooth loading distribution chord-wise. Compared with a linear meanline angle distribution, a design using nonlinear and non-monotonic meanline angle distribution is able to reduce the sonic boom ground loudness by over 20dBP1. The design achieves sonic boom ground loudness less than 70dBP1 and aerodynamic dynamic efficiency 1/D of 8.4. 2) Decreasing sweep angle within the Mach cone will increase 1/D as well as sonic boom. A design with variable sweep from 84° at the very leading edge to 68° at the tip achieves an extraordinarily high 1/D of 10.4 at Mach number 1.6 due to the low wave drag. If no sonic boom constraint is attached, SBiDir-FW concept still has a lot of room to increase the 1/D. 3) The round leading edge and trailing edge under high sweep angle are beneficial to improve aerodynamic performance, sonic boom, and to increase volume of the airplane. 4) Subsonic performance is benefited greatly from the high slenderness of supersonic configuration after rotating goo. A design with excellent supersonic aspect ratio of 0.44, 1/D of 8.g, gives an extraordinary subsonic aspect ration of 10 and 1/D of 1g.7. Two configurations are designed in details to install internal seats, landing gears, and engine installation to demonstrate the feasibility of SBiDir-FW configuration to accommodate all the required volume for realistic airplane. Here we emphasize that the qualitative findings in Phase I are very encouraging, more important than the quantitative results. Qualitative findings give the understanding of physics and provide the path to achieve the ultimate high performance design. The promising quantitative results achieved in Phase I need to be confirmed by wind tunnel testing in Phase II and ultimately proved by flight test. The other important step forward will be made to study the rotation transition from both CFD unsteady simulation and wind tunnel testing.

Zha, Gecheng↗
NASA NTRS1979

Fluid dynamic problems associated with air-breathing propulsive systems

A brief account of research activities on problems related to air-breathing propulsion is made in this final report for the step funded research grant NASA NGL 14-005-140. Problems include the aircraft ejector-nozzle propulsive system, nonconstant pressure jet mixing process, recompression and reattachment of turbulent free shear layer, supersonic turbulent base pressure, low speed separated flows, transonic boattail flow with and without small angle of attack, transonic base pressures, Mach reflection of shocks, and numerical solution of potential equation through hodograph transformation.

Chow, W. L.↗
NASA NTRS1981

International Conference on Numerical Methods in Fluid Dynamics, 7th, Stanford University, Stanford and Moffett Field, CA, June 23-27, 1980, Proceedings

Topics discussed include polygon transformations in fluid mechanics, computation of three-dimensional horseshoe vortex flow using the Navier-Stokes equations, an improved surface velocity method for transonic finite-volume solutions, transonic flow calculations with higher order finite elements, the numerical calculation of transonic axial turbomachinery flows, and the simultaneous solutions of inviscid flow and boundary layer at transonic speeds. Also considered are analytical solutions for the reflection of unsteady shock waves and relevant numerical tests, reformulation of the method of characteristics for multidimensional flows, direct numerical simulations of turbulent shear flows, the stability and separation of freely interacting boundary layers, computational models of convective motions at fluid interfaces, viscous transonic flow over airfoils, and mixed spectral/finite difference approximations for slightly viscous flows.

Reynolds, W. C.↗
DOE OSTI2020

30 cm Drop Modeling

This report documents modeling and analysis of the 30 cm cask drop scenario for spent nuclear fuel (SNF) casks. The 30 cm drop is an important hypothetical dynamic loading scenario for SNF because 10 CFR 71.71 defines testing requirements for normal conditions of transport (NCT) of SNF packages, and free drop (of 30 cm for the heaviest packages) tests are required. Other NCT conditions specified in 10 CFR 71.71 include vibration normally incident to transport, and DOE’s multimodal transportation test (MMTT) and related analysis concluded that mechanical shock and vibration incident to transportation is negligible for SNF structural integrity because the loads are so small. The SNF cask 30 cm drop has a lot more potential to cause damage to SNF because the kinetic energy and momentum of the system is stopped abruptly during the cask impact event. Test data and model predictions both agree that the SNF assemblies are subjected to significantly higher loads in a cask drop event than normal transportation shock and vibration, but currently there is no reason to expect gross structural damage of the fuel assembly or widespread SNF cladding failure as a result of a 30 cm cask drop. The ultimate goal of this work is to close the knowledge gap of stress profiles on SNF rods and fuel assembly components during NCT, and this report describes the progress toward reaching that goal as well as a plan for finishing this work over the next year or two. Key remaining tasks include refining the fuel assembly model, calculating a selection of specific fuel assembly cases to populate the damage model, and finally running the damage model to predict the full general response range of 30 cm drop loads on SNF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1979

Relativistic blast waves in two dimensions. I - The adiabatic case

Approximate solutions are presented for the dynamical evolution of strong adiabatic relativistic blast waves which result from a point explosion in an ambient gas in which the density varies both with distance from the explosion center and with polar angle in axisymmetry. Solutions are analytical or quasi-analytical for the extreme relativistic case and numerical for the arbitrarily relativistic case. Some general properties of nonplanar relativistic shocks are also discussed, including the incoherence of spherical ultrarelativistic blast-wave fronts on angular scales greater than the reciprocal of the shock Lorentz factor, as well as the conditions for producing blast-wave acceleration.

Shapiro, P. R.↗
NASA NTRS1990

Aerobrake plasmadynamics - Macroscopic effects

The flow around an aerobraking spacecraft (such as the Aeroassist Flight Experiment reentry vehicle) will contain a region of partially ionized gas, that is, a plasma. It is shown here by numerical simulation that macroscopic plasmadynamic effects (which are not included in standard aerothermodynamic simulations) will have a noticeable effect on the reentry flow field. In particular, there are thermoelectric phenomena which can have a major influence in flow dynamics at the front of an ionizing bowshock. These thermoelectric phenomena arise because of the presence of large density and temperature gradients at the front of a reentry bowshock, and they include strong local magnetic fields, electric currents, and ohmic heating. One important result is the dramatic increase in temperature (over the case where plasma effects are neglected) at a reentry shock front; the implication is that macroscopic plasmadynamic effects can no longer be neglected in simulations of hypersonic reentry flow fields.

Shebalin, John V.↗