直方管内充分发展湍流的大涡模拟 第一部分

本文用大涡模拟的方法对直方管内Ｒｅ＝４．９×１０￣４的充分发展湍流运动进行数值模拟，所得结果清楚地显示了二次流的存在。统计平均量的数值结果也与实验结果吻合甚好，这表明大涡模拟是高雷诺数复杂湍流运动数值模拟的一个有效手段。     

洁净室流场大涡模拟

应用湍流大涡模拟方法对矢流洁净室和全顶棚送风下部两侧回风式洁净室进行了数值模拟。计算中使用了隐含流线抑风耗散作用的具有三阶部分展开的Taylor-Galerkin离散格式，导出了经Gauss滤波后的基本方程组的有限元列式，应用涡粘性亚格子模式的湍流大涡模拟方法进行数值模拟。数值计算结果与实验结果符合较好，表明湍流大涡模拟能够较准确地捕捉到洁净室流场中旋涡的位置和尺度，在洁净室的数值仿真方面，与标准的K－ε湍流模式相比，湍流大涡模拟结果具有较高的可靠性。     

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

平直槽道内湍流运动的大涡模拟

为在中型机上实现对湍流运动的大涡模拟,采用Moin的基本模型,在数值计算中采用了诸如高斯滤波函数、速度压力交错网格、谱方法以及一种简单的压力计算的修正等方法和技巧,计算结果与Moin结果比较,表明在我国现有条件下也有可能进行大涡模拟方面的工作。     

用于可压缩壁面湍流DNS及LES入口湍流边界条件的回收-调节方法综述

很多湍流流动问题的大涡模拟和直接数值模拟的入口湍流脉动边界条件对于计算结果的影响很大,本文叙述了用于可压缩壁面湍流入口脉动生成的一种方法——“回收－调节”方法（recycling-rescaling method）,叙述了这种方法的理论基础及简要推导过程,介绍了几种常用的“回收－调节”方法并指出了实现中的几个关键问题,分析了这种方法中存在的不足并对其应用进行了总结评述。本文可为超声速流动、燃烧问题的大涡模拟和直接数值模拟方法的发展提供参考。     

低速粘性分离流的数值模拟

本文介绍对一种采用涡控制技术的新概念兵器--环翼以及升力体航天飞机的低速气动力数值模拟.对这类外形复杂、流动状态也十分复杂的流场,需正确模拟出流场不同部分的流态及湍流运动参数.本文采用一种主要与流场平均运动涡量相联系的新的湍流模型,较好地解决了这一问题.数值计算给出了大量满足工程要求的气动参数,并解决了一些关键性气动问题.     

圆截面U形管流有限元模拟

弯曲圆管内流动具有很强的工业应用背景,管内二次涡的大小、位置及发生的频率会对流动系统的阻力及热质交换产生重大作用,研究弯曲圆管内的流动具有重要的理论意义.本文应用湍流大涡模拟方法对弯道流动进行了数值模拟.用时间推进法,给出了隐含流线迎风耗散作用的二阶完全展开Taylor-Galerkin有限元离散格式,数值计算的结果与实验结果符合较好.表明大涡模拟方法适用于边界形状复杂,存在各向异性的大尺度涡的流动仿真.     

固体火箭内流场的时均和大涡模拟技术对比研究

分别采用时均RNG湍流模拟和大涡模拟(LES)方法对某型固体火箭发动机内流场进行了瞬态数值模拟研究,建模过程中针对推进剂燃烧以及燃面边界采用了修正壁面通量的方法,并使用了动网格技术。结果表明:大涡模拟技术可以得到发动机内湍流流场中详细拟序结构,时均湍流模拟对于火箭发动机内流场时均量的模拟仍有较好的工程精度。     

三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

非均匀网格湍流大涡模拟高精度有限体积解法

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

风沙两相流数值模拟研究进展

大气表面层中的风沙运动是典型的高雷诺数湍流气固两相流动,同时具有可蚀地表上的冲击溅起过程、风-沙-电多场耦合、多尺度与跨尺度等复杂特性,使其定量预测极具挑战性。回顾了近30年风沙运动数值模拟的研究进展以及取得的成果,主要包括基于湍流雷诺平均Navier-Stokes （RANS）方法的风沙运动多场耦合模拟与风成地貌耦合尺度模型建立,基于湍流大涡模拟（LES）与高雷诺数湍流理论的风沙运动模拟方法。最后明确了风沙运动现有数值模拟方法的不足和值得进一步研究的方向。     

空间发展混合流的大涡模拟

采用大涡模拟方法数值模拟了空间发展混合流，通过在平均流中引入小扰动，研究了混合流大涡结构的产生和演化过程，捕捉了展向涡的卷起、配对、合并，以及二次流向涡的出现等大尺度的三维拟序结构。与实验及直接数值模拟的结果相比，吻合程度较好，表明本文所采用的LES方法对湍流的数值模拟是切实可行的。     

可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟(英文)

阐述了建立在一方程湍流模型基础上的分离涡模型(DES)和延迟脱体涡模型(DDES),它们很好地模拟了壁面附近的小尺度流动。该方法在大攻角分离流区域,使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型,研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。研究结果说明了方法对模拟分离流动是可行的。     

颗粒两相流自由剪切层大涡模拟

采用大涡模拟(LES)方法数值仿真了可压缩流自由剪切层内颗粒两相流流场结构.气相控制方程采用二维盒式滤波法滤除小于网格尺度的小涡以得到LES控制方程,而流场中小涡对大涡的影响采用K-ωSST湍流模型模拟.在Lagrange坐标系下用颗粒轨道模型仿真颗粒运动,用双线性插值方法实现Euler坐标系到肌Lagrange坐标系的转换,并考虑了颗粒速度和温度的改变对气相流场的反作用.分析了仿真结果中气相自由剪切层的涡结构特性,不同St数颗粒在流场中的运动特性,颗粒对流场的反作用效果,得到了合理的结论.     

舰船空气流场数值模拟方法研究

分离涡数值模拟方法(DES)在物面附近求解雷诺平均Navier-Stokes方程,在其他区域采用大涡模拟方法,兼具前者计算量小的优点和后者能模拟大分离湍流流动、计算精度高的优势。论文分别采用RANS和DES方法对舰船空气流场进行数值模拟,通过所得流场结果与实测数据的对比分析,说明DES方法对大分离流动具有更强的模拟能力,更加适用于舰船空气流场预报。     

平板湍流转捩过程中Λ涡和环状涡的周围流场研究

采用高精度湍流直接数值模拟技术,对平板流动转捩过程中的涡系结构,尤其是Λ涡和环状涡附近的流场进行了详细的研究.高分辨率的数值模拟给出了实验手段无法得到的全流场信息,以及雷诺平均纳维斯托克斯方程方法(RANS)无法给出的精细涡结构信息.数值模拟在精确捕捉Λ涡和环状涡等转捩过程中典型涡系结构、上喷和下扫等典型流动现象的同时,发现了一些与传统认识所不同的现象:1)Λ涡并不是当地涡量集中的地方;2)环状涡的附近存在一个速度显著大于远场来流速度的区域.这些现象对传统湍流转捩理论提出了挑战,并引发人们对湍流生成与维持机理的进一步研究.     

考虑湍流影响的前缘分离涡运动与破碎特性的研究

本文以Ｎ－Ｓ方程和Ｈａｌｌ涡核模型假设为基础，导出了描述湍流的涡核运动方程。利用差分计算方法，对三角翼前缘分离涡运动及其破碎特性进行了数值求解，分析了涡核流场的结构和各物理量的变化特性，反映了三角翼前缘分离涡运动及其破碎特性的实质，说明了湍流涡核方程能更有地模拟三角翼前缘分离涡运动和准确地确定涡核破碎位置，并得出了攻钐雷诺数对破碎位置的影响，计算结果与实验数据十分吻合。     

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

本文采用直接过滤的Navier-Stokes方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似,但仅保留级数的一阶导数项。这样产生的误差相当于丢失了模型的部分耗散作用,本文用一种动力学非线性亚格子模型来补偿丢失的耗散影响。本文根据Caylay-Hamilton定理导出了一个非线性亚格子模型,模型中包含的系数由动力学模式确定。与传统的动力学Smagorinsky模式相比,这种动力学非线性模型(DNM)稳定性更好,且不需要在统计均匀方向进行统计平均来计算模型系数,因此减少了计算开销。本文用这种非线性亚格子模型对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验值、计算值及理论值进行了对比。结果表明,本文模型可以有效地模拟可压缩湍流流场。     

支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。