Rj. Pearson et De. Winterbone, CALCULATING THE EFFECTS OF VARIATIONS IN COMPOSITION ON WAVE-PROPAGATION IN GASES, International journal of mechanical sciences, 35(6), 1993, pp. 517-537
The equations of gas dynamics in one space dimension are solved whilst
considering variations in gas composition and properties. A species c
ontinuity equation is derived and coupled in a single vector equation
with the three usual conservation equations; the resulting system is s
olved directly in conservation law form using the two-step Lax-Wendrof
f technique. The effects of variable gas properties and composition on
wave propagation are illustrated via the shock-tube, or Riemann, prob
lem. Significant errors are introduced, in the cases investigated, by
assuming that the fluid is a perfect gas of constant composition. Flux
-corrected transport (FCT) and nonupwind total variation diminishing (
TVD) approaches are evaluated as a means of mitigating the spurious os
cillations produced at discontinuities by the classical Lax-Wendroff s
cheme; these oscillations can give rise to mass fraction values from t
he species equation which are greater than unity or less than zero. It
is found that the FCT algorithm does. not suppress completely the osc
illations at discontinuities and this corrupts the species transport c
alculations. The TVD algorithm eliminates the oscillations at shock wa
ves and contact surfaces in the cases tested, thus maintaining the int
egrity of the species calculations. This is achieved at the expense of
a 65% increase in computational effort over the constant gas property
, constant composition case. The resolution of both the shock wave and
contact surface can be improved significantly by the use of the artif
icial compression technique, but at an additional computational cost.