三种不同对流结构的行波斑图

利用二维数值分析,探讨了长高比Γ=20、分离比ψ=-0.4的三种行波对流斑图。结果表明:在r?(1.67,2.0]范围内出现了具有两个间歇性缺陷的行波斑图,第一缺陷和第二缺陷发生的位置固定;第一缺陷的出现周期随着相对瑞利数r的增加而增加。当相对瑞利数r较小时,第二缺陷的出现周期不确定;当相对瑞利数r较大时,第二缺陷的出现周期随着相对瑞利数r的增加而增加。在r?(2.0,2.59]范围内出现了具有一个间歇性缺陷的行波斑图,缺陷发生的位置不固定;缺陷的出现周期随着相对瑞利数r的增加而增加。在r?(2.59,4.6]范围内出现无缺陷的行波斑图,这说明随着相对瑞利数r的增加,行波对流结构变得简单化;同时发现不同的行波对流结构有不同的对流振幅变化过程。     

具有间歇性缺陷的混合流体行进波对流斑图

本文通过流体力学基本方程组的数值模拟,探讨了具有中等Soret效应的混合流体行进波斑图的动力学特性.当分离比Ψ=-0.3时,首次发现一种没有源缺陷的左右相对传播的CPW(Counter propagating waves)状态向行进波状态的过渡形式,并且在r=1.50-1.60的范围内,行进波对流斑图中存在着间歇性缺陷结构.这种缺陷出现的周期随瑞利数r增大而增加.在缺陷出现的周期内,对流振幅也以行进波的周期在周期的变化,对流振幅的振动次数或行进波的周围数也随相对瑞利数r增大而增加.当r增加到1.65时,行进波对流斑图中的缺陷结构消失.由于缺陷引起的对流振幅的周期性变化也随之消失,而以行进波的周期在整个时间段上周期的振动.     

具有水平流动的Rayleigh-Benard对流斑图

本文利用数值模拟,探讨了普朗特数 时有水平流动的Rayleigh-Benard对流结构．当水平流动强度 时,发现定常对流的多重稳定性．当 时,Rayleigh-Benard(RB)对流中存在三种对流斑图．它们的出现依赖于水平流动强度 和相对瑞利数 ．与 相比,在 时的行波具有不同的动力学特性．     

具有强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行进波过渡到定常对流状态。     

大长高比腔体中的摆动行波

通过数值求解流体力学方程组,探讨了大长高比Γ(28)20的腔体中的摆动行波。研究结果表明:对于较小的相对瑞利数r,两端壁处有滚动产生且摆动行波消失,波长在空间上变化较大;对于较大的r,两端壁处不再有滚动产生且摆动行波消失,平均波数不再随着时间变化;随着r的增大,摆动幅度减小,摆动周期变小;随着长高比Γ的增加,摆动行波存在的范围增大;相同相对瑞利数情况下,长高比Γ较大的腔体,摆动行波存在的周期较大。     

小长高比腔体内的摆动行波

通过流体力学基本方程组的数值模拟,探讨了具有Soret效应(分离比ψ=?0.47)和小长高比(Γ=8)腔体中混合流体摆动行波对流的动力学特性。研究表明:在相对瑞利数r3.467时系统出现了行波状态;在r=3.647～6.227的范围内,发现了摆动行波对流;且对流振幅随着时间的变化存在两种不同特性,其摆动周期随瑞利数r增大而减小,对流振幅和努塞尔数随瑞利数r增大而增加;当r增大到r=6.228时,摆动行波过渡到定常对流状态。因此,在行波对流向定常对流过渡的过程中存在摆动行波对流.     

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较大时,存在行波与定常波共存的现象。     

缺陷源摆动的对传波

基于流体力学方程组，对长高比 腔体内混合流体对流中缺陷源摆动的对传波的特性进行了数值模拟．结果发现，当分离比 时，缺陷源摆动的对传波的缺陷源由两个对流圈同时被拉长同时被分裂后形成．对传波的存在区间为 ，对传波的摆动周期较小并且基本都稳定在 ．对传波的摆动幅度较小并且几乎不随着相对瑞利数 变化．当 时，对传波的缺陷源下方対流圈轴线与缺陷源移动方向保持一致，没有对传波分支产生．对传波的缺陷源上方存在多个与缺陷源移动方向几乎垂直的対流圈的轴线，存在对传波分支，在对传波分支上出现间歇性的缺陷．从而形成具有单侧缺陷的缺陷源摆动的对传波．对传波的存在区间为 ．对传波的摆动周期为 并且随着相对瑞利数的增加而增加．对传波的摆动幅度较大，并且随着相对瑞利数的增加而减小．     

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

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

矩形倾斜腔体中的两种对流斑图和分区

为了研究矩形倾斜腔体中普朗特数Pr=0.72的流体对流斑图和分区,本文基于流体力学方程组进行了数值模拟。在相对瑞利数r=6.0的情况下,观察了倾角θ=10°和θ=60°时对流斑图随着时间的发展,发现系统存在单圈型对流和多圈型对流两种斑图。流线随着倾角的变化说明:随着倾角增加,对流圈数逐渐减少,对流波长逐渐增加,对流波数减小;然后,随着对流圈数显著地减少,系统由多圈型对流过渡到单圈型对流。根据模拟计算结果,给出了多圈型对流过渡到单圈型对流的临界倾角θc随着相对瑞利数r变化的关系曲线。对流在θ-r平面上分为两个区域:θ<θc时系统是单圈型对流,θ>θc时系统是多圈型对流。对流最大振幅A和努塞尔数Nu随着倾角θ的变化曲线被临界倾角θc分成两段,它们有不同的变化规律。因此,临界倾角也可以由对流最大振幅A或努塞尔数Nu的变化曲线来确定。     

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

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

腔体高度对Rayleigh-Bénard对流的影响的二维数值分析

腔体高度的改变会对Rayleigh-Bénard对流产生影响.本文利用二维流体力学基本方程组进行数值模拟,在相对瑞利数r(28)5,腔体高度0 0 0 0 0d(28)1d,1.4d,1.8d,2d,3d(cm10d(28))时,研究了对流运动的动力学特性,发现随着腔体高度d的增加,稳定对流的对流圈数与之呈线性反比例关系,而对流波数k随之基本保持恒定;对流稳定后,对流振幅和努塞尔数随腔体高度d增大而减小,且得出对流振幅和努塞尔数随着相对腔体高度0/dd变化的函数关系式.在腔体高度0 0 0d(28)1d,1.4d,1.8d情况下,对流振幅和努塞尔数都随相对瑞利数r的增大而增大;且腔体高度d越大,其增长率反而越小.     

格拉晓夫数Gr对侧向局部加热腔体内对流结构的影响

通过二维流体力学基本方程组数值模拟,研究了Pr=6.949时侧向局部加热腔体内不同格拉晓夫数Gr对对流结构的影响.发现随着格拉晓夫数Gr增加出现以下四种情况:准定常流→单局部周期滚动→双局部周期滚动→定常流.还分析了周期对流阶段的对流特性.随着格拉晓夫数Gr增大,对流卷随时间移动的周期减小.随Gr增大,在0≤≤Y 5.2和0.5≤≤Y 5.7范围内对流卷逐渐消失,在5.2≤≤Y5和5.7≤≤Y10范围内对流卷由单对流卷变为多对流卷.     

摆动行波的时空结构

基于流体力学方程组,对长高比Γ=30腔体内混合流体对流中摆动行波的时空结构进行了数值模拟。结果发现:当分离比Ψ=-0.6,-0.4时,在摆动行波存在的下临界附近,摆动行波的对流滚动有消失也有产生,对流平均波数在周期内变化;在摆动行波存在的上临界附近,摆动行波的对流滚动既无消失也无产生,对流平均波数保持为常数。随着相对瑞利数r的增加,对流滚动的摆动幅度和摆动周期明显减小。当分离比Ψ=-0.6时,摆动行波是准周期的;当分离比Ψ=-0.4时,摆动行波是周期的;当分离比Ψ=-0.2时,在摆动行波存在的下临界附近,观察到了一种沿着行波的对称轴线两侧发生对称的摆动行波的新型对流结构。在摆动行波存在的上临界附近,摆动行波是无周期的。随着分离比负值减小,摆动行波存在的上、下限下移,摆动行波存在的稳定区间Δr减小。     

侧向局部加热对流斑图的演化

通过二维流体力学方程组的数值模拟,探讨了普朗特数Pr(28)703.0时高宽比A(28)10矩形腔体内侧向局部加热自然对流在初始小扰动下的成长情况和格拉晓夫数(Gr)对对流的影响规律.研究表明,流速场线性成长阶段扰动的成长率mγ是格拉晓夫数的函数,并给出了成长率随着格拉晓夫数的变化关系式.随着Gr数的增加,对流斑图按定常流—周期—准周期变化,在周期性阶段观察到了2种新的对流斑图,即相位错位对流斑图和同相位对流斑图.最后,探讨了腔体尺寸(d)对对流振幅和热壁面传热能力的影响.     

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

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

瑞利数对行波对流缺陷特性的影响

通过流体力学方程的数值模拟,研究了瑞利数对分离比ψ=-0.6的混合流体行波对流缺陷结构的影响。结果表明:对于给定的相对瑞利数r,缺陷的出现是间歇性的,缺陷出现的位置固定,缺陷出现的周期保持为常数;对于不同的r,缺陷形成时滚动的分裂是不对称的,可以在原滚动的上方也可以在原滚动的下方形成一个新的滚动,从而形成缺陷,缺陷出现的位置基本稳定在腔体的中部,缺陷出现周期随着r的增加而增加;在具有缺陷的行波存在的下限附近,缺陷出现周期减小得较快;在具有缺陷的行波存在的上限附近,缺陷出现周期增加得较快。有缺陷和无缺陷的行波的对流垂直流速最大值δwmax都随着时间在周期变化,但规律是不同的。具有缺陷的行波中,垂直流速最大值δwmax的周期代表缺陷出现的周期;在无缺陷的行波中,垂直流速最大值δwmax的周期代表行波的周期。     

长矩形截面腔体内具有缺陷的对传行波斑图

在长高比Γ=40、分离比Ψ=-0.6情况下,通过流体力学基本方程组的数值模拟,得到了一种新的有趣斑图,即单侧缺陷摆动对传行波.通过分析单侧缺陷摆动对传行波各物理场随时间变化的等值线图、波形图,讨论了其形成过程及时空动力学特性.进一步对比不同相对瑞利数下的单侧缺陷摆动对传行波,探讨了缺陷数、缺陷周期以及缺陷方向等缺陷特征对相对瑞利数的依赖性.     

具有大负分离比的混合流体局部行进波对流

从流体层底部加热引起的对流运动是研究非平衡对流的时空结构或斑图(Pattern)及非线性动力学特性的典型模型之一.本文通过流体力学基本方程组的数值模拟,探讨了具有强Soret效应的混合流体局部行进波的形成过程,发现当分离比ψ=-0.6时,在局部行进波的存在范围内,向局部行进波过渡的不同过程依赖于相对瑞利数r.进一步,讨论了具有强Soret效应的混合流体局部行进波流速场,温度场, 浓度场的结构和特性,分析了局部行进波的存在区间对分离比ψ的依赖性.发现随着Soret效应的增强或负分离比ψ的绝对值的增加,局部行进波稳定存在的区间Δr也在增加.     

混合流体局部行波对流斑图选择的初值依赖性

利用Simple算法对流体力学基本方程组进行了数值模拟,初步研究了局部行波对流斑图选择的初值依赖性问题。分离比ψ-(28)6.0、相对瑞利数r(28)2.1时依赖于初值的有间歇性缺陷的行波,位于腔体右端的局部行波和位于腔体左端的局部行波的多重稳定性;分离比ψ(28)-0.6、相对瑞利数r在1.855~2.118范围内依赖于初值的位于腔体右端的局部行波和位于腔体左端的局部行波的多重稳定性等。虽然在不同初值下,局部行波存在的区间有所不同,局部行波的空间位置有所不同,但局部行波的特性参数变化规律基本一致。结果说明混合流体局部行波对流斑图选择的初值依赖性是存在的。