Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock jumps”

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 199 records · Page 11

Dynamic fields near a crack tip growing in an elastic-perfectly-plastic solid

A full asymptotic solution is presented for the fields in the neighborhood of the tip of a steadily advancing crack in an incompressible elastic-perfectly-plastic solid. There are four findings for mode I crack growth in the plane strain condition. The first is that the entire crack tip in steady crack growth is surrounded by a plastic region and that no elastic unloading is predicted by the complete dynamic asymptotic solution. The second is that, in contrast to the quasi-static solution, the dynamic solution yields strain fields with a logarithmic singularity everywhere near the crack tip. The third is that whereas the stress field varies throughout the entire crack tip neighborhood, it does not exhibit behavior that can be approximated by a constant field followed by an essentially centered-fan field and then by another constant field, especially for small crack growth speeds. The fourth finding is that there are two shock fronts emanating from the crack tip across which certain stress and strain components undergo jump discontinuities. After reviewing the mode III steady-state crack growth, it is concluded that ductile fracture criteria for nonstationary cracks must be based on solutions that include the inertia effects and that for this purpose quasi-static solutions may be inadequate.

Nemat-Nasser, S.↗

Relaxation and approximate factorization methods for the unsteady full potential equation

The unsteady form of the full potential equation is solved in conservation form, using implicit methods based on approximate factorization and relaxation schemes. A local time linearization for density is introduced to enable solution to the equation in terms of phi, the velocity potential. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity, to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi obtained from requirements of density continuity. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. Results are presented for flows over airfoils, cylinders, and spheres. Comparisons are made with available Euler and full potential results.

Shankar, V.↗

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗

Intermittent Flow Regimes in a Transonic Fan Airfoil Cascade

A study was conducted in the NASA Glenn Research Center linear cascade on the intermittent flow on the suction surface of an airfoil section from the tip region of a modern low aspect ratio fan blade. Experimental results revealed that, at a large incidence angle, a range of transonic inlet Mach numbers exist where the leading-edge shock-wave pattern was unstable. Flush mounted high frequency response pressure transducers indicated large local jumps in the pressure in the leading edge area, which generates large intermittent loading on the blade leading edge. These measurements suggest that for an inlet Mach number between 0.9 and 1.0 the flow is bi-stable, randomly switching between subsonic and supersonic flows. Hence, it appears that the change in overall flow conditions in the transonic region is based on the frequency of switching between two stable flow states rather than on the continuous increase of the flow velocity. To date, this flow behavior has only been observed in a linear transonic cascade. Further research is necessary to confirm this phenomenon occurs in actual transonic fans and is not the byproduct of an endwall restricted linear cascade.

Lepicovsky, J.↗

Computational procedures for mixed equations with shock waves

This paper discusses the procedures we have developed to treat a canonical problem involving a mixed nonlinear equation with boundary data that imply a discontinuous solution. This equation arises in various physical contexts and is basic to the description of the nonlinear acoustic behavior of a shock wave near a caustic. The numerical scheme developed is of second order, treats discontinuities as such by applying the appropriate jump conditions across them, and eliminates the numerical dissipation and dispersion associated with large gradients. Our results are compared with the results of a first-order scheme and with those of a second-order scheme we have developed. The algorithm used here can easily be generalized to more complicated problems, including transonic flows with imbedded shocks.

Yu, N. J.↗

Interaction of interplanetary shocks with nonuniform ambient solar wind

Three interplanetary shock wave events are selected from the plasma and magnetic field data of Helios 1 and 2, IMP-8, and Voyagers 1 and 2 for study of the interactions of a weak interplantary shock with a nonuniform ambient solar wind. These events occurred during the periods 22-26 November 1977, 1-7 January 1978, and 2-5 April 1979, respectively. It is found that the shock surfaces of these events are highly distorted. In addition, a portion of the shock surface may be degenerated into a disturbance which does not satisfy the Rankine-Hugoniot jump conditions.

Chao, J. K.↗

Experimental and Theoretical Investigations of Cavitation in Water

The cavitation in nozzles on airfoils of various shape and on a sphere are experimentally investigated. The limits of cavitation and the extension of the zone of the bubbles in different stages of cavitation are photographically established. The pressure in the bubble area is constant and very low, jumping to high values at the end of the area. The analogy with the gas compression shock is adduced and discussed. The collapse of the bubbles under compression shock produces very high pressures internally, which must be contributory factors to corrosion. The pressure required for purely mechanical corrosion is also discussed.

Ackeret, J.↗

Hybrid interferometric diagnostic for the refractive index measurement across nearly discrete shock waves

The measurement of high temperature gas properties is key for characterizing high speed flows. Nearly discrete changes in density across shock waves, in particular, are difficult to resolve through traditional fringe-counting interferometric methods. Existing techniques for estimating large fringe jumps are either resolution limited or require specialized window configurations. In this Letter, we describe a unique hybrid interferometric technique that combines narrowband fringes for high resolution and broadband fringes as an absolute reference to measure changes in refractive index with a resolution of up to 7 × 10 −8 across nearly discrete index changes of up to 1.5 × 10 −4 . By capturing fringes with an ultrahigh-speed camera, the refractive index changes across discrete shock fronts can be estimated inside a shock tube with high accuracy and time resolution. First, a novel hybrid calibration method for tracking finite fringes is discussed. Next, this technique is used to measure the post-initial-shock refractive indices for Mach 2.7 to 4.2 flows (pressures from 90.4 to 228.4 kPa). Results are then compared with theoretical values showing agreement within 2%.

Wang, Gwendolyn T. (ORCID:000000016271893X)↗

Interface structure of co-rotating interaction regions

Plasma and particle observations on Ulysses during its passes through the southern and northern heliosphere have revealed that, inside the streamer belt, the large-scale structure of the quiet global heliosphere is dominated by corotating interaction regions (CIRs). Therefore, considerable attention is now being given to the internal plasma structure of CIRs, and in particular, to the manifestations of the stream interfaces that should mark their origins as interactions between low speed solar wind (in the low-latitude streamer belt) and high speed solar wind (from the equatorial extensions of the high latitude polar coronal holes). The SWICS and HI-SCALE experiments on Ulysses combine plasma and energetic particle measurements that are of considerable utility for such studies because, between them, they cover the proton energy range from 10 eV to 5 MeV. These measurements are used, together with magnetic field data, to study the remarkable series of CIRs that occurred during the period beginning July 1992 and the end of 1993 as Ulysses rose from the ecliptic to a southern heliographic latitude of 48 deg. The structure of the regions between the forward and reverse shocks were previously analyzed in terms of the proton specific entropy argument log that should exhibit a discontinuous jump at the stream interface. It was claimed that the stream interface, defined with respect to specific entropy, is also associated with a discontinuity in energetic proton intensities. The energetic particle data (greater than 60 keV) and how they were ordered with respect to interfaces and with respect to the magnetic field were examined.

Ogilvie, K. W.↗

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗

Chandra and XMM-Newton Observations of A2256: Cold Fronts, Merger Shocks, and Constraint on the IC Emission

We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic (RR). The temperature and metallicity maps show clear evidence of a merger between the western subcluster (SC) and the primary cluster (PC). We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1 is associated with the infalling SC; CF2 is located in the east of the PC; and CF3 is to the west of the PC core. The other two edges at cluster outskirts are shock fronts(SFs): SF1 near the RR in the NW has Mach numbers derived from the temperature and the density jumps, respectively, of MT=1.62±0.12 and Mρ=1.23±0.06; SF2 in the SE has MT=1.54±0.05 and Mρ=1.16±0.13. In the region of the RR, there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse-Compton emission, and therefore set a lower limit on the magnetic field (∼450 kpc from PC center) of B>1.0μG for a single power-law electron spectrum or B>0.4μG for a broken power-law electron spectrum. We propose a merger scenario including a PC, an SC, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects.

Chong Ge↗

Chandra and XMM-Newton Observations of A2256: Cold Fronts, Merger Shocks, and Constraint on the IC Emission

We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic (RR).The temperature and metallicity maps show clear evidence of a merger between the western subcluster (SC) and the primary cluster (PC). We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1is associated with the infalling SC; CF2 is located in the east of the PC; and CF3 is to the west of the PC core. The other two edges at cluster outskirts are shock fronts(SFs): SF1 near the RR in the NW has Mach numbers derived from the temperature and the density jumps, respectively, of MT=1.62±0.12andMρ=1.23±0.06; SF2 in the SE has MT=1.54±0.05andMρ=1.16±0.13. In the region of the RR, there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse-Compton emission, and therefore set a lower limit on the magnetic field (∼450 kpc from PC center) of B>1.0μG for a single power-law electron spectrum or B>0.4μG for a broken power-law electron spectrum. We propose a merger scenario including a PC, an SC, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects.

Chong Ge↗

Sound velocities in shock-compressed soda lime glass: Melting and liquid-state response

Longitudinal sound velocities (and elastic moduli) were determined in soda lime glass (SLG) at peak shock stresses ranging between 40 and 90 GPa. Laser interferometry was used to obtain particle velocity histories and sound velocities by impacting SLG samples on lithium fluoride (LiF) optical windows. In all experiments, the SLG response consists of a sharp jump to a constant state followed by a release wave. The measured longitudinal sound velocities and moduli showed a marked decrease between 52 and 58 GPa, providing experimental evidence for the transformation from an amorphous solid to a liquid in shock-compressed SLG. Furthermore, the stress threshold (~55 GPa) for melting in SLG is considerably lower than the threshold reported in shock-compressed fused silica (~72 GPa), showing the effect of network-modifying cations on the onset of melting. The relative values of sound velocities, shock velocities, and the Hugoniot slopes—between 58 and 90 GPa—are fully consistent with the thermodynamic response of a shock-compressed liquid. Using the experimental results, the Grüneisen parameter (Γ) values were determined for liquid SLG to 90 GPa and then used to provide the Mie-Grüneisen equation of state for liquid SLG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The axisymmetric jet counterflow problem

A numerical method for computing the flow field of an underexpanded axisymmetric supersonic jet flowing counter to a uniform supersonic stream is described. The strong shock waves and shear layers which occur in the flow are treated as normal and tangential discontinuities bounding regions of inviscid flow. Within the regions, the unsteady gasdynamic equations are solved in integral form on a moving computational mesh. As the numerical flow field evolves, the mesh is deformed according to the jump relations applied at the discontinuities, thus determining their positions as part of the solution. Also presented are computed results obtained for a wide range of jet total pressure. The standoff distances of the bow shock and jet shock are in good agreement with those measured experimentally, thus demonstrating the accuracy of the numerical method.

Schiff, L. B.↗

Cylindrical sound wave generated by shock-vortex interaction

The passage of a columnar vortex broadside through a shock is investigated. This has been suggested as a crude, but deterministic, model of the generation of 'shock noise' by the turbulence in supersonic jets. The vortex is decomposed by Fourier transform into plane sinusoidal shear waves disposed with radial symmetry. The plane sound waves produced by each shear wave/shock interaction are recombined in the Fourier integral. The waves possess an envelope that is essentially a growing cylindrical sound wave centered at the transmitted vortex. The pressure jump across the nominal radius R = ct attenuates with time as 1/(square root of R) and varies around the arc in an antisymmetric fashion resembling a quadrupole field. Very good agreement, except near the shock, is found with the antisymmetric component of reported interferometric measurements in a shock tube. Beyond the front r approximately equals R is a precursor of opposite sign, that decays like 1/R, generated by the 1/r potential flow around the vortex core. The present work is essentially an extension and update of an early approximate study at M = 1.25. It covers the range (R/core radius) = 10, 100, 1000, and 10,000 for M = 1.25 and (in part) for M = 1.29 and, for fixed (R/core radius) = 1000, the range M = 1.01 to infinity.

Ribner, H. S.↗

Computation of unsteady transonic flows through rotating and stationary cascades. 1: Method of analysis

A numerical method of solution of the inviscid, compressible, two-dimensional unsteady flow on a blade-to-blade stream surface through a stage (rotor and stator) or a single blade row of an axial flow compressor or fan is described. A cyclic procedure has been developed for representation of adjacent blade-to-blade passages which asymptotically achieves the correct phase between all passages of a stage. A shock-capturing finite difference method is employed in the interior of the passage, and a method of characteristics technique is used at the boundaries. The blade slipstreams form two of the passage boundaries and are treated as moving contact surfaces capable of supporting jumps in entropy and tangential velocity. The Kutta condition is imposed by requiring the slipstreams to originate at the trailing edges, which are assumed to be sharp. Results are presented for several transonic fan rotors and compared with available experimental data, consisting of holographic observations of shock structure and pressure contour maps. A subcritical stator solution is also compared with results from a relaxation method. Finally, a periodic solution for a stage consisting of 44 rotor blades and 46 stator blades is discussed.

Erdos, J. I.↗

Numerical solution of periodic transonic flow through a fan stage

A numerical method of solution of the inviscid, compressible, two-dimensional unsteady flow on a blade-to-blade stream surface through a stage (rotor and stator), or a single blade row, of an axial flow compressor or fan is described. A cyclic procedure has been developed for representation of adjacent blade-to-blade passages, which asymptotically achieves the correct phase between all passages of a stage. A shock-capturing finite-difference method is employed in the interior of the passage, and a method-of-characteristics technique is used at the boundaries. The blade slipstreams form two of the passage boundaries, and are treated as moving contact surfaces capable of supporting jumps in entropy and tangential velocity. The Kutta condition is imposed by requiring the slipstreams to originate at the trailing edges, which are assumed to be sharp. Results are presented for several transonic fan rotors, and compared with available experimental data consisting of holographic observations of shock structure and pressure contour maps. A subcritical stator solution is also compared with results from a relaxation method. Finally, a periodic solution for a stage consisting of 44 rotor blades and 46 stator blades is discussed.

Erdos, J. I.↗

Viscous reacting flows with wall slip and catalysis applied to spheres in arc jets and flight

The influence of wall slip and catalytic atom-recombination on the flow field and wall heat flux are calculated for high altitude flight and arc jet flow conditions. Boundary equations, which include velocity slip, temperature jump, and wall catalytic atom recombination, are coupled to the viscous reacting multicomponent Navier-Stokes equations. These equations are solved using a time-dependent finite difference technique applied to spheres in an arc jet flow (Reynolds number of 550) and a high altitude flight case representative of the Space Shuttle Orbiter (Reynolds number of 450). The results indicate that catalysis strongly influences the temperature jump, but not the velocity slip. Slip increases the atom fraction and temperature at both the wall and the flow field. Likewise, the shock stand-off distance, the wall heat flux, and friction coefficient are increased over the nonslip cases. The reacting gas calculations confirm the chemically frozen nature of the shock layer in arc jet flows.

Scott, C. D.↗