湍流的数值模拟方法进展
湍流研究的现状和进展
湍流研究的现状和进展
湍流研究是一个广泛的领域,其中有许多不同的话题,需要不同的测量、数值建模和理论研究。
近年来,由于现代计算技术的进步,以及先进的测量设备的出现,湍流研究的发展取得了巨大的进展。
首先,科学家们开发出更为精细的流场数据,提高了对于湍流流场中流动结构的认识,例如,现在可以观测到湍流中的微小力矩和能量分布,并通过数值模拟得出更为详细的解释。
其次,科学家提出了一系列新的数值模型,更深入地探讨湍流的细节,以更合理地描述湍流的结构和性质。
结合传统的经验方法,这些数值模型使得对湍流的模拟更加准确,并且可以用于多种应用领域,例如自由混沌运动模拟、叶片流动计算、空气动力学分析、重力和磁场力分析等。
此外,近年来,学者们开发出了一系列新的控制策略,以改善湍流的流动性能。
这些策略的准备方法涉及精细的数值模拟,基于大量的实验测量数据,有效地改进湍流流场的特性,从而提高流体动力学的效率。
总的来说,近年来,湍流研究的发展取得了巨大的进展。
精细的数据和模型,以及新的控制策略,有助于更好地认识和控制湍流流场,进一步提高流体动力学性能。
湍流对水力旋流器内颗粒运动行为影响的数值模拟研究
湍流对水力旋流器内颗粒运动行为影响的数值模拟研究章节一:引言- 研究背景和意义- 国内外研究现状和进展- 研究目的和内容章节二:数值模拟方法- 水力旋流器模型建立- 流场计算方法- 颗粒运动模拟方法章节三:湍流对颗粒运动的影响- 颗粒运动的基本规律- 湍流对颗粒运动的影响机理- 湍流对颗粒运动的数值模拟方法章节四:数值模拟结果与分析- 颗粒运动规律的数值模拟结果- 湍流对颗粒运动的影响分析- 影响因素对颗粒运动的影响分析章节五:结论与展望- 研究结论和贡献- 研究不足和展望- 实践应用前景和意义第一章节:引言水力旋流器是一种用于液-固分离的设备,可在多个行业中应用,如矿业、环保、化工等。
旋流器通过在旋流腔内制造高速旋转的涡流,使颗粒物沉积到旋流器底部,从而实现液-固分离。
水力旋流器具有结构简单、操作方便、效率高等优点。
然而,水力旋流器的颗粒物分离效果和流量稳定性可以受到湍流的影响。
湍流是一种流体运动状态,具有涡旋混合和能量损失等特性。
水力旋流器在使用过程中,液体将以高速旋转的方式进入旋流腔内,此时液体会形成湍流流动,并影响颗粒物的沉降过程。
近年来,数值模拟方法的发展使得对于液-固分离设备内部颗粒物分布状况的研究有了更深入的了解,也为对于水力旋流器内颗粒运动的研究提供了可靠的工具。
本论文旨在通过数值模拟方法,研究湍流对水力旋流器内颗粒运动行为的影响。
首先,本章节将介绍研究背景和意义。
随后,将从国内外研究现状和进展以及研究目的和内容两方面阐述本研究的意义和贡献。
一、研究背景和意义水力旋流器是一种有效的液-固分离设备,在工业生产过程中得到了广泛的应用。
然而,水力旋流器内部液流具有高速旋转的特点,这种特殊的流动状态会对颗粒物的沉降产生影响,导致颗粒物分离效果不佳,影响设备的稳定运行。
因此,对于水力旋流器内部流场结构和颗粒物分布规律进行深入研究,有助于进一步提高旋流器的分离效果和流量稳定性。
二、国内外研究现状和进展已有研究表明,水力旋流器内的颗粒物分离效果与湍流的存在程度有关。
大气湍流模拟与方法研究
大气湍流模拟与方法研究大气湍流是指在自然界中,由于空气分子的热运动和流体不均匀性造成的气流乱流现象。
它对于气象学、环境科学等领域具有重要意义。
为了更好地理解和预测大气湍流,科研人员们进行了大量的模拟与方法研究。
本文将介绍大气湍流模拟的方法和相关研究进展。
一、大气湍流模拟方法1. 直接数值模拟(DNS)直接数值模拟是指通过求解流体动力学方程,对湍流进行精确的数值模拟。
这种方法能够提供精确的湍流数据,但由于计算量巨大,目前只适用于小尺度的湍流问题。
2. 大涡模拟(LES)大涡模拟是指模拟并求解大尺度涡旋,而忽略小尺度湍流的准确表示。
它通过将湍流分解成大涡和小涡来减少计算量,更适用于中等和大尺度的湍流研究。
3. 湍流统计模拟(TST)湍流统计模拟是一种基于概率和统计的模拟方法,通过对湍流的统计特性进行建模,推导出湍流的各种物理参量。
虽然它无法提供湍流的详细结构信息,但能够在计算成本较低的情况下估计湍流的平均性质。
二、大气湍流模拟方法的应用1. 大气环流模拟大气环流是指大尺度的大气运动模式,是全球气候变化和天气预报的重要基础。
通过模拟大气环流中的湍流现象,可以更准确地预测天气变化和气候变化趋势。
2. 污染物扩散模拟污染物扩散是大气科学中的重要研究内容。
将湍流模拟方法应用于污染物扩散模拟中,可以帮助科研人员分析城市污染物的来源、传输路径和浓度分布,为环境保护和污染治理提供科学依据。
3. 风能资源评估利用风能发电是一种清洁、可再生的能源利用方式。
通过模拟大气湍流,可以评估风能资源的分布和利用潜力,为风电场选址和设计提供技术支持。
三、大气湍流模拟方法的挑战与展望1. 精度提升当前的大气湍流模拟方法仍然存在精度不高的问题,特别是对于小尺度湍流的模拟。
因此,需要进一步改进模拟算法和数值计算技术,提高模拟结果的准确性。
2. 计算成本降低目前的大气湍流模拟方法需要耗费大量的计算资源和时间。
在提高精度的前提下,需要寻求更高效的计算方法,降低计算成本。
高温高速湍流流动理论解释及数值模拟方法
高温高速湍流流动理论解释及数值模拟方法一、引言高温高速湍流流动是指介质在高温和高速条件下出现湍流现象的流动过程。
这种流动现象在航空航天、能源转化、燃烧和化学反应等领域中具有重要的应用,但也面临着诸多挑战。
理解高温高速湍流流动的机理以及采用合适的数值模拟方法来模拟这些流动过程对于相关领域的研究和工程应用至关重要。
二、高温高速湍流流动理论解释高温高速湍流流动的理论解释涉及到湍流的起源、湍流的统计性质以及湍流的能量传递过程。
目前存在多种湍流理论模型,如光照湍流理论、能量谱理论和动能传输理论等。
湍流的起源主要是由于流动中的速度梯度引发的不稳定性产生的。
在高温高速条件下,速度梯度的变化更加剧烈,从而使得流动更易产生湍流现象。
湍流的统计性质包括宏观观测量和微观观测量。
宏观观测量可以通过湍流动能的分布和统计参数来描述,而微观观测量则需要了解湍流中的小尺度结构和统计信息。
湍流的能量传递过程是湍流能量从大尺度向小尺度传递的过程。
这一过程可以通过湍流的能谱来描述,能谱将湍流能量和空间尺度联系起来。
对于高温高速湍流流动,湍流能量的传递过程往往更为复杂,需要采用适当的数值模拟方法来加以分析。
三、数值模拟方法在研究高温高速湍流流动时,数值模拟方法是一种非常有效的手段。
常用的数值模拟方法包括直接数值模拟(DNS)、大涡模拟(LES)和雷诺平均纳维-斯托克斯方程(RANS)等。
直接数值模拟(DNS)是一种精确模拟湍流流动过程的方法,可以完全解析湍流中的小尺度结构。
然而,由于高温高速湍流流动的复杂性,直接数值模拟的计算量较大,需要高性能计算设备的支持。
因此,直接数值模拟在实际工程应用中的局限性较大。
大涡模拟(LES)是一种介于直接数值模拟(DNS)和雷诺平均纳维-斯托克斯方程(RANS)之间的方法。
在LES中,较大尺度的湍流结构直接求解,而较小尺度的结构则通过子网格模型来模拟。
这种方法可以在更小的计算量下获得湍流流动的一些重要特征。
湍流燃烧数值模拟的研究与进展
湍流燃烧数值模拟的研究与进展湍流燃烧是指在燃烧过程中,燃料与氧化剂在湍流的条件下相遇和反应。
湍流燃烧数值模拟是一种通过计算机模拟湍流燃烧过程的方法,可以提供燃烧器内部的流场和温度分布等信息,对于燃烧器的设计和优化具有重要的意义。
本文将对湍流燃烧数值模拟的研究与进展进行探讨。
首先,湍流模型的选择是湍流燃烧数值模拟的一个关键问题。
湍流现象十分复杂,需要选择适当的湍流模型来模拟湍流流动。
常用的湍流模型有雷诺平均应力模型(RANS)和大涡模拟(LES)。
RANS是一种将湍流场分为均匀部分和涡旋部分的统计方法,适用于模拟湍流较为稳定的情况;而LES则能模拟较为精细的湍流结构,但计算量较大。
根据具体问题的复杂程度和计算资源的限制,选择适当的湍流模型具有重要意义。
其次,化学反应模型的建立是湍流燃烧数值模拟的另一个关键问题。
燃烧过程中涉及到多种化学反应,需要建立合适的化学反应模型来描述燃烧反应。
常见的化学反应模型有简化化学反应模型和详细化学反应模型。
简化化学反应模型基于简化的反应机理,计算速度较快;而详细化学反应模型则基于包含大量反应步骤的反应机理,计算速度较慢但结果更精确。
根据具体问题的要求和计算资源的限制,选择适合的化学反应模型具有重要意义。
此外,边界条件的设定也是湍流燃烧数值模拟的一个关键问题。
边界条件的合理设定可以保证计算结果的准确性。
常用的边界条件有Inflow Boundary Condition、Outflow Boundary Condition、Wall Boundary Condition等。
对于湍流燃烧数值模拟,还需要考虑湍流场的边界条件,例如由湍流脉动引起的湍流输运方程中的涡粘性项的边界条件等。
最后,计算方法的选择也对湍流燃烧数值模拟的结果和计算速度有着重要的影响。
常用的计算方法有有限差分法(FDM)、有限元法(FEM)和有限体积法(FVM)等。
这些方法在计算精度和计算速度方面各有优势,需要根据具体问题的要求选择适当的方法。
湍流的几种数值模拟方法
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
数值模拟技术的最新进展
数值模拟技术的最新进展近年来,随着计算机技术的飞速发展,数值模拟技术也越来越得到人们的关注。
数值模拟技术是指利用数值方法对物理过程进行仿真和计算的技术,它在许多领域中都有应用,如机械工程、航空航天、地质勘探等等。
本文主要介绍数值模拟技术在目前的最新进展。
1. 基于深度学习的数值模拟深度学习技术是近年来飞速发展的人工智能技术之一,它通过模仿人类神经网络的结构和学习方式,从大量数据中学习并识别模式。
在数值模拟方面,深度学习技术可以用于建立高效、准确的模型。
例如,科学家可以将深度学习技术用于流体力学模拟中,这使得模拟能够更快速地进行,并且可以更加准确地预测流体力学现象,例如湍流、流体分离等等。
2. 并行计算技术随着计算机硬件技术的迅速发展,现在的计算机系统已经能够支持大规模并行计算。
并行计算技术是指将计算任务分成多个子任务来同时执行,以提高计算速度。
在数值模拟中,大规模并行计算技术可以极大地提高计算效率,从而使得更加复杂的模拟成为可能。
例如,人们可以用并行计算技术来模拟地震波传播过程,这将有助于更好地理解地震现象的本质。
3. 多物理场耦合模拟多物理场耦合模拟是指通过数值算法将不同物理学领域(如流体力学、电磁学、结构力学等)中的方程同时求解,以模拟多物理场耦合的物理现象。
例如,在航空航天工程中,飞机的设计要求同时考虑结构力学、燃烧、流体力学等多个因素,这时就需要用到多物理场耦合模拟技术。
目前,多物理场耦合模拟技术已经成为数值模拟领域中的一个重要方向。
4. 高性能计算技术高性能计算技术是指计算机系统通过优化计算资源的配置,以提高计算效率和性能。
在数值模拟中,高性能计算技术尤其重要,因为模拟过程中需要进行大量的计算,需要在有限的时间内完成计算任务。
例如,在气象预报领域中,高性能计算技术可以帮助气象预报模型更加准确地预测天气情况,提高天气预报的准确度。
5. 高精度数值算法高精度数值算法是指通过提高数值计算方法的精度,以提高模拟效果的技术。
湍流的数值模拟方法进展
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 τ即称为次网格剪切应力张量(亦称为亚格子应力)。
流体流动中的湍流弛豫时间研究
流体流动中的湍流弛豫时间研究引言流体是一种非常常见的物质,在自然界和工业生产中都起着重要的作用。
流体的流动特性对于许多工程和科学问题的研究具有重要意义。
其中,湍流是流体流动中最为复杂和困难的问题之一。
湍流的出现常常伴随着复杂的流动模式,这使得湍流的研究成为科学家们一直以来的挑战。
湍流弛豫时间是湍流研究中的一个重要概念。
它描述了湍流流动中的能量耗散过程,对于湍流现象的理解和控制具有重要的意义。
本文将重点研究流体流动中的湍流弛豫时间,并探讨其对湍流流动行为的影响,为湍流研究提供新的思路和方法。
湍流的基本特性湍流是指流体流动中的一种混乱而不规则的状态。
相比于层流流动,湍流流动表现出波动幅度更大、流速方向更混乱的特点。
湍流的出现往往伴随着能量的耗散和能量的重新分配过程,这使得湍流具有一种自维持的能力。
湍流的基本特性可以通过流体的速度分布和涡旋结构来揭示。
在湍流流动中,流体的速度分布呈现出不规则的形态,存在着大小不同的涡旋结构。
这些涡旋结构的产生与湍流中的流体动量交换密切相关,通过分析涡旋的分布和运动规律,可以研究湍流流动的特性和行为。
湍流的产生有多种原因,包括流体的速度梯度、湍流中的涡旋相互作用等。
在具体的流体流动问题中,湍流的产生往往是复杂的多因素综合作用的结果。
因此,对湍流的研究需要深入理解湍流的基本特性和产生机制。
湍流弛豫时间的概念湍流弛豫时间可以看作是湍流能量耗散的时间尺度。
在湍流流动中,湍流能量以某种形式存在,并随着时间的推移逐渐耗散。
湍流弛豫时间描述了湍流能量耗散的速率和过程,反映了湍流流动的动态特性。
湍流弛豫时间的具体定义与湍流的能量耗散存在密切的关系。
湍流能量耗散可以通过湍流流动的涡旋结构和流体动量交换来理解。
涡旋结构的生成和消亡是湍流能量耗散的关键过程之一,而流体动量交换则是湍流能量转移和重新分配的重要途径。
湍流弛豫时间可以通过实验和数值模拟等方法来进行研究。
实验方法包括测量湍流流动的速度分布和涡旋结构,通过分析涡旋的运动规律和能量耗散过程来确定湍流弛豫时间。
湍流燃烧及其数值模拟
湍流燃烧及其数值模拟研究1. 湍流燃烧1.1湍流燃烧基本概念当流动雷诺数数较小时,由于流体粘性的作用,流体呈层流流态。
当流动的特征雷诺数超过相应的临界值,流动从层流转捩到湍流。
湍流燃烧是指湍流流动中可燃气的燃烧,在能源、动力、航空和航天等工程领域,经常遇到的实际燃烧过程几乎全部都是湍流燃烧过程。
湍流燃烧实质是湍流,化学反应和传热传质等过程相耦合的结果。
湍流对燃烧的影响与湍流强度和湍流涡旋尺度有关。
小尺度湍流通过湍流扩散使火焰区内的输运效应增加,从而使化学反应速率增加。
但气流脉动不会火焰面产生皱褶,只能把火焰变成波纹状。
大尺度湍流对火焰内部结构没有影响,但使火焰阵面出现皱褶,增加其燃烧面积,造成火焰表现传播速度增加。
当湍流强度及湍流尺度均较大时,火焰前沿不再连续而分裂成四分五裂。
燃烧对湍流的影响主要表现在燃烧释放的热流流团膨胀,影响气体的密度和运动速度,从而影响当地的涡旋,湍流强度和湍流结构。
1.2湍流燃烧分类湍流燃烧按其燃料和氧化剂的初始混合状态可以分类为:湍流非预混燃烧、预混燃烧和部分预混燃烧。
在湍流非预混燃烧燃料和氧化剂事先是分离的,燃料和氧化剂一边混合一边燃烧,燃烧速率主要受湍流混合过程控制,而在湍流预混燃烧中,燃料和氧化剂在进入核心燃烧区以前已经充分混合,化学反应的速率由火焰前缘从炽热的燃烧区向冷态无反应区的传播所控制。
上面两种燃烧方式是湍流燃烧的两个极限情形,很多情况下两种燃烧模式是并存的,称为部分预混燃烧。
部分预混燃烧可出现在下列情形中叫:(1)在一个完全以非预混燃烧为配置的燃烧装置发牛了局部熄火;(2)当预混火焰前缘穿过非均匀的混气时;(3)射流非预混火焰发生抬举,其根部是一个典型的部分预混火焰。
这三种部分预混燃烧情形涉及了经常受到关注的燃烧研究话题如局部熄火、火焰稳定等,它们对研究湍流燃烧过程的机理有很大意义。
在湍流燃烧中,湍流流动过程和化学反应过程有强烈的相互关联和相互影响.湍流通过强化混合而影响着时平均化学反应速率,同时化学反应放热过程又影响着湍流,如何定量地来描述和确定这种相互作用是湍流燃烧研究的一个重要内容.湍流是非常复杂的,它包括湍流问题,湍流与燃烧的相互作用,流动参数与化学动力参数之间的耦合机理等问题。
流体流动中的湍流动态行为研究
流体流动中的湍流动态行为研究引言湍流是流体力学中的一个重要研究课题,广泛应用于工程、气象、地质等领域。
湍流流动具有复杂的动态行为,其特征包括涡旋结构、乱流能量耗散和湍流尺度等。
理解和控制湍流的动态行为对于提高流动系统的效率、降低能量损耗具有重要意义。
本文将介绍湍流的基本概念、湍流动态行为的研究方法以及最新的研究进展。
湍流的基本概念湍流是流体流动中的一种不规则流动形式,其速度、压力等参量在时间和空间上存在剧烈的波动和涨落。
湍流的发生是由于流体流动速度超过一定的临界速度,引起流动层的不稳定性而产生的。
湍流的形成过程可以分为湍流发展阶段和湍流稳定阶段。
湍流的基本特征包括: - 涡旋结构:湍流中存在着一系列的涡旋,其大小和形态各异,相互之间的交互作用引起了湍流的混沌行为。
- 乱流能量耗散:湍流流动中存在着能量的不断转换和耗散过程,这种能量的转换从而维持了湍流的持续存在。
- 湍流尺度:湍流具有多个尺度,从小尺度到大尺度,相应地存在着不同的湍流结构和运动规律。
湍流动态行为的研究方法湍流动态行为的研究方法主要包括实验观测、数值模拟和理论分析。
实验观测实验观测是研究湍流动态行为的直接方法,可以通过流场的测量和数据处理获得湍流的相关信息。
常用的实验观测方法包括: - 激光多普勒测速技术:通过测量流体中散射的激光来确定流体的速度和速度梯度等信息,得到湍流的速度场和相关统计量。
- 热线方法:利用热线感应的电阻变化来测量流体的速度,可以得到湍流的速度场和能量谱等信息。
- 高速摄影技术:通过高速摄影记录湍流流动的图像序列,可以研究湍流的结构演化和动态行为。
数值模拟数值模拟是研究湍流动态行为的常用方法,通过数值方法求解湍流流动的方程组来模拟湍流的发展和演化。
常用的数值模拟方法包括:- 直接数值模拟(DNS):通过求解流体动力学方程来模拟湍流流动的每个尺度和涡旋结构,适用于小尺度湍流的研究。
- 大涡模拟(LES):通过分解流动的涡旋结构,只模拟大尺度的涡旋,通过子网模型来参数化小尺度湍流的作用。
湍流普朗特数 rans-定义说明解析
湍流普朗特数rans-概述说明以及解释1.引言1.1 概述湍流是一种流动状态,其中流体的运动变得混乱且不规则。
湍流普朗特数RANS(Reynolds Averaged Navier-Stokes)模型是一种常见的湍流模拟方法,用于描述湍流流动的统计平均行为。
在实际工程应用中,湍流模型可以帮助工程师们预测和分析复杂流动情况下的压力变化、速度分布、热传导和物质输运等问题。
为了更好地理解湍流普朗特数RANS模型,首先需要了解湍流的定义和特征。
湍流具有随机性和无规律性,流体中各个点的速度和压力会发生剧烈变化。
相比于层流流动,湍流流动具有更高的能量耗散,表现出不同长度尺度的涡旋结构。
湍流具有非常广泛的应用领域,包括但不限于气象学、航空航天工程、石油化工等。
普朗特数是描述流动性质的重要参数之一,它衡量了动量输送和能量耗散之间的关系。
普朗特数越大,表示动量输送能力越强,能量耗散相对较小。
反之,普朗特数越小,表示动量输送较弱,能量耗散相对较大。
在湍流模拟中,普朗特数在计算动量和能量传输方程中起着重要的作用,与湍流的流动特性密切相关。
RANS模型是基于对湍流进行统计平均的一种模型。
它假设湍流中的速度和压力可以分解为平均分量和涡旋分量,通过解析平均分量方程和湍流涡旋分量方程来计算流动的宏观行为。
RANS模型的应用广泛,可以有效地模拟复杂的湍流流动,例如气流在飞机翼上的流动、液体在流体管道中的传输等。
本文的目的是对湍流普朗特数RANS模型进行深入研究和探讨。
通过对湍流的定义、特征以及普朗特数的概念和意义进行剖析,我们将全面了解湍流模拟方法的原理和应用。
同时,我们还将探讨RANS模型的优缺点以及未来的发展方向。
通过本文的研究,我们可以更好地理解湍流普朗特数RANS模型在工程领域中的应用,并为进一步提高湍流模拟的准确性和可靠性提供参考和指导。
1.2 文章结构本文将按照以下结构进行讨论和分析湍流普朗特数(Reynolds Average Navier-Stokes, RANS)的相关内容:首先,在引言部分对本文的目的和重点进行简要介绍,为读者提供背景信息和整体把握文章的框架。
大气湍流畸变相位屏的数值模拟方法研究_张慧敏
⎪ Θm (θ ) = ⎨ 2(n + 1) sin(mθ )
m≠0
odd term
(10)
⎪ ⎩
(n + 1) Rn0 (r)
m=0
式中 (r,θ)为单位圆域内的极坐标。m,n 分别是多项式角向级次和径向级次。 大气波前的 Zernike 系数向量 a 的协方差矩阵Γa,其元素为 E{apap′}(p,p′=1,…,P),P 为采用的 Zernike
∑ φ(r) =
a ∞
j =1 j
⋅
z j (r)(7) Nhomakorabea其中 zj(r)为各项 Zernike 多项式,aj 为第 j 项 Zernike 多项式系数。 圆域 Zernike 多项式的解析形式不唯一,因为存在幅度因子的选择不同和模式序号的不同排列方式。根
据 Noll[3]的定义,圆域 Zernike 多项式可以表示为径向多项式和角向多项式的乘积
f
(m', n') exp( j2π3− p
( mm' N
+
nn' )) N
(6)
式中 p 为次谐波级数, f (m', n' ) = C ⋅ 3−2 p r0 −5 / 6 ( flx 2 + fly 2 ) , flx = 3− p m' ∆f x , fly = 3− p n' ∆f y
(6)式和(1)式两部分的合并就构成了相位屏的总相位,从而补偿了相位屏的低频部分。
式中 z 为光束传播方向;Cn2(z)代表传播路径上湍流强度的结构常数;κ为三维空间波数,此时,与传播方
向(z 轴方向)垂直的任意薄层切片上的大气相位功率谱为
∫ Fφ (κ r ) = 2πk 2 0.033κ r −11/3
(2021年整理)湍流的数值模拟综述
湍流的数值模拟综述编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(湍流的数值模拟综述)的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为湍流的数值模拟综述的全部内容。
湍流的数值模拟一、引语流体的流动形态分为湍流与层流。
而层流是流体的最简单的一种流动状态。
流体在管内流动时,其质点沿着与管轴平行的方向作平滑直线运动。
此种流动称为层流或滞流,亦有称为直线流动的.流体的流速在管中心处最大,其近壁处最小。
管内流体的平均流速与最大流速之比等于0.5,根据雷诺实验,当雷诺准数引Re<2320时,流体的流动状态为层流.当雷诺数Re〉2320时,流体流动状态开始向湍流态转变,湍流是一种很复杂的流动状态,是流体力学中公认的难题。
自从19世纪末O.Reynolds提出湍流的统计理论以来,已经有一个多世纪了,经过几代科学家的努力,湍流研究取得很大进展,但是仍然不能满足工程应用的需要,以至于经常有悲观的论调侵袭湍流研究。
为什么湍流问题没有圆满地解决会受到如此关注呢?因为湍流是自然界和工程中十分普遍的流动现豫,对于湍流问题的正确认识和模化直接影响到对自然环境的预测和工程的质量.例如,当前影响航天器气动力和气动热预测准确度的主要障碍是缺乏可靠的湍流模型。
和其他一些自然科学的准题不同,解决湍流问题具有迫切性。
湍流运动的最主要特征是不规则性,这是大家公认的。
对于湍流不规则性的深入认识,是一百多年来湍流研究的上要成就之一。
早期的科学家认为,像分子运动一样,湍流是完全不规则运动。
类似于分子运动产生黏性,湍流的耗散可以用涡黏系数来表述。
20世纪初,一些杰出的流体力学家,相继对涡黏系数提出各种流体力学的模型,如Taylor(1921年)的涡模型,Praudtl(1925年)的混合长模型和von Karman (1930年)相似模型等。
流体力学的发展趋势
流体力学的发展趋势21162P21吕鹏2012.3定义流体力学,是研究流体(液体和气体)的力学运动规律及其应用的学科。
主要研究在各种力的作用下,流体本身的状态,以及流体和固体壁面、流体和流体间、流体与其他运动形态之间的相互作用的力学分支。
流体力学是力学的一个重要分支,它主要研究流体本身的静止状态和运动状态,以及流体和固体界壁间有相对运动时的相互作用和流动的规律。
在生活、环保、科学技术及工程中具有重要的应用价值。
重要性上上个世纪在运河河道中发现的孤立波在60年代得到了彻底的解决,既推动了力学和数学的发展,也迅速导致在其它学科如光学、声学中发现类似的现象。
现在孤立波(光学中称孤立子)已成了光通信的基石。
上世纪60年代,为探索为何基于流体力学方程的数值天气预报只能准确到很少几天,通过简化这组方程之后,得到了现在已十分著名的Lorenz方程。
数值计算表明,它的解对初值十分敏感,以致一定时间之后,其值变得几乎完全不可预测的了。
这一发现开辟了混沌研究新领域,奠定了非线性科学的基础。
这一事实还说明,流体力学方程(NS方程)的内涵十分深邃,对它的了解还远不是充分的。
水波中各种波的非线性作用的研究,也丰富了非线性科学的内容。
凡此种种,显示出了本世纪流体力学在科学发展中的作用。
流体力学在工程技术中的作用,更是有目共睹的。
飞机的飞行速度得以超过声速,是空气动力学发展的结果。
人类登月的成功,大型火箭和航天飞机的实现,需要解决成千上万个前所未有的难题,而力学问题往往首当其冲。
为此形成了高超声速气动力学,物理化学流体力学,稀薄气体力学等一系列新的分支学科,并极大地推动了计算科学的发展。
为解决喷气机的噪声问题,提出了流体噪声理论,它完全不同于经典的声学理论。
各种高速、高机动性和高敏捷性的军用飞机和安全、舒适的大型民航机的研制成功,同样需要流体力学提供的新思想和新成果。
70年代兴起的海上采油工业,若没有流体力学的研究成果为依据,设计、建造单台价值超过10亿美元的海上采油平台是不可能的。
流体力学发展概况和趋势
流体力学发展概况和趋势作为物理的一部分,流体力学在很早以前就得到发展。
在19世纪,流体力学沿着两个方面发展,一方面,将流体视为无粘性的,有一大批有名的力学数学家从事理论研究,对数学物理方法和复变函数的发展,起了相当重要的作用;另一方面,由于灌溉、给排水、造船,及各种工业中管道流体输运的需要,使得工程流体力学,特别是水力学得到高度发展。
将二者统一起来的关键是本世纪初边界层理论的提出,其中心思想是在大部分区域,因流体粘性起的作用很小,流体确实可以看成是无粘的。
这样,很多理想流体力学理论就有了应用的地方。
但在邻近物体表面附近的一薄层中,粘性起着重要的作用而不能忽略。
边界层理论则提供了一个将这两个区域结合起来的理论框架。
边界层这样一个现在看来是显而易见的现象,是德国的普朗特在水槽中直接观察到的。
这虽也是很多人可以观察到的,却未引起重视,普朗特的重大贡献就在于他提出了处理这种把两个物理机制不同的区域结合起来的理论方法。
这一理论提出后,在经过约10年的时间,奠定了近代流体力学的基础。
有意思的是在流体力学中发现的这种边界层现象,很快地在别的科学领域得到了响应,因为这里面包含了更广泛和深刻的内容。
由此又大大促进了应用数学的发展,从而形成了现在在很多科学中广泛应用的“渐近匹配法”。
在流体力学中首先发现的现象及为此提出的理论,在一段时间以后被发现在其他学科领域中同样存在和有用,这样的例子并不是唯一的一个。
例如,100年前在水波中观察到的孤立波及其理论到本世纪60年代被发现在声波、光波中同样存在和有用,从而迅速形成了系统的理论。
目前具有重要应用前景的光通讯,正是建立在孤立子(孤立波)理论基础上的。
又如在上个世纪发现的流体从下部加热从而引起对流并能形成有规则图形的现象,以及本世纪20年代发现的两旋转圆筒间所充满的流体在一定条件下能形成有规则的二次流的现象,成了近代在各个学科领域中普遍关注的分岔现象及理论的经典例子。
而且也是最近逐步形成的图形(pattern)动力学的典型例子及实验对象。
湍流的数值模拟
湍流的数值模拟编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(湍流的数值模拟)的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为湍流的数值模拟的全部内容。
2012年秋季学期研究生课程考核湍流的数值模拟一、流体力学概述流体力学是研究流体的力学运动规律及其应用的学科。
主要研究在各种力的作用下,流体本身的状态,以及流体和固体壁面、流体和流体间、流体与其他运动形态之间的相互作用的力学分支。
除水和空气之外,这里的流体还指作为汽轮机工作介质的水蒸气、润滑油、地下石油、含泥沙的江水、血液、超高压作用下的金属和燃烧后产生成分复杂的气体、高温条件下的等离子体等等。
它的主要基础是牛顿运动定律和质量守恒定律,常常还要用到热力学知识,有时还用到宏观电动力学的基本定律、本构方程和高等数学、物理学、化学的基础知识。
气象、水利的研究,船舶、飞行器、叶轮机械和核电站的设计及其运行,可燃气体或炸药的爆炸,汽车制造,以及天体物理的若干问题等等,都广泛地用到流体力学知识.许多现代科学技术所关心的问题既受流体力学的指导,同时也促进了它不断地发展。
二、数值计算在流体力学研究中的应用数值计算是研究流体力学的重要方法。
它是针对流体运动的特点,用数学语言将质量守恒、动量守恒、能量守恒等定律表达出来,从而得到连续性方程、动量方程和能量方程。
此外,还要加上某些联系流动参量的关系式(例如状态方程),或者其他方程。
这些方程合在一起称为流体力学基本方程组。
求出方程组的解后,结合具体流动,解释这些解的物理含义和流动机理。
通常还要将这些理论结果同实验结果进行比较,以确定所得解的准确程度和力学模型的适用范围.从基本概念到基本方程的一系列定量研究,都涉及到很深的数学问题,所以流体力学的发展是以数学的发展为前提.反过来,那些经过了实验和工程实践考验过的流体力学理论,又检验和丰富了数学理论,它所提出的一些未解决的难题,也是进行数学研究、发展数学理论的好课题.按目前数学发展的水平看,有不少题目将是在今后几十年以内难于从纯数学角度完善解决的。
湍流理论发展概述
湍流理论发展概述一、湍流模型的研究背景自然环境和工程装置中的流动常常是湍流流动,模拟任何实际过程首先遇到的就是湍流问题,而湍流问题本身又是流体力学理论上的难题。
对于某些简单的均匀时均流场,如果湍流脉动是各向均匀及各向同性的,可以用经典的统计理论来分析,但实际上的湍流往往是不均匀的,这就给理论分析带来了极大地困难。
这也就引发了对湍流过程进行模拟的想法。
对湍流最根本的模拟方法是在湍流尺度的网格尺寸内求解瞬态的三维N-S 方程的全模拟方法,此时无需引进任何模型。
然而由于计算方法及计算机运算水平的限制,该种方法不易实现。
另一种要求稍低的方法是亚网格尺寸度模拟即大涡模拟(LES),也是由N-S方程出发,其网格尺寸比湍流尺度大,可以模拟湍流发展过程的一些细节,但由于计算量仍然很大,只能模拟一些简单的情况,直接应用于实际的工程问题也存在很多问题[1]。
目前数值模拟主要有三种方法:1.平均N-S方程的求解,2.大涡模拟(LES),3.直接数值模拟(DNS),而模拟的前提是建立合适的湍流模型。
所谓的湍流模型,就是以雷诺平均运动方程与脉动运动方程为基础,依靠理论与经验的结合,引进一系列模型假设,而建立起的一组描写湍流平均量的封闭方程组。
目前常用的湍流模型可根据所采用的微分方程数进行分类为:零方程模型、一方程模型、两方程模型、四方程模型、七方程模型等。
对于简单流动而言,一般随着方程数的增多,精度也越高,计算量也越大、收敛性也越差。
但是,对于复杂的湍流运动,则不一定。
湍流模型可根据微分方程的个数分为零方程模型、一方程模型、二方程模型和多方程模型。
这里所说的微分方程是指除了时均N-S 方程外,还要增加其他方程才能是方程封闭,增加多少个方程,则该模型就被成为多少个模型。
二、基本湍流模型常用的湍流模型有:零方程模型:C-S模型,由Cebeci-Smith给出;B-L模型,由Baldwin-Lomax 给出。
一方程模型:来源由两种,一种从经验和量纲分析出发,针对简单流动逐步发展起来,如Spalart-Allmaras(S-A)模型;另一种由二方程模型简化而来,如Baldwin-Barth(B-B)模型。
(2021年整理)湍流模拟的数值方法介绍
湍流模拟的数值方法介绍编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(湍流模拟的数值方法介绍)的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为湍流模拟的数值方法介绍的全部内容。
湍流模拟的数值方法介绍湍流流动是自然界常见的流动现象,是一种高度非线性的复杂流动,但人们已经能够通过某些数值方法对湍流进行模拟,取得与实际比较吻合的结果。
对于湍流运动,已经采用的数值计算方法主要可以分为三类:直接数值模拟、大涡模拟和雷诺时均方程法。
1。
直接数值模拟(Direct Numerical Simulation,简称DNS)方法就是直接用瞬时的N-S方程对湍流进行计算.DNS的最大好处是无需对湍流流动作任何简化或近似,理论上可以得到相对准确的计算结果.DNS对内存空间及计算速度的要求非常高,目前还无法用于真正意义上的工程计算,但大量的探索性工作正在进行之中.2. 大涡模拟法(large eddy simulation,简称LES)为了模拟湍流流动,一方面要求计算区域的尺寸应大到足以包含湍流运动中出现的最大的涡,另一方面要求计算网格的尺度应小到足以分辨最小涡的运动。
然而,就目前的计算机能力来讲,能够采用的计算网格的最小尺度仍比最小涡的尺度大许多。
因此,目前只能放弃对全尺度范围上涡的运动的模拟,而只将比网格尺度大的湍流运动通过N—S方程直接计算出来,对于小尺度的涡对大尺度运动的影响则通过建立模型来模拟,从而形成目前的大涡模拟法。
LES方法的基本思想可以概括为:用瞬时的N-S方程直接模拟湍流中的大尺度涡,不直接模拟小尺度涡,而小涡对大涡的影响通过近似的模型来考虑。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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)满足 3.3 控制方程将过滤函数作用与N-S 方程的各项,得到过滤后的湍流控制方程组:由于无法同时求解出变量i u 和j i u u ,所以将j i u u 分解成i j i j ij u u u u τ=⋅+,ij τ即称为次网格剪切应力张量(亦称为亚格子应力)。
由此动量方程又可写成:式中ij τ代表了小涡对大涡的影响。
3.4亚格子模型通常在大涡模拟应用最广泛的亚格子模型主要有:动态混合模型、标准的Smagorinsky 模型、尺度相似模型、动态涡粘性模型、选择函数模型、梯度模型、选择函数模型等。
而应用最广泛的模型则是Smagorinsky 模型。
3.4.1亚格子涡扩散和涡粘模型不可压缩湍流的亚格子涡扩散和涡粘模型采用分子热扩散和分子粘性的形式,即 it i x T ∂∂=θκ kk ij ij t ij S τδντ312+= x d x u x x G x u i i '''-=⎰+∞∞-)()()(⎰+∞∞-=1)(dx x G ()(2)1i j ij i j i j u u S u P t x x x γρ∂∂⋅∂∂+=-+∂∂∂∂()(2)1i j ij ij i j i j j u u S u P t x x x x τγρ∂⋅∂∂∂∂+=-+-∂∂∂∂∂式中t κ和t ν分别为亚格子涡扩散系数和亚格子涡粘系数;(1/2)[(/)(/)]ij i j j i S u x u x =⋅∂∂+∂∂为可接尺度的变形率张量。
上式第2项的作用是满足不可压缩连续方程,当ij S 收缩时(ij S =0)公式等号左右两边相等。
把亚格子应力的涡粘模型公式代入到(13)式中,即得)])([()3(ij i i t i kk i i i j i x u x u x p x x u u t u ∂∂+∂∂+∂∂++∂∂-=∂∂+∂∂νντρ 0=∂∂ix u i 3.4.2 Smagorinsky 模型1963年,Smagorinsky 提出了Smagorinsky 模型,这是第一个亚格子模型。
Smagorinsky 模型认为亚格子应力的表达式如下:式中(1/2)[(/)(/)]ij i j j i S u x u x =⋅∂∂+∂∂是可接尺度的变形率张量,T ν是涡粘系数。
1963年Smagorinsky 定义了涡粘系数:) 式中2/1)2(ij ij S S S =是变形率张量,∆是过滤尺度,C S 无量纲参数,称为Smagorinsky 系数。
3.4.3 动态亚格子模式动态亚格子模式 由Germano 在1991年提出的,这种模式是以Smagorinsky ij T kk ij ij S ντδτ231-=-S C S T2)(∆=ν模式作为基本点模型,但却弥补了Smagorinsky 模式的一些不足。
动态亚格子模式其实是使亚格子涡粘模型的系数得到动态确定。
动态亚格子模式要两次对湍流场做进行过滤,第一次为细过滤,然后再进行粗过滤。
应力模型的系数是通过检验滤波器尺度条件和网格尺度计算得到的应力求差而得到的,使模型系数表示成时间和空间的函数,在模拟时就不需要对系数做调节。
在这点上比Smagorinsky 模式定系数计算更合理。
3.4.4 相似性模式尺度相似模式由Bardina于1980年提出。
此模式是假设主要是由大尺度脉动中的最小尺度脉动而产生的由大尺度脉动到小尺度脉动的动量输运,而且过滤掉的小尺度脉动速度和过滤后的最小尺度脉度速度相近。
由相似性假定和二次过滤能导出亚格子应力表达式。
应用该模式可以准确的预测墙壁面附近的渐近特性,但不适于模拟各向不均匀的的复杂流动,比如模拟室内空气的流动结果准确性较差。
3.4.5 混合模式混合模式中的亚格子应力是通过Smagorinsky 模式和相似性模式的叠加来而确定的。
这种模式湍动能耗散充足,而且和实际亚格子应力有紧密的联系,。
3.5大涡模拟的特点(1)可以描述小尺度湍流的流动,而且计算量要小于DNS方法,无论在工程应用还是科学研究都得到广泛的应用。
有很大的发展空间。
(2)用非均匀网格能够使网格数达到最少,在节省计算资源的同时又能保证计算的精度。
(3)湍流尺度比网格尺度小,可以模拟湍流发展过程的一些细节。
(4)相较于RANS方法,LES可以模拟更多的湍流大尺度运动,LES所用的湍流亚网格应力模型受边界的几何形状和流动类别的影响小,比RANS方法所用的Reynolds应力更具普适性。
其不足之处在于:(1)小涡模型网格节点的划分极密集,需要庞大的计算机存储能力;(2) 大量数据处理和非线性偏微分方程的求解需要高速数值处理能力;(3)仅用于比较简单的剪切流运动及管流。
(4) 由于实际湍流极其复杂,数值模拟仍需要非常可观的计算时间和实验经费。
4 直接数值模拟 (DNS)直接数值模拟(DNS)是不用任何湍流模型,直接求解完整的三维非定常的N - S 方程组,计算包括脉动在内的湍流所有瞬时运动量在三维流场中的时间演变。
4.1 控制方程用非定常的N - S 方程对湍流进行直接计算, 控制方程以张量形式给出:21i i i j i j i j i u u u p u f v t x x x x ρ∂∂∂∂+=-+∂∂∂∂∂=0i ju x ∂∂ 4.2 常用数值方法由于最小尺度的涡在时间与空间上都变化很快,为能模拟湍流中的小尺度结构,具有非常高精度的数值方法是必不可少的。
4.2.1谱方法或伪谱方法谱方法或伪谱方法是目前直接数值模拟用得最多的方法,其主要思路为,将所有未知函数在空间上用特征函数展开,成为以下形式:()()()()123(,)m n p m n p m n p V x t a t x x x ψϕχ=∑∑∑ (3) 其中m ψ,n ϕ与p χ,都是已知的正交完备的特征函数族。
在统计均匀性的空间方向或具有周期性一般都采用Fourier 级数展开,该特征函数族精度与效率最高。
除此之外,较多选用Chebyshev 多项式展开,它实质上是在非均匀网格上的Fourier 展开。
此外,也有用Legendre , Jacobi , Hermite 或Laguerre 等函数展开,但它们没有可用的快速变换算法。
若把该式代到N-S 方程组,一组()mnp a t 所满足的常微分方程组就能得到,可用通常的有限差分法来求解对时间的微分。
在用谱方法求解非线性项时,常用伪谱法来代替而直接求卷积。
伪谱法其实是配置法和谱方法相结合的一种方法,就是先把两量用Fourier 反变换回物理空间,然后在物理空间离散的配置点上去计算二者的乘积,再离散Fourier 变换回谱空间。
应用快速Fourier 变换(FFT)算法,伪谱法的计算速度比直接求两Fourier 级数卷积要快。
不过会存在“混淆误差”,就是在卷积计算时,会混进本应落在截断范围以外的高波数分量,引起数值误差。
严重时可使整个计算不正确甚至不稳定,但在多数情形下并不严重,且有一些标准的办法可用来减少混淆误差,但这将使计算工作量增加。
4.2.2高阶有限差分法高阶有限差分法的基本思想是利用离散点上函数值i f 的线性组合来逼近离散点上的导数值。
设i F 为函数()j f x ∂∂的差分逼近式,则j j jF f α=∑ 式中系数j α 由差分逼近式的精度确定,将导数的逼近式代入控制流动的N - S 方程,就得到流动数值模拟的差分方程。
差分离散方程必须满足相容性和稳定性。
4.3 直接数值模拟的特点(1)直接数值求解N-S 方程组,不需要任何湍流模型,因此不包含任何人为假设或经验常数。
(2)由于直接对N - S 方程模拟,故不存在封闭性问题,原则上可以求解所有湍流问题。
(3)能提供每一瞬时三维流场内任何物理量(如速度和压力)的时间和空间演变过程,其中包括许多迄今还无法用实验测量的量。
(4)采用数量巨大的计算网格和高精度流体力学计算方法,完全模拟湍流流场中从最大尺度到最小尺度的流动结构,描写湍流中各种尺度的涡结构的时间演变,辅以计算机图形显示,可获得湍流结构的清晰与生动的流动显示。
DNS 的主要不足之处在于:要求用非常大的计算机内存容量与机时耗费。
据Kim ,Moin &Moser 研究,即使模拟Re 仅为3300 的槽流,所用的网点数N 就约达到了6210⨯ ,在向量计算机上进行了250 h 。
5 不同方法的比较LES, DNS,RANS三种方法中DNS的计算量最大,LES的计算量介于另外两者之间,而RANS的计算量最小.计算量的多少主要取决于网格的数量、流场的时间积分长度、和最小旋涡的时间积分长度,网格的数量是计算量最重要的因素。