CALCULATING THE EFFECTS OF VARIATIONS IN COMPOSITION ON WAVE-PROPAGATION IN GASES

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.