SEARCH · Engineering Papers
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 the present time the efforts on this project are organized into two areas: the development of an approximate stagnation point solution and approximate flowfield studies which can be used to develop and investigate shock jump, electron temperature, radiation, vibration-dissociation coupling, and chemistry models. Progress in each area is discussed.
The breakdown of the conservative potential approximation occurs for all the airfoils tested. It develops as soon as the shock waves appear in the flow field. Since shock waves are not properly represented by the potential approximation, it is conjectured that the breakdown is due to the isentropic shock jump condition of the potential approximation.
NASA's ASPIRE (Advanced Supersonic Parachute Inflation Research and Experiments) project is investigating the supersonic deployment and inflation of full-scale parachutes. To achieve Mars-relevant conditions, the parachutes are deployed at high altitudes over Earth on a sounding rocket platform. During the flight test, Disk-Gap-Band parachutes of 21.5 meter diameter are deployed behind a slender payload 1/6th the diameter of the blunt Mars2020 capsule. Due to the differences in leading body geometry between the test flight and a parachute deployment at Mars, high fidelity numerical simulations of slender and blunt bodywakes, and of rigid parachutes behind them, were used to understand differences and similarities in the flow and the effect on parachute drag. The slender body wake is thinner, closes earlier, and presents a smaller wake deficit. Thus, a parachute deployed in the wake of a slender body is more likely to see a higher dynamic pressure than a parachute deployed behind a blunt body. In the presence of a parachute, the interaction of the unsteady wake with the parachute bow shock is stronger behind the blunt body. Simulations yield highly unsteady forces on the parachute, which was modeled as a rigid body. The mean parachute force behind a slender body is between 3 and 12 percent higher than behind a blunt body, depending on the angle of the parachute with the flow. As the angle of incidence increases, more of the parachute moves out of the leading body wakes, decreasing the sensitivity to leading body shape. To compare the flow past parachutes in Earth's and Mars' atmospheres, simulations were also performed in CO2. At the Mach number considered (1.75), the shock standoff distance ahead of the parachute, post-shock jump conditions, and the resulting parachute forces were found to be very similar in both air and CO2, indicating that a high altitude test is a good proxy for a Mars descent. The results of these numerical simulations and available data on past flight and wind tunnel tests of supersonic Disk-Gap-Band parachutes behind slender bodies were used to generate a parachute drag model for ASPIRE, which in turn was used to help design the flight test. The first flight test occurred in October 2017. The parachute was successfully deployed at Mach 1.77 and an altitude of 42 kilometers. Test instrumentation provided the atmospheric conditions, test vehicle trajectory, and the loads on the parachute along with detailed high-resolution imagery of the inflation process. Reconstruction of the flight test indicated that the measured forces on the parachute were within the model's bounds, although the model over-predicted the parachute force during the first few seconds. The parachute forces during the long subsonic period were well-predicted by the ASPIRE drag model.
The electrostatic potential jump across fast-mode collisionless shocks is examined by comparing published observations, hybrid simulations, and a simple model, in order to better characterize its dependence on the various shock parameters. In all three, it is assumed that the electrons can be described by an isotropic power-law equation of state. The observations show that the cross-shock potential jump correlates well with the shock strength but shows very little correlation with other shock parameters. Assuming that the electrons obey an isotropic power law equation of state, the correlation of the potential jump with the shock strength follows naturally from the increased shock compression and an apparent dependence of the power law exponent on the Mach number which the observations indicate. It is found that including a Mach number dependence for the power law exponent in the electron equation of state in the simple model produces a potential jump which better fits the observations. On the basis of the simulation results and theoretical estimates of the cross-shock potential, it is discussed how the cross-shock potential might be expected to depend on the other shock parameters.
The jump in plasma parameters exhibited by the intense interplanetary shock event of Nov. 12, 1978 is analyzed using ISEE 1, 2 and 3 data. Magnetic and electric field measurements indicated that the shock magnetic field profile was similar to the earth bow shock profile. Data on the electron and proton densities, temperatures, bulk velocities and alpha particles showed a steady electron temperature increase across the shock on a 12 earth radii scale. The upstream and downstream flow parameters are found to be within 10 percent of Rankin-Hugoniot jump conditions. The shock moved at 614 km/sec and had three dissipative scales, one a few Larmor radii determined by the magnetic field jump, a second 10 earth radii correlated with the electron equilibrium and the other 30 earth radii connected to the energetic proton foreshock.
Nonstandard analysis is an area of modern mathematics which studies abstract number systems containing both infinitesimal and infinite numbers. This article applies nonstandard analysis to derive jump conditions for one-dimensional, converging shock waves in a compressible, inviscid, perfect gas. It is assumed that the shock thickness occurs on an infinitesimal interval and the jump functions in the thermodynamic and fluid dynamic parameters occur smoothly across this interval. Predistributions of the Heaviside function and the Dirac delta measure are introduced to model the flow parameters across a shock wave. The equations of motion expressed in nonconservative form are then applied to derive unambiguous relationships between the jump functions for the flow parameters.
The intensities of several lines of molecular hydrogen have been measured from two regions of the supernova-remnant/molecular-cloud shock in IC 443. The lines measured have upper-state energies ranging from 7000 K to 23,000 K. Their relative intensities differ in the two regions, but are consistent with those predicted from the post-shock regions of simple jump-type shocks of different pressure. The pressures so derived are far higher than the pressure in the supernova remnant itself, and a possible reason for this discrepancy is discussed.
Three types of collisionless electrostatic ion acoustic shocks are investigated using a double plasma (DP) device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distrubution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.
Three types of collisionless electrostatic ion-acoustic shocks are investigated using the University of California, Los Angeles, double plasma device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distribution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.
Earth bow wave temperature jump difference between dawn and dusk sides as measured by Vela 3 satellites compared to gasdynamic theory prediction
In the de Hoffmann-Teller reference frame, the cross-shock electric field is simply the thermoelectric field responsible for preserving charge neutrality. As such, it gives information regarding the heating and dissipation occurring within the shock. The total cross-shock potential can be determined by integrating a weighted electron pressure gradient through the shock, but this requires knowledge of the density and temperature profiles. Here, a recently proposed alternative approach relying on particle dynamics is exploited to provide an independent estimate of this potential. Both determinations are applied to slow mode shocks which form the plasma sheet boundary in the deep geomagnetic tail as observed by ISEE 3. The two methods correlate well. There is no indication of the expected transition from resistive to viscous shocks, although the highest Mach number shocks show the highest potentials. The implications of these results for the electron dissipation mechanisms and turbulence at the shock are discussed.
The presence of water/ice/brine in upper layers of Martian crust affects many processes of impact cratering. Modeling of these effects promises better understanding of Martian cratering records. We present here the new ANEOS-based multiphase equation of state for water/ice constructed for usage in hydrocodes and first numerical experiments on permafrost shock melting. Preliminary results show that due to multiple shock compression of ice inclusions in rocks the entropy jump in shocked ice is smaller than in pure ice for the same shock pressure. Hence previous estimates of ice melting during impact cratering on Mars should be re-evaluated. Additional information is included in the original extended abstract.
Jump conditions across MHD shock waves, reducing relations to identities applicable to nonequilibrium situations
This report covers activities funded from October 1, 1998 through September 30, 2002. Two yearly status reports have been filed on this grant, and they are included as Appendix 1. The purpose of this grant was to compare ACE and Wind solar wind parameters when the two spacecraft were near to one another and then to use the intercalibrated parameters to carry out scientific investigations. In September, 2001 a request for a one-year, no-cost extension until September 30, 2002 was submitted and approved. The statement of work for that extension included adjustment of ACE densities below wind speeds of 350 km/s, a study of shock normal orientations using travel time delays between the two spacecraft, comparison of density jumps at shocks, and a study of temperature anisotropies and double streaming to see if such features evolved between the spacecraft.
An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.
An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.
Global, three-dimensional, ideal MHD simulations of Earth's bow shock are reported for low Alfven Mach numbers M(sub A) and quasi-perpendicular magnetic field orientations. The simulations use a hard, infinitely conducting magnetopauause obstacle, with axisymmetric three-dimensional location given by scaled standard model, to directly address previous gasdynamic (GD) and field-aligned MHD (FA-MHD) work. Tests of the simulated shocks' density jumps X for 1.4 approx. less than MA approx. less than 10 and the high M(sub A) shock location, and reproduction of the GD relation between magnetosheath thickness and X for quasi-gasdynamic MHD runs with M(sub A) much greater than M(sub s), confirm that the MHD code is working correctly. The MHD simulations show the standoff distance a(sub s), increasing monotonically with decreasing M(sub A). Significantly larger a(sub s), are found at low M(sub A) than predicted by GD and phenomenological MHD models and FA-MHD simulations, as required qualitatively by observations. The GD and FA-MHD predictions err qualitatively, predicting either constant or decreasing a(sub s), with decreasing M(sub A). This qualitative difference between quasi- perpendicular MHD and FA-MHD simulations is direct evidence for a(sub s), depending on the magnetic field orientation Theta. The enhancement factor over the phenomenological MHD predictions at MA approx. 2.4 agrees quantitatively with one observatiorial estimate. A linear relationship is found between the magnetosheath thickness and X, modified both quantitatively and intrinsically by MHD effects from the GD result. The MHD and GD results agree in the high M(sub A) limit. An MHD theory is developed for a(sub s), restricted to sufficiently perpendicular Theta and high sonic Mach numbers M(sub s). It explains the simulation results with excellent accuracy. Observational and further simulation testing of this MHD theory, and of its predicted M(sub A), Theta, and M(sub s) effects, is desirable.
The electron heating required if protons scatter elastically in a parallel, collisionless shock is calculated. Near-elastic proton scattering off large amplitude background magnetic field fluctuations might be expected if the waves responsible for the shock dissipation are generated by the fire hose instability. The effects of an electrostatic potential jump in the shock layer are included by assuming that the energy lost by protons in traversing the potential jump is converted into electron thermal pressure. It is found that the electron temperature increase is a strong function of the potential jump. Comparison is made to the parallel shock plasma simulation of Quest (1987).