CN103983365A - 多测头瞬态辐射热流计及热辐射热流密度的测定方法 - Google Patents

多测头瞬态辐射热流计及热辐射热流密度的测定方法 Download PDF

Info

Publication number
CN103983365A
CN103983365A CN201410234217.9A CN201410234217A CN103983365A CN 103983365 A CN103983365 A CN 103983365A CN 201410234217 A CN201410234217 A CN 201410234217A CN 103983365 A CN103983365 A CN 103983365A
Authority
CN
China
Prior art keywords
thermo
gauge head
responsive
heat
responsive gauge
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410234217.9A
Other languages
English (en)
Other versions
CN103983365B (zh
Inventor
黄勇
施国栋
朱克勇
王浚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201410234217.9A priority Critical patent/CN103983365B/zh
Publication of CN103983365A publication Critical patent/CN103983365A/zh
Application granted granted Critical
Publication of CN103983365B publication Critical patent/CN103983365B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radiation Pyrometers (AREA)

Abstract

本发明公开了一种快速响应非接触测量热流密度的热流计,包括多敏感面测头、热沉体、热沉体测温热电偶、数据采集装置、数据处理与显示装置,根据由于不同吸收率的两个相同材料,在受到相同强度辐射热流,它们的温度不一样,建立起辐射热流密度与两个敏感面之间温度差关系,同时还避免了***误差,开发一种既可满足空间技术真空环境下使用,也可以在地面对流环境中使用,结构更简单,可持续工作,并可以串联形成热电堆,以增大温差热电势,减少测量中误差,的高精度瞬态辐射热流计。

Description

多测头瞬态辐射热流计及热辐射热流密度的测定方法
技术领域
本发明涉及一种用于测定蒸汽锅炉、均热炉、输热管道、太空热辐射等各种热辐射热流密度的多测头瞬态辐射热流计及利用该热流计进行热辐射热流密度的测定方法。
背景技术
辐射热流计在太阳能、空间技术、气象、工业、冶金、能源、动力、空调等很多领域中都有重要的应用,热流检测理论和技术越来越受到人们的重视。目前使用的全辐射热流计分为热传导型稳态热流计和集总热容式热流计,都是依据热平衡原理设计的。为了提高对变动辐射流的测试精确度,需要高精度瞬态辐射热流计的研究。
主要几种常见瞬态辐射热流计都是由单个敏感面构成。第一种,通过建立起敏感面温度变化与辐射热流量关系来求得辐射热流,但是敏感面温度的变化并不是由敏感面本身直接测量得到,还是要利用其它热电偶来进行,这样热电偶与铜片之间存在着一个导热的过程,为了补偿在辐射热流突变的时候热电偶测得的温度与铜片真实温度因延时产生的动态误差,数据处理通常还要进行了补偿,会影响测量的误差和响应时间。第二种,圆箔式瞬态热流计,通过单个康铜敏感面与铜热沉体接触传递热量,建立热流与敏感面温度梯度关系,但是接触面积较大,热沉体的温度上升会影响到测量结果,不利于长时间工作。第三种,绝热型瞬态热流计,设计了跟踪敏感面温度变化的结构,使敏感面与跟踪面保持相同温度,建立辐射热流量与跟踪加热量的关系,这种方法利于长时间测量,但是结构较为复杂,受干扰影响的因素也较多。
发明内容
针对上述问题,本发明依据热平衡原理设计,结合圆箔式热流计和绝热型瞬态热流计全部优点,巧妙设计了瞬态辐射热流计的模型,通过双测头差分热电偶直接测量到温度,建立起辐射热流密度与两个敏感面之间温度差的关系,同时还避免了***误差,开发一种可满足空间技术真空环境下使用,也可以在地面对流环境中使用的高精度瞬态辐射热流计。
本发明的技术方案为:
一种多敏感面测头的辐射热流计,包括第一热敏感面测头、第二热敏感面测头、热沉体、热沉体测温热电偶、数据采集仪、控制***,与数据处理与显示装置部分,
所述的第一热敏感面测头后设有第一绝热材料,所述的第二热敏感面测头后设有第二绝热材料,第一热敏感面测头和第二热敏感面测头的热吸收率不同,第一热敏感面测头和第二热敏感面测头之间相互绝热,彼此从热源所接收到的辐射热流密在相同条件下是相同的。
所述的热沉体与热沉体测温热电偶相连接,所述第一热敏感面测头、第二热敏感面测头分别与热沉体连接,第一热敏感面测头、第二热敏感面测头、热沉体热电偶分别与数据采集仪连接,并将第一热敏感面测头、第二热敏感面测头,热沉体热电偶接收到的热信号传送到数据采集仪中,数据采集仪产生的模拟信号经放大、A/D转换处理转换成数字信号送入控制***,通过控制***控制数据采集仪,将温度信号传入数据处理与显示装置部分。
优选的,第一热敏感面测头和第二热敏感面测头为两片尺寸材料一致、并经过不同的表面粗糙工艺处理和上色处理的薄康铜片。
优选的,不同的表面粗糙工艺处理和上色处理的薄康铜片,其中之一是经过表面涂深色油漆粗糙表面处理,另一片是经过表面涂浅色油漆光洁表面处理。
优选的,所述的控制***为带有控制程序的计算机,
优选的,所述的辐射热流计的热敏感面测头外不加石英玻璃罩,
优选的,所述的热沉体与热敏感面测头通过铜导线连接。
本发明相对于现有技术所取得的有益技术效果为:
1、本发明多测头瞬态辐射热流计,两个热敏感面与热沉体构成多组铜-康铜热电偶,可瞬间测得各敏感面之间温度,获得辐射热流密度;
2、消除了***的部分误差,测量结果更接近真实值;
3、本发明的多测头瞬态辐射热流计,既可以在真空中测量,也可以在大气中或复杂条件下测量;并且传统单测头模型在大气环境下为减少对流热波动,需要在敏感测头外加石英玻璃罩,由于玻璃罩会对部分辐射能量进行吸收,这样就增大了实验的误差。双测头设计由于环境相同,在计算时可以通过相减消去环境对流换热对测量的影响,因此在大气环境下可以不加石英玻璃罩,这样消除石英玻璃本身对一定波段热辐射的吸收。
4、本发明多测头瞬态辐射热流计,沉体与热敏感面测头通过铜导线连接,使得热沉体与敏感面距离较远,温度不易快速上升,能长期的稳定工作。
具体实施方式
下面将结合附图和具体实施例来对本发明做进一步的详细说明。
图1是本发明多测头瞬态辐射热流计的传感装置结构图;
图中:1为深色粗糙薄康铜片测头;2为浅色光洁薄康铜片2测头;3为铜热沉体;a、b、c均为连接铜导线。
本发明的多测头瞬态辐射热流计,包括两片尺寸材料一致、质量均为m,面积均为A的薄康铜片,深色粗糙薄康铜片1表面涂深色油漆粗糙表面处理,温度为T1,浅色光洁薄康铜片2表面涂浅色油漆光洁表面处理,温度为T2,两片薄康铜片分别作为热流计的热敏感面测头,薄康铜片背后设有绝热材料,两个薄康铜片之间没有热交换。它们彼此从热源处所接收到的辐射热流密在相同条件下是相同的。以a1和a2分别表示深色粗糙薄康铜片1和浅色光洁薄康铜片2的吸收率。
如图1所示,深色粗糙薄康铜片1、浅色光洁薄康铜片2和铜热沉体3组成3个热电偶,铜热沉体温度为T0,由埋入热沉体的热电偶测得,深色粗糙薄康铜片1的温度T1由ab两点电压得到,浅色光洁薄康铜片2的温度T2由bc两点电压得到,深色粗糙薄康铜片1和浅色光洁薄康铜片2之间的温度差T1-T2由ac两点电压得到,这样热敏感面受到的热流发生变化时,将直接测量到两个热敏感面的温度变化。
将敏感面理想化为灰体,光谱吸收率与投入辐射无关,热敏感面受到来自热源的辐射热流密度为q的热辐射后,两个薄康铜片测头的能量平衡方程分别为:
a 1 qA = mc d T 1 dt + q L 1 - - - ( 1 )
a 2 qA = mc d T 2 dt + q L 2 - - - ( 2 )
上式中,a1为深色粗糙薄康铜片1的吸收率,a2为浅色光洁薄康铜片2的吸收率,q为辐射热流密度,A为薄康铜片面积,m为薄康铜片质量,c为薄康铜片比热容,qL1为深色粗糙薄康铜片1的换热损失,qL2为浅色光洁薄康铜片2的换热损失。
qL1和qL2包括三部分的换热损失。分别为热敏感面测头发出辐射损失,薄康铜片与周围温度为T的空气对流换热损失,薄康铜片与温度为Tb的背面隔热层底面的传导换热损失,该三部分的热损失计算方法如下:
1)在辐射背景温度很低情况下,深色粗糙薄康铜片1和浅色光洁薄康铜片2发出辐射损失qL11和qL12分别为:
qL11=ε1σAT1 4
q L11 =ε2σAT2 4
式中σ为斯蒂芬--波尔兹曼常数,σ=5.670×10-8W/(m2·K4),ε1、ε2分别为深色粗糙薄康铜片1和浅色光洁薄康铜片2发射率。物体吸收能量与向外辐射热量发射电磁波两者是方向相反的两个过程,物体吸收热量是落在物体上的辐射能,影响和改变了电子的激发,使得一部分的辐射能转变为热能而被物体所吸收。物体发射热量是发射电磁波辐射式由于电子的激发,使热能转变为辐射能发射出去。当物体被看为灰体时候,根据基尔霍夫定律,物体辐射能力与吸收能力可以从两个表面之间的辐射换热的研究得到:
ϵ = E E 0 = a
即任何可以看为灰体物体的发射率与吸收率在数值上是相等的,这说明了任何物体其吸收能力越强,其辐射能力也越强,即善于吸收者必善于辐射。对于一定波长下的黑度和吸收率也存在上述的关系。
2)深色粗糙薄康铜片1和浅色光洁薄康铜片2与周围温度为T的空气的对流换热qL21和qL22分别为:
qL21=hA(T1-T)
qL22=hA(T2-T)
上式中,h为对流换热系数,预先进行标定,两个测头位置比较接近,上方空气是连通的,这里认为T是一致。
对流换热系数h的大致量级为:
空气自然对流5~25;
气体强制对流20~100;
水的自然对流200~1000;
水的强制对流1000~15000;
油类的强制对流50~1500;
水蒸气的冷凝5000~15000;
有机蒸汽的冷凝500~2000;
水的沸腾2500~25000。
3)深色粗糙薄康铜片1和浅色光洁薄康铜片2通过背面隔热层对温度为Tb的底面的传导换热qL31和qL32分别为:
qL31=λA(T1-Tb)
qL32=λA(T2-Tb)
上式中,λ为隔热层导热系数,由于背部材料性质条件相同,这里显然Tb是一致。因此得到深色和浅色康铜片的总热损失为:
qL1=ε1σAT1 4+hA(T1-T)+λA(T1-Tb)  (3)
qL2=ε2σAT2 4+hA(T2-T)+λA(T2-Tb)  (4)
将(3)、(4)式分别代入(1)、(2)相减得到在大气中的瞬态辐射热流密度测量式为:
( a 1 - a 2 ) q = mc A d ( T 1 - T 2 ) dt + σ ( a 1 T 1 4 - a 2 T 2 4 ) + h ( T 1 - T 2 ) + λ ( T 1 - T 2 ) - - - ( 5 )
如果在真空条件中没有对流换热损失,λ为隔热层导热系数较小,当两个敏感面温差不是特别大情况,背部材料传热引起温差很小,空间的瞬态辐射热流密度式可进一步简化为:
( a 1 - a 2 ) q = mc A d ( T 1 - T 2 ) dt + σ ( ϵ 1 T 1 4 - ϵ 2 T 2 4 ) - - - ( 6 )
双测头深色粗糙薄康铜片1和浅色光洁薄康铜片2接收的热信号传到数据采集仪中,模拟信号经放大,A/D转换等处理转换成数字信号送入计算机,通过计算机软件编程控制数据采集仪,将温度信号传入并进行相应的数据处理与显示装置部分。

Claims (7)

1.一种多敏感面测头的辐射热流计,其特征在于,包括第一热敏感面测头、第二热敏感面测头、热沉体、热沉体测温热电偶、数据采集仪、控制***,与数据处理与显示装置部分;
所述的第一热敏感面测头后设有第一绝热材料,所述的第二热敏感面测头后设有第二绝热材料,第一热敏感面测头和第二热敏感面测头的热吸收率不同,第一热敏感面测头和第二热敏感面测头之间相互绝热,从热源所接收到的辐射热流密在相同条件下是相同的;
所述的热沉体与热沉体测温热电偶相连接,所述第一热敏感面测头、第二热敏感面测头分别与热沉体连接,第一热敏感面测头、第二热敏感面测头、热沉体热电偶分别与数据采集仪连接,并将第一热敏感面测头、第二热敏感面测头,热沉体热电偶接收到的热信号传送到数据采集仪中,数据采集仪产生的模拟信号经放大、A/D转换处理转换成数字信号送入控制***,通过控制***控制数据采集仪,将温度信号传入数据处理与显示装置部分。
2.如权利要求1所述的一种多敏感面测头的辐射热流计,其特征在于,第一热敏感面测头和第二热敏感面测头为两片尺寸材料一致、并经过不同的表面粗糙工艺处理和上色处理的薄康铜片。
3.如权利要求2所述的一种多敏感面测头的辐射热流计,其特征在于,所述的经过不同的表面粗糙工艺处理和上色处理的薄康铜片,其中之一是经过表面涂深色油漆粗糙表面处理,另一片是经过表面涂浅色油漆光洁表面处理。
4.如权利要求1所述的一种多敏感面测头的辐射热流计,其特征在于,所述的控制***为带有控制程序的计算机。
5.如权利要求1所述的一种多敏感面测头的辐射热流计,其特征在于,所述的辐射热流计的热敏感面测头外不加石英玻璃罩。
6.如权利要求1所述的一种多敏感面测头的辐射热流计,其特征在于,所述的热沉体与热敏感面测头通过铜导线连接。
7.一种利用权利要求1-6任一项所述的多敏感面测头的辐射热流计进行热辐射热流密度的测定方法,其特征在于,分别得到第一热敏感面测头的温度T1,第二热敏感面测头的温度T2,热沉体的温度T0,第一热敏感面测头和第二热敏感面测头之间的温度差T1-T2,热敏感面受到的热流发生变化时,将直接测量到两个热敏感面的温度变化;
将敏感面理想化为灰体,热敏感面受到来自热源的辐射热流密度为q的热辐射后,两个热敏感面测头的能量平衡方程分别为:
a 1 qA = mc d T 1 dt + q L 1 - - - ( 1 )
a 2 qA = mc d T 2 dt + q L 2 - - - ( 2 )
上式中,a1为第一热敏感面测头的吸收率,a2为第二热敏感面测头的吸收率,q为辐射热流密度,A为热敏感面测头面积,m为热敏感面测头质量,c为热敏感面测头比热容,qL1为第一热敏感面测头的换热损失,qL2为第二热敏感面测头的换热损失;
qL1和qL2包括三部分的换热损失;分别为热敏感面测头发出辐射损失,热敏感而测头与周围温度为T的空气对流换热损失,热敏感面测头与温度为Tb的背面隔热层底面的传导换热损失,该三部分的热损失计算方法如下:
1)在辐射背景温度很低情况下,第一热敏感面测头和第二热敏感面测头发出辐射损失qL11和qL11分别为:
qL11=ε1σAT1 4
qL12=ε2σAT2 4
2)式中σ为斯蒂芬一波尔兹曼常数,σ=5.670×10-8W/(m2·K4),物体视为灰体,ε1、ε2分别为第一热敏感面测头和第二热敏感面测头发射率;第一热敏感面测头和第二热敏感面测头与周围温度为T的空气的对流换热qL21和qL22分别为:
qL21=hA(T1-T)
qL22=hA(T2-L)
上式中,h为对流换热系数,预先进行标定,两个测头位置比较接近,上方空气是连通的,T是一致的;
3)第一热敏感面测头和第二热敏感面测头通过背面隔热层对温度为Tb的底面的传导换热qL31和qL32分别为:
qL31=λA(T1-Tb)
qL32=λA(T2-Tb)
上式中,λ为隔热层导热系数,由于背部材料性质条件相同,Tb是一致的;
得到第一热敏感面测头和第二热敏感面测头的总热损失分别为:
qL1=ε1σAT1 4+hA(T1-T)+λA(T1-Tb)                (3)
qL2=ε2σAT2 4+hA(T2-T)+λA(T2-Tb)                 (4)
将(3)、(4)式分别代入(1)、(2)相减得实现消去T和Tb得到在大气中的瞬态辐射热流密度测量式为:
( a 1 - a 2 ) q = mc A d ( T 1 - T 2 ) dt + σ ( a 1 T 1 4 - a 2 T 2 4 ) + h ( T 1 - T 2 ) + λ ( T 1 - T 2 ) - - - ( 5 )
在真空条件时,没有对流换热损失,λ为隔热层导热系数较小,当两个敏感而温差不是特别大情况,背部材料传热引起温差很小,空间的瞬态辐射热流密度式可进一步简化为:
( a 1 - a 2 ) q = mc A d ( T 1 - T 2 ) dt + σ ( ϵ 1 T 1 4 - ϵ 2 T 2 4 ) - - - ( 6 )
第一热敏感面测头和第二热敏感面测头接收的热信号传到数据采集仪中,模拟信号经放大,A/D转换成数字信号送入计算机,通过计算机软件编程控制数据采集仪,将温度信号传入并进行相应的数据处理与显示装置部分。
CN201410234217.9A 2014-05-29 2014-05-29 多测头瞬态辐射热流计及热辐射热流密度的测定方法 Expired - Fee Related CN103983365B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410234217.9A CN103983365B (zh) 2014-05-29 2014-05-29 多测头瞬态辐射热流计及热辐射热流密度的测定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410234217.9A CN103983365B (zh) 2014-05-29 2014-05-29 多测头瞬态辐射热流计及热辐射热流密度的测定方法

Publications (2)

Publication Number Publication Date
CN103983365A true CN103983365A (zh) 2014-08-13
CN103983365B CN103983365B (zh) 2017-02-15

Family

ID=51275425

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410234217.9A Expired - Fee Related CN103983365B (zh) 2014-05-29 2014-05-29 多测头瞬态辐射热流计及热辐射热流密度的测定方法

Country Status (1)

Country Link
CN (1) CN103983365B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106525288A (zh) * 2016-10-12 2017-03-22 北京遥测技术研究所 一种纯辐射热流传感器
CN109196318A (zh) * 2016-05-25 2019-01-11 株式会社电装 热流测定装置
CN110807280A (zh) * 2019-10-26 2020-02-18 复旦大学 一种同时适用于传导、对流和辐射的热聚集装置
CN110836725A (zh) * 2018-08-16 2020-02-25 哈尔滨工业大学 高热流耦合环境下双探头热流计及其热流密度的测定方法
US10608269B2 (en) 2015-09-04 2020-03-31 Ceres Intellectual Property Company Limited Temperature measuring method
CN117451217A (zh) * 2023-12-25 2024-01-26 中国空气动力研究与发展中心计算空气动力研究所 一种航天热流传感器及基于双温差补偿的热流修正方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101403638B (zh) * 2008-11-21 2011-08-17 北京航天计量测试技术研究所 带温度补偿的热流传感器
CN102175339B (zh) * 2011-02-18 2013-07-24 北京工业大学 燃料电池内部瞬态薄膜热流传感器
CN202393503U (zh) * 2011-12-01 2012-08-22 北京遥测技术研究所 圆箔热流传感器
CN103033291B (zh) * 2012-12-11 2015-11-25 北京遥测技术研究所 一种圆箔热电堆热流传感器

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10608269B2 (en) 2015-09-04 2020-03-31 Ceres Intellectual Property Company Limited Temperature measuring method
CN109196318A (zh) * 2016-05-25 2019-01-11 株式会社电装 热流测定装置
CN109196318B (zh) * 2016-05-25 2020-06-05 株式会社电装 热流测定装置
CN106525288A (zh) * 2016-10-12 2017-03-22 北京遥测技术研究所 一种纯辐射热流传感器
CN106525288B (zh) * 2016-10-12 2018-10-09 北京遥测技术研究所 一种纯辐射热流传感器
CN110836725A (zh) * 2018-08-16 2020-02-25 哈尔滨工业大学 高热流耦合环境下双探头热流计及其热流密度的测定方法
CN110807280A (zh) * 2019-10-26 2020-02-18 复旦大学 一种同时适用于传导、对流和辐射的热聚集装置
CN110807280B (zh) * 2019-10-26 2023-05-30 复旦大学 一种同时适用于传导、对流和辐射的热聚集装置
CN117451217A (zh) * 2023-12-25 2024-01-26 中国空气动力研究与发展中心计算空气动力研究所 一种航天热流传感器及基于双温差补偿的热流修正方法
CN117451217B (zh) * 2023-12-25 2024-03-12 中国空气动力研究与发展中心计算空气动力研究所 一种航天热流传感器及基于双温差补偿的热流修正方法

Also Published As

Publication number Publication date
CN103983365B (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN103983365A (zh) 多测头瞬态辐射热流计及热辐射热流密度的测定方法
Segrè et al. Rayleigh scattering in a liquid far from thermal equilibrium
CN106066208B (zh) 一种燃煤电站锅炉高温过热器壁面温度在线监测的装置及方法
Yu et al. Thermal conductivity, specific heat, thermal diffusivity, and emissivity of stored canola seeds with their temperature and moisture content
CN105628208B (zh) 一种基于红外成像***的温度测量方法
CN101419095B (zh) 灰体辐射率的测定方法
Manjhi et al. Performance assessment of K-type, E-type and J-type coaxial thermocouples on the solar light beam for short duration transient measurements
CN101144739A (zh) 一种高温材料黑体辐射发射率的测试原理及方法
CN102103085B (zh) 荧光光谱测量玻尔兹曼常数的方法
Diller Heat flux
CN103267773B (zh) 一种双环热保护瞬态辐射热流计及测量方法
CN105606640A (zh) 一种光电式金属线胀系数测定装置
CN101126731A (zh) 导温系数测量方法
Qi et al. Optimal design of E-type coaxial thermocouples for transient heat measurements in shock tunnels
Sapozhnikov et al. Bismuth-based gradient heat-flux sensors in thermal experiment
CN102095507A (zh) 利用对联热电偶测量内燃机热平衡中较小冷却液温差方法
Hubble et al. Development and evaluation of the time-resolved heat and temperature array
Latour et al. Local convective heat transfer identification by infrared thermography from a disc mounted on a cylinder in air crossflow
Katra et al. A method for estimating the spatial distribution of soil moisture of arid microenvironments by close range thermal infrared imaging
CN202631115U (zh) 热红外测温信号处理***
CN102944323A (zh) 一种基于真均方根转换器的温度脉动仪
Jia et al. The study of fluorescence thermal measurement based on DSP
CN105444908A (zh) 热电偶温度信号测量装置及测量方法
Goswami et al. Dynamic calibration of temperature sensors from light rays for transient measurement
Andò et al. A measurement tool for investigating cooling lava properties

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170215

Termination date: 20200529

CF01 Termination of patent right due to non-payment of annual fee