采用快速动网格技术的时空同步流固耦合算法.docVIP

采用快速动网格技术的时空同步流固耦合算法 　　摘要：为了减少流固耦合计算时间，发展了一种时空同步流固耦合算法。在每一次耦合迭代中，首先求解RANS（Reynolds Averaged Navier Stokes）方程，然后采用本课题组所提出的快速动网格技术计算结构及流场网格节点位移，并更新流场网格，实现流场与结构振动的空间同步求解。在每一时间步，通过多次耦合迭代，使流场计算收敛，同时保证结构振动计算的收敛，实现流场与结构振动的时间同步求解。采用该算法对弹性梁流固耦合振动及Wing 445.6颤振问题进行了研究，计算结果与文献中的结果一致。与已有文献的时间同步算法相比，此算法可以减少计算时间81.2%。 　　关键词：流固耦合；动网格；时间同步；空间同步；时空同步 　　中图分类号：0327；TBl23 　　文献标志码：A 　　文章编号：1004-4523（2017）01-0041-08 　　DOI：10.16385/j.cnki.issn.1004-4523.2017.01.006 　　引言 　　结构与流体之问的耦合现象在交通运输、船舶、能源、建筑、机械制造、航天航空等工程领域普遍存在，有时会造成结构的损坏。1940年美国TacomaNarrows bridge在风中发生流固耦合振动，最终坍塌。2006年，某轴流压气机叶片发生流固耦合振动，诱发高周疲劳，造成叶片断裂。2010年，美国洛克希德一马丁公司的验证机在飞行实验中，机翼与空气之问的流固耦合效应诱发的颤振使机翼折断。因此，流固耦合已成为研究人员重点关注的问题。 　　基于CFD/CSD的紧耦合方法可以考虑结构与流场的相互影响，贴近物理实际。因此，研究结构的流固耦合问题通常采用基于CFD/CSD的紧耦合方法。李田采用基于CFD/CSD的流固耦合方法研究了横风下高速列车运行稳定性问题。毛国栋研究了流固耦合效应对建筑中膜结构风振响应的影响。王征基于CFD/CSD技术计算了压气机叶片流固耦合响应并预测了其颤振边界。流固耦合计算主要分为流场分析、结构振动分析及流场网格更新等3个部分。为了考虑边界变动对流场的影响，需要在ALE格式下采用动网格算法更新流场网格节点坐标。目前所发展的动网格算法主要有弹簧法、弹性体方法、温度体方法。及径向基函数方法。上述动网格方法的网格变形效率普遍偏低，采用这些方法更新流场网格的计算时间占流固耦合计算总时问的比重较大。时间同步流固耦合算法包含2个层次的迭代过程，即外部迭代（流场分析和结构振动分析之间的迭代）、流场分析的内部迭代。首先通过内部迭代求解流场，然后采用外部迭代耦合流场和结构振动。流体域和结构分别属于2个不同的空间区域。可以认为，上述时间同步算法在空间域上是不同步的。这种空间不同步的算法在流场计算收敛之后更新流场网格，可以减少流场网格更新的次数，从而减少流固耦合计算的总时间。 　　本课题组针对当前动网格方法网格变形效率低的问题，在弹性体方法的基础上发展了一种快速动网格技术，能够显著提高网格变形效率。将该快速动网格技术用于时间同步流固耦合研究后发现，对流场多次迭代求解会增加流场计算时间，使得流固耦合计算的总时间增加。本文在时间同步流固耦合算法的基础上，结合本课题组所提出的快速动网格技术，发展了一种时空同步流固耦合算法。采用该算法对弹性梁流固耦合振动问题及Wing445.6颤振进行了研究，计算得到的弹性梁振动的位移时间曲线与已有文献的结果一致，计算得到的Wing445.6颤振边界也与实验数据吻合，说明了本文算法的正确性。 　　1.时空同步流固耦合算法 　　流固耦合计算主要分为流场分析、结构振动分析及流场网格更新等3个部分。在问同步流固耦合算法中，每一时间步，都首先进行流场分析，通过迭代求解RANS（Reynolds Averaged NavierStokes）方程得到收敛的流场结果，然后分析结构振动，最后采用动网格算法更新流场网格，重复上述过程直到结构振动计算收敛，具体流程如图1所示。在每一时间步末，流场分析和结构振动分析同时满足收敛，即实现了流场和结构振动的时间同步求解。 　　通常情况下，在时间同步流固耦合计算的每一时间步，往往在流场多次（10次左右）迭代达到收敛之后进行一次流场网格更新。此时，网格变形的时问基本上比流场收敛的计算时间小一个量级。但是，每一次的网格变形的时间与流场计算收敛的过程中的每一次的迭代计算的时间基本相等。如果网格更新次数过多，就会导致网格变形的总时间赶上甚至超过流场计算的总时间。所以，在时间同步流固耦合算法中，通过流场计算收敛之后更新流场网格可以减少网格更新次数，进而减少流固耦合计算的总时间。然而，当采用本课题组所提出的快速动网格技术更新流场网格时，流场网格更新所需的计算时间与流场分析所需计算时间相比可以