航空学报 > 2016, Vol. 37 Issue (7): 2091-2101   doi: 10.7527/S1000-6893.2015.0207

1/4波长驻波型热声发动机的非线性模型

侯薇1, 王晓宇2, 景晓东2, 孙晓峰2   

  1. 1. 中航空天发动机研究院有限公司, 北京 101304;
    2. 北京航空航天大学 能源与动力工程学院, 北京 100083
  • 收稿日期:2015-06-03 修回日期:2015-07-20 出版日期:2016-07-15 发布日期:2015-08-03
  • 通讯作者: 景晓东 男,博士,教授,博士生导师。主要研究方向:气动声学。Tel:010-82338085,E-mail:jingxd@buaa.edu.cn E-mail:jingxd@buaa.edu.cn
  • 作者简介:侯薇 女,博士,工程师。主要研究方向:热声非线性模型。Tel:010-56680655,E-mail:buaahouwei@163.com;王晓宇 女,博士,讲师。主要研究方向:气动声学。Tel:010-82338262,E-mail:bhwxy@buaa.edu.cn;孙晓峰 男,博士,教授,博士生导师。主要研究方向:气动声学和叶轮机稳定性。Tel:010-82317408,E-mail:sunxf@buaa.edu.cn
  • 基金资助:

    国家973计划(2012CB720200);国家自然科学基金(51236001)

A quasi-one-dimensional nonlinear model of an open-closed standing-wave thermoacoustic engine

HOU Wei1, WANG Xiaoyu2, JING Xiaodong2, SUN Xiaofeng2   

  1. 1. AVIC Academy of Aeronautic Propulsion Technology, Beijing 101304, China;
    2. School of Energy and Power Engineering, Beihang University, Beijing 100083, China
  • Received:2015-06-03 Revised:2015-07-20 Online:2016-07-15 Published:2015-08-03
  • Supported by:

    National Basic Research Program of China (2012CB720200); National Natural Science Foundation of China (51236001)

摘要:

为研究声波在热声系统开口端和负载构成的阻抗间断面上的声反射问题,针对一端开口的1/4波长驻波型热声发动机建立了一种基于截面平均的准一维非线性热声模型,并采用具有频散保持特性的高精度计算格式进行数值求解。在发动机开口端采用了对于系统共振频率具有一定选择自适应性的宽频时域声阻抗边界条件。宽频阻抗模型写成共轭复数对的留数和极点分式之和的形式,其时域形式可用递推卷积方法进行快速高效的求解,共轭复数系数采用优化拟合方法近似得出。模型可对一端开口的1/4波长驻波型发动机进行临界起振温度的预测,能够描述热声系统的瞬态压力波由初始小扰动经过非线性增长,最终达到饱和的全过程,计算结果与某微小型驻波热声发动机的实验结果相符。最后分析讨论了板叠的形状对热声系统的影响。

关键词: 热声发动机, 非线性, 模型, 宽频, 时域声阻抗

Abstract:

A quasi-one-dimensional nonlinear model of a standing-wave thermoacoustic engine with a quarter-wavelength resonator is presented. The model is based on cross-sectional averaged equations and is solved by high-order low-dispersion numerical schemes. Considering the acoustic reflection and radiation at the open end, broadband time-domain impedance boundary conditions are employed so that oscillation frequency can be self-adaptive. The impedance is expressed in the mathematical form of partial fraction expansion with complex-conjugate residues and poles, so that the convolution of the impedance with the velocity can be calculated by efficient and causal recursive convolution. Pole-residue pairs of the impedance can be optimized in the frequency domain. Critical temperature difference of an open-ended standing-wave thermoacoustic engine can be predicted and the time-domain simulation shows that the pressure oscillation undergoes a nonlinear amplifying process and eventually reaches saturation amplitude. Numerical results agree very well with the experiment for a small-scale thermoacoustic engine. The simulation results for different stack geometries are included.

Key words: thermoacoustic engines, nonlinear, model, broadband, time-domain impedance boundary condition

中图分类号: