紊流积分尺度实用识别算法及其数值验证
- 格式:docx
- 大小:45.17 KB
- 文档页数:15
一、概述空调系统在现代建筑中扮演着重要的角色,确保室内空气质量和舒适度。
而空调房间气流组织的三维紊流计算则是评估空调系统效能和设计最佳气流分布的关键步骤。
本文将着重探讨空调房间气流组织的三维紊流计算的相关理论、方法和应用。
二、空调房间气流组织的相关理论1. 三维气流模型空调房间气流组织是一个复杂的三维问题,需要建立相应的三维气流模型。
这个模型需要考虑室内外气流的相互作用、温度差异、墙面和家具对气流的阻挡等因素。
2. 紊流模拟方法在三维气流模型中,紊流模拟是非常重要的一步。
常用的方法有LES (大涡模拟)和RANS(雷诺平均Navier-Stokes方程)。
选择合适的模拟方法对于准确地模拟室内气流非常关键。
三、空调房间气流组织的三维紊流计算方法1. 网格划分在进行三维紊流计算之前,首先需要对空间进行网格划分。
通常情况下,使用结构化或非结构化网格来划分室内空间,确保在整个计算空间内都有足够的网格密度。
2. 初始和边界条件设定确定好初始条件和边界条件对于三维紊流计算非常关键。
初始条件包括初始速度场和温度分布,边界条件包括入口和出口的气流速度、温度和湿度等参数。
3. 紊流模拟软件的选择目前市面上有很多用于三维紊流计算的软件,例如ANSYS Fluent、OpenFOAM等。
选择合适的软件对于三维紊流计算的准确性和效率都有很大的影响。
4. 寻找最佳气流分布三维紊流计算的最终目的是寻找最佳的气流分布,以确保室内空气的均匀性和舒适度。
通过对计算结果进行分析和比较,可以找到最佳的室内气流组织方案。
四、空调房间气流组织的三维紊流计算的应用1. 空调系统设计与优化通过三维紊流计算,可以对空调系统的设计和布局进行优化。
确保室内空气的流动均匀性和热舒适度,提高空调系统的效能和节能性能。
2. 室内空气质量评估室内空气质量对于人员的健康与舒适度有着重要的影响。
通过三维紊流计算,可以评估室内空气的CO2浓度和PM2.5等污染物的扩散情况,确保室内空气的新鲜度和清洁度。
空气滤清器内三维紊流数值计算与分析的开题报告一、研究背景空气滤清器是用于净化空气的设备,在汽车、工业设备、建筑物等领域广泛应用。
滤清器内部极易积聚灰尘,降低净化效率,并可能造成阻力增加、耗能增大等问题。
因此,空气滤清器内部的流场特性和清洁性能研究,对于提高其净化效率和降低能耗具有重要意义。
二、研究内容与目标本研究旨在通过数值模拟方法,对空气滤清器内部的三维流场特性进行研究与分析,探究其对滤清器性能的影响,并提出改进措施。
具体研究内容包括:1. 建立空气滤清器的三维几何模型,并对其流场特性进行数值模拟计算。
2. 分析空气滤清器内部的流场特性,如温度、速度、压力等。
3. 探究滤芯构造、过滤介质等因素对滤清器性能的影响。
4. 提出改进空气滤清器性能的措施。
三、研究方法与技术路线本研究主要采用数值模拟方法和CFD(Computational Fluid Dynamics,计算流体动力学)软件进行模拟分析。
具体技术路线如下:1. 建立三维几何模型:通过CAD(Computer-Aided Design,计算机辅助设计)软件建立空气滤清器几何模型,并导入CFD仿真软件中进行后续模拟计算。
2. 数值模拟计算:使用CFD软件对滤清器内部的三维流场进行模拟计算,并得出相关的流场参数,如温度、速度、压力等。
3. 数据分析与处理:对模拟计算得到的数据进行分析处理,并通过图表等方式展示研究结果。
4. 结果验证与改进:根据模拟计算得到的结果,进行实验验证,并提出改进性能的措施。
四、预期成果通过数值模拟计算和实验验证,本研究将得到以下预期成果:1. 获得空气滤清器内部的三维流场特性,如温度、速度、压力等参数。
2. 探究各种因素对空气滤清器性能的影响,提出空气滤清器设计与改进的建议。
3. 建立数值模拟计算的方法与技术,为类似领域的研究提供一定的参考价值。
五、研究难点1. 空气滤清器内部流场特性的复杂性,需要提出合理的数值模拟方案。
三维圆管紊流流动状况的数值模拟分析在工程和生活中,圆管内的流动是最常见也是最简单的一种流动,圆管流动有层流和紊流两种流动状况。
层流,即液体质点作有序的线状运动,彼此互不混掺的流动;紊流,即液体质点流动的轨迹极为紊乱,质点相互掺混、碰撞的流动。
雷诺数是判别流体流动状态的准则数。
本研究用CFD 软件来模拟研究三维圆管的紊流流动状况,主要对流速分布和压强分布作出分析。
1 物理模型三维圆管长2000mm l =,直径100mm d =。
流体介质:水,其运动粘度系数62110m /s ν-=⨯。
Inlet :流速入口,10.005m /s υ=,20.1m /s υ= Outlet :压强出口Wall :光滑壁面,无滑移2 在ICEM CFD 中建立模型2.1 首先建立三维圆管的几何模型Geometry2.2 做Blocking因为截面为圆形,故需做“O ”型网格。
2.3 划分网格mesh注意检查网格质量。
在未加密的情况下,网格质量不是很好,如下图因管流存在边界层,故需对边界进行加密,网格质量有所提升,如下图2.4 生成非结构化网格,输出fluent.msh等相关文件3 数值模拟原理紊流流动当以水流以流速20.1m /s υ=,从Inlet 方向流入圆管,可计算出雷诺数10000υdRe ν==,故圆管内流动为紊流。
假设水的粘性为常数(运动粘度系数62110m /s ν-=⨯)、不可压流体,圆管光滑,则流动的控制方程如下:①质量守恒方程:()()()0u v w t x y zρρρρ∂∂∂∂+++=∂∂∂∂ (0-1)②动量守恒方程:2()()()()()()()()()()[]u uu uv uw u u ut x y z x x y y z z u u v u w p x y z xρρρρμμμρρρ∂∂∂∂∂∂∂∂∂∂+++=++∂∂∂∂∂∂∂∂∂∂'''''∂∂∂∂+----∂∂∂∂ (0-2)2()()()()()()()()()()[]v vu vv vw v v v t x y z x x y y z z u v v v w px y z yρρρρμμμρρρ∂∂∂∂∂∂∂∂∂∂+++=++∂∂∂∂∂∂∂∂∂∂'''''∂∂∂∂+----∂∂∂∂ (0-3)2()()()()()()()()()()[]w wu wv ww w w w t x y z x x y y z z u w v w w px y z zρρρρμμμρρρ∂∂∂∂∂∂∂∂∂∂+++=++∂∂∂∂∂∂∂∂∂∂'''''∂∂∂∂+----∂∂∂∂ (0-4)③湍动能方程:()()()()[())][())][())]t t k k t k k k ku kv kw k k t x y z x x y yk G z zμμρρρρμμσσμμρεσ∂∂∂∂∂∂∂∂+++=+++∂∂∂∂∂∂∂∂∂∂+++-∂∂ (0-5)④湍能耗散率方程:212()()()()[())][())][())]t t k k t k k u v w t x y z x x y y C G C z z k kεεμμρερερερεεεμμσσμεεεμρσ∂∂∂∂∂∂∂∂+++=+++∂∂∂∂∂∂∂∂∂∂+++-∂∂ (0-6)式中,ρ为密度,u 、ν、w 是流速矢量在x 、y 和z 方向的分量,p 为流体微元体上的压强。
流体力学紊流知识点总结一、流体力学紊流的基本概念1. 流体流体是一种可以流动的物质,包括液体和气体两种形态。
在流体力学中,流体的流动规律被广泛研究,紊流的产生和演化也与流体的特性密切相关。
2. 紊流紊流是一种混乱而不规则的流动现象,它具有高度的不可预测性和随机性,是流体中的激动运动状态。
紊流现象可在各种条件下产生,如在管道中的水流、空气中的湍流等,它的产生和演化是非常复杂的,需要借助数学模型和实验研究来揭示其规律。
3. 紊流的特征紊流的主要特征包括混沌性、不可预测性、漩涡结构和能量传递等。
混沌性是指紊流的运动轨迹是不规则的、无序的,不可预测性则是指紊流的演化是随机的、无法准确预测的。
而漩涡结构和能量传递则是紊流内部的重要特征,漩涡结构是形成和维持紊流的基本元素,而能量传递则是紊流中的主要演化机制。
4. 紊流的产生紊流的产生有多种方式,主要包括不稳定性机制、非线性机制和随机机制等。
不稳定性机制是指流体在特定条件下出现的不稳定现象而形成紊流,如雷诺数超过临界值时出现的湍流现象;非线性机制是指流体在非线性条件下产生紊流,如高速运动的流体中的雷诺应力非线性效应等;随机机制则是指由于流体粘性的随机性而导致的紊流现象。
这些产生机制同时也是研究紊流的重要方面。
二、流体力学紊流的数学描述1. 紊流的描述方法紊流的描述通常采用Navier-Stokes方程组和湍流模型,其中Navier-Stokes方程组描述了流体的运动规律,湍流模型则用来描述流体中的湍流演化。
这两种描述方法结合起来可以比较全面地揭示紊流的演化规律。
2. Navier-Stokes方程组Navier-Stokes方程组是描述流体的基本动力学方程,它包括质量守恒方程、动量守恒方程和能量守恒方程三个方程。
这些方程描述了流体的压力、速度、密度等物理量随时间和空间的变化规律,是研究流体运动的重要工具。
3. 湍流模型湍流模型是描述流体中的湍流演化的数学模型,主要包括Reynolds平均方程、湍流能谱方程和湍流能量方程等。
计算流体力学常用数值方法简介李志印 熊小辉 吴家鸣(华南理工大学交通学院)关键词 计算流体力学 数值计算一 前 言任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。
利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。
计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。
一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。
随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。
经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。
现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。
此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。
随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。
目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。
紊流积分尺度实用识别算法及其数值验证王峰峰;赵林;曹曙阳;葛耀君【摘要】对比现有的几种紊流积分尺度算法,针对其在某些来流条件下存在的系统偏差,提出了基于泰勒假定修正的紊流积分尺度识别算法.采用谐波合成法数值化地再现多组基于Von-Karman谱的宽频紊流及窄带单频的风速时程序列,利用该算法对紊流积分尺度进行识别,并将识别结果与预期理论值进行比较,提出基于时间尺度修正系数项的实用紊流积分尺度识别算法并进行了验证,得出对于宽带紊流,时间尺度修正系数适合取为2/3,对于窄带紊流,时间尺度修正系数适合取为6.结果表明:本文方法大大提高了紊流积分尺度的计算精度,具有很好的工程应用价值.%Through the comparison of multiple algorithms of turbulence integral scale, the correctional recognition algorithm based on Taylor assumption was proposed. The harmonic synthesis method was used to generate several broadband time-series based on Von-Karman spectrum and several single-frequency time-series. Through the turbulence integral scale calculation of these signals using the correctional algorithm, the recognition results were compared with the expected theoretical values and the correctness of the algorithm was validated. The correction factor obtained for the broadband turbulence was 2/3, and for the narrowband turbulence was 6. The results show that this method greatly improves the calculation accuracy of the turbulence integral scale, which has a good value in engineering application.【期刊名称】《建筑科学与工程学报》【年(卷),期】2012(029)004【总页数】7页(P67-73)【关键词】紊流积分尺度;泰勒假定;谐波合成法;识别算法【作者】王峰峰;赵林;曹曙阳;葛耀君【作者单位】同济大学土木工程防灾国家重点实验室,上海200092;同济大学土木工程防灾国家重点实验室,上海200092;同济大学土木工程防灾国家重点实验室,上海200092;同济大学土木工程防灾国家重点实验室,上海200092【正文语种】中文【中图分类】TU279.70 引言自然风实际上包含平均风和脉动风2个部分。
脉动风是由于风的不规则性引起的,描述风速脉动特性的参数主要有紊流强度、紊流积分尺度、功率谱密度函数等。
桥梁抖振分析中,理论分析、风洞试验和工程实测之间存在一定的误差,该误差的产生原因主要来自于空气动力学中脉动风特性参数的确定和模拟。
紊流积分尺度是气流中湍流涡旋平均尺寸的量度,其大小决定了脉动风对结构的影响范围,因而在结构风荷载分析中具有不可忽略的意义。
通常认为当来流脉动涡旋的尺寸大于结构物的尺度时,脉动风在结构各个部位引起的动荷载会相互叠加;反之,若脉动涡旋的尺寸比结构物的尺寸小,脉动风在各个部位引起的动荷载可能会抵消。
紊流积分尺度的分析结果与数据记录的长度和平稳程度有关,不同的实测结果也有一定差别,而且,由于紊流结构的差别,即使是同一数据记录采用不同的分析方法,也可能得到不同的结果。
为了进行建筑结构随机振动精细化研究,对脉动风特性参数合理取值至关重要。
本文中笔者列举了多种紊流积分尺度的计算方法,考虑到现有算法存在一定的系统误差,基于谐波合成法数值模拟及验证过程,提出了一种修正的紊流积分尺度的实用识别算法,并对算法的适用性进行了数值验证。
1 紊流积分尺度1.1 紊流积分尺度的定义通过空间某一点(x1,y1,z1)的气流中的速度脉动可以被认为是由平均风所输运的各种尺度的涡旋在该点造成的、按各自不同周期脉动的速度分量的迭加。
大气边界层湍流中的每个涡旋尺度可以看作在那一点引起了频率n的周期脉动,与波相似,因此可以定义涡旋的波长λ=U/n,U 为平均风速,这个波长就是涡旋大小的尺度,涡旋的尺度及湍流脉动能量在不同尺度水平上的分布决定了湍流的结构特征。
紊流积分尺度就是脉动风中湍流涡旋平均尺寸的量度。
由于涡旋具有三维特性,因此对应脉动风速和空间各包含3个方向,u,w,v是脉动风速的纵向、横向和垂直方向,x,y,z是空间的3个方向。
共有9个紊流积分尺度:Lxu,Lyu,Lzu,Lxw,Lyw,Lzw,Lxv,Lyv,Lzv,分别量度与纵向脉动速度有关的涡旋在纵向、横向和垂直方向的平均尺寸。
紊流积分尺度Lra的数学表达式定义为式中:a=u,v,w;r=x,y,z;σ2a 为脉动风速分量a的方差;Ca1a2(r)为相距为r的两点上的脉动风之间的互协方差函数。
纵向脉动风速u在x方向上的紊流积分尺度Lxu为式中:Cu1u2(x)为2个纵向脉动速度u1,u2 的互协方差函数,u1=u(x1,y1,z1,t),u2=u(x1+x,y1,z1,t),t为时间,σ2u为纵向脉动风速的方差。
1.2 泰勒假设按照式(2)来计算紊流积分尺度就必须在大尺度空间内多点同步测量脉动风速,这在实际操作中是非常难实现的。
在实际应用中,一般假设紊流中的涡旋是不衰减地以平均风速向下游传输,则脉动风速u(x1,t+τ)可以定义为u(x1-x/U,τ),这就是泰勒假设。
同时,由于脉动风速本身的平均值为0,自协方差函数可以用自相关函数代替。
根据这一假设,式(2)可以改写为式中:Ru(τ)为来流网速的自相关函数。
泰勒假定的引入将多点测量转化为单点测量,并且用自相关函数代替了空间互相关函数,使紊流积分尺度的求解得到简化。
1.3 经验公式对于脉动风紊流积分尺度的计算,目前已经有一定的研究成果。
Shiotani[1]对开阔地形的脉动风紊流积分尺度进行了测量,结果表明,不同高度脉动风紊流积分尺度有明显差异。
在接近地面的区域,风速和紊流特性往往取决于z/z0,同时受风速影响较小,其中z0为地表粗糙度,是影响近地边界层风特性的最主要的参数,z为距离地表的高度。
对于远离地面的区域,紊流积分尺度是随着风速的增加而增加的。
因此,对于一些低矮的建筑物,在缺少实测数据时,可以用一些简单的经验公式来估算紊流积分尺度。
Flay等[2]测量了垂直风向的长度尺度,并总结了积分尺度的分析方法。
Counihan[3]指出,根据1972年以前的紊流积分尺度实测结果,紊流积分尺度可以描述为地面粗糙度的递减函数,对于高度z=10~240m,认为纵向紊流积分尺度的经验公式为式中:系数C和指数m均为地表粗糙度的函数。
Simiu等[4]指出,根据式(4)计算得到的Lxu 值是实测值的2倍。
对于Lyu,文献[5]中建议取Lxu的0.2倍更好,即对于Lzu,Bietry提出了以下不考虑地表粗糙度的经验公式Teunissen[6-7]提出了在乡村地貌观测的纵向和垂直方向的紊流积分尺度Lxu,Lxw,得到的经验公式分别为Dyrbye等[8]给出的紊流积分尺度计算公式为1.4 计算方法由于紊流积分尺度是与湍流空间相关性关联的参数,理想的分析方法是在空间上多点同时测量,然后由式(1)得到各个方向的紊流积分尺度。
实际中空间多点同时测量较难实现,往往利用泰勒假定将多点测量简化为单点测量。
各国学者对紊流积分尺度提出了各种不同的方法,常用的方法有以下6种:(1)根据紊流积分尺度的定义式(1)进行各个方向的紊流积分尺度计算。
这是紊流积分尺度最直接、最理想的计算方法,但由于该方法需要对空间多点进行同步测量,实际中较少采用。
(2)引入泰勒假定,将空间多点测量变为单点测量,用自相关函数代替空间相关函数,即根据式(3)进行计算。
文献[5]中认为,式(3)的积分上限取Ru(τ)=0.05σ2u 为最佳,这一结论得到了多数学者的认可。
(3)Reed等[9]采用AR线性滤波模型计算紊流积分尺度。
将大气湍流作为理想的随机平稳过程处理,此脉动速度可以由一系列白噪声信号通过线性滤波系统产生,从而建立自拟合模型[10]。
经过比较确定了合适的延迟时间和线性阶次,当平均风速大于10m·s-1时,建议采用二阶线性滤波模型,当平均风速小于10m·s-1时,一阶线性滤波比较合适。
对脉动风速的时间序列进行相关分析可以得到紊流积分尺度。
对于一阶和二阶AR模型,可以表示为:一阶AR模型二阶AR模型式中:φ1,φ2 均为 AR模型系数;ut,ut-1,ut-2分别为t,t-1,t-2时刻的脉动风速;at 为白噪声序列。
应用泰勒假定式(3)计算紊流积分尺度,则有建立AR模型系数φ1,φ2与Ψ的关系,得到:一阶AR模型二阶AR模型式中:p为任意常数;R1 为方程式R21-R1φ1-φ2=0的大根。
将式(12),(13)代入式(11),可以得到紊流积分尺度的计算公式分别为:一阶AR模型二阶AR模型(4)Davenport[11]建议空间相关函数服从指数衰减率,根据泰勒假设,自相关函数也应服从指数衰减率,因此式(3)可以改写为式中:α为拟合的指数衰减系数。
(5)文献[2]和文献[12]中指出,如果脉动风谱符合Von-Karman谱,紊流积分尺度可以从功率谱密度函数直接求出,即式中:Si(0)为对应脉动风速的功率谱在f=0时的值,f为脉动风速对应的频率;σi为对应脉动风速的均方根。
(6)脉动风谱拟合法。
目前被广为接受的脉动风谱有 Von-Karman谱、Davenport谱、Simiu谱和Harris谱等。
假设脉动风谱符合其中的Von-Karman谱,将风谱模型中的紊流积分尺度当作待定参数,根据实测脉动风谱模型进行最小二乘拟合,可以确定出它的紊流积分尺度。
文献[13]中指出,在求某段数据的自相关函数时,实际上利用了矩形窗函数对数据进行截断。
矩形窗函数不仅影响了原来信号在时域的形状,也影响了其在频域的形状,所以采用不同记录长度的数据计算自相关函数时会导致计算结果的波动。