压电谐振风扇激励非定常流动特性数值研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123558 压电谐振风扇激励非定常流动特性 数值研究 谭蕾1 谭晓茗2 张靖周2 1 南京航空航天大学航天学院，南京210016 2 南京航空航天大学能源与动力学院，南京210016 Tel: 025 E-mail: tanlei@nuaa.edu.cn ： 摘要 利用动网格技术对压电谐振风扇产生的非定常流场进行了数值模拟。研究表明：不同谐振频率 驱动下的压电风扇呈现出完全不同的振型，振型符合一端固支悬臂梁的固有振动形态；不同振型下压 电风扇诱导的非定常流场均出现涡流，流场分布随着时间均呈现周期性变化规律，但是涡对的尺度、 位置、扰动范围各异，压电风扇工作在一阶谐振状态下产生的风速和扰动范围相对较优。 关键词：压电风扇；非定常流动；动网格；数值模拟 0前言 压电风扇是一种利用压电振子的谐振作用带动弹性膜片发生弯曲从而产生风力的 装置。由于压电陶瓷的逆压电效应，压电片产生周期性的收缩和膨胀，驱动弹性膜片发 生弯曲谐振，在膜片周围产生一系列的涡串。这些涡串在发展、运动的过程中与周围流 体发生干涉融合，从而聚合成一股自膜片尾端向下游输运的连续射流。虽然压电风扇的 扇面很小，但因其振动频率高、风力集中而不扩散、噪音小和耗能低等优点，同时，在 弹性膜片诱导的流动伴随涡流的特征，可以起到强化对流换热的作用，因此在电子元件 有效散热方面具有广泛的应用前景[1-4] 。 压电风扇诱导的气流流动具有非定常的流动特征，国内外研究人员针对压电风扇的 流动特征、压电风扇的影响因素和压电风扇的对流换热特性[5-9] 已经开展了一系列的研 究工作，但对压电风扇在不同谐振频率驱动下诱导的涡系发展的非定常流动特征的认识 还很不充分。本文运用数值计算的方法，结合膜片振动位移方程，用户自定义函数(UDF) 将膜片运动规律嵌入压电风扇边界，并利用动网格对压电风扇产生的非定常流场进行研 究，以揭示其在不同振型下的流动特性和相关参数影响规律。 1 计算模型 压电风扇可简化为一端固支的悬臂膜 y 片，如图1 所示。计算区域包括压电风扇及 x 外部流域。本文中设定压电风扇长度方向为 l x 坐标，厚度方向为 坐标，膜片固定端中 y 心位置为坐标原点，压电膜片的长度 1 图 压电风扇示意图 l 20mm ，厚度h 0.1mm ，振幅A=5mm ， No No 2011ZB52020 基金项目：国家自然科学基金资助项目（ ： ）；航空基金资助项目（ ： ） 外部流域计算范围为50mm×50mm 。计算区域采用非结构化网格，局部区域采用了网格 加密处理以提高计算精度和准确性。 由于压电风扇产生的风速较低，其流动过程由连续性方程和不可压缩的非稳态N-S 方程所控制。考虑压电风扇在振动过程中产生了湍流涡对，本文采用了不可压缩非稳态 基于重整化群RNG k −ε湍流模型，在近壁区仍然采用壁面函数法来进行处理。 压电风扇简化为一端固支的悬臂梁，将其定义为运动壁面，运动轨迹y (x , t) 既是时 间的函数，同时又是空间的函数，解的形式可以写成时间函数和振型函数的乘积。在已