一种适合迭代求解的反馈力浸入边界法

对Goldstein提出的反馈力浸入边界法进行了新的思考，改进了其对力源项的计算，拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时，含有对速度误差的时间积分项，只能用于含时间项的Navier-Stokes （N-S）方程的求解，且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项，避免了时间相关的参数，使其不仅能适合非定常隐式时间推进，还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性，对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算，计算结果均与文献结果符合较好，表明本文改进的方法是有效的。得出的结论为：可以直接基于迭代次数进行反馈力源项的计算，改进的反馈力浸入边界法不仅可与非定常N-S方程结合，进行隐式求解，还可以与定常N-S方程结合用于定常流动的模拟，可将改进的方法运用到更多的流动问题当中。     

侧向多喷流流场数值研究

为研究侧向喷流对高超声速拦截弹喷流直接力控制性能的影响，发展了在有限体积离散方法框架下的三维非定常N—S方程，采用LU—SSOR隐式时间推进格式和残值项空间迎风型格式求解技术，对大气层内高超声速拦截弹多个姿态控制喷流发动机同时工作时的干扰流场进行了数值求解。计算给出了自由流与侧喷流、侧喷流与侧喷流之间复杂干扰流场结构图(等压云图、等马赫云图)以及喷流附近的压强分布，其结果与单一侧向喷流流场情况进行了比较和分析。结果显示，喷流前的压强场基本不受后喷流的影响，处于前喷流尾流区后的喷流的干扰效率将受到较大影响。     

火箭冲压发动机掺混流场数值方法研究

鉴于火箭冲压发动机补燃室内呈现出复杂的三维化学反应流,传统方法对这一流场的数值分析所遇到的主要困难是求解椭圆型Ｎ－Ｓ方程组极为耗时而不利于工程应用,采用先进的块隐式法求解恒温掺混流场得到了很快的收敛速度,这一尝试为快速计算化学反应流打下了很好的基础。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

用N-S方程有限体积法计算弹体绕流/底喷流

用单元中心有限体积法在物理空间曲线网格离散积分形式N-S方程，用坐标变换计算粘性通量项中的各偏导数值，用Runge-Kutta三步显式格式进行时间推进求解。采用了人工粘性、当地时间步长及隐式残值光顺。计算了旋成体大迎角涡流和带底部喷流的绕流场。计算中发现，各向异性隐式残值光顺较之各向同性隐式残值光顺能更快地加速迭代过程的收敛。计算的物面压力分布与实验基本符合，计算的旋成体大迎角涡流图谱合理。计算结果说明底部喷流对旋成体绕流及空气动力载荷有显著影响。     

超声速底部喷流干扰流场数值模拟

数值模拟了不同马赫数,不同喷流压比下的轴对称超声速底部喷流干扰流场,采用LU隐式算法进行数值求解并引入了Baldwin-Lomax代数湍流模型.采用分区网格将弹身与底部区域合为一个整体进行计算,得到了清晰的流场结构和弹体表面及底部的压力分布,并与试验结果进行了比较.数值模拟结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同, 喷流对底压分布有明显影响,进而对轴向力系数影响显著.     

基于分区的隐式求解二维不可压NS方程的并行实现

本文通过对目前常用并行算法的分析,基于区域分解的思想,给出了适用于隐式求解原始变量二维不可压NS方程压力修正方法的并行算法,并在网络异构编程PVM环境下进行了数值验证,讨论了计算结果,结果显示对任务级的并行,并行效率也同样随粒度的增加而提高.另外此方法还可推广到三维复杂流场的计算问题.     

轴对称欠膨胀羽流流场的分离系数矩阵解法和实验研究

本文对轴对称超音速羽流流场进行了数值计算和实验研究。数值计算方法是空间步进的分离系数矩阵差分法,并用粘性和无粘分裂法处理粘性问题。同时对同轴超音速喷流进行了流场观察和测量,所得纹影照片清晰地显示了尾喷流场的详细结构。计算结果在近场与文献[1]一致,在远场与本文的实验结果基本一致。计算实践表明,与中心差分格式相比,分离系数矩阵方法的稳定性好,精度高,捕捉的激波耗散小,且无数值振荡。     

复杂几何条件下伴随格林函数的数值求解

为简化气动噪声预测中复杂边界的网格处理,采用计算气动声学(CAA)方法与浸入式边界方法(IBM)相结合数值求解气动噪声预测所需的伴随格林函数.根据伴随格林函数的基本形式,文中设计了特定的圆柱声散射算例,并将数值结果与解析解对比,分析该方法的计算精度,验证了该方法在复杂几何边界条件下伴随格林函数求解中的适用性.最后应用该方法求解了考虑吊挂安装效应的锯齿型喷管喷流剪切层内的伴随格林函数.结果表明:由于喷流剪切层的散射效应,伴随格林函数在喷流内分布不均匀,最大值比最小值高出3倍以上,反映了喷流内不同区域声源对远场噪声贡献的差别,降噪设计可以重点考虑降低格林函数较大区域的声源强度以达到降低噪声的目的.      

横向喷流对鸭式导弹滚转特性影响研究

通过数值方法求解三维可压缩雷诺平均N-S方程,对导弹横向喷流的干扰流场进行了数值模拟,计算结果与实验数据吻合较好,基本验证了计算方法在横向喷流复杂流场数值模拟方面的有效性。在此基础上,对固定尾翼鸭式布局导弹亚、跨、超声速流场进行了数值模拟,并且计算分析了横向喷流对鸭式布局导弹滚转控制特性的影响,计算结果表明,横向喷流可以有效提高鸭式布局导弹的滚转控制能力。     

多重网格技术在复杂粘性流场计算中的应用及研究

流体力学控制方程的数值求解过程中,当网格加密、粘性效应强或来流马赫数低的时候,流场收敛非常缓慢。本文在多块网格、LUSGS隐式格式基础上采用多重网格技术,计算了NLR7301case#1两段翼型的低速绕流及DLRF4跨声速流场,在复杂二维、三维粘性流场求解中获得了满意的加速收敛效果。文中还比较了V循环、W循环的加速效率,为多重网格的工程应用奠下基础。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。     

冲压发动机进气道掺混段与弹体内外流场一体化数值模拟

采用二阶迎风隐式TVD格式分区求解可压缩N-S方程,数值模拟了导弹、进气道和掺混段内外流一体化通气模型的复杂流场.本文给出了掺混段出口压力与来流压力比为P/P∞=3.5,4.2和4.9情况的数值模拟结果.并分析了流场特性及激波波系结构,总压恢复系数σ和流量系数ψ.同风洞试验结果进行比较,表明计算结果是合理的.文中还对该CFD分区计算内边界处理所采用的数据结构进行了阐述.     

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

COIL化学非平衡流动的一种解耦方法研究

提出了一种计算氧碘化学激光器化学非平衡流动的解耦方法.求解二维可压缩非定常N-S方程及组分连续方程,计算氧碘化学激光器中的组分混合反应流场.空间离散格式为AUSM up格式,用四步龙格-库塔方法作显式时间推进.假设混合气体为热力学完全而热值非完全气体,化学反应模型采用有限速率反应模型.使用一种"点隐式"方法来处理化学源项的刚性问题.计算结果表明本文方法可以有效解决计算中的刚性问题,与现有隐式算法相比减少了大量的矩阵运算,公式推导简单,易于实现和改进.     

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

超燃冲压发动机三维喷流干扰流场的数值研究

用隐式有限体积法求解三维N－S方程数值研究超燃冲压发动机三维平行、垂直喷流干扰流场。计算中采用了具有当地极值递减（LED）性质的JST和SLIP中心差分格式。为了减缓计算网格数很多时对计算机内存的压力以及提高计算效率，运用了通量守恒的分区计算方法。     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

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

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