Numerical computation for parallel plate thermoacoustic heat exchangers in standing wave oscillatory flow