壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。     

用一种非线性k-ε模型计算突扩台阶后的流动

对一种新的非线性k-ε模型的无量纲系数进行了重新拟合.用新的系数计算了槽道流动、方截面管流和等应变率,平均流场无旋流动的正应力分布.计算结果与实验资料符合很好.将该模型并入常规的k-ε模型的计算程序中,并对突扩台阶后的湍流流动进行了数值计算,同线性k-ε模型的计算结果和实验结果进行了比较.结果表明,这种非线性k-ε模型同线性模型相比有较大的改进,同已有的其它非线性模型相比,也有一定的优点.     

典型城市街道内汽车排放污染物扩散规律的模拟研究

采用自主开发的街道峡谷内污染物扩散模式,对城市中几种典型的街道峡谷内流场和污染物的浓度场进行了模拟研究。模拟结果表明:在一定气象条件下,街道中机动车排放污染物的扩散受街道宽高比、街道两侧建筑物对称性、街道两侧建筑物高度分布及街道两侧的建筑物形状等因素的影响;街道峡谷宽高比接近1时、递升型峡谷以及宽阔街道有利于污染物的扩散。     

建筑物对污染物扩散影响的数值与风洞模拟研究

采用k-ε（RNG）与RSM湍流模型对处于方形建筑物不同位置污染源所排放污染物的扩散规律以及建筑物对流场的影响进行模拟,并且与相应的风洞试验结果进行了比较。流场分析结果表明：数值模拟能够较好地模拟建筑物前方迎风侧停滞回流、顶部回流以及后方空腔区等。浓度场分析结果表明：建筑物前方迎风侧以及顶部回流区污染物扩散的数值模拟结果与风洞试验结果基本相等,而在建筑物后方空腔区污染物的数值模拟结果略高于风洞试验结果。综合分析并与风洞试验结果相比较,RSM模型能够较好地模拟污染物浓度场以及建筑物周围流场的变化规律。     

预处理方法在低速流场数值模拟中的应用

自20世纪70年代以来，随着高分辨率差分格式和高计算效率时间推进格式的建立和发展，可压缩流动数值模拟取得了很大进展。并广泛采用时间推进法。但是。对于低速不可压流动，采用时间推进法效果并不太好。目前对低速不可压流动的数值模拟主要有三种途径：（1）基于压力的算法，如著名的SIMPLE算法（压力偶合方程组的半隐式方法）；     

高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

基于数值模拟的复杂地形风场风资源评估方法

与采用线性模型的WAsP软件相比,CFD方法具有可以逼真模拟复杂地形三维风场大气流动的优势,因此成为了复杂地形风场风资源评估的发展方向。为了把CFD软件包FINE/TURBO和带有壁面函数的k-ε湍流模型用于复杂地形地貌风电场的风资源评估和微观选址中去,本文以南澳岛风场作为研究对象,在中性大气条件假设下,对风场区域以30°风向为间隔进行数值模拟。研究采用不同入口边界条件对数值计算结果的影响。比对模拟结果发现,风场中的风加速与湍流强度基本上不随来流条件发生变化。基于这一规律,结合全年测风数据评估风场的风资源分布状况,建立了一套风资源评估的数值模拟及后处理方法,为工程实际提供了有力参考。     

混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

基于浸入式边界方法的二维结冰机翼的数值模拟

应用浸入式边界方法模拟干净/结冰NACA0012翼形绕流问题.采用SIMPLEC方法离散求解二维不可压雷诺平均Navier-Stokes方程, 湍流模型采用一方程Spalart-Allmaras模型.计算结果与以往的实验结果吻合得很好, 证明了浸入式边界方法和所发展的程序的可靠性.      

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

建筑物顶部形状对街道峡谷内污染物扩散影响的研究

应用数值模拟方法研究建筑物顶部形状对于街道峡谷内汽车排放污染物扩散的影响.研究发现,街道峡谷内漩涡的结构会随街道两侧建筑物顶部形状的不同发生变化,进而影响到街道峡谷内的污染物扩散.数值模拟结果与Rafailidis(1996)等在德国汉堡大学的气象学院进行的风洞实验结果进行了对比,与风洞实验结果相符合.街道两侧建筑物顶部形状不同,峡谷内的污染物扩散规律不同,选择合适的建筑物顶部形状会对街道峡谷内的污染物扩散产生有益的效应,这对于建立有效的机动车污染物控制策略具有重要指导意义.     

重力作用下的粒子扩散模型及数值模拟

建立了粒子在重力、流体阻力和随机力作用下的烟幕扩散模型,测试了某种纳米碳纤维材料在烟幕箱里的扩散,利用线性回归分析估计了测试光路上的质量浓度,计算了烟幕的沉降速率。分别用本文建立的模型和随机游走模型模拟了该材料在烟幕箱里的扩散过程,模拟计算得到的质量浓度和沉降速率表明,本文建立的模型所模拟的结果更接近于实际测量结果。     

隧道内危险性介质泄漏扩散的数值模拟

化学危险品的意外泄漏已经造成很多事故,尤其是在隧道这样相对封闭的空间内,事故发生的概率更大。为了有效地控制和预防事故的发生,本文以RANS方程为控制方程,采用基于simple算法的有限体积法,对隧道内危险性介质泄漏扩散过程进行了数值模拟。讨论隧道内危险品(乙炔)运载车在行驶过程中发生小孔泄漏扩散的情况,计算得到了危险性介质浓度的时间和空间分布规律。同时,分析了运载车不同行驶速度和危险品泄漏释放方式(水平释放和垂直释放)对其扩散过程的影响,为进一步了解危险性介质的扩散行为和控制方式提供了理论基础。     

复杂建筑物对近场扩散影响的数值与风洞模拟的比较分析

采用CFD数值模拟技术与风洞试验相结合的手段模拟了复杂建筑物群对周围流动与污染物扩散的影响。CFD技术采用k-ε（RNG）湍流模型模拟了建筑物群对流动与扩散的影响,风洞试验通过采用多种探测手段,分析建筑物群周围的流动和弥散,并将统计学方法应用于风洞试验结果,验证了CFD数值模拟结果的合理性。研究结果表明,CFD技术能较好地模拟建筑物群对污染物弥散的影响,并与风洞试验吻合,有、无建筑物B的情况下,吻合因子FAC2与FAC5均大于50％,归一化均方误差（NMSE）均小于4,部分偏差（｜FB｜）均小于0．3;建筑物对污染物弥散的影响非常复杂,建筑物尾流对污染物浓度分布影响较为显著,建筑物后方空腔区污染物的数值模拟结果略高于风洞试验结果。综合分析表明,风洞试验与数值模拟相结合的方法是研究这类微小尺度湍流扩散问题的有效手段。     

数值模拟中的边界处理方法研究

为适应于模拟超声速流场中的一些特殊物理现象,研究了不同物理现象对边界处理的要求。针对激波、旋涡和声波等流场中典型的流动现象,采用常用的外推、特征和缓冲三类边界处理方法分别进行数值模拟。计算结果表明超声速流场边界采用外推技术,只适于计算带有激波的边界。特征方法还能适于旋涡运动的计算。声波的计算只能采用缓冲层技术。     

二维爆震波数值模拟

针对爆震发动机的多维效应,利用CFD计算软件对单管脉冲爆震发动机进行了二维的数值模拟。以氢气和氧气混合物为例,仿真了爆震波在管内的传播过程以及传出管外后的流场情况。模拟所得的爆震波参数与文献中的计算结果和实验值基本一致,验证了模拟的正确性。爆震波传出爆震管后退化成一道激波,对该激波的流动情况进行了分析。     

裂纹扩展的无网格数值模拟方法

疲劳断裂是航空材料的重要失效形式 ,由于裂纹尖端应力存在奇异性 ,传统有限元方法模拟裂纹沿任意路径扩展存在很多不足。作为一种新兴的数值模拟方法 ,无网格计算只需将求解问题离散为独立的节点 ,计算过程中可以实时跟踪裂纹尖端区域进行局布细化。将连续的裂纹扩展过程看作多个线性增量 ,每一个增量内裂纹扩展角根据应力强度因子确定 ,通过在裂纹尖端细化节点和引入外部基函数提高了计算精度。本文给出了应用无网格方法模拟裂纹扩展过程的关键技术和计算流程 ,通过对带有中心斜裂纹的 Ti-6 Al-4 V合金平板进行分析 ,预测得到的裂纹扩展路径与实验值吻合的较好。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用.与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要.因此,基于两步四阶时间离散和高精度WE-NO重构,发展了具有四阶时间精度的气体动理学格式.该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力.目前...