刘训臣燃烧诊断实验室主页
  • Alumni
  • People
  • Post
  • Research

Post

科研笔记和一些想法

三月 17, 2026

光谱模拟

基于本实验室之前发布的对于测量燃烧气体的吸收率以及光谱透射率的界面,我在其基础上再次做出了相应的更新: 如有需要请下载下方新的压缩包 flame_spectrum.zip

  1. 气体测量时浓度的设定之前不是很精确,因此更新出用浓度定量来设定气体含量,采用ppm,更加精确。
  2. 对于透射率的设定,有些时候我们仅需要某种气体透射率或者不需要透射率,只看透射系数,我们也可设定某气体的或者不看透射率,只看吸收系数。
  3. 可以手动调节生成图网格划分的精细程度或者图形尺寸的大小。

接下来我详细说明一下如何操作新更新的功能

可以看到手动调节参数组后面跟着气体的比例以及其浓度,因此我们可以调节其浓度来设置其比例。

光谱模拟

看到如图所示除了之前可以设置波数范围,我们也可勾选,可以看到波长的范围,也可以看到分辨率的设置,其次透射率以及吸收系数我们也可以进行选择,对于透射率也可选择某种气体或者全部的透射率。

透射率显示

对于某些我们只想看大体走势或者趋势,不在意某条谱线的情况下,我们可以选择低的网格密度或者图像尺寸,以此来节约时间。或者反之,我们也可以看到的更加精确。

图形显示

最后的模拟结果如下图

加入光谱模拟结果

以上就是更新的内容,所有需求都是自己在实验过程中所设计需要改进的地方,希望能够为同伴提供帮助。

继续阅读
三月 8, 2026

声激励喷嘴的结构与参数设计

今天和大家聊聊声激励脉动系统的结构与参数设计,这也是我的相关实验的核心内容,声激励脉动系统的核心是通过特定结构设计,选择合适的将声能与流体动能结合,观察声流扰动对火焰熄火边界的影响,其结构的设计和参数阈值的选择,直接决定最终激励效果,所以设计环节需要重点把控。

声激励喷嘴

通过阅读相关文献,可以看到前人研究声激励对火焰形态影响的实验装置,因此,本次声激励喷嘴的结构设计,将参考前人的研究,进行声激励喷嘴的设计,如图1为前人所采用的实验装置示意图。如图我们可以看到,前人是在腔体底部安装一个扬声器,在腔体部分通入甲烷和空气通过预混后在出口处直接形成火焰,通过高速相机拍摄火焰形态,对火焰形态进行研究。本研究的声激励实验装置是用于横向射流火焰的研究,因此需要设计适用于本实验的实验装置。

前人实验装置图

需要优先对扬声器进行选择,因为腔体尺寸需要匹配扬声器的参数进行选择,扬声器参数如图2所示。

扬声器参数

由于扬声器需要固定在腔体当中,因此,腔体的直径选择需要大一些,并且在腔体内部安装一个法兰,用于固定扬声器,因此选择直径为325mm的不锈钢腔体,并且在内部高度合适的部位焊接法兰。如图3所示为需要焊接在腔体内部的法兰,其上面分别有八个孔用来固定扬声器,焊接高度距离底部为140mm,适应扬声器的高度。

腔体内部法兰

另外需要在腔体的地步进行密封处理,所以需要在尾部焊接一个法兰并且用一块大小相同的不锈钢板进行尾部的封闭,并且在法兰除使用一个规格相同的密封垫片,用于保证装置的密封性。尾部法兰与密封不锈钢片的尺寸如图4、图5所示,密封垫片与法兰的规格相同,待扬声器连接好后进行密封。

尾部法兰

尾部密封不锈钢片

最终需要将腔体转到10mm的细管上接入横向射流火焰实验装置,但市面上没有合适的漏斗转接,因此通过两次转接最终将235mm的腔体转到10mm的细管上,用于接到横向射流火焰的实验装置上。

另外需要在腔体部位开孔通入空气,并且在细管出口处连接一个三通,一端连接细管,一段连接实验装置,剩余的一端用于连接甲烷。并且需要为扬声器供电,在腔体部分开孔安装航插并且进行密封,即可为扬声器供电使其正常工作。

声激励喷嘴的装配图如图6所示。

声激励喷嘴实物图

功率放大器

根据扬声器的参数,对功率放大器的规格进行选择,选择功率适用于扬声器的功率放大器,确保扬声器能正常工作。

信号发生器

关于信号发生器的选择,需要适配于扬声器的频率范围,该设备中选用10MHz的信号发生器作为输入源。

经过该三个部分,构成了完整的声激励脉动系统,用于实验当中,后续可以采用该系统进行实验,并且可以采用不同相位的激励来观察火焰的燃烧情况,待后续有实验结果还会继续与大家分享。

继续阅读
三月 6, 2026

2026年3月课题组文献讨论


by 张洪民

Wang, Shaojie, Shengming Yin, Weifan Hu, 等. 《MHz-rate laser absorption spectroscopy for temperature and N2O concentration measurements in laser-induced ADN propellant decomposition》. Applications in Energy and Combustion Science 26 (2026年6月): 100469.


创新点

MHz速率激光吸收光谱(LAS)与高速摄影结合

文章主要是将MHz速率激光吸收光谱(LAS)与高速摄影结合,借助 10.6 μm CO₂激光引发 ADN 推进剂快速热分解,通过带偏置三通电路的带间级联激光器(ICL) 提升扫描深度。

主要装置如下图

装置图

加热腔采用不锈钢材质,两侧安装硒化锌(ZnSe)窗口,供红外激光穿过;加热腔前端设置蓝宝石窗口,便于高速摄影记录热分解过程;后端设有进、出气管路,用于通入环境气体和排出废气;加热腔正上方还安装有一个硒化锌窗口,供 10.6 μm 二氧化碳激光器对 ADN 样品进行加热。

问题点

但是我们需要注意到,在其对于N2O测量的数据图,如下面两个图中,其纵坐标温度为470k多一点,而对于一个作为高能量密度的新型含能材料这显然是不符合实际情况的。

测量光谱(1)

测量光谱(2)

究其原因,是因为测量N2O这一个中间产物的环境是密封的采用了吸收池,因为本实验的研究对象需要一个加热的温度较高的环境,因此采用的密封吸收池导致实验测量产物N2O大量聚集,而我们采用激光测量谱线的时候,通常是固定的波数测量范围,因此在激光测量时测量的N2O遍布整个吸收池,这就导致测量数据图并不是真实的分解温度而是激光光程内的温度的平均拟合值,因此出现这种情况。

具体的解决方法也很简单,只需要保证环境不密封或者其中N2O可以及时快速的流出即可,快速的吹扫积累的N2O,保证测量环境的准确性减少对于真实实验环境温度的影响。

点评:本质上,这还是 TDLAS的空间分辨率问题。 在腔体里面用TDLAS测燃烧产生的N2O,很容易就测到的是扩散到环境中的,这是常见的误区。


Dharmaputra B , Nagpure P , Impagnatiello M ,et al.Flame transfer function measurement of a sequential combustor fuelled with natural gas and hydrogen[J].Combustion and Flame, 2025, 274(000).

继续阅读
二月 28, 2026

PID实现流量控制

在实验使用空压机时,总是会出现上游压力不稳导致流量波动幅度大,压缩空气流量大,大量程的质量流量控制器价格过于昂贵,ALICAT-2000slm量程的需要将近五万元,于是考虑用流量计,PID控制器和电控开度阀自制一个控流装置,以实现稳定流量。

流量计可以输出4-20ma电流信号,将输出的电流信号接到PID控制器,pid设置上下限量程,由pid控制器输出控制信号,控制电控阀的开度。 管路连接顺序为,空压机气罐→流量计→电控球阀,流量计的电流信号作为pid的信号输入端,pid接收电流信号后表显流量数值,将控制指令以电流信号输出给电控阀,电控球阀根据指令实现调流稳压。

设备如图设备图

其中,流量计及电控阀需要接DC24V电源,pid控制器需要接85-265VAC,pid接线如图, 接线图

将pid控制器的上下限设置成与流量计量程一致,pid控制器的pv为实时信号换算值,sv为预设值,可根据需要设置成对应数值。 刚开始先使用25s电控阀进行实验。结果流量计数显正常,pid接收电流信号后表显数值与流量计相符,开度阀能正常工作,但开度阀会因流量超过预设流量而全部关闭然后打开,以此往复,未能达到预期结果。 猜测导致实验未达目的的原因是:电控开度阀的动作时间有滞后特性,pid的调节周期(采样频率)短,导致: 1.控制器频繁检测到流量波动,反复输出开度指令,阀门还没完成上一次动作,新的指令又到来,形成“指令打架”; 2.将控制器上下限设为0-1000(与流量计量程一致),阀门的开度指令范围应该是0-100%(对应4-20mA),控制器输出的“0”被阀门识别为“全关”,“1000”被识别为“全开”;当偏差出现时,PID输出容易直接触达0或1000的边界,导致阀门极限动作(全关/全开),而非平滑调节,加剧震荡。

在更换1s快调阀之后,再次进行第二次实验。与之前慢阀门的对比,阀门开关明显变快,开始工作后,上游流量充足,阀门部分开启,流量没有准确达到预设流量,会全部关闭后再重新增大或缩小阀门开度,以此往复,直至准确调整至预设流量,但这种全关的调节,需要消耗大量时间才能使流量达到预设流量,而且每次更改预设流量,都要重新调节,这显然是存在问题的。

次日,进行了第三次实验。通过修改pid指令限幅,阀门全关问题解决,试了几个流量,流量基本能稳稳的控制在设定值,如图,

实验结果图 达到预期成果。有时流量会在区间内波动,可能与上游充气引起压力变化导致的。

数据结果图

此后,进行多次实验总结了使用这套稳流装置的顺序:先开阀给充足的上游流量,然后调限幅,进行PID自整定(之前自整定之所以时间太长是因自整定时流量要超过几次预设值,流量未达到导致自整定失败,只有在实际使用环境才需要自整定。)可知道流量大概上下限,再微调限幅(OPI=FREE,OPL最小电流,OPH最大电流),然后就可以低幅度调流,实现稳流。在使用这个装置以后,上游的流量明显比之前稳定。

继续阅读
一月 12, 2026

Hencken炉的最新进展

自2025年6月12日在连云港实验过第一次设计的Hencken炉后,出现了一些问题,后面又进行了改进,现在再更新一次最新的进展。

在连云港实验过程中发现的主要问题:由于氢氧燃料口排布间距过远,导致氢氧气的混合不均匀,没有达到理想的温度;进气腔没有做一个类似于静压室的结构导致气体流速在不同燃料口不均匀,进而导致燃烧不均匀不稳定。 当时主要的实验条件:氢气流量:50slm-500slm、氧气流量:25slm-150slm;前面几组是在当量比为1的条件下进行试验,之后进行大流量时,将氢气调大。燃烧情况如图: 氢氧流量分别为200-100时的燃烧情况 氢氧流量分别为300-150时的燃烧情况 后面针对遇到的问题进行改进,首先尽量减少燃料口的间距,之前的气管尺寸为1.3mm、间距为3.5mm,现在气管尺寸为0.7mm、间距减小为1.5mm,这样带来的问题就是水冷只能在外围覆盖,不能打水冷孔到内部;对于燃烧不均匀要加入静压室结构的改进,主要是在进气管进去炉体内部后在前方增加一个隔板,让气体先在隔板隔离出来的空间内稳定,然后再流入燃料管。 现在的结构如图: 更改后的燃烧图 同时连接更改前与更改后的炉子查看燃烧情况,可以看出,修改后的炉子燃烧更加均匀且在相同压力下流速更快 对比图 对比动图

继续阅读
十二月 25, 2025

I2-cell

按照如下的图片制作了碘吸收池 流程图

碘的升华温度为30度,蒸汽压力与温度的关系为 蒸汽压公式

其中P为碘的饱和蒸汽压,原来的压力单位是有问题的,不是pa,而为mmHg。 t为冷端温度,单位(℃) 以上部公式可得出,按照常规的制作方法,把碘密封在玻璃瓶中抽真空,因碘在室温下也会挥发,因此抽真空高度不高,同时会造成系统的污染。

制作流程:

  1. 按照图示接好排气系统
  2. 装固态碘:先在液氮筒②内灌入液氮,然后从 d 处把金属膜装入管内,这种方法能使碘蒸汽冷却在管内,不至于污染系统。碘的多少按需要放,装完后用火焰把 a 处封死。
  3. 排气。排气将液氮灌满液氮筒①,然后将液氮筒②往下移,只冷却住底部的碘,然后开启机械泵,扩散泵进行排气。液氮筒①的作用是进一步防止碘蒸汽进入排气系统。在抽气过程中不断地灌满液筒①、②,还可以用煤气火烤碘管的上面部分以保证吸收池的真空度。
  4. 封碘吸收池。封碘吸收池排气完成后,在收口 b 处把吸收池与碘管一起封下,然后用火焰蒸馏碘,把碘蒸汽赶入吸收池内,再在 c 处封下吸收池,这样一个吸收池就制成了。 对于充入碘蒸汽压有严格要求的,将 b 处封离后的吸收池及碘管将放在可调温度的水浴加热容器内,并按碘蒸汽压与温度关系式调节相应的所需温度,使吸收池内的蒸汽压达到饱和值;然后用细小的火焰快速在 c 处封下吸收池。通过以上工艺,我们制成的碘吸收池已用于激光研究工作中,解决了真空系统被碘蒸汽污染的技术难题。

需要注意的问题:

  1. 如想要碘蒸馏瓶,在与吸收池连接时不行使用橡胶或硅胶管连接碘蒸汽会与其发生反应
  2. 在进行吸收池密封后退火时,要先将表面水分或其他杂志清除干净,避免退火后在外壁产生污点

制作的碘吸收池如图: 碘吸收池

继续阅读
十二月 15, 2025

干涉器制作

使用传统的法布里-珀罗标准具如图所示 激光器 在使用过程中需要将激光垂直入射标准具,那使用前就需要进行方向的调节,这个过程就会显得繁琐。因此我们制作了一个更为简单方便的能够让光进行干涉的器件: 将光分为两路,这两路光的光程不相同在稳定的光程差的条件下,光即可产生稳定的干涉。 此外单单只进行分光,会存在一个问题,即分光器不会对激光的传输方向产生显示,在我们将光一分为二再合二为一后,光在光纤内部会正向反向混合传输,导致信号被淹没,因此需要在每一个支路上添加光纤环形器,让光只沿一个方向传播,到最后产生稳定的干涉光。 原理图如图所示: 原理图 实物图如图所示: 实物图 最终所测得的光的干涉如图所示: 测量结果 测量结果 没加环形器的测试结果: 测量结果 在这其中还有几个问题需要注意,首先是要保证光程差的稳定,否则干涉光会产生一定的波动,因此我们在测试没问题后,用胶水将所有光路进行固定;其次是由于加入了环形器因此所制作的干涉器也就有了明确的输入端和输出端;最后是光程差的选择我们首先测试了20m,测量结果并不理想,但初步猜测是由于数据采集卡采集频率不够高的问题,然后我们又是试验了在没有明显的光程差的情况下(不单独加光纤,将原来的一条支路进行微小的拉长)、在2m光程下的结果,最终2m的产生会更稳定一些。20m的测试结果如下: 20m测量结果 20m测量结果

继续阅读
十二月 11, 2025

PIV 流场分析软件介绍

在流体力学、空气动力学、燃烧诊断等科研领域,流场数据的精准计算与深度分析是核心环节。传统数据处理流程繁琐,算法选择受限,往往耗费大量时间精力。今天给大家介绍的软件,集成四大核心功能模块,无需复杂编程,就能实现从原始图像到流场动态特性分析的一站式解决方案,让科研效率直接翻倍!

该软件基于 PyQt6 打造直观友好的图形界面,整合了先进的计算机视觉与数据模态分解算法,专为科研场景设计。无论是单组图像的速度场计算,还是多序列图像的批量处理,亦或是流场数据的深度模态分析,都能轻松应对。无需编写一行代码,只需通过简单的参数设置和鼠标操作,就能完成复杂的流场计算任务。

该软件有四大核心功能模块,覆盖流场分析的全流程。

软件主界面

双图片速度场计算

作为软件的基础核心模块,该功能支持从两幅连续时刻的流场图像中,计算高精度二维速度场(U、V 分量),提供两种主流算法供选择:

  • 光流法(Farneback算法):基于 OpenCV 实现,计算速度快、对硬件资源要求低,适合粒子密度大、有大片连续灰度值的流场场景,比如常规流体流动观测。
  • 神经网络法(PIV-LiteFlowNet-en 模型):依托预训练深度学习模型,精度更高,对复杂流动模式和低信噪比图像的鲁棒性更强,适合复杂流场研究,但需 GPU 支持以保证运算速度。

计算完成后,会实时生成 U 分量云图、V 分量云图、速度矢量场图和速度分布直方图,还能通过滑块调节矢量箭头长度和速度阈值,让结果可视化更清晰。

双图片速度场计算模块界面

多图像序列批处理

针对按时间序列排列的多幅流场图像,该模块可实现自动化批量处理,无需手动逐组计算。它会依次对(img1,img2)、(img2,img3)…进行运算,最终生成 n-1 个速度场数据。

同样支持光流法和神经网络法,内置进度监控和内存管理功能,即使是长时间批量处理或使用耗资源的神经网络法,也能保证稳定性。生成的速度场文件会自动按序号命名保存,大幅提升数据处理效率,特别适合长时间序列流场观测实验。

多图像序列批处理模块界

在选用以上两种方法时,可以根据以下场景选择:

  • 优先选光流法:实时性优先、场景简单稳定、资源有限、无标注数据。
  • 优先选神经网络法:复杂场景、高精度需求、可接受高计算成本、有充足训练数据。

DMD 分析模块

动态模态分解(DMD)是流场分析的关键技术,该模块采用先进的 BOP-DMD 算法,能从流场时间序列数据中,提取主导的时空演化模态、频率及增长率 / 衰减率,助力研究人员挖掘流场潜在规律。

支持加载多种格式的流场文件(.npy、.csv、.dat、.txt),分析结果通过三大可视化图表呈现:

  1. 复平面上的特征值图(直观展示模态稳定性)
  2. 前 9 阶时间模态演化曲线
  3. 前 9 阶空间模态分布云图

还能对比原始数据与 DMD 重构数据,甚至进行短期流场预测。

DMD 分析模块界面

DMD 分析模块特征值

DMD 分析模块时间模态

DMD 分析模块空间模态

文件格式转换

流场数据处理中常面临格式不兼容问题,该模块提供 .npy 与 .csv 格式的批量转换功能:

  • .npy 格式:二进制文件,体积小、读写速度快,适合软件内部计算和数据存储。
  • .csv 格式:文本文件,可读性强,可用 Excel、记事本打开,方便数据交换、人工检查或导入其他不支持 .npy 格式的软件。

操作简单,只需选择目标文件夹和转换方向,一键即可完成批量转换,且完全保留原始数据的 U、V 分量和位置信息。

继续阅读
十二月 11, 2025

光谱模拟

基于本实验室之前发布的对于测量燃烧气体的吸收率以及光谱透射率的界面,我在其基础上做出了相应的更新:

  1. 在原本只有CO2,H2O分子的基础上,添加了CO,NO,N2O以及NO2气体的测量。
  2. 更新了压力调节范围现在其可调压力范围处于10pa-10Mpa
  3. 现在界面可以自主选择使用哪个气体,如果不需要某种气体,可任意删除,如果需要哪种气体,也可自主添加。

压缩包如下,请自行打包下载使用。 flame_spectrum.zip

接下来我会详细的说明一下该界面如何使用: 首先,需要在终端中运行py文件,建议是在anaconda的运行终端中运行运行出来的界面如下图所示。

平衡界面

可以看到左侧栏目中有反应机理,初始条件以及组分设置的选择。我们首先点击初始化气体对象的选项,加载我们的反应机理,确定后续反应的原理。然后我们选择组分设置中的输入方式以及预设,输入方式一般采用当量比,预设中有多种选项,后续我们以甲烷-空气反应为例。选择完毕后,就可以点击计算平衡,我们就可以看到右侧热力学状态对比以及主要组分中分内容,如下图所示。

热力学状态及组分

接下来我们就可以点击传输到光谱模拟,我们的界面就会跳转到光谱模拟,其界面如下图所示:

光谱模拟

可以看到现在只有H2O跟CO2的文件,我们可以点击添加分子来选择加入哪个分子,点击完毕以后就会弹出选择界面,我们下面以CO为例,当我们加入CO以后,我们需要在手动设置参数组设置我们需要的温度,压强,气体浓度等等,并且在下面波数范围以及光程长度,我们同样手动设置一下。在设置完成之后,点击加载分子数据以及开始计算按钮,我这里自己随意设置一下参数,可以看到右图数据如下图所示。可以看到在下图中不同气体利用不同颜色已经标注,非常明显看出气体在不同波数下吸收程度以及透射率等参数变化情况,我们可以尝试多加入几种气体观察一下,如下图所示。

加入CO光谱模拟

混合气体光谱模拟

最后我们也可以删除不需要的气体浓度,总之该方法能为从事气体光谱模拟检测研究工作的研究者提供一定便利,我们会继续改进,后续如有问题欢迎留言。

继续阅读
十二月 9, 2025

反应流体光谱仿真程序 flame spectrum

我们的反应流体光谱仿真程序分两步。

第一步,气体的化学平衡计算,用Cantera 计算。输入温度压力组分,选择化学平衡的约束条件。得到平衡后的热力学状态,比如火焰温度。

化学平衡计算截图

第二步,光谱仿真,用HITRAN计算。化学平衡后得到的温度压力主要组分浓度参数,用来做光谱仿真。主要组分是H2O和CO2 。

当然也可以直接在这个光谱仿真页面输入组分浓度和温度压力参数,直接仿真光谱。

下面图是2um附近仿真,可以看到水的吸收比较强,大概比CO2高了一个数量级。

2um附近水、二氧化碳混合气体光谱仿真截图

下面图是我们之前4um激光测同轴扩散火焰的仿真。可见此处二氧化碳吸收强,而且没有水的干扰。

4.2um附近水、二氧化碳混合气体光谱仿真截图

这个软件框架很灵活,只需要选择不同的数据库就可以模拟不同波段的光谱。

请打包下载几个 py 文件和示例中的par 文件。 flame_spectrum.py 调用了 voigt_simulation 文件夹中 hitran_spectrum_dual.py (可以单独运行对H2O和CO2混合气体仿真)。

flame_spectrum.zip

有问题和建议,请随时留言或者发邮件咨询。

未来新版本软件将包含更多的分子,特别是HITRAN数据库中没有的分子或者波段。比如紫外波段很多自由基还没有高分辨率光谱工作以支持HITRAN数据库条目。完全可以用 NASA JPL数据库 产生HITRAN数据来做。

第一部分反应流体的仿真用Cantera 就可以了,准备主要把第二部分做丰富。 燃烧仿真,湍流燃烧还是很难算的,Fluent 可以大概算算。用 Cantera 和 Chemkin 算一般的反应动力学计算,比如火焰温度。 分子光谱仿真就很小众了,图形化的有 pgophor,最全的应该是 SPFIT/SPCAT 。LIFbase是专门的几个双原子分子的LIF光谱数据库仿真。 几何光学仿真用 zemax,但感觉一般也不需要。

通过一段时间的积累,我们的程序就可以自动化仿真出不同火焰的不同波段光谱了。 Gaydon 在 《The Spectroscopy of Flames》中就给出了很多不同火焰的光谱plate。 希望我们能有一天可以给出这些光谱的仿真,达到对火焰比较好的理解。

继续阅读
  • ««
  • «
  • 1
  • 2
  • 3
  • 4
  • »
  • »»
© 刘训臣燃烧诊断实验室主页 2026