K-∈湍流模型下门桥水阻力的数值计算

采用RANS方程和RNGκ-ε二方程湍流模型及VOF算法对绕门桥的自由表面流动进行了数值计算。数值计算时，使用有限体积法离散控制方程，对流项采用二阶迎风差分格式，扩散项采用中心差分格式。对门桥模型进行了一系列速度下的数值模拟；计算结果与试验测量数据相比较，该计算方法具有较好的计算精度，可用于门桥性能分析和优化。     

带端板地效应翼性能的数值研究

通过求解不可压缩流体ＲＡＮＳ方程,数值模拟带端板三维地效应翼的性能及周围流场。数值方法引进了Ｃｈｏｒｉｎ的人工可压缩性概念,应用近似因式分解技术同时求解速度和压力场,动量方程对流项用二阶迎风差分格式离散,其余空间导数项均采用四阶精度的中心差分格式离散,时间离散采用欧拉隐式格式,计算在非交错网格上进行,为了避免压力场的振荡,在连续方程中隐式地加入了压力的四阶数值耗散项,湍流计算采用了Ｂａｌｄｗｉｎ－     

网格Peclet数和网格尺寸对对流扩散方程差分格式精度的影响

用泰勒级数展开法对对流扩散方程中的扩散项进行了理论分析,从而证明对于给定的差分格式,不仅网格Ｐｅｃｌｅｔ数不能任意选取,而且网格尺寸也不能随意设定。用四种格式对一维有源对流扩散方程进行了计算,结果表明网格尺寸对差分格式精度的影响比网格Ｐｅｃｌｅｔ数更明显。为了得到真实可靠的结果,所用差分格式的阶数愈高,相应的网格尺寸就必须愈小。二阶以上的高精度迎风差分格式与二阶迎风格式相比,并无明显的优势,建议一般不必要采用二阶以上的高精度差分格式。     

网格Peclet数和网格尺寸对流—扩散方程差分格式精度的影响

用泰勒级数展开法对流－扩散方程中的扩散项进行了理论分析,从而证明对于给定的差分格式,不仅网格Ｐｅｃｌｅｔ数不能任意选取,而且网格尺寸也不能随着设定,用四种格式对一维有源对流－扩散方程进行了计算,结果表明网格尺寸对差分格式精度的影响比网格Ｐｅｃｌｅｔ数更明显,为了得到真实可靠的结果,所用差分格式的阶数愈高,相应的网格尺寸就必须愈小,二阶以上的高精度迎风差分格式与二阶迎风格式相比,并无明显的优势,建议     

PWCD 格式及其在湍流流场计算中的应用

针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,提出了一种含有压力梯度项的对流扩散方程的压力加权中心差分格式（ＰＷＣＤ）,且给出了通用的系数表达式．通过对叶栅流道中二维湍流流动的计算结果表明：ＰＷＣＤ格式与ＥＤ格式相比,具有二阶以上精度和很快的收敛速度并且能够得到更合理的压力场和速度场．     

阻尼梁振动方程的高阶数值格式

本文研究了带有阻尼项的四阶梁振动方程初边值问题,基于紧致差分方法,给出了数值求解该问题的四种高阶紧致差分格式.对方程中的一阶和二阶时间导数项采用中心差分离散,对四阶空间导数项分别采用五点、七点和带紧致的五点、七点四种方法进行离散,得到四种高阶紧致差分格式,这四种格式均在时间方向达到二阶精度,在空间方向分别达到二阶、四阶...     

一种求解时间关联型缓坡方程的数值方法

以含底摩阻能量耗散项的时间关联型缓坡方程为控制方程,建立了一种具有二阶精度的数值离散格式.该格式对时间导数使用Euler预测-校正格式离散;对空间导数使用中心差分格式离散.基于统一边界条件表达式,对边界条件进行处理.数值解与物理模型实验值吻合较好,表明数值模拟模型可以有效地模拟波面随时间和空间的变化以及波高的空间分布.     

使用PAFV湍流模型对叶栅流动的数值研究

为进一步发展更有效的湍流模型,提高对高负荷跨声速压气机流动的数值计算精度,以侧偏平均脉动速度和应变率张量为基础构造了新的侧偏平均脉动速度(partial average fluctuation velocity, PAFV)湍流模型。在使用Python计算机语言进行编程计算过程中,将流动控制方程在一般曲线坐标系下变换后,采用有限差分方法进行空间离散。方程对流项采用Steger-Warming通量向量分裂法,扩散项采用中心差分格式离散,时间导数项采用具有TVD(total variation diminishing)性质的二阶Runge-Kutta法离散。使用该湍流模型对德国航空航天研究试验院的L030-4叶栅亚声速和超声速流动进行了计算,在进行网格无关性研究基础上,获得了流场压力等值线云图、马赫数等值线云图,压力系数分布等结果。计算显示侧偏平均脉动速度在激波附近区域受到明显抑制,避免了激波附近区域出现过大的湍流黏性。数值计算结果与实验符合比较好,初步验证了PAFV湍流模型可以有效地模拟该类型的叶栅流动。     

一维四阶双曲方程的紧致差分格式

本文主要研究一维四阶双曲方程初边值问题.首先通过引入一个中间函数将其转化为二阶方程组,然后对方程中的空间导数项采用四阶紧致差分格式离散,时间导数项采用二阶中心差分格式离散,构造出问题的隐式紧致差分格式.数值算例表明该格式具有较好的计算效果.     

Kdv-Burgers方程的两层线性化隐式差分格式

应用非线性对流项和反应项的两层线性化技巧,对非线性Kdv-Burgers方程周期边界问题构建了一类具有二阶截断误差的两层线性化隐式差分格式.用数学归纳原理和离散能量法建立了差分格式的唯一可解性、在最大模意义下的收敛性和稳定性.数值计算表明,该格式在时间和空间上都是二阶收敛的.     

离心泵叶轮内部清水湍流的动态大涡模拟

提出旋转坐标下的单相液体湍流的二阶双系数动态亚格子应力大涡模拟模型并加以验证。考虑亚格子应力的对称性和量纲一致性 ,在基于应变率张量和旋转率张量的不可约量及 Smagorinsky模型和亥姆霍兹定理的基础上 ,提出亚格子应力应表示为可解的应变率张量和旋转率张量的函数。在贴体坐标系下 ,利用交错网格系统和有限体积法离散湍流控制方程 ,采用 SIMPL EC方法求解方程。水泵叶轮中湍流的流速分布规律和压力分布等计算结果都与试验结果基本吻合 ,从而证明了此模型的可行性 ,以及计算方法和程序的可靠性     

基于任意多面体网格的Navier-Stokes方程并行求解器

发展了一种基于任意多面体网格的Navier-Stokes(NS)方程并行求解器,基于积分守恒形式NS方程组,发展了支持任意多面体网格的中心有限体积方法。采用基于面的连接关系,对不同拓扑类型的网格[如结构网格、非结构混合网格、多面体网格及剪切六面体(TRIMMER)网格]进行统一处理。对于无黏通量空间离散格式采用HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)格式,湍流模型采用k-ω两方程模型,时间推进采用适合并行计算的DP-LUR(data-parallel lower-upper relaxation)格式的隐式算法。对RAE2822翼型和ONERA M6机翼的结构网格、混合网格、多面体网格及TRIMMER网格进行了对比验证,结果表明发展的求解器具有较好的网格普适性,在不同拓扑类型的网格上均能求得较为接近实验值的结果;对多面体网格测试了加速比和并行效率,并行计算大大提高了计算效率。证明了求解器具有宽广的网格适应性,能够较为稳定、快速、准确地模拟定常绕流问题。     

双级旋风分离器特性的计算机预报

利用计算流体力学ＣＦＸ４．２ 软件对双级旋风分离器的气固两相流进行了三维数值模拟,采用贴体坐标系统描述复杂的几何边界和基于交错网格的网格系统;对于气相流场,采用雷诺应力模型计算旋转流动,在拉格朗日坐标中计算粒子运动．为优化求解,采用多块结构适体网格,利用ＳＩＭＰＬＥ 算法解离散方程;此外,给出了双级旋风分离器内不同截面的速度分布和５ 种不同尺寸共６０ 个粒子轨迹．研究结果表明,双级旋风分离器比常规的旋风分离器分离效果要好．     

离心风机紊流次声理论与控制的研究

分析了离心风机紊流次声形成的机制，探讨了利用叶轮网罩整流来抑制紊流次声的原理和方法，推导了离心风机紊流诱发的次声模型．通过对叶轮装网罩后叶轮叶道内气动特性的分析，给出了装网罩后叶轮紊流次声辐射功率、叶片尾缘涡流脱落次声辐射功率及网罩次声降噪量的数学模型，并以某窗式空调器离心风机为例进行了次声测试与理论计算，结果表明，理论与实际是相符的．     

叶栅梢部与固壁间隙流动的数值模拟

通过求解不可压缩流体RANS方程，对叶栅梢部与壁面之间的粘性流动进行数值模拟，计算模型基于Chorin的人工可压缩性方法，采用Baldwin-Lomax零方程湍流模式封闭控制方程组，针对间隙流动的特殊性，采用流场分块的方法以提高网络生成质量和计算稳定性，数值计算结构和试验数据的比较表明，该方法能够准确地模拟间隙流动。     

Zonal Embedded Grids for LES of Self-Sustaining Structures of Wall Turbulence

<正>Introduction With the development of computation fluid dynamic (CFD) and the need for complicated flow simulation, more and more research has been conducted using com-putational grids. For the multi-block structured grids of large eddy simulation (LES) or direct numerical simu-lation (DNS) frames, the local grid refinement method     

三维非定常不可压N-S方程一个高精度求解格式

对三维、非定常、不可压Navier-Stokes方程提出了一种新的数值计算方法.在时间离散上应用3阶精度混合显-隐相结合的分裂格式,空间离散在x及y轴向采用非等间距网格的紧致有限差分格式与z向应用Fourier谱展开相杂交的数值方法逼近.经平板边界层流的验证表明,该算法具有计算精度高、稳定性好、收敛速度快等特点.同时也研究了三维、非定常流体运动下游边界问题,提出了无反射出流边界条件,以减少在有限计算区域内人工出流边界反射引起的数值误差,保证直接数值模拟的精度和准确性.该算法的提出对于求解边界层、射流及混合层等流动中的转捩与演化问题具有重要的理论意义.     

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

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

离心压气机实际气体三维粘性流场分析

应用时间推进法求解连续性方程、时均N-S方程、能量方程、实际气体状态方程及Yang＆Shih紊流模型,对压缩重介质实际气体(丙稀)的大流量离心压气机整级流场(闭式叶轮、扩压器、弯道及回流器)进行了数值分析,其中实际气体的热力学属性均表示成温度的函数.分析了压缩实际气体丙稀时模型级内部尤其是弯道及回流器内部的流动规律,将计算结果与丙稀作为理想气体时模拟的结果进行对比分析,并与实验值进行了比较.结果表明:采用实际气体模型进行模拟时压气机性能更接近于实际性能.压缩实际气体及理想气体时叶轮内部速度分布及模型级出口参数是不同的,叶轮出口流场参数分布差别主要集中在轮盘、轮盖两侧.沿弯道及回流器中央流线的压力系数差别显著,模型级出口理想气体的压力系数明显高于实际气体.为了加快收敛,计算时采用了隐式残差平均法及完全多重网格技术.     

Mixed Model for Silt-Laden Solid-Liquid Two-Phase Flows

The kinetic theory of molecular gases was used to derive the governing equations for dense solid-liquid two-phase flows from a microscopic flow characteristics viewpoint by multiplying the Boltzmann equation for each phase by property parameters and integrating over the velocity space. The particle collision term was derived from microscopic terms by comparison with dilute two-phase flow but with consideration of the collisions between particles for dense two-phase flow conditions and by assuming that the particle-phase velocity distribution obeys the Maxwell equations. Appropriate terms from the dilute two-phase governing equations were combined with the dense particle collision term to develop the governing equations for dense solid-liquid turbulent flows. The SIMPLEC algorithm and a staggered grid system were used to solve the discretized two-phase governing equations with a Reynolds averaged turbulence model. Dense solid-liquid turbulent two-phase flows were simulated for flow in a duct. The simulation results agree well with experimental data.