双时间隐式方法求解非定常绕流的相关问题

空间采用三阶精度的WeightedNND(WNND)格式,时间方向则采用具有二阶精度的双时间步隐式方法求解二维非定常Euler方程,数值模拟了NACA0012翼型俯仰强迫振动流场,讨论了内迭代步数、内迭代收敛判据以及物理时间步长对计算结果、计算效率与计算CPU时间的影响,并将计算结果与实验值进行了比较。     

非定常Euler方程数值计算中高精度格式应用

采用五阶精度加权紧致非线性格式(WCNS)和非定常“双时间步”方法求解非定常Euler方程,模拟NACA0012翼型强迫俯仰振动流场,研究了高精度格式应用到非定常计算时“双时间步”方法物理时间步长、子迭代收敛判据、子迭代步数以及物理时间导数离散方法对计算精度和计算效率的影响.     

采用DES模型对6∶1椭球体进行水动力数值计算

采用基于非结构化网格的有限体积法离散控制方程和湍流模式,应用SIMPLE(semi-implicit method for pressure linked equations)方法处理速度-压力耦合,对流项和非定常项分别采用二阶迎风格式和二阶精度离散,隐式选代求解方程组.基于SA和SST的DES模型,得到了不同攻角下粘...     

数值模拟类升力体外形高超声速粘性绕流

采用空间二阶精度的交替方向隐式分解的NND格式求解完全气体假定下的非定常薄层近似Navier Stokes方程 ,并采用抛物化的椭圆型方程生成复杂带翼外形的空间网格。最后给出了类升力体外形在M∞ =8 0 ,α =2 0°下的流场计算结果。     

高超声速尖双锥流动高精度数值模拟

以25°/55°尖双锥外形的高超声速低焓层流流动模拟为例,对高阶加权紧致非线性格式模拟激波/边界层干扰流动的能力进行验证和确认。空间离散采用二阶MUSCL和三阶、五阶加权紧致非线性格式,时间离散采用二阶精度双时间步方法,通量函数采用混合Roe-Rusanov,AUSMPW+,Van Leer等,对比了不同精度空间离散格式对时间、网格收敛特性和通量函数耗散特性的影响。数值模拟结果表明采用高精度空间离散格式能在较疏的网格上获得收敛解,并能消除结果对通量函数的敏感性,但收敛需要推进更久的计算时间。数值模拟结果与实验测量结果吻合良好,满足工程精度要求。     

测试投稿，大家不要管哈~

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

高超声速钝锥体俯仰阻尼导数影响因素分析

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

火箭发动机燃烧室内气体非定常流动数值模拟

应用新的ＰＩＳＯ算法对液体火箭发动机内非定常流动过程进行了数值模拟计算。算法采用一步隐式预测、两步显式校正完成每一时间层计算,而不是通常的多次迭代计算,因而大大缩短计算时间。     

带控制舵机动弹头粘性绕流数值模拟

采用空间二阶精度的交替方向隐式分解的NND格式求解完全气体假定下的非定常薄层近似Navier Stokes方程 ,并采用抛物化的椭圆型方程生成复杂带翼外形的空间网格。最后给出带控制舵机动弹头在M∞ =7 3,α =2 0°下所计算的流场等值线、表面压强分布及表面极限流线     

原始变量守恒形式控制方程的时间准确性分析

基于压力、速度和温度的原始变量为自变量的守恒形式的控制方程可应用于定常流动问题，但是在求解非定常问题，例如某一典型激波管问题时，激波后温度出现过冲现象，即使通过细化网格、提高空间格式精度或者换用其他通量格式仍不能消除，这表明误差可能来自该方法本身。采用一维Euler方程对该方法进行数值分析。分析结果表明，数值误差来自时间项。通过构造相应的双时间步方程，虚拟时间项采用原始变量，而物理时间项采用守恒变量，并在两个相邻物理时间步内作为定常问题求解，可以收敛到相应的守恒形式，消除上述误差，得到准确的非定常数值解。     

二元混压式高超声速进气道流场数值模拟

使用有限差分法,TVD二阶迎风格式,Baldwin-Lomax湍流模型,对二元混压式高超声速进气道设计状态和非设计状态的粘性流场和性能进行了数值模拟,并分析了出口反压及隔离段长度对进气道流场的影响。结果表明隔离段长高比对进气道出口参数影响较大;隔离段出口反压不能太高,否则会引起进气道不起动。     

用动网格法计算理想平板的颤振导数

提出了一种数值模拟振动的理想平板绕流场 ,并提取其气动导数的方法。时间相关的不可压N S方程采用Projetion 2格式解耦 ,关于中间速度的动量方程的时间和空间离散采用二阶半隐格式 ,压力Possion方程迭代用多层网格法加速收敛。分别计算平板作竖向强迫振动或扭转强迫振动的气动力 ,用动网格法考虑平板和气流的耦合作用。由计算得到的气动力用最小二乘法确定 8个气动导数。计算结果和理想平板的Theodorsen理论值有较好的一致性     

球锥体高超声速绕流的俯仰阻尼导数的数值计算

采用Navier Stokes方程描述球锥体外形非定常振动流场 ,在Etkin理论下给出俯仰阻尼导数的计算公式。对定常流场的计算采用ADI形式的NND格式 ,对非定常流场的计算采用四步Runge Kutta方法 ,并引入变系数残值光顺技术加速收敛。将俯仰阻尼导数计算结果与实验及工程计算进行了比较 ,并数值研究了其随振动中心的变化规律。     

根轨迹法解高阶代数方程 (3-5 阶方程的解)

本文采用根轨迹法解高阶代数方程。其基本思想是沿着相应的开环极点和根轨迹走向用迭代法求出部分根,从而降低方程阶次。方法本身对方程性态不提出任何要求,对重根和复根不需作专门处理,是一种比较简便而实用的方法。     

稀薄气体动理论中介观尺度数值模拟的加速方法

针对滑移流、早期过渡流区域采用离散速度法(discrete velocity method, DVM)求解Boltzmann方程收敛速度极慢、计算资源消耗大的难题,提出全流场耦合介观/宏观方程加速方法。在介观层面基于有限差分的DVM求解Boltzmann方程,在宏观层面基于有限体积的压力耦合方程组半隐式方法求解矩方程,充分利用纳维-斯托克斯-傅里叶/R26矩方程在低克努森数下的快速收敛特性,对介观方程进行加速。基于高阶Hermite多项式重构分布函数,完成宏观方程与介观方程之间的数据传输。仿真结果表明:全流场耦合介观/宏观方程的加速方法在滑移流、早期过渡流区域具有显著加速性能,最大降低了95.28%的计算时长;但对中、大克努森数流域,加速性能大幅度下降。     

钝体高超声速气动加热与结构热传递耦合的数值计算

气动加热与结构热传递耦合问题在航天和工程应用领域非常重要。分别采用松耦合与紧耦合方法,数值模拟了高超声速二维圆管绕流的流场与结构传热耦合的非定常过程。在紧耦合方法中,流场部分采用基于Navier Stokes方程的有限体积法,将AUSM+格式与时间方向的显式多步Runge Kutta法结合;结构传热部分采用基于二维热传导方程的Galerkin有限元法。流场与结构区通过交界面的热流和温度边界条件实现耦合。计算结果分别与实验、文献做了对比,结构内部温度变化关系以及壁面的热流分布均较好地吻合。两种耦合方式的计算结果对比表明,对于流场特征时间远小于结构传热特征时间的问题,松耦合方法计算效率高,精度与紧耦合方法接近。     

相空间有限元方法在解二维中子输运方程中的应用

本文采用相空间有限元方法求解了柱形临界多群中子输运问题。其中对于方程中的坐标变量用分片连续线性多项式作为试探函数,对于方程中的角度变量用分片连续双线性多项式作为试探函数。整个求解空间区域和角度区域分别采用三角形和矩形单元划分,然后利用迦辽金方法得到一个以网格点处角通量为未知数的线性联立代数方程组,方程组中的系数矩阵的存储采用了压缩存储技术。最后用高斯消元法解此有限元方程组,表明相空间有限元方法计算收敛性较好、计算精度高。     

前缘和转轴影响翼摇滚特性的数值模拟

采用耦合求解N-S方程和Euler刚体动力学方程的数值模拟方法研究80°后掠三角翼摇滚特性,讨论机翼前缘构型和滚转轴位置对三角翼摇滚特性的影响。N-S方程的离散采用Roe格式和含双时间步的LU-SGS方法,刚体动力学方程的离散采用二阶精度的差分方法,通过交替求解流动和运动控制方程组的耦合策略,模拟80°后掠三角翼自激摇滚的非定常过程。针对转轴安装于上表面的细长三角翼,研究前缘上削尖、下削尖和双面削尖三种构型的翼摇滚特性,分析前缘构型对摇滚振幅及分岔攻角的影响。考虑到上削尖和下削尖构型在翻转后相当于前缘构型互换但转轴位置不同,进一步考察了滚转轴位置对三角翼摇滚特性的影响。结果表明:在一定的攻角下,三种不同前缘构型的三角翼均能形成大幅自维持的摇滚现象,其中上削尖前缘的三角翼振幅最大,双面削尖前缘次之,下削尖前缘三角翼振幅最小;攻角增大到30°以后,双面削尖和上削尖前缘的三角翼会发生翻转现象,翻转180°以后维持等幅振荡。     

基于流场特性的管内检测机器人的外形优选研究

管道内检测机器人的结构形状是影响其运动学、动力学特性的关键因素,在机器人设计中占有重要地位。在合理提出一些基本假设后,用二阶迎风差分格式离散管道内检测机器人附近流场的控制方程,并用SIMPLE算法求解,得到了7种不同形状的管道内检测机器人所受的差压驱动力及其周围的流场信息,并对这几种结构的CFD特性作了比较。最后综合空化条件、驱动力、腔内有效体积和加工难度等因素,确定了上游半球形结构为设计形状。