CN111159857A - 一种音速喷嘴管壁二维瞬态温度场重构方法 - Google Patents

一种音速喷嘴管壁二维瞬态温度场重构方法 Download PDF

Info

Publication number
CN111159857A
CN111159857A CN201911277469.9A CN201911277469A CN111159857A CN 111159857 A CN111159857 A CN 111159857A CN 201911277469 A CN201911277469 A CN 201911277469A CN 111159857 A CN111159857 A CN 111159857A
Authority
CN
China
Prior art keywords
temperature field
pipe wall
temperature
theoretical
function
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
CN201911277469.9A
Other languages
English (en)
Other versions
CN111159857B (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN201911277469.9A priority Critical patent/CN111159857B/zh
Publication of CN111159857A publication Critical patent/CN111159857A/zh
Application granted granted Critical
Publication of CN111159857B publication Critical patent/CN111159857B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

本发明涉及一种音速喷嘴管壁二维瞬态温度场重构方法，包括以下步骤：由经过优化布局后的若干温度传感器得到音速喷嘴管壁监测点的温度值；根据内外壁面结构参数和管壁内部传热特点，将管壁平面划分为m个子区域，基于克里金方程，计算子区域的温度场；对重构温度场进行交叉验证，交叉验证结果用来调整子区域的划分，使区域划分更加合理。

Description

一种音速喷嘴管壁二维瞬态温度场重构方法
技术领域
本发明属于气体流量检测领域，涉及一种音速喷嘴管壁二维瞬态温度场重构算法。
背景技术
临界流文丘里喷嘴(音速喷嘴)作为一种重要的流量测量标准，具有性能稳定、结构简单、坚固耐用、维护方便、准确度等级高等特点，因而受到了广泛的关注[1]。随着航天科技、医学研究的不断发展，高精尖领域对微小气体流量测量的需求促使各种微小气体流量计迅速发展。鉴于音速喷嘴本身的诸多优点，人们开始将其应用到微小流量测量领域。然而由于微小流量气体的特殊性、复杂性，这些微小气体流量计的精度及质量并不高[2]。故对音速喷嘴的性能及其影响因素的研究成为一项重要课题。
目前对音速喷嘴性能影响因素的研究主要有湿度[3]、入口段[4]、扩散段形状和温度[5]等，而温度直接影响音速喷嘴的准确度。由于音速喷嘴自身特殊的结构，气体在喷嘴中不断加速膨胀，并伴有较大的温降，低温气体与喷嘴固体内部结构的传热过程会对喷嘴固体内部结构和热边界层造成一系列复杂影响[6]，故音速喷嘴温度场重建对研究温度对喷嘴特性影响具有重要意义。
现阶段的温度场重建方法大多根据有限的监测点数据进行空间插值，专利[7]采用克里金插值法，但克里金插值法以最小二乘为标准，插值的结果不可避免的存在着平滑效应，即较小的值常常被夸大，而较大值往往被低估[8]。同时，温度场重建过程中往往会忽略实际物理过程。这就导致温度场重建结果出现“牛眼”、“趋中效应”等不符合实际情况的现象[9]，不能真实反映温度场的变化特性。
参考文献
[1]王超,王刚,丁红兵.小喉径音速喷嘴热效应对流量影响的热边界层分析[J].机械工程学报,2015,(16):164-170.
[2]李佳琪,李霞,戚珊珊,等.小流量音速喷嘴装置及测量特性研究[J].仪器仪表学报,2015,36(Z1):317-321.
[3]李春辉,陈宇航,王池,等.湿度对音速喷嘴流出系数影响的实验研究[J].计量技术,2009,(10):3-6.
[4]李春辉,王池.音速喷嘴入口段形状对流出系数的影响[J].计量学报,2008,29(z1):211-214.
[5]郑哈,朱云,蔡晴,等.温度对气体标定装置中喷嘴的特性影响[J].仪表技术与传感器,2015,(7):101-104.
[6]王超,王孝通,丁红兵,等.基于音速喷嘴固体内部温度分布的“热效应”分析[J].天津大学学报,2018,51(8):777-785.
[7]丁红兵，王刚，王超.管壁内部动态温度分布实时监测方法[P].中国专利：CN106052891A，2016-10-26.
[8]杨雨亭,尚松浩,李超.土壤水分空间插值的克里金平滑效应修正方法[J].水科学进展,2010,21(2):208-213.
[9]张世文,宁汇荣,高会议,等.基于各向异性的区域土壤有机碳三维模拟与空间特征分析[J].农业工程学报,2016,32(16):115-124.
发明内容
本发明在克里金插值基础上提出一种新的音速喷嘴管壁二维瞬态温度场重构方法，用以提高重构温度场的准确性，技术方案如下：
一种音速喷嘴管壁二维瞬态温度场重构方法，包括以下步骤：
(1)由经过优化布局后的若干温度传感器得到音速喷嘴管壁n个监测点si(xi,yi)，i＝1,2,…,n，的温度值Ti，对于两个监测点si，sj间x和y方向的滞后距分别为hx＝|xi–xj|和hy＝|yi–yj|。
(2)根据内外壁面结构参数和管壁内部传热特点，将管壁平面划分为m个子区域，第m个子区域的温度场计算步骤如下：
(a)求出x和y方向的实验变差函数γx *(hx)和γy *(hy)；
(b)选择高斯模型作为理论变差函数(h)
Figure BDA0002315937420000021
其中h为滞后距；c0为块金常数；c为拱高；a为变程，采用粒子群优化算法优化理论变差函数，将求解出最优的上述三个参数作为目标，得到x和y方向的理论变差函数γx(hx)和γy(hy)；
Figure BDA0002315937420000022
Figure BDA0002315937420000023
其中cx为x方向理论变差函数的拱高，ax为x方向理论变差函数的变程，cy为y方向理论变差函数的拱高，ay为y方向理论变差函数的变程；
(c)对x和y两方向的理论变差函数进行结构套合，得到套合后的理论变差函数(h)；
Figure BDA0002315937420000024
其中，
Figure BDA0002315937420000025
(d)求解下列克里金方程组得到加权系数λi，并计算预测点温度值：
Figure BDA0002315937420000026
Figure BDA0002315937420000027
其中，μ是拉格朗日乘数，(si,sj)为监测点si与sj间的理论变差函数值，(si,s0)为监测点si与预测点s0间的理论变差函数值，T*(x0,y0)为预测点温度；
(3)对重构温度场进行交叉验证，交叉验证结果用来调整子区域的划分，使区域划分更加合理。
步骤(3)执行如下的步骤：
(a)从重建温度场中选取若干个有限点Ti，与数值实验得到的温度场中相同位置处的温度值Ti′作比较，计算二者的最大误差Δemax和平均误差
Figure BDA0002315937420000031
Figure BDA0002315937420000032
Figure BDA0002315937420000033
(b)考虑温度场二维泊松方程
Figure BDA0002315937420000034
使方程右边f(x,y)值最小，即温度场内热源之和最小。
附图说明
图1为本发明实施例的管壁二维瞬态温度场重构算法流程图。
图2为本发明实施例的管壁二维瞬态温度场重构结果图。
具体实施方式
现结合附图和实例对本发明做进一步说明。
如图1所示为本实施例的管壁二维瞬态温度场重构算法流程图。以下分步骤详细介绍：
1.重构温度场需先获取温度传感器监测到的采样点温度值T(xi,yi)(i＝1,2,…,n)，并计算任意两个采样点si，sj间x和y方向的滞后距hx＝|xi–xj|，hy＝|yi–yj|。
2.考虑到内外壁面结构参数和管壁内部传热特点，为提高重构温度场的准确度，将管壁平面划分为m个子区域，求出各子区域x和y方向的实验变差函数γx *(hx)和γy *(hy)。
Figure BDA0002315937420000035
Figure BDA0002315937420000036
3.选择高斯变差函数模型作为理论变差函数。
Figure BDA0002315937420000037
其中h为滞后距；c0为块金常数；c为拱高；a为变程。现以重建其中某个子区域温度场举例，采用粒子群优化算法寻找该子区域x和y方向最优的理论变差函数，将求解出最优的上述三个参数作为目标，得到x和y方向的理论变差函数γx(hx)和γy(hy)，粒子群优化算法详细步骤如下：
(a)算法开始前需先设定所优化参数的取值范围：0<a<am，0<c<cm，c0≥0。其中am表示两点之间距离最大值，cm表示实际变差函数值最大值。还需设定粒子分量最大速度
Figure BDA0002315937420000041
与算法最大迭代次数n。
(b)在所设定好的参数取值范围内，随机生成M个值作为粒子初值，并作为当前个体最优位置。粒子的初速度在
Figure BDA0002315937420000042
之间随机选取，d为每个粒子的第d个分量。
(c)计算适应度函数的权值系数λi。
Figure BDA0002315937420000043
其中hi为滞后距；Ni是对应滞后距处的样本对数；
Figure BDA0002315937420000044
为实验变差函数的平均值。
(d)计算各个粒子的适应度函数值。
Figure BDA0002315937420000045
(e)确定单个粒子历史最优位置和全局最优位置。
Figure BDA0002315937420000046
pg k+1∈{p1 k,p2 k,...pm k}＝min{F(p1 k),F(p2 k),...F(pm k)}
其中，
Figure BDA0002315937420000047
为第i个粒子经过k次迭代后的当前最优位置，
Figure BDA0002315937420000048
为经过k次迭代后全局最优位置。
(f)更新粒子的速度和位置，得到下一次迭代时的粒子状态。
vid k+1＝ωvid k+c1ξ(pid k-xid k)+c2η(pgd k-xid k)
xid k+1＝xid k+αvid k
其中，k为当前迭代次数；ω为惯性权重因子，表示粒子维持原来速度的程度；c1＝c2＝2；ξ和η是区间[0,1]内均匀分布的随机数；为收敛因子。
(g)判断更新后的粒子分量是否还在参数设定范围内，若某一分量超出设定范围，则在该分量设定范围内重新随机取值；若粒子速度超过最大速度
Figure BDA0002315937420000049
则将速度设为最大速度
Figure BDA00023159374200000410
(h)直到迭代次数到达最大值n，算法结束。全局最优位置
Figure BDA00023159374200000411
的三个分量即为最优拟合的函数参数值。
4.将该子区域x和y两方向的理论变差函数进行结构套合，得到套合后的理论变差函数(h)。
Figure BDA0002315937420000051
其中，
Figure BDA0002315937420000052
5.利用套合后的理论变差函数求解克里金方程组并求解，得到权值系数λi。
Figure BDA0002315937420000053
其中，为拉格朗日乘数，(si,sj)为监测点si与sj间的理论变差函数值，(si,s0)为监测点si与预测点s0间的理论变差函数值，
6.计算预测点温度T*。
Figure BDA0002315937420000054
7.直到该子区域内所有预测点温度值计算完毕，得到该子区域二维温度场。
8.重复步骤3至步骤8，直到所有子区域温度场重构完毕，最终得到喷嘴管壁二维温度场。
9.对重构温度场进行交叉验证。
(a)从重建温度场中选取若干个有限点Ti，与数值实验得到的温度场中相同位置处的温度值Ti′作比较，计算二者的最大误差Δemax和平均误差
Figure BDA0002315937420000055
Figure BDA0002315937420000056
Figure BDA0002315937420000057
(b)考虑温度场二维泊松方程
Figure BDA0002315937420000058
使方程右边f(x,y)值最小，即温度场内热源之和最小。
交叉验证结果用来调整子区域的划分，使区域划分更加合理。若交叉验证结果显示重构温度场与实际情况存在较大偏差，则需对音速喷嘴管壁重新划分子区域并重建新划分区域后的温度场。达到合理划分子区域，且重建出的温度场能在一定程度上反映真实情况的目的。

Claims (2)

1.一种音速喷嘴管壁二维瞬态温度场重构方法，包括以下步骤：
(1)由经过优化布局后的若干温度传感器得到音速喷嘴管壁n个监测点si(xi,yi)，i＝1,2,…,n，的温度值Ti，对于两个监测点si，sj间x和y方向的滞后距分别为hx＝|xi–xj|和hy＝|yi–yj|。
(2)根据内外壁面结构参数和管壁内部传热特点，将管壁平面划分为m个子区域，第m个子区域的温度场计算步骤如下：
(a)求出x和y方向的实验变差函数γx *(hx)和γy *(hy)；
(b)选择高斯模型作为理论变差函数(h)
Figure FDA0002315937410000011
其中h为滞后距；c0为块金常数；c为拱高；a为变程，采用粒子群优化算法优化理论变差函数，将求解出最优的上述三个参数作为目标，得到x和y方向的理论变差函数γx(hx)和γy(hy)；
Figure FDA0002315937410000012
Figure FDA0002315937410000013
其中cx为x方向理论变差函数的拱高，ax为x方向理论变差函数的变程，cy为y方向理论变差函数的拱高，ay为y方向理论变差函数的变程；
(c)对x和y两方向的理论变差函数进行结构套合，得到套合后的理论变差函数(h)；
Figure FDA0002315937410000014
其中，
Figure FDA0002315937410000015
(d)求解下列克里金方程组得到加权系数λi，并计算预测点温度值：
Figure FDA0002315937410000016
Figure FDA0002315937410000017
其中，μ是拉格朗日乘数，(si,sj)为监测点si与sj间的理论变差函数值，(si,s0)为监测点si与预测点s0间的理论变差函数值，T*(x0,y0)为预测点温度；
(3)对重构温度场进行交叉验证，交叉验证结果用来调整子区域的划分，使区域划分更加合理。
2.根据权利要求1所述的方法，其特征在于，步骤(3)执行如下的步骤：
(a)从重建温度场中选取若干个有限点Ti，与数值实验得到的温度场中相同位置处的温度值Ti′作比较，计算二者的最大误差Δemax和平均误差
Figure FDA0002315937410000021
Figure FDA0002315937410000022
Figure FDA0002315937410000023
(b)考虑温度场二维泊松方程
1)
Figure FDA0002315937410000024
使方程右边f(x,y)值最小，即温度场内热源之和最小。
CN201911277469.9A 2019-12-13 2019-12-13 一种音速喷嘴管壁二维瞬态温度场重构方法 Active CN111159857B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911277469.9A CN111159857B (zh) 2019-12-13 2019-12-13 一种音速喷嘴管壁二维瞬态温度场重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911277469.9A CN111159857B (zh) 2019-12-13 2019-12-13 一种音速喷嘴管壁二维瞬态温度场重构方法

Publications (2)

Publication Number Publication Date
CN111159857A true CN111159857A (zh) 2020-05-15
CN111159857B CN111159857B (zh) 2024-02-13

Family

ID=70556801

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911277469.9A Active CN111159857B (zh) 2019-12-13 2019-12-13 一种音速喷嘴管壁二维瞬态温度场重构方法

Country Status (1)

Country Link
CN (1) CN111159857B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112577671A (zh) * 2020-11-27 2021-03-30 武汉工程大学 一种利用克里金法的井盖监测方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006105009A (ja) * 2004-10-04 2006-04-20 Japan Aerospace Exploration Agency 流体機械の流体振動若しくは流体騒音抑制装置
CN105928592A (zh) * 2016-04-27 2016-09-07 中国计量大学 快速检定pVTt法气体流量标准装置与检定方法
CN106052891A (zh) * 2016-07-26 2016-10-26 天津大学 管壁内部动态温度分布实时监测方法
CN106527544A (zh) * 2016-09-16 2017-03-22 天津大学 一种温度、湿度、流量、压力可控的音速喷嘴气体实验装置
CN106714336A (zh) * 2016-10-25 2017-05-24 南京邮电大学 一种基于改进克里金算法的无线传感网温度监测方法
WO2018072351A1 (zh) * 2016-10-20 2018-04-26 北京工业大学 一种基于粒子群优化算法对支持向量机的优化方法
CN109945954A (zh) * 2019-04-02 2019-06-28 上海理工大学 绝热壁面的pVTt标准容器
CN110346654A (zh) * 2019-07-19 2019-10-18 中国人民解放军国防科技大学 基于普通克里金插值的电磁频谱地图构建方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006105009A (ja) * 2004-10-04 2006-04-20 Japan Aerospace Exploration Agency 流体機械の流体振動若しくは流体騒音抑制装置
CN105928592A (zh) * 2016-04-27 2016-09-07 中国计量大学 快速检定pVTt法气体流量标准装置与检定方法
CN106052891A (zh) * 2016-07-26 2016-10-26 天津大学 管壁内部动态温度分布实时监测方法
CN106527544A (zh) * 2016-09-16 2017-03-22 天津大学 一种温度、湿度、流量、压力可控的音速喷嘴气体实验装置
WO2018072351A1 (zh) * 2016-10-20 2018-04-26 北京工业大学 一种基于粒子群优化算法对支持向量机的优化方法
CN106714336A (zh) * 2016-10-25 2017-05-24 南京邮电大学 一种基于改进克里金算法的无线传感网温度监测方法
CN109945954A (zh) * 2019-04-02 2019-06-28 上海理工大学 绝热壁面的pVTt标准容器
CN110346654A (zh) * 2019-07-19 2019-10-18 中国人民解放军国防科技大学 基于普通克里金插值的电磁频谱地图构建方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
徐志鹏;洪育仙;樊奇;谢代梁;: "多传感器pVTt法气体流量标准装置测温方案研究", 传感技术学报, no. 05 *
王超;张哲晓;丁红兵;李金霞;: "湿空气超音速凝结实验平台设计", 仪器仪表学报, no. 06 *
王超;王孝通;丁红兵;安海骄;: "基于音速喷嘴固体内部温度分布的"热效应"分析", 天津大学学报(自然科学与工程技术版), no. 08 *
高全杰;汤红军;汪朝晖;贺勇;: "基于Fluent的超音速喷嘴的数值模拟及结构优化", 制造业自动化, no. 04 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112577671A (zh) * 2020-11-27 2021-03-30 武汉工程大学 一种利用克里金法的井盖监测方法及系统

Also Published As

Publication number Publication date
CN111159857B (zh) 2024-02-13

Similar Documents

Publication Publication Date Title
EP3769091B1 (en) Wind flow sensing system for determining velocity fields of wind flow
CN112513896B (zh) 一种大气污染预测的方法
CN103430030B (zh) 用于通过远程测量径向速度确定流体移动的方法和装置
Sukhinov et al. A mathematical model of pollutant propagation in near-ground atmospheric layer of a coastal region and its software implementation
CN108062595B (zh) 基于wrf/cfd/sahde-rvm耦合的复杂地貌区域短时风能预测方法
CN107677997A (zh) 基于GLMB滤波和Gibbs采样的扩展目标跟踪方法
CN107545100A (zh) 一种venlo型温室环境的高分辨率降阶建模方法
CN114036861B (zh) 基于infoGAN的三维温度场预测方法
EP4127456A1 (en) System and method for fast wind flow measurement by lidar in a complex terrain
CN111753250B (zh) 一种一维非稳态导热反问题方法
CN115943255A (zh) 在复杂地形中借助LiDAR测量风流的湍流的系统和方法
CN111967203B (zh) 一种半解析半数值的大气边界层三维台风风场建模方法
CN111159857A (zh) 一种音速喷嘴管壁二维瞬态温度场重构方法
CN110068399A (zh) 基于径向基函数和正则化的温度场重建算法
CN115659760A (zh) 一种基于Fluent的热熔胶胶枪分析和结构优化方法
CN115470726A (zh) 一种基于深度学习的高超声速进气道流场快速预测方法
Liu et al. Numerical simulation and experimental study on a new type of variable-rate fluidic sprinkler
CN116680994B (zh) 一种基于激光雷达的气溶胶追踪和风场反演方法及系统
CN114462336A (zh) 一种核反应堆主管道冷却剂平均温度计算方法
Lu et al. A method for measuring the inclination of forgings based on an improved optimization algorithm for fitting ellipses
CN107179206A (zh) 一种横流湿式冷却塔阻力特性试验装置及测试方法
Chen et al. Effect of transducer shadowing of ultrasonic anemometers on wind velocity measurement
Cenedese et al. Lagrangian statistics and transilient matrix measurements by PTV in a convective boundary layer
CN109164416A (zh) 一种三平面五元传声器阵列的声源定位方法
Tang et al. Gas flow-rate measurement using a transit-time multi-path ultrasonic flow meter based on PSO-SVM

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