全隐LU-SGS算法在高超声速热化学非平衡流刚性问题中的应用

在高超声速热化学非平衡流动计算中,当地气体能量松弛时间、化学反应特征时间与流动时间推进步长量级差异过大会带来严重数值刚性问题,且在高雷诺数条件下,壁面、拐角等强干扰区网格加密使得该问题加剧,导致初始最大CFL数极小,收敛速度缓慢。原始LU-SGS算法仅考虑化学反应源项和对流项的隐式处理,通过推导黏性项Jacobian矩阵谱半径并采用对角近似处理,发展了热化学非平衡FLU-SGS和BLU-SGS两种全隐LU-SGS算法;针对高焓二维圆柱和轴对称返回舱算例,对比改进前后三种算法的收敛特性。结果表明,FLU-SGS及BLU-SGS算法能够快速建立强黏性干扰和大分离流场、解决热化学非平衡复杂流计算中的刚性问题,实现初始最大CFL数3至5个量级的提升,加速收敛效果明显。

二维非结构网格Euler方程的LU-SGS算法

研究了将二维非结构三角形网格重新排序后,用LU SGS隐式算法计算Euler方程的方法。分析了网格不平衡情况下该方法的可行性和应用效果。计算了不同翼型在不同流动条件下的二维Euler方程,并与四步Runge Kutta显式方法进行了比较,结果验证了此方法具有良好的计算和收敛效果,以及在粘性计算方面的潜力。运用此方法对某预研型号的头部外形进行了设计和选择。该方法克服了以往隐式方法大量耗费内存的弱点,达到了计算耗时短和占用内存少的统一。

非结构网格上浅水方程的LU-SGS隐式算法

针对浅水数值模拟的特点,建立了可以在非结构网格上求解浅水方程的LU-SGS隐式算法.对复杂地形的底坡项进行了修正,并利用淹没节点法处理移动边界问题.利用算例对此隐式算法进行了模拟验证,并与Roe显式算法进行了比较.结果表明,隐式算法可以有效地缩短计算时间,具有良好的计算和收敛效果.长江口实际潮流模拟结果也证明了此隐式算法处理实际浅水问题的能力.

高超声速流动计算中LU-SGS隐式算法的应用

在高超声速条件下,对原始LU-SGS格式及其改进方法的收敛速度做了深入的比较分析,目的是进一步更好地将LU-SGS算法用于工程上复杂外形的计算模拟中。二维圆柱,三维钝锥及空天飞机算例的结果表明:对于高超声速粘性流动的计算,粘性项应进行隐式处理;BLU-SGS方法给出的内迭代方式的收敛性优于DP-LUR方法所给出的内迭代方式;LU-SGS算法中雅克比系数矩阵的计算方式对计算量及收敛性影响较大,若采用精确的矩阵形式,则在流动无分离情况下能取得快速收敛的效果,而在含有流动分离的情况,因受稳定性的影响,精确的矩阵形式的收敛表现不及对角近似形式。

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

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

求解Euler方程的隐式无网格算法

研究了求解Euler方程的隐式无网格算法.用点云离散计算区域,代替通常的网格划分;在当地点云上,引入二次平方极小曲面逼近计算空间导数;用Roe的近似Riemann解确定通量;并用LU SGS算法求解离散得到的Euler方程隐式时间后差联立方程组,数值模拟了二维翼型跨音速绕流.由于无网格算法区域离散只涉及点云,具有灵活性,适合处理复杂的气动外形.

欧拉方程的隐式间断有限元算法研究

针对Euler方程,设计了适合间断Galerkin有限元方法的LU-SGS、GMRES以及修正LU-SGS隐式算法。采用Roe通量以及Van Albada限制器技术实现了经典LU-SGS、GMRES算法,引入高阶项误差补偿,发展了修正LU-SGS算法。以NACA0012、RAE2822翼型为例验证分析了算法的可靠性和高效性。结果表明修正LU-SGS算法存储量较少,程序实现方便,而且计算效率是LU-SGS算法的2.5倍以上,接近于循环GMRES算法。

非结构隐式方法在空心弹丸流场模拟中的应用

采用非结构网格的LU-SGS隐式算法计算三维Euler方程,数值模拟了不同马赫数以及不同攻角下某空心弹丸绕流流场,分析了流场的波系结构及其升阻力特性,计算结果表明空心弹丸的阻力系数比同口径的普通弹丸的阻力系数大约小30％,空心弹丸的阻力系数以及升力系数随攻角的变化规律与普通弹丸一致．