On the solution of the Boltzmann equation for a two-dimensional gas.
Boltzmann equation and statistical properties for two-dimensional gas, analyzing integral iteration for shock wave flow
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.
Boltzmann equation and statistical properties for two-dimensional gas, analyzing integral iteration for shock wave flow
Pulsed gas dynamic laser for radiation generation at high power levels, using shock tube
Pioneer Venus magnetometer observations in the sunlit ionosphere indicate that the ionosphere is dynamic and very responsive to external solar wind conditions. Bow shock location, ionosphere location, the strength of the magnetic field just outside the ionopause, and the field strength in the ionosphere are found to be variable, and the properties of flux ropes in the ionospheric magnetic field are considered. Data on magnetic energy density and on magnetic field strength are presented.
The DAST Aeroelastic Research Wing had been previously in the NASA Langley TDT and an unusual instability boundary was predicted based upon supercritical response data. Contrary to the predictions, no instability was found during the present test. Instead a region of high dynamic wing response was observed which reached a maximum value between Mach numbers 0.92 and 0.93. The amplitude of the dynamic response increased directly with dynamic pressure. The reponse appears to be related to chordwise shock movement in conjunction with flow separation and reattachment on the upper and lower wing surfaces. The onset of flow separation coincided with the occurrence of strong shocks on a surface. A controller was designed to suppress the wing response. The control law attenuated the response as compared with the uncontrolled case and added a small but significant amount of damping for the lower density condition.
A thermo-chemical nonequilibrium ionized nitrogen flow in a shock tube is examined. A one-dimensional computational fluid dynamics code has been developed to study the flowfield incorporating a novel treatment of the electron energy equation. In the present approach, the electron pressure is included in the electron energy flux. In contrast, previous work has removed this term from the flux to facilitate the flux splitting. The code uses a five-species gas model for ionized nitrogen characterized by translational-rotational, vibrational, and electron-electronic temperatures. The results give good agreement with experimental data except for an inadequacy in the electronic source terms. A small difference is observed between the new and old splittings which increases with higher levels of ionization.
Acceleration of interstellar pickup H(+) and He(+) as well as of solar wind protons and alpha particles has been observed on Ulysses during the passage of a corotating interaction region (CIR) at approximately 4.5 AU. Injection efficiencies for both the high thermal speed interstellar pickup ions (H(+) and He(+)) and the low thermal speed solar wind ions (H(+) and He(++) are derived using velocity distribution functions of protons, pickup He(+) and alpha particles from less than 1 to 60 keV/e and of ions (principally protons) above approximately 60 keV. The observed spatial variations of the few keV and the few hundred keV accelerated pickup protons across the forward shock of CIR indicate a two stage acceleration mechanism. Thermal ions are first accelerated to speeds of 3 to 4 times the solar wind speed inside the CIR, presumably by some statistical mechanism, before reaching higher energies by a shock acceleration process. Our results also indicate that (1) the injection efficiencies for pickup ions are almost 100 times higher than they are for solar wind ions, (2) pickup H(+) and He(+) are the two most abundant suprathermal ion species and they carry a large fraction of the particle thermal pressure, (3) the injection efficiency is highest for protons, lowest for He(+), and intermediate for alpha particles, (4) both H(+) and He(+) have identical spectral shapes above the cutoff speed for pickup ions, and (5) the solar wind frame velocity distribution function of protons has the form F(w) = F(sub o)w(sup -4) for 1 is less than w is less than approximately 5, where w is the ion speed divided by the solar wind speed. Above w approximately 5-10 the proton spectrum becomes steeper. These results have important implications concerning acceleration of ions by shocks and CIRs, acceleration of anomalous cosmic rays, and particle dynamics in the outer heliosphere.
In this paper, magnetic and plasma measurements are used to analyze 17 interplanetary coronal mass ejections (CMEs) identified by Ulysses during its in-ecliptic passage to Jupiter. We focus on the expansion characteristics of these CMEs (as inferred from the time rate of change of the velocity profiles through the CMEs) and the properties of 14 forward shocks unambiguously associated with these CMEs. We highlight radial trends from 1 to 5.4 AU. Our results indicate that the CMEs are generally expanding at all heliocentric distances. With regard to the shocks preceding these ejecta, we note the following: (1) There is a clear tendency for the shock speed (in the upstream frame of reference) to decrease with increasing heliocentric distance as the CMEs transfer momentum to the ambient solar wind and slow down; (2) 86% of the shock fronts are oriented in the ecliptic plane such that their normals point westward (i.e., in the direction of planetary motion about the Sun), (3) 86% of the shocks are propagating toward the heliographic equator; and (4) no clear trend was found in the strength of the shocks versus heliocentric distance. These results are interpreted using simple dynamical arguments and are supported by fluid and magnetohydrodynamic (MHD) simulations.
The Ares I-X flight test is the first development flight of the Ares I crew launch vehicle. This mission, first conceived in 2006, will be launching later in 2009. Its primary mission objectives will be to demonstrate flight and roll control of a dynamically similar vehicle, perform a separation event and measure its shock effects, stack and recover a first stage booster, and demonstrate ground operations. All of the primary flight test vehicle s hardware is at Kennedy Space Center, and is being stacked in the Vehicle Assembly Building for a liftoff at Launch Complex 39B. Mission hardware specific to Ares I-X also is being installed at Launch Complex 39, which has been supporting Space Shuttle operations. This presentation will provide a status and preview of the upcoming mission.
Large-scale fluctuations in the solar wind upstream of the termination shock will cause inward and outward motions of the shock. In earlier work, Barnes analyzed such motion by calculating of the response of a planar gasdynamic shock to upstream disturbances. We now generalize this analysis to the case of a spherically symmetric shock. Our procedure is first to solve numerically the set of gasdynamic equations describing the interaction between the solar wind and the interstellar medium to establish a dynamic equilibrium. The next step is to impose upstream fluctuations of the solar wind dynamical pressure on this equilibrium state at an inner boundary, and then to follow the subsequent shock motion.
A study is presented on the effect of embedded supersonic flows and the resulting emerging shock waves on phenomena associated with MHD turbulence, including reconnection, the formation of current sheets and vortex structures, and the evolution of spatial and temporal correlations among physical variables. A two-dimensional model problem, the Orszag-Tang (1979) vortex system, is chosen, which involves decay from nonrandom initial conditions. The system is doubly periodic, and the initial conditions consist of single-mode solenoidal velocity and magnetic fields, each containing X points and O points. The initial mass density is flat, and the initial pressure fluctuations are incompressible, balancing the local forces for a magnetofluid of unit mass density. Results on the evolution of the local structure of the flow field, the global properties of the system, and spectral correlations are presented. The important dynamical properties and observational consequences of embedded supersonic regions and emerging shocks in the Orszag-Tang model of an MHD system undergoing reconnection are discussed. Conclusions are drawn regarding the effects of local supersonic regions on MHD turbulence.
An analytic approach is demonstrated to reveal potential pyroshock-driven dynamic effects causing power losses in the Thermo-Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi-Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high-fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large-amplitude suddenly-applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees-of-freedom (DOF). First, a quasi-static preload was applied on top of the TE module bar, then transient tri-axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high-fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.
Radiative temperatures and electrical conductivities were measured for fluid nitrogen compressed dynamically to pressures of 18-90 GPa, temperatures of 4000-14,000 K, and densities of 2-3 g/cu cm. The data show a continuous phase transition above 30 GPa shock pressure and confirm that (delta-P/delta-T)v is less than 0, as indicated previously by Hugoniot equation-of-state experiments. The first observation of shock-induced cooling is also reported. The data are interpreted in terms of molecular dissociation, and the concentration of dissociated molecules is calculated as a function of density and temperature.
Nonstandard analysis is a relatively new area of mathematics in which infinitesimal numbers can be defined and manipulated rigorously like real numbers. This report presents a fairly comprehensive tutorial on nonstandard analysis for physicists and engineers with many examples applicable to generalized functions. To demonstrate the power of the subject, the problem of shock wave jump conditions is studied for a one-dimensional compressible 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. To use conservations laws, smooth pre-distributions of the Dirac delta measure are applied whose supports are contained within the shock thickness. Furthermore, smooth pre-distributions of the Heaviside function are applied which vary from zero to one across the shock wave. It is shown that if the equations of motion are expressed in nonconservative form then the relationships between the jump functions for the flow parameters may be found unambiguously. The analysis yields the classical Rankine-Hugoniot jump conditions for an inviscid shock wave. Moreover, non-monotonic entropy jump conditions are obtained for both inviscid and viscous flows. The report shows that products of generalized functions may be defined consistently using nonstandard analysis; however, physically meaningful products of generalized functions must be determined from the physics of the problem and not the mathematical form of the governing equations.
Research on the topic of shock wave/turbulent boundary-layer interaction was carried out during the past three years at the Penn State Gas Dynamics Laboratory. This report describes the experimental research program which provides basic knowledge and establishes new data on heat transfer in swept shock wave/boundary-layer interactions. An equilibrium turbulent boundary-layer on a flat plate is subjected to impingement by swept planar shock waves generated by a sharp fin. Five different interactions with fin angle ranging from 10 deg to 20 deg at freestream Mach numbers of 3.0 and 4.0 produce a variety of interaction strengths from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface, and miniature thin-film-resistance sensors mounted on it are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding a total uncertainty of about +/- 10 percent. These experimental data are compared with the results of numerical Navier-Stokes solutions which employ a k-epsilon turbulence model. Finally, a simplified form of the peak heat transfer correlation for fin interactions is suggested.
Understanding the large-scale dynamics of the magnetospheric boundary is an important step towards achieving the ISTP mission's broad objective of assessing the global transport of plasma and energy through the geospace environment. Our approach is based on three-dimensional global magnetohydrodynamic (MHD) simulations of the solar wind-magnetosphere- ionosphere system, and consists of using interplanetary magnetic field (IMF) and plasma parameters measured by solar wind monitors upstream of the bow shock as input to the simulations for predicting the large-scale dynamics of the magnetospheric boundary. The validity of these predictions is tested by comparing local data streams with time series measured by downstream spacecraft crossing the magnetospheric boundary. In this paper, we review results from several case studies which confirm that our MHD model reproduces very well the large-scale motion of the magnetospheric boundary. The first case illustrates the complexity of the magnetic field topology that can occur at the dayside magnetospheric boundary for periods of northward IMF with strong Bx and By components. The second comparison reviewed combines dynamic and topological aspects in an investigation of the evolution of the distant tail at 200 R(sub E) from the Earth.
The proceedings of the 57th Shock and Vibration Symposium are presented. In part 1 is presented the welcoming and the keynote address along with invited papers such as, Dynamic Testing - Seven Years Later. Papers from the Nondevelopment Items workshop such as, Guidelines for Qualifying Nondevelopment Equipment to Shock and Vibration, are presented. And from the Pyrotechnic Shock Workshop the paper entitled, Designing Electronics for Pyrotechnic Shock, is presented.
During the past two decades, there has been an increased reliance on the use of computational fluid dynamics methods for modeling rotors in high speed forward flight. Computational methods are being developed for modeling the shock induced loads on the advancing side, first-principles based modeling of the trailing wake evolution, and for retreating blade stall. The retreating blade dynamic stall problem has received particular attention, because the large variations in lift and pitching moments encountered in dynamic stall can lead to blade vibrations and pitch link fatigue. Restricting to aerodynamics, the numerical prediction of dynamic stall is still a complex and challenging CFD problem, that, even in two dimensions at low speed, gathers the major difficulties of aerodynamics, such as the grid resolution requirements for the viscous phenomena at leading-edge bubbles or in mixing-layers, the bias of the numerical viscosity, and the major difficulties of the physical modeling, such as the turbulence models, the transition models, whose both determinant influences, already present in static maximal-lift or stall computations, are emphasized by the dynamic aspect of the phenomena.
A report is presented regarding the dynamic properties of a rock indigenous to the mare basins of the moon. The reported data were obtained in a study of sample 70215, a very titanium-rich basalt (58% pyroxene, 18% ilmenite, 15% plagioclase, 6% olivine, and 3% quartz by weight). This rock is probably representative of a class of the earliest mare-filling extrusive rocks which are exposed on the present lunar surface. Two series of experiments were performed. One set of experiments involved the measuring of Hugoniot and release adiabats to 15.7 GPa with a propellant gun apparatus. In the second set of experiments, a light-gas gun was employed to yield Hugoniot data at about 120 GPa and release states at about 90 GPa. Lunar basalt 70215 appears to be among the densest rocks in the present lunar sample collection, having a crystal density of 3.38 g/cu cm and a porosity of about 1.3%. The results of the experiments have important implications for both the degree of shock metamorphism expected for impact processes and the extent of ejecta transport on mare surfaces with high-titanium basalt composition.