封闭水平矩形腔内流体自然对流第一次逆叉形分岔的数值研究

本文采用二阶全展开ETG(Euler-Taylor-Galerkin)分裂步有限元方法,对长宽比为3.5(L/B=3.5,如图 1所示)的封闭矩形腔体内,三种Pr数条件下,定常层流范围内,流体自然对流叉形分岔随Rayleigh数的演化过程进行了数值模拟。研究结果表明,该矩形腔内对应三种Pr数条件下,流体的叉形分岔的演化过程中,在第二次模态Ⅱ型叉形分岔之后,均会出现两个较小尺度涡旋合并,突变为一个较大尺度涡旋的全新叉形分岔模态。即在某临界Ra数两侧,存在定常四涡结构和定常三涡结构两个定常解支,当系统控制参数Ra越过临界值,前者被后者突发性取代,这是完全不同于传统叉形分岔的逆叉形分岔。其数值预报,则采用分半法结合流动拓扑结构及典型截面处速度扩线上鞍点的变化来确定。计算结果表明,在计算的Pr数条件下,随Pr数的增加逆叉形分岔对应临界Ra数的取值也会提高。     

不同属性肋片对复杂方腔自然对流影响的数值分析

采用二阶全展开Euler-Taylor-Galerkin分裂步有限元方法,在指定的网格密度条件下,在流动对应的普朗特数取为0.71,雷诺数取为104的情况下,数值分析了热肋、冷肋、上绝热肋、下绝热肋等四种不同属性肋片对封闭方腔内典型自然对流流动的影响.计算结果表明,肋片的存在对封闭方腔内的自然对流及相应的传热效率具有较强的影响,对流流动结构以及平均Nusselt数随肋片的属性发生较大的改变.     

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

长方腔自然对流第一次分岔突变现象的数值分析

在对网格密度作用进行详细分析的基础上,采用二阶全展开 Euler-Taylor-Galerkin分裂步有限元方法,对封闭水平矩形腔体内流体 (Pr=0.71)自然对流的第一次分岔过程进行了数值预报. 计算结果表明,第一次分岔相应的流动拓扑及临界Rayleigh ($Ra$)数随矩形腔体长宽比(W/B)取值的不同时会发生较大变化. 在所计算的长宽比取值范围内,封闭矩形腔内,流体自然对流第一次分岔拓扑的变化对应两种大的类型: 在较小的长宽比取值范围内(W/B\le 2.5),临界Ra数两侧,流动从单一涡核的定常流动突变成为具有不对称结构的定常双涡核运动, 在此范围内临界Ra数的取值随W/B取值的增加而减小;当对应长宽比取值2.6 \le W/B \le 4.0时,第一次分岔拓扑结构的变化呈现出更加复杂的特性,临界Ra数两侧流动从定常双涡核突变为定常非对称的三涡核流动,相应的临界Ra 数也随W/B的增加而减少. 而在区间[2.5,2.6]两端,临界Ra数的取值发生一次阶跃式突增,将该长宽比取值的区间定义为长方腔内该流体第一次分岔的突变区间.      

侧加热腔内的自然对流

开展侧加热腔内自然对流的研究具有重大的环境及工业应用背景. 总结侧加热腔内水平温差驱动的自然对流的最新研究进展, 并概述相应的流动性质、动力机制和传热特性以及对不同无量纲控制参数的依赖也有重要的科学价值. 已取得的研究结果显示突然侧加热的腔内自然对流的发展可包括初始阶段、过渡阶段和定常或准定常阶段. 不同发展阶段的流动依赖于瑞利数、普朗特数及腔体的高宽比, 且定常或准定常阶段的流态可以是定常层流流动、非定常周期性流动或者湍流流动. 此外, 回顾了对流流动失稳机制的研究成果以及湍流自然对流方面的新进展. 最后, 展望了侧加热腔内的自然对流研究的前景.      

热肋片宽度对复杂方腔自然对流影响的数值分析

采用笔者等发展的二阶全展开Euler-Taylor-Galerkin分裂步有限元方法,经过对算法精度验证,在确定不同网格密度对计算结果影响的基础上,确定出计算所需网格密度的必要条件,并在指定的网格密度条件下,分析了热肋片宽度对封闭方腔内典型自然对流流动的影响,流动对应的Pr=0.71,Ra=104。计算结果表明,随热肋...     

基于热格子Boltzmann方法的自然对流数值模拟

对Eggels和Somers提出的热格子Boltzmann格式进行了改进. 在不可压缩流动的假设下,提出了一种新的温度平衡分布函数,可以克服压缩性对温度统计的影响,并且相应地修正了统计宏观温度的方法. Eggels和Somers的方法对速度和温度均采用半步长反弹格式边界条件,适合无滑移的速度边界条件.但是对温度采用该边界条件在物理本质上显得不够准确,所以在边界上对二者统一采取算法既简单又容易实现的非平衡态外推格式,同时可以与Boltzmann格式的整体二阶精度保持一致. 最后,利用改进的热格子Boltzmann方法(TLBM)模拟了Ra=10^6和Pr=0.71(空气)的方腔中的自然对流,模拟得到的流动参数与其它数值方法的结果吻合得很好,表明改进的热格子Boltzmann方法可以有效准确地模拟非等温流动.      

三维方腔介电液体电对流的数值模拟研究

本文对封闭方腔内介电液体电对流进行了三维数值模拟研究.方腔的6个边界为固壁;4个侧边界为电绝缘边界;上下界面为两个电极.直流电场作用在从底部电极注入的自由电荷上,从而对液体施加库伦体积力并驱动流体流动形成电对流.为了求解这一物理问题,发展了一种二阶精度的有限体积法来求解完整的控制方程,包括Navier-Stokes方程和一组简化的Maxwell方程.考虑到电荷密度方程的强对流占优特性,采用了全逆差递减格式来求解该方程,获得了准确有界的解.通过研究发现,该流动在有限振幅区内的分叉类型为亚临界,即系统存在一个线性和非线性临界值,分别对应流动的开始和终止.由于非线性临界值比线性值小,因此两个临界值之间有一个迟滞回线.与无限大域中的自由对流相比,侧壁施加的额外约束改变了流场结构,使这两个临界值均有所增大.此外,还讨论了电荷密度和速度场的空间分布特征,发现电荷密度分布中存在电荷空白区.最后对更小空间尺寸情况计算结果表明,流动的线性分叉类型为超临界.本文的结果拓展了已有的二维有限空间内电对流的研究,并为三维电对流的线性和弱非线性理论分析提供参考.      

正弦波温度边界下多孔介质方腔内非热平衡对流传热数值模拟

采用局部非热平衡模型,在方腔左侧壁面温度正弦波变化、右侧壁面温度均一的边界条件下,通过SIM-PLER算法数值研究了固体骨架发热多孔介质方腔内的稳态非达西自然对流,主要探讨了不同正弦波波动参数N及方腔的高宽比M/L对方腔内自然对流与传热的影响规律。计算结果表明:正弦波温度边界使得方腔内的流场出现了复杂的变化,流体及固体区域左侧壁面附近出现了周期性的正负变化的温度场分布,左侧壁面局部Nusselt数出现了周期性的震荡现象;存在一个最佳温度波动参数N=1,此时多孔介质方腔内的整体散热量达到最大值;增加方腔高宽比会显著地削弱方腔内的自然对流传热过程,小高宽比也会在一定的程度上削弱多孔介质方腔内的对流传热。     

线圈绕Y轴倾斜时方腔内空气磁力与重力耦合对流的数值模拟

数值分析了微重力与重力环境下圆形载流线圈绕Y轴倾斜时方腔内空气热磁对流。磁场计算采用毕奥萨伐定律;控制方程基本变量采用控制容积法离散,求解采用SIMPLE算法。计算过程中Ra数的变化范围为104~105,线圈倾斜角yeuler的变化范围为-90°~90°,磁场力数γ数的变化范围为0~200。获得了空气热磁对流的流场和温...     

Transient natural convection in a cavity with heat input and a constant temperature wall on opposite sides

The transient natural convection of a fluid with Prandtl number of order 200 in a two-dimensional square cavity has been numerically studied. One of the vertical walls of the cavity is kept at a constant (ambient) temperature and a constant heat flux is applied on the opposite wall. The other walls are adiabatic. Initially, a boundary layer is formed near the heated wall; subsequently, a large vortical structure is generated, together with an upper intrusion layer. As time progresses, the average temperature in the cavity increases, and a descending boundary layer is formed near the constant temperature wall. During the transition to the steady-state regime, a thermal stratification pattern is formed. The results are compared with the scale analysis presented by Patterson and Imberger (1980).     

Stability and slightly supercritical oscillatory regimes of natural convection in a 8:1 cavity: solution of the benchmark problem by a global Galerkin method

The global Galerkin method is applied to the benchmark problem that considers an oscillatory regime of convection of air in a tall two‐dimensional rectangular cavity. The three most unstable modes of the linearized system of the Boussinesq equations are studied. The converged values of the critical Rayleigh numbers together with the corresponding oscillation frequencies are calculated for each mode. The oscillatory flow regimes corresponding to each of the three modes are approximated asymptotically. No direct time integration is applied. Good agreement with the previously published results obtained by solution of the time‐dependent Boussinesq equations is reported. Copyright © 2004 John Wiley & Sons, Ltd.     

Large-scale flow structures and heat transport of turbulent forced and mixed convection in a closed rectangular cavity

Results of an experimental study of flow structure formation and heat transport in turbulent forced and mixed convection are presented. The experiments were conducted in a rectangular cavity with a square cross section, which has an aspect ratio between length and height of Γ_xz = 5. Air at atmospheric pressure was used as working fluid. The air inflow was supplied through a slot below the ceiling, while exhausting was provided by another slot, which is located directly above the floor. Both vents extend over the whole length of the cell. In order to induce thermal convection the bottom of the cell is heated while the ceiling is maintained at a constant temperature. This configuration allows to generate and study mixed convection under well defined conditions. Results of forced convection at Re = 1.07 × 10⁴ as well as mixed convection at 1.01 × 10⁴ ? Re ? 3.4 × 10⁴ and Ra = 2.4 × 10⁸ (3.3 ? Ar ? 0.3), which were obtained by means of Particle Image Velocimetry and local temperature measurements, are presented. For purely forced convection a 2D mean wind, which can be approximated by a solid body rotation, is found. With increasing Archimedes number this structure becomes unstable, leading to a transition of the solid body rotation into additional smaller convection rolls. Proper orthogonal decomposition of the instantaneous velocity fields has been performed for further analysis of these coherent large-scale structures. Their fingerprint is found in the spatial temperature distribution of the out flowing air at the end of the outlet channel, which reveals a temporally stable profile with two maxima over the length of the outlet. Moreover a maximum in the global enthalpy transport by the fluid is found at Ar ≈ 0.6.     

LES-based filter-matrix lattice Boltzmann model for simulating turbulent natural convection in a square cavity

In this paper, a novel thermal filter-matrix lattice Boltzmann model based on large eddy simulation (LES) is proposed for simulating turbulent natural convection. In this study, the Vreman subgrid-scale eddy-viscosity model is introduced into the present framework of LES to accurately predict the flow in near-wall region. Two dimensional numerical simulations of natural convection in a square cavity were performed at high Rayleigh number varying from 10⁷ to 10¹⁰ with a fixed Prandtl number of Pr = 0.71. The influences of the higher-order terms upon the present results at high Rayleigh numbers are examined, taking Ra = 10⁷ and 10⁸ as the example, revealing that the proper minimization of the higher-order terms can improve numerical accuracy of present model for high Rayleigh convective flow. For the turbulent convective flow, the time-averaged quantities in the median lines are presented and compared against those available results from previous studies. The general structure of turbulent boundary layers is well predicted. All numerical results exhibit good agreement with the benchmark solutions available in the previous literatures.     

Experimental observations of the thermal flow around a square obstruction on a vertical wall in a differentially heated cavity

The transient thermal boundary layer flow around a square obstruction placed at the middle of the hot wall in a differentially heated cavity is visualized using a shadowgraph technique. The results show that the thermal boundary layer flow, which is blocked by the obstruction, firstly forms an intrusion head under the obstruction (the lower intrusion head). Subsequently, the lower intrusion head bypasses the obstruction and reattaches to the down-stream boundary. During the reattachment process, a more complicated flow is induced, and eventually both the lower intrusion head and the thermal boundary layer destabilize. After the lower intrusion head is convected away, the thermal boundary layer flow re-stabilizes. At the quasi-steady state, the thermal boundary layer forms a double-layer structure, which is split into two sections by the obstruction. It is demonstrated that both the transient processes and the quasi-steady state flow structures of the thermal boundary layer are significantly altered by the obstruction in comparison with the case without the obstruction.     

Experiment and prediction of a cavity type separated flow

A detail study involving flow visualization, Laser Doppler Velocimeter (LDV) measurements and numerical prediction is presented. The visualization experiments revealed striking results of a pulsatile motion in the separated flow region associated with the formation and passage of large eddy structures. Measurements of mean velocities and turbulence intensity profiles across the separated flow field, provided information about the separated shear layer development and the recirculating flow pattern. The numerical predictions, obtained with a two-layer turbulence model in conjunction with the SIMPLE algorithm, failed to reproduce the coherent eddies and the pulsatile motion, but the mean velocities are reasonably reproduced.     

Numerical study of buoyancy- and thermocapillary-driven flows in a cavity

Thermocapillary- and buoyancy-driven convection in open cavities with differentially heated endwalls is investigated by numerical solutions of the two-dimensional Navier-Stokes equations coupled with the energy equation. We studied the thermocapillary and buoyancy convection in the cavities, filled with low-Prandtl-number fluids, with two aspect-ratiosA=1 and 4, Grashof number up to 10⁵ and Reynolds number ⋎Re⋎≤10⁴. Our results show that thermocapillary can have a quite significant effect on the stability of a primarily buoyancy-driven flow, as well as on the flow structures and dynamic behavior for both additive effect (i.e., positiveRe) and opposing effect (i.e., negativeRe).