Hostname: page-component-cd9895bd7-dk4vv Total loading time: 0 Render date: 2024-12-26T07:58:34.728Z Has data issue: false hasContentIssue false
  • English
  • Français

The internal structure of shock waves

Published online by Cambridge University Press:  29 March 2006

Bruce L. Hicks ,
Shee-Mang Yen  and
Barbara J. Reilly
Bruce L. Hicks
Affiliation:
Co-ordinated Science Laboratory, University of Illinois
Shee-Mang Yen
Affiliation:
Co-ordinated Science Laboratory, University of Illinois
Barbara J. Reilly
Affiliation:
Co-ordinated Science Laboratory, University of Illinois
Get access Rights & Permissions [Opens in a new window]

Abstract

The nonlinear Boltzmann equation has been solved for shock waves in a gas of elastic spheres. The solutions were made possible by the use of Nordsieck's Monte Carlo method of evaluation of the collision integral in the equation. Accurate solutions were obtained by the same numerical procedure for eight values of the upstream Mach numbers M1 ranging from 1.1 to 10, even though the corresponding degree of departure from equilibrium varies by a factor greater than 100. Many more characteristics of the internal structure of the shock waves have been calculated from the solutions than have hitherto been available. Each solution of the Boltzmann equation requires about 108 multiplications to obtain statistical errors of 3% in values of the velocity distribution function and collision integral and much smaller errors in the moments of these functions. The reciprocal shock thickness is in agreement with that of the Mott-Smith shock (u2 moment) from M1 = 2.5-8. The density profile is asymmetric with an upstream relaxation rate (measured as density change per mean free path) approximately twice as large as the downstream value for weak shocks and equal to the downstream value for strong shocks. The temperature density relation is in agreement with that of the Navier-Stokes shocks for Mach numbers in the range 1·1-1-56. The Boltzmann reciprocal shock thickness is smaller than the Navier-Stokes value in this range of Mach number because the viscosity-temperature relation computed is not constant as predicted by the linearized theory.

The velocity moments of the distribution function are, like the Mott-Smith shock, approximately linear with respect to the number density; however, the deviations from linearity are statistically significant. Four functionals of the distribution function that are discussed show maxima within the shock. The entropy is a good approximation to the Boltzmann function for all M1. The solutions obtained satisfy the Boltzmann theorem for all Mach numbers. The ratio of total heat flux q to qx (that associated with the longitudinal degree of freedom) correlates well with local Mach number for all Ml in accordance with a relation derived by Baganoff & Nathenson (1970). The Chapman-Enskog linearized theory predicts that this ratio is constant. The (effective) transport coefficients are larger than the Chapman-Enskog equivalents by as much as a factor of three at the mid-shock position.

At M1 = 4, and for 40% of the velocity bins, the distribution function is different from the corresponding Mott-Smith value by more than three times the 90% confidence limit. The r.m.8. value of the percentage difference in distribution functions is 15% for this Mach number. At MI = 1-59, the half width and several other characteristics of the function \[ \int f\,dv_y dv_z \] differ from that of the Chapman-Enskog first iterate, and many of the deviations are in agreement with an experiment by Muntz & Harnett (1970).

Type
Research Article
Information
Journal of Fluid Mechanics , Volume 53 , Issue 1 , 9 May 1972 , pp. 85 - 111
Copyright
© 1972 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Bauanoff, D. & Natbenson, M. 1970 Phys. Fluids, 13, 596.
Bird, G. A. 1965 Proc. 4th Int. Symp. on Rarefied Gaa Dynamics, 1, 216.
Bird, G. A. 1967 J. Fluid Mech. 30, 479.
Bird, G. A. 1970a Phys. Fluids, 13, 1172.
Bird, G. A. 1970b Phys. Fluids, 13, 2676.
Grad, H. 1952 Comm. Pure Appl. Math. 5, 257.
Hicks, B. L. 1965 C.S.L. Rep. University of Illinois, R-236.
Hicks, B. L. & Smith, M. A. 1967 C.S.L. Rep. University of Illinois, R-347.
Hicks, B. L. & Smith, M. A. 1968 J. Comp. Phys. 3, 58.
Hicks, B. L. & YEN, S. M. 1967 Phys. Fluids, 10, 458.
Hicks, B. L. & YEN, S. M. 1969 Proc. 6th Int. Symp. on Rarefied Gas Dynamics, 1, 313.
Hicks, B. L. & YEN, S. M. 1971 Proc. 7th Int. Symp. on Rarefied Glas Dynamics (to be published).
Hicks, B. L., YEN, S. M. & Reilly, B. 1969 CSL Rep. University of Illinois, R-412.
Holtz, T., Mwtz, E. P. & YEN, S. M. 1971 Phys. Fluids, 14, 545.
Mordtjceow, M. & Libby, P. A. 1962 Pibol Rep. Brooklyn Institute of Technology, no. 749.
Mott-Smith, H. M. 1951 Phys. Rev. 82, 885.
Muntz, E. P. & Harnett, L. N. 1970 Phys. Fluids, 12, 2027.
Nordsieck, A. & Hicks, B. L. 1967 Proc. 5th Inst. Symp. on Rarefied Gaa Dynamics, 1 675.
Schmidt, B. 1969 J. Fluid Mech. 39, 361.
Schmidt, H. J. 1965 M.S. thesis, University of Illinois.
Talbot, L. & Sherman, F. S. 1959 N.A.S.A. Memo. 12-14-58W.
Wang-Chang, C. B. 1948 Univers & of Michigan Rep. APL/JHO CM 504.
Yen, S.M. 1966 Phys. Fluids, 9, 1417.
Yen, S. M. 1971 Proc. 7th Int. Symp. on Rarefied Gas Dynamics (to be published).
Yen, S. M. & Hicks, B. L. 1967a CSL Rep. University of Illinois, R-350.
Yen, S. M. & Hicks, B. L. 1967b Proc. 5th Int. Symp. on Rarefied Gas Dynamics, 1, 785.
Yen, S. M. & Schmidt, H. J. 1969 Proc. 6th Int. Symp. on Rarefied Gas Dynamics, 1, 205.