Loading [MathJax]/jax/element/mml/optable/SuppMathOperators.js

搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

倾斜Poiseuille-Rayleigh-Bénard流动的对流分区与动力学特性

宁利中 张珂 宁碧波 刘爽 田伟利

宁利中, 张珂, 宁碧波, 刘爽, 田伟利. 倾斜Poiseuille-Rayleigh-Bénard流动的对流分区与动力学特性. 物理学报, 2020, 69(12): 124401. doi: 10.7498/aps.69.20191941
引用本文: 宁利中, 张珂, 宁碧波, 刘爽, 田伟利. 倾斜Poiseuille-Rayleigh-Bénard流动的对流分区与动力学特性. 物理学报, 2020, 69(12): 124401. doi: 10.7498/aps.69.20191941
Ning Li-Zhong, Zhang Ke, Ning Bi-Bo, Liu Shuang, Tian Wei-Li. Convection partition and dynamics in inclined Poiseuille-Rayleigh-Bénard flow. Acta Phys. Sin., 2020, 69(12): 124401. doi: 10.7498/aps.69.20191941
Citation: Ning Li-Zhong, Zhang Ke, Ning Bi-Bo, Liu Shuang, Tian Wei-Li. Convection partition and dynamics in inclined Poiseuille-Rayleigh-Bénard flow. Acta Phys. Sin., 2020, 69(12): 124401. doi: 10.7498/aps.69.20191941

倾斜Poiseuille-Rayleigh-Bénard流动的对流分区与动力学特性

宁利中, 张珂, 宁碧波, 刘爽, 田伟利

Convection partition and dynamics in inclined Poiseuille-Rayleigh-Bénard flow

Ning Li-Zhong, Zhang Ke, Ning Bi-Bo, Liu Shuang, Tian Wei-Li
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 利用流体力学基本方程组的数值模拟, 探讨了具有通过流动的倾斜腔体中普朗特数Pr = 6.99的流体的对流分区与动力学特性. 结果表明, 对于相对瑞利数r = 9, 在通过流动雷诺数Re = 1.5时, 随着腔体倾斜角θ的增加, 系统出现均匀行波对流、非均匀行波对流以及单对流圈斑图; 在通过流动雷诺数Re = 12.5时, 随着腔体倾斜角θ的增加, 系统出现局部行波对流、平行流及局部单对流圈斑图; 进一步, 对通过流动雷诺数Re和腔体倾斜角θ的不同组合情况的数值模拟, 发现在通过流动雷诺数Re和腔体倾斜角θ构成的平面上, 具有通过流动的倾斜腔体中的对流可以分成前述六种斑图区域, 即均匀行波对流区、非均匀行波对流区、单对流圈区、局部行波对流区、平行流区及局部单对流圈区. 研究了不同对流区域对流最大垂直流速wmax和努塞尔数Nu随着时间的变化特性. 探讨了不同对流区域对流振幅A和努塞尔数Nu随着腔体倾斜角θ变化的动力学特性.
    Since Henri Bénard first carried out experiments on convection in the fluid layer heated from below at the beginning of last century, and Lord Rayleigh first analyzed small disturbance theoretically, Rayleigh-Bénard convection has received much attention from many researchers, and has become one of the models to study the spatiotemporal structure, flow stability and dynamic characteristics of convection. The methods of studying the Rayleigh-Bénard convection are divided into experimental research, theoretical analysis and numerical simulation. With the development of computer, the research of numerical simulation has made great progress. Because the Rayleigh-Bénard convection can be accurately described by continuity equation, momentum equation and energy equation of hydrodynamics. Therefore, the numerical simulation based on hydrodynamics equations has aroused a lot of research interest. Based on the classical Rayleigh-Bénard convection, the influence of horizontal flow on the Rayleigh-Bénard convection can be studied by applying horizontal flow to one end of the horizontal cavity. On the other hand, the influence of cavity inclination on Rayleigh-Bénard convection can be studied by considering the variation of inclined angles in the cavity. Some valuable convective properties have been obtained. In order to reveal some new convection structures or phenomena, the effects of cavity inclination and through-flow on Rayleigh-Bénard convection are considered at the same time in this paper.By using the numerical simulation of the basic equations of hydrodynamics, the convection partition and dynamic characteristics of the fluid with Prandtl number Pr = 6.99 in the inclined cavity with through-flows are discussed. The results show that for the reduced Rayleigh number r = 9, the system presents uniform traveling wave convection, non-uniform traveling wave convection and single roll convection pattern at the through-flow Reynolds number Re = 1.5 with the increase of the inclined angle θ in the cavity, that for the through-flow Reynolds number Re = 12.5, the system presents the localized traveling wave convection, parallel flow and localized single roll convection pattern with the increase of the inclined angle θ in the cavity, that furthermore, the numerical simulation of different values of through-flow Reynolds number Re and inclined angle θ in the cavity shows that on the plane composed of through flow Reynolds number Re and inclined angle θ in the cavity, the convection in the inclined cavity with through-flow can be divided into six kinds of pattern regions, namely, uniform traveling wave convection region, non-uniform traveling wave convection region, single roll convection region, localized traveling wave convection region, parallel flow region, and localized single roll convection region. The characteristics of the maximum vertical velocity wmax and Nusselt number Nu of convection varying with time in different convection regions are studied. The dynamic properties of convective amplitude A and Nusselt number Nu in different convective regions varying with inclined angle θ in the cavity are discussed.
      PACS:
      44.25.+f(Natural convection)
      44.15.+a(Channel and internal heat flow)
      05.70.Ln(Nonequilibrium and irreversible thermodynamics)
      通信作者: 宁碧波, ning_bibo@163.com
    • 基金项目: 国家级-国家自然科学基金项目(10872164)
      Corresponding author: Ning Bi-Bo, ning_bibo@163.com

    自从20世纪初Henri Bénard首先对底部加热流体层对流进行实验, Lord Rayleigh首先进行小扰动理论分析以来, Rayleigh-Bénard对流问题一直受到许多研究者的广泛关注[1, 2]. 文献[1, 2]就Rayleigh-Bénard对流问题的研究进展进行过综述. 关于Rayleigh-Bénard对流问题的研究方法分为实验研究、理论分析以及数值模拟. 最初, 人们通过实验及理论分析的方法进行研究. 随着计算机的发展, 数值模拟分析研究获得了很大进展. 基于对流振幅方程组的模拟具有计算时间少又能定性揭示对流特性的特点[3,4]. 由于Rayleigh-Bénard对流现象可以精确地由流体力学的连续方程、动量方程及能量方程来描述, 因此, 基于流体力学方程组的数值模拟引起了广泛的研究兴趣[5-16]. 已经揭示了行波[5-7]、局部行波[8-12]、振动对流[9]、行波中的缺陷结构[6,7,12,13]、摆动行波[14,15]及对传波[16]等对流结构或者现象, 取得了丰富的成果. 在对经典Rayleigh-Bénard对流问题研究的基础上, 对水平腔体的一侧施加水平流动[17-19], 研究水平流动对Rayleigh-Bénard对流问题的影响. 已经发现了对流发生的附近的线性特性[17], 局部行波的周期性[18], 对流斑图的分区特性及斑图成长特性[19]等一些新的对流结构. 如果对于经典Rayleigh-Bénard对流问题, 考虑腔体倾斜时[20-25], 研究获得了一些新的对流特性, 如缺陷湍流, 对流层的双稳定性与竞争, 不同普朗特数情况下的对流斑图、分区及临界条件等. 本文进一步在考虑腔体倾斜的同时, 也考虑通过流动对Rayleigh-Bénard对流问题的影响, 试图获得一些新的对流结构或者对流现象.

    本文利用流体力学基本方程组的数值模拟, 探讨了具有通过流动的倾斜腔体中普朗特数Pr = 6.99的流体的对流分区与动力学特性. 对不同通过流动雷诺数Re和腔体倾斜角θ的组合情况的数值模拟, 发现在通过流动雷诺数Re和腔体倾斜角θ构成的平面上, 具有通过流动的倾斜腔体中的对流可以分成六种斑图区域, 即均匀行波对流区、非均匀行波对流区、单对流圈区、局部行波对流区、平行流区及局部单对流圈区. 这是本文模型不同于以前模型的最大特点. 不同对流区域对流最大垂直流速wmax和努塞尔数Nu随着时间的变化具有不同特性. 不同对流区域对流振幅A和努塞尔数Nu随着腔体倾斜角θ的变化也具有不同动力学规律.

    对于图1所示的倾斜角为θ, 坐标为x,y的腔体, 当下热壁面与上冷壁面之间的温度差ΔT特征的瑞利数Ra(Ra=gαΔTd3/κν)(α是体积膨胀系数, d是腔体垂直高度, g是重力加速度, κ为热扩散系数κ=λ/(ρ0cp), cp为定压比热容, ρ是流体密度, λ为热传导系数, ν是运动黏性系数. 下标0表示对应物理参数的参考值, 按传导状态下沿腔体垂直高度的平均值计算.)超过相应的临界值时, 对流就发生在腔体中. 对于这一问题, 采用Boussinesq假设, 流体密度ρ随着流体温度T的变化关系式为ρ=ρ0[1α(TT0)]. 这时描述倾斜腔体中流体对流的流体力学基本方程组为

    图 1 倾斜腔体的流动示意图\r\nFig. 1. Flow diagram in inclined cavity.
    图 1  倾斜腔体的流动示意图
    Fig. 1.  Flow diagram in inclined cavity.
    ux+wy=0, (1)
    ut+(u2)x+(uw)y=peffρ0x+ν(2ux2+2uy2)gα(TT0)sinθ, (2)
    wt+(uw)x+(w2)y=peffρ0y+ν(2wx2+2wy2)+gα(TT0)cosθ, (3)
    Tt+(uT)x+(wT)y=κ(2Tx2+2Ty2), (4)

    式中有效压强peff=p+ρ0gycosθρ0gxsinθ, u表示x方向的速度, w表示y方向的速度, T表示温度, p是流体压强, t是时间.

    在腔体中上下壁面都是固定的, 速度u, w为0, 上下壁面等温. 左端进口速度为平面Poiseuille流的速度分布, 温度上下壁面之间线性变化. 右端出口为自由出流. 因此, 上下壁面边界条件为

    y=0时, u=w=0, T=T0+0.5ΔT;

    y = d时, u=w=0, T=T00.5ΔT.

    左端进口条件为

    x=0时, u=4umaxzd(1zd), w = 0.

    x=0时, T=T0+ΔT(0.5zd).

    x=Lx时, 右端出口为自由出流.

    这里umax为进口速度最大值, Lx为计算腔体的长度.

    初始条件流速为u=w=0, 温度为平均温度T0.

    在数值计算中, 利用有限容积法离散了控制力学方程组, Simple算法用于求解速度-压力耦合方程. 扩散项采用中心差分格式, 对流项采用迎风格式, 时间项采用一阶隐式格式离散. 采用均匀交错网格系统. 数值模拟的可靠性和正确性可以由文献[26]的实例证明. 由于文献[27]对d/20d/30的网格进行了计算, 发现计算结果几乎相同, 误差很小. 所以采用d/20的网格进行模拟, 时间步长采用Δt=0.01 s.

    讨论中采用量纲归一化参数描述腔体体形, 坐标为X=x/d, Y=y/d, 长高比为Γ=Lx/d. 对流系统可以由瑞利数Ra(为了方便, 使用相对瑞利数r=Ra/Rac, 其中Rac=1708), 普朗特数Pr=ν/κ, 表征进口流动的雷诺数Re=umaxdν及倾斜腔体的倾角θ来控制.

    本文采用T0=293.15K, Pr=6.99的流体. 长高比Γ=20.

    对于具有通过流动的倾斜腔体中Pr=6.99的对流, 当给定相对瑞利数r=9时, 不同的通过流动雷诺数和不同的倾斜角度组合时, 系统出现六种不同的对流斑图.

    3.1.1   均匀行波对流

    Re=1.5并且θ=10时, 出现均匀行波对流结构, 如图2所示. 图2(a)为腔体中的温度场, 可以看出, 温度场上升羽流和下降羽流间隔出现, 羽流以垂直线为对称; 在腔体长度方向上均匀分布. 图2(b)为腔体中的流线分布, 在向上弯曲的通过流线的下部或者向下弯曲的通过流线的上部出现封闭的流线圈, 流线圈均匀分布. 温度场下降羽流的位置和两个流线圈之间的下降流位置相对应; 温度场上升羽流的位置和两个流线圈之间的上升流位置相对应. 腔体中共有23个流线圈, 它的波数为 k=23π/20=3.61. 图2(c)为腔体中011范围内的速度矢量场. 在向上弯曲的通过流线的下部的矢量场即下壁面附近的矢量场是顺时针方向; 在向下弯曲的通过流线的上部的矢量场即上壁面附近的矢量场是逆时针方向. 这一对具有不同旋转方向的对流圈沿腔体长度方向上交替出现. 在对流圈的左侧如果是箭头向下的下降流, 在对流圈的右侧必然是箭头向上的上升流. 沿腔体长度方向上, 可以清楚地看出, 上升流和下降流交替地出现. 在上升流和下降流处垂直速度较大, 对流圈的上下部水平速度较大. 对流圈的中心部位速度较小. 由图2可以看出, 温度场、流线及速度矢量场都是均匀分布并整体上向一个方向流动, 类似于Rayleigh-Bénard-Poiseuille 流动中出现的行波, 因此称其为均匀行波.

    图 2 ${{Re}} = 1.5$且${{\theta }} = 10^\circ $时的均匀行波对流 (a) 温度; (b)流线; (c)速度矢量\r\nFig. 2. Uniform traveling wave convection at ${{Re}} = 1.5$ and ${{\theta }} = 10^\circ $: (a) Temperature; (b) streamlines; (c) velocity vector
    图 2  {{Re}} = 1.5{{\theta }} = 10^\circ 时的均匀行波对流 (a) 温度; (b)流线; (c)速度矢量
    Fig. 2.  Uniform traveling wave convection at {{Re}} = 1.5 and {{\theta }} = 10^\circ : (a) Temperature; (b) streamlines; (c) velocity vector
    3.1.2   非均匀行波对流

    保持雷诺数为{{Re}} = 1.5, 增加倾斜角到{{\theta }} = 30^\circ 时, 系统出现非均匀行波对流结构, 如图3所示. 由图3(a)的温度场可知, 温度场上升羽流和下降羽流间隔出现, 但在某些部位的羽流不再关于垂直线对称, 在腔体长度方向上为非均匀分布. 图3(b)为腔体中的流线分布, 可以看出, 在向上弯曲的通过流线的下部或者向下弯曲的通过流线的上部出现封闭的流线圈. 在向上弯曲的通过流线的下部出现封闭的流线圈在某些部位被拉长, 变大. 流线圈在空间上的分布不均匀, 没有规律. 由于某些流线圈的变大, 腔体空间上对流圈数减少到13个, 它的空间平均波数也随之减小, 变为{{k}} = 13{\text{π}}/20 = 2.04. 图3(c)为腔体中0 \leqslant {{X}} \leqslant 13范围内的速度矢量场. 可以看出, 顺时针方向旋转的对流圈和逆时针方向旋转的对流圈沿腔体长度方向交替出现, 上升流和下降流也沿腔体长度方向交替出现. 但对流圈的长度分布不均匀, 大的对流圈是小对流圈的2—3倍长. 由图3可以看出, 温度场、流线及速度矢量场都是非均匀分布但整体上向一个方向流动, 因此称其为非均匀行波.

    图 3 ${{Re}} = 1.5$和${{\theta }} = 30^\circ $时的非均匀行波对流 (a) 温度; (b)流线; (c)速度矢量\r\nFig. 3. Non-uniform traveling wave convection at ${{Re}} = 1.5$ and ${{\theta }} = 30^\circ $: (a) Temperature; (b) streamlines; (c) velocity vector
    图 3  {{Re}} = 1.5{{\theta }} = 30^\circ 时的非均匀行波对流 (a) 温度; (b)流线; (c)速度矢量
    Fig. 3.  Non-uniform traveling wave convection at {{Re}} = 1.5 and {{\theta }} = 30^\circ : (a) Temperature; (b) streamlines; (c) velocity vector
    3.1.3   单对流圈斑图

    保持雷诺数为{{Re}} = 1.5, 继续增加倾斜角到{{\theta }} = 60^\circ 时, 系统出现单对流圈斑图, 如图4所示. 图4(a)的温度场除腔体的两端附近外, 是近似于平行上下壁面的直线. 图4(b)的流线是由光滑曲线构成的一个大对流圈. 图4(c)的速度矢量图是一个顺时针方向旋转的对流圈. 速度矢量显示在对流圈的外周速度较大, 对流圈的中心速度较小. 由于图4(b)图4(c)的腔体中具有一个对流圈, 因此, 称其为单对流圈.

    图 4 ${{Re}} = 1.5$和${{\theta }} = 60^\circ $时的单对流圈 (a) 温度; (b)流线; (c)速度矢量\r\nFig. 4. Single roll convection at ${{Re}} = 1.5$ and ${{\theta }} = 60^\circ $: (a) Temperature; (b) streamlines; (c) velocity vector.
    图 4  {{Re}} = 1.5{{\theta }} = 60^\circ 时的单对流圈 (a) 温度; (b)流线; (c)速度矢量
    Fig. 4.  Single roll convection at {{Re}} = 1.5 and {{\theta }} = 60^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.
    3.1.4   局部行波对流斑图

    当雷诺数增加到{{Re}} = 12.5, 倾斜角减小到{{\theta }} = 1^\circ 时, 系统出现局部行波对流斑图, 如图4所示. 可以看出, 在腔体左端部0 \leqslant {{X}} \leqslant 16的范围内, 温度等值线和流线是一些平行于腔体上下壁面的直线. 仅存在平行的通过流动. 在腔体右端部16 \leqslant {{X}} \leqslant 20的范围内, 温度等值线形成了上升羽和下降羽, 但沿垂直线不对称, 流线形成了上下弯曲的流动, 在向上弯曲和在向下弯曲的部位垂直尺度较小的对流圈. 10 \leqslant {{X}} \leqslant 20范围内的速度矢量如图5(c)所示, 可以看出, 左端部附近平行流和右端部附近对流结构的矢量场的特性. 由于平行的通过流动和对流共存于腔体中, 故称之为局部行波对流斑图, 类似于Rayleigh-Bénard-Poiseuille 流动中出现的局部行波.

    图 5 ${{Re}} = 12.5$和${{\theta }} = 1^\circ $时的局部行波对流 (a) 温度; (b)流线; (c)速度矢量\r\nFig. 5. Localized traveling wave convection at ${{Re}} = 12.5$ and ${{\theta }} = 1^\circ $: (a) Temperature; (b) streamlines; (c) velocity vector.
    图 5  {{Re}} = 12.5{{\theta }} = 1^\circ 时的局部行波对流 (a) 温度; (b)流线; (c)速度矢量
    Fig. 5.  Localized traveling wave convection at {{Re}} = 12.5 and {{\theta }} = 1^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.
    3.1.5   平行流斑图

    保持雷诺数为{{Re}} = 12.5, 增加倾斜角到{{\theta }} = 10^\circ 时, 系统出现平行流斑图, 如图6所示. 可以看出, 0 \leqslant {{X}} \leqslant 20范围内的温度等值线和流线与10 \leqslant {{X}} \leqslant 15范围内的速度矢量图都是一些平行于上下壁面的直线, 故称其为平行流.

    图 6 ${{Re}} = 12.5$和${{\theta }} = 10^\circ $时的平行流动 (a) 温度; (b)流线; (c)速度矢量\r\nFig. 6. Parallel flows at ${{Re}} = 12.5$ and ${{\theta }} = 10^\circ $: (a) Temperature; (b) streamlines; (c) velocity vector.
    图 6  {{Re}} = 12.5{{\theta }} = 10^\circ 时的平行流动 (a) 温度; (b)流线; (c)速度矢量
    Fig. 6.  Parallel flows at {{Re}} = 12.5 and {{\theta }} = 10^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.
    3.1.6   局部单对流圈斑图

    保持雷诺数为{{Re}} = 12.5, 增加倾斜角到{{\theta }} = 45^\circ 时, 系统出现局部单对流圈斑图, 如图7所示. 可以看出, 图7(a)的温度等值线和图7(b)的流线在腔体的上部是一些平行于上下壁面的直线. 流线在腔体的下部是一个单对流圈. 从图7(c)的速度矢量图可以更明显地看出, 腔体的上部是平行上壁面的流动, 下部是一个顺时针的滚动. 将这种上部具有平行流并且下部具有一个滚动的流动斑图称为局部单对流圈斑图.

    图 7 $R{{e}} = 12.5$和${{\theta }} = 45^\circ $时的局部单对流圈 (a) 温度; (b)流线; (c)速度矢量\r\nFig. 7. Localized single roll convection at $R{{e}} = 12.5$ and ${{\theta }} = 45^\circ $: (a) Temperature; (b) streamlines; (c) velocity vector.
    图 7  R{{e}} = 12.5{{\theta }} = 45^\circ 时的局部单对流圈 (a) 温度; (b)流线; (c)速度矢量
    Fig. 7.  Localized single roll convection at R{{e}} = 12.5 and {{\theta }} = 45^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.
    3.1.7   对流垂直最大流速随着时间的变化特性

    图8(a)是相对瑞利数{{r}} = 9, 雷诺数R{{e}} = 1.5及倾斜角分别为{{\theta }} = 10^\circ , 30^\circ , 60^\circ 时对应的均匀行波对流、非均匀行波对流以及单对流圈斑图的最大垂直流速{{{w}}_{\max}}随着时间的变化. 可以看出, 当{{\theta }} = 10^\circ 30^\circ 时对应的均匀行波对流、非均匀行波对流的最大垂直流速{{{w}}_{\max}}随着时间增加而增长, 在t = 150 s左右达到稳定, 然后不再随着时间变化. {{\theta }} = 60^\circ 时对应的单对流圈斑图的最大垂直流速{{{w}}_{\max}}, 首先随着时间增加而增长, 然后随着时间增加而上下振荡, 在t = 200 s以后达到稳定, 不再随着时间变化. 显然, 后者稳定下来所需要时间较前者长. 最大垂直流速{{{w}}_{\max}}随着倾斜角的增加而变大.

    图8(b)是相对瑞利数{{r}} = 9, 雷诺数R{{e}} = 12.5及倾斜角分别为{{\theta }} = 1^\circ , 10^\circ , 45^\circ 时对应的局部行波对流、平行流及局部单对流圈斑图的最大垂直流速{{{w}}_{\max}}随着时间的变化. 可以看出, 当{{\theta }} = 45^\circ 时局部单对流圈斑图的最大垂直流速{{{w}}_{\max}}随着时间增加而增长, 在t = 110 s左右达到稳定, 然后不再随着时间变化. {{\theta }} = 1^\circ , 10^\circ 时对应的局部行波对流与平行流的最大垂直流速{{{w}}_{\max}}, 首先随着时间增加而增长, 然后随着时间增加而上下振荡, 不再增加. 在t = 340 s以后{{\theta }} = 10^\circ 对应的平行流的最大垂直流速{{{w}}_{\max}}稳定下来; 在t = 290 s以后达到{{\theta }} = 1^\circ 对应的平行流的最大垂直流速{{{w}}_{\max}}达到稳定. 最大垂直流速{{{w}}_{\max}}随着倾斜角的增加而变大. 比较图8(a)图8(b)可以明显地看出, 在大倾斜角情况下, 最大垂直流速{{{w}}_{\max}}随着时间单调增加后达到稳定; 在小倾斜角情况下, 最大垂直流速{{{w}}_{\max}}随着时间单调增加后变成上下振荡, 然后达到稳定. 大雷诺数情况下最大垂直流速{{{w}}_{\max}}上下振荡的时间较长. 当系统的参数和初始条件在某种组合条件下, 系统出现瞬态振荡现象.

    图 8 不同对流结构的最大垂直流速的时间演化 (a) Re = 1.5; (b) Re = 12.5\r\nFig. 8. Time evolution of maximum vertical velocity in different convection structures: (a) Re = 1.5; (b) Re = 12.5.
    图 8  不同对流结构的最大垂直流速的时间演化 (a) Re = 1.5; (b) Re = 12.5
    Fig. 8.  Time evolution of maximum vertical velocity in different convection structures: (a) Re = 1.5; (b) Re = 12.5.
    3.1.8   努塞尔数Nu随着时间的变化特性

    努塞尔数Nu是表示对流换热强烈程度的一个无量纲数, 定义为Nu = {{hd}}/{\lambda }, 其中h为流体的对流传热系数. 图9(a)是相对瑞利数{{r}} = 9, 雷诺数R{{e}} = 1.5及倾斜角分别为{{\theta }} = 10^\circ , 30^\circ , 60^\circ 时对应的均匀行波对流、非均匀行波对流以及单对流圈斑图的努塞尔数{{Nu}}随着时间的变化. 可以看出, 当{{\theta }} = 10^\circ 时对应的均匀行波对流努塞尔数{{Nu}}随着时间增加而减小, 然后出现振荡, 在t = 300 s左右达到稳定. 当{{\theta }} = 30^\circ 时对应的非均匀行波对流努塞尔数{{Nu}}随着时间增加而减小, 在t = 240 s左右达到稳定. {{\theta }} = 60^\circ 时对应的单对流圈斑图的努塞尔数{{Nu}}首先随着时间增加而减小, 然后随着时间增加而增加, 在t = 190 s以后达到稳定. 努塞尔数{{Nu}}随着倾斜角的增加而变大.

    图 9 不同对流结构的努塞尔数的时间演化 (a) Re = 1.5; (b) Re = 12.5\r\nFig. 9. Time evolution of Nusselt number in different convection structures: (a) Re = 1.5; (b) Re = 12.5.
    图 9  不同对流结构的努塞尔数的时间演化 (a) Re = 1.5; (b) Re = 12.5
    Fig. 9.  Time evolution of Nusselt number in different convection structures: (a) Re = 1.5; (b) Re = 12.5.

    图9(b)是相对瑞利数{{r}} = 9, 雷诺数R{{e}} = 12.5及倾斜角分别为{{\theta }} = 1^\circ , 10^\circ , 45^\circ 时对应的局部行波对流、平行流及局部单对流圈斑图的努塞尔数{{Nu}}随着时间的变化. 三种情况下, 努塞尔数{{Nu}}随着时间增加而减小, 在t = 10 s左右达到稳定. 与小雷诺数情况下相反, 努塞尔数{{Nu}}随着倾斜角的增加而减小.

    对前述的六种对流斑图进行了大量的数值模拟后, 发现当{{r}} = 9时在倾斜角\theta 与雷诺数{{Re}} 平面上可以划分为6个不同的区域, 如图10所示. 对于较小的雷诺数{{Re}} , 当倾斜角\theta 较小时, 由于倾斜角\theta 对对流的影响微乎其微, 在腔体长度方向冷面上的重力分力和热面上浮力都很小, 可以忽略, 系统出现的是类似水平腔体情况下的均匀行波对流; 随着倾斜角\theta 的增加, 在腔体长度方向冷面上的重力分力和热面上浮力都已经发挥作用, 随着两个方向浮力作用的竞争, 系统中的均匀行波对流转变成非均匀行波对流; 随着倾斜角\theta 的进一步增加, 在腔体长度方向冷面上的重力分力和热面上浮力明显加强, 起到控制作用, 系统过渡到单对流圈斑图. 对于雷诺数{{Re}} = 10附近的情况, 当倾斜角\theta 较小时, 由于倾斜角\theta 对对流的影响微不足道, 在腔体长度方向冷面上的重力分力和热面上浮力很小, 其影响可以忽略, 系统出现的是类似水平腔体情况下的局部行波对流; 随着倾斜角\theta 的增加, 在腔体长度方向冷面上的重力分力和热面上浮力变得明显, 与雷诺数{{Re}} 的综合作用结果使系统转化成局部单对流圈斑图. 对于雷诺数Re \approx 15的情况, 雷诺数的影响明显, 起了控制作用, 对于倾斜角θ约小于{{50}}^\circ 的情况, 系统出现的都是平行流斑图. 在倾斜角θ约大于{{50}}^\circ 时, 形成平行流动的临界雷诺数几乎线性变大.

    图 10 不同对流斑图在倾斜角度$\theta $-雷诺数${{Re}}$平面上的分区\r\nFig. 10. Partition of different convective patterns in the plane $\theta $-${{Re}}$.
    图 10  不同对流斑图在倾斜角度\theta -雷诺数{{Re}}平面上的分区
    Fig. 10.  Partition of different convective patterns in the plane \theta -{{Re}}.
    3.3.1   对流振幅的特性

    {{r}} = 9, 考虑{{Re}} = 1.5{{Re}} = 12.5时不同倾斜角度情况下的对流特征物理量. 对最大垂直流速进行量纲归一化处理, 得到对流最大振幅, 即A = {w_{\max}}/\left({\kappa /d} \right). 图11是最大振幅A随着倾斜角的变化. 可以看出, 图11(a)中的均匀行波对流、非均匀行波对流以及单对流圈斑图的最大振幅A随着倾斜角的增加而增加. 图11(b)中的局部行波对流、平行流及局部单对流圈斑图的最大振幅A也随着倾斜角的增加而增加. 相应倾斜角情况下, 图11(b)中的最大振幅A大于图11(a)中的最大振幅A.

    图 11 不同对流结构时A随着${{\theta }}$的变化 (a) Re = 1.5; (b) Re = 12.5\r\nFig. 11. Variation of A in different convection structures with ${{\theta }}$: (a) Re = 1.5; (b) Re = 12.5.
    图 11  不同对流结构时A随着{{\theta }}的变化 (a) Re = 1.5; (b) Re = 12.5
    Fig. 11.  Variation of A in different convection structures with {{\theta }}: (a) Re = 1.5; (b) Re = 12.5.
    3.3.2   努塞尔数的特性

    图12是努塞尔数{{Nu}}随着倾斜角的变化. 从图12(a)可以看出, 均匀行波对流与非均匀行波对流时努塞尔数{{Nu}}随着倾斜角的变化几乎不发生变化, 是一个常数. 对于单对流圈斑图, 努塞尔数{{Nu}}随着倾斜角的增加而增加. 图12(b)中的局部行波对流、平行流及局部单对流圈斑图的努塞尔数{{Nu}}随着倾斜角的增加而减小. 在不同倾斜角情况下, 图12(b)中的努塞尔数{{Nu}}变化范围与图12(a)中的努塞尔数{{Nu}}类似.

    图 12 不同对流结构时${{Nu}}$随着${{\theta }}$的变化 (a) Re = 1.5; (b) Re = 12.5\r\nFig. 12. Variation of ${{Nu}}$ in different convection structures with ${{\theta }}$: (a) Re = 1.5; (b) Re = 12.5.
    图 12  不同对流结构时{{Nu}}随着{{\theta }}的变化 (a) Re = 1.5; (b) Re = 12.5
    Fig. 12.  Variation of {{Nu}} in different convection structures with {{\theta }}: (a) Re = 1.5; (b) Re = 12.5.

    利用流体力学基本方程组的数值模拟, 探讨了具有通过流动的倾斜腔体中普朗特数Pr = 6.99的流体的对流分区与动力学特性. 可以得出以下结论.

    1)相对瑞利数r = 9时, 对通过流动雷诺数Re和腔体倾斜角\theta 的不同组合情况进行了数值模拟, 发现在通过流动雷诺数Re和腔体倾斜角\theta 构成的平面上, 具有通过流动的倾斜腔体中的对流可以分成六种斑图区域, 即均匀行波对流区、非均匀行波对流区、单对流圈区、局部行波对流区、平行流区及局部单对流圈区.

    2)最大垂直流速{{{w}}_{\max}}随着时间的变化特性说明, 在大倾斜角情况下, 最大垂直流速{{{w}}_{\max}}随着时间单调增加后达到稳定; 在小倾斜角情况下, 最大垂直流速{{{w}}_{\max}}随着时间单调增加后变成上下振荡, 然后达到稳定. 大雷诺数情况下最大垂直流速{{{w}}_{\max}}上下振荡的时间较长.

    3)最大振幅A随着倾斜角的增加而增加, 不同对流斑图区变化规律不同. 努塞尔数{{Nu}}随着倾斜角\theta 的变化特性说明, 均匀行波对流与非均匀行波对流时努塞尔数{{Nu}}随着倾斜角的变化几乎不发生变化. 单对流圈斑图, 努塞尔数{{Nu}}随着倾斜角的增加而增加. 局部行波对流、平行流及局部单对流圈斑图的努塞尔数{{Nu}}随着倾斜角的增加而减小.

    [1]

    Cross M C, Hohenberg P C 1993 Rev. Mod. Phys. 65 851

    [2]

    Bodenschatz E, Pesch W, Ahlers G 2000 Annu. Rev. Fluid Mech. 32 709Google Scholar

    [3]

    Watanabe T, Iima M, Nishiura Y 2012 J. Fluid Mech. 712 219Google Scholar

    [4]

    Mercader I, Batiste O, Alonso A,Knobloch E 2010 Fluid Dyn. Res. 42 025505Google Scholar

    [5]

    Mercader I, Batiste O, Alonso A, Knobloch E 2013 J. Fluid Mech. 722 240Google Scholar

    [6]

    宁利中, 余荔, 袁喆, 周洋 2009 中国科学: 物理学 力学 天文学 39 746

    Ning L Z, Yu L, Yuan Z, Zhou Y 2009 Sci. Sin.: Phys. Mech. Astron. 39 746

    [7]

    宁利中, 王娜, 袁喆, 李开继, 王卓运 2014 物理学报 63 104401Google Scholar

    Ning L Z, Wang N, Yuan Z, Li K J, Wang Z Y 2014 Acta Phys. Sin. 63 104401Google Scholar

    [8]

    Barten W, Lucke M, Kamps M, Schmitz R 1995 Phys. Rev. E 51 5662Google Scholar

    [9]

    Jung D, Lucke M 2002 Phys. Rev. Lett. 89 054502Google Scholar

    [10]

    Taraut A V, Smorodin B L, Lucke M 2012 New J. Phys. 14 093055Google Scholar

    [11]

    Mercader I, Batiste O, Alonso A, Knobloch E 2011 J. Fluid Mech. 667 586Google Scholar

    [12]

    宁利中, 王永起, 袁喆, 李开继, 胡彪 2016 科学通报 61 872Google Scholar

    Ning L Z, Wang Y Q, Yuan Z, Li K, Hu B 2016 Chin. Sci. Bull. 61 872Google Scholar

    [13]

    宁利中, 齐昕, 周洋, 余荔 2009 物理学报 58 2528Google Scholar

    Ning L Z, Qi X, Zhou Y, Yu L 2009 Acta Phys. Sin. 58 2528Google Scholar

    [14]

    Zhao B X, Tian Z F 2015 Phys. Fluids 27 074102Google Scholar

    [15]

    宁利中, 刘爽, 宁碧波, 袁喆, 田伟利, 渠亚伟 2017 应用力学学报 34 1086Google Scholar

    Ning L Z, Liu S, Ning B B, Yuan Z, Tian W L, Qu Y W 2017 Chin. J. Appl. Mech. 34 1086Google Scholar

    [16]

    Ning L Z, Qi X, Yuan Z, Shi F 2008 J. Hydrodyn. 20 567Google Scholar

    [17]

    Jung C, Lucke M, Buchel P 1996 Phys. Rev. E 54 1510Google Scholar

    [18]

    胡彪, 宁利中, 宁碧波, 田伟利 2017 科学通报 62 4278Google Scholar

    Hu B, Ning L Z, Ning B B, Tian W L 2017 Chin. Sci. Bull. 62 4278Google Scholar

    [19]

    宁利中, 胡彪, 宁碧波, 田伟利 2016 物理学报 65 214401Google Scholar

    Ning L Z, Hu B, Ning B B, Tian W L 2016 Acta Phys. Sin. 65 214401Google Scholar

    [20]

    Daniels K E, Bodenschatz E 2002 Phys. Rev. Lett. 88 034501Google Scholar

    [21]

    Daniels K E, Plappb B, Bodenschatz E 2000 Phys. Rev. Lett. 84 5320Google Scholar

    [22]

    Daniels K E, Brausch O, Pesch W, Bodenschatz E 2008 J. Fluid Mech. 597 261Google Scholar

    [23]

    Busse F H, Clever R M 1992 J. Eng. Math 26 1Google Scholar

    [24]

    宁利中, 吴昊, 宁碧波, 田伟利, 宁景昊 2019 应用数学和力学 40 398Google Scholar

    Ning L Z, Wu H, Ning B B, Tian W L, Ning J H 2019 Appl. Math. Mech. 40 398Google Scholar

    [25]

    宁利中, 张珂, 宁碧波, 吴昊, 田伟利 2020 应用力学学报 37 737Google Scholar

    Ning L Z, Zhang K, Ning B B, Wu H, Tian W L 2020 Chin. J. Appl. Mech. 37 737Google Scholar

    [26]

    Patankar S V 1982 Computation of Heat Transfer and Fluid Flow (Minneapolis: University of Minnesota) pp20−30

    [27]

    胡彪, 宁利中, 宁碧波, 田伟利, 吴昊, 宁景昊 2017 水动力学研究与进展 32 336Google Scholar

    Hu B, Ning L Z, Ning B B, Tian W L, Wu H, Ning J H 2017 Chin. J. Hydrodyn. 32 336Google Scholar

    期刊类型引用(2)

    1. 朱玉枝,李友荣. 矩形流道内密度极值流体Poiseuille-Rayleigh-Bénard对流传热研究. 工程热物理学报. 2024(11): 3480-3486 . 百度学术
    2. 宁利中,宁碧波,郝建武,田伟利. 长矩形腔体中混合流体的对称对传波. 应用力学学报. 2021(06): 2392-2397 . 百度学术

    其他类型引用(3)

  • 图 1  倾斜腔体的流动示意图

    Fig. 1.  Flow diagram in inclined cavity.

    图 2  {{Re}} = 1.5{{\theta }} = 10^\circ 时的均匀行波对流 (a) 温度; (b)流线; (c)速度矢量

    Fig. 2.  Uniform traveling wave convection at {{Re}} = 1.5 and {{\theta }} = 10^\circ : (a) Temperature; (b) streamlines; (c) velocity vector

    图 3  {{Re}} = 1.5{{\theta }} = 30^\circ 时的非均匀行波对流 (a) 温度; (b)流线; (c)速度矢量

    Fig. 3.  Non-uniform traveling wave convection at {{Re}} = 1.5 and {{\theta }} = 30^\circ : (a) Temperature; (b) streamlines; (c) velocity vector

    图 4  {{Re}} = 1.5{{\theta }} = 60^\circ 时的单对流圈 (a) 温度; (b)流线; (c)速度矢量

    Fig. 4.  Single roll convection at {{Re}} = 1.5 and {{\theta }} = 60^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.

    图 5  {{Re}} = 12.5{{\theta }} = 1^\circ 时的局部行波对流 (a) 温度; (b)流线; (c)速度矢量

    Fig. 5.  Localized traveling wave convection at {{Re}} = 12.5 and {{\theta }} = 1^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.

    图 6  {{Re}} = 12.5{{\theta }} = 10^\circ 时的平行流动 (a) 温度; (b)流线; (c)速度矢量

    Fig. 6.  Parallel flows at {{Re}} = 12.5 and {{\theta }} = 10^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.

    图 7  R{{e}} = 12.5{{\theta }} = 45^\circ 时的局部单对流圈 (a) 温度; (b)流线; (c)速度矢量

    Fig. 7.  Localized single roll convection at R{{e}} = 12.5 and {{\theta }} = 45^\circ : (a) Temperature; (b) streamlines; (c) velocity vector.

    图 8  不同对流结构的最大垂直流速的时间演化 (a) Re = 1.5; (b) Re = 12.5

    Fig. 8.  Time evolution of maximum vertical velocity in different convection structures: (a) Re = 1.5; (b) Re = 12.5.

    图 9  不同对流结构的努塞尔数的时间演化 (a) Re = 1.5; (b) Re = 12.5

    Fig. 9.  Time evolution of Nusselt number in different convection structures: (a) Re = 1.5; (b) Re = 12.5.

    图 10  不同对流斑图在倾斜角度\theta -雷诺数{{Re}}平面上的分区

    Fig. 10.  Partition of different convective patterns in the plane \theta -{{Re}}.

    图 11  不同对流结构时A随着{{\theta }}的变化 (a) Re = 1.5; (b) Re = 12.5

    Fig. 11.  Variation of A in different convection structures with {{\theta }}: (a) Re = 1.5; (b) Re = 12.5.

    图 12  不同对流结构时{{Nu}}随着{{\theta }}的变化 (a) Re = 1.5; (b) Re = 12.5

    Fig. 12.  Variation of {{Nu}} in different convection structures with {{\theta }}: (a) Re = 1.5; (b) Re = 12.5.

  • [1]

    Cross M C, Hohenberg P C 1993 Rev. Mod. Phys. 65 851

    [2]

    Bodenschatz E, Pesch W, Ahlers G 2000 Annu. Rev. Fluid Mech. 32 709Google Scholar

    [3]

    Watanabe T, Iima M, Nishiura Y 2012 J. Fluid Mech. 712 219Google Scholar

    [4]

    Mercader I, Batiste O, Alonso A,Knobloch E 2010 Fluid Dyn. Res. 42 025505Google Scholar

    [5]

    Mercader I, Batiste O, Alonso A, Knobloch E 2013 J. Fluid Mech. 722 240Google Scholar

    [6]

    宁利中, 余荔, 袁喆, 周洋 2009 中国科学: 物理学 力学 天文学 39 746

    Ning L Z, Yu L, Yuan Z, Zhou Y 2009 Sci. Sin.: Phys. Mech. Astron. 39 746

    [7]

    宁利中, 王娜, 袁喆, 李开继, 王卓运 2014 物理学报 63 104401Google Scholar

    Ning L Z, Wang N, Yuan Z, Li K J, Wang Z Y 2014 Acta Phys. Sin. 63 104401Google Scholar

    [8]

    Barten W, Lucke M, Kamps M, Schmitz R 1995 Phys. Rev. E 51 5662Google Scholar

    [9]

    Jung D, Lucke M 2002 Phys. Rev. Lett. 89 054502Google Scholar

    [10]

    Taraut A V, Smorodin B L, Lucke M 2012 New J. Phys. 14 093055Google Scholar

    [11]

    Mercader I, Batiste O, Alonso A, Knobloch E 2011 J. Fluid Mech. 667 586Google Scholar

    [12]

    宁利中, 王永起, 袁喆, 李开继, 胡彪 2016 科学通报 61 872Google Scholar

    Ning L Z, Wang Y Q, Yuan Z, Li K, Hu B 2016 Chin. Sci. Bull. 61 872Google Scholar

    [13]

    宁利中, 齐昕, 周洋, 余荔 2009 物理学报 58 2528Google Scholar

    Ning L Z, Qi X, Zhou Y, Yu L 2009 Acta Phys. Sin. 58 2528Google Scholar

    [14]

    Zhao B X, Tian Z F 2015 Phys. Fluids 27 074102Google Scholar

    [15]

    宁利中, 刘爽, 宁碧波, 袁喆, 田伟利, 渠亚伟 2017 应用力学学报 34 1086Google Scholar

    Ning L Z, Liu S, Ning B B, Yuan Z, Tian W L, Qu Y W 2017 Chin. J. Appl. Mech. 34 1086Google Scholar

    [16]

    Ning L Z, Qi X, Yuan Z, Shi F 2008 J. Hydrodyn. 20 567Google Scholar

    [17]

    Jung C, Lucke M, Buchel P 1996 Phys. Rev. E 54 1510Google Scholar

    [18]

    胡彪, 宁利中, 宁碧波, 田伟利 2017 科学通报 62 4278Google Scholar

    Hu B, Ning L Z, Ning B B, Tian W L 2017 Chin. Sci. Bull. 62 4278Google Scholar

    [19]

    宁利中, 胡彪, 宁碧波, 田伟利 2016 物理学报 65 214401Google Scholar

    Ning L Z, Hu B, Ning B B, Tian W L 2016 Acta Phys. Sin. 65 214401Google Scholar

    [20]

    Daniels K E, Bodenschatz E 2002 Phys. Rev. Lett. 88 034501Google Scholar

    [21]

    Daniels K E, Plappb B, Bodenschatz E 2000 Phys. Rev. Lett. 84 5320Google Scholar

    [22]

    Daniels K E, Brausch O, Pesch W, Bodenschatz E 2008 J. Fluid Mech. 597 261Google Scholar

    [23]

    Busse F H, Clever R M 1992 J. Eng. Math 26 1Google Scholar

    [24]

    宁利中, 吴昊, 宁碧波, 田伟利, 宁景昊 2019 应用数学和力学 40 398Google Scholar

    Ning L Z, Wu H, Ning B B, Tian W L, Ning J H 2019 Appl. Math. Mech. 40 398Google Scholar

    [25]

    宁利中, 张珂, 宁碧波, 吴昊, 田伟利 2020 应用力学学报 37 737Google Scholar

    Ning L Z, Zhang K, Ning B B, Wu H, Tian W L 2020 Chin. J. Appl. Mech. 37 737Google Scholar

    [26]

    Patankar S V 1982 Computation of Heat Transfer and Fluid Flow (Minneapolis: University of Minnesota) pp20−30

    [27]

    胡彪, 宁利中, 宁碧波, 田伟利, 吴昊, 宁景昊 2017 水动力学研究与进展 32 336Google Scholar

    Hu B, Ning L Z, Ning B B, Tian W L, Wu H, Ning J H 2017 Chin. J. Hydrodyn. 32 336Google Scholar

  • [1] 秦铭宏, 赖强, 吴永红. 具有无穷共存吸引子的简单忆阻混沌系统的分析与实现. 物理学报, 2022, 71(16): 160502. doi: 10.7498/aps.71.20220593
    [2] 扶龙香, 贺少波, 王会海, 孙克辉. 离散忆阻混沌系统的Simulink建模及其动力学特性分析. 物理学报, 2022, 71(3): 030501. doi: 10.7498/aps.71.20211549
    [3] 刘春杰, 赵新军, 高志福, 蒋中英. 高分子混合刷吸附/脱附蛋白质的模型化研究. 物理学报, 2021, 70(22): 224701. doi: 10.7498/aps.70.20211219
    [4] 李春曦, 程冉, 叶学民. 接触角迟滞和气-液界面张力温度敏感性对液滴蒸发动态特性的影响. 物理学报, 2021, 70(20): 204701. doi: 10.7498/aps.70.20210294
    [5] 徐珂, 许龙, 周光平. 考虑水蒸气蒸发和冷凝的球状泡群中泡的动力学特性. 物理学报, 2021, 70(19): 194301. doi: 10.7498/aps.70.20210045
    [6] 扶龙香, 贺少波, 王会海, 孙克辉. 离散忆阻混沌系统的Simulink建模及其动力学特性分析. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211549
    [7] 肖利全, 段书凯, 王丽丹. 基于Julia分形的多涡卷忆阻混沌系统. 物理学报, 2018, 67(9): 090502. doi: 10.7498/aps.67.20172761
    [8] 宁利中, 胡彪, 宁碧波, 田伟利. Poiseuille-Rayleigh-Bnard流动中对流斑图的分区和成长. 物理学报, 2016, 65(21): 214401. doi: 10.7498/aps.65.214401
    [9] 包涵, 包伯成, 林毅, 王将, 武花干. 忆阻自激振荡系统的隐藏吸引子及其动力学特性. 物理学报, 2016, 65(18): 180501. doi: 10.7498/aps.65.180501
    [10] 王小发, 李骏. 短外腔偏振旋转光反馈下1550 nm垂直腔面发射激光器的动力学特性研究. 物理学报, 2014, 63(1): 014203. doi: 10.7498/aps.63.014203
    [11] 李春曦, 裴建军, 叶学民. 倾斜粗糙壁面上含不溶性活性剂溶液的动力学特性. 物理学报, 2013, 62(21): 214704. doi: 10.7498/aps.62.214704
    [12] 蒋树庆, 甯家敏, 陈法新, 叶繁, 薛飞彪, 李林波, 杨建伦, 陈进川, 周林, 秦义, 李正宏, 徐荣昆, 许泽平. Z箍缩动态黑腔动力学及辐射特性初步实验研究. 物理学报, 2013, 62(15): 155203. doi: 10.7498/aps.62.155203
    [13] 李洪伟, 周云龙, 刘旭, 孙斌. 基于随机子空间结合稳定图的气液两相流型分析. 物理学报, 2012, 61(3): 030508. doi: 10.7498/aps.61.030508
    [14] 黄丽莲, 辛方, 王霖郁. 新分数阶超混沌系统的研究与控制及其电路实现. 物理学报, 2011, 60(1): 010505. doi: 10.7498/aps.60.010505
    [15] 吴渝, 杨晶晶, 王利. Swarm模型突现行为的动力学特性分析. 物理学报, 2011, 60(10): 108902. doi: 10.7498/aps.60.108902
    [16] 和兴锁, 李雪华, 邓峰岩. 平面柔性梁的刚-柔耦合动力学特性分析与仿真. 物理学报, 2011, 60(2): 024502. doi: 10.7498/aps.60.024502
    [17] 包伯成, 周国华, 许建平, 刘中. 斜坡补偿电流模式控制开关变换器的动力学建模与分析. 物理学报, 2010, 59(6): 3769-3777. doi: 10.7498/aps.59.3769
    [18] 蒋飞, 刘中, 胡文, 包伯成. 连续混沌调频信号的动力学设计与分析. 物理学报, 2010, 59(1): 116-122. doi: 10.7498/aps.59.116
    [19] 郑桂波, 金宁德. 两相流流型多尺度熵及动力学特性分析. 物理学报, 2009, 58(7): 4485-4492. doi: 10.7498/aps.58.4485
    [20] 蔡国梁, 谭振梅, 周维怀, 涂文桃. 一个新的混沌系统的动力学分析及混沌控制. 物理学报, 2007, 56(11): 6230-6237. doi: 10.7498/aps.56.6230
  • 期刊类型引用(2)

    1. 朱玉枝,李友荣. 矩形流道内密度极值流体Poiseuille-Rayleigh-Bénard对流传热研究. 工程热物理学报. 2024(11): 3480-3486 . 百度学术
    2. 宁利中,宁碧波,郝建武,田伟利. 长矩形腔体中混合流体的对称对传波. 应用力学学报. 2021(06): 2392-2397 . 百度学术

    其他类型引用(3)

计量
  • 文章访问数:  6862
  • PDF下载量:  73
  • 被引次数: 5
出版历程
  • 收稿日期:  2019-12-21
  • 修回日期:  2020-03-25
  • 刊出日期:  2020-06-20

/

返回文章
返回