CN108007668B - 一种后掠翼三维边界层摩擦阻力测量装置及测量方法 - Google Patents
一种后掠翼三维边界层摩擦阻力测量装置及测量方法 Download PDFInfo
- Publication number
- CN108007668B CN108007668B CN201711196440.9A CN201711196440A CN108007668B CN 108007668 B CN108007668 B CN 108007668B CN 201711196440 A CN201711196440 A CN 201711196440A CN 108007668 B CN108007668 B CN 108007668B
- Authority
- CN
- China
- Prior art keywords
- oil
- model
- vapor lamp
- sodium vapor
- frictional resistance
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 65
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 claims abstract description 38
- 229910052708 sodium Inorganic materials 0.000 claims abstract description 38
- 239000011734 sodium Substances 0.000 claims abstract description 38
- 238000002474 experimental method Methods 0.000 claims abstract description 31
- 229920002545 silicone oil Polymers 0.000 claims abstract description 30
- 238000009826 distribution Methods 0.000 claims abstract description 22
- 241000209140 Triticum Species 0.000 claims abstract description 10
- 235000021307 Triticum Nutrition 0.000 claims abstract description 10
- 239000012528 membrane Substances 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000007664 blowing Methods 0.000 claims description 8
- 238000000691 measurement method Methods 0.000 claims description 7
- 230000009471 action Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000010008 shearing Methods 0.000 claims description 4
- 239000010410 layer Substances 0.000 claims 4
- 239000011229 interlayer Substances 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 14
- 239000003921 oil Substances 0.000 description 62
- 238000009434 installation Methods 0.000 description 5
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 4
- 239000004744 fabric Substances 0.000 description 4
- 229910052710 silicon Inorganic materials 0.000 description 4
- 239000010703 silicon Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005305 interferometry Methods 0.000 description 2
- 238000007639 printing Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000003102 growth factor Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种油膜干涉后掠翼三维边界层摩擦阻力测量装置及测量方法，该测量包括CCD摄像机、单色钠灯光源、稳压器、反光板、油滴涂抹工具、黑色麦拉膜、计算机及试验用风洞。CCD摄像机、单色钠灯光源、反光板设在风洞试验区域外，单色钠灯对准反光板。单色钠灯光源通过稳压器与外接电源相连。CCD摄像机的信号输出端与计算机对应接口相连。麦拉膜紧密贴合在模型试验段。按照曲率大小，将模型试验段分块，调节单色钠灯光源和反光板的照射区域，使之覆盖试验段全部区域；利用油滴涂抹工具，将硅油按模型展向和弦向等比例均匀点涂在模型试验件表面的观察区域。本发明能有效的测出三维边界层表面摩擦阻力分布，并且能减少大曲率模型测量时的误差，也能根据摩擦阻力的分布推断出后掠翼三维边界层不同展向的转捩位置。
Description
技术领域
本发明涉及实验流体力学领域，本发明公开了一种后掠翼三维边界层摩擦阻力系数测量方法。
背景技术
转捩位置一直是流体力学研究中十分关心的一个问题，因为层流区和湍流区的分布直接影响到飞行器表面摩擦阻力的大小。特别是超声速飞行时，从层流边界层到为湍流边界层，超声速飞行器表面对流换热率以及粘性阻力都将提高一个量级。因此，利用风洞实验对模型表面摩擦阻力系数和转捩位置进行测量就显得十分重要。
传统油膜干涉方法往往是试验中在模型表面划一条油线，此种方法对于翼型二维边界层适用，但对于三维边界层，模型表面展向和弦向的流动情况具有较大差异，油线的方法缺少对这种差异的捕捉，不能获取模型表面全局的摩擦阻力系数，所以应用到三维边界层会有较大的误差；同样的，也缺少对于三维边界层转捩位置的判断和分析。对于具有较大曲率的模型，传统方法主要是通过增加光源来增大试验模型测量的有效区域。这种测量方法虽然可以实现单次成像的全局测量，但是试验模型表面大曲率也造成了不同位置处油膜对应的入射光折射角有很大的差别，因此需要对不同位置出的测量点都要做曲率修正，这种修正工作量繁重，且对于不同大曲率表面缺乏通用的修正公式，在工程上往往不适用，测量往往存在较大误差。
发明内容
本发明要解决的技术问题是提供一种油膜干涉后掠翼三维边界层摩擦阻力测量装置，并且提供一种油膜干涉后掠翼三维边界层摩擦阻力测量方法及其对应的转捩预测方法。
一种油膜干涉后掠翼三维边界层摩擦阻力测量设备装置，包括CCD摄像机、单色钠灯光源、稳压器、反光板、油滴涂抹工具、黑色麦拉膜、计算机及试验用风洞。CCD摄像机、单色钠灯光源、反光板设在风洞试验区域外，单色钠灯对准反光板，反光板反射的光要能照射到试验模型上，CCD相机能够对试验模型做有效的拍摄。单色钠灯光源通过稳压器与外接电源相连。CCD摄像机的信号输出端与计算机对应接口相连。油滴涂抹工具用细毛笔或者注射针头。麦拉膜紧密贴合在模型试验段。
一种油膜干涉后掠翼三维边界层摩擦阻力测量方法及其对应的转捩预测方法，包括下列步骤：
1.试验前，将模型试验件安装在试验用风洞中，并将黑色麦拉膜紧密贴合在模型试验件表面；
2.安装，固定CCD摄像机、单色钠灯光源和反光板在风洞试验区域外的位置；
3.根据后掠翼表面曲率的分布特点，将模型试验件按照弦向位置分成若干个区域，使划分的各个测量区域表面的曲率近似为零；同时调节单色钠灯光源、反光板、CCD摄像机的摆放位置和角度，使对于划分后的测量区域里的油膜，入射光的入射角都为一个固定的值30°，然后固定单色钠灯光源、反光板、CCD摄像机的相对位置和角度，改变测量区域，最终获得全部测量区域内的油膜图像；
4.利用油滴涂抹工具，将硅油按模型展向和弦向等间距均匀点涂在模型试验件表面的测量区域，油滴形状应尽量为圆形，直径为油滴间距的七分之一，使油滴能覆盖模型试验件表面的全部测量区域，并且通过展向和弦向的等间距布置，便于分析展向和弦向的摩擦阻力变化趋势，进一步分析转捩的三维特性；
5.涂抹完成后，启动风洞控制系统，开始吹风试验，液滴状的硅油在风洞来流剪切力的作用下，在模型试验件表面连续移动，逐渐被吹成楔形，在单色钠灯光源的照射下，利用CCD摄像机可以记录下明暗相间的干涉条纹，并将摄取的信息送给计算机存储起来，记录二级条纹出现后的帧数时刻以及条纹移动明显时的帧数时刻；
6.根据已知的油膜干涉条纹光学流动算法公式，
式中n0表示硅油的折射率，θ表示折射角，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ表示试验中选用的单色钠灯的波长，q∞(t)表示不同时刻的来流动压，μ(t)表示试验所用硅油不同时刻的粘性，t1和t2分别代表油膜图案拍摄的二级条纹出现后的帧数时刻和条纹移动明显时的帧数时刻。根据本发明涉及试验的特点，可以对上述算法进行简化，本发明涉及试验为低速、定常试验，则来流动压q∞(t)和硅油粘度μ(t)随时间的变化可以忽略，则公式可以简化为：
其中Δt表示计算选用的两帧照片拍摄的时间差。通过读取油膜条纹图像的以上参数，带入上述简化公式得到模型表面观察区域表面的摩擦力分布。
7.综合所测得的模型表面全部的摩擦阻力系数分布，进一步分析出模型表面三维边界层各个弦向和展向的摩擦阻力系数分布；同时通过湍流边界层摩擦阻力系数明显大于层流边界层的特点，判断出模型表面各个展向截面的转捩位置，分析各个展向转捩位置的差异。
本发明有益效果是，
效果一：改进了传统的油膜干涉法通过测表面摩擦阻力的方法。本方法中沿着展向和弦向等间距涂抹油滴，使得分布的油点，通过规律性的排列，最大程度上代表了后掠翼三维边界层表面上全局的摩擦阻力特性，通过对各点位置处摩擦阻力的测量，即可获得模型表面全局的摩擦阻力系数特性，使得本测试方法适用于三维边界层表面摩擦阻力的测量。
效果二：改进了传统油膜干涉法针对大曲率试验模型测量数据误差大的情况，发展了分区域测量的方法，通过分区域测量的手段，改善了测量精度，避免了传统方法的繁重修正计算，通过分区域测量、固定入射角的手段简化了试验操作，提高了测量精度。
效果三：利用试验模型表面三维边界层摩擦阻力的分布情况，可以判断出各个展向截面转捩位置的不同，分析其中的差异。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明测量装置布置图。
图2是本发明硅油点涂分布示意图。
图3是本发明大曲率模型分区域布光图。
图4是本发明模型试验件实物及在风洞中的安装位置。
图5是本发明后掠翼模型试验件翼型图。
图6是本发明模型试验件示意图。
图7是本发明三个布光区域油膜图。
图8是本发明初次布光区域硅油点涂分布图
图9是本发明中摩擦阻力系数在展向和弦向的大小。
附图中
1、CCD摄像机 2、单色钠灯光源 3、稳压器
4、反光板 5、麦拉膜 6、计算机
7、风洞 8、模型试验件 81、点状硅油分布
82、初次布光区域 83、第二次布光区域 84、第三次布光区域
85、后掠翼模型试验件翼
型图
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例和附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
实施例一：
一种油膜干涉后掠翼三维边界层摩擦阻力测量方法，包括以下步骤：
第一步：如图1和图4，试验前，将模型试验件8安装在试验用风洞7中，并将黑色麦拉膜5紧密贴合在模型试验件8表面。
第二步：如图1，安装，固定CCD摄像机1、单色钠灯光源2和反光板4在风洞试验区域外的位置。
第三步：如图1和图3，根据后掠翼表面曲率的分布特点，将模型试验段按照弦向位置分成3个区域，分别为初次布光区域82，第二次布光区域83以及第三次布光区域84，使划分的各个测量区域表面的曲率近似为零。同时调节单色钠灯光源2、反光板3、CCD摄像机1的摆放位置和角度，使对于划分后的测量区域里的油膜，入射光的入射角都为一个固定的值30°。然后固定单色钠灯光源、反光板、CCD摄像机的相对位置和角度，改变测量区域，最终获得全部测量区域内的油膜图像。如图7，为来流风速20m/s,模型试验件8攻角2°工况下，三次布光区域下的油膜图案。
第四步：如图2，利用油滴涂抹工具，将硅油按模型展向和弦向等比例均匀点涂在模型试验件8表面的观察区域；油滴形状应尽量为圆形，直径为油滴间距的七分之一，使油滴能覆盖模型试验件8的测量表面的全部区域。如图8，为上述工况下，模型吹风前，初次布光区域82下的油滴图案。
第五步：涂抹完成后，通知风洞控制系统，开始吹风试验，液滴状的硅油在风洞来流剪切力的作用下，在模型试验件8表面连续移动，逐渐被吹成楔形，在单色钠灯光源2的照射下，利用CCD摄像机1可以记录下明暗相间的干涉条纹。并将摄取的信息送给计算机6存储起来，记录二级条纹出现后的帧数时刻以及条纹移动明显时的帧数时刻。
第六步：根据已知的油膜干涉条纹光学流动算法公式:
式中n0表示硅油的折射率，θ表示折射角，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ表示试验中选用的单色钠灯的波长，q∞(t)表示不同时刻的来流动压，μ(t)表示试验所用硅油不同时刻的粘性，t1和t2分别代表油膜图案拍摄的初始时刻和截至时刻。根据本发明涉及试验的特点，可以对上述算法进行简化，本发明涉及试验为低速、定常试验，则来流动压q∞(t)和硅油粘度μ(t)随时间的变化可以忽略，则公式可以简化为：
其中Δt表示计算选用的两帧照片拍摄的时间差。通过读取油膜条纹图像的以上参数，带入上述简化公式得到模型表面观察区域表面的摩擦力分布。本实施例中，硅油的折射率n0＝1.4，cosθ＝0.92，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ单色钠灯的波长λ＝589nm，来流动压q∞＝245kg/(m·s2)，硅油粘性μ＝20mpa·s，选取的前后油膜图案时间间隔Δt＝18.9s。综合上述信息，得到了来流风速20m/s,攻角2°条件下，试验模型表面三维边界层摩擦阻力系数分布，如下表：
最终可以得到后掠翼三维边界层表面摩擦阻力系数分布。
实施例二：
一种油膜干涉后掠翼三维边界层摩擦阻力测量方法，包括以下步骤：
第一步：如图1和图4，试验前，将模型试验件8安装在试验用风洞7中，并将黑色麦拉膜5紧密贴合在模型试验件8表面。
第二步：如图1，安装，固定CCD摄像机1、单色钠灯光源2和反光板4在风洞试验区域外的位置。
第三步：如图1和图3，根据后掠翼表面曲率的分布特点，将模型试验段按照弦向位置分成3个区域，分别为初次布光区域82，第二次布光区域83以及第三次布光区域84，使划分的各个测量区域表面的曲率近似为零。同时调节单色钠灯光源2、反光板3、CCD摄像机1的摆放位置和角度，使对于划分后的测量区域里的油膜，入射光的入射角都为一个固定的值30°。然后固定单色钠灯光源、反光板、CCD摄像机的相对位置和角度，改变测量区域，最终获得全部测量区域内的油膜图像。如图7，为来流风速20m/s,模型试验件8攻角2°工况下，三次布光区域下的油膜图案。
第四步：如图2，利用油滴涂抹工具，将硅油按模型展向和弦向等比例均匀点涂在模型试验件8表面的观察区域；油滴形状应尽量为圆形，直径为油滴间距的七分之一，使油滴能覆盖模型试验件8的测量表面的全部区域。如图8，为上述工况下，模型吹风前，初次布光区域82下的油滴图案。
第五步：涂抹完成后，通知风洞控制系统，开始吹风试验，液滴状的硅油在风洞来流剪切力的作用下，在模型试验件8表面连续移动，逐渐被吹成楔形，在单色钠灯光源2的照射下，利用CCD摄像机1可以记录下明暗相间的干涉条纹。并将摄取的信息送给计算机6存储起来，记录二级条纹出现后的帧数时刻以及条纹移动明显时的帧数时刻。
第六步：根据已知的油膜干涉条纹光学流动算法公式:
式中n0表示硅油的折射率，θ表示折射角，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ表示试验中选用的单色钠灯的波长，q∞(t)表示不同时刻的来流动压，μ(t)表示试验所用硅油不同时刻的粘性，t1和t2分别代表油膜图案拍摄的初始时刻和截至时刻。根据本发明涉及试验的特点，可以对上述算法进行简化，本发明涉及试验为低速、定常试验，则来流动压q∞(t)和硅油粘度μ(t)随时间的变化可以忽略，则公式可以简化为：
其中Δt表示计算选用的两帧照片拍摄的时间差。通过读取油膜条纹图像的以上参数，带入上述简化公式得到模型表面观察区域表面的摩擦力分布。本实施例中，硅油的折射率n0＝1.4，cosθ＝0.92，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ单色钠灯的波长λ＝589nm，来流动压q∞＝245kg/(m·s2)，硅油粘性μ＝20mpa·s，选取的前后油膜图案时间间隔Δt＝18.9s。综合上述信息，得到了来流风速20m/s,攻角2°条件下，试验模型表面三维边界层摩擦阻力系数分布，如下表：(W1-W5依次代表展向翼尖位置到翼根位置)
最终可以得到大曲率模型表面摩擦阻力系数分布。
实施例三：
一种油膜干涉后掠翼三维边界层摩擦阻力测量方法，包括以下步骤：
第一步：如图1和图4，试验前，将模型试验件8安装在试验用风洞7中，并将黑色麦拉膜5紧密贴合在模型试验件8表面。
第二步：如图1，安装，固定CCD摄像机1、单色钠灯光源2和反光板4在风洞试验区域外的位置。
第三步：如图1和图3，根据后掠翼表面曲率的分布特点，将模型试验段按照弦向位置分成3个区域，分别为初次布光区域82，第二次布光区域83以及第三次布光区域84，使划分的各个测量区域表面的曲率近似为零。同时调节单色钠灯光源2、反光板3、CCD摄像机1的摆放位置和角度，使对于划分后的测量区域里的油膜，入射光的入射角都为一个固定的值30°。然后固定单色钠灯光源、反光板、CCD摄像机的相对位置和角度，改变测量区域，最终获得全部测量区域内的油膜图像。如图7，为来流风速20m/s,模型试验件8攻角2°工况下，三次布光区域下的油膜图案。
第四步：如图2，利用油滴涂抹工具，将硅油按模型展向和弦向等比例均匀点涂在模型试验件8表面的观察区域；油滴形状应尽量为圆形，直径为油滴间距的七分之一，使油滴能覆盖模型试验件8的测量表面的全部区域。如图8，为上述工况下，模型吹风前，初次布光区域82下的油滴图案。
第五步：“涂抹完成后，通知风洞控制系统，开始吹风试验，液滴状的硅油在风洞来流剪切力的作用下，在模型试验件8表面连续移动，逐渐被吹成楔形，在单色钠灯光源2的照射下，利用CCD摄像机1可以记录下明暗相间的干涉条纹。并将摄取的信息送给计算机6存储起来，记录二级条纹出现后的帧数时刻以及条纹移动明显时的帧数时刻。
第六步：根据已知的油膜干涉条纹光学流动算法公式:
式中n0表示硅油的折射率，θ表示折射角，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ表示试验中选用的单色钠灯的波长，q∞(t)表示不同时刻的来流动压，μ(t)表示试验所用硅油不同时刻的粘性，t1和t2分别代表油膜图案拍摄的初始时刻和截至时刻。根据本发明涉及试验的特点，可以对上述算法进行简化，本发明涉及试验为低速、定常试验，则来流动压q∞(t)和硅油粘度μ(t)随时间的变化可以忽略，则公式可以简化为：
其中Δt表示计算选用的两帧照片拍摄的时间差。通过读取油膜条纹图像的以上参数，带入上述简化公式得到模型表面观察区域表面的摩擦力分布。本实施例中，硅油的折射率n0＝1.4，cosθ＝0.92，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ单色钠灯的波长λ＝589nm，来流动压q∞＝245kg/(m·s2)，硅油粘性μ＝20mpa·s，选取的前后油膜图案时间间隔Δt＝18.9s。综合上述信息，得到了来流风速20m/s,攻角2°条件下，试验模型表面三维边界层摩擦阻力系数分布，如下表：(W1-W5依次代表展向翼尖位置到翼根位置)
将上述数据绘制成散点图，如图9，结合图8，分析后掠翼三维边界层转捩过程中展向的变化:试验中三维边界层在各个展向截面均发生了转捩，其中翼根截面的转捩点约在0.8x/c位置处，其余截面约在0.9x/c位置处。通过分析，在翼尖方向的实验区域由于展向压力梯度的影响，分离区要比翼根方向更长。从翼根方向到翼尖方向，转捩分界线基本相同，但是湍流区条纹间距逐渐变小，说明摩擦阻力系数随着顺展向逐渐减小。实验结果与文献中得到的结论“转捩分界线呈等百分比弦向的特点”相符。由于扰动幅值增长因子翼根处更强，使得实验段靠近翼根方向的湍流区区域内具有更强的扰动，条纹间距更大，摩阻系数也更强。
实例三中利用试验模型表面三维边界层摩擦阻力的分布情况，可以判断出各个展向截面转捩位置的不同，分析了其中的差异。
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员能够理解本发明，但是本发明不仅限于具体实施方式的范围，对本技术领域的普通技术人员而言，只要各种变化只要在所附的权利要求限定和确定的本发明精神和范围内，一切利用本发明构思的发明创造均在保护之列。
Claims (1)
1.一种油膜干涉后掠翼三维边界层摩擦阻力测量方法，其特征在于，包括下列步骤：
(1)试验前，将模型试验件安装在试验用风洞中，并将黑色麦拉膜紧密贴合在模型试验件表面；
(2)安装，固定CCD摄像机、单色钠灯光源和反光板在风洞试验区域外的位置；
(3)根据后掠翼表面曲率的分布特点，将模型试验件按照弦向位置分成若干个区域，使划分的各个测量区域表面的曲率近似为零；同时调节单色钠灯光源、反光板、CCD摄像机的摆放位置和角度，使对于划分后的测量区域里的油膜，入射光的入射角都为一个固定的值30°，然后固定单色钠灯光源、反光板、CCD摄像机的相对位置和角度，改变测量区域，最终获得全部测量区域内的油膜图像；
(4)利用油滴涂抹工具，将硅油按模型展向和弦向等间距均匀点涂在模型试验件表面的测量区域，油滴形状应尽量为圆形，直径为油滴间距的七分之一，使油滴能覆盖模型试验件表面的全部测量区域，并且通过展向和弦向的等间距布置，便于分析展向和弦向的摩擦阻力变化趋势，进一步分析转捩的三维特性；
(5)涂抹完成后，启动风洞控制系统，开始吹风试验，液滴状的硅油在风洞来流剪切力的作用下，在模型试验件表面连续移动，逐渐被吹成楔形，在单色钠灯光源的照射下，利用CCD摄像机可以记录下明暗相间的干涉条纹，并将摄取的信息送给计算机存储起来，记录二级条纹出现后的帧数时刻以及条纹移动明显时的帧数时刻；
(6)根据已知的油膜干涉条纹光学流动算法公式，计算表面摩擦阻力系数Cf：
式中n0表示硅油的折射率，θ表示折射角，Δx表示N级条纹移动的总距离，N代表计算用的条纹级数，λ表示试验中选用的单色钠灯的波长，q∞(t)表示不同时刻的来流动压，μ(t)表示试验所用硅油不同时刻的粘性，t1和t2分别代表油膜图案拍摄的二级条纹出现后的帧数时刻和条纹移动明显时的帧数时刻；
(7)综合所测得的模型表面全部的摩擦阻力系数分布，进一步分析出模型表面三维边界层各个展向和弦向的摩擦阻力系数分布；同时通过湍流边界层摩擦阻力系数明显大于层流边界层的特点，判断出模型表面各个展向截面的转捩位置，分析各个展向转捩位置的差异。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711196440.9A CN108007668B (zh) | 2017-11-25 | 2017-11-25 | 一种后掠翼三维边界层摩擦阻力测量装置及测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711196440.9A CN108007668B (zh) | 2017-11-25 | 2017-11-25 | 一种后掠翼三维边界层摩擦阻力测量装置及测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108007668A CN108007668A (zh) | 2018-05-08 |
| CN108007668B true CN108007668B (zh) | 2018-09-21 |
Family
ID=62053686
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711196440.9A Active CN108007668B (zh) | 2017-11-25 | 2017-11-25 | 一种后掠翼三维边界层摩擦阻力测量装置及测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108007668B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108918078A (zh) * | 2018-05-14 | 2018-11-30 | 南京航空航天大学 | 一种高超声速油膜干涉法摩阻测量实验装置及实验方法 |
| CN112268684B (zh) * | 2020-12-14 | 2021-03-19 | 中国空气动力研究与发展中心低速空气动力研究所 | 低速风洞旋翼模型可变方位角表面压力测量系统及方法 |
| CN114858404B (zh) * | 2022-04-15 | 2024-06-25 | 西华大学 | 基于光流算法的全局表面摩擦阻力测试方法 |
| CN116952525B (zh) * | 2023-09-20 | 2023-12-01 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种用于风洞实验的翼型壁面摩阻非接触测量方法及系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1822252A1 (ru) * | 1990-01-02 | 1996-06-27 | Центральный аэрогидродинамический институт им.проф.Н.Е.Жуковского | Способ определения напряжения трения на поверхности аэродинамических моделей |
| CN103267622A (zh) * | 2013-05-13 | 2013-08-28 | 中国航空工业集团公司沈阳空气动力研究所 | 一种荧光油流摩擦力场测量装置及测量方法 |
| CN104697705A (zh) * | 2015-03-19 | 2015-06-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种油膜干涉表面摩擦应力测量方法 |
| EP2950071A1 (en) * | 2014-05-30 | 2015-12-02 | Airbus Operations S.L. | Aerodynamic pressure sensing system for an airfoil-shaped body |
| CN204988587U (zh) * | 2015-08-05 | 2016-01-20 | 北京航空航天大学 | 一种高精度非接触气动摩擦阻力测量装置 |
| CN106226199A (zh) * | 2016-09-13 | 2016-12-14 | 中海石油环保服务（天津）有限公司 | 一种油滴‑海水间界面张力测试装置及测试方法 |
-
2017
- 2017-11-25 CN CN201711196440.9A patent/CN108007668B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1822252A1 (ru) * | 1990-01-02 | 1996-06-27 | Центральный аэрогидродинамический институт им.проф.Н.Е.Жуковского | Способ определения напряжения трения на поверхности аэродинамических моделей |
| CN103267622A (zh) * | 2013-05-13 | 2013-08-28 | 中国航空工业集团公司沈阳空气动力研究所 | 一种荧光油流摩擦力场测量装置及测量方法 |
| EP2950071A1 (en) * | 2014-05-30 | 2015-12-02 | Airbus Operations S.L. | Aerodynamic pressure sensing system for an airfoil-shaped body |
| CN104697705A (zh) * | 2015-03-19 | 2015-06-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种油膜干涉表面摩擦应力测量方法 |
| CN204988587U (zh) * | 2015-08-05 | 2016-01-20 | 北京航空航天大学 | 一种高精度非接触气动摩擦阻力测量装置 |
| CN106226199A (zh) * | 2016-09-13 | 2016-12-14 | 中海石油环保服务（天津）有限公司 | 一种油滴‑海水间界面张力测试装置及测试方法 |
Non-Patent Citations (1)
| Title |
|---|
| 翼型边界层转捩热/油膜及红外测量技术的对比;董昊等;《南京航空航天大学学报》;20131215(第06期);第792-796页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108007668A (zh) | 2018-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108007668B (zh) | 一种后掠翼三维边界层摩擦阻力测量装置及测量方法 | |
| Zhang et al. | An experimental investigation on the surface water transport process over an airfoil by using a digital image projection technique | |
| Meier | Computerized background-oriented schlieren | |
| Driver | Application of oil-film interferometry skin-friction measurement to large wind tunnels | |
| Ota et al. | Improvement in spatial resolution of background-oriented schlieren technique by introducing a telecentric optical system and its application to supersonic flow | |
| Lang et al. | Boundary-layer transition on a rotor blade measured by temperature-sensitive paint, thermal imaging and image derotation | |
| Disotell et al. | Global surface pressure measurements of static and dynamic stall on a wind turbine airfoil at low Reynolds number | |
| CN104697705B (zh) | 一种油膜干涉表面摩擦应力测量方法 | |
| Lee et al. | Three-component velocity field measurements of propeller wake using a stereoscopic PIV technique | |
| Lunte et al. | Wall shear stress measurements by white-light oil-film interferometry | |
| Vinuesa et al. | Measurement of wall-shear stress | |
| Baldwin et al. | Skin friction measurements using oil film interferometry in a 3-D supersonic flowfield | |
| CN108918078A (zh) | 一种高超声速油膜干涉法摩阻测量实验装置及实验方法 | |
| Wolf et al. | Time-resolved volumetric particle tracking velocimetry of large-scale vortex structures from the reattachment region of a laminar separation bubble to the wake | |
| Weiss et al. | Dynamic-stall measurements using time-resolved pressure-sensitive paint on double-swept rotor blades | |
| Yorita et al. | Time-resolved pressure-sensitive paint measurements for cryogenic wind tunnel tests | |
| Gluzman et al. | A simplified photogrammetry procedure in oil-film interferometry for accurate skin-friction measurement over arbitrary geometries | |
| CN206378107U (zh) | 轮廓测量装置 | |
| Konrath et al. | Combined application of pressure sensitive paint and particle image velocimetry to the flow above a delta wing | |
| Garimella | Flow visualization methods and their application in electronic systems | |
| Costantini et al. | Experimental analysis of a wind-turbine rotor blade airfoil by means of temperature-sensitive paint | |
| CN109489934B (zh) | 基于剪切敏感液晶涂层的流道内壁面摩擦力场测量方法 | |
| CN112229820B (zh) | 一种染料池折射率测量方法 | |
| Husen | Skin friction measurements using luminescent oil films | |
| Kahler et al. | High resolution velocity profile measurements in turbulent boundary layers |
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 |