基于非错列网格的三维椭圆型流动问题的算法

采用虚拟错列网格和动量插值技术，在非错列网格上发展了求解三维随圆型（有回流）流动问题的差分算法和通用程度，用于典型层流问题的计算结果与实验数据和近似分析解吻合；用于三维湍流旋流问题也得到了合理的结果。     

NS方程的非结构化网格方法及其差分格式

采用同位非结构化网格上的有限体积法,对ＮａｖｉｅｒＳｔｏｋｅｓ 方程的ＳＩＭＰＬＥ 算法及差分格式进行了研究．提出了适用于非结构化网格的不用求解单元顶点变量值的二阶混合差分格式．该格式的优点在于：１） 减少了计算工作量;２） 避免了普通混合差分格式因使用简单的一阶迎风差分格式所引起的网格界面方向相关性的问题．最后,采用三角形网格,利用提出的方法及差分格式,对方腔内的驱动层流及绕翼型湍流进行了数值计算,计算结果与基准解或实验值的符合程度优于普通混合差分格式     

基于同位网格的一种新的定压边界处理方法

将SIMPLE-P算法中处理压力边界的思想推广到同位网格中,综合考虑压力和速度边界条件的数学特性,重新推导边界控制体的压力方程,借鉴IDEAL算法处理压力和速度耦合的思想,讨论定压边界流动传热问题的求解步骤。通过Kovasznay流、压力驱圆管层流和浮力驱动流等几个典型算例对新方法和SIMPLE-P算法进行考核并将计算结果与文献中的数值解和实验值进行对比。结果表明:新方法是准确有效的,对于定压边界问题的求解较之SIM-PLE-P算法具有更好的收敛性能;对于大空间自然对流问题使用新方法在较小的计算区域就可获得高精度的结果,从而节省了大量的存储空间和运算时间。     

圆柱绕流水动力学特性的数值模拟

应用计算流体动力学软件FLUENT对低雷诺数(Re)(Re≤300)下圆柱绕流问题进行数值模拟。采用分区结构化网格及层流模型求解不可压缩Navier-Stokes方程,在未施加任何扰动的条件下获得了圆柱绕流2种特殊的绕流场,系统分析了壁面压力系数、平均阻力系数、脉动阻力、升力系数及旋涡脱落频率等水动力学特性参数随Re数变化的规律,结果与实验研究吻合较好。     

通风分裂算法在有限元非结构网格中的推广

将在结构网格中发展起来的通风分裂算法推广到有限体积非结构网格中,并成功地进行了定常跨声速二维欧拉方程的数值计算．此格式具有流通量矢量分裂格式的简单性及通量差分分裂格式的高阶精度．最后给出的ＮＡＣＡ００１２翼型绕流的结果与实验结果符合的较好．     

火灾问题中湍流模型和数值方法的研究

本文用求解化学流体力学基本方程组的方法(场模化)研究火灾中常见的浮升流在室内的运动规律,着重比较了三种湍流模型和两种差分格式对计算结果的影响,提出了能改善收敛性的非均匀网格的指数型分布。运用本文推荐的湍流模型和差分格式,计算得到了室内气流速度和温度在空间的分布及其随时间的变化,与实验观测一致。     

基于AMG法及多块非正交网格的三维方腔流场数值模拟

为了将AMG法与多块非正交结构网格结合起来,应用于求解基于多块格离散的三维N-S方程离散系统,针对经典算例三维顶盖驱动直方腔和45°倾斜方腔流动问题,构建了三维多块非正交结构网格系统,设计了三维N-S方程块结构网格离散系统系数矩阵CSR存贮格式;根据非正交同位网格SIMPLE耦合算法,建立了N-S方程的AMG数值求解方法,编制了用AMG法计算的三维数值仿真程序,并采用AMG法求解控制方程中所有流动参数;同已有的结果做了比较,给出了三维方腔流场内复杂的涡系图像,并对算法的有效性进行了验证.结果表明:AMG法及相关算法可用于复杂区域流动问题的数值模拟,为大规模的CFD数值模拟问题奠定了基础.     

旋流泵流动研究

文章以分块非结构六面体网格划分计算模型,运用雷诺平均N-S方程和标准的k-ε双方程湍流模型结合SIMPLEC算法,数值模拟了旋流泵内部三维不可压湍流场,得出旋流泵内速度和静压分布,以此来更好的了解旋流泵的内部流动。     

基于FLUENT的离心泵内部流动的数值模拟

本文基于三维雷诺时均的Navier-Stokes方程和标准的κ-ε湍流模型,采用三维无结构网格及压强连接的隐式修正SIMPLEC算法,利用Fluent中提供的多重参考系(MRF)模型,对IB型离心泵进行内部流动的数值模拟,根据计算结果分析了离心泵内的流动规律,并揭示出了离心泵内存在漩涡、二次流等不良的流动现象。本文还将离心泵性能的预测值与实验值作了比较以验证计算结果的正确性。     

基于VOF模型的导弹低速入水数值模拟方法

基于商用软件Fluent框架,利用动网格技术和UDF方法实现了回转体低速0°攻角直入水过程的数值模拟。采用水、气2相VOF方法捕捉界面,分别采用湍流与层流建立了导弹入水冲击的数值模拟方法,并与实验获得的照片进行了对比分析。分析结果表明:层流模型下导弹入水冲击的数值模拟和试验结果吻合较好,湍流模型下的数值模拟结果误差较大,因此,在低速条件下,湍流模型不适用,应该用层流模型。     

求解椭圆问题的一类外推三层网格法

使用线性拉格朗日有限元离散一类二维椭圆问题,选择合适剖分尺度形成最粗网格、次粗网格和最细网格和对应的方程组。在最粗网格和次粗网格上使用外推法（新外推法或经典外推法）得到次粗网格上高精度近似解,然后使用三次样条插值为细网格提供初始值,结合磨光算子,构造了经典外推三层网格法和新外推三层网格法,并给出相应的数值实验。与通常的瀑布型多重网格法相比,数值实验表明了两种新算法计算精度更高,细层上迭代步数非常少,计算时间更短,具有较强的稳健性。     

二维机翼振动非定常流场的数值模拟

准确数值模拟机翼动态失速时非定常流场主要有湍流模型和数值算法两方面的难点。在研究湍流模型的基础上 ,采用了 q-ω低 Reynolds数双方程湍流模型 ,并采用了一种高精度、高分辨率的数值计算方法 (L U- SGS- GE隐式算法和四阶 MUSCL TVD格式 )数值模拟了振动机翼周围的附着流动、轻度失速和深度失速非定常流场。结果表明 :采用本方法计算附着流动和轻微失速状态流场时 ,其升力系数迟滞曲线和阻力系数迟滞曲线与实验结果吻合得令人满意。深度失速的计算精度比前人工作有明显改善     

单锥式油水分离旋流器内流场的数值模拟

为了研究水力旋流器用于油水分离的复杂情况,使用FLUENT软件中的多相流欧拉分析方法,结合雷诺应力湍流模型对单锥式旋流器的内部流场进行了数值模拟.分析了旋流器内部的体积浓度分布、压力分布,以及切向、轴向和径向速度分布的规律,揭示了油水两相流的分离特性.在不同流量下,计算出了旋流器的流量-效率曲线,计算结果与实验数据吻合较好,从而证明了该湍流模型和数值算法的可靠性.     

失速翼型气动性能及噪声特性计算

基于计算流体力学、数值研究风力机叶片翼型失速条件下的流动转变。目前大多数CFD商业软件采用的湍流模型假设翼型表面边界层完全湍流,而实际上存在流动状态由层流到湍流过渡;因此导致了数值研究结果与实验偏差较大。为预测不可压缩流动下翼型表面层流到湍流的过渡,采用K-w SST模型分析NACA0063翼型失速流动过渡行为,并与S-A湍流模型计算结果和风洞实验数据比对。为评价翼型气动噪声水平,同时对翼型自身噪声进行计算和研究。     

建筑结构非定常绕流风场的数值模拟方法

给出了一种适用于计算建筑结构非定常绕流风场的大涡数值模拟算法.该算法基于有限差分法,采用曲线坐标结构网格,能精确描述形状复杂的物面边界,为下一步准确模拟含有因结构受风变形所致运动边界的绕流场奠定了基础.该算法采用投影法解耦纳维-斯托克斯方程中的压力和速度,对非定常流场的时间步进采用二阶Adams-Bashforth方法.采用同位网格以减少计算所需内存,为了平抑同位网格下中心差分格式导致的固有压力波动现象,计算对流速度时采用Rhie-Chow动量插值方法,利用编制的曲线坐标系下大涡模拟数值计算程序,对德州理工大学(TTU)建筑足尺模型绕流风场进行了模拟,所得结果与现场实测和风洞试验结果进行了比较.结果表明,本文算法是建筑结构非定常绕流风场数值模拟的有效方法.     

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

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

基于计算流体动力学法的风场模拟和流/固耦合问题

介绍了计算流体动力学（CFD）方法中的各种湍流模型和数值求解方法,阐述了流/固耦合问题的描述方法与求解策略,分析了强耦合算法和弱耦合算法、流场和流/固耦合研究中的网格生成与更新技术,以及耦合界面的信息传递方法,指出了其网格建模中存在的一些难题.     

基于网格数值耗散修正的圆湍射流模型

基于民用机舱内个性通风口,以圆湍射流作为研究对象,以标准κ-ε湍流模型较精细数值模拟与实验对比结果为出发点,通过对粗网格下数值结果的分析,针对网格尺度增大导致计算中数值耗散相应增大、射流速度沿轴向和径向衰减过快问题,从湍流耗散率方程出发,对影响湍流耗散的生成项和耗散项进行修正,平衡因网格尺度产生流动衰减过快的因素,进而提出了耗散率修正的κ-ε湍流模型.采用该修正模型验证了不同雷诺数和非等温条件下射流的适用性.应用修正的湍流模型可以避免机舱内大尺度空间流动对网格的苛刻要求,明显减少网格总数,同时达到对舱内复杂结构和边界条件下流场特征的正确模拟.     

平面管道内钝体后湍流预混火焰的数值分析

本文基于运用SCASM 燃烧模型的“k—ε—gr”湍流模型,借助于GM80通用程序对平面管道内钝体后的湍流预混火焰进行了数值分析.计算得到的火焰扩展角大约为2°,和与来流状态无关的结果与普遍承认的实验结果一致。本文给出并求解了控制反应度脉动均方值gr 的微分方程,从而构造了“k—ε—gr”湍流模型;提出了以最大温度脉动值的轨迹作为火焰面的数学表征的建议;考虑了SCASM 燃烧模型中局部温度对层流火焰传播速度的影响;给出了计算近壁区gr 值的壁面函数.     

旋风分离器数值计算及性能实验分析

借助计算流体力学软件Fluent,采用三维贴体坐标网格,基于非稳态雷诺应力湍流模型,对旋风分离器内部流场进行数值计算.研究不同粒径固相颗粒的运动轨迹,揭示颗粒在分离器中的运动机理,得到旋风分离器内部气流切向速度、轴向速度及切面旋转矢量速度的分布规律,并与实验测试值进行比对.结果表明:在相同条件下,数值计算与实验测试结果非常接近,能很好地预测切向速度的"驼峰"结构及轴向速度分布的上行流和下行流;随着颗粒粒径的增加,分离器外壁呈螺旋流分布,内部流夹带随粒径的增加而逐渐减小.