搜索

x
中国物理学会期刊

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

CSTR: 32037.14.aps.72.20230430

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

CSTR: 32037.14.aps.72.20230430
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.

     

    目录

    /

    返回文章
    返回