Search NASA⌕ Search

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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Comparison of jet plume shape predictions and plume influence on sonic boom signature

An Euler shock-fitting marching code yields good agreement with semiempirically determined plume shapes, although the agreement decreases somewhat with increasing nozzle angle and the attendant increase in the nonisentropic nature of the flow. Some calculations for the low boom configuration with a simple engine indicated that, for flight at altitudes above 60,000 feet, the plume effect is dominant. This negates the advantages of a low boom design. At lower altitudes, plume effects are significant, but of the order that can be incorporated into the low boom design process.

Barger, Raymond L.↗
NASA NTRS1992

Spectral solution of inviscid supersonic flows over wedges and axisymmetric cones

A shock-fitted multidomain spectral collocation method is used to solve both steady and unsteady inviscid supersonic flows over bodies. New aspects of the method include two subdomain interface types and a zonal solution procedure to get efficient convergence to steady-state in supersonic regions. Steady-state examples include flow over a sharp cone, a hyperbolic cone and a hyperbolic wedge. For the cone, the exact flow solution is used to show that the method is spectrally accurate. As an example of an unsteady flow, a calculation of the interaction of a free-stream hot ring with the flow over a sharp cone is presented.

Kopriva, David A.↗
NASA NTRS2000

Astrophysical Constraints on Chondrule Formation Models

Three petrologic constraints, heating and cooling timescales and the frequency of compound chondrules, are applied to the shock wave and lightning models of chondrule formation. Shocks fit the constraints well; lightning requires high gas densities.

Desch, S. J.↗
NASA NTRS1971

Influence of thermal anisotropy on best-fit estimates of shock normals

The influence of thermal anisotropy on the estimates of interplanetary shock parameters and the associated normals is discussed. A practical theorem is presented for quantitatively correcting for anisotropic effects by weighting the before and after magnetic fields by the same anisotropy parameter h. The quantity h depends only on the thermal anisotropies before and after the shock and on the angles between the magnetic fields and the shock normal. The theorem can be applied to most slow shocks, but in those cases h usually should be lower, and sometimes markedly lower, than unity. For the extreme values of h, little change results in the shock parameters or in the shock normal.

Lepping, R. P.↗
NASA NTRS1972

Influence of thermal anisotropy on best-fit estimates of shock normals.

This paper deals with the influence of thermal anisotropy on least-squares estimates of interplanetary shock parameters and the associated normals by using the Rankine-Hugoniot equations. A practical theorem is given for quantitatively correcting for anisotropic effects by weighting the before and after magnetic fields by the same 'anisotropy parameter' h. The quantity h depends only on the thermal anisotropies before and after the shock and on the angles between the magnetic fields and the shock normal. It is shown that, for fast shocks and for a liberal range of realistic conditions, the quantity h lies in the range from 0.90 to 1.22. The theorem can also be applied to most slow shocks, but in those instances h usually should be lower and sometimes markedly lower than unity.

Lepping, R. P.↗
NASA NTRS1983

Transonic small disturbance calculations including entropy corrections

Murman's fully conservative mixed type finite-difference operators are first modified. A special sonic point operator with an iterative damping term is introduced which helps the convergence and does not affect the spatial conservative differences. Reliable calculations with second order supersonic schemes are obtained using two sonic operators, the regular sonic point operator followed by a first order supersonic scheme. Also, shock point operator is shown to be equivalent to fitting a locally normal shock terminating the supersonic region. The potential calculations are then modified to account for the non-isentropic jump conditions using a simple shock fitting procedure based on Prandtl relation. The entropy increase across the shock is calculated in terms of the Mach number upstream of the shock and the effect of the generated vorticity is estimated via Crocco relation. Different examples are calculated and extensions to the full potential equation are discussed.

Hafez, M.↗
NASA NTRS1982

Computation of nonlinear supersonic potential flow over three dimensional surfaces

A fully implicit three dimensional marching technique is presented for the computation of inviscid supersonic flows. The nonconservative full potential equation is solved on spherical cross flow planes using standard transonic SLOR techniques. The three dimensional bow shock is fitted as a boundary with isentropic jump conditions and embedded shocks are captured. Applications of these techniques are presented for three dimensional bodies, wings, and wing-body configurations and the computed solutions are compared to experimental pressure data. Zero lift wave drag calculations are also presented and compared with measurements to evaluate the accuracy of the present techniques.

Siclari, M. J.↗
NASA NTRS1999

Multi-Spacecraft Observations of Interplanetary Shocks

Using multi-spacecraft observations primarily from ACE and WIND, and from IMP 8 and Geotail when available, the 3-dimensional structure of interplanetary shocks on the hundred Earth radii scale will be discussed. The complete magnetic field, and solar wind ion and electron data sets were used to fit the shocks with a full non-linear least squares fit "Rankine-Hugoniot" technique yielding the local shock surface normals and speeds with associated uncertainties. Multi-spacecraft results reveal that on the distance scale of ACE's L1 halo orbit the shocks deviate significantly from a simple planar geometry. This result has important consequences for the prediction of the exact arrival times of interplanetary shocks at the Earth's magnetosphere, and hence, on the reliability of space weather predictions. It also has implications on the coherence scale of solar wind structures and their evolution from the Sun to Earth.

Szabo, A.↗
NASA NTRS1999

Multi-Spacecraft Observations of Interplanetary Shocks

Using multi-spacecraft observations primarily from ACE and WIND and from IMP 8 and Geotail when available, the 3-dimensional structure of interplanetary shocks on the hundred Earth radii scale will be discussed. The complete magnetic field, and solar wind ion and electron data sets were used to fit the shocks with a full non-linear least squares fitting "Rankine-Hugoniot" technique yielding the local shock surface normals and speeds with associated uncertainties. Multi-spacecraft results reveal that on the distance scale of ACE's L1 halo orbit the shocks deviate from a simple planar geometry. This result has important consequences for the prediction of the exact arrival times of interplanetary shocks at the Earth's magnetosphere, and hence, on the reliability of space weather predictions. It also has implications on the coherence scale of solar wind structures and their evolution from the Sun to Earth.

Smith, C. W.↗
NASA NTRS1982

The observation of a coronal transient directed at earth

The paper reports the observation of a large coronal transient that can only be interpreted as a three-dimensional structure. Its form is one which has not been observed before: a gradually expanding, sun-centered disk of excess brightness, whose projected radius increased from 4 to 8 solar radii during 0832-0958 UT on November 27, 1979. This earth-directed transient was the source of an interplanetary shock wave that reached ISEE 3 at 0649 UT, November 30, and earth at 0738 UT, November 30. Fitting the shock speed at ISEE 3 and the average transit speed from the sun to ISEE 3 to a power law of the form V = (V sub 0)(r exp -n), it is found that V sub 0 = 1980 km/s and n = 0.294, in good agreement with shock wave models. The shock speed predicted by the power law at 10 solar radii is 1000 km/s, in good agreement with the estimated frontal speed of the transient.

Howard, R. A.↗
NASA NTRS1992

Tests of shock chemistry in IC 443G

Eight molecular species, in the hot dense clump IC 443G, believed to be impacted by the shock wave from the SNR IC 443, are investigated. The clump consists of two distinct regions, one relatively cool, and one hotter and denser. Region 1 contains CO, HCO(+), HCN, and CN, whose abundances may be explained either by ion-molecule chemistry, or by a D shock of 60-90 km/s, passing through a clump of about 100,000/cu cm. Region 2 gives rise to SiO, CS, SO, and H2CO, and requires an ND shock of 5-15 km/s passing through a region of about 1,000,000/cu cm. Observed fractional abundances fit ND shock models if L is about 6.6 x 10 exp 15 cm. In general, observed line widths vary inversely with derived excitation density, while centroid velocities of all species are essentially identical.

Turner, B. E.↗
NASA NTRS1985

Numerical studies of 2-dimensional flows

A formulation of the lambda scheme for the analysis of two dimensional inviscid, compressible, unsteady transonic flows is presented. The scheme uses generalized Riemann variables to determine the appropriate two point, one sided finite difference approximation for each derivative in the unsteady Euler equations. These finite differences are applied at the predictor and corrector levels with shock updating at each level. The weaker oblique shocks are captured, but strong near normal shocks are fitted into the flow using the Rankine-Hugoniot relations. This code is demonstrated with a numerical example of a duct flow problem with developing normal and oblique shock waves. The technique is implemented in a code which has been made efficient by streamlining to a minimal number of operations and by eliminating branch statements. The scheme is shown to provide an accurate analysis of the flow, including formation, motions, and interactions of shocks; the results obtained on a relatively coarse mesh are comparable to those obtained by other methods on much finer meshes.

Moretti, G.↗
NASA NTRS1985

The energy spectra of solar flare electrons

A survey of 50 electron energy spectra from .1 to 100 MeV originating from solar flares was made by the combination of data from two spectrometers onboard the International Sun Earth Explorer-3 spacecraft. The observed spectral shapes of flare events can be divided into two classes through the criteria of fit to an acceleration model. This standard two step acceleration model, which fits the spectral shape of the first class of flares, involves an impulsive step that accelerates particles up to 100 keV and a second step that further accelerates these particles up to 100 MeV by a single shock. This fit fails for the second class of flares that can be characterized as having excessively hard spectra above 1 MeV relative to the predictions of the model. Correlations with soft X-ray and meter radio observations imply that the acceleration of the high energy particles in the second class of flares is dominated by the impulsive phase of the flares.

Evenson, P. A.↗
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↗
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.↗
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.↗