文档库 最新最全的文档下载
当前位置:文档库 › 近10年晴空湍流的研究进展

近10年晴空湍流的研究进展

近10年晴空湍流的研究进展
近10年晴空湍流的研究进展

第三章,湍流模型

第三章,湍流模型 第一节, 前言 湍流流动模型很多,但大致可以归纳为以下三类: 第一类是湍流输运系数模型,是Boussinesq 于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即: 2 1 21 x u u u t ??=-μρ 3-1 推广到三维问题,若用笛卡儿张量表示,即有: ij i j j i t j i k x u x u u u δρμρ32 -??? ? ????+ ??=''- 3-2 模型的任务就是给出计算湍流粘性系数t μ的方法。根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。 第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。 第三类是大涡模拟。前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes 方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。 实际求解中,选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。 FLUENT 提供的湍流模型包括:单方程(Spalart-Allmaras )模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。 湍流模型种类示意图 第二节,平均量输运方程 包含更多 物理机理 每次迭代 计算量增加 提的模型选 RANS-based models

雷诺平均就是把Navier-Stokes 方程中的瞬时变量分解成平均量和脉动量两部分。对于速度,有: i i i u u u '+= 3-3 其中,i u 和i u '分别是平均速度和脉动速度(i=1,2,3) 类似地,对于压力等其它标量,我们也有: φφφ'+= 3-4 其中,φ表示标量,如压力、能量、组分浓度等。 把上面的表达式代入瞬时的连续与动量方程,并取平均(去掉平均速度i u 上的横线),我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式: 0)(=?? +??i i u x t ρρ 3-5 () j i j l l ij i j j i j i i u u x x u x u x u x x p Dt Du -?? +???????????? ????-??+????+??-=ρδμρ32 3-6 上面两个方程称为雷诺平均的Navier-Stokes (RANS )方程。他们和瞬时Navier-Stokes 方程有相同的形式,只是速度或其它求解变量变成了时间平均量。额外多出来的项j i u u ''-ρ是雷诺应力,表示湍流的影响。如果要求解该方程,必须模拟该项以封闭方程。 如果密度是变化的流动过程如燃烧问题,我们可以用法夫雷(Favre )平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外,所有变量都用密度加权平均。变量的密度加权平均定义为: ρρ/~ Φ=Φ 3-7 符号~表示密度加权平均;对应于密度加权平均值的脉动值用Φ''表示,即有: Φ''+Φ=Φ~ 。很显然,这种脉动值的简单平均值不为零,但它的密度加权平均值等于零,即: 0≠Φ'', 0=Φ''ρ Boussinesq 近似与雷诺应力输运模型 为了封闭方程,必须对额外项雷诺应力j i u u -ρ进行模拟。一个通常的方法是应用Boussinesq 假设,认为雷诺应力与平均速度梯度成正比,即: ij i i t i j j i t j i x u k x u x u u u δμρμρ)(32 ??+-??? ? ????+??=''- 3-8 Boussinesq 假设被用于Spalart-Allmaras 单方程模型和ε-k 双方程模型。Boussinesq 近似 的好处是与求解湍流粘性系数有关的计算时间比较少,例如在Spalart-Allmaras 单方程模型中,只多求解一个表示湍流粘性的输运方程;在ε-k 双方程模型中,只需多求解湍动能k 和耗散率ε两个方程,湍流粘性系数用湍动能k 和耗散率ε的函数。Boussinesq 假设的缺点是认为湍流粘性系数t μ是各向同性标量,对一些复杂流动该条件并不是严格成立,所以具有其应用限制性。

湍流流动的近壁处理详解

壁面对湍流有明显影响。在很靠近壁面的地方,粘性阻尼减少了切向速度脉动,壁面也阻止了法向的速度脉动。离开壁面稍微远点的地方,由于平均速度梯度的增加,湍动能产生迅速变大,因而湍流增强。因此近壁的处理明显影响数值模拟的结果,因为壁面是涡量和湍流的主要来源。 实验研究表明,近壁区域可以分为三层,最近壁面的地方被称为粘性底层,流动是层流状态,分子粘性对于动量、热量和质量输运起到决定作用。外区域成为完全湍流层,湍流起决定作用。在完全湍流与层流底层之间底区域为混合区域(Blending region),该区域内分子粘性与湍流都起着相当的作用。近壁区域划分见图4-1。 图4-1,边界层结构 第一节,壁面函数与近壁模型 近壁处理方法有两类:第一类是不求解层流底层和混合区,采用半经验公式(壁面函数)来求解层流底层与完全湍流之间的区域。采用壁面函数的方法可以避免改进模型就可以直接模拟壁面存在对湍流的影响。第二类是改进湍流模型,粘性影响的近壁区域,包括层流底层都可以求解。 对于多数高雷诺数流动问题,采用壁面函数的方法可以节约计算资源。这是因为在近壁区域,求解的变量变化梯度较大,改进模型的方法计算量比较大。由于可以减少计算量并具有一定的精度,壁面函数得到了比较多的应用。对于许多的工程实际流动问题,采用壁面函数处理近壁区域是很好的选择。 如果我们研究的问题是低雷诺数的流动问题,那么采用壁面函数方法处理近壁区域就不合适了,而且壁面函数处理的前提假设条件也不满足。这就需要一个合适的模型,可以一直求解到壁面。FLUENT提供了壁面函数和近壁模型两种方法,以便供用户根据自己的计算问题选择。

4.1.1壁面函数 FLUENT 提供的壁面函数包括:1,标准壁面函数;2,非平衡壁面函数两类。标准壁面函数是采用Launder and Spalding [L93]的近壁处理方法。该方法在很多工程实际流动中有较好的模拟效果。 4.1.1.1 标准壁面函数 根据平均速度壁面法则,有: **1 ln()U Ey k = 4-1 其中,1/41/2 * /p p w U C k U μτρ ≡ ,1/41/2 * p p C k y y μρμ≡,并且 k =0.42,是V on Karman 常数;E =9.81,是实验常数;p U 是P 点的流体平均速度;p k 是P 点的湍动能;p y 是P 点到壁面的距离;μ是流体的动力粘性系数。 通常,在*30~60y >区域,平均速度满足对数率分布。在FLUENT 程序中,这一条件改变为*11.225y >。 当网格出来*11.225y <的区域时候,FLUENT 中采用层流应力应变关系,即:**U y =。这里需要指出的是FLUENT 中采用针对平均速度和温度的壁面法则中,采用了*y ,而不是y +(/u y τρμ≡)。对于平衡湍流边界层流动问题,这两个量几乎相等。 根据雷诺相似,我们可以根据平均速度的对数分布,同样给出平均温度的类似分布。FLUENT 提供的平均温度壁面法则有两种:1,导热占据主要地位的热导子层的线性率分布;2,湍流影响超过导热影响的湍流区域的对数分布。 温度边界层中的热导子层厚度与动量边界层中的层流底层厚度通常都不相同,并且随流体介质种类变化而变化。例如,高普朗特数流体(油)的热导子层厚度比其粘性底层厚度小很多;对于低普朗特数的流体(液态金属)相反,热导子层厚度比粘性底层厚度大很多。 1/41/2 * ()w p p P T T c C k T q μρ-≡ '' 4-2 =()1/41/2 *2*1/41/222 1Pr Pr 21Pr ln()1Pr Pr Pr 2p p t p t p t c C k y U q Ey P k C k U U q μμρρ?+?''? ????++???? ??????+-??''?? ** **()()T T y y y y <> 4-3

湍流的研究进展

湍流的研究进展 XXX (XXX大学化工学院,青岛 266042) 摘要:本文对一百多年来湍流研究的进展作了简要回顾,并概述了湍流产生的原因及湍流对流体造成的影响,从不同的方向阐述了当今流体湍流的研究成果,展现了湍流研究的深入对于科学技术与社会发展产生的重要作用,展望了对于湍流研究的前景,并对湍流研究的发展提出了一些建议和设想。 关键词:湍流;湍流模式;流体湍流;湍流强度; The Turbulence of Research Progress XXXXX (Qingdao University of Science and Technology, Qingdao 266042) Abstract: Stupid hundred years Turbulence Research progress made brief review and an overview of the the turbulence causes and turbulent fluid caused today's fluid turbulence research, elaborated from a different direction, to show the turbulentdepth study of the important role of science and technology and social development, the future prospects for turbulence research, development and turbulence research has made some suggestions and ideas. Keywords: Turbulence; Turbulence models; Fluid turbulence; Turbulence intensity; 一、湍流研究的历史进程 人类很久前就已经观察到湍流运动了,但对它系统地进行研究则仅仅有一百多年的历史。经过一百多年的研究工作,人们的认识日益深化, 预测方法不断改进。回顾一下湍流研究取得进展的历程对于进一步揭示这一十分复杂流动现象是有益的。 涡团粘度概念首先是由波希尼斯克(Boussinesq)于1877年提出的,他的观点是湍流是一团杂乱无章的涡团。而现代湍流理论的创始人O.Reynolds则认为,湍流是由层流不稳定性发展起来的。这两位湍流研究的先驱者对湍流的认识有所不同。 本世纪二十年代湍流研究取得了巨大进展,有电子管补偿线路的热线风速计为湍流实验研究提供了有效的手段。 从四十年代到六十年代末湍流研究在理论和实验两方面都没有很大的突破。但是应用热线风速计测量各种湍流特性的资料大大充实了湍流的数据库。 六十年代末以后, 湍流研究又出现了一个新高潮,切变湍流中拟序结构的发现,复杂的湍流模式的建立和发展。湍流的直接数值模拟的尝试以及在方程中发现奇异吸引子或其它混沌现象的探索是近二十多年来湍流研究中的重大突破。

湍流研究简史-温景嵩

湍流研究简史-温景嵩 长春实验所发现的湍流不连续性及其对柯尔莫果洛夫理论基础的冲击具有十分重要的意义。(长春实验是指作者1972年9月在长春郊区采用类似热线风速仪的仪器测量大气湍流的温度脉动,也称温度脉动仪,然后通过频谱分析仪进行各谱段频谱分析。作者从中发现了湍流不连续性,也称间歇性。)因为湍流不仅是流体运动中的一个重大的世纪性的前沿课题,不仅它普遍存在于自然界,也普遍地存在于工程界,它是基础科学中一个重大的前沿分支---20世纪下半叶兴起的非线性科学的先驱和归宿。正由于以上两个原因,所以湍流问题的研究不仅吸引了众多的流体力学家,力学家的兴趣,而且也吸引了众多的数学家,物理学家,大气科学家,甚至包括了众多的工程技术界的专家学者的兴趣,大家都想在这一领域里一显身手。可以说湍流这一领域真正是“江山如此多娇,引无数英雄竞折腰”。自1883年英国曼彻斯特大学著名流体力学大师雷诺发表他的现代湍流开创性工作以来,一百二十多年里在湍流领域中已积累起浩如烟海的文献,发表了成百上千种的学说和理论,尽管如此,由于湍流这一课题固有的十分严重的困难,一百二十多年的众多科学家的奋斗结果,真正成功的理论并不多,算起来也就四个。 1. 普朗特的半经验混合长理论 第一个是1925年普朗特发表的半经验混合长理论,以及由此而导出的平板平均流速与所在高度的对数成正比的对数分布律。(冯. 卡尔曼1930,普朗特1933)这个对数分布律已由大量实验所证明。在工程上有很好的应用,可以用以计算平板表面所受的摩擦阻力,经过推广后,现在还可以用以计算飞船模型表面所受摩擦阻力。应该承认普朗特的半经验混合长理论解决了工程应用上的一大难题。后来前苏联学者莫宁(Monin)和奥布霍夫又把它成功地推广到近地面边界层大气风速的分布问题中去,为解决大气物理中的大气扩散等难题开辟了道路。然而普朗特的混合长理论并不是在工程应用中产生,也不是在大气中应用产生,也不是由实验带出来的结果。相反,它是在解决湍流这一学科发展中所面对的难题而产生的。它产生了以后,才有了工程的应用,才有了在大气中的应用,并且也才有了实验的证实。普朗特的半经验混合长理论是为解决雷诺方程的不闭合难题而创造出。1895年,也就是雷诺用实验证明湍流发生规律工作后的十二年,同样是由他研制成著名的雷诺方程。该方程从支配黏性流体运动的基本方程---纳维-斯托克斯方程出发,然后把瞬时流场分解为平均流场和湍流脉动速度流场的和,把这个和式代入到纳维-斯托克斯方程再取平均就形成了雷诺方程,这是一个支配湍流场中平均流场变化的方程,不幸方程不闭合。因为除了待求的平均流场外,又多了一个未知数,即同一点上湍流脉动速度的两个分量相关矩,它具有应力的量纲,又叫雷诺应力。它表征了湍流脉动场对平均场的影响,相关矩肯定不为0 ,即雷诺应力不是0。否则有湍流发生后的平均流场分布规律就应和没有湍流发生时的层流流场规律相同。而实验已证实,两者确实不同,这就证实湍流场的雷诺应力对平均场确有重要影响。可惜这是未知的。于是一个雷诺方程无法同时解出平均场和雷诺应力两个未知数,形成湍流研究中著名的不闭合难题,这个难题是由纳维-斯托克斯的非线性,以及湍流特有的随机性,在对方程求取平均值过程中必然产生。所以是湍流研究中固有的一个难点。用同样的雷诺方法,原则上可以求出湍流脉动速度两个分量相关矩的方程,这样方程就多了一个,此时和原来的雷诺方程一起现在有了两个方程,两个未知数,似乎可以闭合,其实不然。从纳维-斯托克斯方程的非线性特点,可以断定在建立两个分量的二阶相关矩方程时,必然又会增加一个新的未知的三阶相关矩,方程仍然不闭合,依此类推,若建立三阶相关矩方程,则同样还会多出一个未知的四阶相关矩,可以断言,沿着这条路线下去,未知数永远要比方程多一个,方程不可能闭合。这样下去,湍流问题就无法严格在数学上求解。雷诺方程建立后又过了三十年,即1925年由普朗特用混合长理论解决了这个难题。他的解决办法就是用物理模型方法来切断雷诺方程在数学上的不封闭链条,在雷诺方程那里就打住,引入混合长的物理模型,使雷诺

中国湍流研究的发展史_中国科学家早期湍流研究的回顾

中国湍流研究的发展史 I 中国科学家早期湍流研究的回顾 黄永念 北京大学力学与工程科学系,湍流与复杂系统国家重点实验室,北京,100871 摘要总结了二十世纪三十年代到六十年代中国老一辈科学家(包括物理学家,力学家)周培源、王竹溪、张国藩、林家翘、谢毓章、张守廉、黄授书、胡宁、柏实义、陈善模、庄逢甘、陆祖荫、李政道、蔡树棠、是勋刚、李松年、谈镐生、包亦和等诸位先生的湍流研究工作。介绍他们对流体力学中最为困难的湍流问题所作出的努力和贡献。 关键词湍流统计理论,能量衰变规律,均匀各向同性湍流,剪切湍流。 引言 湍流一直被认为是物理学中最难而又久未解决的基础理论研究的一个课题。从1883年Reynolds圆管湍流实验研究算起已经跨越了两个世纪,湍流问题仍未得到解决。在跨入二十一世纪时,很多从事湍流研究工作的科学家都在思考这样的问题:二十世纪的湍流研究留给我们哪些宝贵财富?二十一世纪又应该如何面对这个老大难问题?Yaglom在2000年法国举行的一次湍流讲习班上回顾了二十世纪的湍流理论发展过程[1],指出了其中两个最重要的成就:一个是Kolmogorov的局部均匀各向同性湍流理论,另一个是von Karman的湍流平均速度的对数分布律。同时又一次向世人介绍著名科学家Lamb在临终前对解决湍流问题的悲观看法。由于中国与世界各国在文字和语言上的差异和长期缺乏国际间的交流,历次湍流研究工作的总结和回顾中,人们往往忽略了中国科学家的作用。只有周培源教授在1995年流体力学年鉴上发表了“中国湍流研究50年”才打破了这种隔阂[2]。但是这篇文章也只局限于周培源教授率领的北京大学研究组所做的系列研究工作。实际上有很多中国科学家在上一世纪中做了非常出色的工作。本文仅就半个世纪前的三十年代到六十年代他们的湍流研究工作做一个简单的介绍,目的是要引起大家关注中国科学家的湍流研究和对湍流研究所做的贡献。 中国科学家的湍流研究工作可以分成两个方面,一是在国内极其困难的条件下坚持开展的研究工作,这方面的工作国际上鲜为人知。另一方面是在国外开展的研究工作,这部分工作国内也不很熟悉。因此,本文将把他们的不懈努力介绍给大家。 胡非在1995年发表的专著《湍流,间隙性与大气边界层》中曾专门介绍了中国学者的湍流研究工作[3],但他的介绍还不够全面,特别是缺少对早期工作的报道。本文可以弥补其中的不足。 1 三十年代的研究工作 在我国最早发表湍流论文的是当时在清华大学的王竹溪先生。他在周培源先生的指导下

湍流与层流_湍流研究概述

第一篇 大气的组成与物理特性 第一章 第二章 第三章 第四章 第五章 大气的气体成份 大气中的粒子群 大气的运动、能量与构造 大气的光学特性 大气的电学特性
1

第二篇 大气湍流
粘性流体的两种形态: 层流和湍流。 层流是流体运动中较简单的状态, 普遍的却是湍流。
2

湍流研究的意义
湍流的研究与国防建设和国民经济中 的航空、船运、环境保护、气象、化工、 冶金、水利、医学等学科密切相关,如果 能掌握它的运动规律,对它进行合理的应 用和有效的控制,那么对基础研究与实际 应用将有重大的意义。
3

湍流研究的成果
人们对湍流结构、湍流边界层、湍流 剪切流、湍流的传热传质、湍流扩散、湍 流统计模型、大气湍流、晴空湍流、等离 子湍流、湍流测量等问题进行了广泛的研 究,并取得了丰硕的成果。
4

本节的内容
湍流的一般定义和描述; 湍流与层流的区别; 湍流理论发展的历史; 湍流理论简介; 湍流的特点; 大气湍流的复杂性; 湍流研究技术的发展。
5

湍流的一般定义和描述
1. 湍流是随机的(Reynolds,Taylor,Von Karman ,Hinze等),又具有拟序结 构。 2. 流体的湍流运动是由各种大小和涡量 不同的涡旋叠加而成的,其中最大涡 尺度与流动环境密切相关,最小涡尺 度则由粘性确定;流体在运动过程中, 涡旋不断破碎、合并,流体质点轨迹 不断变化。
6

关于湍流理论研究进展精品资料

关于湍流理论研究进展 摘要本文对近年来湍流理论在某些方面的研究进展作了概要介绍,对具有代表性的理论假设的思想方法,进行了扼要阐述,指出了相应的实用价值和局限性。 关键词湍流湍流统计理论混沌理论湍流拟序结构湍流剪切流动 1 无处不在的湍流现象 湍流是自然界中流体的一种最普遍的运动现象,它广泛的存在于我们生活周围。在大风吹过地面障碍物的旁边,在湍急的河水流过桥墩的后面,在烟囱中冒出的浓烟随风渐渐扩散等地方,都能观察到湍流运动现象。简单地说,湍流运动就是流体的一种看起来很不规则的运动。由于湍流现象广泛存在于自然界和工程技术的各个领域,因此湍流基础理论研究取得的进展就可能为经济建设和国防建设的广泛领域带来巨大的效益。例如,提高各种运输工具的速度以大量节约能源,提高各种流体机械的效益;改善大气和水体的环境质量,降低流体动力噪声,防止流体相互作用引发的结构振动乃至破坏;加强反应器内部物质的热交换与化学反应的速度等等。 然而像湍流这样,虽经包括许多著名科学家在内长达一个世纪多的顽强努力,正确反映客观规律的系统的湍流理论至今还没有建立,在整个科学研究史上也是不多见的。因此,可以说湍流是力学中没有解决的最困难的难题之一。因此,世界上许多国家一直坚持把湍流研究列为需要最优先发展的若干重大基础研究课题之一。 2 湍流理论的发展历史 湍流理论从它的思路来说大体可分为两类[1]。一类是先把流体动力学方程组平均以后,然后再设法使方程组封闭,求解后再和实验结果比较,看封闭办法是否正确。湍流中绝大部分理论是属于这一类型。另一类是先求解,取特殊模型,再引进平均,得到要求的物理量,和相应的实验结果进行比较。 2.1 Reynolds方程和混合长度理论 十九世纪70年代是Maxwell-Boltzmann分子运动理论取得辉煌成果的时代。它成功地解释了气体状态方程、气体粘性、气体热传导和气体扩散等

湍流理论若干问题研究进展

第15卷第4期水利水电科技进展1995年8月 湍流理论若干问题研究进展 刘兆存 金忠青 (河海大学 南京 210098) 摘要 本文对近年来湍流理论在某些方面的研究进展作了概要介绍,对拟序结构发现后人们对湍流内部结构的新认识和近年来发展很快的从微分方程分析角度出发对湍流机理新的探索进行了评价,说明引入混沌后在时、空演化方面对湍流机理的模拟,最后阐述了流动稳定性和层流向湍流的转捩。 关键词 湍流 N-S方程 流动结构 流动机理 封闭性 近年来,在围绕湍流结构和统计两条主线的研究工作中出现了新观点和新趋势,虽然从历史的观点来看有些可能是错的——在科学容忍的范围内,但在现阶段却是研究的主流。 1 简要回顾及发展 1.1 半经验理论和模式理论 湍流的控制方程是N-S方程,但和层流相比,方程不封闭。为满足工程需要,发展了一系列的以普朗特混合长理论为代表的湍流半经验理论或早期模式理论。这种理论虽然对于增进对湍流机理的了解没有提供更多的贡献,但对解决工程实际问题却起了重大的作用[1]。半经验理论是一种唯像理论,并不涉及湍流内部机理。以速度分布公式为例,半经验理论的速度分布公式大致有对数型和指数型。对数型速度分布得到的假定是充分发展的剪切湍流中主流区(不含边界层的)的流速梯度和分子粘性无关,指数型(或渐近指数型)则假定分子粘性不能忽略[2],两种类型的流速分布公式在工程实践中都获得了非常广泛的应用。半经验理论的一个发展方向是吸收统计理论的成果,用统计理论的精细成果丰富半经验理论不足并保留便于应用的优点,如文[3]所作的工作。 近代的模式理论在封闭湍流基本方程组时特别吸收了统计理论的成果,如二方程模型、应力通量代数模型、应力通量方程模型等。关于这方面的详细论述,将另文给出。 1.2 统计理论 湍流的统计理论的目标则是从最基本的物理守恒定律——N-S方程和连续性方程出发,探讨湍流的机理。理查逊-柯尔莫哥洛夫湍流图像部分被实验所证实。统计理论中湍流的能量传递关系被更符合实际的U. Fr isch等所提出的B-模型所代替。湍流统计理论历时半个多世纪的发展,经泰勒、陶森德等人的努力,取得丰硕的成果,但仍不能绕过封闭性的困难,所得成果都还是很不完善的。湍流统计理论的重要性目前已有所下降[1]。我国周培源等提出了均匀各向同性湍流的准相似性条件以及相应均匀各向同性湍流的涡旋结构统计理论并得到实验的验证[4],进一步将在均匀各向同性湍流中得到的准相似性条件推广到一般的剪切湍流中,然后对关联方程的耗散项作出假定,利用逐级近似方法发展了湍流的统计理论[5],所得结果部分经实验证实。文[6]采用逐级迭代法对湍流平均运动方程和脉动速度关联方程 · 12·

近壁面函数的简单理解

一个成功的湍流计算离不开好的网格。在许多的湍流中,空间的有效粘性系数不同,是平均动量和其它标量输运的主要决定因素。因此,如果需要有足够的精度,这就需要保证湍流量要比较精确求解。由于湍流与平均流动有较强的相互作用,因此求解湍流问题比求解层流时候更依赖网格。对于近壁网格而言,不同的近壁处理对网格要求也不同。下面对常见的几种近壁处理的网格要求做个说明。采用壁面函数时候的近壁网格:第一网格到壁面距离要在对数区内。对数区的y+ >30~60。FLUENT在y+ <时候采用层流(线性)准则,因此网格不必要太密,因为壁面函数在粘性底层更本不起作用。对数区与完全湍流的交界点随压力梯度和雷诺数变化。如果雷诺数增加,该点远离壁面。但在边界层里,必须有几个网格点。壁面函数处理时网格划分采用双层模型时近壁网格要求当采用双层模型时,网格衡量参数是y+ ,并非y* 。最理想的网格划分是需要第一网格在y+ =1位置。如果稍微大点,比如=4~5,只要位于粘性底层内,都是可以接收的。理想的网格划分需要在粘性影响的区域内(Rey<200 )至少有十个网格,以便可以计算粘性区域内的平均速度和湍流量。采用双层区模型时网格划分采用Spalart-Allmaras 模型时的近壁网格要求该模型属于低雷诺数模型。这就要求网格能满足求解粘性影响区域内的流动,引入了阻尼函数,用以削弱粘性底层的湍流粘性影响。因此,理想的近壁网格要求和采用双层模型时候的网格要求一致。采用大涡模拟的近壁网格要求对于大涡模拟,壁面条件采用了壁面法则,因此对近壁网格划分没有太多限制。但是,如果要得到比较好的结果,最好网格要细,最近网格距离壁面在 y+=1的量级上。 for Hexa mesh, ==>Y+是第一层高度一半和 viscous length scale 的比值 for Tetra mesh==>Y+是第一层高度1/3和 viscous length scale 的比值 y+就是Yplus,它跟你在湍流模型里采用的近壁面函数选取有关,若Yplus为个位数,选增强型壁面函数,若在两位数以上,选标准或非平衡的壁面函数。 y+的意思是底层网格必须划分在对数率成立的区域内。 一般应使y+的值为15~300,但是y+是模拟完成后才知道的。 而且同一个模型不同地方不同流速y+不一样,所以不是很精确。如果模拟传热应注意y+对结果的影响。

湍流理论发展概述

湍流理论发展概述 一、湍流模型的研究背景 自然环境和工程装置中的流动常常是湍流流动,模拟任何实际过程首先遇到的就是湍流问题,而湍流问题本身又是流体力学理论上的难题。对于某些简单的均匀时均流场,如果湍流脉动是各向均匀及各向同性的,可以用经典的统计理论来分析,但实际上的湍流往往是不均匀的,这就给理论分析带来了极大地困难。这也就引发了对湍流过程进行模拟的想法。 对湍流最根本的模拟方法是在湍流尺度的网格尺寸内求解瞬态的三维N-S 方程的全模拟方法,此时无需引进任何模型。然而由于计算方法及计算机运算水平的限制,该种方法不易实现。另一种要求稍低的方法是亚网格尺寸度模拟即大涡模拟(LES),也是由N-S 方程出发,其网格尺寸比湍流尺度大,可以模拟湍流发展过程的一些细节,但由于计算量仍然很大,只能模拟一些简单的情况,直接应用于实际的工程问题也存在很多问题[1]。目前数值模拟主要有三种方法:1. 平均N-S方程的求解,2.大涡模拟(LES),3.直接数值模拟(DNS),而模拟的前提是建立合适的湍流模型。 所谓的湍流模型,就是以雷诺平均运动方程与脉动运动方程为基础,依靠理论与经验的结合,引进一系列模型假设,而建立起的一组描写湍流平均量的封闭方程组。目前常用的湍流模型可根据所采用的微分方程数进行分类为:零方程模型、一方程模型、两方程模型、四方程模型、七方程模型等。对于简单流动而言,一般随着方程数的增多,精度也越高,计算量也越大、收敛性也越差。但是,对于复杂的湍流运动,则不一定。湍流模型可根据微分方程的个数分为零方程模型、一方程模型、二方程模型和多方程模型。这里所说的微分方程是指除了时均N-S 方程外,还要增加其他方程才能是方程封闭,增加多少个方程,则该模型就被成为多少个模型。

湍流的研究进展论文

湍流的研究进展 丁立新 (青岛科技大学) 摘要本文重点就湍流的理论研究进展作一阐述,从湍流的相干结构、表征及发展由来,到上世纪末湍流研究进展的雷诺方程,本世纪湍流的统计理论和半经验理论发展,湍流的模式理论,湍流的高级数值模拟分别论述,并为主要的工程应用做简要的介绍。 关键词湍流理论研究工程应用 Research process of turbulence Dinglixin Qingdao University of Science & technology Abstract This article focuses on the turbulence of research process as elaborated. From coherent structure of turbulence, characterization and development of turbulence to Reynolds equation about research process of turbulence on the end of the century, the development of semi-empirical theory and statistical theory of turbulence of this century, mode theory of turbulence, advanced numerical simulation of turbulence. Finally, brief description of turbulence industrial applications is suggested. Keywords Turbulence, Theoretical research of turbulence, Engineering applications 湍流是自然界和工程中最常出现的流动形态,湍流的出现将使动量、质量、能量的输送速率极大地加快,一方面造成能量消耗加快,污染物加快扩散等严重消极

珠峰北坡地区近地层大气湍流与地气能量交换特征

第21卷第12期2006年12月 地球科学进展 A DVANCE S I N E AR TH S C I ENC E V o l.21 N o.12 D e c.,2006 文章编号:1001-8166(2006)12-1293-11 珠峰北坡地区近地层大气湍流与地气能量交换特征 仲雷1,2,3,马耀明1,4,苏中波5,刘 新1,李茂善3,4,马伟强3,4,王永杰1,3 (1.中国科学院青藏高原研究所,北京 100085;2.中国气象局成都高原气象研究所高原气象开放实验室,四川 成都 610071;3.中国科学院研究生院,北京 100039;4.中国科学院寒区旱区环境与工程研究所,甘肃 兰州 730000;5. I n t er na ti o nal I n s tit u t e f o r G e o-I n f o r m a t i o n S c i e n c e a nd E a r t h O b serv a ti o n,E n sc h e de7500,t he N e t h er l ands) 摘 要:利用珠峰北坡曲宗地区连续一年的大气观测资料(2005年4月至2006年3月),分析了珠峰北坡地区近地层大气湍流宏观统计特征和西南季风爆发前后地气能量交换特征。研究表明在珠峰北坡地区M o n i n-O bukh o v相似定律同样适用。拟合得到了珠峰北坡曲宗地区近地层无因次风速分量方差以及温度和湿度归一化标准差和静力学稳定度的函数关系。研究得出曲宗地区能量平衡各分量(净辐射通量、感热通量、潜热通量和土壤热通量)以及地面加热场具有明显的季节变化和日变化规律。尤其是在西南季风的影响下,曲宗地区感热通量和潜热通量在季风爆发前后具有明显相反的变化趋势。其它特征参数(波文比和地表反射率)在西南季风爆发前后的变化规律也十分明显。 关 键 词:近地层;大气湍流;能量通量;珠峰;曲宗 中图分类号:P425.2 文献标识码:A 1 引 言 青藏高原地处我国西部,约占国土面积的四分之一,平均海拔在4000m以上。高原热力、动力作用以及地—气间的物质能量交换过程对我国、亚洲乃至全球的气候变化均有重大影响。青藏高原对大气的动力和热力作用主要是通过下垫面与大气的相互作用,并以湍流方式进行物质和能量交换而实现的。开展高原上地气系统物理过程的观测,分析确定热量、水汽等湍流参数的变化特征,将有助于改进全球气候模式和区域天气、气候模式在该地区的参数化方案。因此,对青藏高原的研究愈来愈受到中外学者的关注。20世纪50年代以来科学家进行了多次关于青藏高原的气象科学试验,如第一、第二次青藏高原大气科学实验(Q X P M E X,1979年5~8月,T I P E X,1998年5~8月),全球能量水分循环亚洲季风之青藏高原试验研究(G E W E X/G A M E-T i b e t,1996—2000年),“全球协调加强观测计划之亚澳季风青藏高原试验”(C EO P/C A M P-T i b e t,2001—2005年),积累了大量的宝贵资料,并且取得了丰硕的科研成果[1~10]。但是青藏高原特殊的地理条件和恶劣的气候环境给野外观测试验造成极大困难,使得很多试验只能在现有城市附近和交通相对便利的地区展开,不足以了解像珠穆朗玛峰(以下简称珠峰)这样的大地形对大气环流的影响。喜马拉雅山脉山体是北半球地表与对流层大气物质交换的重要通  收稿日期:2006-10-11;修回日期:2006-10-24. *基金项目:科技部社会公益研究专项“珠穆朗玛峰地区对全球变化的响应”(编号:2005D I A3J106);中国气象局成都高原气象研究所高原气象开放基金课题“青藏高原地表特征参数卫星遥感反演研究”(编号:L P M2006011);中国科学院知识创新工程重要方向项目“喜马拉雅山北坡地区地面大气与对流层大气交换研究”(编号:K ZCX3-S W-231);国家自然科学基金项目“西藏高原能量水循环降雨共同观测研究”(编号:40520140126)资助.  作者简介:仲雷(1979-),男,安徽蚌埠人,博士生,主要从事大气边界层观测和卫星遥感应用研究.E-m a i l:z hongl@i t pca s. a c. c n *通讯作者:马耀明(1964-),男,山西夏县人,研究员,博导,主要从事陆面过程和遥感应用研究.E-m a i l:y m m a @i t pcas. a c. c n

第三章_湍流模型

第三章 湍流模型 第一节 前言 湍流流动模型很多,但大致可以归纳为以下三类: 第一类是湍流输运系数模型,是Boussinesq 于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即: 2 1 21 x u u u t ??=''-μρ 3-1 推广到三维问题,若用笛卡儿张量表示,即有: ij i j j i t j i k x u x u u u δρμρ32 -??? ? ????+ ??=''- 3-2 模型的任务就是给出计算湍流粘性系数t μ的方法。根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。 第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。 第三类是大涡模拟。前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes 方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。 实际求解中,选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。 FLUENT 提供的湍流模型包括:单方程(Spalart-Allmaras )模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。 湍流模型种类示意图 Direct Numerical Simulation 包含更多 物理机理 每次迭代 计算量增加 提的模型选 RANS-based models

浅谈湍流的认识与发展

浅谈湍流的认识与发展 摘要:本文结合流体力学课程的学习以及对湍流相关书籍的阅读,阐述个人对湍流运动的发展、特点、性质的理解。湍流作为“经典物理学最后的疑团”,人们不断地进行探索,建立湍流模型对其进行研究理论分析。近年来,对于湍流这一不规则运动,人们提出了并且倾向于应用混沌理论进行分析,并取得了一些成果。对湍流的认识在不断深入。 关键字:湍流概念湍流性质湍流强度模型建立混沌理论 在流体力学的学习过程中, 湍流一度被称为“经典物理学最后的疑团”,我对湍流这一流体的状态极其相关的力学性质进行了更深入的了解与学习,结合课堂上老师的讲解以及课后对相关参考文献的阅读理解,在此我想浅谈一下这一阶段我对湍流的学习与认识。 从湍流的定义出发,初识湍流,湍流是流体的一种流动状态。对于流体,大家都知道,当流速很小时,流体分层流动,互不混合,称为层流,也称为稳流或片流;逐渐增加流速,流体的流线开始出现波浪状的摆动,摆动的频率及振幅随流速的增加而增加,此种流况称为过渡流;当流速增加到很大时,流线不再清楚可辨,流场中有许多小漩涡,层流被破坏。这时的流体作不规则运动,有垂直于流管轴线方向的分速度产生,这种运动称为湍流。流体作湍流时,阻力大流量小,能量耗损增加。能量耗损E与速度的关系为△ E= kv2(k是比例系数,它与管道的形状、大小以及管道的材料有关。v是平均流速)。所有流体都存在湍流现象。 我们可以用雷诺数的范围量化湍流。在直径为d的直管中,若流体的平均流速为v,由流体运动粘度v组成的雷诺数有一个临界值(大约为2300~2800),若Re小于该范围则流动是层流,在这种情况下,一旦发生小的随机扰动,随着时间的增长这扰动会逐渐衰减下去;若Re大于该范围,层流就不可能存在了,一旦有小扰动,扰动会增长而转变成湍流。雷诺在1883年用玻璃管做试验,区别出发生层流或湍流的条件。把试验的流体染色,可以看到染上颜色的质点在层流时都走直线。当雷诺数超过临界值时,可以看到质点有随机性的混合,在对时间和空间来说都有脉动时,这便是湍流。不用统计、概率论的方法引进某种量的平均值就难于描述这一流动。除直管中湍流外还有多种多样各具特点的湍流,虽经大量实验和理论研究,但至今对湍流尚未建立起一套统一而完整的理论。在流

智能材料研究进展(毕业论文doc)

智能材料研究进展 摘要 智能材料是一门多门类、多学科交叉的科学,与物理学、材料力学、电子学、化学、仿生学、生命科学、控制理论、人工智能、信息技术、生物技术、计算机技术、材料合成与加工等诸多的前沿科学及高新技术戚戚相关、紧紧相连。因此,它一旦有所突破,便会导致众多学科的理论创新和许多领域的技术变革,大大地推动国家科学技术的进步和综合实力的提高。智能材料具有十分重要的现实用途和极为广阔的应用前景。从高精尖的宇宙探索,到普通人的日常生活,智能材料都起着重要的作用。 未来社会发展的趋势是智能化。智能化的首要问题是大力发展智能材料,智能材料的研究是材料科学研究的重要方向。智能材料的本质特征是材料具有仿生功能,即材料能根据感受到的信息而自动判断、控制和调整以适应外界条件变化。 本文介绍了智能材料的概念、定义及智能材料的特征,阐述和评价了智能材料——形状记忆合金、电流变体材料、光致变色材料、电致变色材料、形状记忆复合材料和智能型药物释放体系等的种类、组成、特点、用途、研究现状与市场前景。重点论述了压电陶瓷材料的制造工艺、特点、性质、研究现状及市场前景等。论述了发展智能材料的战略意义,展望了它的发展前景。 关键词:智能材料,研究,应用,发展

DEVELOPMENT PROGRESS OF SMART MATERIALS ABSTRACT Smart materials is more than one categories, interdisciplinary science, and physics, mechanics, electronics, chemistry, bionics, life sciences, control theory, artificial intelligence, information technology, biotechnology, computer technology, materials synthesis and processing and many other leading edge science and very much related to high-tech, tightly linked. Therefore, once it has been a breakthrough, it will lead to many disciplines in many areas of theoretical innovation and technological change; greatly promote national scientific and technological progress and the improvement of overall strength. Smart materials is of great practical use and very broad application prospects. Explore the universe from the sophisticated to the daily lives of ordinary people, smart materials play an important role. The trends of coming society are intellectualization. The essential issue of intellectualization is to develop intelligent materials vigorously. The study of intelligent materials is a crucial direction of material science.The main characteristic of intelligent materials is bionics functions. That is, it can judge, control and adjust it automatically to adapt the change of the external environment according to accepting information. In this paper, the concept of smart materials, definitions, describes the characteristics of smart materials, intelligent materials described and evaluated - shape memory alloys, electrorheological materials, photochromic materials, electrochromic materials, shape memory composites and smart based drug delivery system, and the type, composition, characteristics, uses, current situation and market prospects. Focuses on the manufacturing process of piezoelectric ceramic materials, characteristics, nature, current situation and

相关文档