用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

有相对运动多体动力学系统流动计算方法若干问题讨论

对应用非结构动网格技术数值模拟有相对运动多体动力学系统的流动现象时发现的若干问题进行了讨论.应用动网格技术求解均匀流场,验证了几何守恒律方程离散形式对流场计算精度的重要性,还发现格式时间精度越高、网格运动速度越接近声速,几何守恒律的作用越明显;采用类似于可控轨迹实验的准定常计算方法和非结构动网格技术,分别对超声速气流中物体运动的动力学过程进行了数值模拟,通过流场和运动特性比较和理论分析,指出准定常计算方法本质上无法反映非定常流动特征,即使时间步长趋向零也不会趋向于真解;还讨论了数值模拟非定常流动的双时间法,分析时间推进步长对解的影响,随着时间步长减小解的误差降低,但是这一方法与动网格技术结合解决多体分离问题有难度.     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

基于黏性涡域法的二维钝体绕流模拟

针对传统计算流体力学方法需要划分网格且计算耗时长等问题，一种确定性涡方法(黏性涡域法)逐渐发展并日臻成熟。本文利用该方法计算模拟了二维钝体的流场、气动系数及斯特劳哈尔数，并与其他成熟CFD软件计算结果对比来验证模拟精度。黏性涡域法的模拟过程为生成涡量、涡量演变及气动计算。文章首先推导了该方法各过程的基本方程，然后基于Lua语言编制了圆柱和方柱绕流计算程序，最后利用Gnuplot实现流场可视化。利用黏性涡域法模拟不同雷诺数下圆柱绕流、不同无量纲频率和无量纲振幅的振动圆柱气动力时程、不同风向角下的方柱绕流，结果表明：雷诺数为2×10²～1×10⁵时，基于黏性涡域法的圆柱绕流气动计算结果与其他方法的计算结果基本一致；基于该方法的不同风向角下方柱绕流气动计算结果与其他文献结果趋势相同。     

高超声速连续/稀薄流自适应界面推进重叠网格方法及应用研究

根据高超声速飞行器在过渡流域飞行时绕流流场具有连续流与稀薄流并存的特点,提出了一种适用于混合流场自适应界面推进的多尺度重叠网格方法。用两套独立生成不同尺度的结构网格和非结构网格覆盖整个流场。根据连续/稀薄流界面失效函数的截止值在结构网格内进行自适应"挖洞"。在挖洞后的稀薄流区利用非结构网格DSMC方法进行计算;在连续流区的结构网格内采用基于N-S方程的CFD算法。两种算法进行耦合计算和求解,在连续/稀薄流分界面进行流场信息的插值、交换和处理。根据流场实时耦合计算结果,利用界面失效函数自适应改变两套网格的重叠部分。最后,通过对超声速圆柱和高超声速带扩张角圆管绕流的数值模拟,验证了方法的可行性和有效性。     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

薄片分离过程流场特性的数值仿真研究

利用弹簧近似法和局部网格重构相结合的非结构动网格方法,耦合求解Navier-Stokes方程和六自由度刚体动力学方程,数值模拟了三维薄片非定常分离过程。数值仿真结果有效反映出薄片分离过程中运动速度、轨迹、受力变化规律及流场特性的发展与变化。根据数值仿真结果对薄片分离机理、多体气动干扰作用进行了分析。分离过程中单个薄片与主体段周围形成复杂流动,并伴随着尾迹涡的运动、发展等非定常过程。气动载荷的非线性变化是薄片分离的主要原因。     

非结构动网格在三维可动边界问题中的应用研究

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法并通过耦合求解基于(Arb itrary Lagrangian-Eu lerian ALE)描述的Eu ler方程,模拟了做谐和振动的刚性无后掠机翼及柔性M6机翼的跨音速绕流,计算结果与实验结果吻合良好。     

高压燃气涡轮径向内冷叶片气热耦合的数值分析

采用气热耦合方法对采用径向内冷方式的MarkⅡ型跨声速高压燃气涡轮金属导叶进行数值模拟,通过分析叶片通道内的传热和流动过程发现叶片表面附面层内流动非常复杂,包含层流流动、转捩和湍流流动状态,所以只有使用转捩模型计算的叶片附面层内流动与实际情况相符,叶片壁面温度和换热系数分布与实验结果吻合的较好,使用其他湍流模型由于不能准确描述附面层内流动而使得计算结果误差相对较大,但是所有的湍流模型都能很好的模拟附面层以外流动.      

基于自由变形技术的分流叶片形状优化设计

基于自由变形技术(FFD)和计算流体力学(CFD)建立了离心压缩机叶片形状优化设计的方法,从而避免了对复杂的实体造型本身进行参数化,提高了优化效率。首先基于自由变形技术建立分流叶片外形的参数化表示,根据优化拉丁超立方试验设计方法构建控制点变化的样本空间,接着通过CFD数值仿真获得各样本的性能参数并建立响应面分析法(RSM)模型,以压缩比与等熵效率最大化目标构建多目标优化模型,进一步采用最小偏差法转化为单目标优化模型并进行求解,最后分析比较了形状优化对压缩机性能和流道内流体流动的影响。结果表明形状优化后压缩机的压缩比提升0.46%,等熵效率提升0.84%,同时减少了流道内低速区域,降低了流动损失。      

间断Galerkin方法计算叶型转捩流动

使用间断Galerkin方法研究叶型转捩流动,进行了大涡模拟(LES)和转捩模型的求解,对T106低压透平和Zierke-Deutsch压气机叶栅内的流动进行了计算.通过T106低压透平的计算对LES和转捩模型进行了比较,结果表明两种方法得到的压力分布和分离泡位置均与实验值吻合良好.LES得到的分离泡的轴向位置为0.145~0.165m,转捩模型得到的分离泡的轴向位置为0.150~0.156m.LES可以再现复杂的瞬时流动细节,对于深入研究流动机理很有意义,而转捩模型尽管无法获得足够的流动细节,但是也能预测边界层的分离和转捩现象,并且结果与LES时均结果相差不大,对于工程应用很有价值.使用转捩模型计算Zierke-Deutsch压气机叶栅内的流动也获得了与实验值符合的结果.      

基于图像处理思想的自适应网格技术

为了提高激波和旋涡等复杂流动现象的捕捉精度和计算效率,基于图像处理的思想,研究了非结构网格的自适应加密和稀疏算法.将流场计算网格与图像像素类比,借鉴图像处理的方法,提出了一种流场自适应区域的分割算法.该方法能自动对自适应区进行分层,中间区域网格可以达到原来多次自适应的效果,而且提高了生成网格的质量.将算法嵌入到流场求解器中,根据流场求解结果对网格进行自适应,并使用算例进行了验证.     

MLP高阶重构格式应用

在有限体积框架下,利用MLP(multi-dimensional limiting process)系列重构格式结合HLLHLLC(Harten-Lax-van Leer with contact)近似黎曼求解器,同时引入激波探测函数进一步降低MLP在光滑流动区域的数值耗散,数值模拟了超声速前台阶流动、结冰翼型的非定常流动、高超声速双楔流动和DLR F6-WB跨声速流场,研究了MLP系列格式在可压缩复杂流场中的表现.结果表明:在多维空间中,MLP格式能够在如强斜激波与网格线不重合等复杂流场数值模拟中保持严格的流场单调性;具有和传统MUSCL(monotone upstream-centered schemes for conservation laws)格式类似的计算效率,可以实现5阶,甚至更高阶重构;数值耗散更低,捕获更准确的激波位置,对航空工程数值模拟具有重要意义.     

A methodology for simulating 2D shock-induced dynamic stall at flight test-based fluctuating freestream

A comprehensive methodology for simulating 2 D dynamic stall at fluctuating freestream is proposed in this paper.2 D CFD simulation of a SC1095 airfoil exposed to a fluctuating freestream of Mach number 0.537 ± 0.205 and Reynolds number 6.1 × 10~6(based on the mean Mach number) and undergoing a 10° ± 10° pitch oscillation with a frequency of 4.25 Hz was conducted.These conditions were selected to be representative of the flow experienced by a helicopter rotor airfoil section in a real-life fast forward flight.Both constant freestream dynamic stall as well as fluctuating freestream dynamic stall simulations were conducted and compared.The methodology was carefully validated with experimental data for both transonic flow and dynamic stall under fluctuating freestream.Overall, the results suggest that the fluctuating freestream alters the dynamic stall mechanism documented for constant freestream in a major way, emphasizing that inclusion of this effect in the prediction of dynamic stall related rotor loads is imperative for rotor performance analysis and blades design.     

一种旋流式喷嘴的实验和数值研究

对一种小尺寸复杂结构的旋流式喷嘴进行了实验和数值研究，实验研究中确定了喷嘴前后压差与流量间的对应关系，数值研究中应用VOF方法对喷嘴内三维气液两相流动进行了计算，并与实验结果做了比较，二者吻合较好，而后应用该算法对喷嘴在不同结构尺寸下的流动过程进行了计算，对比分析后发现，旋流室及旋流器的结构和尺寸均会对喷嘴出口速度产生影响，但只有旋流器螺旋升角和槽道数目才会对喷嘴雾化角产生显著影响。     

涡轴发动机加速过程中的寿命延长优化控制

在涡轴发动机加速控制中,同时考虑发动机性能和寿命,优化加速控制策略,设计寿命延长控制.分析了加速过程中的气动稳定性、强度、燃烧稳定性以及功率等限制条件,采用序列二次规划法(SQP)优化算法进行加速控制中的燃油流量优化;根据燃气涡轮第1级静子冷却叶片温度最大值及其叶片型面与端壁的最大温差是影响发动机寿命的主要因素,采用两种寿命延长优化算法,一种是改变涡轮前温度的限制值,另一种是将涡轮前温度和动力涡轮转速同时作为优化目标,根据加权法进行折中.针对某涡轴发动机加载控制过程分别对两种寿命延长控制优化算法进行了仿真,结果表明两种方法均能有效降低涡轮前温度,并对动态性能影响较小,因此优化后的控制策略能有效延长发动机寿命.      

反向旋转盘腔内部流动特性

对反向旋转盘腔系统内部流动特性进行了数值模拟,计算结果与已有文献中的实验结果对比表明,采用的数学模型和计算方法是准确的.数值模拟结果表明:反向旋转盘腔内流动结构较为复杂,转速比对流动结构有重要影响.根据流动结构特点研究了旋转雷诺数和转速比对盘腔内压力分布以及转盘壁面摩擦力矩的影响,旋转雷诺数越大,盘腔内的压力越低,转盘壁面的摩擦力矩越大.反向同速旋转有助于增大盘腔内的压力,但与其他转速比相比其摩擦力矩最大.      

迎风面三维积冰过程中水膜流动的计算方法

对影响迎风面三维积冰过程的冰层表面水膜流动进行分析,建立描述其流动的数学模型,引入在流动前锋下游预设薄水膜的方法处理其在表面干区的流动推进过程,提出超过临界厚度的水膜被吹除的方法处理角状冰下游局部水膜被吹离表面的现象.通过对翼型-平板结构中翼型上的冰形对比和水膜流动规律分析,验证了该水膜流动数学模型的合理性和计算方法的可行性.研究表明:与Messinger模型相比,积冰模型中引入水膜流动后,可以较好地模拟三维明冰积冰,同时明冰极限的模拟结果也更接近试验值;冰层表面的水膜厚度约为10-5m,流动速度约为10-2m/s.