期刊文献+
共找到15篇文章
< 1 >
每页显示 20 50 100
Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method 被引量:4
1
作者 柴振华 施保昌 郑林 《Chinese Physics B》 SCIE EI CAS CSCD 2006年第8期1855-1863,共9页
By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improv... By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination. 展开更多
关键词 multi-relaxation-time lattice Boltzmann method non-equilibrium extrapolation scheme high Reynolds number lid-driven cavity flow
在线阅读 下载PDF
The second Hopf bifurcation in lid-driven square cavity 被引量:1
2
作者 Tao Wang Tiegang Liu Zheng Wang 《Chinese Physics B》 SCIE EI CAS CSCD 2020年第3期119-125,共7页
To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of... To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of theoretical analyses and numerical simulations. In this paper, we study the characteristics of the second Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme recently developed by us. We numerically identify the critical Reynolds number of the second Hopf bifurcation located in the interval of(11093.75, 11094.3604) by bisection. In addition, we find that there are two dominant frequencies in its spectral diagram when the flow is in the status of the second Hopf bifurcation, while only one dominant frequency is identified if the flow is in the first Hopf bifurcation via the Fourier analysis. More interestingly, the flow phase portrait of velocity components is found to make transition from a regular elliptical closed form for the first Hopf bifurcation to a non-elliptical closed form with self-intersection for the second Hopf bifurcation. Such characteristics disclose flow in a quasi-periodic state when the second Hopf bifurcation occurs. 展开更多
关键词 unsteady lid-driven cavity flows SECOND HOPF bifurcation critical REYNOLDS number numerical simulation
在线阅读 下载PDF
类蜂窝结构对腔体流动分岔的影响研究
3
作者 管永科 郭世鹏 +1 位作者 桑为民 安博 《力学学报》 北大核心 2025年第3期569-577,共9页
流动分岔对于厘清流场物理特性意义非凡,对于工程实际中遇到的流动问题分析和相关应用有着非常重要的指导意义.作为流体力学中的经典问题,腔体流动在很多重要领域都被广泛应用,具备极高的研究价值.对于腔体内流流动而言,腔体几何构型对... 流动分岔对于厘清流场物理特性意义非凡,对于工程实际中遇到的流动问题分析和相关应用有着非常重要的指导意义.作为流体力学中的经典问题,腔体流动在很多重要领域都被广泛应用,具备极高的研究价值.对于腔体内流流动而言,腔体几何构型对于流动分岔的类型和临界雷诺数影响较大.针对目前国内外研究现状,部分几何构型对于流动分岔的影响规律尚未明晰.文章通过基于直角网格的LBM算法和流场稳定性分析理论,针对类蜂窝构型的腔体内流开展数值模拟和流动现象机理分析研究,揭示了类蜂窝结构对于内流流场流动分岔的影响机制.对比之前的研究,着重探究类蜂窝构型对于Hopf、Neimark-Sacker以及Turbulence-triggering流动分岔的影响以及该构型内流动演化模式和涡系发展特征,从物理层面分析流动现象的形成机理.研究结果表明,虽然类蜂窝构型不会影响流动分岔的类型以及流场演化模式,但是对于流动分岔的临界值影响较大.可以有效推迟流动失稳和湍流出现,极大地提升了内流流场的稳定性.此外,该流场内的涡系结构及演化规律异常稳定. 展开更多
关键词 顶盖驱动腔体内流 类蜂窝结构 流场稳定性 格子玻尔兹曼方法 流动分岔
在线阅读 下载PDF
任意精度的三点紧致显格式及其在CFD中的应用 被引量:6
4
作者 林建国 谢志华 +1 位作者 周俊陶 林建忠(推荐) 《应用数学和力学》 CSCD 北大核心 2007年第7期843-852,共10页
通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间任意精度的三点紧致显格式的表达式,又由Fourier分析法得到了格式的数值弥散和耗散特性.与以往的高精度紧致差分格式不同,提出的格式不用隐式求解代数方程组并且可以达到任意精度.... 通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间任意精度的三点紧致显格式的表达式,又由Fourier分析法得到了格式的数值弥散和耗散特性.与以往的高精度紧致差分格式不同,提出的格式不用隐式求解代数方程组并且可以达到任意精度.通过方波问题和顶盖方腔流的算例表明,格式在稀疏网格下可以得到很高的精度,不仅能节省计算量,而且易于编程,有很高的计算效率. 展开更多
关键词 任意精度 紧致 三节点 显格式 方腔流
在线阅读 下载PDF
不可压缩粘性流动的CBS有限元解法 被引量:3
5
作者 孙旭 张家忠 +1 位作者 周志宏 徐忠 《计算力学学报》 EI CAS CSCD 北大核心 2010年第5期862-867,共6页
对于二维不可压缩粘性流动,首先通过坐标变换的方式得到了的不含对流项的NS方程,并给出了CBS有限元方法求解的一般过程。结合一类同时含有压力和速度的出口边界条件,对方腔顶盖驱动流、后向台阶绕流和圆柱绕流进行了计算。所得结果与基... 对于二维不可压缩粘性流动,首先通过坐标变换的方式得到了的不含对流项的NS方程,并给出了CBS有限元方法求解的一般过程。结合一类同时含有压力和速度的出口边界条件,对方腔顶盖驱动流、后向台阶绕流和圆柱绕流进行了计算。所得结果与基准解符合良好,验证了CBS算法对于定常、非定常粘性不可压缩流动问题的可行性和所用出口边界条件的无反射特性。特别的,对于圆柱绕流,Re=100时非定常升、阻力系数及漩涡脱落等非定常都得到了较好地模拟,为一进步研究自激振动等更加复杂的非定常流动问题奠定了基础。 展开更多
关键词 CBS 有限元 方腔顶盖驱动流 后向台阶绕流 圆柱绕流
在线阅读 下载PDF
非稳态流动的隐式最小二乘等几何方法 被引量:2
6
作者 陈德祥 刘帅 +1 位作者 徐自力 胡哺松 《计算力学学报》 CAS CSCD 北大核心 2015年第5期639-643,共5页
对非稳态粘性不可压流动问题提出了一种隐式最小二乘等几何计算方法。该方法先用隐式的向后多步差分格式对Navier-Stokes方程进行时间离散,再用Newton法线性化对流项,最后在每个时间步上用最小二乘等几何方法进行求解。根据该算法编制... 对非稳态粘性不可压流动问题提出了一种隐式最小二乘等几何计算方法。该方法先用隐式的向后多步差分格式对Navier-Stokes方程进行时间离散,再用Newton法线性化对流项,最后在每个时间步上用最小二乘等几何方法进行求解。根据该算法编制了计算程序,通过构造解析解的方法验证了程序的正确性,用该程序求解了雷诺数为5000时的非稳态二维顶盖驱动流问题,计算结果捕捉到了流动过程中涡的演化过程,表明本文方法可用于非稳态流动的求解。 展开更多
关键词 非稳态流动 NAVIER-STOKES 方程 最小二乘等几何方法 隐式差分格式 顶盖驱动流
在线阅读 下载PDF
镜像对称顶盖驱动方腔内流过渡流临界特性研究 被引量:2
7
作者 安博 孟欣雨 桑为民 《力学学报》 EI CAS CSCD 北大核心 2022年第9期2409-2418,共10页
流场过渡流临界特性是指流场因流动状态改变而引起的流场物理特性变化.如流动从定常演化为非定常周期性时,流动处于过渡状态的物理性质.它从根本上决定了流动演化模式和流场特性等物理规律,对认清流动现象的形成机理有重要意义.本文在... 流场过渡流临界特性是指流场因流动状态改变而引起的流场物理特性变化.如流动从定常演化为非定常周期性时,流动处于过渡状态的物理性质.它从根本上决定了流动演化模式和流场特性等物理规律,对认清流动现象的形成机理有重要意义.本文在之前腔体内流流场过渡流临界特性研究的基础上,针对镜像对称顶盖驱动方腔内流开展数值模拟和流场稳定性分析研究,捕捉各流动分岔点,如Hopf流动分岔点和Neimark-Sacker流动分岔点等,并揭示其对流场特性的影响;分析流场演化模式,随着雷诺数增大从定常状态依次演化为非定常周期性流动、准周期性流动和湍流;揭示各种流动现象的形成机理,如流动滞后、对称性破坏、能量级串等;分析流场拓扑结构,阐明流场镜像对称性和流场稳定性的关系.本文研究成果有助于揭示该流场的物理特性,进一步完善了内流流场特性的研究.研究发现,针对本文镜像对称方腔顶盖驱动内流,流场稳定性的破坏总是以Hopf流动分岔点的出现而发生并且伴随着流场对称性的破坏;流场演化模式符合经典的Ruelle-Takens模式;流动从定常状态演化至非定常周期性流动时存在流动滞后现象. 展开更多
关键词 镜像对称顶盖驱动方腔 过渡流临界特性 流场稳定性 流场对称性 流动分岔点
在线阅读 下载PDF
改进的基于特征线的N-S方程算子分裂有限元法
8
作者 水庆象 王大国 沈连山 《船舶力学》 EI CSCD 北大核心 2016年第4期381-392,共12页
基于特征线的Navier-Stokes(N-S)方程算子分裂有限元法(CBOS有限元法)的核心是在每一个时间层上,采用算子分裂法将N-S方程的对流项和扩散项分开求解;对流项的求解过程借鉴了简单显式特征线时间离散,显式求解。该文在原CBOS有限元法基础... 基于特征线的Navier-Stokes(N-S)方程算子分裂有限元法(CBOS有限元法)的核心是在每一个时间层上,采用算子分裂法将N-S方程的对流项和扩散项分开求解;对流项的求解过程借鉴了简单显式特征线时间离散,显式求解。该文在原CBOS有限元法基础上推导了一种更加精确的对流项显式离散方法。通过自编程序对方腔流进行数值模拟,表明该算法具有更高的计算精度。低雷诺数下圆柱绕流计算所得的阻力系数、升力系数、斯特劳哈数等与已有数据较为接近,表明该文中算法能较准确地模拟圆柱绕流的流场特性;最后,文中分析了Re=200时圆柱绕流在一个周期内所受升力变化与对应流场中压力和流线演化的关系。 展开更多
关键词 N-S方程 CBOS有限元法 方腔流 圆柱绕流
在线阅读 下载PDF
深纵比对方腔过渡流临界特性的影响研究
9
作者 安博 孟欣雨 +1 位作者 郭世鹏 桑为民 《力学学报》 EI CAS CSCD 北大核心 2023年第6期1247-1256,共10页
流场过渡流临界特性是指流场因流动分岔而引起的流动状态和流场物理特性变化.它从根本上决定了流动演化模式和流场特性等物理规律,对解释复杂流动现象意义重大.文章针对不同深纵比(R∈[0.1,2.0])的顶盖驱动方腔内流开展数值模拟和流场... 流场过渡流临界特性是指流场因流动分岔而引起的流动状态和流场物理特性变化.它从根本上决定了流动演化模式和流场特性等物理规律,对解释复杂流动现象意义重大.文章针对不同深纵比(R∈[0.1,2.0])的顶盖驱动方腔内流开展数值模拟和流场稳定性分析研究.预测Hopf,Neimark-Sacker和period-doubling流动分岔及湍流始现的临界雷诺数;分析流场演化模式,发现对应不同的深纵比,有些流动遵循经典的Ruelle-Takens模式,有些流动则会由周期性流动跃变至湍流;捕捉和分析各种流动现象,如流场稳定性丧失、能量级串、流场拓扑结构变化规律等.研究成果对于揭示深纵比这一几何参数对腔体内流过渡流临界特性的影响规律意义重大,进一步完善了内流流场特性的研究.研究发现,Moffatt效应不仅存在于拥有尖锐夹角的内部流动中,也出现于挤压拉伸的狭长空间;无论是深腔还是浅腔,流场稳定性最初的破坏总是以Hopf流动分岔的出现而开始;就浅腔(R<1.0)而言,随着深纵比逐渐增加,Hopf流动分岔的临界雷诺数越来越小,流动更容易变为非定常状态,说明流场稳定性变得越来越容易被破坏;就深腔(R>1.0)而言,相较于经典方腔驱动内流(R=1.0),流场稳定性更容易丧失;沿纵向的几何外形拉伸并不是提升流场稳定性的强制约束. 展开更多
关键词 顶盖驱动方腔流动 过渡流临界特性 流场稳定性 深纵比 流动分岔
在线阅读 下载PDF
SIMPLE算法的一种改进格式 被引量:3
10
作者 秦剑 潘华辰 +1 位作者 田晓庆 RAHMAN M M 《应用力学学报》 CAS CSCD 北大核心 2021年第6期2398-2404,共7页
在连续性方程中人为加入压力的时间导数项,将人工可压缩方法与压力泊松方程相结合起来,有效抑制了由于守恒变量和非守恒变量转换引起的密度波动;运用改进的非线性动量插值格式对网格界面处的速度进行离散,消除了交错网格下棋盘式的压力... 在连续性方程中人为加入压力的时间导数项,将人工可压缩方法与压力泊松方程相结合起来,有效抑制了由于守恒变量和非守恒变量转换引起的密度波动;运用改进的非线性动量插值格式对网格界面处的速度进行离散,消除了交错网格下棋盘式的压力波动,减小了网格界面处的数值耗散限。为验证该方法的通用性,依次对方腔自然对流、顶盖驱动流和后向阶梯流进行了数值分析,证实该算法在收敛速度、稳定性方面均有较大提升。 展开更多
关键词 SIMPLE算法 人工压缩方法 方腔自然对流 顶盖驱动流 后向阶梯流
在线阅读 下载PDF
惯量松弛因子在格子Boltzmann 方法中的应用 被引量:2
11
作者 王颖娟 龚光彩 +2 位作者 石星 龚子彻 刘永超 《工程数学学报》 CSCD 北大核心 2021年第2期195-206,共12页
格子Boltzmann方法作为一种较成熟的数值模拟方法被广泛应用到了各个领域,尤其在解决多孔介质问题时有其独特的优越性,但当流动问题过于复杂时计算效率较低.因此本文将惯量松弛因子引入到格子Boltzmann方法中,对二维、三维顶盖驱动方腔... 格子Boltzmann方法作为一种较成熟的数值模拟方法被广泛应用到了各个领域,尤其在解决多孔介质问题时有其独特的优越性,但当流动问题过于复杂时计算效率较低.因此本文将惯量松弛因子引入到格子Boltzmann方法中,对二维、三维顶盖驱动方腔流动进行了数值模拟.模拟分别从计算效率、计算精度、以及计算稳定性等方面将使用不同惯量松弛因子所得的结果与基准解进行比较,并进行讨论和分析.结果显示当惯量松弛因子取0.03到0.05之间时能使模拟结果在保持较高精度的同时提高计算效率,而且随着惯量松弛因子的增大计算效率提高得也越快,在工程材料与能源环境领域将有着重要应用. 展开更多
关键词 惯量松弛因子 格子BOLTZMANN方法 顶盖驱动流 计算效率 收敛速度
在线阅读 下载PDF
配置点谱方法-人工压缩法(SCM-ACM)求解不可压缩流体流动
12
作者 牟永强 郝建超 +2 位作者 张敬奎 董华 林欢 《计算力学学报》 EI CAS CSCD 北大核心 2018年第3期304-309,共6页
开发了配置点谱方法SCM(spectral collocation method)与人工压缩法ACM(artificial compressibility method)相结合的方法 SCM-ACM,用于求解不可压缩粘性流动问题。选取典型的方腔顶盖驱动流为研究测试对象,首先建立人工压缩格式的控制... 开发了配置点谱方法SCM(spectral collocation method)与人工压缩法ACM(artificial compressibility method)相结合的方法 SCM-ACM,用于求解不可压缩粘性流动问题。选取典型的方腔顶盖驱动流为研究测试对象,首先建立人工压缩格式的控制方程,其次采用SCM离散控制方程的空间偏微分项,推导出矩阵形式的代数方程,最后测试了SCM-ACM代码的有效性。结果显示,SCM-ACM能够有效求解不可压缩流动问题,并继承了谱方法的指数收敛特性,且具有ACM求解过程简单及易于实施的特点。 展开更多
关键词 SCM-ACM 配置点谱方法 人工压缩法 顶盖驱动流 不可压缩流动
在线阅读 下载PDF
相邻双侧边盖驱动方腔流动的三维线性整体稳定性 被引量:1
13
作者 胡军 李杰权 《气体物理》 2018年第3期1-8,共8页
文章考察了相邻双侧边盖驱动方腔流动(即上壁面向右运动和左侧壁面向下运动)的三维线性整体稳定性.首先,采用Taylor-Hood有限元方法并经由Newton迭代过程计算得到双侧边盖驱动方腔流动的二维稳态基本流.其次, Taylor-Hood有限元在Chebys... 文章考察了相邻双侧边盖驱动方腔流动(即上壁面向右运动和左侧壁面向下运动)的三维线性整体稳定性.首先,采用Taylor-Hood有限元方法并经由Newton迭代过程计算得到双侧边盖驱动方腔流动的二维稳态基本流.其次, Taylor-Hood有限元在Chebyshev Gauss配置点上进行离散,同时Gauss配置点也可以用于线性稳定性方程的高阶有限差分格式离散.然后,离散得到的矩阵形式的广义特征值问题可以结合shift-and-invert算法采用隐式重启Arnoldi方法计算.最后,通过对线性稳定性方程特征值的计算,发现了一个最不稳定的驻定模态和两对对称行波模态.最不稳定的三维驻定模态的临界Reynolds数为Re_c=261.5,远远小于二维不稳定的临界Reynolds数Re_c^(2d)=1 061.7.通过画出这3类三维不稳定模态的流向扰动速度和扰动涡量的空间等值面图像,可以发现不稳定扰动位于稳态基本流的两个主涡区域,因此可以认为主涡区域是三维扰动失稳的主要能量来源地. 展开更多
关键词 边盖驱动方腔流动 整体稳定性 临界Reynolds数 行波模态 Taylor—Hood有限元
在线阅读 下载PDF
一种基于非结构化网格的无压力型投影法
14
作者 李国杰 李杨 +6 位作者 王勇 李建峰 张晓 王憨鹰 黄萌 周丹丹 胡广涛 《榆林学院学报》 2022年第4期8-11,共4页
考虑到流动区域的复杂性及投影法求解非定常问题的高效性,本文发展了一种基于非结构化网格的无压力型投影法,使其能够完成复杂流动区域内非定常流动问题的高效求解。算法实施中,首先基于不包含压力项的预测方程完成中间速度求解,然后根... 考虑到流动区域的复杂性及投影法求解非定常问题的高效性,本文发展了一种基于非结构化网格的无压力型投影法,使其能够完成复杂流动区域内非定常流动问题的高效求解。算法实施中,首先基于不包含压力项的预测方程完成中间速度求解,然后根据包含压力项的修正方程完成压力场求解,最后基于压力场对速度场进行修正。基于本文算法,对方腔顶盖驱动流典型算例进行了模拟分析,通过与前人结果对比,验证了本文算法的准确性。 展开更多
关键词 非结构化网格 投影法 无压力型 顶盖驱动
在线阅读 下载PDF
单个中性球形颗粒在三维方腔中的运动 被引量:1
15
作者 李洋 夏振华 《空气动力学学报》 CSCD 北大核心 2022年第3期109-119,I0002,共12页
采用多松弛格子Boltzmann方法模拟了单个中性球形颗粒在三维顶盖驱动方腔流中的运动。考虑了展向弱受限的对称边界和强受限的固壁边界两种情况下初始位置、颗粒大小以及雷诺数的影响。对于展向弱受限的情形,发现颗粒的初始位置显著影响... 采用多松弛格子Boltzmann方法模拟了单个中性球形颗粒在三维顶盖驱动方腔流中的运动。考虑了展向弱受限的对称边界和强受限的固壁边界两种情况下初始位置、颗粒大小以及雷诺数的影响。对于展向弱受限的情形,发现颗粒的初始位置显著影响着最终的运动轨迹。根据相图被划分为三个区域:外层稳定区、内层稳定区以及涡中心区域。通过对颗粒受力的分解解释了其在极限环上运动的机理。此外,还详细介绍了球形颗粒在极限环上的顺时针旋转运动。随着雷诺数的增加,颗粒逐渐向外围靠近,不断旋转达到相应的极限环轨迹。对于选定的初始位置,观察到在高雷诺数时大颗粒向外迁移,而在低雷诺数时大颗粒的极限环靠近涡中心。对于展向强受限的情形,极限环与颗粒的初始位置无关。随着雷诺数的增加,除方腔左上角外,极限环轨迹有向外迁移的趋势。最后,随着颗粒尺寸的增大,极限环向方腔内部收缩。 展开更多
关键词 格子BOLTZMANN方法 顶盖驱动流 方腔 中性悬浮 球形颗粒
在线阅读 下载PDF
上一页 1 下一页 到第
使用帮助 返回顶部