Numerical simulation on onset characteristics of traveling-wave thermoacoustic engines based on a time-domain network model

• Main Authors: Wang, Kai, Zhang, Jie, Zou, Jiang, Wu, Ke, Qiu, Li Min, Huang, Z. Y., Sun, Da Ming
• Other Authors: Energy Research Institute @ NTU (ERI@N)
• Format: Article
• Language: English
• Published: 2015
• Subjects: DRNTU::Science::Physics::Heat and thermodynamics
• Online Access: https://hdl.handle.net/10356/103665; http://hdl.handle.net/10220/25842
• Institution: Nanyang Technological University

Onset characteristics of thermoacoustic engines are of great importance for understanding the internal working mechanisms of thermoacoustic conversion. A one-dimensional time-domain network model for predicting the onset characteristics of traveling-wave thermoacoustic engines with helium as working gas is built. The acoustic resistance, inertance, compliance, and thermal-relaxation effects of all the acoustic components are included. The viscous and heat transfer terms in the time-domain governing equations of the acoustic tubes and the heat exchangers are deduced from the frequency-domain linear thermoacoustic theory. Combining the time-domain governing equations of the regenerator, numerical simulations of the whole onset process are then conducted in a wide operating condition range. The complete dynamic pressure wave evolution processes are simulated successfully. It is shown that a steady standing-wave acoustic field forms in almost all parts of the traveling-wave thermoacoustic engine except for the regenerator area. Onset temperature, operating frequency, and quality factor are calculated with a relatively high accuracy. The thermal relaxation effects in the regenerator are found to have a remarkable impact on the onset characteristics, especially at high mean pressures. It is also shown that the experimental damping temperature is closer to the calculated onset temperature than the experimental onset temperature. Furthermore, the reasonable distributions of the pressure and volume flow rate and the phase relationship between them in the whole system are obtained and analyzed.