用大涡模拟计算流致振动的流体激励力 large eddy simulation of turbulent buffet forces in flow induced vibration.pdfVIP

第28卷第5期 核动力工程 Vjl．28．No．5 20O 7年lO月 NuclearPower Oct．2007 En西neering 文章编号：0258-0926(2007)05．0014．04 用大涡模拟计算流致振动的流体激励力 席志德，陈炳德，李朋洲 (中国棱动力研究设计院反应堆结构力学实验室，成都．6l004” 摘要：本文对湍流大涡模拟理论的脉动过滤和sma90一ns时亚格子模型进行了介绍，并用大涡模拟对秦 山II期反应堆l：5模型压力容器和吊篮的环腔内流场进行了数值模拟，计算了吊篮表面压力的时间分布和 空间分布，其压力脉动的功率谱密度计算值和实验值在同一个量级内。因此，可以用计算所得的压力作为吊 篮振动的激励载荷。 关键词：流致振动；流体激励力；大涡模拟；cFx 2 中图分类号；0351 文献标识码：A 1前言 模拟模型进行了介绍，并以秦山Ⅱ期反应堆l：5 反应堆中存在着许多由于流体激发而导致的 模型为算例，对压力容器和吊篮的环腔内流场进 振动，其研究内容主要包括振动激励源、结构的 行了数值模拟，计算了作用在吊篮表面的压力， 动态特性、结构的响应等。当反应堆结构受到一 这些压力将作为流体激振力来计算结构的响应。 回路冷却剂的冲刷时，产生振动的激励源有：湍 流、旋涡脱落、流弹性失稳和泵压力波动。其中 2湍流的CFD模型 主要是湍流激励。由于反应堆结构的复杂性及湍 湍流是一种涡旋运动。在高雷诺数情况下， 流本身的复杂性给理论研究带来了一定的困难， 这种涡旋运动占优势，涡旋的尺寸和相应的脉动 因此，人们开展了实验研究、数值研究以及实验 频率谱域很宽。大涡旋与低频脉动相联系，由粘 和数值相结合的研究方法。在流体中振动的结构 性力所决定，小涡旋与高频脉动相联系。从原则 受到流体的作用力可以分成两部分来考虑：①与 上说，求解湍流问题并无理论上的困难。因为， 流体运动有关，与结构运动无关，称为流体激振 一方面，描述湍流运动的精确的微分方程已经得 力；②与结构运动有关的力，称为运动相关流体 出，即N．s方程；另一方面，数值计算方法的发 力。对于运动相关流体力的推导见文献【1】，而流 展，已足以直接求解N．s方程，目前已有研究机 体激振力可以通过cFD技术计算流场的压力来 构对N-S方程不作任何形式的模化和简化，利用 获得。 极细密的网格直接数值求解N．s方程，即DNs。 随着计算机及计算方法的发展。直接求解流 但是，现代计算机的储存能力和运算速度尚不足 场并把流场参数转化成结构激励载荷的方法成为 以求解任何问题。因此，人们提出了许多湍流的 可能。但是，目前受到湍流某些特性的限制，完 计算模型。对于工程中的湍流时均效应，可以用 全求解流场参数的所有尺度还存在一些困难。因 基于雷诺应力的统计方程来加以计算，但是方程 此，人们提出了许多计算湍流的模型，如雷诺应 中会出现雷诺应力和湍动扩散通量，使得方程不 力模型bs方程、大涡模拟等。雷诺应力模型抹 再封闭。众多的湍流模型就是对这两项的封闭而 平了湍流的著干细节，而直接求解N．S方程又受 提出来的，如零方程模型、单方程模型、双方程 到计算条件的限制，因而介于二者之间的大涡模 模型等。对于流体激励力，要求知道压力沿时间 拟便成了解决这些问题的有力工具。本文对大涡