CFD／CSD紧耦合及新型动网格方法在气动弹性模拟中的应用

研发出一套基于紧耦合的CFD／CSD耦合方法和程序。非定常流场求解采用混合网格有限体积方法,时间离散采用基于LU-SGS隐式格式的双时间步长法。通过求解雷诺平均Navier—Stokes方程模拟了三维机翼的跨音速气动弹性现象。得到了其颤振边界,与风洞实验结果吻合较好,验证了方法和程序的有效性和实用意义。同时将Delau...     

多重网格法在非结构网格中的应用

将多重网格法运用于非结构网格.网格是通过聚合法得到的,网格之间是相互关联的.方程的求解采用Jamson的有限体积法.给出了二维、三维情况的数值算例.     

耦合辐射的三维热化学非平衡流数值模拟

发展耦合辐射的三维热化学非平衡流场计算方法,可用于非结构网格.采用Jameson有限体积法求解耦合辐射源项的三维N-S方程.辐射源项通过求解辐射输运方程(Radiative Transport Equation RTE)获得.在空间方向上离散后,采用有限体积法求解辐射输运方程.化学模型包含11个组元,20个化学反应.采用该数值方法计算MUSES-C模型在速度为11.6 km·s^-1时的绕流流场及前驻点处的辐射热流密度.并通过对比,分析热辐射对流场的影响.     

二维自适应非结构网格DSMC并行算法研究

研究了二维自适应非结构网格DSMC并行算法实现的过程.首先提出了一类非结构网格自适应策略,有效降低了网格尺度对计算结果的影响,提高了流场的分辨率;然后基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,利用分区并行思想,设计了非结构网格DSMC并行算法,节约了计算时间.利用For-tran90的动态分配内存技术编制了通用计算程序;最后对过渡流域高超声绕流进行了数值模拟,计算结果初步验证了算法的可行性与有效性.     

轴对称再入舱模型气动热特性数值模拟研究

应用基于块结构网格的有限体积求解方法,对热化学非平衡环境下轴对称再入舱模型的气动热特性进行了数值模拟。控制方程为带化学反应的多组元轴对称N-S方程,空间离散采用VanLeer迎风格式,时间推进为隐式LU-SGS格式;采用7组元7化学反应模型及Park双温模型模拟再入流场的热化学非平衡效应。对Hollis MP-1模型的气动热特性进行了数值模拟,分别就网格效应、湍流模型、流场的热力学性质对流场的气动力、热环境的影响进行了深入研究。研究结果表明:SST模型与k-w1998模型能更准确地计算再入流场热流峰值的位置与大小;在再入舱模型的局部区域,采用热力学非平衡模型计算的物面压强与热流结果要明显低于热力学平衡模型的结果。     

EGO方法的训练算法及应用

针对高效全局优化(Efficient Global Optimization,简称EGO)方法的训练问题,选择元启发式(Meta-heuristic)算法、随机取样算法以及低频序列算法,并选用三个无约束、两个带约束解析优化算例以及两个气动优化算例,对这三类训练算法进行详细地比较研究,发现在元启发式算法中差分进化算法最具应用潜力,而低频序列算法可以有效降低EGO方法的随机性,其中Faure序列平均性能最优.     

进化算法与确定性算法在优化控制问题中的收敛性对比

对比了进化算法(基因算法)与确定性算法(共轭梯度法)在优化控制问题中的优化效率．两种方法都与分散武优化策略-Nash对策进行了结合，并成功地应用于优化控制问题。计算模型采用绕NACA0012翼型的位流流场．区域分裂技术的引用使得全局流场被分裂为多个带有重叠区的子流场，使用4种不同的方法进行当地流场解的耦合，这些算法可以通过当地的流场解求得全局流场解。数值计算结果的对比表明．进化算法可以得到与共轭梯度法相同的计算结果．并且进化算法的不依赖梯度信息的特性使其在复杂问题及非线性问题中具有广泛的应用前景。     

用混合网格及各向异性多重网格法求解三维可压紊流流动

采用混合网格求解紊流Navier Stokes方程。在物面附近采用柱状网格 ,其他区域则采用完全非结构网格。方程的求解采用Jamson的有限体积法 ,紊流模型采用两层Baldwin Lomax代数紊流模型。用各向异性多重网格法来加速解的收敛。数值算例表明 ,用混合网格及各向异性多重网格求解紊流流动是非常有效的     

三维非结构聚合多重网格法数值模拟研究

在三维非结构网格上应用聚合式多重网格技术来加速Euler方程的收敛过程．自行设计了一种高效率的网格聚合方法．采用四重三维非结构网格，在每一层网格上采用有限体积法进行计算．通过对M6翼型的数值求解验证了多重网格加速收敛的高效性．     

非结构网格热化学非平衡辐射流场数值计算

针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法. 根据原子分子光谱理论逐条计算了100$\sim$1\,500\,nm间N, O, N$^{ + }$, O$^{ + }$的谱线以及N计算流体力学; 辐射; 热化学非平衡; 非结构网格; 有限体积法针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法.根据原子分子光谱理论逐条计算了100～1 500nm间N,O,N+,O+的谱线以及N2,O2,NO,N+2等分子的10个谱带,特别分析了NO的β'带,γ'带,δ带和ε带的辐射特性.采用耦合辐射的双温模型计算热化学非平衡流场,辐射源项通过直接求解辐射输运方程RTE(radiative transport equation)获得.在空间和方向上分别离散后,利用有限体积法求解不同方向上的辐射输运方程.计算得出了再入飞行器前驻点的辐射强度分布.采用该数值方法计算了MUSES-C模型在速度为11.6km/s时的绕流流场及前驻点处的辐射热流密度.并通过对比分析了热辐射对流场的影响.