A simplified computational method for studying the heat transfer characteristics of parallel plate thermoacoustic heat exchangers is presented. The model integrates the thermoacoustic equations of the standard linear theory into an energy balance-based numerical calculus scheme. Details of the time-averaged temperature and heat flux density distributions within a representative domain of the heat exchangers and adjoining stack are given. The effect of operation conditions and geometrical parameters on the heat exchanger performance is investigated and main conclusions relevant for HX design are drawn as far as fin length, fin spacing, blockage ratio, gas and secondary fluid-side heat transfer coefficients are concerned. Most relevant is that the fin length and spacing affect in conjunction the heat exchanger behavior and have to be simultaneously optimized to minimize thermal losses localized at the HX-stack junctions. Model predictions fit experimental data found in literature within 36% and 49% respectively at moderate and high acoustic Reynolds numbers. (C) 2011 Elsevier Ltd. All rights reserved.