一种求解低速跨流域流固耦合问题的离散统一气体动理学格式

运动边界及流固耦合问题是低速及高速、连续流域及稀薄流域流动中常面临的多场耦合问题。本文在原始离散统一气体动理学格式基础上，采用任意拉格朗日-欧拉方法及动网格技术，并耦合计算结构动力学，实现适用于低速流动运动边界及流固耦合问题的计算方法。在方法验证方面，进行连续流域圆柱绕流强迫振动和自由振动数值模拟，计算结果与参考文献结果吻合良好，准确捕捉圆柱在锁频区和非锁频区的流动特性。在稀薄流域范围，将方法中的连续流D2Q9离散速度模型替换成跨流域离散速度模型，对比低速圆柱绕流强迫振动算例在两种离散速度模型下的数值结果。虽主要关注低速流动问题，但发展的算法仍具备模拟高超声速运动边界问题的能力，对此进行了Ma=5.0稀薄流振动圆柱研究。低速和高速运动边界模拟能力显示了方法在考虑稀薄气体效应流固耦合问题方面的潜在应用价值。     

高空飞艇的流固耦合数值研究

发展了一套用于研究高空飞艇在来流作用下的流固耦合数值方法.利用经过验证的自编制的有限元程序来求解飞艇的非线性大变形,利用商业软件FLUENT求解飞艇的外流场,通过非动态交替的方法实现飞艇的流固耦合计算.结果表明:飞艇的柔性变形改变了飞艇周围的流场结构,飞艇表面气动载荷分布发生了较大的改变;由于飞艇表面气动载荷的作用,导致了囊体内外压差的变化,进而导致囊体内的应力随之发生了较大的变化.     

任意运动变形边界流场测量PIV算法的研究：薄膜涡激振动

提出了一种精确测量任意运动变形边界流场的粒子图像速度场PIV算法,能智能识别强流固耦合问题流场测量中任意运动变形的边界信息,然后生成贴体自适应的图像互相关计算窗口,计算获取运动变形边界附近的速度场数据。这一算法可大大提高任意运动边界附近流场的测量精度,为流固耦合问题的理论分析和数值计算验证提供可靠的高质量实验数据。针对数字合成给定流场的粒子图像序列,采用所发展的PIV算法对粒子图像区域中的运动变形边界进行了精确识别,高质量地复现了原始流场信息。最后,对低速闭式循环水洞中的钝体尾流柔性薄膜涡激振动现象进行了PIV实验测量,获得了柔性薄膜大变形运动状态下的瞬态流场特征。     

大展弦比机翼流固耦合迭代求解模拟

基于某型变掠翼气动布局概念构型,结合适应于大展弦比机翼的复合材料结构设计,构建基于数据通讯的流动数值方法与结构有限元分析相互迭代的流固耦合分析平台,讨论布局在不同飞行条件下的机翼流固特性,解决大尺度结构气动弹性的耦合作用方式及准确计算的方法问题.算例结果表明流固耦合迭代求解计算方法高效、合理,较为准确地获得了大尺度结构变形后的布局气动特性改变特征,完善了特殊布局的气动特性仿真方法.     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

考虑流固耦合效应的航空发动机风扇叶片应力数值模拟与试验测量研究

航空发动机风扇转子在高压比、高转速、高负荷的级环境中工作时,存在叶片固体域与流体域之间强烈的耦合作用。针对风扇工作中的流固耦合问题,采用基于流固耦合的数值模拟方法对风扇叶片的结构特性进行模拟,研究考虑流固耦合效应前后叶片结构特性的变化。通过风扇转子加减速试验测量叶片表面测点应力变化,并将数值模拟与试验测量结果进行了对比分析。分析结果表明：考虑流固耦合效应后叶片表面的受力情况变化较大,导致叶片表面的应力与变形分布产生较大的变化;仅考虑离心力作用的计算方法得到的应力值与试验测量值误差最大达到50%,而考虑流固耦合效应的计算值误差在10%左右;考虑叶片流固耦合效应得到的应力分布更满足实际工程应力与强度分析要求。     

压气机叶片流固耦合数值计算

采用流固耦合方法对压气机单排跨声叶片进行气动弹性数值模拟.首先对压气机流场进行三维定常求解,并与试验结果进行比对.在定常流场分析的基础上,进行流固耦合计算,采用非定常方法分别求解流体域Navier-Stokes(N-S)方程和固体域结构动力学方程.耦合过程中,流体域和固体域在物理边界上进行能量传递.计算结果表明:所采用的流固耦合数值方法用于判断叶轮机械叶片气弹稳定性是可行的,可以得到叶片的瞬态响应,从而判断叶片是否发生颤振.      

基于非结构自适应网格和嵌入边界法的残差格式研究

嵌入边界法由于在求解NS方程时能够简化网格生成问题而在计算流体领域受到越来越广泛的关注.简言之,嵌入边界法能够简化大变形和运动条件下多物理流动模拟、流固相互作用耦合问题,然而壁面边界条件的精确处理仍旧是该方法需要解决的问题.在本文工作中,为考虑壁面边界条件而在NS方程中增加了补偿项,同时采用非结构网格自适应技术保持了壁面边界条件的精度.     

高速飞行器流场和固体结构温度场一体化计算新方法的初步研究

本文探索和研究了一种不需作流／固耦合迭代的一体化计算方法，通过将气体流动和固体结构传热用统一的方程组进行描述，并用统一的方法进行求解，避免了流／固之间的耦合迭代问题。为了验证该数值方法的可靠性，本文对高焓高超声速气流绕二维不锈钢圆管的流动、气动加热以及结构热响应问题进行了计算，结果与实验符合良好。     

柔性扑翼全机气动结构耦合特性及稳定性研究

由于微型扑翼飞行器具有体积小、重量轻等特点,在飞行过程中会产生明显的柔性变形,因此在对扑翼气动特性进行数值模拟时,有必要考虑柔性变形的影响。基于结构动力学理论,发展了一种适用于扑翼全机气动结构耦合特性的数值模拟方法与一种适用于扑翼的结构动力学求解方法。方法利用哈密顿原理,对动能和应变能进行变分,进而得到扑翼动力学方程,采用结构有限元方法对运动方程进行离散并求解;采用基于嵌套网格的CFD求解器对微型扑翼全机非定常绕流进行数值模拟。采用CFD/CSD耦合求解器对微型扑翼飞行器全机气动结构耦合特性进行数值模拟,分析了结构变形对柔性扑翼气动特性的影响,并分别对扑动过程中刚性扑翼和柔性扑翼的压心变化范围进行研究,分析了结构变形对柔性扑翼稳定性的影响。