一种基于微纳米多尺度表面结构的流动转捩控制方法与流程

本发明涉及一种流动转捩控制方法，尤其涉及一种基于微纳米多尺度表面结构的流动转捩控制方法，属于流动控制。

背景技术：

1、长期以来流动从层流到湍流的转捩问题一直是空气动力学领域的热点问题，层流减阻技术通过维持机翼表面大范围的层流流动，能够大幅降低摩擦阻力，是提高飞行器气动效率的一项核心技术。随着新一代超声速民机研究热潮的兴起，超声速层流减阻技术开始得到广泛关注，对于超声速飞机，为了获得较低的波阻，采用大后掠翼气动构型已成为气动界的共识，研究表明，大后掠翼构型的流动不稳定机制相比传统飞机构型变得更为复杂和敏感，机翼边界层可能包含前缘线、横流c-f波、流向t-s波、görtler涡等多种不稳定机制。

2、流动转捩控制方法是指根据边界层内不稳定扰动波的时间和空间尺度，采用相应的控制方法，以抑制不稳定扰动波进一步发展，从而达到延迟转捩的目的，目前使用较多的为表面结构控制和脉冲电弧等离子体激励控制，其中表面结构控制主要根据边界层内不稳定扰动波的空间尺度，采用相应的表面结构，以达到抑制不稳定扰动波沿程放大的目的。

3、由于超声速机翼构型复杂，其边界层内的不稳定波的发生机理是多样的，因而其扰动波的空间尺度也是多样的，现有技术中传统的单一尺度的表面结构仅能针对同一种不稳定机制，无法兼顾其它尺度的不稳定扰动波，无法适用于多种转捩机制共存的情形，鲁棒性和工程实用性较差，因此，需要一种宽适用性流动转捩控制方法。

技术实现思路

1、在下文中给出了关于本发明的简要概述，以便提供关于本发明的某些方面的基本理解。应当理解，这个概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关键或重要部分，也不是意图限定本发明的范围。其目的仅仅是以简化的形式给出某些概念，以此作为稍后论述的更详细描述的前序。

2、鉴于此，为解决现有技术中传统的单一尺度的表面结构无法适用于多种转捩机制共存情形的问题，本发明提供一种基于微纳米多尺度表面结构的流动转捩控制方法。

3、技术方案如下：一种基于微纳米多尺度表面结构的流动转捩控制方法，包括以下步骤：

4、s1.根据线性稳定性理论的计算数据设定微纳米多尺度表面结构设计参数及布置准则；

5、具体的：

6、微米级结构单元阵列布置在机翼表面前缘后方，与前缘距离为机翼弦长的2%-5%；

7、微米级结构单元为凸起半球体结构，其直径为0.2-0.3，单位为mm，微米级结构单元凸起高度为边界层厚度的40%-60%；

8、微米级结构单元沿平行于机翼前缘方向等间距排列，微米级结构单元展向间距为最不稳定低频扰动的展向波长的60%-80%；

9、纳米级结构单元为直径不等的凸起球体结构，其直径为20-50，单位为nm，纳米级结构单元凸起高度为其直径的一半；

10、纳米级结构单元环绕式布置在微米级结构单元上，纳米级结构单元间距为其直径的4倍-5倍；

11、s2.利用cfd结构网格求解器，计算出微纳米多尺度表面结构布置所需的流动参数；

12、具体的：

13、s21.建立直角坐标系下的n-s方程以及理想气体状态方程，得到流场速度信息；

14、s22.根据流场速度信息，结合边界层厚度的数学定义得到距离机翼前缘后方2%-5%处的边界层厚度；

15、s3.基于线性稳定性分析，计算出微纳米多尺度表面结构布置所需的最不稳定低频扰动的展向波长；

16、具体的：

17、s31.对可压缩边界层内流场进行参数的无量纲化，构建瞬时流场，给定方向波数和扰动频率的扰动波，通过求解特征值得到方向波数；

18、s32.构造多种扰动波，给定不同数值的方向波数和扰动频率，根据得到的最大扰动增长率σ，计算出最不稳定低频扰动的展向波长；

19、s4.根据计算出的边界层厚度和最不稳定低频扰动的展向波长，按照微纳米多尺度表面结构布置准则进行现场施工对流动转捩进行控制。

20、进一步地，所述s21中，n-s方程表示为：

21、；

22、；

23、；

24、；

25、；

26、其中，、、分别为笛卡尔坐标系下的流向、法向和展向，为密度，为流向速度，为法向速度，为展向速度，为温度，为粘性系数，为气体压强，为体膨张粘性系数，为耗散函数，为热传导系数，为定压比热容，为时间变量；

27、理想气体状态方程表示为：

28、；

29、其中，为气体常数；

30、耗散函数表达式为：

31、；

32、根据stokes假设，设，求解出基本流场速度信息，即基本流场的流向速度、基本流场的法向速度和基本流场的展向速度的值；

33、所述s22中，边界层厚度的数学定义为：边界层内从壁面开始，沿着壁面切向的流动速度达到自由来流速度u的99%的位置，即流向速度等于0.99u时所在位置垂直于壁面的高度为边界层厚度。

34、进一步地，所述s31中，对可压缩边界层内流场进行参数的无量纲化，采用长度l作为参考量，，其中，为边界层外缘的运动粘性系数，为边界层外缘的自由来流速度，采用边界层外缘的自由来流速度、边界层外缘的温度、边界层外缘的密度、边界层外缘的粘性系数、边界层外缘的热传导系数和边界层外缘的动压作为参考量；

35、引入平行流假设，忽略基本流垂直速度分量，在基本流上叠加三维小扰动，得到瞬时流场；

36、瞬时流场表示为：

37、；

38、其中，为基本流场信息，即通过n-s方程求解得到（），为基本流场的密度，为基本流场的温度，为基本流场的粘性系数，为基本流场的热传导系数；

39、针对三维超声速转捩问题，三维小扰动表示为：

40、；

41、其中，为扰动的特征向量，， t为矩阵的转置符号；

42、整理上述方程得到八阶的常微分方程组表示为：

43、；

44、；

45、；

46、；

47、；

48、其中，，为边界层外缘的特征雷诺数，，为虚数，为比热比，为边界层外缘马赫数的平方，为普朗特数；

49、根据给定方向波数和扰动频率的扰动波，通过求解特征值得到方向波数；

50、方向波数表示为：

51、；

52、其中，为方向波数的实部，表示扰动波的展向波数，为方向波数的虚部，表示扰动波的展向增长率；

53、扰动频率表示为：

54、；

55、方向波数表示为：

56、；

57、其中，为方向波数的实部表示扰动波的流向波数，为方向波数的虚部，表示扰动波的流向增长率；

58、所述s32中，构造多种扰动波，给定不同数值的方向波数和扰动频率，计算出最大扰动增长率σ；

59、最大扰动增长率σ表示为：

60、；

61、根据最大扰动增长率σ，得到最不稳定低频扰动的展向波数，计算出最不稳定低频扰动的展向波长；

62、最不稳定低频扰动的展向波长表示为：

63、。

64、本发明的有益效果如下：本发明可以根据风洞来流条件以及机翼参数计算出边界层厚度和最不稳定低频扰动的展向波长对机翼进行微纳米多尺度表面结构布置，微米级结构利用低频扰动模态发展的竞争机制，通过引入合适展向波长的扰动，来抑制最不稳定模态的增长，重点针对展向波长尺度较大的低频扰动；纳米级结构通过调控高频扰动波的偏振、相位、振幅、频率等特性来抑制或破坏其增长和传播，重点针对展向波长尺度较小的高频扰动；本发明改善了传统单尺度表面结构控制模式无法适于多种转捩机制共存的情形，大幅扩展了表面结构控制的鲁棒性和工程实用性，且本发明无需额外的系统装置，也无需主动消耗能量。