基于自适应笛卡尔网格的可压缩粘性流动数值模拟

基于自适应笛卡尔网格方法求解Navier-Stokes方程,网格以四叉树数据结构存储,固壁边界条件通过一种虚拟单元体方法引入。在几何外形确定的前提下自动完成网格生成、加密和流场的求解任务。对含有激波的NACA0012翼型的超音速绕流工况和含有回流区的双NACA0012翼型亚音速绕流工况进行了数值模拟,并与现有的非等距笛卡尔网格解和非结构网格解进行了对比。结果表明：基于自适应笛卡尔网格能够准确模拟可压缩粘性流动,同非等距笛卡尔网格相比,自适应技术的使用显著降低了网格量,但是同非结构网格相比,现有的自适应笛卡尔网格技术在边界层的分辨上效率较低,有待进一步发展。     

基于自适应笛卡尔网格的可压缩黏性流动数值模拟

基于自适应笛卡尔网格方法求解Navier-Stokes方程,网格以四叉树数据结构存储,固壁边界条件通过一种虚拟单元体方法引入。在几何外形确定的前提下自动完成网格生成、加密和流场的求解任务。对含有激波的NACA0012翼型的超音速绕流工况和含有回流区的双NACA0012翼型亚音速绕流工况进行了数值模拟,并与现有的非等距笛卡尔网格解和非结构网格解进行了对比。结果表明:基于自适应笛卡尔网格能够准确模拟可压缩黏性流动,同非等距笛卡尔网格相比,自适应技术的使用显著降低了网格量,但是同非结构网格相比,现有的自适应笛卡尔网格技术在边界层的分辨上效率较低,有待进一步发展。     

基于自动网格生成技术的精细油藏数值模拟器

针对笛卡尔网格油藏数值模拟软件不能精细模拟油藏特殊区域的问题，在综合一些自动网格生成算法优点的基础上，提出了一套新的适合于有限体积差分的网格自动生成算法，用面向对象语言C^ 开发出能对油藏特殊区域生成柔性网格的全隐式油藏数值模拟软件－VGRNS，并通过实例将此软件与笛卡尔网格划分的油藏数模软件进行了对比。     

基于自动网格生成技术的精细油藏数值模拟器

针对笛卡尔网格油藏数值模拟软件不能精细模拟油藏特殊区域的问题 ,在综合一些自动网格生成算法优点的基础上 ,提出了一套新的适合于有限体积差分的网格自动生成算法 ,用面向对象语言 C+ +开发出能对油藏特殊区域生成柔性网格的全隐式油藏数值模拟软件—VGRNS,并通过实例将此软件与笛卡尔网格划分的油藏数模软件进行了对比 .     

笛卡尔切割单元法在湖泊二维流场模拟中的应用

针对湖泊几何边界复杂的特点,为了高效并准确地计算流场,采用了基于有限体积法的笛卡尔切割单元法.该方法在流场主要区域采用笛卡尔网格,而对与固体边界相交处的单元采用特殊的处理.这样既保留了笛卡尔网格求解速度快的优点,同时在边界处又保证了精度,避免了传统差分法处理该类问题需采用阶梯形网格逼近而导致的失真问题.通过与一个标准算例的贴体网格计算结果对比,显示计算结果一致.最后运用笛卡尔切割单元法,求解了杭州西湖的二维流场,结果表明该方法适用于具有复杂几何边界的湖泊流场求解,且精度较高.     

带浸入边界法的新型五阶WENO格式求解双曲守恒律方程

采用一种带浸入边界法的新型五阶有限差分WENO(weighted essentially non-oscillatory)格式在笛卡尔网格上求解含有复杂物面的双曲型守恒律方程。这种结构网格上的新型WENO格式因对计算网格质量依赖性较高,故一般不能直接应用于上述问题的数值模拟。而浸入边界法是一种能较好处理复杂物面边界的方法。将两种方法结合起来,可在笛卡尔网格上数值解决跨音速复杂流动问题,并用四个经典算例验证新型五阶WENO方法的有效性。     

自适应网格在水平集方法中的优化与应用

在水平集方法对运动界面进行追踪的过程中,利用自适应笛卡尔网格可以在水平集演化的关键区域对网格进行自动细化.水平集方程每个时间步求解后,细化区域和粗化区域的选择是决定迭代求解效率的一个重要因素,由于不同细化级别边界处的节点部分缺少同级别的邻接点,加大了求解水平集演化方程中偏微分近似解的复杂度.提出一种优化后的笛卡尔自适应...     

自动网格体系在柱体绕流大涡模拟中的适用性评估

为了评估基于snappyHexMesh方法生成的自动网格体系在二维柱体绕流大涡模拟中的适用性,比较了该自动网格体系与人工网格体系对于Re为3 900圆柱绕流和Re为22 000方柱绕流的数值模拟结果。通过设置合理的计算域以及数值格式,采用snappyHexMesh自动网格以及人工网格的算例都表现出良好的数值稳定性。将不同网格体系的数值模拟结果与物理试验结果进行对比,结果表明,采用snappyHexMesh网格可以提高数值求解效率;圆柱绕流对网格体系的变化比较敏感,不同密度的snappyHexMesh网格会显著影响圆柱气动力特征以及尾流区域的流场结果;snappyHexMesh网格体系可以准确预测方柱绕流,在方柱绕流大涡模拟中具有相较于圆柱绕流更好的适用性。     

三维笛卡尔网格生成前处理软件测试

针对自主开发的三维笛卡尔网格生成前处理软件,基于CT扫描思想,提出了一种通过对三维网格进行二维切片/三维剖切显示来检测网格生成质量的测试方法. 通过该测试方法发现部分测试实例中实体内部存在异常空网格的问题,该问题是由于三维实体的三角面片化近似引起的,文中给出了消除异常空网格的实现方法. 网格生成测试算例验证了本文所给出的软件测试方法的有效性.     

水力压裂裂缝模拟研究

针对低渗油气藏水力压裂裂缝通常采用直角网格加密模拟准确性不高的问题,采用了油藏数值模拟中的新型网格结构——PEBI网格来模拟水力裂缝.在充分认识水力裂缝渗流机理的基础上,深入分析了数值模拟中处理水力裂缝的常规直角网格加密方法和PEBI网格方法的原理;并以某火山岩油藏为例运用两种方法对裂缝方位、导流能力和裂缝穿透比进行优化研究.研究结果表明,PEBI网格模拟水力裂缝的方法相对于常规直角网格加密方法具有较高的准确性和优越性.运用该方法可以为现场压裂施工设计和效果评价提供一种有效的研究手段,对于提高低渗油气藏压裂效果具有实际意义.     

二维非结构网格Euler方程高效解法

在非结构网格上应用多重网格技术加速 Euler 方程的收敛,在多重网格中通过聚合法进行粗网格生成,并对粗网格中的多边形网格做了等价面处理.在空间离散上采用 Roe 格式,在时间推进上分别采用了显式和隐式算法.通过对 NA-CA0012 翼型和 RAE2822 翼型的流场模拟,比较了显式多重网格法和隐式多重网格法的计算效率.     

二维非定常有涡流动的数值模拟方法

为有效地模拟二维有分离现象的粘性流动,发展了求解涡量流函数方程的数值方法。对涡量输运方程的时间项运用四阶Runge-Kutta法离散,从而将方程分解为4个计算步分别求解。在每一计算步对方程进行Steger-Warming近似因式分解,从而使涡量的计算在空间的两个方向上分别进行。对流项采用Chakravaythy-Oscher总变差减少(TVD)格式离散。流函数的Poisson方程运用Tschebyscheff SLOR方法交替方向迭代求解。运用该方法对突然起动圆柱和高雷诺数时弯曲薄翼的非定常有涡流动进行了数值模拟,并将计算结果与其它方法的计算结果和试验结果进行了对比,结果表明本方法具有精度高、数值稳定性好和计算效率高的优点。     

多段翼型绕流数值分析对比及结冰影响研究

飞机结冰是影响飞机飞行安全的重要因素之一,研究结冰对多段翼型的气动特性影响具有重要意义.应用多块结构化网格生成方法,采用两种不同的方法求解雷诺平均N-S方程分别对多段翼型绕流流场进行数值对比分析,然后选用一种较合适的方法进一步针对三种不同的冰型:钝头体、双角体和尖头体冰型,分析不同形状的冰型对多段翼型绕流流场及气动特性的影响.     

基于伴随方法的动态非定常气动外形优化设计

为实现动态非定常条件下的气动外形优化设计,基于非定常形式的欧拉方程及伴随方程,建立了一套完整的针对非定常动态问题的气动外形优化设计方法,并采用所建立的方法,对俯仰振荡的NACA0012翼型在跨声速(Ma=0.8)及限制条件下进行了翼型时均减阻外形优化. 对优化设计的结果进行了分析. 研究结果表明:该动态非定常气动外形优化设计方法具有较高的可靠性和实用性.      

基于二维非结构网格的全隐式,迎风格式欧拉方程研究

对二维非结构网格进行了网格重构,使用格点格式有限体积法离散流场方程,采用Roe矢通量分裂方法计算通量,时间推进采用Lu-SGS方法.通过对NACA0012翼型和RAE2822翼型的流场模拟,验证了本算法的可行性和高效性.     

低雷诺数下翼型粘性绕流主动流动控制的数值模拟

通过数值求解Navier-Stokes方程模拟零质量射流对低雷诺数翼型绕流的控制作用.空间离散采用中心有限体积格式,时间推进为双时间推进方法.计算结果表明,采用零质量射流进行主动流动控制能有效抑制大攻角下大尺度的流动分离,改善翼型的气动特性 ,起到显著的增升减阻作用.     

低雷诺数翼型分离流动的大涡模拟研究

采用高精度大涡模拟算法,对低雷诺数下的孤立翼型分离流动问题进行研究,计算了雷诺数为55000、马赫数0.2、来流5°攻角下的NACA-0025翼型,生成数值数据库,从时均流场、瞬态流场、频谱和高阶统计量等多个角度进行分析.研究结果表明:大涡模拟方法能够很好的描述低雷诺数翼型分离流动,其瞬态流场图画与实验结果吻合的很好;翼型上表面出现大尺度的开放式分离区,在Kelvin-Helmholtz(K-H)不稳定性作用下,自由剪切层失稳卷起展向涡,展向涡二次失稳发生旋涡配对现象;分离区流场的演化受大尺度涡结构控制,流场中高阶统计量的分布也与涡结构密切相关.     

扇翼飞行器绕翼型流动数值研究

扇翼飞行器以其低速大载荷、低噪音特性在军用及民用领域有着广泛的应用前景．基于CFD软件采用滑移网格技术对简化为二维的绕翼型流动进行数值模拟,得到低速大迎角下绕翼型的流场,分析并比较了不同迎角下的流场结构,给出了升阻力随迎角的变化关系．在此基础上通过比对分析翼型局部几何形状改变后得到的数值模拟结果,得到了飞行器机翼前缘高度对绕翼型流场的影响以及相应的升阻力变化情况,并讨论了引起这些影响和变化的可能原因．最终给出可用于扇翼飞行器设计的方案．     

欧拉两相流法在翼型积冰中的应用

应用欧拉两相流理论对各种结冰气象条件下翼型积冰进行了数值模拟,采用格心型有限体积方法求解N-S方程来获得翼型的流场.对4°攻角下NACA0012翼型前缘的霜冰,明冰,混合冰进行了预测,最后将计算的结果和文献中实验数据进行对比,结果吻合良好,说明该方法是有效的.同时对积冰翼型的气动特性进行了分析,并比较了不同冰形对升、阻力的影响,计算结果显示明冰对气动特性的破坏最大.