基于多尺度数值模拟的粗糙元对边界层影响机制及应用研究.docxVIP

基于多尺度数值模拟的粗糙元对边界层影响机制及应用研究

一、引言

1.1研究背景与意义

在航空航天领域，飞行器在高速飞行时，其表面会形成边界层，边界层的状态对飞行器的气动性能有着至关重要的影响。例如，当飞行器以高超声速飞行时，湍流边界层的热负荷与摩擦阻力可达层流边界层的3-5倍，严重影响飞行器的性能和安全。而表面的粗糙元会显著改变边界层的流动特性，进而影响飞行器的气动力、热传递以及飞行稳定性。研究粗糙元对边界层的影响，有助于深入理解高超声速流动的物理机制，为飞行器的设计和优化提供理论依据，从而提高飞行器的性能和可靠性，降低飞行风险。

在能源领域，风力发电机的叶片在运行过程中，表面的粗糙度会影响边界层的流动，进而影响叶片的气动效率和能量捕获能力。通过研究粗糙元对边界层的作用，能够为风力发电机叶片的设计和维护提供指导，提高风力发电的效率，降低能源成本，促进可再生能源的发展。此外，在其他涉及流体流动的工业领域，如船舶航行、管道输送等，了解粗糙元对边界层的影响也有助于优化设备性能，提高能源利用效率，减少运行成本。

传统的实验研究方法虽然能够直观地观察到粗糙元对边界层的影响，但受到实验条件的限制，如实验设备的精度、实验环境的可控性等，难以深入探究其内在的物理机制。而数值模拟方法能够克服这些局限性，通过精确的数学模型和计算算法，对复杂的流动现象进行详细的模拟和分析，为研究粗糙元对边界层的影响提供了有力的工具。因此，开展粗糙元对边界层影响的数值模拟研究具有重要的理论和实际意义。

1.2国内外研究现状

国外对粗糙元影响边界层的研究开展较早。Liepmann率先进行了垂直于流向的半圆柱对转捩影响的详细研究，发现当粗糙元高度与当地边界层位移厚度相当时（k0.92δ*，其中k为粗糙元高度，δ为当地边界层位移厚度），转捩会在粗糙元下游立即发生；而当k=0.76δ时，流动仍然保持层流流动。Fage在同类实验中得到了不同形状粗糙元使得转捩立即发生的临界高度。Tani、Hama和Mituisi应用半圆形丝线作为粗糙元进行实验，得出了粗糙元高度、位置和转捩发生位置之间的相互关系。Dryden根据量纲分析选取k/δ作为独立无量纲参数，对实验数据重新分析得到转捩位置Ret与k/δ之间的关系。在三维粗糙元研究方面，Khebanoff、Scrubauer和Tidstrom的实验得到了Ret与Ukk/ν之间的关系。近年来，一些研究从粗糙元对边界层内局部湍流度的影响出发，分析转捩提前的机制。

国内的相关研究也取得了一定进展。周云使用直接数值模拟方法研究了粗糙元3种外形对流场的影响，发现矩形粗糙元抑制效果最好。成江逸等人借助华中科技大学Φ0.25mMa6Ludwieg管风洞，研究了不同位置单粗糙元、多粗糙元与在不同来流情况下粗糙元对边界层转捩以及第二模态不稳定波发展的影响，发现当第二模态不稳定波经过粗糙元时，不稳定波幅值会明显减小，但其频段不稳定波的增长率会增大并最终导致转捩的结束点被提前，同时粗糙元对第二模态不稳定波传播速度无明显影响。

然而，当前研究仍存在不足。在数值模拟方面，加入的扰动频率往往单一，与自由来流情况下扰动频率呈离散式分布的实际情况不符，且数值模拟的边界条件设定与风洞实验存在差异，导致仅依靠数值模拟方法探究粗糙元对边界层转捩的影响存在较大难度。在实验研究中，对于复杂工况下粗糙元对边界层影响的系统性研究还不够完善，不同研究结果之间的对比和整合也有待加强。

1.3研究目标与内容

本研究旨在通过数值模拟，深入揭示粗糙元对边界层的影响规律及其物理机制。具体研究内容包括：首先，建立高精度的数值模拟模型，对不同类型、尺寸和布置方式的粗糙元在不同流动条件下对边界层的影响进行模拟，分析边界层的速度分布、温度分布、压力分布等参数的变化规律。其次，探究粗糙元诱导边界层转捩的物理机制，通过分析流场中的涡结构、扰动传播等现象，揭示粗糙元如何影响边界层的稳定性和转捩过程。最后，基于模拟结果，总结粗糙元对边界层影响的一般性规律，为工程应用提供理论支持和指导。

1.4研究方法与技术路线

本研究主要运用直接数值模拟（DNS）和大涡模拟（LES）方法。DNS方法能够直接求解Navier-Stokes方程，精确模拟流场中的所有尺度的湍流运动，对于研究粗糙元对边界层的精细影响具有重要意义。LES方法则通过对大尺度涡进行直接模拟，对小尺度涡采用亚网格模型进行模拟，在保证一定精度的前提下，能够大大降低计算成本，适用于模拟复杂的流动问题。

在技术路线上，首先根据研究需求，建立合理的物理模型和数值计算模型，确定计算区域、边界条件和初始条件。然后，利用DNS和LES方法进行数值模拟计算，对模拟