广义一维势热声发动机微循环的性能分析
鄂青1, 2，吴锋1, 2，雷霆1
(1. 武汉工程大学 光电信息与能源工程学院，湖北 武汉，430205；
2. 武汉工程大学 热科学与动力工程研究所，湖北 武汉，430205)
摘要：对热声发动机中的工质粒子进行量子力学行为分析，得出粒子处在不同一维势场条件下时的广义量子热声发动机微循环分析模型；运用此模型推导出广义量子热声发动机的性能参数表达式。研究结果表明：要使热声发动机性能达最优，必须对声场进行控制，使其能够在回热器中建立起一维无限深的或谐振的势阱条件。
关键词：广义一维势；量子热声发动机；循环性能
中图分类号：TB65 文献标志码：A
文章编号：1672-7207（2021）06-1883-09
Performance analysis of micro-cycle in quantum thermoacoustic engine with generalized 1D potential
E Qing1, 2, WU Feng1, 2, LEI Ting1
(1. School of Optical Information and Energy Engineering, Wuhan Institute of Technology, Wuhan 430205, China;
2. Institute of Thermal Science and Power Engineering, Wuhan Institute of Technology, Wuhan 430205, China)
Abstract: The quantum mechanical behaviour of working fluid particles in thermoacoustic engine was analyzed, and a set of analysis model for generalized quantum thermoacoustic micro-cycle working in various one-dimensional(1D) potential wells was established. The expression of performance parameters of generalized quantum thermoacoustic engine was derived from this model. The results show that to achieve the optimal performance of the thermoacoustic engine, the acoustic field must be controlled so that 1D infinite or resonant potential well condition can be established in the regenerator.
Key words: generalized 1D potential; quantum thermoacoustic engine; performance of micro-cycle
有限时间热力学(finite time thermodynamics，FTT)[1-5]作为经典热力学的延伸，主要研究非平衡系统中能量流和熵流的规律。在低温系统和微纳米能量转换装置的进一步研究和工程实践中[6-7]，系统的量子特性变得越来越重要。在有限时间热力学与量子热力学的基础上，人们建立并分析了不同种类的量子能量转换系统，这代表了热力学的一个重要的发展方向。
人们在研究量子卡诺[8]、奥托[9]、斯特林[10-11]、迪塞尔[12]和热声[13]循环时，将其工质视为一系列粒子的集合。这些粒子有的是囚禁于一个一维无限深势阱的，也可以是被囚禁于谐振势或者四次势的，甚至可以是箱势中的极端相对论粒子[8-17]。王建辉等[18-19]验证了循环性能计算中粒子的能量谱简化式，并据此得出了广义量子奥托循环的效率及做功量的表达式。本文作者进一步对由2个等压过程和2个等熵过程组成的热声发动机微循环展开分析。
1 广义量子热声发动机微循环模型的建立
1.1 系统的量子力学基础
普朗克的量子理论给出了粒子的能量表达式为
。可将处于1个广义一维势场内的粒子所具有的能量谱公式简化[18]为
(1)
式中：
为粒子的能量；
为约化普朗克数；
，是为了满足量子化的能量表达形式而构建出的频率；L为势阱宽度；
为特定势场所固有的常数；
和
为由势阱种类决定的指数；n为量子数。，和在不同势阱条件下的取值见表1。表1中，c和m分别为粒子的速度与质量；L1为系统处于所示循环的宏观状态1时粒子所处势阱的宽度；
为系统处于宏观状态1时其中处于基态(n=1)的粒子能量。这里选用了4种特殊势：一维无限深势阱、谐振势阱、包含相对论粒子的一维势阱以及四次势阱，它们代表了广义一维势的典型情况。系统在宏观状态i的内能
计算式为
(2)
式中：
为处于状态i时n级本征态的粒子的能量；
为状态i时处于n级本征态的粒子的概率。pin应满足如下归一化条件：
(3)
表1 不同势阱的参数比较
Table 1 Comparisons of different potential parameters
在经典热声热力过程中，外界通过振动膜片的运动对系统做功。与振动膜片的经典运动类似，可以假设广义一维势阱的壁面在有限的速度下运动，这样，当系统消耗外界功量或对外界输出功量时，其施加于势阱壁面的合力Fi就可以写为
(4)
式中：Fin为处于状态i时n级本征态的粒子施加于势阱壁面的力；Li为处于状态i时的势阱宽度。
将式(1)代入式(4)可得
(5)
1.2 理想广义量子热声发动机微循环
图1 理想热声发动机微循环示意图
Fig. 1 Schematic diagrams of ideal thermoacoustic engine micro-cycle
图1所示为在理想情况下，热声回热器间气体微团与固体工质的热作用微循环过程。其中，热声板叠两端存在明显的温度梯度，板间距小于或等于气体工质热渗透深度的2倍[20]。设Tx为气团在平衡位置处的平均温度(下标x表示微团的当地坐标)，Tm为温度幅值，p0为系统静压，pm为压力波动幅值，ξm为振荡的位移幅值，
为热声回热器固体壁的纵向温度梯度，
为气团处于状态i(i=1，3)时的温度与平均温度偏差的绝对值。当气体微团处于最右端状态1时，气团的温度
低于与之相接触的器壁温度为
，因而有热量Q1从固体壁面流向气团使其温度升高，变为Tx+Tm，并因膨胀而对邻近气团做功。这一过程可近似为定压吸热过程，即图1中的过程1—2。当气团由于振荡和热弛豫快速从右端向左端运动时，此气团无法与器壁进行热交换，因而可近似认为该过程是绝热的。当气团到达左端时，温度变为
，压力从右端的p0+pm变为p0-pm，与气团相接触的当地壁温为
，气团由于膨胀而对邻近气团做功，这就是图1中的过程2—3。当气团到达左端时，由于气团温度比当地壁温高，因而有热量Q2流出气团，使其温度降低为Tx-Tm，并因被压缩而消耗邻近气团的功。这一过程可近似为定压放热过程，即图1中的过程3—4。图1(b)中的过程4—1为近似的绝热压缩过程，当气团因振荡而从左端向右端快速运动时，由于热弛豫，它无法与器壁交换热量。当到达右端时，气团的温度变为，压力变为p0+pm，此时，气团的状态回复到初始状态，完成了1个热力学循环。根据循环命名传统，可将热声发动机微循环简称为热声微循环。由热力学第一定律可知，在循环的每个过程中，气体微团作为闭口系统应满足能量守恒方程
(其中，Q和W分别为在任意过程中热力系与外界交换的热量与功，
为系统内能)。当系统自外界吸热时，Q>0；当向外界放热时，Q<0；当系统消耗外功时，W<0；对外界做功时，W>0。
本文从量子力学的角度分析，可将上述气体微团视为1种被限制在广义一维势阱中的粒子。为了简单，在分析量子热声微循环的过程中，只考虑由出现概率较高的2个特征态粒子构成的二能级系统。1台真正的发动机中的工作介质是由无数这样的粒子组成的。从粒子的量子行为角度分析，每个微循环都可归纳为由2个量子绝热过程和2个量子等压过程环绕而成，由此得到本文的主要研究对象即理想广义量子热声微循环(Ideal generalized quantum thermoacoustic micro-cycle, IGQTAC)，如图2所示。
图2 理想广义量子热声微循环示意图
Fig. 2 Schematic diagram of ideal generalized quantum thermoacoustic micro-cycle
在过程1—2中，系统与1个功源及1个高温热源(温度为TH)相偶联。随势阱壁外推，系统在相同作用力下膨胀，此时，它会从高温热源吸收大量热量，同时向功源输出体积功。通常势阱壁的运动耗时比高温热源的量子弛豫时间(约
)长得多，即系统能够保持与高温热源的充分热接触，以确保其始终处于热平衡状态。在系统中，处于激发态(n=2)的粒子的占有概率pie与处于基态(n=1)的粒子的占有概率pig=1-pie之间应满足玻尔兹曼分布规律：
(6)
式中：k为玻尔兹曼常数；Ti为系统平衡温度。
由式(6)可求得系统处于激发态的概率为
(7)
式中：
，为系统处于宏观状态i时的能级宽度。
由式(7)可得
(8)
联立式(1)，(3)和(5)可得出在此过程中系统做功W12的计算式为
(9)
由于系统始终处于热平衡状态，所以，气团从高温热源的吸收的热量Q12为
(10)
式中：(i=1，2)为图2所示状态1和状态2的系统内能。由式(1)~(3)可得
(11)
将式(9)和式(11)代入式(10)可得
(12)
在过程2—3中，系统与外界没有热交换，系统对外界做的功将全部来自系统自身内能。在量子绝热过程中，系统内部粒子不会发生能级跃迁，即各能级粒子的占有概率维持不变，即
。由式(1)~(3)和式(5)可得交换的功W23为
(13)
在过程3—4中，系统在等作用力状态下向低温热源(温度为TL)释放热量。同时，气团自身体积随势阱壁的内移被压缩，从而消耗外界功。采用类似过程1—2的分析方法可得此过程中交换的功W34和热量Q34分别为：
＜0(14)
(15)
在过程4—1中，由系统绝热压缩有
，故消耗的功W41为
(16)
根据热力学第一定律，在理想广义量子热声微循环过程中，系统从高温热源吸收的热量QH为
(17)
系统输出净功Wnet为
＞0(18)
1.3 循环周期
为了计算完整的循环周期，可假设势阱宽度变化的平均速度为
，这里，势阱宽度的平均变化速度很慢，因此其变化时间远大于系统及热源的弛豫时间(约为)，从而使系统能够始终处于热平衡状态下。则系统经历循环过程1—2—3—4—1所需的时间，即循环周期
[21]为
(19)
1.4 循环性能参数
由式(1)~(3)可得1个二能级系统的内能为
(20)
由式(5)可知：1个二能级系统施加在势阱壁上的力
可写为
(21)
由
和
，可得
(22)
(23)
由式(7)可知当系统处在宏观状态2和状态4时，其内部激发态粒子的占有概率分别为：
(24)
(25)
式中：
>1，为势阱宽度比；T2和T4分别为状态点2和状态点4的系统温度。
在理想情况下，可取状态点4和状态点2的温度分别为冷、热端温度(图2)，即
显然，对于特定势场，p4e=p1e，它们是由L4，Tx，Tm及决定的。对于p4e，以及已给定的情况，仅由
决定，因此，是1个重要的控制参数。
将式(22)和式(23)代入式(17)和式(18)可将循环从高温热源获得的热量及做功计算式改写为：
即可得出广义量子热声微循环的效率η为
(26)
将式(22)和式(23)代入式(19)得
循环的输出功率为
(27)
令
，可得量纲一输出功率
为
(28)
对于任意特定的势场，和都是常数。当冷热端温度T4和T2已知时，由式(25)可知，式(26)和式(28)所描述的效率与量纲一功率都仅仅是p4e和势阱宽度比的函数。
2 广义量子热声微循环性能分析
式(26)和式(28)指出了广义量子热声微循环的性能与系统参数间的关系，在给定了一部分参数的情况下，可用图线的形式对其性能与重要参数间的关系进行研究。
对于工作介质气体微团被束缚于各种一维势阱(如一维无限深势阱、包含相对论粒子的一维势阱、谐振势阱、四次势阱等)中的量子热声微循环，可统称为一维量子热声微循环(1D quantum thermoacoustic cycle, 1DQTAC)，即一维量子热声微循环包括工作于一维无限深势阱的量子热声微循环(1D infinite potential quantum thermoacoustic cycle, 1DIQTAC)、相对论粒子系统量子热声微循环(relativistic particles quantum thermoacoustic cycle, RQTAC)、谐振系统量子热声微循环(harmonic potential quantum thermoacoustic cycle, HQTAC)和四次势系统量子热声微循环(quartic potential quantum thermoacoustic cycle, QQTAC)。
由表1可见：一维无限深势阱系统与谐振系统具有相同的和，但与基态能量本征值
是有差别的，因此，当L4和T4给定时，这2种系统唯一的区别在于具有不同的
。又由式(7)可推得
，其中
(
)为图2中宏观状态4时一维无限深势阱系统(谐振系统)中激发态粒子的占有概率。换言之，当一维无限深势阱系统和谐振系统的p4e相同时，由式(26)和式(28)可知，这2类量子热声微循环的性能也相同。因此，下面的分析中合并考虑这2种情况。为了便于比较，在对各种不同势场系统进行分析时，采用以下条件：Tx=600 K，Tm/Tx=0.1。
由式(26)可绘出图3所示的粒子在特定量子力学势场中运动时的热声微循环效率
与p4e及之间的关系曲面。由图3可见：当p4e给定时，会随的增大呈单调递增的变化；而当取值较大时，3种不同系统的都会随p4e的取值变化而出现较明显的变化；但当不断减小时，p4e对的影响会变得非常微弱。在任意取值下，由
可求出ηmax和ηmin，此时对应的p4e分别为
和
。由图3可见：在相同参数条件下，相对论粒子系统中的ηmax和ηmin明显比谐振系统的低；而谐振系统中的ηmax和ηmin比四次势系统的略高，即粒子处于谐振势中时，热声发动机效率较高。
图3 
与p4e及
的关系
Fig. 3 Relationship among , p4e and
由式(28)可绘出图4所示的气体微团在特定势场中运动时的P*与p4e及之间的关系曲面。由图4可见：当p4e一定时，P*与呈类抛物线关系；而当一定时，P*与p4e亦呈类抛物线关系；并且P*-曲线及P*-p4e曲线对于3种势场条件都具有相似的形状；此外，存在最大值
以及相应的参数值
和
(可根据2个极值条件方程
，
联立得出)。在相同参数下，谐振系统的输出功率远比四次势系统的高；而四次势系统的输出功率又比相对论粒子系统的略高，即粒子处于谐振势中时热声发动机输出功率较高。此外，由于热声发动机的目的是对外输出有效功量，因此，若要保证
，则p4e的取值应满足
。
图4 P*与p4e及的关系
Fig. 4 Relationship among P*, p4e and
图5所示为p4e取不同数值时P*与之间的关系。由图5可知：在3种情况下，同一所对应的P*会先随p4e增加而增大，达到极限后继续增大p4e反而会使P*下降；在相同参数下，3种系统中循环的输出功率最大值满足
>
>
，其中下标H，R及Q分别代表谐振系统、相对论粒子系统和四次势系统；且3种系统中对应于的虽然很接近但依然满足
>
>
。可见，当粒子气体微团在谐振势中工作时，循环的输出功率与循环效率均比其他2种势场系统的高。
图5 P*与的关系
Fig. 5 Relationship between P* and
3 结论
1) 当p4e给定时，工作于一维无限深势阱的量子热声微循环与工作于谐振势的量子热声微循环的性能是相同的。而当L4和T4给定时，这2种系统中的p4e会因为不同而不同。
2) 当p4e给定时，3种势场条件下循环的都会随的增大呈单调递增趋势；而当给定时，3种势场条件下循环的都会随p4e的变化而分别达到最大值ηmax和最小值ηmin。通常，ηmax出现在p4e>0.5的区间中，而ηmin出现在p4e<0.5的区间。在相同参数下，处于3种不同势场的量子热声微循环满足
。
3) 当一定时，与呈类抛物线关系；而当一定时，与亦呈类抛物线关系；P*-曲线及P*-p4e曲线对于3种势场条件下的微循环都具有相似的形状。在相同参数下，处于3种不同势场的量子热声微循环满足
>
>
，并且P*达最大值时对应的p4e大于0.5。输出功率达最大值时循环效率也接近最大值，此时热机性能达最好。
4) 在相同参数下，当工质粒子工作于一维无限深势阱或谐振势场时，循环的效率和输出功率均比其他2种势场系统的高。
参考文献：
[1] 吴锋, 陈林根, 孙丰瑞. 斯特林机的有限时间热力学优化[M]. 北京: 化学工业出版社, 2008: 4-17.
WU Feng, CHEN Lingen, SUN Fengrui. Finite time thermodynamic optimization of Stirling engine[J]. Beijing: Chemical Industry Press, 2008: 4-17.
[2] CHEN Lingen, SUN Fengrui. Advances in finite time thermodynamics: analysis and optimization[M]. New York: Nova Science Publishers, 2004: 53-67.
[3] ANDRESEN B. Current trends in finite-time thermodynamics[J]. Angewandte Chemie International Edition, 2011, 50(12): 2690-2704.
[4] WANG Hao, WU Guoxing. Performance analysis and optimum criteria of a quantum dot engine with two discrete energy levels[J]. Physics Letters A, 2012, 376(33): 2209-2216.
[5] CHEN Lingen, LIU Xiaowei, GE Yanlin, et al. Ecological optimization of irreversible harmonic oscillator Carnot refrigerator [J]. Journal of the Energy Institute, 2016, 86(2): 85-96.
[6] MA Xiaomeng, LI Min, ZHOU Yueming, et al. Nonsequential double ionization of Xe by mid-infrared laser pulses[J]. Optical and Quantum Electronics, 2017, 49(4): 1-14.
[7] QIN Meiyan, ZHU Xiaosong. Molecular orbital imaging for partially aligned molecules[J]. Optics & Laser Technology, 2017, 87: 79-86.
[8] LIU Xiaowei, CHEN Lingen, GE Yanlin, et al. Fundamental optimal relation of a generalised irreversible quantum Carnot heat pump with harmonic oscillators[J]. International Journal of Ambient Energy, 2012, 33(3): 118-129.
[9] WU Feng, YANG Zhichun, CHEN Lingen, et al. Work output and efficiency of a reversible quantum Otto cycle[J]. Thermal Science, 2010, 14(4): 879-886.
[10] WU Feng, CHEN Lingen, SUN Fengrui, et al. Finite-time exergoeconomic performance bound for a quantum Stirling engine[J]. International Journal of Engineering Science, 2000, 38(2): 239-247.
[11] 吴锋, 汪拓, 陈林根, 等. 量子斯特林热机的输出功和热效率[J]. 机械工程学报, 2014, 50(4): 150-154.
WU Feng, WANG Tuo, CHEN Lingen, et al. Work output and efficiency of quantum stirling engine[J]. Journal of Mechanical Engineering, 2014, 50(4): 150-154.
[12] ACIKKALP E, CANER N. Application of exergetic sustainable index to the quantum irreversible diesel refrigerator cycles for 1D box system[J]. The European Physical Journal Plus, 2015, 130(4): 73-1-8.
[13] 鄂青, 吴锋. 广义一维势中热声制冷微循环的性能分析[J]. 中南大学学报(自然科学版), 2019, 50(3): 726-733.
E Qing, WU Feng. Performance analysis for quantum thermoacoustic refrigeration micro-cycle working in generalized 1D potential[J]. Journal of Central South University(Science and Technology), 2019, 50(3): 726-733.
[14] WANG Jianhui, HE Jizhou, MAO Zhiyuan. Performance of a quantum heat engine cycle working with harmonic oscillator systems[J]. Science in China Series G: Physics, Mechanics and Astronomy, 2007, 50(2): 163-176.
[15] WANG Jianhui, XIONG Shuangquan, HE Jizhou, et al. Performance analysis of a quantum heat engine working with a particle in a one-dimensional harmonic trap[J]. Acta Physica Sinica, 2012, 61(8): 080509.
[16] E Qing WU Feng, YIN Yong, et al. Optimal power and efficiency of quantum thermoacoustic micro-cycle working in 1D harmonic trap[J]. Journal of Low Temperature Physics, 2017, 189(1/2): 84-97.
[17] E Qing, WU Feng, CHEN Lingen, et al. Thermodynamic optimization for a quantum thermoacoustic refrigeration micro-cycle[J]. Journal of Central South University, 2020, 27(9): 2754-2762.
[18] WANG Jianhui, WU Zhaoqi, HE Jizhou. Quantum Otto engine of a two-level atom with single-mode fields[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 2012, 85(4): 041148.
[19] WANG Jianhui, MA Yongli, HE Jizhou. Quantum-mechanical engines working with an ideal gas with a finite number of particles confined inapower-law trap[J]. Euro Physics Letters, 2015, 111(20): 20006.
[20] 阚绪献, 吴锋, 张晓青, 等. 热声热机微热力学循环的最优性能[J]. 武汉理工大学学报, 2009, 31(14): 130-133.
KAN Xuxian, WU Feng, ZHANG Xiaoqing, et al. Performance optimization of a thermoacoustic engine micro-cycle[J]. Journal of Wuhan University of Technology, 2009, 31(14): 130-133.
[21] ANDRESEN B, BERRY R S, ONDRECHEN M J, et al. Thermodynamics for processes in finite time[J]. Accounts of Chemical Research, 1984, 17(8): 266-271.
(编辑 刘锦伟)
收稿日期： 2020 -11 -21； 修回日期： 2021 -01 -05
基金项目(Foundation item)：国家留学基金资助项目(202008420030)；武汉工程大学科学研究基金资助项目(K201919) (Project(202008420030) supported by the China Scholarship Council; Project(K201919) supported by the Scientific Research Foundation of Wuhan Institute of Technology)
通信作者：鄂青，博士，副教授，从事热声热机系统研究；E-mail: eqinghust@163.com
DOI: 10.11817/j.issn.1672-7207.2021.06.018
引用格式：鄂青, 吴锋, 雷霆. 广义一维势热声发动机微循环的性能分析[J].中南大学学报(自然科学版), 2021, 52(6): 1883-1891.
Citation:E Qing, WU Feng, LEI Ting. Performance analysis of micro-cycle in quantum thermoacoustic engine with generalized 1D potential[J]. Journal of Central South University(Science and Technology), 2021, 52(6): 1883-1891.
