非结构网格在周期性充分发展流动数值模拟中的应用

为了解决周期性充分发展流动数值模拟的计算复杂性，该文提出了非结构网格上周期性边界进出口条件相互回代的方法。应用此方法对横掠叉排管束周期性充分发展层流流动进行了数值模拟，并将其结果与流动进口内8个管束层流流动的充分发展区域的数值模拟结果进行了比较，得到了较一致的结果。实验表明：单周期上进出口条件相互回代的方法可以替代多周期的数值计算，从而简化了计算过程，节省了计算时间。     

基于单流道和结构网格的湿式潜水泵数值模拟

采用流场模拟软件ANSYS-CFX对BQ1100-425型湿式潜水泵进行数值模拟,在完成网格无关性验证的基础上,分别研究了全流道与单流道模型,结构网格与非结构网格模型的数值模拟结果,如压力、速度等内流场和扬程、效率等外特性。结果表明:全流道数值模拟相比于单流道能更准确地反映真实流场的流动情况,但计算时间较长,计算效率低。单流道数值模拟在一定范围内能较好地反映真实流场的运行情况,且计算时间较短。结构网格与非结构网格都能准确反映实际流场的运行变化情况,但结构网格计算准确性更高,计算时间短,效率高。鉴于湿式潜水泵结构的复杂性,基于单流道结构网格的数值模拟方法能有效地模拟内流场并进行扬程、效率等性能预测。     

副油箱从机翼分离流场的数值模拟

采用非结构网格的方法对含有动边界的流场进行数值模拟，非结构动网格的生成采用了改进的弹簧方法和局部网格重新生成相结合的方法，在网格局部重新生成时，采用了阵面推进方法.以基于格心的有限体积法为基础，空间二阶精度，网格单元边界处守恒量通量的计算采用Hanel方法.分别采用四阶Runge-Kutta方法、GMRES方法求解基于ALE形式的Euler方程对振动的NACA0012翼型绕流进行了计算，两种方法都是可行的，其中GMRES隐式方法计算的结果与文献中结果以及实验结果吻合较好.对副油箱从机翼分离的过程进行模拟，给出的副油箱运动趋势和风洞实验结果相吻合.     

不同湍流模型数值模拟三维轴对称凸体分离流动的比较

采用基于非结构化网格的有限体积法对三维轴对称凸体分离流动进行了数值模拟.分别采用标准k-ε模型、可实现(Realizable)k-ε模型、低Re数k-ε模型、剪切应力输运(Shear Stress Transport,SST)k-ω模型和低Re数SST k-ω模型,进行了数值模拟和计算.给出了计算模型、流场网格以及计算结果,并与实验结果进行了对比,结果发现:低Re数k-ε模型模拟效果最佳,与实验结果最为接近.     

通道中有稳定器的湍流流场的数值计算

采用非正交曲线坐标系下的ＳＩＭＰＬＥ方法对含有稳定器的通道湍流流场进行了数值计算,网格采用分区方法生成,计算时对整个流场进行迭代计算而不需要分区迭代,避免了分区计算中界面上的插值计算．文中还对后台阶流动进行了计算,结果和实验比较吻合．     

斜流式喷水推进器内部三维紊流流动的数值模拟与性能分析

应用数值模拟方法研究斜流式喷水推进器内部紊流流场,基于紊流流场计算分析了其特性.采用三维贴体坐标网格对计算区域进行剖分离散,选择k-ε双方程紊流模型,应用SIMPLEC数值方法计算内流场各流动参数,进而根据流场的计算结果预测了该推进器内部泵在不同工况下的性能.计算结果与实验结果进行了对比,二者有较好的一致性.     

导弹构形横向喷流干扰流场数值模拟

利用对称TVD格式和重叠网格技术数值模拟了翼面位置不同的两种导弹构形的横向喷流干扰流场,研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性及翼面位置对干扰效应及气动性能的影响,并将计算得到的干扰力放大因子和干扰力矩系数的结果与实验结果进行比较,吻合良好.     

飞行器操纵面嗡鸣的非结构网格并行计算方法

为了数值模拟飞行器操纵面的嗡鸣现象,在集群计算机MPI并行计算环境下建立基于三维非定常欧拉方程耦合结构运动方程的嗡鸣计算方法.气动流场求解采用基于非结构网格的中心有限体积法进行空间离散,时间推进采用双时间方法,结构运动方程采用Adams预估校正方法求解.针对翼面与操纵面缝隙间存在的网格运动问题,在非结构网格系统上采用Delaunay图映射方法实现网格的运动变形.最后,使用飞行器操纵面标准嗡鸣计算模型对计算方法进行验证,结果表明:所建立的并行计算方法正确,程序具有很好的计算效率,能够对飞行器操纵面嗡鸣进行高效的数值分析.     

多面体网格在船体绕流计算中的应用

为了验证多面体网格在船体绕流、特别是存在自由面的高速船体绕流计算中的有效性,将多面体网格应用于动网格中,基于多面体网格对DARPA Suboff潜艇进行了数值计算;基于混合网格和动区域法对一种高速单体穿浪船进行了数值计算,并将数值计算结果同水池实验结果进行了对比.研究表明:同四面体网格相比,多面体网格能有效减少非结构网格的单元数量,并能用更少的网格单元获取相似甚至更好的计算结果;多面体网格在潜艇绕流计算上的应用效果非常好;结合Robust(Octree)网格生成方法可改善多面体网格的正交性和分层性,进而能够更好地捕捉波形和飞溅;混合网格结合动区域法可以将多面体网格应用于包含船体运动的数值计算.     

新的多重网格算法在压气机流场计算中的应用

根据在流场数值计算中使用多重网格方法可以加速收敛的原理,在叶栅流场控制方程使用有限体积差分法进行离散的前提下,提出了一种新的4重网格方法来加速求解压气机叶栅可压缩流动,该方法可以与显式时间推进法有效配合使用。通过对NASA67超音速轴流压气机叶栅流场的数值计算表明,该方法明显加快了流场计算的收敛速度,并且计算结果与试验结果吻合较好。     

裂缝性油藏非结构四边形网格数值模拟研究

在裂缝性油藏数值模拟中,引入非结构四边形网格系统是目前高精度数值模拟研究的热点,它能对复杂油藏边界、交叉裂缝、多条裂缝、多方向分布裂缝、复杂分支结构井等条件下任意连通区域进行网格剖分,相对其他非结构网格,它需要的网格数量少,且精度更高。建立连续介质模型和离散裂缝模型相结合的混合渗流模型并数值离散求解发现,非结构四边形网格的引入导致数值模拟计算效率严重下降。针对此问题,提出了3种解决方案：改进Paving算法的裂缝末端形态变形法、二维和三维网格编号重新排序方法和MPI并行计算方法,计算效率平均提高70%。     

湍流射流火焰的PDF模拟及误差分析

将有限容积法和MonteCarlo法结合起来求解概率密度函数 (PDF)方程的混合算法 ,是湍流燃烧数值模拟的PDF方法中目前最先进的算法 .论文将此算法发展到无结构网格中 ,并以湍流射流火焰为算例 ,将数值模拟的结果和实验数据作了比较 .文中分析了该混合算法数值误差的来源 ,指出影响误差大小的参数 ,并通过数值实验研究了统计偏差受样本数目影响的规律     

非结构网格RKDG方法求解多介质界面问题

在多介质流体力学方程数值模拟中,界面的计算是一个非常重要的问题.采用非结构网格上RKDG(Runge-Kutta Discontinuous Galerkin)有限元方法,给出了捕捉多介质界面的Lev-el Set方程和重新初始化方程的计算方法.数值实验表明多介质界面在非结构网格计算上,RKDG有限元方法具有很大的优点.     

无结构三角形网格下运动界面追踪的WENO有限体积法

在处理运动界面追踪问题的流体体积函数(VOF)法的基础上,给出了一种无结构三角形网格下的高分辨率的运动界面捕捉方法.该方法采用高精度的加权本质无振荡(WENO)有限体积格式离散VOF函数的空间导数,采用三阶TVD Runge-kutta方法离散时间导数,采用Lax-Friedrichs通量作为数值流通量.数值试验结果表明,用该方法来进行旋转速度场和剪切速度场的运动界面追踪,可以得到与理论解非常一致的追踪结果.     

激波诱导圆柱运动数值模拟

该文采用基于含有运动边界自适应非结构网格点的显式有限体积方法求解Euler方程，提出了一种运动边界的数值处理方法，对激波诱导的圆柱运动进行了数值模拟。由于物体的运动，物体边界附近的网格将发生变形。对运动边界附近严重变形的网格采用局部网格重新生成的方法进行处理。采用HLLC方法来处理运动边界问题，计算并比较了气流流过静止圆柱和圆柱在静止流场中运动2种情形。结果表明，应用HLLC方法处理运动边界是可行的。     

采用非结构化网格求解颗粒流的一种数值方法

离散颗粒模型是计算稠密颗粒两相流的一种重要模型,该模型通常采用硬球模式描述颗粒间的作用。提出一种在二维非结构化网格下结合硬球模式求解颗粒群运动的数值方法,对颗粒间的碰撞检测策略进行了优化,通过对颗粒采用局部移动的方法大大提高了计算的效率。数值算例显示在颗粒粒径小于网格尺度时,该方法能适应多种类型的网格,这为进一步实现不规则区域下稠密颗粒两相流的数值模拟奠定了基础。     

非结构网格下两步压力校正算法的潮流模拟

针对非结构网格下的潮流数值计算,采用有限体积法离散基本方程,对流-扩散项离散采用幂率格式,交叉扩散项引入超松弛校正方法,水位校正方程应用Rhie-Chow动量插值思想和SIMPLEC类算法导出.采用两步压力校正法将交叉扩散项单独考虑,以提高模型对高度畸变网格的适应能力.计算结果表明,该处理方法可提高压力欠松弛系数,并能有效改善复杂区域潮流模拟的健壮性.通过对黄浦江感潮河段潮流的模拟,表明该模型的预测结果与实测资料吻合良好,模型对潮流的模拟结果是可靠的.     

HLLC方法在二维非结构动网格上的应用

以自适应非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的Godunov一阶精度方法求解Euler方程，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，空间和时间都是一阶精度，对含有运动边界的多体干扰流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步进行多体分离数值模拟打下了基础。     

非结构网格上动能守恒的低耗散数值方法

非结构网格下的大涡模拟方法是解决叶轮机械内部复杂几何流道中强剪切湍流问题的有效途径,而降低算法的数值耗散是提高模拟精度和可靠性的关键之一。该文基于动能守恒特性推导了非结构网格下的低耗散数值方法,使其在维持格式稳定的前提下达到数值耗散最小,并在此基础上发展了与其相适应的无反射边界条件。进一步选取了5个具有代表性的基础算例进行模拟,结果表明:该文方法较传统迎风格式具有低耗散的特性且计算量减少了近75%。     

非结构四面体网格上三维非线性扩散方程的有限体积法

为了有效地数值模拟科学和工程中有广泛应用的非线性扩散方程,在三维线性扩散方程非结构四面体网格的有限体积法的基础上,提出了一个计算非结构四面体网格上非线性扩散方程的有限体积法。方法采用网格单元中心作为计算节点,相对于网格点的方法,计算量减少了一半。用L agrange因子法得到网格点上的值,考虑了网格中心点和网格点的相对位置,更适应大变形的网格。利用算子分裂,使计算更加简单。用N ew ton-B iCG STAB法来求解得到非线性方程组。数值结果表明:该方法具有二阶精度、保持通量守恒、对大变形的网格适应性强。