槽道湍流主动控制直接数值模拟.pdfVIP

中国力学学会学术大会2009 （CCTAM2009 ） 槽道湍流主动控制的直接数值模拟 孙振生 *，2) Cédric Larricq 任玉新 * （北京清华大学航天航空学院，北京 100084 ） 摘要：本文对可压缩槽道湍流进行了直接数值模拟，并研究了利用槽道下壁面的运动来减小槽道湍流壁面阻力的 情况。为了达到减阻的目的，需要对槽道壁面的运动进行主动控制。本文详细的研究了壁面垂向反向控制和次优控 制两种控制方案，数值计算结果表明这两种控制方案分别取得了 13％和 4％的减阻效果。仔细分析了这两种方法对 近壁区流场结构以及流场统计信息的影响。研究结果表明，这两种方案都在不同程度上减小了湍流脉动，雷诺应力 并且削弱了近壁区的涡结构；再次验证了近壁涡结构是流场阻力主要来源这一观点。分析了壁面在展向和流向的变 形结构，发现展向变形宽度与近壁条带结构宽度相当，而流向变形长度与条带结构的长度相当；这些都对设计使用 的微执行机构具有一定的指导意义。 关键词： 直接数值模拟，主动控制，槽道湍流减阻，壁面运动 引言 以减阻为目的的湍流主动控制技术因其广阔的应用前景而成为现代湍流研究的热点之一。近代 微电子机械系统（MEMS）的研制成功，特别是微传感器、微控制器和微执行机构集成的 M3 芯片的问 世，为壁湍流的主动控制提供了有力的工具。其工作原理是通过微传感器采集合适的壁面参数，交 由微控制器根据控制算法得到壁面变形率，使得微执行机构变形，削弱湍流中高摩察阻力产生机制， 达到湍流减阻的目的。但是何种壁面参数作为微控制器的输入信号，何种执行动作能够控制壁面湍 流，微传感器和微执行机构在壁面如何分布等都是亟待解决的问题。因此，对于实际的微执行机构 来说，需要考虑执行机构的形状、尺度、变形规律和阵列分布等许多参数，而这些参数对流场的影 响程度仍然未知。因此，本文不考虑微传感器和微执行机构的大小，假定在壁面可以任意变形，主 要研究控制算法对减阻效果的影响，进行减阻机理的分析，通过对壁面变形情况的分析，可以对实 际执行机构的设计提供有益的参考。 近代研究表明，近壁区的流向涡结构在湍流输运中扮演着重要角色[1,2], 近壁区准流向涡引发的 高速流体的下扫运动是湍流壁面高摩察力的重要原因，所以通过控制近壁区流向涡结构可以达到减 小阻力的目的。智能壁面正是这样一种控制手段，通过寻求有效的控制方案来达到控制流向涡的目 的。本文用直接数值模拟的方法实验了两种方案，一种是Choi等提出的壁面垂向反向控制[3]，但是 为了计算的稳定性，本文添加了粘性阻尼项；另一种是Lee等人提出的以壁面切应力作为输入信号的 次优控制[4]。仔细研究对比了这两种方案的减阻效果以及对近壁区涡结构的影响，进而对减阻机理 进行了进一步的探讨。 1 数值方法 本文所求解的控制方程为三维非定常可压缩Navier-Stokes方程，空间方向采用二阶精度的有限 体积格式，无粘通量的计算首先使用MUSCL插值得到控制体界面左右两侧的状态，然后根据Roe格式 得到界面上的无粘通量，粘性通量的计算采用二阶中心差分，时间方向上采用了双时间步LU-SGS方 法[5]，为二阶隐式格式，与显式格式相比，隐式算法具有较好的稳定性，可以使用较大的物理时间 步长，同时，LU-SGS算法采用显式扫描，避免了直接求解线性方程组带来的巨大计算开销。动网格 的计算采用TFI插值。为了提高计算效率，本文采用了基于消息传递模式MPI的并行编程，利用多个 2) E-mail: szs07@mails.tsinghua.edu.cn 中国力学学会学术大会2009 （CCTAM2009 ） CPU进行并行计算，显著提高了计算效率。 x x 本文所用的计算域和坐标系如图 1 所示，其中 ，y ，z 分别表示流向，展向和垂向，计算域在 ， 4πδ 2πδ δ y 方向的长度分别为 ， ， 为半槽道高度。流向和展向是统计均匀的，采用均匀网格，垂 直壁面方向采用余弦函