一种基于磁粒子跟踪的流化床内流场测速方法技术

本发明专利技术公开了一种基于磁粒子跟踪技术的流化床内流场测速方法，在磁粒子跟踪测速时，将具有磁性的粒子在理论上等效为磁偶极子，通过实时测量磁粒子运动时其周围产生的磁矢量场，估算得到磁偶极子的位置矢量和磁矩矢量，单个采样间隔内磁粒子的位置矢量差，即可计算出该磁粒子的瞬时平移速度，获取单个采样间隔内磁粒子的磁矩矢量旋转角度，即可计算出该磁粒子的瞬时旋转速度,所述瞬时平移速度、瞬时旋转速度即代表流化床内流场当前的局部流速；本发明专利技术的方法可以不依赖于光学，通过磁性信息跟踪粒子，可以在完全不透明的环境中使用，仅需要三轴磁传感器阵列和高精度数据同步采集卡，试验成本较低。试验成本较低。试验成本较低。

【技术实现步骤摘要】
一种基于磁粒子跟踪的流化床内流场测速方法


[0001]本专利技术属于流场测量领域，具体涉及到一种基于磁粒子跟踪的流化床内流场测速方法。

技术介绍

[0002]流动显示是认知流场内流动介质间相互作用机理、分析流动现象、改进流动设计的重要试验手段，通过各种流动显示实验，可以了解复杂的流动现象，探索其物理机制，为人们发现新的流动现象、建立新的概念和物理模型提供依据。
[0003]随着光学、激光技术、电子技术以及计算机技术的快速发展，复杂流动内部结构的显示能力及流动信息定量提取和分析处理能力得到了长足的进步，多种流动显示方法已用于科研实验和工程实际，如水洞中常用的有染色线、空气泡、氢气泡、激光-荧光灯方法；低速风洞中常用的有烟线、氦气泡、荧光微丝、油流、激光片光等方法；高速风洞中常用的荧光油流、彩色纹影、全息干涉等方法；流化床中常用的数字图像分析、激光多普勒测速、放射性粒子自动跟踪以及磁共振测量等方法。近年来，由于新一代飞行器研制、新型化学工艺、医学及石油等工业对复杂流动认知的迫切需要，以散班技术、粒子图像测速(Particle Image Velocimetry,PIV)、激光诱发荧光和层析技术为代表的新一代非接触式流动测量与显示技术脱颖而出。该类技术因其非侵入式、高精度、可视化、测量信息丰富等优点得到快速发展，兼有定性显示和定量测量的能力，在显示流动结构的同时，可测量密度、温度、速度、压力和浓度等参数，以PIV为例的部分测量方法在功能和测量适应性上优势显著，逐步脱颖而出，技术不断发展完善，发展了二维/二分量、二维/三分量、全息/三分量、解析时间、微型等多种PIV测量形式，形成了多种品牌的PIV商业化产品。
[0004]随着流动机理研究的深入以及流动细节试验测量的需要，人们开始不仅仅关注流场内粒子的平移速度，也同时关注粒子的旋转速度，以期望更精细地研究流动细节。而以PIV为主的图像处理方式在测量粒子的旋转速度时，无法对颗粒球形度较高或者颗粒表面特性比较均匀的颗粒进行旋转运动观察，只能对具有表面局部特征的颗粒或者形状不规则的颗粒进行旋转速度的测量，且测得的旋转角度也仅仅为某一个平面映射下的角度。

技术实现思路

[0005]本专利技术的目的是提出一种能够同时测量粒子平移速度和旋转速度的方法，为流化床内复杂的流动现象提供新的流动显示测量方法。
[0006]为了实现上述目的，本专利技术采用如下技术方案：
[0007]一种基于磁粒子跟踪的流化床内流场测速方法，包括以下步骤：
[0008]S1:在流化床外壁上设置三轴磁传感器阵列，该阵列中三轴磁传感器的数量不少于两个；
[0009]S2:制备含一个含有磁性的颗粒，使其与其他无磁性颗粒的外部形状保持一致，并将磁性颗粒与若干个无磁性颗粒混合后设置在流化床内部；
[0010]S3：同步采集三轴磁传感器阵列的信号，测量得到磁粒子在三轴磁传感器位置处产生的磁矢量场；
[0011]S4：利用磁粒子在空间任一点产生的磁矢量场的理论计算公式和实际测得的磁传感器阵列的磁矢量场，建立目标函数，所述目标函数为理论值与实测值差量的二范数最小构建目标函数；
[0012]S5：利用Levenberg-marquart算法，使得目标函数值最小，求解得到磁粒子的位置矢量和磁矩矢量；
[0013]S6:获取在设定时间范围内磁性粒子的最终位置变化的位置矢量差，即可计算出该磁粒子的平移速度，获取在设定时间范围内磁性粒子的最终位置变化的磁矩矢量的旋转角度，即可计算出该磁性粒子的旋转速度。
[0014]在上述技术中，包括以下步骤：
[0015]A1:在流化床的外部设置有三轴磁传感器，磁传感器的位置坐标为(x
k
,y
k
,z
k
)
T
，k＝1
…
N,其中：N为三轴磁传感器的数量；
[0016]A2：采集无磁粒子情况下的流化床所处位置的背景磁场
[0017]A3：在流化床内设置一个磁粒子；
[0018]A4：开启流化床的流场控制设备使得磁粒子随流场运动，同时启动磁场采集设备，实时采集磁粒子运动时在磁传感器位置处产生的磁矢量场，每个磁传感器测得的实际磁矢量值分别记为
[0019]A5：磁粒子在每个磁传感器位置处产生的理论磁矢量场定义为
[0020]构建目标函数：利用Levenberg-marquart算法获取磁粒子的位置矢量r＝(a,b,c)
T
和磁矩矢量m＝(m,n,p)；
[0021]A6：定义r
i
＝(a
i
,b
i
,c
i
)
T
为t
i
时刻计算得到的磁粒子的位置矢量，
[0022]m
i
＝(m
i
,n
i
,p
i
)为t
i
时刻计算得到的磁粒子的磁矩矢量，
[0023]r
i+1
＝(a
i+1
,b
i+1
,c
i+1
)
T
为t
i+1
时刻计算得到的磁粒子的位置矢量，
[0024]m
i+1
＝(m
i+1
,n
i+1
,p
i+1
)为t
i+1
时刻计算得到的磁粒子的磁矩矢量，
[0025]则磁粒子的平移速度计算公式为：
[0026]磁粒子旋转单位向量为：
[0027]A7：采用矢量图示法将计算得到的磁粒子在长方体内部的平移速度矢量和旋转速度矢量可视化显示；
[0028]其中i是自然数。
[0029]在上述技术方案中，所述磁粒子的旋转向量为三维向量，其方向与旋转轴一致，长
度等于旋转角，所述旋转角为：
[0030]在上述技术方案中，任一一个三轴磁传感器的误差函数定义为
[0031]综上所述，由于采用了上述技术方案，本专利技术的有益效果是：
[0032]本专利技术的方法可以不依赖于光学，通过磁性信息跟踪粒子，可以在完全不透明的环境中使用，仅需要三轴磁传感器阵列和高精度数据同步采集卡，试验成本较低。
附图说明
[0033]本专利技术将通过例子并参照附图的方式说明，其中：
[0034]图1基于磁粒子跟踪的内流场测速传感器安装示意图；
[0035]图2基于磁粒子跟踪的内流场测速方法时序图；
[0036]图3为第一个磁通门传感器测得的磁矢量信号；
[0037]图4为估计得到的磁粒子位置矢量信号；
[0038]图5为磁矩矢量估计结果；
[0039]图6为磁粒子平移速度；
[0040]图7为磁粒子一个旋转周期的旋转向量三维示意图；
[0041]图8为整个测试周期内的旋转角；...

【技术保护点】

【技术特征摘要】
1.一种基于磁粒子跟踪的流化床内流场测速方法，其特征在于包括以下步骤：S1:在流化床外壁上设置三轴磁传感器阵列，该阵列中三轴磁传感器的数量不少于两个；S2:制备含一个含有磁性的颗粒，使其与其他无磁性颗粒的外部形状保持一致，并将磁性颗粒与若干个无磁性颗粒混合后设置在流化床内部；S3：同步采集三轴磁传感器阵列的信号，测量得到磁粒子在三轴磁传感器位置处产生的磁矢量场；S4：利用磁粒子在空间任一点产生的磁矢量场的理论计算公式和实际测得的磁传感器阵列的磁矢量场，建立理论值与实测值差量的二范数最小的目标函数；S5：利用非线性最小二乘算法求解得到磁粒子的位置矢量和磁矩矢量；S6:获取单个采样间隔内磁粒子的位置矢量差，即可计算出该磁粒子的瞬时平移速度，获取单个采样间隔内磁粒子的磁矩矢量旋转角度，即可计算出该磁粒子的瞬时旋转速度；所述瞬时平移速度、瞬时旋转速度即代表流化床内流场当前的局部流速。2.根据权利要求1所述的一种基于磁粒子跟踪的流化床内流场测速方法，其特征在于包括以下步骤：A1:在流化床的外壁上设置有若干个三轴磁传感器阵列，磁传感器的位置坐标为(x
k
,y
k
,z
k
)
T
，k＝1
…
N,其中：N为三轴磁传感器的数量；A2：采集无磁粒子情况下的流化床所处位置的背景磁场A3：在流化床内设置一个磁粒子；A4：开启流化床的流场控制设备使得磁粒子随流场运动，同时启动磁场采集设备，实时采集磁粒子运动时在磁传感器位置处产生的磁矢量场；A5：将每个磁传感器测得的实际磁矢量值分别记为磁粒子在每个传感器位置处产生的理论磁矢量场定义为构建目标函数：利用Levenberg-m...

【专利技术属性】
技术研发人员：尹刚，郑杰匀，魏志，刘刚，张仲志，王武，邓吉龙，廖晓林，杨振华，蒋鸿，赵莉，谢艳，乔志远，邓路军，何登，
申请(专利权)人：中国空气动力研究与发展中心高速空气动力研究所，
类型：发明
国别省市：