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 847 records · Page 47

Synoptic observations of coronal transients and their interplanetary consequences

Conclusions reached to date after 5 yr of pooled, extensive monitoring of the solar coronal mass ejection (CME) at the Hawaiian High Altitude Observatory, Skylab and with the SMM are reported. Additional white light data have been gathered with the OSO-7 and P78-1 spacecraft. CME provides 5 percent of the solar wind mass flux and was the dominant driving force in interplanetary shocks in the last solar maximum phase. Both bubble and cloud shapes have been observed in CME events, which are nearly ubiquitous in proton events. Each of the CME shapes possesses distinctive dynamical characteristics. Finally, steady emissions of soft X-rays have been identified as precursors to CMEs, which display some latitudinal confinement.

Michels, D. J.↗

A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes

A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.

Batina, John T.↗

A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes

A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.

Batina, John T.↗

Internal Designs Application for Inlet and Nozzle Aeroperformance Improvement

The following research results are based on development of an approach previously proposed by the authors for optimum nozzle design to obtain maximum thrust. The design was denoted a Telescope nozzle. A Telescope nozzle contains one or several internal designs of certain location, which are inserted at certain locations into a divergent conical or planar main nozzle near its exit. Such a design provides additional thrust augmentation over 20% by comparison with the optimum single nozzle of equivalent lateral area. What is more, recent experimental acoustic tests have discovered an essential noise reduction due to Telescope nozzles application. In this paper, some additional theoretical results are presented for Telescope nozzles and a similar approach is applied for aeroperformance improvement of a supersonic inlet. In addition, a classic gas dynamics problem of a similar supersonic flow into a plate has been analyzed. In some particular cases, new exact analytical solutions are obtained for a flow into a wedge with an oblique shock wave. Numerical simulations were conducted for supersonic flow into a divergent portion of a 2D or axisymmetric nozzle with several plane or conuical designs as well as into a 2D or axisymmetric supersonic inlet with a forebody. The 1st order Kryko-Godunov marching numerical scheme for inviscid supersonic flows was used. Several cases were tested using the NASA CFL3d code based on full Navier-Stokes equations. Numerical simulation results have confirmed essential benefits of Telescope design applications in propulsion systems.

Gilinsky, M.↗

Computational Study of Axisymmetric Off-Design Nozzle Flows

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles operating off-design at transonic Mach numbers have been completed. These computations span the very difficult transonic flight regime with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined, including the Explicit Algebraic Stress model. Computations have been completed at freestream Mach numbers of 0.9 and 1.2, and nozzle pressure ratios (NPR) of 4 and 6. Calculations completed with variable time-stepping (steady-state) did not converge to a true steady-state solution. Calculations obtained using constant timestepping (timeaccurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was the result of using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show reasonable agreement with experimental data. The SST turbulence model demonstrates the best overall agreement with experimental data.

DalBello, Teryn↗

Internal Designs Application for Inlet and Nozzle Aeroperformance Improvement

The following research results are based on development of an approach previously proposed by the authors for optimum nozzle design to obtain maximum thrust. The design was denoted a Telescope nozzle. A Telescope nozzle contains one or several internal designs of certain location, which are inserted at certain locations into a divergent conical or planar main nozzle near its exit. Such a design provides additional thrust augmentation over 20% by comparison with the optimum single nozzle of equivalent lateral area. What is more, recent experimental acoustic tests have discovered an essential noise reduction due to Telescope nozzles application. In this paper, some additional theoretical results are presented for Telescope nozzles and a similar approach is applied for aeroperformance improvement of a supersonic inlet. In addition, a classic gas dynamics problem of a similar supersonic flow into a plate has been analyzed. In some particular cases, new exact analytical solutions are obtained for a flow into a wedge with an oblique shock wave. Numerical simulations were conducted for supersonic flow into a divergent portion of a 2D or axisymmetric nozzle with several plane or conical designs as well as into a 2D or axisymmetric supersonic inlet with a forebody. The 1st order Kryko-Godunov marching numerical scheme for inviscid supersonic flows was used. Several cases were tested using the NASA CFL3d code based on full Navier-Stokes equations. Numerical simulation results have confirmed essential benefits of Telescope design applications in propulsion systems.

Gilinsky, M.↗

Internal Designs Application for Inlet and Nozzle Aeroperformance Improvement

The following research results are based on development of an approach previously proposed by the authors for optimum nozzle design to obtain maximum thrust. The design was denoted a Telescope nozzle. A Telescope nozzle contains one or several internal designs of certain location, which are inserted at certain locations into a divergent conical or planar main nozzle near its exit. Such a design provides additional thrust augmentation over 20% by comparison with the optimum single nozzle of equivalent lateral area. What is more, recent experimental acoustic tests have discovered an essential noise reduction due to Telescope nozzles application. In this paper, some additional theoretical results are presented for Telescope nozzles and a similar approach is applied for aeroperformance improvement of a supersonic inlet. In addition, a classic gas dynamics problem of a similar supersonic flow into a plate has been analyzed. In some particular cases, new exact analytical solutions are obtained for a flow into a wedge with an oblique shock wave. Numerical simulations were conducted for supersonic flow into a divergent portion of a 2D or axisymmetric nozzle with several plane or conical designs as well as into a 2D or axisymmetric supersonic inlet with a forebody. The 1st order Kryko-Godunov march- ing numerical scheme for inviscid supersonic flows was used. Several cases were tested using the NASA CFL3d code based on full Navier-Stokes equations. Numerical simulation results have confirmed essential benefits of Telescope design applications in propulsion systems.

Gilinsky, M.↗

Time-Frequency Analysis of Rocket Nozzle Wall Pressures During Start-up Transients

Surveys of the fluctuating wall pressure were conducted on a sub-scale, thrust- optimized parabolic nozzle in order to develop a physical intuition for its Fourier-azimuthal mode behavior during fixed and transient start-up conditions. These unsteady signatures are driven by shock wave turbulent boundary layer interactions which depend on the nozzle pressure ratio and nozzle geometry. The focus however, is on the degree of similarity between the spectral footprints of these modes obtained from transient start-ups as opposed to a sequence of fixed nozzle pressure ratio conditions. For the latter, statistically converged spectra are computed using conventional Fourier analyses techniques, whereas the former are investigated by way of time-frequency analysis. The findings suggest that at low nozzle pressure ratios -- where the flow resides in a Free Shock Separation state -- strong spectral similarities occur between fixed and transient conditions. Conversely, at higher nozzle pressure ratios -- where the flow resides in Restricted Shock Separation -- stark differences are observed between the fixed and transient conditions and depends greatly on the ramping rate of the transient period. And so, it appears that an understanding of the dynamics during transient start-up conditions cannot be furnished by a way of fixed flow analysis.

Baars, Woutijn J.↗

Development of a Numerical Model of Hypervelocity Impact into a Pressurized Composite Overwrapped Pressure Vessel

As the outlook for space exploration becomes more ambitious and spacecraft travel deeper into space than ever before, it is increasingly important that propulsion systems perform reliably within the space environment. The increased reliability compels designers to increase design margin at the expense of system mass, which contrasts with the need to limit vehicle mass to maximize payload. Such are the factors that motivate the integration of high specific strength composite materials in the construction of pressure vessels commonly referred to as composite overwrapped pressure vessels (COPV). The COPV consists of a metallic liner for the inner shell of the COPV that is stiff, negates fluid permeation and serves as the anchor for composite laminates or filaments, but the liner itself cannot contain the stresses from the pressurant it contains. The compo-site-fiber reinforced polymer (CFRP) is wound around the liner using a combination of hoop (circumferential) and helical orientations. Careful consideration of wrap orientation allows the composite to evenly bear structural loading and creates the COPV's characteristic high strength to weight ratio. As the CFRP overwrap carries most of the stresses induced by pressurization, damage to the overwrap can affect mission duration, mission success and potentially cause loss-of-vehicle/loss-of-crew. For this reason, it is critical to establish a fundamental understanding of the mechanisms involved in the failure of a stressed composite such as that of the COPV. One of the greatest external threats to the integrity of a spacecraft's COPV is an impact from the meteoroid and orbital debris environments (MMOD). These impacts, even from submillimeter particles, generate extremely high stress states in the CFRP that can damage numerous fibers. As a result of this possibility, initial assumptions in survivability analysis for some human-rated NASA space-craft have assumed that any alteration of the vessel due to impact is considered a catastrophic failure. This assumption is conservative and made due to lack of knowledge on the level of allow-able damage to the composite overwrap that can be sustained and still allow successful completion of the mission. To quantify the allowable damage level to the composite overwrap involves assessing stress redistribution following damage as well as evaluating possible time-dependent mechanisms involved in the COPV response to an impact event. Limited published work in this subject has shown that COPV can withstand at least some level of damage due to high energy impacts. These observations have been confirmed and expanded upon in recent experimental research performed by NASA. This research has demonstrated that there is not only robustness in a COPV to compensate for CFRP damage, but has also identified two significant failure modes for pressurized COPV. The lowest threshold failure mode involves the perforation of the vessel, and the highest threshold failure mode is the catastrophic rupture. While both of these failure modes mean a loss of the COPV, system robustness affords some tolerance to the venting as opposed to the more catastrophic rupture. As a consequence, it is necessary to understand the conditions that result in the transition between these failure modes. The aforementioned experimental research has been performed in both the unpressurized and pressurized condition to identify the damage level that triggered the failure thresh-old. This COPV test program was sponsored by the NASA Engineering and Safety Center (NESC), and tests were performed at NASA White Sands Test Facility (WSTF). Planning and coordination were provided by NASA JSC Hypervelocity Impact Technology (HVIT) group, and the COPVs were provided by the ISS Program. Unpressurized testing has been conducted at the pressure of the vacuum test chamber, while, the pressurized testing has been conducted at 290 +/- 10 bar (4,200  100 psi) using nitrogen as the pressurizing gas, which corresponds to the design pressure for the target COPV. In this research, spherical aluminum projectiles with varying diameter has been chosen as the impactor. For the unpressurized COPV, the dependence of penetration up to the dependence of hole size in the liner has been obtained as a function of impact conditions. For the pressurized research, the dependence of penetration up to rupture has been obtained as a function of im-pact conditions. Two representative post-test photographs of the failed COPV's from a nor-mal impact into the COPV surface are shown in Fig. 1. These images display the dramatic difference between failure modes, venting (Fig. 1a) and rupture (Fig. 1b). For venting, liner perforation, severed composite fibers/tows and ply delamination are commonly observed damage characteristics of this COPV failure mode. In the case of rupture, the COPV typically experienced a separation of its domed regions and severe break-up of the cylindrical region. Fully understanding the transition from venting to rupture experimentally is costly and potentially unachievable for conditions that cannot be generated in the laboratory. These shortcomings have motivated the performance of three-dimensional numerical simulations to expand the existing experimental database. These simulations have been carried out with the nonlinear-structural-dynamics, analysis-tool, CTH. A typical pressure contour plot from an impact simulation of an entire COPV is shown in Fig. 2. To generate the COPV stress state without initiating a shock wave, the pressure in the simulated COPV is ramped up to the final pressure over a millisecond prior to impact of the projectile with nitrogen gas. Figure 2a shows the system in this initial condition. After one millisecond, a projectile is initiated into the simulation and impacts the COPV. Figure 2b shows the system after this impact. In the figure, the onset of venting is represented as the change in pressure (μbar), red to green, at the perforation site. Also seen in the figure is the eroded projectile that had passed into the COPV vessel with the generated shock wave in the pressurant propagating just ahead of the material. In this paper, pertinent experimental details and the development of the material constitutive models necessary for this work along with the efforts to validate their use are dis-cussed. The simulation results are presented and compared with the NASA experimental observations. While work is on-going from this effort, early observations pertinent to the failure threshold are presented.

Garcia, M. A.↗

Budget of Turbulent Kinetic Energy in a Shock Wave Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Quantities present in the exact equation of the turbulent kinetic energy transport were accumulated and used to calculate terms like production, dissipation, molecular diffusion, and turbulent transport. The present results for a turbulent boundary layer were validated by comparison with direct numerical simulation data. It was found that a longer development domain was necessary for the boundary layer to reach an equilibrium state and a finer mesh resolution would improve the predictions. In spite of these findings, trends of the present budget match closely with that of the direct numerical simulation. Budgets for the SBLI region are presented at key axial stations. These budgets showed interesting dynamics as the incoming boundary layer transforms and the terms of the turbulent kinetic energy budget change behavior within the interaction region.

shock wave/boundary-layer interaction↗

Dynamics of the Uranian Rings

Some of the problems of the shepherding satellite model of Goldreich ant tremaine are discussed. The following topics are studied: (1) optical depths of the all the observed narrow rings; (2) satellite and ring separation timescales; (3) ring edge sharpness; (4) shock formation in narrow rings; (5) the existence of small satellites near the Uranian rings; and (6) the apse and node alignments of the eccentric and inclined rings.

Dermott, S. F.↗

NASA LaRC Contribution to the High Angle Working Group of the Third Aeroelastic Prediction Workshop: BSCW Shock Buffet

FUN3D Core Capabilities - Established as a research code in late 1980s; now supports numerous internal and external efforts across the speed range - Solves 2D/3D steady and unsteady Euler and RANS equations on node-based mixed element grids for compressible and incompressible flows - General dynamic mesh capability: any combination of rigid / overset / morphing grids, including 6-DOF effects - Aeroelastic modeling using mode shapes, full FEM, etc. - Constrained / multipoint adjoint-based design and mesh adaptation - Distributed development team using agile/extreme software practices including 24/7 regression, performance testing - Capabilities fully integrated, online documentation, training videos, tutorials

Pawel Chwalowski↗

Laser-Induced Thermal Acoustics Theory and Expected Experimental Errors when Applied to a Scramjet Isolator Model

A scramjet isolator model test apparatus is being assembled in the Isolator Dynamics Research Lab (IDRL) at the NASA Langley Research Center in Hampton, Virginia. The test apparatus is designed to support multiple measurement techniques for investigating the flow field in a scramjet isolator model. The test section is 1-inch high by 2-inch wide by 24-inch long and simulates a scramjet isolator with an aspect ratio of two. Unheated, dry air at a constant stagnation pressure and temperature is delivered to the isolator test section through a Mach 2.5 planar nozzle. The isolator test section is mechanically back-pressured to contain the resulting shock train within the 24-inch isolator length and supports temperature, static pressure, and high frequency pressure measurements at the wall. Additionally, nonintrusive methods including laser-induced thermal acoustics (LITA), spontaneous Raman scattering, particle image velocimetry, and schlieren imaging are being incorporated to measure off-wall fluid dynamic, thermodynamic, and transport properties of the flow field. Interchangeable glass and metallic sidewalls and optical access appendages permit making multiple measurements simultaneously. The measurements will be used to calibrate computational fluid dynamics turbulence models and characterize the back-pressured flow of a scramjet isolator. This paper describes the test apparatus, including the optical access appendages; the physics of the LITA method; and estimates of LITA measurement uncertainty for measurements of the speed of sound and temperature.

Middleton, Troy F.↗

Approximate and exact numerical integration of the gas dynamic equations

A highly accurate approximation and a rapidly convergent numerical procedure are developed for two dimensional steady supersonic flow over an airfoil. Examples are given for a symmetric airfoil over a range of Mach numbers. Several interesting features are found in the calculation of the tail shock and the flow behind the airfoil.

Lewis, T. S.↗

Comet-solar wind interaction - Dynamical length scales and models

ICE magnetometer measurements at Comet Giacobini-Zinner and model simulations of comet-solar wind interactions are analyzed. The magnetometer data reveal the existence of intense hydromagnetic turbulence, a draping of the magnetic field lines to form a magnetotail, a weak shock, and a magnetic barrier region in the magnetosphere. The global models of the comet-solar wind interaction are described. The observed data and models are compared and good correlation is displayed.

Mendis, D. A.↗

Observations of Deep Ionospheric F-Region Density Depletions with FPMU Instrumentation and their Relationship with the Global Dynamics of the June 22-23, 2015 Geomagnetic Storm

The magnetic storm that commenced on June 22, 2015 was one of the largest storms in the current solar cycle, resulting from an active region on the Sun that produced numerous coronal mass ejections (CMEs) and associated interplanetary shock waves. On June 22 at 18:36 UT the magnetosphere was impacted by the leading-edge shock wave and a sheath carrying a large and highly variable interplanetary magnetic field (IMF) Bz with values ranging from +25 to -40 nT. During the subsequent interval from 0000 to 0800 UT, there was a second intensification of the geomagnetic storm resulting from the impact of the CME. We present dramatic responses of simultaneous particle measurements from the high-altitude Magnetospheric Multiscale Mission (MMS) at high altitudes in the magnetosphere (approx. 9-12 Re) and from the low-altitude (F-region) Floating Potential Measurement Unit (FPMU) on board the International Space Station (ISS). We analyze potential causes of these dramatic particle flux dropouts by putting them in the context of storm-time electrodynamics, and support our results with numerical simulations of the global magnetosphere and ionosphere. During the sheath phase of the storm, the MMS spacecraft in the near-earth equatorial plane observed a rapid reconfiguration of the magnetic field near 1923 UT. Initially in the warm plasmasheet, particle flux dropouts were observed as they tracked the plasma-sheet to lobe transitions with the stretching and thinning of the plasmasheet. Anti-sunward flowing O+ ions of ionospheric origin were also measured during this period, confirming that the MMS spacecraft temporarily was in a lobe.

Coffey, Victoria↗

Observations of Deep Ionospheric F-Region Density Depletions with FPMU Instrumentation and Their Relationship with the Global Dynamics of the June 22-23, 2015 Geomagnetic Storm

The magnetic storm that commenced on June 22, 2015 was one of the largest storms in the current solar cycle, resulting from an active region on the Sun that produced numerous coronal mass ejections (CMEs) and associated interplanetary shock waves. On June 22 at 18:36 UT the magnetosphere was impacted by the leading-edge shock wave and a sheath carrying a large and highly variable interplanetary magnetic field (IMF) Bz with values ranging from +25 to -40 nT. During the subsequent interval from 0000 to 0800 UT, there was a second intensification of the geomagnetic storm resulting from the impact of the CME. We present dramatic responses of simultaneous particle measurements from the high-altitude Magnetospheric Multiscale Mission (MMS) at high altitudes in the magnetosphere (approx. 9-12 Re) and from the low-altitude (F-region) Floating Potential Measurement Unit (FPMU) on board the International Space Station (ISS). We analyze potential causes of these dramatic particle flux dropouts by putting them in the context of storm-time electrodynamics, and support our results with numerical simulations of the global magnetosphere and ionosphere. During the sheath phase of the storm, the MMS spacecraft in the near-earth equatorial plane observed a rapid reconfiguration of the magnetic field near 1923 UT. Initially in the warm plasma sheet, particle flux dropouts were observed as they tracked the plasma-sheet to lobe transitions with the stretching and thinning of the plasma sheet. Anti-sunward flowing O+ ions of ionospheric origin were also measured during this period, confirming that the MMS spacecraft temporarily was in a lobe.

Coffey, Victoria↗

On the practical use of high-order methods for hyperbolic systems

The paper tests a number of high order methods on a variety of dynamic problems in one, two, and three space dimensions. The problems covered include wave propagation phenomena as well as an asymptotic approach to a steady state. Consideration is given to both smooth and shocked flows. It is shown that the methods compared require only minor modifications of many existing second-order schemes. Further, the results show that significant gains can be expected from the use of fourth-order methods. Finally, spectral methods are also considered for some of the problems presented.

Turkel, E.↗