In terms of numerical implementation of the algorithm, the time domain pressure travelling waves are stored in vectors with elements referred to by an index (which we shall call ) which runs from 1 to . Care must be taken over what time is represented by the index. On surfaces 1 and 2 we know the pressure at , , , ..., so the th element of the vector in the numerical implementation therefore corresponds to the pressure at time .

Equations (5.9) and (5.15) are
then simply

The forward going pressure wave on surface 3 is sampled
at , , , ..., and the th element refers to
the time . Comparison with equation (5.17) then yields

The pressure waves and are then the forward going and backward going sequences for the left hand side of the junction at surface 3. We have therefore peeled off the first layer in the object. Equations (5.21) to (5.24) can then be used to peel off the next layer. First the primary reflection coefficient for the junction from surface 3 to surface 4 is calculated, then the forward and backward travelling waves on surface 4 are obtained and then the delay for travel along the cylinder to next discontinuity is added. The algorithm repeats as necessary, peeling off successive layers in the object, calculating the new cross-section using equation (5.10) each time.

Back to Kemp Acoustics Home