存在热浮力的剪切湍流中颗粒近壁运动的大涡模拟
LES,DNS,RANS三种模拟模型计算量比较及其原因
LES,DNS,RANS模型计算量比较 摘要：湍流流动是一种非常复杂的流动，数值模拟是研究湍流的主要手段，现有的湍流数值模拟的方法有三种：直接数值模拟（Direct Numerical Simulation: DNS），Reynolds平均方法（Reynolds Average Navier-Stokes: RANS）和大涡模拟(Large Eddy Simulation: LES)。直接数值模拟目前只限于较小Re数的湍流，其结果可以用来探索湍流的一些基本物理机理。RANS方程通过对Navier-Stokes方程进行系综平均得到描述湍流平均量的方程；LES方法通过对Navier-Stokes方程进行低通滤波得到描述湍流大尺度运动的方程，RANS和LES方法的计算量远小于DNS，目前的计算能力均可实现。 关键词：湍流；直接数值模拟；大涡模拟；雷诺平均模型 1 引言 湍流是空间上不规则和时间上无秩序的一种非线性的流体运动，这种运动表现出非常复杂的流动状态，是流体力学中有名的难题，其 性。传统计算复杂性主要表现在湍流流动的随机性、有旋性、统计[]1 流体力学中描述湍流的基础是Navier-Stokes（N-S）方程，根据N-S 方程中对湍流处理尺度的不同，湍流数值模拟方法主要分为三种：直接数值模拟（DNS）、雷诺平均方法（RANS）和大涡模拟（LES）。直接数值模拟可以获得湍流场的精确信息，是研究湍流机理的有效手段，但现有的计算资源往往难以满足对高雷诺数流动模拟的需要，从而限制了它的应用范围。雷诺平均方法可以计算高雷诺数的复杂流动，但给出的是平均运动结果，不能反映流场紊动的细节信息。大涡模拟基于湍动能传输机制，直接计算大尺度涡的运动，小尺度涡运动对大尺度涡的影响则通过建立模型体现出来，既可以得到较雷诺平均方法更多的诸如大尺度涡结构和性质等的动态信息，又比直接数值模拟节省计算量，从而得到了越来越广泛的发展和应用。
管内湍流的数值模拟
管内湍流的数值模拟 摘要：当Reynolds数大于临界值时，平滑流动会出现一系列复杂的变化，最终会导致流动特征的本质变化，流动呈无序的混乱状态，这种状态称为湍流。计算流体力学是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。本文以湍流管流模型为例，借助Fluent软件进行空气动力学分析，对该管内湍流流动进行模拟。 关键词：计算流体力学；Fluent；管内湍流；数值模拟 1 引言 流体试验表明，当Reynolds数大于临界值时，平滑流动会出现一系列复杂的变化，最终会导致流动特征的本质变化，流动呈无序的混乱状态。这时，即使是边界条件保持不变，流动也是不稳定的，速度等流动特性都随机变化，这种状态称为湍流。 随着高速电子计算机的出现，数值模拟越来越多地应用于流场的模拟。计算流体力学(Computational Fluid Dynamics ,简称为CFD)就是其中一种有效的研究流体动力学的数值模拟方法，它是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析；是基于计算机技术的一种数值计算工具, 用于求解流体的流动和传热问题。它能够描述几何体边界的复杂的流动现象，能够在设计的初期快速地评价设计并做出修改；在设计的中期，用来研究设计变化对流动的影响，减少未预料到的负面影响；设计完成后，CFD提供各种数据和图像，证实设计目的。CFD大大减少了费用、时间以及新设计带来的风险。近年来，CFD越来越多地应用于翼型设计和流场的分析中，成为一种重要的设计和计算方法。 Fluent软件是用于模拟和分析在复杂几何区域内的流体流动与热交换问题的专用CFD软件。它用于计算计算流体流动和传热问题的软件，其应用的范围有一般流体的流场、自由表面的问题、紊流、非牛顿流流场、化学反应等。Fluent提供了灵活的网格特性，用户可以方便的使用结构网格和非结构网格对各种复杂区域进行网格划分。本文以湍流管流模型为例，借助Fluent软件进行空气动力学分析，对该管内湍流流动进行模拟，并分析了模型内的中心速度分布、表面摩擦系数和流速剖面。 2 数学及物理模型的建立 2.1 数学模型
大涡模拟
4.6.3大涡模拟LSE 大涡模拟LES 基本思想是：湍流运动是湍流运动是由许多大小不同尺度的涡旋组成，大尺度的涡旋对平均流动影响比较大，各种变量的湍流扩散、热量、质量、动量和能量的交换以及雷诺应力的产生都是通过大尺度涡旋来实现的，而小尺度涡旋主要对耗散起作用，通过耗散脉动来影响各种变量。不同的流场形状和边界条件对大涡旋有较大影响，使它具有明显的各向不均匀性。而小涡旋近似于各向同性，受边界条件的影响小，有较大的共性，因而建立通用的模型比较容易。据此，把湍流中大涡旋(大尺度量)和小涡旋(小尺度量)分开处理，大涡旋通过N-S 方程直接求解，小涡旋通过亚格子尺度模型，建立与大涡旋的关系对其进行模拟，而大小涡旋是通过滤波函数来区分开的。对于大涡旋，LES 方法得到的是其真实结构状态，而对小涡旋虽然采用了亚格子模型，但由于小涡旋具有各向同性的特点，在采用适当的亚格子模式的情况下，LES 结果的准确度很高。 大涡模拟LES 有四个一般的步骤： ①定义一个过滤操作，使速度分解u(x,t)为过滤后的成分(),u x t 和亚网格尺度成分u ’(x,t)，这里要特别指出：过滤操作和Reynolds 分解是两个不同的概念，亚网格尺度SGS 成分u ’(x,t)与Reynolds 分解后的速度脉动值是两个不同的量。过滤后的三维的时间相关的成分()t x u ,表示大尺度的涡旋运动； ②由N-S 方程推导过滤后的速度场进化方程，该方程为一个标准形式，其中包含SGS 应力张量； ③封闭亚网格尺度SGS 应力张量，可采用最简单的涡黏性模型； ④数值求解模化方程，从而获得大尺度流动结构物理量。 （1）过滤操作 LES 方法和一般模式理论不同之处在于对N-S 方程第一步的处理过程不一样。一般模式理论方法是对变量取平均值，LES 方法是通过滤波操作，将变量分成大尺度量和小尺度量。对任一流动变量(),u x t 划分为大尺度量(,)u x t 和小尺度量(),u x t '（亚格尺度）： (,)(,)(,)u x t u x t u x t '=+ 其中大尺度量是通过滤波获得:，过滤操作定义为： ()?-=dr t r x u x r G t x u ),(),(, （4.78） 式中积分遍及整个流动区域，(,)G r x 是空间滤波函数，它决定于小尺度运动的尺寸和结构。 滤波器G 要满足正规化条件 ?=1),(dr x r G （4.79） 亚网格尺度SGS 成分定义为 ),(),(),('t x u t x u t x u -= （4.80） 与Reynolds 分解不同的是，),(t x u 为一个随机的场分布，且 0),('≠t x u
湍流的数值模拟方法进展
3 大涡模拟（LES ） 湍流大涡数值模拟（LES ）是有别于直接数值模拟和雷诺平均模式的一种数值模拟手段。利用次网格尺度模型模拟小尺度湍流运动对大尺度湍流运动的影响即直接数值模拟大尺度湍流运动， 将N-S 方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。 3.1 基本思想 很多尺度不同的旋涡一起组成了湍流运动平均流动主要取决于大漩涡的流动，大尺度运动则受到小旋涡的影响。流动中的大涡实现了动量、能量质量、热量的交换，耗散主要是由于小涡作用的。大旋涡中受到流场形状、阻碍物的影响，，使大漩涡的各向异性更加明显。然而小漩涡之间各项同性，相互没有太大的区别，所以建立统一的模型比较容易一些。综上所述，大涡模拟将湍流瞬时运动量通过滤波将运动分成小尺度和大尺度。大尺度的运动受到小尺度的运动的影响可以通过应力项（类似于雷诺应力项）来表示，即为亚格子雷诺应力，以建立这种模型的方法来模拟。而大尺度则是求解运动微分方程而计算出来的，也就是说大涡模拟，要先过滤掉小尺度的脉动，然后再推出小尺度的运动封闭方程以及大尺度的运动控制方程。 3.2 滤波函数 正如上面提到，大涡模拟要先将流动变量分解成小尺度量和大尺度量，我们把这个作用叫做滤波。滤波运算就是在一区域内按照一定的条件对函数进行加权平均，作用是将高波数滤掉，使低波数保留，滤波函数的特征尺度决定了截断波数的最大波长，下面三种滤波函数是最为常用的主要有以下三种：盒式、富氏截断以及高斯滤波函数。 不可压常粘性系数的湍流运动控制方程为N-S 方程: j ij i j j i i x S x P x u u t u ???+??-=??+??)2(1γρ 式中：S 拉伸率张量，表达式为：2/)//(i j j i ij x u x u S ??+??=；γ分子粘性系数；ρ流体密度。设将变量i u 分解为方程(11)中i u 和次网格变量(模化变量)'i u ，
大涡模拟的FLUENT算例2D
Tutorial:Modeling Aeroacoustics for a Helmholtz Resonator Using the Direct Method(CAA) Introduction The purpose of this tutorial is to provide guidelines and recommendations for the basic setup and solution procedure for a typical aeroacoustic application using computational aeroacoustic(CAA)method. In this tutorial you will learn how to: ?Model a Helmholtz resonator. ?Use the transient k-epsilon model and the large eddy simulation(LES)model for aeroacoustic application. ?Set up,run,and perform postprocessing in FLUENT. Prerequisites This tutorial assumes that you are familiar with the user interface,basic setup and solution procedures in FLUENT.This tutorial does not cover mechanics of using acoustics model,but focuses on setting up the problem for Helmholtz-Resonator and solving it.It also assumes that you have basic understanding of aeroacoustic physics. If you have not used FLUENT before,it would be helpful to?rst review FLUENT6.3User’s Guide and FLUENT6.3Tutorial Guide. Problem Description A Helmholtz resonator consists of a cavity in a rigid structure that communicates through a narrow neck or slit to the outside air.The frequency of resonance is determined by the mass of air in the neck resonating in conjunction with the compliance of the air in the cavity. The physics behind the Helmholtz resonator is similar to wind noise applications like sun roof bu?eting. We assume that out of the two cavities that are present,smaller one is the resonator.The motion of the?uid takes place because of the inlet velocity of27.78m/s(100km/h).The ?ow separates into a highly unsteady motion from the opening to the small cavity.This unsteady motion leads to a pressure?uctuations.Two monitor points(Point-1and Point-2) act as microphone points to record the generated sound.The acoustic signal is calculated within FLUENT.The?ow exits the domain through the pressure outlet.
湍流大涡数值模拟进展
第22卷第2期空气动力学学报Vol.22,No.2 2004年06月ACTA AERODYNAMICA SINICA Jun.,2004 文章编号:0258-1825(2004)02-0121-09 湍流大涡数值模拟进展 崔桂香,许春晓,张兆顺 (清华大学工程力学系,北京100084) 摘要:本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题。文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径。文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向。 关键词:湍流;大涡数值模拟;亚格子模型;近壁模型;标量湍流 中图分类号:V211.3文献标识码:A * 0引言 复杂流动的准确数值预测是当前航空、航天器研究和设计中迫切需要解决的空气动力学前沿问题之一。随着计算空气动力学方法的不断完善,计算机品质的不断提高,湍流的数值模拟方法成为提高数值预测航天器空气动力特性的瓶颈。 目前,数值预测湍流的方法有:直接数值模拟(DNS),大涡数值模拟(LE S)和雷诺平均模拟(RANS)。 直接数值模拟(DNS)是精确数值模拟湍流的方法,它的主要困难在于湍流是一种不规则多尺度运动,无论在空间上或者时间上湍流都有十分宽广谱。准确数值模拟湍流既要精确计算大尺度流动;又要足够准确地计算小尺度运动。在最简单的各向同性湍流中湍流的空间尺度有以下估计:L ma x/l mi n~Re3/4K,同样最大和最小时间尺度之比T max/t min~Re3/4K,它们都和流动的泰勒雷诺数Re K有关。按照上述估计,空间网格数至少应有:N=N x@N y@N z~Re9/4K,运算量超过Re3K,航空航天器复杂绕流计算的网格数和运算时间远远超过上述估计。因此,目前不具备直接数值模拟复杂工程湍流需要的计算机,湍流直接数值模拟只能作为低雷诺数简单湍流的研究工具。 工程中常用的复杂湍流数值模拟方法是求解雷诺平均的控制方程,这种方法只计算大尺度平均流动,而所有湍流脉动对平均流动的作用,即雷诺应力,用模型假设封闭。由于雷诺应力主要由大尺度脉动贡献,而大尺度脉动和流动的几何特性密切相关,因此雷诺平均模式不是普适的,而是和流动有关,就是说,不存在对一切流动都适用的统一模式;对于不同类型的流动,模式的形式或系数需要修正,而这种修正常常带有经验性。所以,雷诺平均模式不是理想的封闭方法。 湍流大涡数值模拟是有别于直接数值模拟和雷诺平均模式的一种数值预测湍流的方法。这种方法是基于对各种尺度湍流脉动在输运和耗散中作用的认识:大尺度湍流脉动具有主要的能量和动量并支配湍流脉动的动量和能量输运;而湍动能的耗散主要发生在小尺度脉动中[1];根据这一认识产生了湍流大涡数值模拟。它的具体实施方法如下:首先,用滤波方法将小尺度脉动从湍流脉动中去掉,假设空间任意一点的滤波函数为G(x-x0),最简单的滤波器是盒式滤波: G(G)=1,|G i|[$/2, G(G)=0,|G i|>$/2(1)利用滤波器对湍流速度场过滤,过滤后的速度脉动中不存在过滤尺度$以下的脉动成分,称为可解湍流: u i(x,t)= 1 $3 m D u(y,t)G(x-y)d y(2) *收稿日期:2003-03-26;修订日期:2003-06-02. 基金项目:国家自然科学基金资助项目(批准号:10272065,10232020). 作者简介:崔桂香(1950-),女,清华大学工程力学系教授,从事湍流大涡数值模拟和标量湍流的研究.
数值模拟中的湍流模型
流体力学是力学的一个重要分支，它是研究流体（包括液体和气体）这样一个连续介质的宏观运动规律以及它与其他运动形态之间的相互作用的学科，在现代科学工程中具有重要的地位。宏观上讲，黏性流体的流动形态有三种：层流、湍流以及从层流到湍流的转捩。从工程应用的角度看，大多数情况下转捩过程对流体流动的影响不大可以忽略，层流在很少情况下才出现，而在自然界和工程中最普遍存在的是湍流，因此湍流是科学家和工程师研究的重点。湍流理论的研究主要集中在两个方面：一是湍流的触发；二是湍流的描述和湍流问题的求解。 对于工程中出现的湍流问题，其求解方法可归纳为四种：理论分析、风洞实验、现场测试和数值模拟。四种方法相互补充，以风洞实验和现场测试为主，理论分析和数值模拟为辅。数值模拟又称数值风洞，它的出现才十几年却取得迅猛发展，是目前数值计算领域的热点之一，它是数值计算方法、计算机软硬件发展的结果。我们知道，描述流体运动（层流）的流体力学基本方程组是封闭的，而描述湍流运动的方程组由于采用了某种平均（时间平均或网格平均等）而不封闭，须对方程组中出现的新未知量采用模型而使其封闭，这就是CFD中的湍流模型。湍流模型的主要作用是将新未知量和平均速度梯度联系起来。目前，工程应用中湍流的数值模拟主要分三大类：直接数值模拟（DNS）；基于雷诺平均N-S方程组（RANS）的模型和大涡模拟（LES）。 DNS是直接数值求解N-S方程组，不需要任何湍流模型，是目前最精确的方法。其优点在于可以得出流场内任何物理量（如速度和压力）的时间和空间演变过程，旋涡的运动学和动力学问题等。由于直接求解N-S方程，其应用也受到诸多方面的限制。第一：计算域形状比较简单，边界条件比较单一；第二：计算量大。影响计算量的因素有三个：网格数量、流场的时间积分长度（与计算时间长度有关）和最小旋涡的时间积分长度（与时间步长有关），其中网格数量是重要因素。为了得到湍流问题足够精确的解，要求能够数值求解所有旋涡的运动，因此要求网格的尺度和最小旋涡的尺度相当，即使采用子域技术，其网格规模也是巨大的。为了求解各个尺度旋涡的运动，要求每个方向上网格节点的数量与Re3/4成比例，考虑一个三维问题，网格节点的数量与Re9/4成比例。目前，DNS能够求解Re(10e4)的范围。 基于RANS的湍流模型采用雷诺平均的概念，将物理量区分为平均量和脉动量，将脉动量对平均量的影响用模型表示出来。目前，基于RANS方程已经发展了许多模型，几乎能对所有雷诺数范围的工程问题求解，并得出一些有用的结果。其缺点在于：第一：不同的模型解决不同类型的问题，甚至对于同一类型的问题，对应于不同的边界条件需要修改模型的常数；第二：由于不区分旋涡的大小和方向性，对旋涡的运动学和动力学问题考虑不足，不能用来对流体流动的机理进行描述。
湍流的数值模拟
2012年秋季学期研究生课程考核 （读书报告、研究报告） 考核科目高等流体力学 学生所在院（系）机电工程学院 学生所在学科机械制造及自动化学生姓名高强 学号12S008123 学生类别工学硕士 考核结果阅卷人
湍流的数值模拟 一、流体力学概述 流体力学是研究流体的力学运动规律及其应用的学科。主要研究在各种力的作用下，流体本身的状态，以及流体和固体壁面、流体和流体间、流体与其他运动形态之间的相互作用的力学分支。除水和空气之外，这里的流体还指作为汽轮机工作介质的水蒸气、润滑油、地下石油、含泥沙的江水、血液、超高压作用下的金属和燃烧后产生成分复杂的气体、高温条件下的等离子体等等。它的主要基础是牛顿运动定律和质量守恒定律，常常还要用到热力学知识，有时还用到宏观电动力学的基本定律、本构方程和高等数学、物理学、化学的基础知识。气象、水利的研究，船舶、飞行器、叶轮机械和核电站的设计及其运行，可燃气体或炸药的爆炸，汽车制造，以及天体物理的若干问题等等，都广泛地用到流体力学知识。许多现代科学技术所关心的问题既受流体力学的指导，同时也促进了它不断地发展。 二、数值计算在流体力学研究中的应用 数值计算是研究流体力学的重要方法。它是针对流体运动的特点，用数学语言将质量守恒、动量守恒、能量守恒等定律表达出来，从而得到连续性方程、动量方程和能量方程。此外，还要加上某些联系流动参量的关系式(例如状态方程)，或者其他方程。这些方程合在一起称为流体力学基本方程组。 求出方程组的解后，结合具体流动，解释这些解的物理含义和流动机理。通常还要将这些理论结果同实验结果进行比较，以确定所得解的准确程度和力学模型的适用范围。 从基本概念到基本方程的一系列定量研究，都涉及到很深的数学问题，所以流体力学的发展是以数学的发展为前提。反过来，那些经过了实验和工程实践考验过的流体力学理论，又检验和丰富了数学理论，它所提出的一些未解决的难题，也是进行数学研究、发展数学理论的好课题。按目前数学发展的水平看，有不少题目将是在今后几十年以内难于从纯数学角度完善解决的。
湍流理论若干问题研究进展
第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·
化工装置中两相流模型的建立
化工装置中的两相流模型的建立 摘要：通过文献调研，本文重点分析了大涡模型在离心泵两相流中的应用。较为详细的概述了模型的建立以及边界条件的确定和求解方法。 关键词：文献调研、大涡模型、边界条件 前言 两相流动是流体力学中一门重要的分支学科，它在很多现代工程技术甚至医学中得到广泛的应用。可以认为，绝大多数的流动都是多相流动，纯粹的单相流动只是个别情况。降雾，下雨、下冰雹、云层流动、流沙、尘暴等是自然界中两相流动的一些例子。各种发动机和窖炉中的喷雾燃烧、核反应堆的冷却、宇航飞行器的两相绕流、含铝推进剂固体火箭发动机中的燃气流动、石油和天然气的开采和输运、热力设备与制冷系统的工作过程、化学工艺中的流态化、吸收、蒸发、凝结和化学反应过程、采矿和冶金过程中的旋流分离和输运、气力和液力输送、煤的气化和液化、煤粉和煤浆燃烧、空气和水的污染、环保、粉尘爆炸、血液的循环与凝固、水利工程中的泥沙运动和高速渗气流等工程实际问题无不与两相流动有关。离心泵是化工生产中最常见的装置之一，泵内流体的运动以及流体对泵的的磨蚀尤为突出，而两相流动的研究就是为设计泵以及如何防止这些机械磨蚀产生的基础和关键性的内容。近几年，两相流动己发展到与可压缩流体力学及边界层理论有同等重要的地位。因此固液两相流动及多相流动的研究不仅对流体力学的发展，而且对解决工程中的实际问题具有重大的理论价值和实际意义。 下面就离心泵叶轮内高浓度液-固两相湍流的大涡模拟为例阐述化工装置中两相流数学模型的建立、边界条件的确定以及求解方法的选择。 湍流大涡数值模拟（LES）是有别于直接数值模拟和雷诺平均模拟的一种数值模拟手段。利用次网格尺度模型模拟小尺度湍流运动对大尺度湍流运动的影响即直接数值模拟大尺度湍流运动，将N-S方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。 1 大涡模拟 1.1 大涡模拟的基本思想 湍流运动是由许多尺度不同的旋涡组成的。那些大旋涡对于平均流动有比较明显的影响，而那些小旋涡通过非线性作用对大尺度运动产生影响。大量的质量、热量、动量、能量交换是通过大涡实现的，而小涡的作用表现为耗散。流场的形状，阻碍物的存在，对大旋涡有比较大的影响，使它具有更明显的各向异性。小旋涡则不然，它们有更多的共性，更接近各向同性，因而较易于建立有普遍意义的模型。基于上述物理基础，LES把包括脉动运动在内的湍流瞬时运动量通过滤
湍流理论发展概述
湍流理论发展概述 一、湍流模型的研究背景 自然环境和工程装置中的流动常常是湍流流动，模拟任何实际过程首先遇到的就是湍流问题，而湍流问题本身又是流体力学理论上的难题。对于某些简单的均匀时均流场，如果湍流脉动是各向均匀及各向同性的，可以用经典的统计理论来分析，但实际上的湍流往往是不均匀的，这就给理论分析带来了极大地困难。这也就引发了对湍流过程进行模拟的想法。 对湍流最根本的模拟方法是在湍流尺度的网格尺寸内求解瞬态的三维N-S 方程的全模拟方法，此时无需引进任何模型。然而由于计算方法及计算机运算水平的限制，该种方法不易实现。另一种要求稍低的方法是亚网格尺寸度模拟即大涡模拟（LES），也是由N-S 方程出发，其网格尺寸比湍流尺度大，可以模拟湍流发展过程的一些细节，但由于计算量仍然很大，只能模拟一些简单的情况，直接应用于实际的工程问题也存在很多问题[1]。目前数值模拟主要有三种方法：1. 平均N-S方程的求解，2.大涡模拟（LES），3.直接数值模拟（DNS），而模拟的前提是建立合适的湍流模型。 所谓的湍流模型，就是以雷诺平均运动方程与脉动运动方程为基础，依靠理论与经验的结合，引进一系列模型假设，而建立起的一组描写湍流平均量的封闭方程组。目前常用的湍流模型可根据所采用的微分方程数进行分类为：零方程模型、一方程模型、两方程模型、四方程模型、七方程模型等。对于简单流动而言，一般随着方程数的增多，精度也越高，计算量也越大、收敛性也越差。但是，对于复杂的湍流运动，则不一定。湍流模型可根据微分方程的个数分为零方程模型、一方程模型、二方程模型和多方程模型。这里所说的微分方程是指除了时均N-S 方程外，还要增加其他方程才能是方程封闭，增加多少个方程，则该模型就被成为多少个模型。
湍流的数值模拟方法进展
《高等计算流体力学》课程作业 湍流的数值模拟方法进展
1概述 自然环境和工程装置中的流动常常是湍流，模拟任何实际过程首先遇到的就是湍流问题，而湍流问题本身又是流体力学理论上的难题。对于某些简单的均匀时均流场，如果湍流脉动是各向均匀及各向同性的，可以用经典的统计理论来分析，但实际上的湍流往往是不均匀的，给理论分析带来了极大困难。 湍流是空间上不规则和时间上无秩序的一种非线性的流体运动，表现出非常复杂的流动状态，主要表现在湍流流动的随机性、有旋性、统计性。传统计算流体力学中描述湍流的基础是Navier-Stokes（N-S）方程，根据N-S方程中对湍流处理尺度的不同，湍流数值模拟方法主要分为：直接数值模拟（DNS）、雷诺平均方法（RANS）和大涡模拟（LES）。 直接数值模拟可以获得湍流场的精确信息，是研究湍流机理的有效手段，但现有的计算资源往往难以满足对高雷诺数流动模拟的需要，从而限制了它的应用范围。雷诺平均方法可以计算高雷诺数的复杂流动，但给出的是平均运动结果，不能反映流场脉动的细节信息。大涡模拟基于湍动能传输机制，直接计算大尺度涡的运动，小尺度涡运动对大尺度涡的影响则通过建立模型体现出来，既可以得到比雷诺平均方法更多的诸如大尺度涡结构和性质等的动态信息，又比直接数值模拟节省计算量，从而得到了越来越广泛的发展和应用。 2 雷诺平均方法(RANS) 雷诺平均模拟（RANS）即应用湍流统计理论，将非定常的N - S方程对时间作平均，求解工程中需要的时均量。利用湍流模式理论，对Reynolds应力做出各种假设，即假设各种经验的和半经验的本构关系，从而使湍流的平均Reynolds方程封闭。 2.1控制方程 对非定常的N - S 方程作时间演算，并采用Boussinesp 假设，得到Reynolds 方程
大涡模拟简单介绍
《粘性流体力学》小论文 题目：浅谈大涡模拟 学生姓名：丁普贤 学生学号：103911018 完成时间：2010/12/16
浅谈大涡模拟 丁普贤 （中南大学,能源科学与工程学院，湖南省长沙市，410083） 摘要：湍流流动是一种非常复杂的流动，数值模拟是研究湍流的主要手段，现有的湍流数值模拟的方法有三种：直接数值模拟、大涡模拟和雷诺平均模型。本文主要是介绍大涡模拟，大涡模拟的思路是：直接数值模拟大尺度紊流运动，而利用亚格子模型模拟小尺度紊流运动对大尺度紊流运动的影响。大涡模拟在计算时间和计算费用方面是优于直接数值模拟的，在信息完整性方面优于雷诺平均模型。本文还介绍了对N-S方程过滤的过滤函数和一些广泛使用的亚格子模型，最后简单对一些大涡模拟的应用进行了阐述。 关键词：计算流体力学；湍流；大涡模拟；亚格子模型
A simple study of Large Eddy Simulation DING Puxian (Central South University, School of Energy Science and Power Engineering, Changsha, Hunan, 410083) Abstract：Turbulent flow is a very complex flow, and numerical simulation is the main means to study it. There are three numerical simulation methods: direct numerical simulation, large eddy simulation,Reynolds averaged Navier-Stokes method. Large eddy simulation (LES) is mainly introduced in this paper. The main idea of LES is that large eddies are resolved directly and the effect of the small eddies on the large eddies is modeled by subgrid scale model. Large eddy simulation calculation in computing time and cost is superior to direct numerical simulation, and obtain more information than Reynolds averaged Navier-Stokes method. The Navier-Stokes equations filtering filter function and some extensive use of the subgrid scale model are simply discussed in this paper. Finally, some simple applications of large eddy simulation are told. Key words：computational fluid dynamics; turbulence; large eddy simulation; subgrid scale model
湍流的数值模拟综述
湍流的数值模拟 一、引语 流体的流动形态分为湍流与层流。而层流是流体的最简单的一种流动状态。流体在管内流动时，其质点沿着与管轴平行的方向作平滑直线运动。此种流动称为层流或滞流，亦有称为直线流动的。流体的流速在管中心处最大，其近壁处最小。管内流体的平均流速与最大流速之比等于0.5，根据雷诺实验，当雷诺准数引Re<2320时，流体的流动状态为层流。当雷诺数Re>2320时，流体流动状态开始向湍流态转变，湍流是一种很复杂的流动状态，是流体力学中公认的难题。 自从19世纪末O.Reynolds提出湍流的统计理论以来，已经有一个多世纪了，经过几代科学家的努力，湍流研究取得很大进展，但是仍然不能满足工程应用的需要，以至于经常有悲观的论调侵袭湍流研究。为什么湍流问题没有圆满地解决会受到如此关注呢?因为湍流是自然界和工程中十分普遍的流动现豫，对于湍流问题的正确认识和模化直接影响到对自然环境的预测和工程的质量。例如，当前影响航天器气动力和气动热预测准确度的主要障碍是缺乏可靠的湍流模型。和其他一些自然科学的准题不同，解决湍流问题具有迫切性。 湍流运动的最主要特征是不规则性，这是大家公认的。对于湍流不规则性的深入认识，是一百多年来湍流研究的上要成就之一。早期的科学家认为，像分子运动一样，湍流是完全不规则运动。类似于分子运动产生黏性，湍流的耗散可以用涡黏系数来表述。20世纪初，一些杰出的流体力学家，相继对涡黏系数提出各种流体力学的模型，如Taylor(1921年)的涡模型，Praudtl(1925年)的混合长模型和von Karman(1930年)相似模型等。当科学家用流体力学观念(不是分子观念)来建立湍流耗散的涡黏模型时，就开始考虑连续介质不规则运动的特点，其中有别于气体分子不规则运动的最主要特点是运动的多尺度性。第一个提出流体湍流运动中多尺度输运特性的科学家mchardson(1922年)曾描述湍动能的多尺度传输过程如下：“大涡包含小涡，并喂予速度；小涡包含更小的涡，如此继续直到黏性耗散”。多尺度的思想导致产生描述多尺度的谱概念和谱分析方法，并最终产生了Kolmogorov(1941年)的局部各向同性的通用谱(即5/3谱)。 湍流不仅是多尺度的而且是有结构的运动。20世纪中叶，大量的湍流实验(包括测量和显示)发现多尺度的湍流运动存在某种特殊的运动状态。Townsend(1951年)，Corrsin(1955年)和Lumley(1965年)等从脉动序列的间歇性和空间相关相继推测湍流结构的可能形态。理论上也提出过各种湍涡的模型：球涡模型，柱涡模型等。早期的湍流结构主要是从运动学上考虑，把旋涡结构作为湍流统计的样本。我国的周培源教授是近代湍流模式的奠基人之一，他首先提出先解方程后平均的统计方法，就是说湍涡必须满足Navier—Stokes方程(Chou and Chou，1995年)。 真实的、可以观察到的湍流结构通过流动显示，以及稍后湍流直接数值模拟所证实。典型的例子是混合层的Brown—Roshko涡(1976年)，图1明显地展示了混合层中存在规则的大涡和分布在大涡周围的细小湍涡。在边界层、槽道和圆管湍流中也存在各式各样的大涡结构。例如，用激光诱导荧光的显示方法，我们可以在圆管湍流中观察到周向（图2a）和流向大涡（图2b）。值得提出的是，不仅在剪切湍流中有大涡结构，简单的均匀各向同性湍流中也存在涡结构。图3展示的是各向同性湍流的直接数值模拟中强涡量等值面，它们是管状结构。仔细分析还可以确定管状涡的平均长度约等于各向同性湍流的积分尺度，它们的平均直径约等于湍流TayLor微尺度，更进一步分析可以算出管状涡内部的平均速度
关于湍流理论研究进展
关于湍流理论研究进展 摘要本文对近年来湍流理论在某些方面的研究进展作了概要介绍，对具有代表性的理论假设的思想方法，进行了扼要阐述，指出了相应的实用价值和局限性。 关键词湍流湍流统计理论混沌理论湍流拟序结构湍流剪切流动 1 无处不在的湍流现象 湍流是自然界中流体的一种最普遍的运动现象，它广泛的存在于我们生活周围。在大风吹过地面障碍物的旁边，在湍急的河水流过桥墩的后面，在烟囱中冒出的浓烟随风渐渐扩散等地方，都能观察到湍流运动现象。简单地说，湍流运动就是流体的一种看起来很不规则的运动。由于湍流现象广泛存在于自然界和工程技术的各个领域，因此湍流基础理论研究取得的进展就可能为经济建设和国防建设的广泛领域带来巨大的效益。例如，提高各种运输工具的速度以大量节约能源，提高各种流体机械的效益；改善大气和水体的环境质量，降低流体动力噪声，防止流体相互作用引发的结构振动乃至破坏；加强反应器内部物质的热交换与化学反应的速度等等。 然而像湍流这样，虽经包括许多著名科学家在内长达一个世纪多的顽强努力，正确反映客观规律的系统的湍流理论至今还没有建立，在整个科学研究史上也是不多见的。因此，可以说湍流是力学中没有解决的最困难的难题之一。因此，世界上许多国家一直坚持把湍流研究列为需要最优先发展的若干重大基础研究课题之一。 2 湍流理论的发展历史 湍流理论从它的思路来说大体可分为两类[1]。一类是先把流体动力学方程组平均以后，然后再设法使方程组封闭，求解后再和实验结果比较，看封闭办法是否正确。湍流中绝大部分理论是属于这一类型。另一类是先求解，取特殊模型，再引进平均，得到要求的物理量，和相应的实验结果进行比较。 2.1 Reynolds方程和混合长度理论 十九世纪70年代是Maxwell-Boltzmann分子运动理论取得辉煌成果的时代。它成功地解释了气体状态方程、气体粘性、气体热传导和气体扩散等一系列现象。湍流理论开始发展的时候，就受着这种思想支配。1877年T.V.Bonssinesq[2]又开始
湍流燃烧数值模拟研究
湍流燃烧及其数值模拟研究 1. 湍流燃烧 1.1 湍流燃烧基本概念 当流动雷诺数数较小时，由于流体粘性的作用，流体呈层流流态。当流动的特征雷诺数超过相应的临界值，流动从层流转捩到湍流。湍流燃烧是指湍流流动 中可燃气的燃烧，在能源、动力、航空和航天等工程领域,经常遇到的实际燃烧过程几乎全部都是湍流燃烧过程。湍流燃烧实质是湍流，化学反应和传热传质等过程相耦合的结果。湍流对燃烧的影响与湍流强度和湍流涡旋尺度有关。小尺度湍流通过湍流扩散使火焰区内的输运效应增加，从而使化学反应速率增加。但气流脉动不会火焰面产生皱褶，只能把火焰变成波纹状。大尺度湍流对火焰内部结构没有影响，但使火焰阵面出现皱褶，增加其燃烧面积，造成火焰表现传播速度增加。当湍流强度及湍流尺度均较大时，火焰前沿不再连续而分裂成四分五裂。 燃烧对湍流的影响主要表现在燃烧释放的热流流团膨胀，影响气体的密度和运动速度，从而影响当地的涡旋，湍流强度和湍流结构。 1.2 湍流燃烧分类湍流燃烧按其燃料和氧化剂的初始混合状态可以分类为：湍流非预混燃烧、预混燃烧和部分预混燃烧。在湍流非预混燃烧燃料和氧化剂事先是分离的，燃料和氧化剂一边混合一边燃烧，燃烧速率主要受湍流混合过程控制，而在湍流预混燃烧中，燃料和氧化剂在进入核心燃烧区以前已经充分混合，化学反应的速率由火焰前缘从炽热的燃烧区向冷态无反应区的传播所控制。上面两种燃烧方式是湍流燃烧的两个极限情形，很多情况下两种燃烧模式是并存的，称为部分预混燃烧。部分预混燃烧可出现在下列情形中叫：(1) 在一个完全以非预混燃烧为配置的燃烧装置发牛了局部熄火；(2) 当预混火焰前缘穿过非均匀的混气时；(3) 射流非预混火焰发生抬举，其根部是一。个典型的部分预混火焰。这三种部分预混燃烧情形涉及了经常受到关注的燃烧研究话题如局部熄火、火焰稳定等，它们对研究湍流燃烧过程的机理有很大意义。 在湍流燃烧中,湍流流动过程和化学反应过程有强烈的相互关联和相互影响. 湍流通过强化混合而影响着时平均化学反应速率,同时化学反应放热过程又影响着湍流,如何定量地来描述和确定这种相互作用是湍流燃烧研究的一个重要内容. 湍流是非常复杂的，它包括湍流问题，湍流与燃烧的相互作用，流动参数与
流体力学的研究进展
流体力学的研究进展 环境工程01班张东元 20106370 摘要： 利用流体力学研究解决不同方面和不同角度的问题，主要是构建数学模型，利用数学模型预测的结果进行分析研究实际情况，从而解决实际问题。 关键词：流体、数学模型、数值模拟、预测 简介： 流体力学是力学的一个分支，它研究流体静止和运动的力学规律，及其在工程技术中的应用。流体力学的研究进展，就是其相关方面的研究的前沿动态。 正文： 流体力学中研究得最多的流体是水和空气。它的主要基础是牛顿运动定律和质量守恒定律，常常还要用到热力学知识，有时还用到宏观电动力学的基本定律、本构方程和物理学、化学的基础知识。 目前来说，流体力学在供热通风和燃气工程中应用的非常广泛。热的供应，空气的调节，燃气的输配，排毒排湿，除尘降温等等，都是一流体作为工作介质，通过流体的各种物理作用，对流体的流动有效地加以组织来实现的。除了上述的两个方面，在地质考察，环境监测，气象观测及海洋评估等多个领域中流体力学都有或多或少的涉猎。下面本文就针对其中几个方面来做个概述。
在环境流体力学方面而言，其本身建立之初，主要是为了研究污染问题，所涉及的尺度小、范围窄，不能在宏观上、整体上把握环境要素的变化规律。近年来，随着流体力学的发展与突破，环境流体力学逐渐趋向于研究多尺度、多科学的综合问题。 在河流水体方面，主要研究水利学中的环境问题，比如像05年的松花江苯类物质污染事件，其中的污染主体是100吨左右的苯类污染物，而根据流体力学的知识，黑龙江省水文局的工程师判断出了污染水团留出哈尔滨市区范围及中国范围的时间，并以此为后续处理工作的计划做了时间上的预算。这里面就涉及了一个叫水质数学模型的概念，是水体中污染物随时间和空间迁移转化规律的描述，表示一个定量关系。从文献中可知，GIS（地理信息系统）对构造这一模型起了十分重要的作用，这也是流体力学中科技的运用。 在大气环境中，流体力学一般研究居住环境的空气质量问题。随着计算机资源的迅速改善，用现代计算流体力学方法数值模拟大气环境效果显著。对于不同的尺度范围，研究对象的物理条件不同，数值模拟方法也随之不同。在城市尺度和居民小区尺度大气环境中，采用大涡数值模拟方法等方法。数值计算的边界条件是实现数值预测城市大气环境的关键，其中包括高空边界条件、计算域侧面的边界条件和下垫面边界条件。在预测中，选择合理的模型就能获得较为接近真实值的预测结果。这种数值模拟能得到所测地的风速风向、颗粒物大小及污染物浓度等方面的预测值，而不需要实地地观测考察，只需要单纯的模拟计算就行了。这样一来，不仅仅是单纯的居住环境的大气研