搜索

x

留言板

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

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

基于GPU的二维梯形空腔流的格子Boltzmann模拟与分析

陈百慧 施保昌 汪垒 柴振华

引用本文:
Citation:

基于GPU的二维梯形空腔流的格子Boltzmann模拟与分析

陈百慧, 施保昌, 汪垒, 柴振华

GPU based lattice Boltzmann simulation and analysis of two-dimensional trapezoidal cavity flow

Chen Bai-Hui, Shi Bao-Chang, Wang Lei, Chai Zhen-Hua
PDF
HTML
导出引用
  • 采用格子Boltzmann方法模拟上下壁面驱动的二维梯形空腔流, 并使用GPU-CUDA程序进行加速计算. 主要采用本征正交分解方法, 分析了4种壁面驱动条件的流场模态, 并探究了雷诺数和驱动速度方向对流场形态的影响. 结果表明: 1)当上壁面单驱动(T1a)时, 若雷诺数为1000—8000, 流场处于稳态流动; 雷诺数为8500时, 流场处于周期性非稳态流动; 雷诺数大于10000时, 流场处于非周期非稳态流动. 2)当下壁面单驱动(T1b)时, 若雷诺数在1000—8000之间, 流场处于稳态流动; 雷诺数增大至11500时, 流场处于周期性非稳态流动; 雷诺数大于12500时, 流场进入非周期非稳态流动. 3)当上下壁面同方向同速度双驱动(T2a)时, 若雷诺数在1000—10000区间, 流场均为稳态流动; 雷诺数为12500—15000时, 流场处于周期性非稳态流动; 当雷诺数大于20000时, 流场为非周期非稳态流动. 4)当上下壁面反方向同速度双驱动(T2b)时, 若雷诺数在1000—5000之间, 流场处于稳态流动; 雷诺数为6000时, 流场处于周期性非稳态流动; 雷诺数大于8000时, 流场为非周期非稳态流动.
    In this study, we utilize the lattice Boltzmann method to investigate the flow behavior in a two-dimensional trapezoidal cavity, which is driven by both sides on the upper wall and lower wall. Our calculations are accelerated through GPU-CUDA software. We conduct an analysis of the flow field mode by using proper orthogonal decomposition. The effects of various parameters, such as Reynolds number (Re) and driving direction, on the flow characteristics are examined through numerical simulations. The results are shown below. 1) For the upper wall drive (T1a), the flow field remains stable, when the Re value varies from 1000 to 8000. However, when Re = 8500, the flow field becomes periodic but unstable. The velocity phase diagram at the monitoring point is a smooth circle, and the energy values of the first two modes dominate the energy of the whole field. Once Re exceeds 10000, the velocity phase diagram turns irregular and the flow field becomes aperiodic and unsteady. 2) For the lower wall drive (T1b), the flow is stable when Re value is in a range of 1000-8000, and it becomes periodic and unsteady when Re = 11500. The energy values of the first three modes appear relatively large. When Re is greater than 12500, the flow field becomes aperiodic and unsteady. At this time, the phase diagram exhibits a smooth circle, with the energy values of the first two modes almost entirely dominating the entire energy. 3) For the case of upper wall and lower wall moving in the same direction at the same speed (T2a), the flow field remains stable when Re changes from 1000 to 10000. When Re varies from 12500 to 15000, the flow becomes periodic and unstable. The velocity phase diagram is still a smooth circle, with the first two modes still occupying a large portion of the energy. Once Re exceeds 20000, the energy proportions of the first three modes significantly decrease, and the flow becomes aperiodic and unsteady. 4) For the case in which the upper wall and lower wall are driven in opposite directions at the same velocity (T2b), the flow field remains stable when Re changes from 1000 to 5000. When Re = 6000, the energy of the first mode accounts for 86%, and the flow field becomes periodic but unstable. When Re exceeds 8000, the energy proportions of the first three modes decrease significantly, and the flow field becomes aperiodic and unsteady.
      通信作者: 施保昌, shibc@hust.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 12072127, 51836003)资助的课题
      Corresponding author: Shi Bao-Chang, shibc@hust.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 12072127, 51836003)
    [1]

    Peng Y F, Shiau Y H, Hwang R R 2003 Comput. Fluids 32 337Google Scholar

    [2]

    Isaev S A, Baranov P A, Sudakov A G, et al. 2008 Thermophys. Aeromech. 15 463Google Scholar

    [3]

    Balachandar S, Eaton J K 2010 Annu. Rev. Fluid Mech. 42 111Google Scholar

    [4]

    Kuhlmann H C, Romanò F 2019 The Lid-Driven Cavity (Austria: Springer International Publishing) p233

    [5]

    Nonino C, Del Giudice S 1988 Int. J. Numer. Methods Eng. 25 313Google Scholar

    [6]

    Alleborn N, Raszillier H, Durst F 1999 Int. J. Heat Mass Tran. 42 833Google Scholar

    [7]

    Talebi F, Mahmoudi A H, Shahi M 2010 Int. Commun. Heat Mass 37 79Google Scholar

    [8]

    Chowdhury M, Kumar B V R 2023 J. Nonnewton. Fluid Mech. 312 104975Google Scholar

    [9]

    Ribbens C J, Watson L T, Wang C Y 1994 J. Comput. Phys. 112 173Google Scholar

    [10]

    McQuain W D, Ribbens C J, Wang C Y, et al. 1994 Comput. Fluids 23 613Google Scholar

    [11]

    Ismael M A, Chamkha A J 2015 Numer. Heat Transfer, Part A 68 312Google Scholar

    [12]

    Mebarek-Oudina F, Laouira H, Hussein A K, et al. 2022 Mathematics 10 929Google Scholar

    [13]

    Mondal P, Mahapatra T R 2021 Int. J. Mech. Sci. 208 106665Google Scholar

    [14]

    Kareem A K, Mohammed H A, Hussein A K, et al. 2016 Int. Commun. Heat Mass 77 195Google Scholar

    [15]

    Rashad A M, Sivasankaran S, Mansour M A, et al. 2017 Numer. Heat Transfer, Part A 71 1223Google Scholar

    [16]

    Darr J H, Vanka S P 1991 Phys. Fluids 3 385Google Scholar

    [17]

    Paramane S B, Sharma A 2008 Numer. Heat Transfer, Part B 54 84Google Scholar

    [18]

    张恒, 任峰, 胡海豹 2021 物理学报 70 184703Google Scholar

    Zhang H, Ren F, Hu H B 2021 Acta Phys. Sin. 70 184703Google Scholar

    [19]

    Yang F, Shi X, Guo X, et al. 2012 Energy Procedia 16 639Google Scholar

    [20]

    Nagendra K, Lakshmisha K N 2009 Int. J. Numer. Methods Heat Fluid Flow 19 790Google Scholar

    [21]

    Perumal D A, Dass A K 2011 Comput. Math. Appl. 61 3711Google Scholar

    [22]

    Nemati H, Farhadi M, Sedighi K, et al. 2010 Int. Commun. Heat Mass 37 1528Google Scholar

    [23]

    Guo X, Zhong C, Zhuo C, et al. 2014 Theor. Comput. Fluid Dyn. 28 215Google Scholar

    [24]

    Patil D V, Lakshmisha K N, Rogg B 2006 Comput. Fluids 35 1116Google Scholar

    [25]

    Hou S, Zou Q, Chen S, et al. 1995 J. Comput. Phys. 118 329Google Scholar

    [26]

    Chai Z H, Shi B C, Zheng L 2006 Chin. Phys. 15 1855Google Scholar

    [27]

    Yuana K A, Budiana E P 2019 IOP Conf. Ser.: Mater. Sci. Eng. 546 052088Google Scholar

    [28]

    Perumal D A, Dass A K 2008 WIT Trans. Eng. Sci. 59 45Google Scholar

    [29]

    Bo A N, Mellibovsky F, Bergada J M, et al. 2020 Appl. Math. Modell. 82 469Google Scholar

    [30]

    Shahid H, Yaqoob I, Khan W A, et al. 2021 Int. Commun. Heat Mass 129 105658Google Scholar

    [31]

    Perumal D A 2018 Therm. Sci. Eng. Prog. 6 48Google Scholar

    [32]

    Guo Y, Bennacer R, Shen S, et al. 2010 Int. J. Numer. Methods Heat Fluid Flow 20 130Google Scholar

    [33]

    Shahid H, Yaqoob I, Khan W A, et al. 2021 Int. Commun. Heat Mass 128 105552Google Scholar

    [34]

    An B, Bergada J M, Mellibovsky F 2019 J. Fluid Mech. 875 476Google Scholar

    [35]

    An B, Guo S, Bergadà J M 2023 Appl. Sci. 13 888Google Scholar

    [36]

    Sidik N A C, Munir F A 2012 Arab. J. Sci. Eng. 37 1723Google Scholar

    [37]

    Zhang T, Shi B, Chai Z 2010 Compu. Fluids 39 1977Google Scholar

    [38]

    Qian Y H, D’Humières D, Lallemand P 1992 EPL 17 479Google Scholar

    [39]

    Ziegler D P 1993 J. Stat. Phys. 71 1171Google Scholar

    [40]

    Guo Z L, Zheng C G, Shi B C 2002 Chin. Phys. 11 366Google Scholar

    [41]

    Guo Z L, Zheng C G, Shi B C 2002 Phys. Fluids 14 2007Google Scholar

  • 图 1  二维等腰梯形空腔示意图

    Fig. 1.  Coordinate of two-dimensional isosceles trapezoidal cavity

    图 2  曲边界上各节点示意图

    Fig. 2.  Geometry condition of curved wall boundary

    图 3  对比Re = 100和500时垂直和水平中线上的速度 (a), (c) 沿垂直中线的速度u和纵坐标y的关系图; (b), (d) 沿水平中线的速度v和横坐标x的关系图

    Fig. 3.  Comparison of the velocity profiles along the vertical and horizontal lines through cavity center at Re = 100, 500: (a) and (c) The y-component of velocity u profiles along the vertical line; (b) and (d) x-component of velocity v profiles along the horizontal line.

    图 4  对比网格数为$ 256\times128 $, $ 512\times256 $, $ 1024\times512 $, Re = 500, 10000时垂直和水平中线上的速度. $ U_1=0.1 $, $ U_2=0 $, Re = 500时(a)沿垂直中线的速度u和纵坐标y的关系图, (b)沿水平中线的速度v和横坐标x的关系图. $ U_1=0.1 $, $ U_2=0.1 $, Re = 10000时(c) 沿垂直中线的速度u和纵坐标y的关系图, (d) 沿水平中线的速度v和横坐标x的关系图

    Fig. 4.  Comparison of the velocity profiles along the vertical and horizontal lines at Re = 500, 10000 for $ 256\times128 $, $ 512\times256 $, $ 1024\times512 $. (a) The y-component of velocity u profiles along the vertical line, (b) x-component of velocity v profiles along the horizontal line at $ U_1=0.1 $, $ U_2=0 $ and Re = 500. (c) The y-component of velocity u profiles along the vertical line, (d) x-component of velocity v profiles along the horizontal line at $ U_1=0.1 $, $ U_2=0.1 $ and Re = 10000.

    图 5  4种不同壁面驱动速度下的梯形空腔示意图 (a) Case T1a, 上壁面向右单驱动; (b) Case T1b, 下壁面向右单驱动; (c) Case T2a, 上壁面和下壁面同时向右双驱动; (d) Case T2b, 上壁面向右驱动, 下壁面向左双驱动

    Fig. 5.  Trapezoidal cavity under four different velocity boundaries: (a) Case T1a, only the upper wall-driven; (b) Case T1b, only the lower wall-driven; (c) Case T2a, the upper and lower wall are driven at the same speed and in the same direction; (d) Case T2b, the upper and lower wall are driven at the same speed and in the opposite direction

    图 6  Case T1a条件下, Re为1000—8000时的流场形态 (a) Re = 1000; (b) Re = 3200; (c) Re = 5000; (d) Re = 6000; (e) Re = 8000

    Fig. 6.  Stream function contours of Reynolds numbers from 1000 to 8000 for Case T1a: (a) Re = 1000; (b) Re = 3200; (c) Re = 5000; (d) Re = 6000; (e) Re = 8000

    图 7  Case T1a条件下, Re为1000—8000时, (a) $ x=1 $的中线上, 纵坐标y和水平速度u的关系, (b) $ y=0.5 $ 的中线上, 横坐标x和竖直速度v的关系

    Fig. 7.  (a) Variation of the velocity u magnitude along $ x=1 $, (b) variation of the velocity v magnitude along $ y=0.5 $ for Re from 1000 to 8000 for Case T1a.

    图 8  Case T1a条件下, Re = 8500时呈现周期性非稳态流动 (a) 时间步$ t=3000000 $时的流场图; (b) 监控点处速度相图; (c) 速度监控点处速度随时间的变化

    Fig. 8.  Stream function contours are periodic and unsteady at Re = 8500 for Case T1a: (a) Stream function contours at time steps $ t=3000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 9  Case T1a条件下, Re = 8500时的各阶模态图 (a)速度平均场模态图; (b) 前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 9.  Each order of mode for Case T1a when Re = 8500: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 10  Case T1a条件下, Re = 10000时呈现非周期非稳态流动 (a)时间步$ t=5000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 10.  Stream function contours are non-periodic and unsteady at Re = 10000 for Case T1a: (a) Stream function contours at time steps $ t=5000000 $; (b) velocity phase diagrams at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 11  Case T1a条件下, Re = 12500时呈现非周期非稳态流动 (a)时间步$ t=5000000 $时的流场图; (b)监控点处速度相图; (c)监控点处速度随时间的变化

    Fig. 11.  Stream function contours are non-periodic and unsteady at Re = 12500 for Case T1a: (a) Stream function contours at time steps $ t=5000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 12  Case T1b条件下, Re为1000—8000时的流场形态 (a) Re = 1000; (b) Re = 3200; (c) Re = 5000; (d) Re = 6000; (e) Re = 8000

    Fig. 12.  Stream function contours of Reynolds numbers from 1000 to 8000 for Case T1b: (a) Re = 1000; (b) Re = 3200; (c) Re = 5000; (d) Re = 6000; (e) Re = 8000

    图 13  Case T1b条件下, Re为1000—8000时, (a) $ x=1 $的中线上, 纵坐标y 和水平速度u的关系; (b) $ y=0.5 $的中线上, 横坐标x和竖直速度v的关系

    Fig. 13.  (a) Variation of the velocity u magnitude along $ x=1 $, (b) variation of the velocity v magnitude along $ y=0.5 $ for Re from 1000 to 8000 for Case T1b.

    图 14  Case T1b条件下, Re = 11500时呈现周期性非稳态流动 (a)时间步$ t=5000000 $时的流场图; (b)监控点处速度相图; (c)监控点处速度随时间变化图

    Fig. 14.  (a) Stream function contours are periodic and unsteady at Re = 11500 for Case T1b: (a) Stream function contours at time steps $ t=5000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 15  Case T1b条件下, Re = 11500时的各阶模态图 (a)速度平均场模态图; (b)前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 15.  Each order of mode for Case T1b when Re = 11500: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 16  Case T1b条件下, Re = 12500时呈非周期非稳态流动 (a)时间步$ t=6000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 16.  Stream function contours are non-periodic and unsteady at Re = 12500 for Case T1b: (a) Stream function contours at time steps $ t=6000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 17  Case T1b条件下, Re = 15000时呈非周期非稳态流动 (a)时间步$ t=6000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 17.  Stream function contours are non-periodic and unsteady at Re = 15000 for Case T1b: (a) Stream function contours at time steps $ t=6000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 18  Case T2a条件下, Re为1000—10000时的流场形态 (a) Re = 1000; (b) Re = 3200; (c) Re = 5000; (d) Re = 6000; (e) Re = 8000; (f) Re = 10000

    Fig. 18.  Stream function contours of Reynolds numbers from 1000 to 10000 for Case T2a: (a) Re = 1000; (b) Re = 3200; (c) Re = 5000; (d) Re = 6000; (e) Re = 8000; (f) Re = 10000

    图 19  Case T2a条件下, Re为1000—10000时, (a) $ x=1 $的中线上, 纵坐标y 和水平速度u的关系; (b) $ y=0.5 $的中线上, 横坐标x和竖直速度v的关系

    Fig. 19.  (a) Variation of the velocity u magnitude along $ x=1 $, (b) variation of the velocity v magnitude along $ y=0.5 $ at Re = 1000–10000 for Case T2a.

    图 20  Case T2a条件下, Re = 12500时呈现周期性非稳态流动 (a)时间步$ t=6000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 20.  Stream function contours are periodic and unsteady at Re = 12500 for Case T2a: (a) Stream function contours at time steps $ t=6000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 21  Case T2a条件下, Re = 12500时的各阶模态图 (a)速度平均场模态图; (b)前n 阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 21.  Each order of mode for Case T2a when Re = 12500: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 22  Case T2a条件下, Re = 15000时呈现周期性非稳态流动 (a)时间步$ t=6000000 $时的流场图; (b) 监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 22.  Stream function contours are periodic and unsteady at Re = 15000 for Case T2a: (a) Stream function contours at time steps $ t=6000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 23  Case T2a条件下, Re = 15000时的各阶模态图 (a)速度平均场模态图; (b)前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 23.  Each order of mode for Case T2a when Re = 15000: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 24  Case T2a条件下, Re = 20000时呈现非周期非稳态流动 (a)时间步$ t=6000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 24.  Stream function contours are non-periodic and unsteady at Re = 20000 for Case T2a: (a) Stream function contours at time steps $ t=6000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 25  Case T2a条件下, Re = 20000时的各阶模态图 (a)速度平均场模态图; (b)前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 25.  Each order of mode for Case T2a when Re = 20000: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity.

    图 26  Case T2a条件下, Re = 25000时呈非周期非稳态流动 (a)时间步$ t=8000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 26.  Stream function contours are non-periodic and unsteady at Re = 25000 for Case T2a: (a) Stream function contours at time steps $ t=8000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 27  Case T2a条件下, Re = 25000时的各阶模态图 (a)速度平均场模态图; (b)前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 27.  Each order of mode for Case T2a when Re = 25000: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 28  Case T2b条件下, Re为1000—5000时的流场形态 (a) Re = 1000; (b) Re = 3200; (c) Re = 5000

    Fig. 28.  Stream function contours of Reynolds numbers from 1000 to 5000 for Case T2b: (a) Re = 1000; (b) Re = 3200; (c) Re = 5000

    图 29  Case T2b条件下, Re为1000—5000时, (a) $ x=1 $的中线上, 纵坐标y 和水平速度u的关系; (b) $ y=0.5 $的中线上, 横坐标x和竖直速度v的关系

    Fig. 29.  (a) Variation of the velocity u magnitude along $ x=1 $, (b) variation of the velocity v magnitude along $ y=0.5 $ at Re = 1000–5000 for Case T2b.

    图 30  Case T2b条件下, Re = 6000时呈现周期性非稳态流动 (a)时间步$ t=3000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 30.  Stream function contours are periodic and unsteady at Re = 6000 for Case T2b: (a) Stream function contours at time steps $ t=3000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 31  Case T2b条件下, Re = 6000时的各阶模态图 (a)速度平均场模态图; (b)前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 31.  Each order of mode for Case T2b when Re = 6000: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 32  Case T2b条件下, Re = 8000时呈现非周期非稳态流动 (a)时间步$ t=5000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间变化图

    Fig. 32.  Stream function contours are non-periodic and unsteady at Re = 8000 for Case T2b: (a) Stream function contours at time steps $ t=5000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 33  Case T2b条件下, Re = 8000时的各阶模态图 (a)速度平均场模态图; (b)前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 33.  Each order of mode for Case T2b when Re = 8000: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 34  Case T2b条件下, Re = 10000时呈现非周期性非稳态流动 (a)时间步$ t=10000000 $时的流场图; (b)监控点处速度相图; (c)速度监控点处速度随时间的变化

    Fig. 34.  Stream function contours are non-periodic and unsteady at Re = 10000 for Case T2b: (a) Stream function contours at time steps $ t=10000000 $; (b) velocity phase diagram at the velocity monitoring point; (c) velocity over time at the velocity monitoring point

    图 35  Case T2b条件下, Re = 10000时的各阶模态图 (a) 速度平均场模态图; (b)前n阶模态的能量占比; (c)各阶模态的能量占比; (d)—(f)一阶、二阶、三阶速度模态图

    Fig. 35.  Each order of mode for Case T2b when Re = 10000: (a) The mean field modal diagrams of velocity; (b) energy share of the first n order of mode for velocity; (c) energy share of each order of mode for velocity; (d)–(f) the first-order, second-order and third-order mode of velocity

    图 36  4种驱动情况下, Re对流场形态的影响图. 黑色圆圈表示流场处于稳态流动状态, 红色三角形表示流场处于周期性非稳态流动状态, 蓝色方框表示流场处于非周期性非稳态流动状态

    Fig. 36.  Effect of Reynolds number on the flow pattern for four different cases. The black circle indicates that the flow field is steady, the red triangle shows that the flow field is periodic but still unsteady, and the blue box indicates that the flow field is non-periodic and unsteady

    表 1  CPU, OpenAcc, CUDA的计算时间和加速比

    Table 1.  Compute time and acceleration ratio of CPU, OpenAcc, CUDA

    计算方式CPUOpenAccCUDA加速比(CPU/OpenAcc)加速比(CPU/CUDA)
    计算时间/s6459.28408.7548.9515.80131.96
    下载: 导出CSV

    表 2  对比Re = 100, 1000, 3200时流场中的涡心位置

    Table 2.  Comparison of the vortex center position in the cavity at Re = 100, 1000, 3200

    Re$x_1$$y_1$$x_2$$y_2$
    $x_{\rm 2 l}$$x_{\rm 2 r}$$y_{\rm 2 l}$$y_{\rm 2 r}$
    100Zhang et al.[37]0.57210.4212
    Present0.57220.4210
    Error/%0.030.05
    1000Zhang et al.[37]0.54730.35580.34230.01800.63510.0450
    Present0.54790.35610.34280.01790.63700.0451
    Error/%0.110.090.150.420.300.12
    3200Zhang et al.[37]0.7140.43920.33780.34910.45040.0788
    Present0.72250.44480.34270.34320.45390.0809
    Error/%1.181.281.461.700.772.66
    下载: 导出CSV

    表 3  Case T1a条件下, Re = 1000, 3200, 5000, 6000, 8000时流场中一级涡、二级涡和三级涡的涡心位置

    Table 3.  The first, second and third primary eddies characteristics at Re = 1000, 3200, 5000, 6000, 8000 for Case T1a.

    Re$x_1$$y_1$$x_2$$y_2$$x_3$$y_3$
    $x_{\rm 2 l}$$x_{\rm 2 r}$$y_{\rm 2 l}$$y_{\rm 2 r}$$x_{\rm 3 l}$$x_{\rm 3 r}$$y_{\rm 3 l}$$y_{\rm 3 r}$
    1000286.89148.44153.65358.669.9311.39
    3200359.03180.97163.84158.09142.9015.17
    5000192.02134.52241.44384.21241.30192.66155.94347.5213.7725.36
    6000200.42132.50261.28393.16242.12196.81155.06353.9812.9335.50
    8000210.51130.69290.17405.75243.39202.60154.62358.0412.4837.51
    下载: 导出CSV

    表 4  Case T1b条件下, Re = 1000, 3200, 5000, 6000, 8000时流场中一级涡, 二级涡和三级涡的涡心位置

    Table 4.  The first, second and third primary eddies characteristics at Re = 1000, 3200, 5000, 6000, 8000 for Case T1b.

    Re$x_1$$y_1$$x_2$$y_2$$x_3$$y_3$
    $x_{\rm 2 l}$$x_{\rm 2 r}$$y_{\rm 2 l}$$y_{\rm 2 r}$$x_{\rm 3 l}$$x_{\rm 3 r}$$y_{\rm 3 l}$$y_{\rm 3 r}$
    1000282.15113.99102.26426.01194.18207.8417.85492.65237.54247.61
    3200280.78118.95110.43412.34158.93218.7225.99487.95229.41245.42
    5000280.79118.95110.43412.34158.95218.7225.97487.81229.78245.57
    6000280.79118.95110.43412.34158.94218.7326.10487.84229.99245.57
    8000283.22121.00121.00408.63117.03222.4661.80469.33210.11233.00
    下载: 导出CSV

    表 5  Case T2a条件下, Re = 1000, 3200, 5000, 6000, 8000, 10000时流场中的涡心位置

    Table 5.  The first, second and third primary eddies characteristics at Re = 1000, 3200, 5000, 6000, 8000, 10000 for Case T2a.

    Re$x_1$$y_1$$x_2$$y_2$$x_3$$y_3$
    $x_{\rm 2 l}$$x_{\rm 2 r}$$y_{\rm 2 l}$$y_{\rm 2 r}$$x_{\rm 3 l}$$x_{\rm 3 r}$$y_{\rm 3 l}$$y_{\rm 3 r}$
    $x_{\rm 3 r_1}$$x_{\rm 3 r_2}$$y_{\rm 3 r_1}$$y_{\rm 3 r_2}$
    1000265.6583.21166.69390.56209.14202.40412.69422.93101.61117.56
    3200276.6490.43139.51404.39190.46205.4934.89410.36421.46216.82115.62134.97
    5000282.1092.07141.76408.64186.10206.4440.16410.52421.23222.49120.92138.91
    6000284.2192.69142.88410.13184.42206.7641.25410.77421.29223.56122.65140.12
    8000287.0793.68144.54412.15182.08207.1442.18411.69422.02224.22124.32141.38
    10000288.9494.48145.67413.49180.45207.3343.32412.89423.35223.99124.48141.80
    下载: 导出CSV

    表 6  Case T2b条件下, Re = 1000, 3200, 5000时流场中的涡心位置

    Table 6.  The first, second and third primary eddies characteristics at Re = 1000, 3200, 5000 for Case T2b.

    Re$x_1$$y_1$$x_2$$y_2$$x_3$$y_3$
    $x_{3 {\rm l}}$$x_{3 {\rm r} }$$y_{3 {\rm l} }$$y_{3 {\rm r}}$
    1000259.99132.4459.83183.94
    3200260.97131.1266.16206.2684.50369.55245.4214.67
    5000261.42130.7465.46208.2288.80365.54246.8614.39
    下载: 导出CSV
  • [1]

    Peng Y F, Shiau Y H, Hwang R R 2003 Comput. Fluids 32 337Google Scholar

    [2]

    Isaev S A, Baranov P A, Sudakov A G, et al. 2008 Thermophys. Aeromech. 15 463Google Scholar

    [3]

    Balachandar S, Eaton J K 2010 Annu. Rev. Fluid Mech. 42 111Google Scholar

    [4]

    Kuhlmann H C, Romanò F 2019 The Lid-Driven Cavity (Austria: Springer International Publishing) p233

    [5]

    Nonino C, Del Giudice S 1988 Int. J. Numer. Methods Eng. 25 313Google Scholar

    [6]

    Alleborn N, Raszillier H, Durst F 1999 Int. J. Heat Mass Tran. 42 833Google Scholar

    [7]

    Talebi F, Mahmoudi A H, Shahi M 2010 Int. Commun. Heat Mass 37 79Google Scholar

    [8]

    Chowdhury M, Kumar B V R 2023 J. Nonnewton. Fluid Mech. 312 104975Google Scholar

    [9]

    Ribbens C J, Watson L T, Wang C Y 1994 J. Comput. Phys. 112 173Google Scholar

    [10]

    McQuain W D, Ribbens C J, Wang C Y, et al. 1994 Comput. Fluids 23 613Google Scholar

    [11]

    Ismael M A, Chamkha A J 2015 Numer. Heat Transfer, Part A 68 312Google Scholar

    [12]

    Mebarek-Oudina F, Laouira H, Hussein A K, et al. 2022 Mathematics 10 929Google Scholar

    [13]

    Mondal P, Mahapatra T R 2021 Int. J. Mech. Sci. 208 106665Google Scholar

    [14]

    Kareem A K, Mohammed H A, Hussein A K, et al. 2016 Int. Commun. Heat Mass 77 195Google Scholar

    [15]

    Rashad A M, Sivasankaran S, Mansour M A, et al. 2017 Numer. Heat Transfer, Part A 71 1223Google Scholar

    [16]

    Darr J H, Vanka S P 1991 Phys. Fluids 3 385Google Scholar

    [17]

    Paramane S B, Sharma A 2008 Numer. Heat Transfer, Part B 54 84Google Scholar

    [18]

    张恒, 任峰, 胡海豹 2021 物理学报 70 184703Google Scholar

    Zhang H, Ren F, Hu H B 2021 Acta Phys. Sin. 70 184703Google Scholar

    [19]

    Yang F, Shi X, Guo X, et al. 2012 Energy Procedia 16 639Google Scholar

    [20]

    Nagendra K, Lakshmisha K N 2009 Int. J. Numer. Methods Heat Fluid Flow 19 790Google Scholar

    [21]

    Perumal D A, Dass A K 2011 Comput. Math. Appl. 61 3711Google Scholar

    [22]

    Nemati H, Farhadi M, Sedighi K, et al. 2010 Int. Commun. Heat Mass 37 1528Google Scholar

    [23]

    Guo X, Zhong C, Zhuo C, et al. 2014 Theor. Comput. Fluid Dyn. 28 215Google Scholar

    [24]

    Patil D V, Lakshmisha K N, Rogg B 2006 Comput. Fluids 35 1116Google Scholar

    [25]

    Hou S, Zou Q, Chen S, et al. 1995 J. Comput. Phys. 118 329Google Scholar

    [26]

    Chai Z H, Shi B C, Zheng L 2006 Chin. Phys. 15 1855Google Scholar

    [27]

    Yuana K A, Budiana E P 2019 IOP Conf. Ser.: Mater. Sci. Eng. 546 052088Google Scholar

    [28]

    Perumal D A, Dass A K 2008 WIT Trans. Eng. Sci. 59 45Google Scholar

    [29]

    Bo A N, Mellibovsky F, Bergada J M, et al. 2020 Appl. Math. Modell. 82 469Google Scholar

    [30]

    Shahid H, Yaqoob I, Khan W A, et al. 2021 Int. Commun. Heat Mass 129 105658Google Scholar

    [31]

    Perumal D A 2018 Therm. Sci. Eng. Prog. 6 48Google Scholar

    [32]

    Guo Y, Bennacer R, Shen S, et al. 2010 Int. J. Numer. Methods Heat Fluid Flow 20 130Google Scholar

    [33]

    Shahid H, Yaqoob I, Khan W A, et al. 2021 Int. Commun. Heat Mass 128 105552Google Scholar

    [34]

    An B, Bergada J M, Mellibovsky F 2019 J. Fluid Mech. 875 476Google Scholar

    [35]

    An B, Guo S, Bergadà J M 2023 Appl. Sci. 13 888Google Scholar

    [36]

    Sidik N A C, Munir F A 2012 Arab. J. Sci. Eng. 37 1723Google Scholar

    [37]

    Zhang T, Shi B, Chai Z 2010 Compu. Fluids 39 1977Google Scholar

    [38]

    Qian Y H, D’Humières D, Lallemand P 1992 EPL 17 479Google Scholar

    [39]

    Ziegler D P 1993 J. Stat. Phys. 71 1171Google Scholar

    [40]

    Guo Z L, Zheng C G, Shi B C 2002 Chin. Phys. 11 366Google Scholar

    [41]

    Guo Z L, Zheng C G, Shi B C 2002 Phys. Fluids 14 2007Google Scholar

  • [1] 赖瑶瑶, 陈鑫梦, 柴振华, 施保昌. 基于格子Boltzmann方法的钉扎螺旋波反馈控制. 物理学报, 2024, 73(4): 040502. doi: 10.7498/aps.73.20231549
    [2] 张恒, 任峰, 胡海豹. 基于格子Boltzmann方法的幂律流体二维顶盖驱动流转捩研究. 物理学报, 2021, 70(18): 184703. doi: 10.7498/aps.70.20210451
    [3] 臧晨强, 娄钦. 复杂微通道内非混相驱替过程的格子Boltzmann方法. 物理学报, 2017, 66(13): 134701. doi: 10.7498/aps.66.134701
    [4] 王佐, 张家忠, 王恒. 非正交多松弛系数轴对称热格子Boltzmann方法. 物理学报, 2017, 66(4): 044701. doi: 10.7498/aps.66.044701
    [5] 周光雨, 陈力, 张鸿雁, 崔海航. 基于格子Boltzmann方法的自驱动Janus颗粒扩散泳力. 物理学报, 2017, 66(8): 084703. doi: 10.7498/aps.66.084703
    [6] 张娅, 潘光, 黄桥高. 疏水表面减阻的格子Boltzmann方法数值模拟. 物理学报, 2015, 64(18): 184702. doi: 10.7498/aps.64.184702
    [7] 黄桥高, 潘光, 宋保维. 疏水表面滑移流动及减阻特性的格子Boltzmann方法模拟. 物理学报, 2014, 63(5): 054701. doi: 10.7498/aps.63.054701
    [8] 任晟, 张家忠, 张亚苗, 卫丁. 零质量射流激励下诱发液体相变及其格子Boltzmann方法模拟. 物理学报, 2014, 63(2): 024702. doi: 10.7498/aps.63.024702
    [9] 解文军, 滕鹏飞. 声悬浮过程的格子Boltzmann方法研究. 物理学报, 2014, 63(16): 164301. doi: 10.7498/aps.63.164301
    [10] 史冬岩, 王志凯, 张阿漫. 任意复杂流-固边界的格子Boltzmann处理方法. 物理学报, 2014, 63(7): 074703. doi: 10.7498/aps.63.074703
    [11] 刘邱祖, 寇子明, 贾月梅, 吴娟, 韩振南, 张倩倩. 改性疏水固壁润湿性反转现象的格子Boltzmann方法模拟. 物理学报, 2014, 63(10): 104701. doi: 10.7498/aps.63.104701
    [12] 郭亚丽, 徐鹤函, 沈胜强, 魏兰. 利用格子Boltzmann方法模拟矩形腔内纳米流体Raleigh-Benard对流. 物理学报, 2013, 62(14): 144704. doi: 10.7498/aps.62.144704
    [13] 曾建邦, 李隆键, 蒋方明. 气泡成核过程的格子Boltzmann方法模拟. 物理学报, 2013, 62(17): 176401. doi: 10.7498/aps.62.176401
    [14] 刘邱祖, 寇子明, 韩振南, 高贵军. 基于格子Boltzmann方法的液滴沿固壁铺展动态过程模拟. 物理学报, 2013, 62(23): 234701. doi: 10.7498/aps.62.234701
    [15] 曾建邦, 李隆键, 廖全, 蒋方明. 池沸腾中气泡生长过程的格子Boltzmann方法模拟. 物理学报, 2011, 60(6): 066401. doi: 10.7498/aps.60.066401
    [16] 曾建邦, 李隆键, 廖全, 陈清华, 崔文智, 潘良明. 格子Boltzmann方法在相变过程中的应用. 物理学报, 2010, 59(1): 178-185. doi: 10.7498/aps.59.178
    [17] 石自媛, 胡国辉, 周哲玮. 润湿性梯度驱动液滴运动的格子Boltzmann模拟. 物理学报, 2010, 59(4): 2595-2600. doi: 10.7498/aps.59.2595
    [18] 卢玉华, 詹杰民. 三维方腔温盐双扩散的格子Boltzmann方法数值模拟. 物理学报, 2006, 55(9): 4774-4782. doi: 10.7498/aps.55.4774
    [19] 李华兵, 黄乒花, 刘慕仁, 孔令江. 用格子Boltzmann方法模拟MKDV方程. 物理学报, 2001, 50(5): 837-840. doi: 10.7498/aps.50.837
    [20] 吕晓阳, 李华兵. 用格子Boltzmann方法模拟高雷诺数下的热空腔黏性流. 物理学报, 2001, 50(3): 422-427. doi: 10.7498/aps.50.422
计量
  • 文章访问数:  2727
  • PDF下载量:  59
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-03-21
  • 修回日期:  2023-04-16
  • 上网日期:  2023-05-25
  • 刊出日期:  2023-08-05

/

返回文章
返回