缺陷源摆动的对传波

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

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

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

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

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

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

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

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

利用二维数值分析,探讨了长高比Γ=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效应的混合流体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的增大,摆动幅度减小,摆动周期变小;随着长高比Γ的增加,摆动行波存在的范围增大;相同相对瑞利数情况下,长高比Γ较大的腔体,摆动行波存在的周期较大。     

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

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

混合流体对流中行波斑图的三重稳定性

本文利用Simple算法对流体力学基本方程组进行了数值模拟,探讨了矩形腔体混合流体中对流行波斑图的多重稳定性问题.当分离比Ψ=-0.4,r=1.8时,首次发现存在三种不同的依赖于初值的稳定对流斑图,即缺陷源摆动的对传波,向左传播的有缺陷行波和向右传播的有缺陷行波的三重稳定性.两种有缺陷的行波传播方向不同,但缺陷出现周期及最大振幅随着时间的变化规律一致;缺陷源摆动的对传波与有缺陷的行波的斑图不同,最大振幅随着时间的变化规律差别较大.结果说明在本文参数情况下混合流体对流行波斑图存在三重稳定性或者初值依赖性.     

摆动行波的时空结构

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

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

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

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

局部行波对流的时空结构

通过流体力学方程的数值模拟,研究了矩形腔体长高比Γ(28)20和混合流体分离比ψ(28)-0.4的无缺陷的局部行波对流和有缺陷的局部行波对流的特性。结果表明:局部行波可以存在于腔体左端壁附近,也可以存在于右端壁附近,但传播方向始终在背离端壁面的方向,对流存在区间长度随着相对瑞利数r的增加而增加;沿着对流分叉曲线增加r,当r超过临界值后,系统从无缺陷的局部行波过渡到有缺陷的局部行波状态;对于流体参数,无缺陷的局部行波存在于r∈[1.519,1.600],有缺陷的局部行波存在于r∈(1.600,1.670];有缺陷的局部行波中对流区间长度随着r增加基本保持为常数;缺陷出现在腔体中的相对位置基本不变,距背离传播方向的端壁5.8左右,缺陷出现的周期随着r增加呈增加趋势,腔体的长高比足够大是出现有缺陷的局部行波的一个必要条件。     

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

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

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

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

无黏流体与固体界面的漏瑞利波

漏瑞利波存在于半无限无黏性流体和半无限固体媒质的界面处. 首先推导流固无限各向同性介质界面处漏瑞利波的特征方程和位移及应力的解析计算公式. 然后结合典型结构通过数值计算研究了漏瑞利波特性以及位移和应力在流体和固体中的分布规律. 数值计算结果表明漏瑞利波的相速度和衰减随流固密度比的增大而增大, 在流固界面上法向位移连续而切向位移不连续. 流固密度比对固体媒质中沿垂直于漏瑞利波的传播方向的位移、正应力和剪应力有比较大的影响,而对沿漏瑞利波的传播方向的正应力几乎没影响. 为利用漏瑞利波的无损检测与评价提供了理论基础.      

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

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

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

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

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

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

腔体高度对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越大,其增长率反而越小.