【摘 要】
：
自由界面流动是自然界中常见的流动现象,传统模拟方法一般都是基于网格划分,操作复杂且计算中容易出现网格扭曲畸变等问题。无网格的光滑粒子流体动力学方法(SPH方法)能够有效克服传统的网格方法在处理类似大变形问题时存在网格畸变或变形界面难以追踪等困难,是计算自由界面流动有效方法。本文即利用该方法对自由界面流动中常见的溃坝、液滴入水、固体入水等问题进行了数值模拟,模拟结果表明该方法简单有效,在处理此类问题
【机 构】
：
中国空气动力研究与发展中心,四川绵阳 621000
论文部分内容阅读
自由界面流动是自然界中常见的流动现象,传统模拟方法一般都是基于网格划分,操作复杂且计算中容易出现网格扭曲畸变等问题。无网格的光滑粒子流体动力学方法(SPH方法)能够有效克服传统的网格方法在处理类似大变形问题时存在网格畸变或变形界面难以追踪等困难,是计算自由界面流动有效方法。本文即利用该方法对自由界面流动中常见的溃坝、液滴入水、固体入水等问题进行了数值模拟,模拟结果表明该方法简单有效,在处理此类问题时具有十分优越的特性。
其他文献
本文是张涵信关于定常旋涡运动的拓扑结构理论的推广。采用微分方程奇点理论,定性研究了非定常旋涡破裂的拓扑结构。研究发现,在垂直于涡轴的截面上,其截面流线形态依赖于参数：λ(x,t)=1/ρ((e)ρ/(e)t+(e)ρu/(e)x)。如果λ(x,t)>0,截面流线向里转,如果λ(x,t)<0,截面流线向外转。如果沿着涡轴λ(x,t)变号,则截面流线出现极限环。在过涡轴的截面上。采用五阶WENO格式,
本文发展了一套新的大涡模拟显式滤波算法,将数值稳定滤波与湍流大涡模拟计算的显式滤波统一起来,通过一次滤波同时达到中心型格式计算的数值稳定要求和大涡模拟显式滤波对数值误差控制的要求。新滤波算法既避免了数值稳定滤波引入的滤波累积效应,又克服了通常采用的、对非线性项滤波的显式滤波方案破坏控制方程的伽利略变换不变性的问题。
本文基于Menter SST模型的RANS/LES混合模式对二维拐角流动进行了数值模拟研究,空间无粘项采用五阶WENO格式,粘性项采用二阶中心差分格式,时间项采用三步三阶Runge-Kutta方法。通过计算给出了壁面压力分布,分离区大小以及拐角干扰区密度等值线,计算结果与实验结果基本相符。研究还表明,对分离区大小和位置的预测还受到RANS/LES混合模式中的混合函数的影响。
研究高阶紧致差分方法构造原理,发展求解Boltzmann模型方程气体运动论高阶格式。使用BGK、Shakhov和ES椭球碰撞模型,以一维激波管、激波内流动和二维翼型、圆柱绕流为算例,开展不同流区一维、二维绕流高阶格式应用研究。将计算结果与实验数据、DSMC模拟值对比分析,验证高阶紧致格式用于求解Boltzmann模型方程可靠性,分析不同精度格式和碰撞模型对气体运动论数值算法效率与精度影响。
通过对速度分布函数的高阶展开和对初值的高阶重构,构造了三阶气体动理论BGK格式,并通过几个典型算例验证了新格式在光滑流场和间断流场的性能。研究表明,相比于传统的基于Riemann解的高阶格式,本文提出的格式不仅考虑了网格单元界面上物理量的高阶重构,而且在初始场的演化阶段耦合了流体的对流和粘性扩散,也能够保证解的高阶精度。
基于直接模拟Monte Carlo(DSMC)法发展了一种处理复杂几何边界时与网格无关的MonteCarlo方法,称之为GI-DSMC(Grid-Independence DSMC)法。该方法既保留了规则网格在处理直边界时简单的优点又引入了分子与复杂几何形状边界作用时的简易处理方案,从而消除在处理复杂几何边界的网格生成时易导致的网格与物理系统的不匹配性。以圆球绕流为例,给出气体分子在圆球表面反射时
针对复杂外形高超声速流场,将平衡粒子模拟方法(EPSM)与蒙特卡罗直接模拟方法(DSMC)结合,构造了可模拟含近连续流区及过渡流区的非结构网格DSMC/EPSM混合算法,采用了适用于非整数自由度振动能激发的EPSM方法,使算法能够模拟高温非平衡效应及化学反应,并通过算例验证了算法的有效性。
本文针对影响热环境计算精度的网格尺度问题,结合分子运动论的成果,提出基于物面局部流动参数的法向网格模拟准则,并通过不同数值方法的适用性和复杂外形的适应性等计算分析,验证了此网格模拟标准的可用性。
研究求解喷管内流场N-S方程数值方法,发展基于N-S方程物面边界滑移流计算技术。发展求解羽流核心区轴对称DSMC方法与远场三维DSMC仿真方案,建立多组元混合物羽流DSMC仿真方法。研究求解卫星姿控发动机内外近场、远场、倒流区和羽流与物面相互作用影响区多流域流场分区N-S/DSMC耦合计算技术,提出一套用于求解发动机燃气混合物羽流及对太阳电池帆板与卫星体表面撞击污染影响数值模拟方法。通过某在轨卫星
在二维非结构网格上,对结合斜率限制器的间断Galerkin有限元方法求解二维欧拉方程进行研究。该斜率限制器为van Albada类型,从有限体积法中演变而来,易于在非结构网格上实现。文中使用Lax-Friedrich和HLLC通量函数计算网格单元边界处的数值通量,时间推进采用经典四步Runge-Kutta方法。运用该方法对NACA0012在跨声速情况下的无粘流动进行数值模拟,结果表明：采用该类型斜