基于非结构自适应网格的二维Euler方程数值求解方法研究

提出了一种基于非结构自适应网格的二维Euler方程的数值解法.采用有限体积法进行空间离散,通量计算采用Jamson中心格式,使得它适用于任意多边形计算单元.为了得到定常解,采用一种显式的四步Runge-Kutta迭代方法对时间进行积分.根据流场参数的变化梯度确定加密边,由加密准则进行自适应网格剖分,然后得到分布合理的加密过后的网格.求解二维Euler方程,对NACA0012翼型进行了数值模拟,通过对自适应前后的数值解的对比,说明所建立的方法是正确的.     

二维喷动床CFD-DEM模拟中曳力模型的影响

CFD-DEM已经广泛应用到喷动床的研究中,其模拟的准确性与用于处理颗粒-流体相互作用的曳力模型密切相关。为了探究不同曳力模型对喷动床CFD-DEM模拟结果的影响,基于非结构化网格的喷动床仿真,使用7个曳力模型分别对锥底喷动床内气固两相运动进行了数值模拟。综合床层压降、喷动高度和颗粒速度特性三个方面,Wen-Yu模型和Gibilaro模型预测的气固两相运动最剧烈,其次是Di Felice模型、Syamlal-O’Brien模型、Gidaspow模型和Huilin-Gidaspow模型,BVK模型预测的气固两相运动最平缓。由于模拟的气固两相体系属于密相体系,Huilin-Gidaspow模型的光滑过渡函数没有产生效果,所以Gidaspow模型和Huilin-Gidaspow模型在各个方面的预测结果基本一致。     

基于FlowStar软件的栅格舵气动特性模拟

栅格舵是由外部边框和内部若干薄栅格组成的一种新型气动力面和控制面。由于其具有尺寸小、重量轻和易折叠等结构特点,以及升力特性好、铰链力矩小和压心位置随马赫数变化小等良好的气动性能,越来越受到重视和广泛应用。但是,由于栅格舵是特殊的蜂窝结构,常规的结构网格生成困难,并且栅格壁之间存在严重的波系干扰,流场结构复杂,给数值模拟带来了挑战。本文针对栅格舵流场结构复杂、网格生成和数值模拟难度大的难题,采用国家数值风洞工程(NNW)项目团队自主开发的非结构混合网格流场软件FlowStar,对栅格舵气动特性开展了数值模拟研究。首先,通过十字栅格舵全弹标模外形,验证了数值模拟方法的可靠性。然后,通过研究栅格不同形状对栅格舵气动特性的影响规律,给出了栅格舵设计时尽量采用气动性能较高和结构强度较好的菱形栅格,以及尽可能少用或不用正三角形栅格的指导建议。最后,开展了类Space X迎风面局部弧形栅格舵气动特性模拟研究,结果表明,弧形后掠可以有效降低阻力,提高升阻比;攻角为0°时,与平直栅格舵相比,弧形后掠栅格舵在亚跨声速阶段,阻力降低约8%,马赫数大于2时,阻力降低约15%。本文研究可为栅格舵的设计提供参考和...     

在动网格上的多介质界面数值处理方法研究

运用动网格上的ALE方法对一维可压缩多介质Riemann问题进行求解,在处理介质界面 上的数值通量时提出了3种不同的数值方法：Lagrange方法、GFM和HLLC方法,并且对 这3种方法的数值结果进行了比较,认为GFM方法和HLLC方法在介质界面上出现大压力 梯度时能够明显消除界面上密度的非物理震荡,提高介质界面上数值解的精度.     

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.