CN107976384A - 单色激光诱导白炽光纳米级碳烟粒径测量装置及方法 - Google Patents
单色激光诱导白炽光纳米级碳烟粒径测量装置及方法 Download PDFInfo
- Publication number
- CN107976384A CN107976384A CN201710984620.7A CN201710984620A CN107976384A CN 107976384 A CN107976384 A CN 107976384A CN 201710984620 A CN201710984620 A CN 201710984620A CN 107976384 A CN107976384 A CN 107976384A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- flame
- incandescence
- 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
Links
- 239000004071 soot Substances 0.000 title claims abstract description 140
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 94
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000002245 particle Substances 0.000 claims abstract description 155
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000004590 computer program Methods 0.000 claims abstract description 8
- 230000005055 memory storage Effects 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 20
- 230000002776 aggregation Effects 0.000 claims description 11
- 238000004220 aggregation Methods 0.000 claims description 11
- 238000001816 cooling Methods 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 10
- 235000019504 cigarettes Nutrition 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000007423 decrease Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000003750 conditioning effect Effects 0.000 claims description 5
- 230000003287 optical effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- 239000000779 smoke Substances 0.000 claims description 2
- 238000002485 combustion reaction Methods 0.000 claims 1
- 238000004891 communication Methods 0.000 claims 1
- 230000006698 induction Effects 0.000 abstract description 14
- 238000005516 engineering process Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 238000004821 distillation Methods 0.000 description 2
- 230000005622 photoelectricity Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
- G01N15/0205—Investigating particle size or size distribution by optical means
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明公开了一种单色激光诱导白炽光纳米级碳烟粒径测量装置及方法,属于碳烟颗粒的粒径测量技术,测量装置包括火焰发生单元、对火焰的碳烟颗粒进行加热的激光器、用于捕获火焰中白炽光信号的捕获单元、显示并记录所捕获的白炽光信号与时间的关系的示波器以及与示波器通信连接的处理单元;处理单元包括存储器和处理器,存储器存储有计算机程序,该计算机程序被处理器执行时能实现以下步骤:定义碳烟颗粒的激光诱导炽光信号强度降低到其峰值的1/e所需时间为信号特征时间,基于包含碳烟粒径与信号特征时间关系的粒径计算模型得出火焰中碳烟颗粒的粒径,并用标定物性参数的方法获得更加准确的粒径值,从而方便快捷地获得碳烟的二维粒径分布。
Description
技术领域
本发明涉及碳烟颗粒的粒径测量技术,具体地说,涉及一种单色激光诱导白炽光纳米级碳烟粒径测量装置及方法。
背景技术
近些年来,大气中微粒的危害受到人们的广泛关注,而碳烟对人体的危害大小与粒径有很大的关系,碳烟粒径越小,进入人体肺部停滞的比例越大,对人体危害也就也越大,因此研究纳米级碳烟的粒径测量更有实际意义。而激光诱导白炽光技术具有高时间和空间分辨率,对碳烟测试是一种很有发展前景的技术。
目前的激光诱导白炽光技术通常为双色激光诱导技术,需要建立双色激光诱导白炽光***,基于激光加热颗粒的温度衰减公式准确计算碳烟粒径。然而,在实验条件下双色激光诱导白炽光***需要同时检测和严格校准,而且因为大颗粒的冷却速度比小颗粒慢,多分散颗粒的测量温度通常只能视为颗粒温度的有效值,不能准确反映每个颗粒的温度,所以双色激光诱导白炽光操作困难且准确率低。
发明内容
本发明的目的为提供一种单色激光诱导白炽光纳米级碳烟粒径测量装置以及单色激光诱导白炽光纳米级碳烟粒径测量方法。该测量装置和测量方法准确率高,操作简单,易于推广。
为了实现上述目的,本发明提供的单色激光诱导白炽光纳米级碳烟粒径测量装置包括火焰发生单元、对火焰的碳烟颗粒进行加热的激光器、用于捕获火焰中白炽光信号的捕获单元、显示并记录所捕获的白炽光信号与时间的关系的示波器以及与示波器通信连接的处理单元;捕获单元包括放置在火焰一侧的光电倍增管以及放置在光电倍增管前用于阻止火焰中除白炽光信号之外的光信号通过的窄带通滤光片;处理单元包括存储器和处理器,存储器存储有计算机程序,该计算机程序被处理器执行时能实现以下步骤:
建立火焰中碳烟颗粒的白炽光信号强度与时间的关系,定义白炽光信号强度降低到其峰值的1/e时所需时间为信号特征时间τλ,基于包含碳烟粒径与信号特征时间τλ关系的粒径计算模型得出火焰中碳烟颗粒的粒径。
上述测量装置中,通过获取碳烟颗粒的信号特征时间,并通过粒径计算模型得出火焰中碳烟颗粒的粒径,与现有技术相比,能有效简化测量过程,提高了测量效率。
具体的方案为纳米级碳烟颗粒粒径计算模型是
其中,C是由聚集体的屏蔽效应参数以及碳烟颗粒与周围气体的热力学参数组成的系数,Tg碳烟颗粒周围气体的温度,Th-max表示碳烟颗粒冷却过程中的最高温度,λ是光电倍增管的设定探测波长,h表示普朗克常数,c表示光速,k表示玻尔兹曼常数。
C的表达式是
其中,αT表示热调节系数,pg表示碳烟颗粒周围气体的压力,Np表示聚集体中基元体颗粒的数量,ρs是碳烟的密度,cs是碳烟的比热,R是通用气体常数,Wg是碳烟颗粒周围气体的平均摩尔质量,γ是周围气体的比热比,kh和Dh均为碳烟颗粒的缩放系数。
作为优选,测量装置还包括对碳烟颗粒进行TEM采样的透镜电镜,粒径计算模型为碳烟粒径与特征时间的二元关系,该二元关系由以下步骤得到:
对于同一工况下的碳烟颗粒,首先在2800-3700K的范围内设定若干个Th-max的值得到不同位置处的若干个信号特征时间τλ,采用透射电镜TEM采样观察并得到这些位置的碳烟粒径Dp,基于实验所设定的探测波长和火焰温度,计算选择最优的Th-max,然后在二维坐标系中标定由τλ和Dp确定的点,并拟合得到碳烟粒径Dp与信号特征时间τλ的二元关系式,最后根据该关系式由不同测量位置处的信号特征时间τλ得到火焰中Dp的二维分布结果;对于不同工况下的碳烟颗粒,按照上述步骤得到不同的二元关系并由各工况下信号特征时间τλ得到Dp的值。
作为优选,激光器为脉冲激光器,其所发射的激光的波长为1064nm,作用在碳烟上的激光能量密度不超过0.1J/cm2。
为减小火焰背景光的影响,光电倍增管的响应波段和窄带滤光片的透过波长应位于可见光范围。为保持信号测量的良好跟随性,光电倍增管的上升时间应小于1ns。
作为优选,激光器与火焰之间设有用于对正火焰中心的光阑以及用于调节激光器发射脉冲能量的衰减器;激光器放置在一移动平台上,该移动平台根据光阑的对正信号校正激光器的位置以使其对准火焰中心。
本发明提供的单色激光诱导白炽光纳米级碳烟粒径测量方法包括:
获取信号特征时间步骤,采用脉冲激光对火焰中碳烟颗粒进行加热,获取火焰中碳烟颗粒的白炽光信号强度降低到其峰值的1/e时所需时间,即信号特征时间;
确定系数C步骤,C的表达式为
其中,αT表示热调节系数,pg表示碳烟颗粒周围气体的压力,Np表示聚集体中基元体颗粒的数量,ρs是碳烟的密度,cs是碳烟的比热,R是通用气体常数,Wg是碳烟颗粒周围气体的平均摩尔质量,γ是周围气体的比热比,kh和Dh均为碳烟颗粒的缩放系数;
计算步骤,基于纳米级碳烟粒径计算模型,依据获取的信号特征时间及确定的系数C计算碳烟粒径,所述粒径计算模型为
其中,C是由聚集体的屏蔽效应参数以及碳烟颗粒与周围气体的热力学参数组成的系数,Tg碳烟颗粒周围气体的温度,Th-max表示碳烟颗粒冷却过程中的最高温度,λ是光电倍增管的设定探测波长,h表示普朗克常数,c表示光速,k表示玻尔兹曼常数。
本发明提供的另一种单色激光诱导白炽光纳米级碳烟粒径测量方法包括:
建模步骤,建立碳烟粒径Dp与信号特征时间τλ的关系模型,定义火焰中碳烟颗粒的白炽光信号强度降低到其原始值的1/e时所需时间为信号特征时间τλ;
获取步骤,测量并计算火焰中激光诱导白炽光信号的信号特征时间τλ,并结合透射电镜采样观测得到的碳烟粒径结果,获取碳烟粒径Dp与信号特征时间τλ之间的关系系数C;
推演步骤,基于所述碳烟粒径Dp与信号特征时间τλ的二元关系系数C,根据火焰不同位置处的激光诱导白炽光信号的信号特征时间τλ,得出火焰中碳烟粒径的二维分布图,对应得出碳烟粒径。
优选地,碳烟颗粒最高温度设定为2800-3700K。
与现有技术相比,本发明的有益效果为:
本发明提供了一种新型的可行性方法测量碳烟粒径,避免了大量计算拟合中的费时费力和双色激光诱导白炽光***中可操作性差等问题,在热传导为主导的冷却过程中准确测量火焰中碳烟的二维粒径分布,相对误差小。
附图说明
图1为本发明实施例的单色激光诱导白炽光纳米级碳烟粒径测量装置示意图;
图2为本发明实施例1的存储器中计算机程序被处理器执行时的步骤图;
图3为本发明实施例在不同激光能量密度下Dp-λ和Dp-T大小及其一致性关系图;
图4为本发明实施例2的存储器中计算机程序被处理器执行时的步骤图;
图5为本发明实施例的激光器发射能量密度为0.09J/cm2时粒径Dp-λ和Dp-T与火焰高度HAB之间的关系曲线图;
图6为本发明实施例的Tg=1600K、λ=500nm时由颗粒实际最高温度Th-max的变化导致的粒径误差曲线;
图7为本发明实施例的碳烟颗粒在放大倍数20万倍后经透射电镜测量的结果图。
具体实施方式
以下结合实施例及其附图对本发明作进一步说明。
实施例1
现有技术中,在典型的Santoro火焰上进行双色激光诱导白炽光实验用如下公式计算粒径
其中,τT表示温度衰减时间,Tg表示周围气体的温度,聚集体的屏蔽效应参数和颗粒及周围气体的热力学参数组成系数C,即
其中,αT表示热调节系数,pg表示碳烟颗粒周围气体的压力,Np表示聚集体中基元体颗粒的数量,ρs是碳烟的密度,cs是碳烟的比热,R是通用气体常数,Wg是碳烟颗粒周围气体的平均摩尔质量,γ是周围气体的比热比,kh和Dh均为碳烟颗粒的缩放系数。
在实验条件下双色激光诱导白炽光***需要同时检测和严格校准,而且因为大颗粒的冷却速度比小颗粒慢,多分散颗粒的测量温度通常只能视为颗粒温度的有效值,不能准确反映每个颗粒的温度,所以双色激光诱导白炽光操作困难且准确率低。
本实施例采用单色激光诱导白炽光进行纳米级碳烟粒径测量,参见图1,其测量装置包括火焰发生单元4、激光器1、捕获单元、示波器8以及处理单元。
火焰发生单元4用于产生火焰,激光器1的激光信号发射到火焰中心对火焰中心的碳烟颗粒进行加热,为了使激光器1的激光信号能够准确地发射到火焰中心,在火焰发生装置4与激光器1之间安装有光阑3,且将激光器1安装在移动平台上,光阑3用于对正火焰中心,移动平台根据光阑的对正信号校正激光器的位置以使其对准火焰中心。在火焰发生装置4与激光器1之间还设有衰减器2,用于调节激光器1发射的脉冲能量。激光器1穿过火焰中心后被光束截止器5截止。
捕获单元包括放置在火焰一侧的光电倍增管7以及放置在光电倍增管7前用于阻止火焰中除白炽光信号之外的光信号通过的且以405nm和600nm为中心波长的窄带通滤光片6(FWHM=10nm)。示波器8接收由光电倍增管7捕获的白炽光信号,显示并记录该信号与时间的关系。同时,示波器8与处理单元通信连接。
处理单元包括存储器和处理器,存储器存储有计算机程序,该计算机程序被处理器执行时能实现图2所示步骤:
获取信号特征时间步骤S1,采用脉冲激光对火焰中碳烟颗粒进行加热,获取火焰中碳烟颗粒的白炽光信号强度降低到其峰值的1/e时所需时间,即为信号特征时间τλ;
确定系数C步骤S2,C的表达式为上述式(2)
根据参考文献将各物性参数具体值代入C中计算,各参数取值如表1所示,
表1:C中各物性参数取值
将上述各取值代入公式(2)中计算得到C1=0.648。
计算步骤S3,基于单色激光诱导白炽光纳米粒径计算模型,依据获取的信号特征时间及确定的系数C计算碳烟粒径,纳米粒径计算模型为
其中,C是由聚集体的屏蔽效应参数以及碳烟颗粒与周围气体的热力学参数组成的系数,Tg碳烟颗粒周围气体的温度,Th-max表示碳烟颗粒冷却过程中的最高温度,λ是激光器发射的激光的波长,h表示普朗克常数,c表示光速,k表示玻尔兹曼常数。
通过控制激光加热颗粒,能够在以热传导为主导的冷却过程中从激光诱导白炽光信号确定粒径,使用典型的C2H4/空气扩散火焰来验证该方法。本实施例的激光器发射的激光能量密度为0.09J/cm2,通过功率计监测能量大小和变化情况。由于测量过程中需要保证碳烟的冷却过程以热传导为主导,因此加热碳烟颗粒的最高温度Th-max不得超过碳烟的升华点3700K,本实施例根据激光能量密度的取值,碳烟最高温度Th-max范围为2800-3700K,此状态下颗粒的升华效应可以忽略不计。
利用公式(2)和公式(3)分别计算纳米级碳烟粒径得到Dp-λ和Dp-T,公式(2)和公式(3)中的C根据各物性参数计算均取0.648。将传统的双色激光诱导白炽光***测量的粒径Dp-T和本次单色激光诱导白炽光***测量的粒径Dp-λ进行一致性比较,得到的关系图,其中Dp-405nm和Dp-600nm分别表示滤光片中心波长为405nm和600nm时的Dp-λ值。由图3显示的结果分析发现,Dp-T和Dp-λ有很好的一致性,波长为405nm和600nm时Dp-λ与Dp-T的总体误差分别为3.6%和5.4%。
实施例2
本实施例同样基于实施例1中的单色激光诱导白炽光纳米级碳烟粒径测量装置实现,相同之处不再赘述,本实施例的测量装置还包括对碳烟颗粒进行TEM采样的透镜电镜,通过透射电镜TEM采样观察碳烟粒径Dp,得到的电镜结果如图7所示。本实施例的存储器中计算机程序被处理器执行时能实现图4所示步骤:
获取信号特征时间步骤S1,获取碳烟颗粒的白炽光信号强度降低到其峰值1/e的所需时间,即为信号特征时间τλ,测量火焰中激光诱导白炽光信号的特征时间τλ;
建模步骤S2,建立碳烟粒径Dp与信号特征时间τλ的关系模型;
获取步骤S3,结合透射电镜采样观测得到的碳烟粒径结果,获取碳烟粒径Dp与特征时间τλ之间的关系系数C;
推演步骤S4,基于碳烟粒径Dp与信号特征时间τλ的二元关系系数C,通过测量火焰不同位置处的激光诱导白炽光信号,得出火焰中碳烟粒径的二维分布图,对应得出碳烟粒径。
本实施例的碳烟颗粒的最高温度Th-max设定为2800-3700K。为了验证该范围内Th-max对Dp-λ的测量影响不大,假设在恒定的最高温度Th-max=3600K条件下尝试用τλ测量粒径。图5表示在最低的能量密度0.09J/cm2条件下,粒径Dp随着中心线上火焰高度HAB的变化曲线。结果显示,在405nm和600nm的检测波长中测量的粒径依然有很好的一致性,比较Dp-λ和Dp-T在405nm和600n波长下的相对平均误差分别是5.6%和5.5%。因此可以说明Th-max对Dp-λ的测量影响很小。可以在给定的温度范围2800-3700K假设一个合适的Th-max值,而无需用双色激光诱导白炽光***测量Th-max。
进一步分析选取的Th-max值对粒径误差的影响。图6显示了在Tg=1600K、λ=500nm时由颗粒实际最高温度Th-max的变化导致的粒径误差曲线,不同曲线代表不同假设最高温度,横坐标代表实际最高温度Th-max,纵坐标代表在实际温度和假设温度下求得的Dp之间的相对误差。由图6分析发现,给定温度范围在2800-3700K条件下平均误差在Th-max为3600K时达到最小,因此在本实施例Tg=1600K、λ=500nm的条件下选用Th-max为3600K,但是3600K并不总是对粒径来说是最适合的温度值,例如,当Tg=1800K,λ=700nm时,在Th-max=3600K下误差达到最大13.5%,而在Th-max=3050K时为7.4%。即Th-max的最佳假定值取决于实际的实验条件,通过调整至适当的检测波长和最佳假设最高温度可以降低测量误差。
对C传统的求解过程是根据条件设定C中各物性参数的值,但由于不同的燃烧条件下参数的取值偏差较大,又Dp和C呈线性关系,因此误差对粒径的求解影响很大。本实施例直接得出碳烟粒径Dp与信号特征时间τλ的二元关系,根据该二元关系以及获取待测碳烟颗粒的信号特征时间,对应得出碳烟粒径的二维分布图,直观且方便,大大提高了测量效率。
Claims (10)
1.一种单色激光诱导白炽光纳米级碳烟粒径测量装置,其特征在于:包括火焰发生单元、对火焰的碳烟颗粒进行加热的激光器、用于捕获火焰中白炽光信号的捕获单元、显示并记录所捕获的白炽光信号与时间关系的示波器以及与所述示波器通信连接的处理单元;
所述的捕获单元包括放置在火焰一侧的光电倍增管以及放置在所述光电倍增管前用于阻止火焰中除白炽光信号之外的光信号通过的窄带通滤光片;
所述处理单元包括存储器和处理器,存储器存储计算机程序,该计算机程序被所述处理器执行时能实现以下步骤:
建立火焰中碳烟颗粒的白炽光信号强度与时间的关系,定义碳烟颗粒的白炽光信号强度降低到其峰值的1/e时所需时间为信号特征时间τλ,基于包含碳烟粒径与信号特征时间τλ关系的粒径计算模型得出火焰中碳烟颗粒的粒径。
2.根据权利要求1所述的单色激光诱导白炽光纳米级碳烟粒径测量装置,其特征在于:
所述的纳米级碳烟颗粒粒径计算模型为
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>p</mi>
<mo>-</mo>
<mi>&lambda;</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>C</mi>
<mfrac>
<msub>
<mi>&tau;</mi>
<mi>&lambda;</mi>
</msub>
<msqrt>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
</msqrt>
</mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>-</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>-</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mi>k</mi>
</mrow>
<mrow>
<mi>h</mi>
<mi>c</mi>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
其中,C是由聚集体的屏蔽效应参数以及碳烟颗粒与周围气体的热力学参数组成的系数,Tg碳烟颗粒周围气体的温度,Th-max表示碳烟颗粒冷却过程中的最高温度,λ是光电倍增管的设定探测波长,h表示普朗克常数,c表示光速,k表示玻尔兹曼常数。
3.根据权利要求2所述的单色激光诱导白炽光纳米级碳烟粒径测量装置,其特征在于,
所述的C的表达式为
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msub>
<mi>&alpha;</mi>
<mi>T</mi>
</msub>
<msub>
<mi>p</mi>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>R</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;W</mi>
<mi>g</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>&gamma;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>k</mi>
<mi>h</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>2</mn>
<mo>/</mo>
<msub>
<mi>D</mi>
<mi>h</mi>
</msub>
</mrow>
</msup>
</mrow>
其中,αT表示热调节系数,pg表示碳烟颗粒周围气体的压力,Np表示聚集体中基元体颗粒的数量,ρs是碳烟的密度,cs是碳烟的比热,R是通用气体常数,Wg是碳烟颗粒周围气体的平均摩尔质量,γ是周围气体的比热比,kh和Dh均为碳烟颗粒的缩放系数。
4.根据权利要求1所述的单色激光诱导白炽光纳米级碳烟粒径测量装置,其特征在于,还包括对碳烟颗粒进行TEM采样的透镜电镜,所述的粒径计算模型为碳烟粒径与特征时间的二元关系图,该二元关系图由以下步骤得到:
对于同一工况下的碳烟颗粒,首先在2800-3700K的范围内设定若干个Th-max的值得到不同位置处的若干个信号特征时间τλ,采用透射电镜TEM采样观察并得到这些位置的碳烟粒径Dp,基于实验所设定的探测波长和火焰温度,计算选择最优的Th-max,然后在二维坐标系中标定由τλ和Dp确定的点,并拟合得到碳烟粒径Dp与信号特征时间τλ的二元关系式,最后根据该二元关系式由不同测量位置处的信号特征时间τλ得到火焰中Dp的二维分布结果;
对于不同工况下的碳烟颗粒,按照上述步骤得到不同的二元关系并由各工况下信号特征时间τλ得到Dp的值。
5.根据权利要求1所述的单色激光诱导白炽光纳米级碳烟粒径测量装置,其特征在于:
所述的激光器为脉冲激光器,在高温燃烧测量环境中,所使用的能量密度不超过0.1J/cm2。
6.根据权利要求1所述的单色激光诱导白炽光纳米级碳烟粒径测量装置,其特征在于:
所述的光电倍增管的响应波段和所述的窄带通滤光片的透过波长位于可见光范围;所述光电倍增管的上升时间小于1ns。
7.根据权利要求1所述的单色激光诱导白炽光纳米级碳烟粒径测量装置,其特征在于:
所述的激光器与火焰之间设有用于对正火焰中心的光阑以及用于调节激光器发射的脉冲能量的衰减器;
所述激光器放置在一移动平台上,该移动平台根据光阑的对正信号校正激光器的位置以使其对准火焰中心。
8.一种单色激光诱导白炽光纳米级碳烟粒径测量方法,其特征在于,包括:
获取信号特征时间步骤,采用脉冲激光对火焰中碳烟颗粒进行加热,获取火焰中碳烟颗粒的白炽光信号强度降低到其峰值的1/e时所需时间,即信号特征时间;
确定系数C步骤,C的表达式为
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msub>
<mi>&alpha;</mi>
<mi>T</mi>
</msub>
<msub>
<mi>p</mi>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>R</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;W</mi>
<mi>g</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>&gamma;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>k</mi>
<mi>h</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>2</mn>
<mo>/</mo>
<msub>
<mi>D</mi>
<mi>h</mi>
</msub>
</mrow>
</msup>
</mrow>
其中,αT表示热调节系数,pg表示碳烟颗粒周围气体的压力,Np表示聚集体中基元体颗粒的数量,ρs是碳烟的密度,cs是碳烟的比热,R是通用气体常数,Wg是碳烟颗粒周围气体的平均摩尔质量,γ是周围气体的比热比,kh和Dh均为碳烟颗粒的缩放系数;
计算步骤,基于纳米级碳烟粒径计算模型,依据获取的信号特征时间及确定的系数C计算碳烟粒径,所述粒径计算模型为
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>p</mi>
<mo>-</mo>
<mi>&lambda;</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>C</mi>
<mfrac>
<msub>
<mi>&tau;</mi>
<mi>&lambda;</mi>
</msub>
<msqrt>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
</msqrt>
</mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>-</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>-</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mi>k</mi>
</mrow>
<mrow>
<mi>h</mi>
<mi>c</mi>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
其中,C是由聚集体的屏蔽效应参数以及碳烟颗粒与周围气体的热力学参数组成的系数,Tg碳烟颗粒周围气体的温度,Th-max表示碳烟颗粒冷却过程中的最高温度,λ是光电倍增管的设定探测波长,h表示普朗克常数,c表示光速,k表示玻尔兹曼常数。
9.一种单色激光诱导白炽光纳米级碳烟粒径测量方法,其特征在于,包括:
建模步骤,建立碳烟粒径Dp与信号特征时间τλ的关系模型,定义火焰中碳烟颗粒的白炽光信号强度降低到其原始值的1/e时所需时间为信号特征时间τλ;
获取步骤,测量并计算火焰中激光诱导白炽光信号的信号特征时间τλ,并结合透射电镜采样观测得到的碳烟粒径结果,获取碳烟粒径Dp与信号特征时间τλ之间的关系系数C;
推演步骤,基于所述碳烟粒径Dp与信号特征时间τλ的二元关系系数C,根据火焰不同位置处的激光诱导白炽光信号的信号特征时间τλ,得出火焰中碳烟粒径的二维分布图,对应得出碳烟粒径。
10.根据权利要求9所述的单色激光诱导白炽光纳米级碳烟粒径测量方法,其特征在于:
所述的碳烟颗粒的最高温度设定为2800-3700K。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710984620.7A CN107976384B (zh) | 2017-10-20 | 2017-10-20 | 单色激光诱导白炽光纳米级碳烟粒径测量装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710984620.7A CN107976384B (zh) | 2017-10-20 | 2017-10-20 | 单色激光诱导白炽光纳米级碳烟粒径测量装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107976384A true CN107976384A (zh) | 2018-05-01 |
CN107976384B CN107976384B (zh) | 2019-07-09 |
Family
ID=62012515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710984620.7A Active CN107976384B (zh) | 2017-10-20 | 2017-10-20 | 单色激光诱导白炽光纳米级碳烟粒径测量装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107976384B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108593530A (zh) * | 2018-06-01 | 2018-09-28 | 大连海事大学 | 微颗粒前向散射光检测装置及其检测方法 |
CN110702568A (zh) * | 2019-10-09 | 2020-01-17 | 浙江大学 | 一种激光诱导炽光微米级含碳微粒质量浓度测量装置及方法 |
CN111896436A (zh) * | 2020-08-11 | 2020-11-06 | 哈尔滨工业大学 | 碳黑团聚体一次粒径分布和热适应系数同时测量方法及装置 |
CN113310857A (zh) * | 2021-05-27 | 2021-08-27 | 上海交通大学 | 湍流燃烧场四维碳烟颗粒初级粒径分布测量***及其方法 |
CN113348359A (zh) * | 2019-01-23 | 2021-09-03 | 罗伯特·博世有限公司 | 用于在使用激光诱导白炽光原理的情况下探测流动流体中的颗粒或气溶胶的传感器装置 |
CN113758840A (zh) * | 2021-09-07 | 2021-12-07 | 安徽工业大学 | 一种基于Image J的碳烟颗粒粒径及面积的测量方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101200005A (zh) * | 2007-09-04 | 2008-06-18 | 广西师范大学 | 简便快速制备Au-Ag复合纳米微粒的方法 |
CN102003936A (zh) * | 2010-09-14 | 2011-04-06 | 浙江大学 | 同时测量液滴位置、粒径和复折射率的方法和装置 |
CN201974157U (zh) * | 2010-12-22 | 2011-09-14 | 浙江大学 | 颗粒场紧凑式数字全息装置 |
CN102248177A (zh) * | 2011-07-29 | 2011-11-23 | 上海龙翔新材料科技有限公司 | 激光诱导球型银粉制备方法 |
CN103698256A (zh) * | 2013-12-25 | 2014-04-02 | 浙江大学 | 一种全场彩虹在线测量液体喷雾的方法和装置 |
CN103760142A (zh) * | 2014-01-14 | 2014-04-30 | 中国科学院工程热物理研究所 | 一种燃油喷嘴液滴空间分布光学测量方法及装置 |
CN103940709A (zh) * | 2014-05-06 | 2014-07-23 | 南京中科神光科技有限公司 | 一种实时微生物粒子计数器 |
-
2017
- 2017-10-20 CN CN201710984620.7A patent/CN107976384B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101200005A (zh) * | 2007-09-04 | 2008-06-18 | 广西师范大学 | 简便快速制备Au-Ag复合纳米微粒的方法 |
CN102003936A (zh) * | 2010-09-14 | 2011-04-06 | 浙江大学 | 同时测量液滴位置、粒径和复折射率的方法和装置 |
CN201974157U (zh) * | 2010-12-22 | 2011-09-14 | 浙江大学 | 颗粒场紧凑式数字全息装置 |
CN102248177A (zh) * | 2011-07-29 | 2011-11-23 | 上海龙翔新材料科技有限公司 | 激光诱导球型银粉制备方法 |
CN103698256A (zh) * | 2013-12-25 | 2014-04-02 | 浙江大学 | 一种全场彩虹在线测量液体喷雾的方法和装置 |
CN103760142A (zh) * | 2014-01-14 | 2014-04-30 | 中国科学院工程热物理研究所 | 一种燃油喷嘴液滴空间分布光学测量方法及装置 |
CN103940709A (zh) * | 2014-05-06 | 2014-07-23 | 南京中科神光科技有限公司 | 一种实时微生物粒子计数器 |
Non-Patent Citations (3)
Title |
---|
吴学成: "基于轨迹图像的煤粉颗粒速度和粒径测量", 《浙江大学学报》 * |
岳宗宇: "激光诱导炽光法定量测量碳烟", 《燃烧科学与技术》 * |
潘振艳: "激光诱导辐射法测定燃煤氢气火焰中碳粒子的研究", 《中国激光》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108593530A (zh) * | 2018-06-01 | 2018-09-28 | 大连海事大学 | 微颗粒前向散射光检测装置及其检测方法 |
CN113348359A (zh) * | 2019-01-23 | 2021-09-03 | 罗伯特·博世有限公司 | 用于在使用激光诱导白炽光原理的情况下探测流动流体中的颗粒或气溶胶的传感器装置 |
CN110702568A (zh) * | 2019-10-09 | 2020-01-17 | 浙江大学 | 一种激光诱导炽光微米级含碳微粒质量浓度测量装置及方法 |
CN110702568B (zh) * | 2019-10-09 | 2020-08-25 | 浙江大学 | 一种激光诱导炽光微米级含碳微粒质量浓度测量装置及方法 |
CN111896436A (zh) * | 2020-08-11 | 2020-11-06 | 哈尔滨工业大学 | 碳黑团聚体一次粒径分布和热适应系数同时测量方法及装置 |
CN113310857A (zh) * | 2021-05-27 | 2021-08-27 | 上海交通大学 | 湍流燃烧场四维碳烟颗粒初级粒径分布测量***及其方法 |
CN113758840A (zh) * | 2021-09-07 | 2021-12-07 | 安徽工业大学 | 一种基于Image J的碳烟颗粒粒径及面积的测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107976384B (zh) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107976384B (zh) | 单色激光诱导白炽光纳米级碳烟粒径测量装置及方法 | |
Sun et al. | Single-shot, Time-Resolved planar Laser-Induced Incandescence (TiRe-LII) for soot primary particle sizing in flames | |
CN105606222B (zh) | 一种火焰三维温度场的测量装置及测量方法 | |
JP4563525B2 (ja) | レーザー誘起白熱光における絶対光度測定 | |
WO2021057860A1 (zh) | 基于置障条件下天然气长输管道喷射火研究的实验装置 | |
CN105372263B (zh) | 一种电涡流加热式管道缺陷红外成像检测装置及检测方法 | |
WO2018120444A1 (zh) | 一种动目标的红外辐射光谱特性仿真分析方法 | |
CN109655161B (zh) | 基于红外热像仪的目标红外积分辐射强度测试方法及设备 | |
Ma et al. | Combined soot optical characterization using 2-D multi-angle light scattering and spectrally resolved line-of-sight attenuation and its implication on soot color-ratio pyrometry | |
JP2009518584A (ja) | リアルタイム粒子状物質計測システム | |
Cenker et al. | Assessment of soot particle-size imaging with LII at Diesel engine conditions | |
US20040141176A1 (en) | Small particle analysis by laser induced incandescence | |
CN109655017A (zh) | 一种管道同轴度测量*** | |
US9196032B1 (en) | Equipment and method for three-dimensional radiance and gas species field estimation | |
Meléndez et al. | Infrared thermography of solid surfaces in a fire | |
Kheirkhah et al. | Development and validation of a multi-angle light scattering method for fast engine soot mass and size measurements | |
CN104181172B (zh) | 一种固体推进剂燃烧喷焰自由电子浓度测试方法 | |
US10012628B2 (en) | Multifunctional particle analysis device and method of calibrating the same | |
CN115639124B (zh) | 用单波长激光实现碳烟浓度测量和原位标定的装置及方法 | |
CN114034405B (zh) | 一种非接触式测温方法及*** | |
CN113310857B (zh) | 湍流燃烧场四维碳烟颗粒初级粒径分布测量***及其方法 | |
CN111366254A (zh) | 一种大气透过率检测方法及装置 | |
Kelly et al. | A novel multi-band plenoptic pyrometer used for temperature measurements of strand burner plumes | |
CN209372158U (zh) | 干熄炉空气导入流量的测量*** | |
CN102253234B (zh) | 本生灯下气体燃料层流火焰传播速度的精确测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |