基于RNG k-ε湍流模型钝体绕流的数值模拟

采用RNG k-ε湍流模型对钝体绕流流场进行了数值模拟．计算采用非交错网格的有限体积法求解二维不可压N-S方程．计算结果表明，RNG k-ε湍流模型可以成功地模拟绕钝体的不稳定、非定常和剧烈分离流动。     

微射流改变绕流圆柱气动性能的数值模拟

采用有限体积法求解Faver平均N-S方程,对微射流改变绕流圆柱的气动性能进行数值模拟。计算结果显示了在微射流作用下绕流圆柱的气动性能发生改变,揭示了微射流产生的低压回流区是绕流圆柱气动性能发生改变的根本原因,该结果和文献[1]中实验观察到的现象比较吻合。     

民用飞机翼身整流粘性绕流数值模拟方法研究

利用CFD方法,针对民用飞机开展翼身整流黏性绕流数值模拟方法研究.采用Roe三阶迎风偏置通量差分分裂方法(FDS)和隐式近似因子分解方法(AF)对三维非定常雷诺平均N-S方程进行数值离散求解,采用多重网格技术加速流场求解收敛过程,并比较了不同湍流对翼身结合处流动分离的模拟能力.计算结果与风洞实验的对比表明,所采用的数值方法可以有效地模拟翼身结合处黏性分离流动,具有一定的工程实用价值.     

求解不可压缩流动的一种新算法

在Jameson有限体积法的基础上发展了一种模拟不可压流动的数值方法。该算法用四步Runge-Kutta时间推进法直接求解非线性的动量方程,充分利用了良好的非线性收敛特性,并从连续方程出发推导了求解压力的方程,解决了压力与速度的耦合问题。最后用该方法在同位网格上,对二维平板边界层流、二维顶盖驱动方腔层流及二维扩压器管道湍流进行了数值模拟。计算结果与理论值、标准解及试验数据进行了比较分析,从而证实了该方法对模拟不可压流动的适用性,并显示出了编程简单、易于推广、计算效率高等优点。     

二维不可压N-S方程二次四边形单元有限元解

粘性不可压流体问题是众多工程中重要的力学问题.数值求解Navier-Stokes方程会遇到两大困难:非线性和不可压性.针对二维不可压Navier-Stokes方程的特点,建立了以流函数为求解变量的四阶微分控制方程,有效地避免了处理涡量边界的难题.采用8节点二次四边形单元,单元基函数为2次非线性高阶函数,建立了求解二维不可压N-S方程的有限元方程,并自主开发了二次四边形单元有限元程序.数值实验结果验证了该方法的精确性和可靠性.因此,该方法在计算流体力学中有较好的应用前景.     

考虑黏性及压缩效应的旋涡运动特性

采用扩散速度法处理黏性作用项的离散涡法,数值模拟了旋涡在黏性不可压流场中的运动;推导了可压缩流中的旋涡运动方程,对不可压缩流中的Biot-Savart方程进行了可压缩修正,进而用离散涡方法求解出非定常、不稳定、可压缩流场中的旋涡运动。分析了胀量及黏性扩散作用对旋涡运动的影响,通过大量长时间行为的数值试验得出了一系列结果。结果表明：为了真实地反映旋涡在实际流体中的运动规律,必须考虑黏性及压缩性的影响。     

四阶紧致格式有限体积法湍流大涡模拟

为准确预测不可压复杂湍流,提出了一种应用于大涡模拟的高精度有限体积法。该方法在非交错网格上数值求解大涡模拟方程,空间离散采用有限体四阶紧致格式,时间推进用四阶Runge-Kutta法,压力速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动和振荡平板上方的流动,证实了该方法具有近四阶精度;并在此基础上,采用动力Smagorinsky亚格子应力模式,成功地完成了充分发展槽道湍流的大涡模拟计算,所得结果与直接数值模拟结果吻合良好。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

冷却模块气动噪声的数值研究

针对冷却模块噪声预测过程中散热器等部件是否可以忽略的问题,首先,建立等效的多孔介质流场模型,采用大涡模拟分析来流湍流及均匀性,并求解风扇周围的湍动能和涡量,通过涡声理论辨识主要噪声源。然后,建立声学边界元模型,对冷却模块的离散及宽频气动噪声进行三维声场预测与声压频谱分析。最后与试验进行对比分析。研究结果表明:散热器等部件对来流湍流及噪声源分布影响显著,在气动噪声计算过程中不可忽略。     

二维时间发展混合层噪声的直接数值模拟

为了考查不同波长的初始扰动对二维时间发展混合层流动与噪声的影响,采用直接数值模拟的方法进行了研究。采用多尺度展开技术从耦合双分布函数玻尔兹曼方程中恢复出质量、动量和能量守恒方程;融合特征边界条件与玻尔兹曼方程来满足边界的无反射特征;运用高精度的时空离散格式求解玻尔兹曼方程并直接模拟了5种不同工况下的时间发展混合层。结果表明:高频扰动主要控制混合层中旋涡的对并和缠绕,低频扰动主要控制旋涡的迁移和旋转;涡对的演化过程类似于一个旋转的四极子源模型;扰动能以平面波的形式向远场辐射并且以当地声速传播。该结果可为混合层湍流机理研究提供参考。     

基于k-ε模型的剪切流中湍流动能衰减规律研究

分别从理论推导和数值模拟的角度,分析了剪切流中湍流动能的衰减规律.通过简化湍流变量输运方程,采用Matlab求解了使用标准k-ε模型计算剪切流时湍动能在不同高度处的数值解,通过fluent数值模拟,验证了理论数值解,并与均匀流中的数值模拟结果进行对比,分析了在剪切流中入口湍流强度和湍流黏性比、来流速度梯度和输运方程扩散项等对湍流变量衰减的影响规律.结果表明:剪切流中存在湍流动能衰减;入口湍流强度对湍流变量输运方程生成项的影响较大,湍流黏性比对湍流变量输运方程生成项的影响较小;当入口湍流强度相同时,入口湍流黏性比越大,湍动能衰减的越快,来流速度梯度对湍动能衰减的影响越小;当入口湍流黏性比相同时,入口湍流强度越大,湍动能衰减的越慢,来流速度梯度对湍动能衰减的影响越大.     

空调室外机气动噪声的数值分析

采用雷诺平均数值模拟方法对空调室外机内部的三维黏性非定常流动进行了计算,获取了用于噪声分析的气动声源。在此基础上,依据Lowson方程对脉动压力产生的气动声场进行了仿真,预测了空调室外机各部位气动声源对声场的贡献,并对锯齿轴流风叶的降噪效果进行了分析。为了识别偶极子声源,对固体壁面处的静压脉动进行了分析。结果表明:静压脉动的相位一致是采用偶极子声源强度进行声源识别的前提条件。室外机噪声数值预测结果与试验在定性上吻合较好,说明噪声数值分析方法具有良好的工程应用价值。     

基于非定常流场的离心风机气动噪声分析

提出了一种不直接求解声场却能为离心风机降噪提供有用信息的分析方法.首先,利用有限容积法对风机内部的非定常流场进行计算.然后,采用时域和频域分析方法对流场内静压脉动的强度和频率进行分析.最后,根据声学基本理论,判定风机内部主要气动噪声源的位置及噪声类型.应用该方法对某离心风机进行了计算,并将分析计算结果与该风机的噪声测量结果进行了对比,证明该方法能够有效地判断气动噪声源的位置和噪声类型.     

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

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

轴向流中平行简支弹性薄板的大挠度流固耦合的数值模拟

分析研究了轴向流中简支弹性薄板大挠度流固耦合系统的振动响应和流场特性.板结构动力学方程采用基于位移的有限元法离散;流场采用二维不可压缩粘性流体N-S方程,并用有限体积法离散;在此基础上结合动网格控制技术,建立模拟双向流固耦合作用下轴向流中简支弹性薄板的二维数值模型.利用该数值模型得到了单块简支板随流速变化流致振动特性,研究了结构大挠度的振动稳定性,分别得到了Pitchfork分岔曲线和非线性系统结构的Hopf分叉曲线.通过轴向流恒定流速下不同间距的平行两块简支弹性薄板流固耦合的数值模拟得到了的流致振动特性.     

基于弱可压缩方法的水中声流数值分析

为研究水中的非线性声波传播问题,本文主要开展了基于弱可压缩方程的理论分析及水中声波动及声流模拟工作。首先,采用正压流体的密度与压力函数,结合粘性流体动力学方程,对方程组采用特征变量的波动分析,证明了该方程能描述水中的涡旋和声波传递。然后,本文结合时空守恒元解元算法求解该方程组,并对一些典型气动声学问题进行了计算和对比研究：通过计算分析了二维顶盖半圆腔驱动流,本文计算结果与其他研究结论进行对比分析,证实了该方法可以准确模拟出流体的粘性;通过对平面声波在不同粘性流体中传播进行了数值模拟,声波衰减趋势与理论解完全一致;最后,本文还通过对Rayleigh声流问题进行非定常数值模拟,得到的模拟结果与其他研究学者结论完全一致。这说明,采用二维弱可压缩方程组及时空守恒元解元方法,能够模拟水中的声波传播及耗散过程。     

喷流与超声速主流粘性干扰流场的数值模拟

采用二维雷诺平均Navier-Stokes(N-S)方程时具有声速喷流的超声速主流干扰流场进行了数值求解。N－S方程的无粘项采用有限体积法在非结构网格上进行空间离散，通量分裂采用二阶精度的Vanleer格式，紊流粘性系数采用Bald-win-Lomax二层代数紊流模型计算。对当地时间步长进行了粘性修正，并发展了一种基于线化流量的逆风隐式格式以提高求解效率。     

采用非结构网格的IDEAL算法实现及定压边界条件处理

在非结构网格系统中采用有限容积法对不可压N-S方程进行了离散,其中对流项采用二阶迎风格式,扩散项采用超松弛校正格式.通过修正的动量插值格式构造了压力方程,利用IDEAL算法处理压力与速度的耦合.利用Taylor展式导出了邻近边界内节点的压力计算公式,进而将边界静压引入到压力方程的迭代计算中,确保了边界信息向内部区域的有效传递.根据质量守恒定律分别导出了压力入口和出口边界的速度计算公式,结合伯努利定律给出了总压边界的处理方法.通过压力驱动管内层流、三通管层流以及浮力驱动流等几个典型算例对文中方法进行了考核,并将计算结果与文献中的基准解和实验值进行了对比.结果表明,不管是静压边界还是总压边界,采用该处理方法均能获得令人满意的结果,尤其对于大空间自然对流问题,可以在保证精度的前提下节省大量的存储空间和运算时间.     

轴向流中两平行弹性薄板大挠度流固耦合系统的数值模拟

研究了轴向流中两平行弹性薄板大挠度流固耦合系统的振动响应和流场特性.板的结构动力学方程采用基于位移的有限元法离散,流场采用完全的二维不可压缩粘性流体N-S方程,用有限体积法离散,结合动网格控制技术,建立了模拟轴向流中两平行悬臂弹性薄板双向流固耦合作用的二维数值模型,实现了弹性薄板—流体—弹性薄板之间的交互作用.利用该数值模型得到了单块板的流致振动特性;发现了随着两板间距大小的不同,两板表现出三种不同的的极限环振动:同相位、异相位以及不确定相位,并得到了同相位与异相位两种振动状态下的板边界层脱落及尾流变化规律;另外,分析还发现在同相位和异相位过渡阶段,两板均表现出拍现象.     

致密油藏有限元数值模拟

 目前数值模拟方法均基于达西渗流,无法对致密油藏非达西渗流进行有效模拟。致密油藏内在复杂,为实现致密油藏非达西渗流的科学模拟,建立了油水2相微可压缩非达西流的油藏数值模拟数学模型。同时考虑到致密油藏渗透率低、非均质性强的特点并兼顾模拟精度采用有限单元法进行离散,从而建立了有限元数值模型,形成了新的致密油藏数值模拟方法,最后将该方法进行实例检验,验证了方法的正确性。     

轴向流中固支弹性薄板的大挠度流固耦合系统的数值模拟

就轴向流中两端固支大挠度弹性薄板的流固耦合振动特性,固支薄板的结构动力学方程用有限元法离散,流场采用不可压缩的二维粘性流体(N-S方程)用有限体积法离散,结合ADINA中的流体单元划分技术,建立了双向流固耦合作用下轴向流中两端固支薄板的二维仿真模型.通过模拟仿真分析研究了给定不同流速下固支板的流固耦合振动特征和大挠度系统的振动稳定性.分别得出了不同流速下固支板中点的挠度—流速曲线、挠度时程曲线及挠曲线图.结果表明:当流速小于固支板的临界流速时,板将处于稳定的直线平衡状态;当流速大于固支板的临界流速时,板将在新的位置达到弯曲平衡状态,以及在弯曲平衡位置附近发生极限环振动.