计算风工程中k-ε模型的一类边界条件

给定适当的湍流来流边界条件是计算风工程中一个重要的问题。然而通常情况下，由于问题的复杂性，很难给定一个既满足边界层自保持要求又与试验相符的边界条件。本文首先从模型方程本身出发，推导出一类近似满足k-ε模型自保持边界条件要求的湍动能表达式。此表达式与有关文献建议的湍流边界条件比较，可以考虑湍动能随高度的变化，因而更符合实测规律；然后基于风洞试验数据，定义了一组k-ε模型中的模型常数；最后通过模拟风洞试验的算例，验证并分析了此类湍动能边界条件配合试验常数，在标准k-ε模型中模拟均匀大气边界层的适用性。     

基于标准k-ε模型的平衡大气边界层模拟

基于标准k-ε模型和FLUENT软件,采用适当的平均风和湍流特征量入口条件,同时在壁面处采用修正壁面函数边界条件,进行空流场的平衡大气边界层模拟.将流场中不同位置处的平均风廓线、湍动能及其粘性耗散等流动特征量的廓线进行对比.结果表明,采用标准k-ε模型结合修正壁面函数对气流流经粗糙度相同的壁面时,不同位置处流动特征量的廓线吻合程度均较高.研究提出对建筑绕流等问题模拟之前需进行空流场的平衡性进行分析,另对文中尚未解决的问题进行讨论并提出解决思路.     

计算风工程中κ-ε模型的一类边界条件

给定适当的湍流来流边界条件是计算风工程中一个重要的问题.然而通常情况下,由于问题的复杂性,很难给定一个既满足边界层自保持要求又与试验相符的边界条件.本文首先从模型方程本身出发,推导出一类近似满足κ-ε模型自保持边界条件要求的湍动能表达式.此表达式与有关文献建议的湍流边界条件比较,可以考虑湍动能随高度的变化,因而更符合实测规律;然后基于风洞试验数据,定义了一组κ-ε模型中的模型常数;最后通过模拟风洞试验的算例,验证并分析了此类湍动能边界条件配合试验常数,在标准κ-ε模型中模拟均匀大气边界层的适用性.     

Sajben扩压器复杂流动中湍流模型的性能评估

用8个常用的湍流模型对Sajben扩压器中跨声速流动进行了数值模拟,评估了Spalart-Allmaras, 标准k-ε, RNG (re-normalization group) k-ε,realizable k-ε,标准k-ω,SST(shear stress transport) k-ω,v2-f,Reynolds stress共8个湍流模型对激波/湍流边界层相互作用的模拟预测能力.通过与实验数据比较发现:SST k-ω模型和v2-f模型比其他模型模拟的更准确,其中SST k-ω模型比v2-f更能准确地预测壁面压力,然而对于分离点、再附点以及分离区长度v2-f比SST k-ω预测得更准确.      

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

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

粗糙壁面表面传热数值计算比较

为了计算冰形模拟中的表面传热系数,给出了三种粗糙壁面表面传热系数计算方法,即传统的边界层积分法、基于k-ε湍流模型及标准壁面函数的计算方法(k-ε法)和基于Spalart-Allmaras湍流模型在粗糙壁面扩展的计算法(S-A扩展法).计算了直径为15cm、表面等效沙粒粗糙高度为1.35mm的圆柱,在Re分别为8.8×...     

应用修正的k-ε模型研究超声速H2/Air燃烧

为了研究超声速燃烧中流体可压缩性的影响,对标准k-ε湍流模型进行可压缩性修正(包括结构可压缩性修正和膨胀可压缩性修正两部分)。分别应用标准k-ε模型、修正的k-ε模型和雷诺应力模型(RSM),考虑氢气/空气详细化学反应机理(GR I-M ech 2.11机理,10组分,28基元反应),数值模拟有壁面限制的超声速混合层冷态及热态流场。结果表明:壁面和燃烧对湍流影响都很大;修正模型对冷态以及燃烧场的预测结果优于其它两个;修正模型预测的混合层厚度更薄,燃烧区域更窄,与实验结果吻合地更好。     

收敛-扩张喷管中运用激波矢量控制法的计算研究

采用二阶TVD格式的有限体积法耦合RNG k-ε湍流模型和非平衡壁面函数求解二维守恒型雷诺平均N-S方程，数值模拟了运用激波矢量控制法的二维收敛一扩张喷管中的流动现象，很好的模拟出了由于次流喷射所产生的激波及上下壁面的压力分布，并将计算结果与国外的实验结果进行了对比。     

风对建筑绕流流动的数值方法

本文利用k-ε湍流模型及SIMPLE方法,在二维条件下,就风对建筑绕流流动进行了数值模拟,建筑及地面的边界条件采用以低Re数k-ε湍流模型中导出的壁函数处理近壁湍流。自由边界采用外推格式。对于计算域中的建筑物则利用控制粘性系数的方法处理。并对计算结果和风洞实验进行了比较。 为了研究建筑物之间的气流及压力分布,本文还对并列两栋建筑的不同情况进行了数值计算和分析。     

障碍物导致火焰变形的数值模拟

工业部大量存在着障碍物。实践证明，障碍物对加速燃烧和诱导爆炸具有非常重要的影响。为了深入研究障碍物与火焰之间的关系本文基于湍流的k-ε模型和改进的EBU－Arrhenius燃烧模型，考虑了障碍物对流动的附加作用，通过修改方程的源项，建立了湍流加速火焰使火焰形状改变的理论模型。本文选用Simple格式，壁面边界层采用壁面函数法处理，模拟了火焰在含有障碍物的三维空间中的传播现象，其数值模拟得到的火焰形状与实验结果做了对比，结果表明本计算是成功的。最后，还分析了导致火焰变形的原因。     

考虑可压缩与热传导的壁面函数边界条件及其应用

考虑到可压缩和热传导效应的壁面函数边界条件,被耦合到了采用k-ω两方程湍流模型、用有限体积法求解N-S方程的程序中。壁面函数基于耦合的速度和温度型,并且在边界层内的粘性子区和对数区内一致有效。引入壁面函数边界条件后,通过算例验证在y <100的范围内,得到的物面压力、摩阻、热流与实验结果比较,结果可靠。而无壁面函数边界条件时,要得到相同精度的结果,要求y ≈1。壁面函数的引入,为工程上准确预测飞行器在湍流流动中表面受力与气动热提供了保障。     

RNG k-ε模型数值模拟油雾燃烧流场

采用RNG(R enorm alization G roup)k-ε紊流模型对某环形燃烧室火焰筒内三维两相燃烧流场进行数值模拟。运用偏微分方程法和区域法生成贴体网格,EBU-A rrhen ius紊流燃烧模型估算化学反应速率,六通量热辐射模型估算辐射通量。在非交错网格下求解气相采用S IM PLE算法,液相采用颗粒轨道模型和PS IC算法。近壁区处理分别采用两层和三层壁面函数。计算数据与试验值比较表明,采用三层壁面函数的RNG k-ε模型更适用于模拟三维两相燃烧流场。      

液体煤油超声速燃烧数值分析的湍流模型

利用计算流体力学软件对煤油在所设计的双模态燃烧室内的喷雾燃烧进行了数值模拟。采用k-ε三种形式、k-ω两种形式和雷诺应力湍流模型,计算了在飞行马赫数为5,煤油与空气的当量比为0.551时的情况。通过将数值模拟结果与实验数据进行对比,证明这些湍流模型都可用于煤油超声速燃烧的数值计算。并且,k-εRNG模型和k-ωSST模型的计算结果与实验数据最为接近,可以认为这两个模型在模拟煤油超声速燃烧流场中是非常有效的;而雷诺应力模型效果较差。     

两种差分格式和两种湍流模型在轴对称冲击射流数值计算中的比较

用延迟修正QUICK和乘方两种差分格式，结合标准k-ε模型和一种基于重正化群（RNG）思想的k-ε模型，对半封闭轴对称涡流冲击射流场进行了数值模拟。通过与实验结果的比较，考察了两种差分格式对计算结果的影响，并讨论了两种湍流模型对冲击射流场的数值预测能力。研究表明：当都采用标准k-ε模型时，QUICK格式得到的结果比乘方格式更接近于实验值，尤其对湍流能分布；当都采用QUICK格式时，RNG k-ε模型的预测值比标准k-ε模型更准确。     

基于气体射流的气液两相流动减阻特性

为了明确射流表面减阻特性,基于气液两相流动原理,建立具有射流功能的气液两相流动减阻计算模型,采用数值计算方法,通过RNG k-ε湍流模型对射流表面不同射流方向角、以及不同气体射流速度情况下的流场进行数值模拟及实验验证,分析气体射流对壁面黏性阻力、减阻率及流场的影响,研究气体射流的气液两相的减阻特性。结果表明:气体射流表面能够减小壁面黏性阻力,具有较好的减阻效果;四种射流方向角下,射流表面减阻率均在气体射流速度为5m/s时达到极大值,并且当射流方向角在xoz平面为-30°时,减阻率最大,为7.18%;数值模拟的减阻率效果略高于实验结果;气体射流的气液两相流动对壁面边界层进行了有效控制,减小了壁面的黏性剪应力和雷诺应力,导致壁面的黏性阻力减小,使得射流表面具有减阻性能。     

液体火箭发动机推力室冷却通道流动与传热数值研究

采用气固耦合算法对液体火箭发动机推力室再生冷却通道的流动与传热过程进行了三维湍流流动与传热数值模拟，冷却工质为氢气，其密度、导热系数、动力粘度随着温度和压力而变化。应用大涡模拟及标准k-ε双方程模型两种湍流模型分别进行数值模拟，详细揭示了再生冷却通道固体区和流体区内的速度场和温度场，并在不同的计算网格数目下对两种湍流模型的计算结果进行了对比。结果表明，在相同的网格条件下，标准k-ε双方程模型与实验数据的吻合精度比大涡模拟模型更好，且满足工程计算精度。随着网格数的增加，大涡模拟的计算精度逐渐得到改善。     

用湍流模型研究Richtmyer-Meshkov不稳定性诱导的湍流混合

将k-ε模型应用到存在激波的Richtmyer-Meshkov不稳定性诱导混合的研究中,采用代数关系式补充模型的封闭假设,给出了合理描述RMI特征的湍动能产生项的表达式。采用该模型对两个激波管实验进行了模拟。通过与实验和文献中的直接模拟结果比较分析,验证了所采用的模型封闭、模型常数、数值算法、程序实现是合适的。     

三级压缩楔前缘半径对壁面静压及热流分布影响研究

为了研究三级压缩楔前缘半径对壁面静压和热流分布的影响,在0.6m激波风洞上开展了试验测量,模型长约0.6m前缘半径为0mm3mm试验名义马赫数为5.98。研究结果表明,试验得到的壁面静压和热流数据重复性很好。采用Fluent软件进行了二维和三维流场参数辅助模拟分析,三种不同湍流模型获得的壁面静压分布差别不大,均与试验结果吻合较好;不同湍流模型获得的壁面热流分布差异较大,采用标准k-ε模型得到的结果与试验吻合较好。试验和数值模拟结果均表明,在第二和第三压缩面上,经过激波后,壁面静压逐渐上升到一个压强平台;壁面热流逐渐上升到一个局部极大值,然后在同一压缩面持续下降。随前缘半径增加,壁面静压和壁面热流整体减小,压强平台值和热流局部极大值也减小,而达到压强平台和热流局部极大值需要的长度增加,显示激波边界层干扰影响区域增大。     

气氧/甲烷涡流冷壁燃烧室流场与壁面耦合传热分析

针对涡流燃烧室的试验模型,耦合计算了其中的燃烧、流动和传热过程。湍流模型采用RNG k-ε两方程模型,辐射传热采用P1辐射模型,化学反应采用有限速率模型。数值模拟了涡流燃烧室在20s内的传热过程,得到了燃烧室和喷管的壁面温度分布随时间的变化。燃烧室侧壁面和头部温度在4秒内就达到较稳定状态,在涡流保护下侧壁面最高温度在650K左右,头部最高温度在785K左右。对于热沉喷管,壁面温度随时间基本呈缓慢线性上升,在一定的时间内可以满足实验要求。     

影响冷气掺混数值模拟精度的若干问题分析

对冷气掺混数值模拟中的计算网格、湍流模型、射流边界条件等影响精度的若干问题进行分析。在吹风比0.5下,对不同网格分布、湍流模型、射流边界条件进行数值模拟,所得结果与Ajersch的实验结果进行对比。结果表明,出口下游的网格分布可适当稀疏以减少计算量;在所研究的几种湍流模型中,k-ε模型所得结果与实验结果吻合得最好;考虑射流通道内流动能提高精度,在射流出口给定流场分布也可保证计算结果的精度。