一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面，设计了局部化的数据组织方式，节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法，该算法仅需少量的逻辑运算与代数运算，不仅可以快速跟踪模拟分子在网格之间的迁移，而且可以准确判别分子与物面是否相互作用，搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置，避免了重新计算。引入碰撞距离的思想，既减少存储又保证正确的模拟结果。在程序编制过程中，我们充分展现了Fortran90高级语言的主要特性，引入动态数组、指针、链表以及派生类型数据，编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

二维流热耦合数值模拟程序及方法研究

本文发展了一套二维流热耦合数值模拟程序,流体域中采用结构化网格,固体域中采用非结构化网格,流体域和固体域耦合面上的网格为非匹配连接。为了使该程序能够同时模拟低速不可压流动,程序中还加入了预处理方法。此外,本文还发展了一种能较好保持通量守恒的界面插值技术,来处理流体域和固体域耦合面上的非匹配网格。数值计算结果与实验值或商业软件结果的对比表明:本程序中的预处理方法有效地改善了低速不可压流动下的计算收敛性和收敛精度,使本程序能够同时模拟高速和低速不可压流动;扩压器和叶栅的流热耦合计算表明本文程序在流热耦合计算方面具有较好的模拟能力,同时也表明本文发展的界面插值技术的可靠性。     

子母弹抛壳过程非定常流动的数值模拟

利用弹簧近似和网格重构相结合的非结构动网格技术,耦合求解Euler方程及弹道方程,模拟飞行在46千米高空、0°攻角、马赫数6.5、壳片分离初速分别为6m/s,12m/s、壳片分离角速度分别为0°/s,50°/s,150°/s的子母弹抛壳动态过程。研究表明,抛壳初期,壳片与母弹周围形成复杂流动,在壳片内表面存在激波的生成、发展、减弱等非定常过程,导致壳片受到气动力随着分离距离增加表现为很强的非线性特征;壳片分离初始角速度对壳片姿态的影响不是线性的,尽管分离初始角速度增加到150°/s,由于转动惯量小、俯仰力矩大,壳片实际飞行姿态达不到设计要求,壳片攻角较小,气动阻力小,导致壳片伴随母弹较长时间飞行;增大分离初速是实现快速分离的有效措施。     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

基于Riemann问题的解在自适应无结构网格上求解Euler方程

给出了一种根据 Riemann问题的解 ,计算网格单元边界处的守恒量通量的方法。该方法不受网格单元形状的限制 ,具有较好的通用性。实际计算表明 ,采用这一方法在自适应无结构网格上求解 Euler方程 ,对于复杂的流动细节具有较高的分辨率 ,并且使用多重网格法能够加速收敛     

非结构网格上Euler方程的区域分裂算法及并行计算

改进了“波阵面”区域分裂算法,并应用于流场区域的划分;对于子区域边界的不光滑现象,为尽量减少通讯消耗,提出了一种边界并行优化策略。利用PVM并行环境,探讨了非结构网格上求解Euler方程的分区并行算法。根据改进的区域分裂算法及优化策略,运用Jameson有限体积法,对二维翼型流场进行了分区并行求解,多区计算的结果与单区计算的结果作了比较,表明了本文研究方法的有效性。     

非结构网格多区域动态加密及分布式并行计算

针对分布式并行系统,提出一种二维非结构网格多区域动态并行策略,可在计算过程中对各计算区域的网格进行局部自适应网格加密,可在不过多消耗计算机资源的情况下,快速,有效地实现多台同构计算机上的计算网格的动态负载平衡.并针对二维Euler方程,利用已有的高精度MmB有限体积格式,应用提出的并行计算策略,进行了多个数值试验,得到令人满意的效果.     

分级型进化算法及其在翼型跨声速气动优化中的应用

提出了变精度模型的分级型进化算法并对初始翼型为NACA0012的二维翼型进行了跨声速流动下的形状增升优化。借鉴自然进化中不同进化阶段个体适应度评估环境不同的机理，构造了分级精度模型以加速优化过程。对翼型进行了给定气动条件下的形状增升优化，给出了优化结果，并与传统基因算法及单精度模型算法结果进行了对比。     

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。     

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

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