SEARCH · Search NASA

Results for “Shock fitting”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗

DOE OSTI2021

PBX9012 Tracker Gauge Data

Five shock-to-detonation transition (SDT) gas gun experiments for PBX 9012 were performed by Burns and Chiquete [2020]. The data from the experiments is available on the small-scale database. Here we focus on the tracker gauge data for the lead shock trajectory, x(t). In order to determine the shock speed U s (t) and the detonation transition point, a least squared fit to the shock trajectory is used based on a pair of ODEs with a fitting form for the shock acceleration developed by Hill and Gustavsen [2002].

42 ENGINEERING↗

NASA NTRS1976

A shock surface geometry - The February 15-16, 1967, event

An estimated shape is presented for the surface of the flare-associated interplanetary shock of February 15-16, 1967, as seen in the ecliptic-plane cross section. The estimate is based on observations by Explorer 33 and Pioneers 6 and 7. The estimated shock normal at the Explorer 33 position is obtained by a least-squares shock parameter-fitting procedure for that satellite's data; the shock normal at the Pioneer 7 position is found by using the magnetic coplanarity theorem and magnetic-field data. The average shock speed from the sun to each spacecraft is determined along with the local speed at Explorer 33 and the relations between these speeds and the position of the initiating solar flare. The Explorer 33 shock normal is found to be severely inclined and not typical of interplanetary shocks. It is shown that the curvature of the shock surface in the ecliptic plane near the earth-Pioneer 7 region is consistent with a radius of not more than 0.4 AU.

Lepping, R. P.↗

DOE OSTI2022

The structure of cluster merger shocks: turbulent width and the electron heating time-scale

We present a new 2MsChandra observation of the cluster merger Abell 2146, which hosts two huge M ~ 2 shock fronts each ~500 kpc across. For the first time, we resolve and measure the width of cluster merger shocks. The best-fitting width for the bow shock is 17 ± 1kpc and for the upstream shock is 10.7 ± 0.3 kpc. A narrow collisionless shock will appear broader in projection if its smooth shape is warped by local gas motions. We show that both shock widths are consistent with collisionless shocks blurred by local gas motions of 290 ± 30 km s –1 . The upstream shock forms later on in the merger than the bow shock and is therefore expected to be significantly narrower. From the electron temperature profile behind the bow shock, we measure the time-scale for the electrons and ions to come back into thermal equilibrium. We rule out rapid thermal equilibration of the electrons with the shock-heated ions at the 6σ level. The observed temperature profile instead favours collisional equilibration. For these cluster merger shocks, which have low sonic Mach numbers and propagate through a high β plasma, we find no evidence for electron heating over that produced by adiabatic compression. Our findings are expected to be valid for collisionless shocks with similar parameters in other environments and support the existing picture from the solar wind and supernova remnants. Furthermore, the upstream shock is consistent with this result but has a more complex structure, including a ~2ke V increase in temperature ~50 kpc ahead of the shock.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2017

Feature Detection and Curve Fitting Using Fast Walsh Transforms for Shock Tracking: Applications

Walsh functions form an orthonormal basis set consisting of square waves. Square waves make the system well suited for detecting and representing functions with discontinuities. Given a uniform distribution of 2p cells on a one-dimensional element, it has been proven that the inner product of the Walsh Root function for group p with every polynomial of degree < or = (p - 1) across the element is identically zero. It has also been proven that the magnitude and location of a discontinuous jump, as represented by a Heaviside function, are explicitly identified by its Fast Walsh Transform (FWT) coefficients. These two proofs enable an algorithm that quickly provides a Weighted Least Squares fit to distributions across the element that include a discontinuity. The detection of a discontinuity enables analytic relations to locally describe its evolution and provide increased accuracy. Time accurate examples are provided for advection, Burgers equation, and Riemann problems (diaphragm burst) in closed tubes and de Laval nozzles. New algorithms to detect up to two C0 and/or C1 discontinuities within a single element are developed for application to the Riemann problem, in which a contact discontinuity and shock wave form after the diaphragm bursts.

Gnoffo, Peter A.↗

NASA NTRS1974

Improved curve fits for the thermodynamic properties of equilibrium air suitable for numerical computation using time-dependent or shock-capturing methods, part 1

Simplified curve fits for the thermodynamic properties of equilibrium air were devised for use in either the time-dependent or shock-capturing computational methods. For the time-dependent method, curve fits were developed for p = p(e, rho), a = a(e, rho), and T = T(e, rho). For the shock-capturing method, curve fits were developed for h = h(p, rho) and T = T(p, rho). The ranges of validity for these curves fits were for temperatures up to 25,000 K and densities from 10 to the minus 7th power to 10 to the 3d power amagats. These approximate curve fits are considered particularly useful when employed on advanced computers such as the Burroughs ILLIAC 4 or the CDC STAR.

Tannehill, J. C.↗

NASA NTRS2008

Relationship of Interplanetary Shock Micro and Macro Characteristics: A Wind Study

The non-linear least squared MHD fitting technique of Szabo 11 9941 has been recently further refined to provide realistic confidence regions for interplanetary shock normal directions and speeds. Analyzing Wind observed interplanetary shocks from 1995 to 200 1, macro characteristics such as shock strength, Theta Bn and Mach numbers can be compared to the details of shock micro or kinetic structures. The now commonly available very high time resolution (1 1 or 22 vectors/sec) Wind magnetic field data allows the precise characterization of shock kinetic structures, such as the size of the foot, ramp, overshoot and the duration of damped oscillations on either side of the shock. Detailed comparison of the shock micro and macro characteristics will be given. This enables the elucidation of shock kinetic features, relevant for particle energization processes, for observations where high time resolution data is not available. Moreover, establishing a quantitative relationship between the shock micro and macro structures will improve the confidence level of shock fitting techniques during disturbed solar wind conditions.

Szabo, Adam↗

DOE OSTI2022

Vibrational response and temperature of shock-compressed Pt: In situ extended x-ray absorption fine structure measurements to 325 GPa

Single pulse (~100 ps duration) extended x-ray absorption fine structure (EXAFS) measurements were obtained in laser shocked Pt to examine atomic vibrational properties and to determine temperatures for shock pressures from 72-325 GPa. Pt x-ray absorption spectra were recorded while a planar shock wave was propagating through the Pt sample. The absorption spectrum corresponding to the shocked Pt was obtained by subtracting the ambient Pt contribution from the measured absorption spectrum. Fits to the shocked state EXAFS data provided Pt lattice parameters consistent with the known Pt Hugoniot and the mean squared relative displacements (MSRDs) for nearest neighbor Pt atoms. Pt temperatures in the shocked state, estimated from the MSRDs using the correlated Debye model and correlated Einstein models, were consistent with both the Hugoniot temperatures calculated by integration along the Hugoniot and with the Hugoniot temperatures from published first-principles calculations. However, the agreement between the Hugoniot temperatures from EXAFS measurements and the calculated Hugoniot temperatures is not as good above 200 GPa, likely due to anharmonic effects. Furthermore, the present results demonstrate that single-pulse synchrotron EXAFS measurements in laser-shocked solids are useful for quantitative temperature determination and for examination of vibrational properties, including anharmonicity, at extreme pressure-temperature conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1991

Shock wave equation of state of muscovite

Shock wave data were obtained between 20 and 140 GPa for natural muscovite obtained from Methuen Township (Ontario), in order to provide a shock-wave equation of state for this crustal hydrous mineral. The shock equation of state data could be fit by a linear shock velocity (Us) versus particle velocity (Up) relation Us = 4.62 + 1.27 Up (km/s). Third-order Birch-Murnaghan equation of state parameters were found to be K(OS) = 52 +/-4 GPa and K-prime(OS) = 3.2 +/-0.3 GPa. These parameters are comparable to those of other hydrous minerals such as brucite, serpentine, and tremolite.

Sekine, Toshimori↗

NASA NTRS1980

Molecular hydrogen emission in NGC 7027

The spatial distribution of the emission in the v = 1-0 S(1) line of H2 in the planetary nebula NGC 7027 is presented. The excited H2 molecules do not fill the same volume as the ionized gas and most likely reside near the outer edge of the nebula. The spatial resolution of the data, 5 arcsec, is insufficient to define the location of the molecules precisely. An upper limit to the strength of the v = 2-1 S(1) line is given; it is low enough to preclude simple ultraviolet fluorescence as the source of excitation. A simple shock model can fit the line ratio data, however, and there is enough energy in the expanding nebula to sustain such a shock. A rough estimate of 1-4 solar masses for the mass in the molecular cloud surrounding NGC 7027 is derived.

Beckwith, S.↗

NASA NTRS1974

Curve fits of predicted inviscid stagnation-point radiative heating rates, cooling factors, and shock standoff distances for hyperbolic earth entry

Curve-fit formulas are presented for the stagnation-point radiative heating rate, cooling factor, and shock standoff distance for inviscid flow over blunt bodies at conditions corresponding to high-speed earth entry. The data which were curve fitted were calculated by using a technique which utilizes a one-strip integral method and a detailed nongray radiation model to generate a radiatively coupled flow-field solution for air in chemical and local thermodynamic equilibrium. The range of free-stream parameters considered were altitudes from about 55 to 70 km and velocities from about 11 to 16 km.sec. Spherical bodies with nose radii from 30 to 450 cm and elliptical bodies with major-to-minor axis ratios of 2, 4, and 6 were treated. Powerlaw formulas are proposed and a least-squares logarithmic fit is used to evaluate the constants. It is shown that the data can be described in this manner with an average deviation of about 3 percent (or less) and a maximum deviation of about 10 percent (or less). The curve-fit formulas provide an effective and economic means for making preliminary design studies for situations involving high-speed earth entry.

Suttles, J. T.↗

NASA NTRS1981

Acceleration theory for 5-40 keV ions at interplanetary shocks

A theory is presented to explain the acceleration of suprathermal ions observed near propagating shocks in the solar wind. The hard power-law spectra of ions of energies 5 to 40 keV are shown to be accounted for by diffusive scattering across a plane shock with a small amount of adiabatic deceleration losses on both sides of the shock. The theory fits the observations of the event of November 25, 1977, which followed a large solar flare, and indicates a mean free path perpendicular to the shock of less than 0.0003 AU behind the shock and 0.01 AU in front of it. The theory also predicts a steepening ion energy spectrum at higher energies.

Forman, M. A.↗

DOE OSTI2021

Re-calibration of PBX9501 SURF model

PBX 9501 is a plastic bonded explosive composed of 95 wt % HMX and a binder; see [Gibbs and Popolato, 1980, pp. 109–119]. SURF is a reactive burn model for shock initiation and propagation of detonation waves. It has previously been calibrated for PBX 9501. Here the SURF model is recalibrated for PBX 9501; specifically, lot 730-010 at ρ = 1.837 g/cc 3 . The new calibration uses the Davis reactants and products EOS calibrated for the AWSD model [Aslam et al., 2020]. The burn rate in the shock initiation regime is fit to the Pop plot from 5 embedded gauge shock-to-detonation transition (SDT) experiments from [Gustavsen et al., 1999, see fig 12 and table 5]. In the propagation regime, the burn rate is fit to curvature effect data (detonation speed as function of front curvature); see [Aslam, 2007]. Also the burn parameters are adjusted to fit the gap-stick experiment [Hill et al., 2018]. Simulating the detonation wave speed in this experiment requires a model that is accurate for initiation with complex shock loading; in particular, a pressure decreasing gradient behind a curved lead shock. This is more difficult than calibrating to the standard SDT experiments which are 1-D and driven by a sustained shock. Simulations of the gap-stick experiment for PBX 9501 with the SURF model will be discussed in a subsequent report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

NASA NTRS1986

High speed viscous flow calculations about complex configurations

Applications of the NASA Ames Parabolized Navier-Stokes (PNS) code to a variety of complex generic configurations is presented. The algorithm, boundary conditions, initial conditions, and grid generators are discussed as applied to these configurations. The PNS code was used as the mainline procedure to numerically simulate the viscous supersonic flow over these generic configurations. The turbulence model that was used in this study is the Baldwin-Lomax model. The boundary conditions are the usual viscous no slip at the wall, and a characteristic procedure is used to fit the bow shock wave which is the outermost boundary. An elliptic grid generator is employed to discretize the flow domain. In addition, an equilibrium air capability has been incorporated into the code. It uses the curve fits of Tannehill, et al. The flow regimes vary from a Mach number of 2 up to 25. Both laminar and turbulent flow are considered. Varying angles of attack have also been computed. Configurations vary from simple cone-type bodies to lifting winged bodies, such as the space shuttle or the generic supersonic cruise fighter.

Chaussee, D. S.↗

DOE OSTI2024

Yttrium solid phase equation of state with uncertainty quantification

We discuss the development of an equation of state (EOS) for solid yttrium up to roughly 80 GPa. The EOS makes use of different experimental datasets, including measurements of the isobaric density, isobaric specific heat, room temperature isotherm, and principal shock Hugoniot. The fitting procedure is carried out using Markov Chain Monte Carlo (MCMC), where we fit model parameters for both the cold curve and ion thermal models simultaneously. The results show close agreement with experimental data and provide meaningful uncertainty estimates on the model parameters. This work serves as a first step towards a multiphase EOS for yttrium, which will include higher pressure solid phases (> 80 GPa), as well as modeling of the liquid phase.

36 MATERIALS SCIENCE↗

NASA NTRS1974

Recent developments on time-dependent calculation of nonequilibrium flows

It is shown that, in a time-dependent nonequilibrium flow analysis, erroneous and unstable solutions can result from the use of temperature equations and flow variables other than pressure in fitting the bow shock. A formulation is presented that has been critically evaluated by comparison with a previous one and found to yield steady and more accurate solutions in conserving total enthalpy.

Li, C. P.↗

NASA NTRS1976

A mathematical model of an active control landing gear for load control during impact and roll-out

A mathematical model of an active control landing gear (ACOLAG) was developed and programmed for operation on a digital computer. The mathematical model includes theoretical subsonic aerodynamics; first-mode wing bending and torsional characteristics; oleo-pneumatic shock strut with fit and binding friction; closed-loop, series-hydraulic control; empirical tire force-deflection characteristics; antiskid braking; and sinusoidal or random runway roughness. The mathematical model was used to compute the loads and motions for a simulated vertical drop test and a simulated landing impact of a conventional (passive) main landing gear designed for a 2268-kg (5000-lbm) class airplane. Computations were also made for a simply modified version of the passive gear including a series-hydraulic active control system. Comparison of computed results for the passive gear with experimental data shows that the active control landing gear analysis is valid for predicting the loads and motions of an airplane during a symmetrical landing. Computed results for the series-hydraulic active control in conjunction with the simply modified passive gear show that 20- to 30-percent reductions in wing force, relative to those occurring with the modified passive gear, can be obtained during the impact phase of the landing. These reductions in wing force could result in substantial increases in fatigue life of the structure.

Mcgehee, J. R.↗

NASA NTRS1989

A computational study on the interaction between a vortex and a shock wave

A computational study of two-dimensional shock vortex interaction is discussed in this paper. A second order upwind finite volume method is used to solve the Euler equations in conservation form. In this method, the shock wave is captured rather than fitted so that the cases where shock vortex interaction may cause secondary shocks can also be investigated. The effects of vortex strength on the computed flow and acoustic field generated by the interaction are qualitatively evaluated.

Meadows, Kristine R.↗