基于圆箔热流计的瞬态辐射热流测量方法
基于圆箔热流计的瞬态辐射热流测量方法
【技术领域】
[0001] 本发明属于瞬态辐射热流测量技术领域，特别涉及一种基于圆箔热流计的瞬态热 流测量方法。
【背景技术】
[0002] 热量传递是一种普遍存在的物理现象，作为一个基本的热物理量，热流密度的测 量在工业技术、健康环境等各领域受到广泛关注。圆箔热流计是一种应用于辐射热流测量 的传感器，具有结构稳定、测量量程大等优点，在燃烧室、火灾实验等领域中被广泛应用。圆 箔式热流计最早由RobertGardon发明，现在已成为一款发展成熟的商用热流传感器，美国 Medtherm、Vatell等公司都是圆箔式热流计的著名生产商。
[0003] 圆箔热流计（也称之为Gardon热流计）的选型主要根据厂家提供的传感器参数， 例如测量灵敏度、时间常数等。圆箔热流计的灵敏度在数值上等于稳态条件下的传感器输 出电动势与测量热流的比值，因此基于灵敏度参数，根据传感器的输出电动势，可以获得热 流密度，这是商用圆箔热流计测量输出的基本原理。圆箔式热流计的时间常数是描述传感 器测量响应时间的参数，即描述热流测量响应的延迟程度。
[0004] 尽管现有商用的圆箔热流计可以较好地应用于稳态辐射热流测量，但对于高瞬态 热流测量仍有一定的局限性：圆箔热流计对辐射热流的测量具有一定延迟，测量延迟使得 测量结果产生误差，尤其是在高瞬变热流的测量中，基于稳态条件下的热流计灵敏度参数 无法精确描述高瞬态热流；为解决瞬态测量问题，一般选用圆箔片半径较小的圆箔热流计， 然而通过减小圆箔片的半径提高测量响应速度的方法会受到机械加工技术的显著制约，热 流计高瞬态响应性能的提高对传感器机械微加工的技术要求变得更为苛刻。
【发明内容】
[0005] 基于上述问题，本发明提供一种基于圆箔热流计的瞬态辐射热流测量方法，在不 改变圆箔热流计硬件结构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出 电动势的测量，计算得到瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精 确测量。
[0006] 针对上述目的，本发明提供一种基于圆箔热流计的瞬态辐射热流测量方法，其特 征在于，包括：
[0007] 建立圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中 心的瞬态温度分布函数；
[0008] 建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函 数；
[0009] 测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和 所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到 所述圆箔热流计的瞬态辐射热流随时间的分布。
[0010] 其中，所述圆箔热流计的瞬态辐射热流测量数学模型通过所述圆箔热流计的圆箔 片热量传递守恒方程建立。
[0011] 其中，所述建立所述圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热 流计的圆箔片中心的瞬态温度分布函数，具体包括：
[0012] 建立所述圆箔热流计的瞬态辐射热流测量数学模型；
[0013] 求解所述数学模型，得到所述圆箔热流计的圆箔片中心的瞬态温度分布函数：
[0014] T(T) =F[T，q(T)，R，6,入，P，C，T0]
[0015] 其中，q(t)为t时刻的投射到所述圆箔片表面的瞬态辐射热流；R为所述圆箔片 的半径；S为所述圆箔片的厚度；X为所述圆箔片的导热系数；p为所述圆箔片的密度；c 为所述圆箔片的比热容；I；为热沉体温度。
[0016] 其中，所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函 数具体为：
[0017] E= 〇[T(t),T0]
[0018] 其中，E为所述圆箔热流计的瞬态输出电动势，T(t)为所述圆箔片中心的瞬态温 度，I；为所述圆箔热流计的热沉体温度。
[0019] 其中，所述基于所述圆箔片的瞬态温度分布函数和所述圆箔片中心的瞬态温度与 所述圆箔热流计的瞬态输出电动势的关联函数，测量所述圆箔热流计的瞬态输出电动势， 计算投射到所述圆箔热流计的瞬态辐射热流随时间的分布q(T)，具体包括：
[0020] 在测量的时间间隔AT内，采集所述圆箔热流计测量的m个瞬态输出电动势E;
[0021 ] 通过所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函 数，得到所述圆箔片中心的m个瞬态温度；
[0022] 基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布函 数，计算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0023] 其中，所述瞬态辐射热流随时间的分布用函数表示为：
[0025] 其中，i为关于时间t的级数序号（i= 1，. . .，n)，Q为待定系数，并且m彡n。
[0026] 其中，所述基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温 度分布函数，计算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布，具体包括：
[0027] 根据所述圆箔片中心的m个瞬态温度值，反演求解所述瞬态辐射热流随时间的分 布函数中的n个待定系数Q，即得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0028] 本发明的基于圆箔热流计的瞬态辐射热流测量方法，在不改变圆箔热流计硬件结 构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出电动势的测量，能够计 算获得该时间间隔内瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精确测 量，解决了商用圆箔热流计无法更好地应用于高瞬态辐射热流测量的技术局限性。另外，本 发明的测量方法简单可行，为圆箔热流计应用于高瞬态辐射热流的测量提供了重要的技术 指导，具有非常好的应用前景。
【附图说明】
[0029] 图1示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的流程图。
[0030] 图2示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的步骤S3的具体 流程图。
【具体实施方式】
[0031] 下面结合附图和实施例，对本发明的【具体实施方式】作进一步详细描述。以下实施 例用于说明本发明，但不用来限制本发明的范围。
[0032] 本发明的一个实施例提供了一种基于圆箔热流计的瞬态辐射热流测量方法。
[0033] 图1示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的流程图。
[0034] 参照图1，本发明的基于圆箔热流计的瞬态辐射热流测量方法，具体包括：
[0035] S1、建立所述圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的 圆箔片中心的瞬态温度分布函数；
[0036] 本实施例中，在入射瞬态辐射热流的作用下，圆箔热流计的圆箔片从中心向边缘 产生径向导热，圆箔片的中心与圆箔片边缘形成温差，圆箔片瞬态温度分布是外部瞬态辐 射热流、内部热传导共同作用的结果，满足热量传递守恒方程，此方程为瞬态辐射热流测量 的数学模型。
[0037] 通过求解圆箔热流计的圆箔片热量传递守恒方程，在t时刻，圆箔片中心的瞬态 温度分布T(t)与t、q(T)、R、S、A、p以及c相关，可以用函数F来表示：
[0038] T(t) =F[t,q(t),R, 8 ,A,p,c,T0]
[0039] 其中，q(t)为t时刻的投射到圆箔片表面的瞬态辐射热流；R为所述圆箔片的半 径；S为所述圆箔片的厚度；X为所述圆箔片的导热系数；p为所述圆箔片的密度；c为所 述圆箔片的比热容；I；为热沉体温度。
[0040] S2、建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联 函数；
[0041] 本实施例中，圆箔片的边缘与热流计的热沉体温度相同，I；是热沉体温度，为已知 值。圆箔片的中心与边缘具有温度梯度，构成了热电偶回路，因此，圆箔热流计瞬态输出电 动势E与圆箔片中心的瞬态温度T(T)、热沉体温度I；相关，热流计圆箔片中心的瞬态温度 与热流计瞬态输出电动势的关联函数为：
[0042]E= 〇[T(t),T0] 〇
[0043] 具体地，在一个实施例中，上述圆箔片中心的瞬态温度与热流计瞬态输出电动势 的关联函数可以为：
[0044] E=kAT(l+gAT)
[0045] 式中，AT二!'"）-!'#为热电偶电动势的输出系数；g为热电偶电动势的输出系 数的非线性修正项，并且在一个实施例中，k可以设置为：k= 0. 0387mV/K，g可以设置为：g =0? 0012K'
[0046] S3、测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函 数和所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投 射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0047] 图2示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的步骤S3的具体 流程图。
[0048] 如图2所示，步骤S3具体包括：
[0049]S31、在测量的时间间隔At内，采集所述圆箔热流计测量的m个瞬态输出电动势 E；
[0050] S32、通过所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关 联函数，得到所述圆箔片中心的m个瞬态温度；
[0051] S33、基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布 函数，计算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0052] 上述实施例中，在测量的任一时间间隔At内，瞬态辐射热流q(t)与时间t的 关系函数为：
[0054] 其中，Q为待定系数，在瞬态辐射热流级数展开表达式中共有n个待定系数（(^，… ，Cn)〇
[0055] 基于圆箔片中心的瞬态温度分布函数可知，热流计圆箔片中心的瞬态温度T(t) 与瞬态辐射热流q(t)相关，即T(t)与待定系数，…，Cn)相关。
[0056] 具体地，在一个实施例中，基于上述瞬态辐射热流与时间的关系函数，圆箔片中心 的的瞬态温度分布函数T(t)具体可以为：
[0058] 式中，T。为热沉体温度；a。为T。温度下圆箔片的热扩散率；a4 |a|是根据 所述圆箔片的测点中心温度确定的常数；JnSn阶bessel函数；0 "为
的第 n个解。
[0059] 在上述过程中，当m多n时，基于瞬态温度分布函数和瞬态辐射热流与时间的关系 函数，通过m个瞬态温度值，结合数学算法反演求解n个待定系数，…，Cn)，即确定了瞬 态辐射热流q(T)随时间的具体分布。
[0060] 本发明的基于圆箔热流计的瞬态辐射热流测量方法，在不改变圆箔热流计硬件结 构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出电动势的测量，能够计 算获得该时间间隔内瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精确测 量，解决了商用圆箔热流计无法更好地应用于高瞬态辐射热流测量的技术局限性。另外，本 申请的测量方法简单可行，为圆箔热流计应用于高瞬态辐射热流的测量提供了重要的技术 指导，具有非常好的应用前景。
[0061] 以上实施方式仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通 技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有 等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定。
【主权项】
1. 一种基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，包括： 建立圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中心的 瞬态温度分布函数； 建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数； 测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和所述 圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到所述 圆箔热流计的瞬态辐射热流随时间的分布。2. 如权利要求1所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 圆箔热流计的瞬态辐射热流测量数学模型通过所述圆箔热流计的圆箔片热量传递守恒方 程建立。3. 如权利要求2所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 建立所述圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中心的 瞬态温度分布函数，具体包括： 建立所述圆箔热流计的瞬态辐射热流测量数学模型； 求解所述数学模型，得到所述圆箔热流计的圆箔片中心的瞬态温度分布函数： Τ( τ )=F[ τ，q( τ )，R，δ，λ，p，c，T0] 其中，q(T)为τ时刻的投射到所述圆箔片表面的瞬态辐射热流；R为所述圆箔片的半 径；S为所述圆箔片的厚度；λ为所述圆箔片的导热系数；ρ为所述圆箔片的密度；c为所 述圆箔片的比热容；TO为热沉体温度。4. 如权利要求1所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数具体为： E = Φ[Τ(τ),Τ〇] 其中，E为所述圆箔热流计的瞬态输出电动势，Τ(τ)为所述圆箔片中心的瞬态温度，?； 为所述圆箔热流计的热沉体温度。5. 如权利要求1所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和所述圆箔 片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到所述圆箔 热流计的瞬态辐射热流随时间的分布q( τ )，具体包括： 在测量的时间间隔△ τ内，采集所述圆箔热流计测量的m个瞬态输出电动势E ; 通过所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，得 到所述圆箔片中心的m个瞬态温度； 基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布函数，计 算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。6. 如权利要求5所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 瞬态辐射热流随时间的分布用函数表示为：其中，i为关于时间τ的级数序号（i = 1，. . .，n)，Ci为待定系数，并且m彡η。7. 如权利要求6所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布函数，计算得 到投射到所述圆箔热流计的瞬态辐射热流随时间的分布，具体包括： 根据所述圆箔片中心的m个瞬态温度值，反演求解所述瞬态辐射热流随时间的分布函 数中的η个待定系数Ci,即得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
【专利摘要】本发明公开了一种基于圆箔热流计的瞬态热流测量方法，包括：建立圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中心的瞬态温度分布函数；建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数；测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到所述圆箔热流计的瞬态辐射热流随时间的分布。本发明在不改变圆箔热流计硬件结构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出电动势的测量，计算得到瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精确测量。
【IPC分类】G01K17/00, G01K7/04
【公开号】CN104897306
【申请号】CN201510266825
【发明人】符泰然, 宗安州, 庞传和
【申请人】清华大学
【公开日】2015年9月9日
【申请日】2015年5月22日
【技术领域】
[0001] 本发明属于瞬态辐射热流测量技术领域，特别涉及一种基于圆箔热流计的瞬态热 流测量方法。
【背景技术】
[0002] 热量传递是一种普遍存在的物理现象，作为一个基本的热物理量，热流密度的测 量在工业技术、健康环境等各领域受到广泛关注。圆箔热流计是一种应用于辐射热流测量 的传感器，具有结构稳定、测量量程大等优点，在燃烧室、火灾实验等领域中被广泛应用。圆 箔式热流计最早由RobertGardon发明，现在已成为一款发展成熟的商用热流传感器，美国 Medtherm、Vatell等公司都是圆箔式热流计的著名生产商。
[0003] 圆箔热流计（也称之为Gardon热流计）的选型主要根据厂家提供的传感器参数， 例如测量灵敏度、时间常数等。圆箔热流计的灵敏度在数值上等于稳态条件下的传感器输 出电动势与测量热流的比值，因此基于灵敏度参数，根据传感器的输出电动势，可以获得热 流密度，这是商用圆箔热流计测量输出的基本原理。圆箔式热流计的时间常数是描述传感 器测量响应时间的参数，即描述热流测量响应的延迟程度。
[0004] 尽管现有商用的圆箔热流计可以较好地应用于稳态辐射热流测量，但对于高瞬态 热流测量仍有一定的局限性：圆箔热流计对辐射热流的测量具有一定延迟，测量延迟使得 测量结果产生误差，尤其是在高瞬变热流的测量中，基于稳态条件下的热流计灵敏度参数 无法精确描述高瞬态热流；为解决瞬态测量问题，一般选用圆箔片半径较小的圆箔热流计， 然而通过减小圆箔片的半径提高测量响应速度的方法会受到机械加工技术的显著制约，热 流计高瞬态响应性能的提高对传感器机械微加工的技术要求变得更为苛刻。
【发明内容】
[0005] 基于上述问题，本发明提供一种基于圆箔热流计的瞬态辐射热流测量方法，在不 改变圆箔热流计硬件结构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出 电动势的测量，计算得到瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精 确测量。
[0006] 针对上述目的，本发明提供一种基于圆箔热流计的瞬态辐射热流测量方法，其特 征在于，包括：
[0007] 建立圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中 心的瞬态温度分布函数；
[0008] 建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函 数；
[0009] 测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和 所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到 所述圆箔热流计的瞬态辐射热流随时间的分布。
[0010] 其中，所述圆箔热流计的瞬态辐射热流测量数学模型通过所述圆箔热流计的圆箔 片热量传递守恒方程建立。
[0011] 其中，所述建立所述圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热 流计的圆箔片中心的瞬态温度分布函数，具体包括：
[0012] 建立所述圆箔热流计的瞬态辐射热流测量数学模型；
[0013] 求解所述数学模型，得到所述圆箔热流计的圆箔片中心的瞬态温度分布函数：
[0014] T(T) =F[T，q(T)，R，6,入，P，C，T0]
[0015] 其中，q(t)为t时刻的投射到所述圆箔片表面的瞬态辐射热流；R为所述圆箔片 的半径；S为所述圆箔片的厚度；X为所述圆箔片的导热系数；p为所述圆箔片的密度；c 为所述圆箔片的比热容；I；为热沉体温度。
[0016] 其中，所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函 数具体为：
[0017] E= 〇[T(t),T0]
[0018] 其中，E为所述圆箔热流计的瞬态输出电动势，T(t)为所述圆箔片中心的瞬态温 度，I；为所述圆箔热流计的热沉体温度。
[0019] 其中，所述基于所述圆箔片的瞬态温度分布函数和所述圆箔片中心的瞬态温度与 所述圆箔热流计的瞬态输出电动势的关联函数，测量所述圆箔热流计的瞬态输出电动势， 计算投射到所述圆箔热流计的瞬态辐射热流随时间的分布q(T)，具体包括：
[0020] 在测量的时间间隔AT内，采集所述圆箔热流计测量的m个瞬态输出电动势E;
[0021 ] 通过所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函 数，得到所述圆箔片中心的m个瞬态温度；
[0022] 基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布函 数，计算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0023] 其中，所述瞬态辐射热流随时间的分布用函数表示为：
[0025] 其中，i为关于时间t的级数序号（i= 1，. . .，n)，Q为待定系数，并且m彡n。
[0026] 其中，所述基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温 度分布函数，计算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布，具体包括：
[0027] 根据所述圆箔片中心的m个瞬态温度值，反演求解所述瞬态辐射热流随时间的分 布函数中的n个待定系数Q，即得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0028] 本发明的基于圆箔热流计的瞬态辐射热流测量方法，在不改变圆箔热流计硬件结 构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出电动势的测量，能够计 算获得该时间间隔内瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精确测 量，解决了商用圆箔热流计无法更好地应用于高瞬态辐射热流测量的技术局限性。另外，本 发明的测量方法简单可行，为圆箔热流计应用于高瞬态辐射热流的测量提供了重要的技术 指导，具有非常好的应用前景。
【附图说明】
[0029] 图1示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的流程图。
[0030] 图2示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的步骤S3的具体 流程图。
【具体实施方式】
[0031] 下面结合附图和实施例，对本发明的【具体实施方式】作进一步详细描述。以下实施 例用于说明本发明，但不用来限制本发明的范围。
[0032] 本发明的一个实施例提供了一种基于圆箔热流计的瞬态辐射热流测量方法。
[0033] 图1示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的流程图。
[0034] 参照图1，本发明的基于圆箔热流计的瞬态辐射热流测量方法，具体包括：
[0035] S1、建立所述圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的 圆箔片中心的瞬态温度分布函数；
[0036] 本实施例中，在入射瞬态辐射热流的作用下，圆箔热流计的圆箔片从中心向边缘 产生径向导热，圆箔片的中心与圆箔片边缘形成温差，圆箔片瞬态温度分布是外部瞬态辐 射热流、内部热传导共同作用的结果，满足热量传递守恒方程，此方程为瞬态辐射热流测量 的数学模型。
[0037] 通过求解圆箔热流计的圆箔片热量传递守恒方程，在t时刻，圆箔片中心的瞬态 温度分布T(t)与t、q(T)、R、S、A、p以及c相关，可以用函数F来表示：
[0038] T(t) =F[t,q(t),R, 8 ,A,p,c,T0]
[0039] 其中，q(t)为t时刻的投射到圆箔片表面的瞬态辐射热流；R为所述圆箔片的半 径；S为所述圆箔片的厚度；X为所述圆箔片的导热系数；p为所述圆箔片的密度；c为所 述圆箔片的比热容；I；为热沉体温度。
[0040] S2、建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联 函数；
[0041] 本实施例中，圆箔片的边缘与热流计的热沉体温度相同，I；是热沉体温度，为已知 值。圆箔片的中心与边缘具有温度梯度，构成了热电偶回路，因此，圆箔热流计瞬态输出电 动势E与圆箔片中心的瞬态温度T(T)、热沉体温度I；相关，热流计圆箔片中心的瞬态温度 与热流计瞬态输出电动势的关联函数为：
[0042]E= 〇[T(t),T0] 〇
[0043] 具体地，在一个实施例中，上述圆箔片中心的瞬态温度与热流计瞬态输出电动势 的关联函数可以为：
[0044] E=kAT(l+gAT)
[0045] 式中，AT二!'"）-!'#为热电偶电动势的输出系数；g为热电偶电动势的输出系 数的非线性修正项，并且在一个实施例中，k可以设置为：k= 0. 0387mV/K，g可以设置为：g =0? 0012K'
[0046] S3、测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函 数和所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投 射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0047] 图2示出了本发明的基于圆箔热流计的瞬态辐射热流测量方法的步骤S3的具体 流程图。
[0048] 如图2所示，步骤S3具体包括：
[0049]S31、在测量的时间间隔At内，采集所述圆箔热流计测量的m个瞬态输出电动势 E；
[0050] S32、通过所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关 联函数，得到所述圆箔片中心的m个瞬态温度；
[0051] S33、基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布 函数，计算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
[0052] 上述实施例中，在测量的任一时间间隔At内，瞬态辐射热流q(t)与时间t的 关系函数为：
[0054] 其中，Q为待定系数，在瞬态辐射热流级数展开表达式中共有n个待定系数（(^，… ，Cn)〇
[0055] 基于圆箔片中心的瞬态温度分布函数可知，热流计圆箔片中心的瞬态温度T(t) 与瞬态辐射热流q(t)相关，即T(t)与待定系数，…，Cn)相关。
[0056] 具体地，在一个实施例中，基于上述瞬态辐射热流与时间的关系函数，圆箔片中心 的的瞬态温度分布函数T(t)具体可以为：
[0058] 式中，T。为热沉体温度；a。为T。温度下圆箔片的热扩散率；a4 |a|是根据 所述圆箔片的测点中心温度确定的常数；JnSn阶bessel函数；0 "为
的第 n个解。
[0059] 在上述过程中，当m多n时，基于瞬态温度分布函数和瞬态辐射热流与时间的关系 函数，通过m个瞬态温度值，结合数学算法反演求解n个待定系数，…，Cn)，即确定了瞬 态辐射热流q(T)随时间的具体分布。
[0060] 本发明的基于圆箔热流计的瞬态辐射热流测量方法，在不改变圆箔热流计硬件结 构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出电动势的测量，能够计 算获得该时间间隔内瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精确测 量，解决了商用圆箔热流计无法更好地应用于高瞬态辐射热流测量的技术局限性。另外，本 申请的测量方法简单可行，为圆箔热流计应用于高瞬态辐射热流的测量提供了重要的技术 指导，具有非常好的应用前景。
[0061] 以上实施方式仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通 技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有 等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定。
【主权项】
1. 一种基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，包括： 建立圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中心的 瞬态温度分布函数； 建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数； 测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和所述 圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到所述 圆箔热流计的瞬态辐射热流随时间的分布。2. 如权利要求1所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 圆箔热流计的瞬态辐射热流测量数学模型通过所述圆箔热流计的圆箔片热量传递守恒方 程建立。3. 如权利要求2所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 建立所述圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中心的 瞬态温度分布函数，具体包括： 建立所述圆箔热流计的瞬态辐射热流测量数学模型； 求解所述数学模型，得到所述圆箔热流计的圆箔片中心的瞬态温度分布函数： Τ( τ )=F[ τ，q( τ )，R，δ，λ，p，c，T0] 其中，q(T)为τ时刻的投射到所述圆箔片表面的瞬态辐射热流；R为所述圆箔片的半 径；S为所述圆箔片的厚度；λ为所述圆箔片的导热系数；ρ为所述圆箔片的密度；c为所 述圆箔片的比热容；TO为热沉体温度。4. 如权利要求1所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数具体为： E = Φ[Τ(τ),Τ〇] 其中，E为所述圆箔热流计的瞬态输出电动势，Τ(τ)为所述圆箔片中心的瞬态温度，?； 为所述圆箔热流计的热沉体温度。5. 如权利要求1所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和所述圆箔 片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到所述圆箔 热流计的瞬态辐射热流随时间的分布q( τ )，具体包括： 在测量的时间间隔△ τ内，采集所述圆箔热流计测量的m个瞬态输出电动势E ; 通过所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，得 到所述圆箔片中心的m个瞬态温度； 基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布函数，计 算得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。6. 如权利要求5所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 瞬态辐射热流随时间的分布用函数表示为：其中，i为关于时间τ的级数序号（i = 1，. . .，n)，Ci为待定系数，并且m彡η。7. 如权利要求6所述的基于圆箔热流计的瞬态辐射热流测量方法，其特征在于，所述 基于上述圆箔片中心的m个瞬态温度，根据所述圆箔片中心的瞬态温度分布函数，计算得 到投射到所述圆箔热流计的瞬态辐射热流随时间的分布，具体包括： 根据所述圆箔片中心的m个瞬态温度值，反演求解所述瞬态辐射热流随时间的分布函 数中的η个待定系数Ci,即得到投射到所述圆箔热流计的瞬态辐射热流随时间的分布。
【专利摘要】本发明公开了一种基于圆箔热流计的瞬态热流测量方法，包括：建立圆箔热流计的瞬态辐射热流测量数学模型，获得所述圆箔热流计的圆箔片中心的瞬态温度分布函数；建立所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数；测量所述圆箔热流计的瞬态输出电动势，基于所述圆箔片的瞬态温度分布函数和所述圆箔片中心的瞬态温度与所述圆箔热流计的瞬态输出电动势的关联函数，计算投射到所述圆箔热流计的瞬态辐射热流随时间的分布。本发明在不改变圆箔热流计硬件结构和使用条件的前提下，通过在任一个时间间隔内的多个瞬态输出电动势的测量，计算得到瞬态辐射热流随时间变化的分布，从而实现高瞬态辐射热流的精确测量。
【IPC分类】G01K17/00, G01K7/04
【公开号】CN104897306
【申请号】CN201510266825
【发明人】符泰然, 宗安州, 庞传和
【申请人】清华大学
【公开日】2015年9月9日
【申请日】2015年5月22日
最新回复(0)