湍流直接数值模拟
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)。
直接数值模拟可以获得湍流场的精确信息,是研究湍流机理的有效手段,但现有的计算资源往往难以满足对高雷诺数流动模拟的需要,从而限制了它的应用范围。
雷诺平均方法可以计算高雷诺数的复杂流动,但给出的是平均运动结果,不能反映流场紊动的细节信息。
大涡模拟基于湍动能传输机制,直接计算大尺度涡的运动,小尺度涡运动对大尺度涡的影响则通过建立模型体现出来,既可以得到较雷诺平均方法更多的诸如大尺度涡结构和性质等的动态信息,又比直接数值模拟节省计算量,从而得到了越来越广泛的发展和应用。
2 直接数值模拟(DNS)湍流直接数值模拟(DNS)就是不用任何湍流模型,直接求解完整的三维非定常的N - S 方程组,计算包括脉动在内的湍流所有瞬时运动量在三维流场中的时间演变。
cfd中湍流模型与控制方程

CFD中湍流模型与控制方程在计算流体动力学(CFD)中,湍流模型和控制方程是非常重要的概念。
湍流模型:湍流是一种高度复杂、非线性的流体运动状态,其特点是流体中的速度、压力等物理量随时间和空间发生随机变化。
为了模拟湍流,需要采用湍流模型。
湍流模型通常分为两类:直接数值模拟(DNS)和非直接数值模拟。
1.直接数值模拟(DNS):DNS直接求解Navier-Stokes方程,不需要对湍流进行任何假设或简化。
然而,由于湍流的多尺度特性,DNS需要极高的计算资源,因此在实际应用中受到限制。
2.非直接数值模拟:为了降低计算成本,非直接数值模拟方法被广泛应用。
这些方法包括雷诺平均法(RANS)、大涡模拟(LES)和统计平均法(SAS)等。
这些方法通过对湍流进行某种程度的平均或滤波,将湍流分解为可解析的大尺度运动和需要模型化的小尺度运动。
控制方程:在CFD中,流体的运动遵循基本的物理定律,如质量守恒定律、动量守恒定律和能量守恒定律。
这些定律在数学上表现为一系列偏微分方程,称为控制方程。
1.质量守恒方程(连续性方程):描述流体微元的质量不随时间变化,即流体微元的质量流入率等于其质量流出率。
在不可压缩流体中,连续性方程简化为速度场的散度为零。
2.动量守恒方程(Navier-Stokes方程):描述流体微元的动量不随时间变化,即流体微元的动量流入率加上外力等于其动量流出率。
Navier-Stokes方程是流体动力学的基本方程,描述了流体运动的基本规律。
3.能量守恒方程:描述流体微元的能量不随时间变化,即流体微元的能量流入率加上外力做功和热源等于其能量流出率。
在不可压缩流体中,能量守恒方程通常简化为温度场的热传导方程。
在求解这些控制方程时,需要选择合适的湍流模型来封闭方程组,以便进行数值求解。
不同的湍流模型和控制方程组合可以适用于不同的流体流动场景,如层流、湍流、可压缩流体、不可压缩流体等。
流体力学中的多尺度湍流模拟与建模

流体力学中的多尺度湍流模拟与建模湍流是流体力学中一个复杂而普遍存在的现象,涉及到多尺度的运动和相互作用。
在实际应用中,对湍流进行准确模拟和有效建模具有极大的重要性。
本文将介绍流体力学中的多尺度湍流模拟与建模方法,并探讨其在工程实践中的应用。
第一部分:湍流模拟方法湍流模拟是通过数值方法模拟湍流流动,以获得流场的详细信息。
在多尺度湍流模拟中,常用的模拟方法包括直接数值模拟(DNS)、大涡模拟(LES)、雷诺平均导数模拟(RANS)等。
直接数值模拟是一种最为精确的模拟方法,通过求解流动的Navier-Stokes方程来模拟湍流现象。
由于湍流流动存在广泛的空间和时间尺度,直接数值模拟的计算成本极高,通常只能用于精细的研究和小规模的流动模拟。
大涡模拟是在直接数值模拟的基础上发展起来的一种方法,通过将大涡的运动精确模拟,而对小涡采用模型进行参数化。
相比于直接数值模拟,大涡模拟的计算成本较低,可以在一定程度上模拟湍流的多尺度特性。
雷诺平均导数模拟是一种更为常用的湍流模拟方法,在工程实践中得到广泛应用。
该方法通过将流场的各个变量进行平均处理,然后引入湍流模型来描述湍流效应。
由于雷诺平均导数模拟只考虑了平均尺度上的湍流特性,无法准确模拟湍流的具体结构,因此在一些对流动细节要求较高的场合,该方法的精度有限。
第二部分:湍流建模方法湍流建模是为了在湍流模拟中描述湍流效应而引入的方法。
这些模型基于湍流的统计性质和物理规律,对湍流的各种参数进行描述和计算。
常用的湍流建模方法包括湍流能量方程、湍流应力传输方程等。
湍流能量方程是湍流建模中的一种重要方法,用于描述湍流的能量传输过程。
该方程通过考虑湍流的产生、消耗和传输等过程,以及湍流能量的耗散来描述湍流的演化规律。
基于湍流能量方程,可以计算湍流的能谱和湍流能量的分布等参数。
湍流应力传输方程是湍流建模中的另一种关键方法,用于描述湍流的动量传输过程。
该方程通过考虑湍流的各向异性和湍流的剪切作用等因素,计算湍流应力的分布和演化规律。
湍流数值模拟及其在工程热力学中的应用

湍流数值模拟及其在工程热力学中的应用湍流是自然界和工程中广泛存在的一种流动状态,其具有不规则、不稳定、非线性等特点。
因此,湍流研究成为了流体力学中的一个重要分支。
湍流数值模拟(Large Eddy Simulation)是目前研究湍流问题的重要手段之一,广泛应用于工程热力学中。
湍流数值模拟技术的发展历程湍流数值模拟技术起源于20世纪50年代,当时主要应用于理论模拟。
20世纪80年代后,随着计算机技术的发展,数值模拟技术应用于实际工程中,并得到广泛应用。
近年来,由于计算机性能的不断提高和算法的不断改进,湍流数值模拟技术越来越成熟,其应用范围也更加广泛。
湍流数值模拟技术的基本原理湍流数值模拟技术的基本原理是将流场分为宏观湍流和微观湍流两部分,并通过不同方法对二者进行模拟。
具体而言,宏观湍流采用平均场方程进行模拟,微观湍流则通过小尺度涡结构之间的相互作用进行模拟。
在湍流数值模拟过程中,关键是要准确地描述湍流的能量转移和钝化机制,以便合理地模拟湍流特性。
目前,湍流数值模拟技术主要有两种方法:直接数值模拟和大涡模拟。
直接数值模拟(Direct Numerical Simulation,DNS)是最为精确的湍流数值模拟方法,它直接求解完整的Navier-Stokes方程,但计算量也是最大的。
而在工程应用中,一般采用次网格模型,采用模型对小尺度湍流进行近似处理,减少计算量。
其中,大涡模拟(Large Eddy Simulation,LES)是一种很有代表性的方法,它将外部湍流场分解为大尺度湍流和小尺度湍流两部分,对大尺度湍流进行直接数值模拟,对小尺度湍流采用模型进行处理。
湍流数值模拟在工程热力学中的应用湍流数值模拟技术在工程热力学中有着广泛的应用。
具体而言,湍流数值模拟可以用来模拟涡流管道的流动、火焰、燃烧室和喷气发动机等复杂流场问题。
下面,我们将从两个方面来介绍湍流数值模拟在工程热力学中的应用:(1)流体力学问题湍流数值模拟技术在流体力学问题中得到了广泛应用,例如现代汽车设计中对车身和车厢空气动力学的研究,对于气动设计、噪声控制和气密性等方面的分析有很大的帮助。
流体力学实验装置的流体流动湍流分析方法

流体力学实验装置的流体流动湍流分析方法流体力学是研究流体在各种条件下的力学性质和运动规律的科学。
在流体力学实验中,对流体流动的湍流现象进行分析是非常重要的一项工作。
本文将探讨流体力学实验装置中流体流动湍流分析的方法。
一、流体流动的湍流现象在流体力学实验中,流体流动的湍流现象是非常常见的。
湍流是一种不规则而复杂的流动现象,其特点是流体中的速度、密度、温度等参数呈现不规则的瞬时变化。
湍流会导致流体的速度场、压力场等出现不规则的分布,使流体运动变得难以预测和分析。
二、湍流分析的基本方法在流体力学实验中,分析流体流动的湍流现象可以采用各种方法。
常用的湍流分析方法包括直接数值模拟(DNS)、大涡模拟(LES)、雷诺平均法(RANS)等。
这些方法各有特点,适用于不同的流动情况和要求。
1. 直接数值模拟(DNS)直接数值模拟是一种计算流体动力学的方法,通过数值模拟计算流体在空间和时间上的变化。
DNS方法需要考虑流场中所有长度尺度的细节,对计算资源要求较高,适用于小尺度流动情况的湍流分析。
2. 大涡模拟(LES)大涡模拟是一种通过滤波方法将湍流分解成大尺度涡和小尺度湍流的方法。
LES方法可以利用数值模拟计算大尺度涡的运动规律,通过模型计算小尺度湍流的效应。
LES方法适用于中等尺度流动情况的湍流分析。
3. 雷诺平均法(RANS)雷诺平均法是一种通过时间平均和空间平均将湍流分解成平均流场和湍流脉动的方法。
RANS方法通过求解平均流动方程和湍流模型方程,可以得到流场的平均参数。
RANS方法适用于大尺度流动情况的湍流分析。
三、实验装置中流体流动湍流分析的应用流体力学实验装置中流体流动湍流分析的方法应用广泛。
通过对流体流动的湍流现象进行分析,可以了解流场中的湍流结构、湍流能量传输规律等重要信息,为实验设计和数据处理提供依据。
例如,在风洞实验中,通过对流体流动的湍流现象进行模拟和分析,可以评估空气动力学性能、气流湍流特性等参数,为风洞实验的设计和优化提供参考。
高精度差分格式及湍流数值模拟(三)

~m; (~10m)
最小尺度涡: 微米量级; 中小尺度涡: 0.1mm~ 1mm量级
直接数值模拟: 分辨出最小尺度涡; 网格量极为巨大: 计算量、存储量大
对高性能计算需求强烈
二、 高精度CFD软件 OpenCFD
OpenCFD: 作者开发的一套高精度、开放的CFD程序
1) OpenCFD-SC : 高精度差分 2) OpenCFD-EC: 有限体积 3) OpenCFD-Comb: 化学反应计算
computational region
200
z 300
outflow boundary
buffer region
400
500
示例2:有攻角小头钝锥边界层转捩的DNS
3) 计算结果验证
a.不同网格之间的比较 (网格收敛性) b. 不同扰动形式的结果比较 (模型正确性) d. 与Horvath静风洞实验比较 (与实验结果比较) c. 与 eN方法比较 (与理论结果比较) e. 与Stenson实验比较 (与实验结果比较)
面向工程计算开放的CFD代码(Open CFD code for Engineering Computing)
算法: 多块结构网格有限体积 + RANS 特点: 适用于复杂外形工程计算
差分-有限体积混合方法
3) 高精度化学反应模拟程序—— OpenCFD-Comb
面向化学反应的开放CFD程序: An OpenCFD code for Combustion
高精度差分格式及 湍流数值模拟 (三)
Part 3 可压缩湍流DNS
1. 背景 2. OpenCFD软件简介 3. 典型可压缩湍流的直接数值模拟示例 4. 湍流模型的评估及改进 5. 湍流燃烧的DNS 6. 小结
湍流模型介绍

湍流模型目前计算流体力学常用的湍流的数值模拟方法主要有以下三种:直接模拟(direct numerical simulation, DNS)直接数值模拟(DNS)特点在湍流尺度下的网格尺寸内不引入任何封闭模型的前提下对Navier-Stokes方程直接求解。
这种方法能对湍流流动中最小尺度涡进行求解,要对高度复杂的湍流运动进行直接的数值计算,必须采用很小的时间与空间步长,才能分辨出湍流中详细的空间结构及变化剧烈的时间特性。
基于这个原因,DNS目前仅限于相对低的雷诺数中湍流流动模型。
另外,利用DNS模型对湍流运动进行直接的数值模拟对计算工具有很高的要求,计算机的内存及计算速度要非常的高,目前DNS模型还无法应用于工程数值计算,还不能解决工程实际问题。
大涡模拟(large eddy simulation, LES)大涡模拟(LES)是基于网格尺度封闭模型及对大尺度涡进行直接求解N-S方程,其网格尺度比湍流尺度大,可以模拟湍流发展过程的一些细节,但其计算量仍很大,也仅用于比较简单的剪切流运动及管流。
大涡模拟的基础是:湍流的脉动与混合主要是由大尺度的涡造成的,大尺度涡是高度的非各向同性,而且随流动的情形而异。
大尺度的涡通过相互作用把能量传递给小尺度的涡,而小尺度的涡旋主要起到耗散能量的作用,几乎是各向同性的。
这些对涡旋的认识基础就导致了大涡模拟方法的产生。
Les大涡模拟采用非稳态的N-S方程直接模拟大尺度涡,但不计算小尺度涡,小涡对大涡的影响通过近似的模拟来考虑,这种影响称为亚格子Reynolds应力模型。
大多数亚格子Reynolds模型都是将湍流脉动所造成的影响用一个湍流粘性系数,既粘涡性来描述。
LES对计算机的容量和CPU的要求虽然仍然很高,但是远远低于DNS方法对计算机的要求,因而近年来的研究与应用日趋广泛。
应用Reynolds时均方程(Reynolds-averaging equations)的模拟方法许多流体力学的研究和数值模拟的结果表明,可用于工程上现实可行的湍流模拟方法仍然是基于求解Reynolds时均方程及关联量输运方程的湍流模拟方法,即湍流的统观模拟方法。
湍流的几种数值模拟方法

LES特点
抓大不放小 非常有利,有力的工具 是最近,可预见未来流体 力学研究和应用的热点 近来又出现了VLES, DES等在LES上发展而 来的工具
Will RANS survive LES? Hanjalic自问自答
会。Journal of Fluids Engineering -V127, 5, pp. 831-839 (Will RANS
Prandtl(1925)混合长度模型
也被称作零方程模型 还在被广泛应用 廉价,易收敛 基本在流场比较简单,或者对计算结果 精度要求不高或者流场形状比较复杂的 行业中,比如暖通空调,流体机械等。
Prandtl混合长度模型 缺点
最明显的缺点是:当速度梯度 为零的 时候, 消失, 这与事实不符
Launder and Li(1994), Craft and Launder (1995)
目前有很多学者在继续此方面的工作
Brian E. Launder
本科Imperial College, London 硕博 MIT 实验流体力学 1964-1976 Imperial College 讲师
涡流粘度
Eddy viscosity or turbulent viscosity
二维流场分子粘性力
为描述雷诺应力,Boussinesq 1887 定义了与之相对应的
RANS模型的核心在于给出 的数 学表达式,要求精度高,适用范围广
涡流粘度,
Prandtl 1925 Prandtl 1945 Bradshaw 1968 Kolmogorov, 1942 Hanjalic 1970 Rotta 1951 Chou 1945 Davidov 1961
湍流的数值模拟方法进展

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 ,即'+=i i i u u u ,i u 可以采用Leonard 提出的算式表示为:(11)式中)(x x G '-称为过滤函数,显然G(x)满足x d x u x x G x u i i '''-=⎰+∞∞-)()()(⎰+∞∞-=1)(dx x G3.3 控制方程将过滤函数作用与N —S 方程的各项,得到过滤后的湍流控制方程组:由于无法同时求解出变量i u 和j i u u ,所以将j i u u 分解成i j i j ij u u u u τ=⋅+,ij τ即称为次网格剪切应力张量(亦称为亚格子应力)。
流体的湍流模型和湍流模拟

流体的湍流模型和湍流模拟流体力学是研究流体的运动规律和性质的学科,其中湍流模型和湍流模拟是其中非常重要的研究方向。
湍流是流体力学中一种复杂而普遍存在的现象,它具有不规则、无序和随机性等特点。
湍流模型和湍流模拟的发展,对于理解和预测真实世界中的湍流现象,以及涉及湍流的工程设计和应用具有重要意义。
一、湍流模型湍流模型是描述湍流现象的数学模型,在流体力学中起着扮演着非常重要的作用。
根据流体力学理论,湍流是由于流体中微小尺度的速度涡旋突然出现和消失所导致的现象。
由于湍流涡旋的尺度范围很广,从而难以直接模拟和计算。
因此,使用湍流模型来近似描述湍流现象,成为了一种常用的方法。
常见的湍流模型包括雷诺平均湍流模型(Reynolds-averaged Navier-Stokes equations, RANS)和大涡模拟(large eddy simulation, LES)等。
雷诺平均湍流模型是基于平均流场的统计性质,通过求解雷诺平均速度和湍流应力来评估湍流效应。
而大涡模拟是将湍流现象分解为不同尺度的涡旋,并通过直接模拟大涡旋来研究湍流运动。
二、湍流模拟湍流模拟是利用计算机来模拟湍流现象的方法,通常基于数值方法对流体力学方程进行求解。
湍流模拟分为直接数值模拟(direct numerical simulation, DNS)、雷诺平均湍流模拟和大涡模拟等。
直接数值模拟是将流场划分为网格,并通过离散化流体力学方程和湍流模型来求解湍流流场的详细信息。
由于该方法需要计算微小尺度的细节,计算量非常大,限制了其在实际工程中的应用。
因此,直接数值模拟主要用于湍流现象的基础研究和理论验证。
相比之下,雷诺平均湍流模拟和大涡模拟能够更有效地模拟湍流现象。
雷诺平均湍流模拟通过对湍流参数进行求解,来描述平均的湍流效应。
而大涡模拟则将湍流现象分为大涡旋和小涡旋,通过模拟大涡旋来捕获湍流流场的主要特征。
三、湍流模型与湍流模拟的应用湍流模型和湍流模拟在工程设计和应用中有着广泛的应用。
工程流体力学中的湍流模型比较与分析

工程流体力学中的湍流模型比较与分析引言:湍流是流体力学中一种复杂的流动现象,它广泛存在于自然界和工程应用中。
研究和模拟湍流流动是工程流体力学中的一个重要课题。
湍流模型是用来描述湍流流动的数学模型,对于工程实践中的湍流模拟有着重要的影响。
本文将比较和分析几种常用的湍流模型,包括雷诺平均Navier-Stokes方程(RANS)模型、大涡模拟(LES)和直接数值模拟(DNS)。
1. 雷诺平均Navier-Stokes方程(RANS)模型雷诺平均Navier-Stokes方程是湍流模拟中最常用的模型之一。
它基于雷诺平均的假设,将流动场分解为平均流动和湍流脉动两部分。
RANS模型通过求解平均流动方程和湍流脉动方程来描述流场的平均状态和湍流效应。
经典的RANS模型包括k-ε模型和k-ω模型,它们通过引入湍流能量和正应力来描述湍流的传输和衰减。
2. 大涡模拟(LES)大涡模拟是一种介于RANS模型和DNS模型之间的模型。
在LES模拟中,较大的湍流涡旋被直接模拟,而较小的涡旋则通过子网格模型(subgrid model)来描述。
LES模型可以较好地模拟湍流的空间变化特性,对于流动中的尺度较大的湍流结构有着较好的描述能力。
然而,由于需要模拟较小的湍流结构,LES模拟通常需要更高的计算资源和更复杂的数值算法。
3. 直接数值模拟(DNS)直接数值模拟是一种最为精确的湍流模拟方法,它通过直接求解包含所有空间和时间尺度的Navier-Stokes方程来模拟湍流流动。
DNS模拟可以精确地捕捉湍流流动中的所有涡旋和尺度结构,提供最为详细的湍流统计信息。
然而,由于湍流流动具有广泛的空间和时间尺度,DNS模拟通常需要巨大的计算资源和较长的计算时间。
4. 模型比较与选择在实际工程应用中,选择合适的湍流模型需要综合考虑计算资源、计算效率和模拟精度。
如果在工程实践中仅关注流场的整体特征和平均效应,RANS模型是一种简便且有效的选择,尤其是k-ε模型和k-ω模型在工程应用中得到了广泛的应用。
流体力学中的流体中的湍流模拟技术

流体力学中的流体中的湍流模拟技术流体力学中的流体湍流模拟技术为了更好地理解和研究流体行为,科学家和工程师在许多领域,如航空航天、海洋工程和能源领域等,依赖于流体力学。
流体力学研究中一个重要的问题是湍流现象的模拟。
本文将探讨流体力学中的流体湍流模拟技术。
一、湍流的概念和特征湍流是流体力学中的一种复杂流动状态,其特点是流速和压力的瞬时变化,无规则的旋涡结构以及尺度的不确定性。
湍流对于流体力学来说是一个挑战,因为湍流过程难以解析地描述。
二、传统的湍流模拟方法1. 直接数值模拟(Direct Numerical Simulation,DNS):DNS是一种通过求解雷诺平均套用的纳维尔-斯托克斯方程来模拟湍流的方法。
然而,DNS需要非常细密的网格以捕捉湍流涡旋的小尺度结构,因此计算量非常大。
2. 大涡模拟(Large Eddy Simulation,LES):LES是一种将湍流流体划分为大尺度涡旋和小尺度涡旋的方法。
大尺度涡旋通过求解时间平均的纳维尔-斯托克斯方程来模拟,小尺度涡旋通过模型来近似。
LES在一定程度上减少了计算量,但仍然需要较精细的网格来解决小尺度涡旋。
三、基于计算流体力学的湍流模拟方法随着计算机技术的发展,计算流体力学(Computational Fluid Dynamics,CFD)成为了模拟湍流的重要工具。
CFD基于数值方法对流体力学方程进行离散求解,可以模拟复杂的湍流流动。
1. 雷诺平均湍流模型(Reynolds-Averaged Navier-Stokes,RANS):RANS是CFD中最常用的湍流模拟方法。
它通过对流体力学方程进行时间平均和空间平均,然后引入湍流模型来描述整体的湍流效应。
RANS方法计算量相对较小,适用于许多工程应用。
2. 湍流模型的改进与发展:针对RANS方法在湍流模拟中存在的局限性,研究人员提出了许多改进的湍流模型。
如雷诺应力输运模型(Reynolds Stress Transport Model,RSTM)和湍动能方程模型(Turbulent Kinetic Energy,TKE)等。
水力学中的湍流流场数值模拟方法

水力学中的湍流流场数值模拟方法湍流是自然界中最常见的流动现象之一,它不仅出现在河道中,也出现在空气中、海洋中等自然环境中。
湍流带有不规则、无序的运动形式,可以将能量从大尺度输送到小尺度。
然而,湍流流场的物理机理十分复杂,难以通过实验和经验来全面理解和研究。
因此,采用数值模拟方法来模拟湍流流场已成为一种重要的研究手段。
本文将介绍目前水力学中常用的湍流流场数值模拟方法,包括雷诺平均NAVIER-STOKES方程模型(RANS)、大涡模拟(LES)、直接数值模拟(DNS)等。
1. 雷诺平均NAVIER-STOKES方程模型(RANS)RANS是目前水力学中常用的湍流流场数值模拟方法,它的基本思想是用平均流动变量来描述湍流流场,从而将部分湍流运动视为均匀的分析。
RANS假设流场中的湍流运动呈现稳定流动形式(平均流动),模拟平均流动状态,再通过额外的方程组描述湍流中的脉动变化,求解平均流动和湍流脉动变化的复合方程。
RANS方法否认任何尺度上的湍流结构,其主要适用于稳态的湍流运动,如河流、管道流动等。
2. 大涡模拟(LES)与RANS不同,LES方法重点关注大尺度上的湍流结构,将湍流流场分解成大尺度流动和小尺度结构,对大尺度结构进行数值模拟,对小尺度结构进行忽略(或近似处理)。
因此,LES适用于小尺度结构对大尺度流动影响较显著的湍流流场,例如紊流冲击波、湍流尾涡等。
在LES中,大尺度上的湍流结构通过方程组求解,而小尺度上的结构则需借助湍流模型的辅助说明。
由此,需要找到适合模拟大尺度流动和小尺度结构的模型参数。
3. 直接数值模拟(DNS)与RANS和LES不同,DNS方法直接模拟所有尺度上的湍流结构,没有任何参数模型的干扰,相比其他两种方法更加精确和准确。
但DNS需要在计算机模拟中处理每个细节,内存和处理能力的要求比较高。
因此DNS目前仅应用于小尺度流动的研究,例如涡街、微小水滴的湍流等。
综上所述,湍流流场数值模拟是研究湍流流场运动机理的重要手段。
fluent 湍流数值模拟方法

fluent 湍流数值模拟方法湍流数值模拟方法是一种在计算流体力学中用于模拟湍流现象的数值方法。
它能够通过数值模拟来预测流体中的湍流运动和湍流相关的现象。
湍流是一种复杂的流动状态,它涉及到大范围的速度和压力变化,以及涡旋的产生和传输。
传统的流体力学模型往往难以处理湍流问题,因此需要采用数值模拟方法。
湍流数值模拟方法通常基于雷诺平均的Navier-Stokes方程,即雷诺平均Navier-Stokes方程(RANS)。
该方程将流场分解为平均部分和脉动部分,并通过对平均部分施加平均性质来减小方程的复杂性。
湍流数值模拟方法的核心是湍流模型。
湍流模型是用来描述湍流运动和湍流相关现象的数学模型,常见的湍流模型有雷诺平均应力模型(Reynolds-Averaged Navier-Stokes model,简称RANS模型)和大涡模拟(Large Eddy Simulation,简称LES)模型。
RANS模型是常用的湍流模型,它基于雷诺应力的概念,通过近似计算湍流运动的影响。
RANS模型通常采用湍流能量方程和湍流运输方程来描述湍流的能量传输和湍流量的传输。
LES模型是一种更为详细的湍流模型,它不仅考虑了雷诺平均流动的影响,还能够模拟湍流中的大尺度湍流涡旋。
LES模型通常通过将流场分解为大尺度和小尺度湍流结构来描述湍流运动。
湍流数值模拟方法的求解过程通常分为网格生成、离散化、求解和后处理四个步骤。
其中网格生成是将流场划分为离散网格的过程,离散化是将Navier-Stokes方程离散成代数方程组的过程,求解是通过迭代计算求解方程组的过程,后处理是对计算结果进行分析和可视化的过程。
湍流数值模拟方法是研究湍流现象和湍流相关问题的重要工具,它在航空航天、汽车工程、能源领域等众多领域具有广泛应用。
流体力学中的湍流模拟方法比较与评估

流体力学中的湍流模拟方法比较与评估引言:湍流是流体力学领域中一个重要且复杂的现象,在自然界和工程应用中都普遍存在。
由于湍流的不稳定性和高度的非线性特性,准确预测和模拟湍流是一个具有挑战性的问题。
因此,为了更好地理解湍流的性质和行为,并预测其对工程应用的影响,研究人员开发了多种湍流模拟方法。
本文将对流体力学中常用的湍流模拟方法进行比较与评估。
一、直接数值模拟(DNS)方法直接数值模拟(DNS)是一种较为精确的湍流模拟方法。
该方法通过解析求解Navier-Stokes方程,将湍流现象的所有空间和时间尺度都考虑在内。
DNS可以提供准确的湍流统计数据,但由于计算量巨大,限制了其在工程领域的应用。
二、雷诺平均纳维-斯托克斯(RANS)方法雷诺平均纳维-斯托克斯(RANS)方法是湍流模拟中最常用的方法之一。
该方法基于统计平均,将湍流视为时间均匀的平均流场。
RANS方法通过引入湍流模型来描述湍流的效应,并求解平均速度和湍流应力的方程。
虽然RANS方法计算相对快速,但由于使用了湍流模型,其预测精度受到模型误差的限制。
三、大涡模拟(LES)方法大涡模拟(LES)方法是介于DNS和RANS之间的一种方法。
该方法通过数值滤波将湍流中的大尺度结构进行直接模拟,而将小尺度结构根据模型进行参数化或直接忽略。
LES方法可以提供较高的模拟精度,并在一定程度上保留了湍流的具体特征。
然而,LES方法的计算成本较高,对网格分辨率的要求也很高。
四、湍流模型比较与评估为了评估湍流模拟方法的准确性和适用性,通常需要进行模型比较和验证。
湍流模型的性能评价通常通过与实验数据或更精确的模拟方法进行对比来完成。
1. 实验验证法:实验验证法是评估湍流模拟方法的常用手段之一。
通过与实验数据进行对比,可以直观地了解模拟结果的准确性。
这样的比较涉及到湍流统计量、湍流能谱、湍流结构等方面的对比。
然而,受限于实验条件和设备,实验数据的获取可能受到局限,也可能存在误差。
大气边界层中湍流运动的模拟与分析

大气边界层中湍流运动的模拟与分析大气边界层中的湍流运动对天气预报、空气质量评估以及风电场的建设等领域具有重要的影响。
因此,对大气边界层中的湍流运动进行模拟与分析,能够为解决相关问题提供有效的支持和参考。
本文将介绍湍流运动的模拟方法以及相关分析技术。
一、湍流模拟方法湍流模拟是通过数值方法对大气边界层中的湍流运动进行数值模拟,从而获取湍流场的详细信息。
目前常用的湍流模拟方法包括直接数值模拟(DNS)、大涡模拟(LES)和雷诺平均湍流模拟(RANS)等。
1. 直接数值模拟(DNS)直接数值模拟是一种以最基本的方程组为基础,对大气边界层中湍流运动进行精确模拟的方法。
它通过离散化时间和空间,使用计算机求解Navier-Stokes方程组,得到湍流场的精确解。
但直接数值模拟的计算量非常大,通常仅适用于小尺度或小时间尺度的模拟。
2. 大涡模拟(LES)大涡模拟是一种介于直接数值模拟和雷诺平均湍流模拟之间的方法。
它通过将流场分解为一个大尺度的结构和一个小尺度的湍动结构,只对小尺度湍动进行模拟,通过模拟大尺度结构来减小计算量。
大涡模拟在模拟大气边界层湍流运动方面具有一定的优势。
3. 雷诺平均湍流模拟(RANS)雷诺平均湍流模拟是一种通过对时间和空间进行平均,将湍流场表示为平均量和脉动量的和的方法。
它通过求解雷诺平均Navier-Stokes方程和湍流能量方程,得到湍流场的平均解。
雷诺平均湍流模拟在计算上相对简单,适用于大尺度湍流的模拟。
二、湍流分析技术湍流模拟得到的湍流场数据需要进行进一步的分析才能得到有用的信息。
下面介绍几种常用的湍流分析技术。
1. 自相关函数自相关函数是一种分析湍流场中各点相关性的方法。
它可以通过计算不同点之间的相关性来获取湍流运动的相关长度。
自相关函数可以用于描述湍流场的时空结构。
2. 能谱分析能谱分析是一种通过计算湍流场不同频率分量的能量来了解湍流场特性的方法。
它可以用于表征湍流场的能量分布情况和主导长度尺度。
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)。
直接数值模拟可以获得湍流场的精确信息,是研究湍流机理的有效手段,但现有的计算资源往往难以满足对高雷诺数流动模拟的需要,从而限制了它的应用范围。
雷诺平均方法可以计算高雷诺数的复杂流动,但给出的是平均运动结果,不能反映流场紊动的细节信息。
大涡模拟基于湍动能传输机制,直接计算大尺度涡的运动,小尺度涡运动对大尺度涡的影响则通过建立模型体现出来,既可以得到较雷诺平均方法更多的诸如大尺度涡结构和性质等的动态信息,又比直接数值模拟节省计算量,从而得到了越来越广泛的发展和应用。
2 直接数值模拟(DNS)湍流直接数值模拟(DNS)就是不用任何湍流模型,直接求解完整的三维非定常的N - S 方程组,计算包括脉动在内的湍流所有瞬时运动量在三维流场中的时间演变。
湍流直接数值模拟

Kolmogorov耗散尺度
动速度均方根值),将以上关系代入式,可得:
1/4 3
ห้องสมุดไป่ตู้
,而 ~ u '3 / l ( u ' 是脉
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率
那么三维总网格数N:
=104
这是一个天文数字的估计 ,假设 Rel =104,就要求网格数为109, 考虑到计算的流动变量数,需要约1010字长的计算机内存。直接数值 模拟实际工程湍流运动时,对网格分辨率的要求更高。 应当指出选定最小的网格长度还和数值方法有关。谱方法的数 值精度最高,差分法的精度和差分格式有关。
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率 6.2.2 湍流直接数值模拟的时间分辨率 6.2.3 初始条件和边界条件
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
湍流是多尺度的不规则运动,湍流直接数值模拟和层 流运动的数值计算有很大区别:
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.3 初始条件和边界条件
边界条件的提法: (1) 固壁采用无滑移条件
(2)周期条件
(3)渐近条件→续 (4)进口条件→续
(5)出口条件→续
(6)可压缩湍流的附加边界条件
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.3 初始条件和边界条件
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.3 初始条件和边界条件
湍流的模拟和建模

湍流的模拟和建模湍流是自然界中普遍存在的现象,其涵盖的规模从大气层中的云团到船舶和管道中的流体,十分广泛。
湍流现象表现为流体的不规则而混乱的流动,其中的旋涡和涡旋不断形成和消失。
湍流的复杂性和不可预测性使其对于物理学家和工程师来说是非常具有挑战性的问题。
然而,通过数字模拟和建模,我们可以更好地理解和控制湍流现象,进而提高生产和人类生活的质量。
湍流的模拟和建模一直是流体力学领域的研究热点,旨在通过计算机模拟来预测复杂流动中的物理性质。
对于湍流的模拟,目前主要有两类方法:直接数值模拟(Direct Numerical Simulation, DNS)和大涡模拟(Large-Eddy Simulation, LES)。
其中DNS方法对于湍流的描述最为详细,可以剖析流场中的每一处涡旋,但计算成本极高,通常只适用于小规模的问题。
LES方法通过简化较小尺度的湍流结构来减少计算量,虽然无法完全描述每个涡旋,但是在较大的尺度下仍能准确预测湍流的行为。
湍流现象的建模通常可以基于Navier-Stokes方程进行,这是一组描述流体本质的偏微分方程。
针对这些方程的求解方法和算法不断更新和优化,使得模拟计算变得更加高效和准确。
其中著名的流体力学软件包,包括ANSYS Fluent、OpenFOAM等已经成为工业和研究界广泛应用的工具。
当然,与模拟和建模相伴的,是精度和计算成本之间的取舍。
对于湍流现象的模拟通常需要对涡旋的尺度、湍流能量转化等参数进行详细定量的计算,因此准确度成为了模拟中一个重要的考量因素。
在确定准确度之余,如何减少计算成本也是一个必须解决的问题。
因此,研究人员通常采用增加计算资源的方式,如改进集群计算机和高性能计算机的配置来提升计算速度,并利用一些优化算法和计算技巧来控制误差和减少计算成本。
在湍流模拟和建模方面,模型验证也是一个很重要的步骤,这也是模拟不能完全取代实验的原因之一。
验证过程通常会与实验数据进行比对,用实验数据的帮助来验证模型的准确性。
湍流模拟的数值方法介绍

湍流模拟的数值方法介绍湍流流动是自然界常见的流动现象,是一种高度非线性的复杂流动,但人们已经能够通过某些数值方法对湍流进行模拟,取得与实际比较吻合的结果。
对于湍流运动,已经采用的数值计算方法主要可以分为三类:直接数值模拟、大涡模拟和雷诺时均方程法。
1.直接数值模拟(Direct Numerical Simulation,简称DNS)方法就是直接用瞬时的N-S方程对湍流进行计算。
DNS的最大好处是无需对湍流流动作任何简化或近似,理论上可以得到相对准确的计算结果。
DNS对内存空间及计算速度的要求非常高,目前还无法用于真正意义上的工程计算,但大量的探索性工作正在进行之中。
2. 大涡模拟法(large eddy simulation, 简称LES)为了模拟湍流流动,一方面要求计算区域的尺寸应大到足以包含湍流运动中出现的最大的涡,另一方面要求计算网格的尺度应小到足以分辨最小涡的运动。
然而,就目前的计算机能力来讲,能够采用的计算网格的最小尺度仍比最小涡的尺度大许多。
因此,目前只能放弃对全尺度范围上涡的运动的模拟,而只将比网格尺度大的湍流运动通过N-S方程直接计算出来,对于小尺度的涡对大尺度运动的影响则通过建立模型来模拟,从而形成目前的大涡模拟法。
LES方法的基本思想可以概括为:用瞬时的N-S方程直接模拟湍流中的大尺度涡,不直接模拟小尺度涡,而小涡对大涡的影响通过近似的模型来考虑。
总体而言,LES方法对计算机内存及CPU速度的要求仍比较高,但低于DNS方法。
3.雷诺平均法(RANS:Reynolds-averaged Navier-Stokes)虽然N-S方程可以用于描述湍流,但N-S方程的非线性使得用解析的方法精确描写三维时间相关的全部细节极端困难,即使能真正得到这些细节,对于解决实际问题也没有太大的意义。
这是因为,从工程应用的观点上看,重要的是湍流所引起的平均流场的变化,是整体的效果。
雷诺平均法(Reynolds-averaged Navier-Stokes,简称RANS)是将非稳态的N-S控制方程组作时间平均运算,湍流的各种瞬时量被表示成时均值和脉动值之和,在所得的时均方程中会出现脉动值的乘积的时均值这一类新未知量,从而使方程组不封闭。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
混淆误差产生原因: 当用伪谱方法计算时,谱空分量用,它等于:
利用三角级数公式
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
在上式中,p、q、k的取值范围都是{一N/2,N/2-1},所 以取值不等于零的波数组合情况只有两种:p+q-k=0,或 p+q-k= N,于是有
Kolmogorov耗散尺度
动速度均方根值),将以上关系代入式,可得:
1/4 3
,而 ~ u '3 / l ( u ' 是脉
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率
那么三维总网格数N:
=104
这是一个天文数字的估计 ,假设 Rel =104,就要求网格数为109, 考虑到计算的流动变量数,需要约1010字长的计算机内存。直接数值 模拟实际工程湍流运动时,对网格分辨率的要求更高。 应当指出选定最小的网格长度还和数值方法有关。谱方法的数 值精度最高,差分法的精度和差分格式有关。
边界条件的提法: (3)渐近条件 对于湍流边界层或其他薄湍流切变层,在远离薄层和物面的渐近区 域,速度场趋近于无旋的均匀场,因此对于不可压缩流体可以采用如下:
刚盖假定:数值方法只能计算有限域内的流动,渐近条件只能采用 近似形式,一种方法是在离开薄层或物体横向一定距离的平面上设 置“虚拟边界”,在虚拟边界y=H满足以下条件,称为刚盖假定:
式中,第一项是乘积w=uv在谱空间中的投影,第二项是 在伪谱运算中产生的误差,故称混淆误差。
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
混淆误差的消除--3/2规则 将函数在傅里叶空间中展开的系数进行延拓到原来的3/2倍, 延拓后函数的傅里叶系数有M=3N/2项。延拓按以下的规则赋 值:
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.3 初始条件和边界条件
边界条件的提法: (3)渐近条件 另一种更好的方法是先做一个指数变换,将无限域变到有限域, 例如,令:
然后,在有限域里数值求解Navier-Stokes方程。如果y=0是 固壁,则在指数变换时,在y=0附近自动加密网格,而在 方向 则是均匀网格。在(x,y,z)坐标系里,原渐近边界条件可写作:
(1)加权余量法
L u f (u)
L表示微分算子,f (u) 是已知函数。将未知函数 u 用一组完备的 线性独立函数族 展开: k
k 0,1
当展开式只取有限项时,上式是原函数 u 的近似。把uN代入原 来的微分方程,将产生误差,并称之为残差或余量,用RN表示:
第6章 湍流直接数值模拟
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率
以均匀各向同性湍流为例,假定各向同性湍流的含能尺度或积分 尺度为 l ,Kolmogorov耗散尺度等于 。为了足够准确地计算湍 流的大尺度运动,立方体的长度L必须大于含能尺度 l ,另一方 面,为了保证准确模拟湍流小尺度运动,网格长度 必须小于耗 散尺度 。因此,一维网格数至少应满足以下不等式:
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.3 平面槽道湍流的直接数值模拟
平面槽道湍流是典型的有固壁的简单湍流运动,如图6-8所示。流体 在平行平板之间流动,假定平板在流向和展向都是无限长。平均定常的 槽道湍流的流向平均压强梯度是常数,展向的平均压强梯度等于零。计 算过程中保持平均流量不变,因此以平均速度为特征长度的流动雷诺数 也保持常数。槽道宽度等于2H,计算域的长度和宽度根据流动雷诺数确 定,原则是计算域应当包含足够多的近壁结构。
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
再根据-M/2≤p≤M/2-1和-M/2≤q≤M/2-1(M=3N/2)的取值限制,满足 上式的p和q只有以下四种可能取值范围:
在以上p,q取值范围里,按照延拓公式有 U p Vq 0 ,因此按3/2延 拓规则,式(6.19)的第二项等于零,就是说,式(6.19)的最后结果 是:
第6章 湍流直接数值模拟
6.1 湍流数值模拟的方法
三者对比与联系: (1)直接数值模拟要求模拟所有尺度的湍流脉动;雷 诺平均方法网格尺度应当大于脉动的积分尺度;大涡数 值模拟的网格分辨率介于DNS和RANS之间。
(2)三种湍流数值模拟方法给出的信息量有很大差别。
(3)必须根据需要来选择数值模拟的方法。
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.3 平面槽道湍流的直接数值模拟
控制方程:
边界条件:槽道的上下壁面应是无滑移条件
结合实际可将速度和脉动压强展开如下:
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.3 平面槽道湍流的直接数值模拟
网格坐标为:
理论上,直接数值模拟要求网格最小长度达到Kolmogorov尺 度的量级,因此在给定算例的雷诺数后,应当估计需要的网格数, 以满足分辨率的要求。 首先,确定槽道湍流的雷诺数,Um是槽道截面的平均流速,H是 槽道的半宽度。然后估算Kolmogorov的耗散尺度和时间尺度
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
再选择另一组完备的线性独立函数族 k 作为权函数.要求余量 的加权积分等于零:
如果微分算子L是线性的,则最后求解的是线性代数方程组; 如果微分算子L是非线性的,则最后求解的是非线性代数方程组。 求出代数方程的解,就得到微分方程的近似解。 ①伽辽金方法(Galerkin法) 加权余量法三种形式 ②Tau方法 ③配置点法
以延拓的谱系数做伪谱运算,得到类似的乘积公式
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
上面式中p,q,k的取值范围都是{-M/2,M/2-1},所以,当时,由 延拓公式(6.18a)~(6.18d)可得式(6.19)的第一项,该项恰好是准确的 卷积,即:
由于M=3N/2,式(6.18)中的第二项中p+q=k 3N/2,当 N / 2 k N / 2 1 时,p+q的取值范围是:
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率 6.2.2 湍流直接数值模拟的时间分辨率 6.2.3 初始条件和边界条件
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
湍流是多尺度的不规则运动,湍流直接数值模拟和层 流运动的数值计算有很大区别:
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理 谱方法的优点是精度高,计算速度快 。但是适应复杂边界的 试探函数十分难找,所以对于复杂边界的湍流问题,特别是对于 流场中存在间断的情况,只能采用差分离散方法。
(2)伪谱法和混淆误差 混淆误差概念: 对于非线性方程,有限项谱展开的非线性项会产生附加的 误差,这种误差在谱方法中称做混淆误差(aliasingerror)。
a. 近似非定常出口条件:
(5)出口条 件 b. 嵌边区(或强粘性区)方法
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理 6.3.2 格栅湍流的直接数值模拟 6.3.3 平面槽道湍流的直接数值模拟
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理 谱方法是一种加权余量的数值计算方法,它把微分方程离散化 为代数方程组,其基本原理如下。设有微分方程:
数值模拟结果:各向同性湍流中 的涡结构算例
强迫各向同性湍流直接数值模 拟的结果表明,在各向同性湍流中 强涡是以细长涡管形式出现,如图 6-5所示。涡管的直径是Kolmogorov 尺度,涡管的平均长度是积分尺度, 涡管的强度(环绕涡管的环量)随雷诺 数的增加而增大。图中管状结构是 涡量的等值面,虽然它们并非真正 的涡管,但由于强涡量集中在很细 的管状结构中,可以推断它们接近 于当地的涡管。
第六章 湍流直接数值模拟
第6章 湍流直接数值模拟
6.1 湍流数值模拟的方法 6.2 湍流数值模拟的基本原理 6.3 湍流直接数值模拟的谱方法 6.4 湍流直接数值模拟的差分法
第6章 湍流直接数值模拟
6.1 湍流数值模拟的方法
直接数值模拟(DNS): 从完全精确的流动控制方程出发,对所有尺度的湍流运动进行 数值模拟,这种最精细的数值模拟称为直接数值模拟。 雷诺平均数值模拟(RANS): 从雷诺平均方程出发,在这一层次上的数值模拟称为雷诺平均 数值模拟。 大涡数值模拟(LES): 介于DNS和RANS之间的数值模拟方法称为大涡数值模拟。思想 是:大尺度脉动(或大尺度湍涡)用数值模拟方法计算,只将小尺度 脉动对大尺度运动的作用做模型假设。LES的理论依据是小尺度脉 动有局部平衡的性质,很可能存在某种局部普适的统计规律,如局 部各向同性或局部相似性等。
第一,由于湍流脉动具有宽带的波数谱和频谱,因此湍流直接 数值模拟要求有很高的时间和空间分辨率。
第二,为了求得湍流统计特性,需要足够多的样本流动;如果 湍流是时间平稳态,就要有足够长的时间序列,通常在充分发展的 湍流中,需要105以上的时间积分步。 由于这些特殊要求,需要有内存大、速度快的计算机才能实现 湍流直接数值模 拟。
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理