周期加热的Rayleigh-Benard对流时空结构

通过二维流体力学基本方程组的数值模拟,研究了普朗特数Pr=6.99时矩形渠槽周期加热对Rayleigh-Benard对流时空结构的影响．当水平流动强度Re=0时,发现稳定的由周期加热引起的局部定常对流．当Re比较小时,对流滚动抑制水平流动,获得了由周期加热引起的局部行波对流．当水平流动强度比较大时,由于周期加热与水平流动相互作用,水平流动抑制部分对流滚动,导致对流区域上游附近出现传导区域,对流区域减小,从而形成一种新的局部行波对流结构．并进一步讨论了Rayleigh-Benard对流时空结构的动力学特性．     

Rayleigh-Benard对流中小扰动的发展及对流斑图

通过二维流体力学的扰动方程组的数值模拟,探讨了分离比ψ=-0.2时,长高比Γ=30的矩形腔体中混合流体Rayleigh-Benard对流发生点附近扰动的成长和斑图的形成。结果表明:温度场线性成长阶段扰动的成长率γ_m是相对瑞利数r的函数,成长率γ_m随着相对瑞利数r的变化关系式为γ_m=0.9351r~(5.2039);在对流发生点附近的瞬态斑图取决于相对瑞利数r。给出了不同的相对瑞利数r(r分别为1.5、1.7、1.8)的情况下从小振幅到大振幅稳定状态的过渡过程中的两种不同的对流斑图,并讨论了其动力学特性。研究发现,当r较大时,存在行波与定常波共存的现象。     

利用格子Boltzmann方法数值模拟Rayleigh-Benard对流

采用格子Boltzmann方法对较大Rayleigh数范围下的二维Rayleigh-Benard对流进行了模拟研究.引入能量分布函数,利用该能量分布函数与粒子速度分布函数耦合来求解一个热流场,能量分布函数与粒子速度分布函数和Boltzmann方程构成了一个新的双分布格子Boltzmann模型.在考虑密度随温度变化的情况下,进行数值模拟,得到了Rayleigh-Benard对流速度、温度随时间的变化规律、系统的流线和等温线分布及平均Nusselt数与Rayleigh数的之间的关系,与相关文献数据进行了对比,模拟结果非常吻合,证明了改进的双分布格子Boltzmann模型的有效性.     

低Pr数流体非稳态自然对流流动的标度分析

施加了恒定热通量的竖直平板置于静止的均匀或线性分层流体中会在板壁上形成自然对流流 动. 它在达到稳态前将先历经起始和过渡两个阶段. 而表征其瞬时流动特征的主要参数是壁温、热边界层厚度、内层和外层速度边界层厚度、边界层内竖直方向的最大速度、以及层流边界层发展达到稳态所需的时间. 利用标度分析得到了均匀和线性分层低普郎特数(Pr<1)流体沿恒定热通量竖直平板的非稳态自然对流流动在各个发展阶段的标度关系.     

极小长高比腔体内混合流体Undulation行波对流

通过二维流体力学基本方程组模拟了具有较强Soret效应(分离比ψ=-0.47)的混合流体在极小长高比(Γ=4)腔体内的Rayleigh-Benard对流运动.研究了极小长高比行波对流的动力学特性,得到了稳定的Undulation行波存在的r值范围,给出了稳定后的Undulation行波摆动周期Tp的变化规律,分析了极小长高比行波对流的r依赖性及稳定性.首次发现极小长高比Γ=4时,与长高比Γ=12和Γ=8时相比,在腔体两端的滚动生成和消失的现象不再出现.讨论了长高比对Undulation行波向行波过渡的影响.     

低Pr数流体自然对流边界层流动的直接数值模拟

在太阳能采暖通风系统中,空气在太阳辐射作用下于吸热壁上所形成的非稳态自然对流边界 层流动是决定系统热性能的关键所在. 对于线性分层的低普郎特数($Pr<1$)流体而言,标度分 析表明在起始阶段边界层的发展是与高度无关而只依赖于时间的,而处于稳态时,各项标度 关系与处于起始阶段终结时刻的大不相同,且都与高度和时间无关,而只与$Pr$数和线性分 层度相关. 直接数值模拟计算结果表明由标度分析所得的各项标度关系揭示了流动的特性参 数与控制参数之间的决定性内在联系. 但同时也表明标度关系并没有充分体现出 对$Pr$数的所有依赖关系,这一不足可利用直接数值模拟计算结果得到弥补.     

具有水平流动的Rayleigh-Bénard对流进口段特性分析

采用SIMPLE算法对二维流体力学基本方程组进行了数值模拟,对普朗特数Pr=0.0272的具有水平流动的Rayleigh-Bénard对流进口段特性进行了研究。结果表明,在某些流体参数下具有水平流动的Rayleigh-Bénard对流由进口段和行波对流段组成,进口段长度取决于雷诺数和相对瑞利数。对于雷诺数Re=150的情况,当相对瑞利数r≤4时腔体内是水平流动;当r≥12时腔体内是行波对流;在40随着雷诺数变化的表达式。进口段长度随相对瑞利数的增大而减小,随着雷诺数的增大而增大。行波对流区的对流波数随着雷诺数的增大而减小。     

混合流体分离比对摆动行波的影响

利用Simple算法对流体力学基本方程组进行了数值模拟,研究了摆动行波的特性。结果表明:在摆动行波存在的范围内,当相对瑞利数r较小时,腔体内沿着空间的平均波数是随着时间周期变化的;当r较大时,腔体内沿着空间的平均波数随着时间增大保持为常数。摆动行波的摆动周期随r的增大而减小;分离比负值越大,变化越平缓,摆动行波出现的范围越大。随着长高比增加,较小分离比时,摆动周期的上、下限明显提高;分离比较大时,摆动行波存在区间对应的r上、下限明显增加。     

底部周期加热的局部对流斑图

运用Simple算法对二维流体力学基本方程组进行了数值模拟,探讨了普朗特数(Pr)为0.0272时矩形腔体底部周期加热对对流时空斑图的影响。当水平流动雷诺数(Re)为0时,发现了由正弦波周期加热引起的稳定的局部定常对流。当Re≠0时,由于正弦波周期加热与水平流动相互作用,获得了由正弦波周期加热和水平流动引起的局部行波对流。进一步比较和讨论了底部正弦波周期加热局部对流和混合流体Rayleigh-Benard局部对流的时空斑图,发现它们存在不同的机理。     

双流体Rayleigh—Benard对流中的局部行波

Rayleigh-Benard对流是研究非平衡对流的斑图（Pattem）及非线性动力学特性的典型模型之一。据此,通过流体力学基本方程组的数值模拟,探讨了分离比沙-0．4时双流体局部行波的形成过程;讨论了双流体局部行波流速场、温度场、浓度场、平均浓度流的结构和特性,分析了局部行波被局部化的原因。研究结果表明：局部行波始终在背离端壁方向传播并被限定在此端,它存在于腔体左端还是右端取决于经过过渡后对传波中最终控制腔体的一支行波的传播方向;在局部行波的存在范围内（1．519≤r≤1．604）,随相对瑞利数增加,衡量对流振幅的最大垂直流速、特征通过流体层的垂直热流量的努塞尔数、局部行波宽度等都在增加,但反映浓度特性的混合参数减小。     

Rayleigh-Benard Convection in a Near-Critical Fluid in the Neighborhood of the Stability Threshold

Steady-state Rayleigh-Benard convection in a medium with parameters close to the thermodynamic critical point is simulated within the framework of the complete Navier-Stokes equations with a two-scale representation of the pressure and the Van-der-Waals equation of state. A calibration relation is obtained for a realistic Rayleigh number in a compressible stratified medium. The parameters of the numerical simulation are determined from experimental data for near-critical helium on the basis of the calibration relation. The threshold Rayleigh numbers are found without and with allowance for stratification and a comparison with the experimental and theoretical data is carried out. The effect of compressibility of the near-critical fluid on steady-state convection flows is investigated beyond the stability threshold and the effect of adiabatic compression of the medium is analyzed.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, 2005, pp. 48–61.Original Russian Text Copyright © 2005 by Polezhaev and Soboleva.     

具有强SORET效应的混合流体 Undulation行进波对流斑图

本文通过流体力学基本方程组的数值模拟,探讨了具有强Soret效应（分离比ψ=-0．6）的混合流体Undulation行进波对流斑图的动力学特性。在相对瑞利数r〈6．436时,首次发现一种没有源缺陷的左右相对传播的CPW（Counter propagating waves）状态向行进波状态的过渡形式。在r=6．436—10.8的范围内,发现了两种不同结构的Undulation行进波对流斑图。当6．436〈r〈10时,出现了腔体内的平均波数在时间上变化且局部波数或当地波数在空间和时间上连续变化的Undulation行进波对流斑图。当r=10—10．8时,出现了腔体内的平均波数在时间上保持为常数而局部波数或当地波数在空间和时间上连续变化的Undulation行进波斑图。在两种状态下,Undulation行进波的摆动周期随瑞利数r增大而减小,它的对流振幅和Nusselt数随瑞利数r增大而增加。在Undulation行进波斑图形成以前,存在以中心为对称的Undulation行进波斑图,它的存活时间依赖于r。当r增加到11．0时,Undulation行进波过渡到定常对流状态。     

Rayleigh-Benard系统不稳定性的DSMC方法研究

本文应用DSMC方法初步研究了稀薄气体状态下Rayleigh-Benard系统的不稳定性特征，着重考察了气体密度梯度对Rayleihg-Benard系统不稳定性特征的影响。研究表明DSMC方法能够成功地模拟Rayleigh-Benard系统内涡的运动，得到了与外力方向相反的密度梯度会增加Rayleigh-Benard系统流动不稳定性的结果。     

具有强SORET效应的充分发展的混合流体行进波对流

混合流体Rayleigh-Benard对流是研究对流稳定性,时空结构和非线性特性的典型模型之一。本文利用流体力学扰动方程组的数值模拟,讨论了偏离传导状态具有强SORET效应的混合流体行进波对流的温度场和浓度场的成长过程,分析了充分发展对流情况下的对流振幅,Nusselt数及混合参数与相对瑞利数的关系。并给出了行进波相速度对相对瑞利数的依赖关系。结果说明混合参数的曲线与行进波相速度的分布曲线是类似的。文末,给出了垂直速度,温度和浓度场的分布并讨论了相对瑞利数对场的分布及不同场之间的相位差的影响。     

Application of GDQ Method for Study of Mixed Convection in Horizontal Eccentric Annuli

A generalized differential-integral quadrature (GDQ) discretization technique was used to solve a mixed heat convection problem in a body-fit coordinate system in its primitive variables form. A special treatment of the boundary condition to satisfy the continuity and momentum equations along the boundaries with the implementation of the GDQ method was investigated. Comparisons with the experimental and numerical results of other investigators are presented and discussed. In contrast with the existing published results, this highly accurate method was able to reveal extremely weak net circulation around the outer cylinder. In the horizontal annulus with the mixed heat convection problem the combination of unbalance of buoyancy and centrifugal forces causes net circulation. The net circulation decreases and approaches to zero with the rise of Rayleigh number, and it reaches its minimum value with high eccentricity when the inclination angle of eccentricity is π.     

Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel

We analyse the convection flow of a viscous fluid through a horizontal channel enclosing a fully saturated porous medium. The Galerkin finite element analysis is used to discuss the flow and heat transfer through the porous medium using serendipity elements. The velocity, the temperature distributions and the rate of heat transfer are analysed for variations in the governing parameters. The profiles at different vertical levels are asymmetric curves, exhibiting reversal flow everywhere except on the midplane. In a given porous medium, for fixed G or N, the temperature in the fluid region at any position in fluids with a higher Prandtl number, is much higher than in fluids with a lower Prandtl number. Likewise, other parameters being fixed, lesser the permeability of the medium, lower the temperature in the flow field. Nu reduces across the flow at all axial positions, while it enhances along the axial direction of the channel. Nu reduces with decrease in the Darcy parameter D, and thus lesser the permeability of the medium, lesser the rate of heat transfer across the boundary at any axial position of the channel.     

Forced and mixed convection boundary layer flow along a flat plate with variable viscosity and variable Prandtl number: new results

The steady laminar boundary layer flow along a flat plate is studied taking into account the variation of fluid viscosity and fluid Prandtl number with temperature. In the forced convection case the plate moves with constant velocity and its temperature varies in power law with x. In the mixed convection case the plate temperature is constant and the fluid moves upwards due to an external free stream and due to buoyancy forces. The results are obtained with the direct numerical solution of the boundary layer equations. The study concerns the wall heat transfer, the wall shear stress and velocity and temperature profiles across the boundary layer. The results of the present work are different from those existing in the literature, which have been obtained with the assumption of constant Pr number.     

Effect of a Standing Acoustic Wave on the Development of Large-Scale Convection in a Horizontal Layer

The effect of a standing acoustic wave on the development of long-wave convective perturbations in a horizontal layer with thermally insulated boundaries is investigated. The main two-dimensional flow is determined. A nonlinear amplitude equation with spatially-periodic coefficients is derived for investigating the stability of the main flow and secondary convection flows in the neighborhood of the stability threshold. The intensity of the acoustic field is assumed to be low. It is shown that the acoustic action leads to destabilization of the layer. Plane and three-dimensional perturbations are critical at large and small Prandtl numbers, respectively. Nonlinear one-dimensional steady-state solutions of the amplitude equation are obtained and their stability is investigated.