共查询到20条相似文献,搜索用时 531 毫秒
1.
2.
3.
4.
挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题.本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象,采用颜色模型模拟出口处粘弹流体和空气的两相流动,通过重新标色获得两种流体的界面,并最终获得胀大的形状.Navier-Stokes方程和本构方程的求解采用双分布函数模型.将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较,发现格子Boltzmann两相模型结果与解析解和实验结果相吻合,相比于单相模型,收敛速度更快,解的稳定性更高.研究了流道尺寸对胀大率的影响,并对挤出胀大的内在机理进行了分析. 相似文献
5.
多组份流体质量扩散的格子Boltzmann方法 总被引:1,自引:0,他引:1
建立了多组份流体、多速度格子Boltzmann模型.利用Chapman-Enskog渐近展开法推导了流体质量扩散方程和运动方程,并且给出了相应的扩散系数和粘性系数表达式.通过正弦波的衰减过程测量了两流体间的扩散系数.测量值与理论预测值相吻合.作为本模型的一个应用实例,对圆形区域内的静止流体在具有均匀来流速度的另外一种流体中的对流扩散问题进行了数值模拟. 相似文献
6.
基于格子 Bhatnagar-Gross-Krook模型的地震压力波模拟 总被引:1,自引:0,他引:1
给出一种新的用于模拟地震压力波的格子Boltzmann模型. 通过使
用Chapman-Enskog展开和多重尺度技术,得到了一系列的格子Boltzmann方程和时间
尺度$t_0$上守恒律,给出了满足地震压力波方程所要求的高阶矩以及简单的平衡态分布
函数表达式. 数值结果表明这种方法可以用来模拟地震压力波. 相似文献
7.
基于BGK碰撞模型,通过在迁移方程中引入作用力项,建立了粘弹流体的轴对称格子Boltzmann模型.通过Chapman-Enskog展开,获得了准确的柱坐标下轴对称宏观流动方程.采用双分布函数对运动方程和本构方程进行迭代求解,模拟分析了粘弹流体管道流动,获得了流场中的速度和构型张量的分布,通过与解析解进行比较,验证了模型的准确性.研究了作为粘弹流体流动基准问题的收敛流动,对涡旋位置进行了定量分析,将回转长度的计算结果与有限体积法进行了比较,两种数值结果十分吻合.研究结果表明,模型能够准确表征粘弹流体的轴对称流动,具有较广阔的应用前景. 相似文献
8.
9.
运用一种改进的非结构化四边形格子法,对含孔正交各向异性板条受面内冲击拉伸时弹性应力波的传播过程和孔边的动应力集中进行了研究.非结构化格子法采用与有限元类似的网格剖分方法,并基于围绕每个节点的积分平衡方程,并自然满足复杂边界的自由边界条件.计算中不需存储刚度矩阵,因而计算速度快、效率高、节省内存,在解决应力波传播问题中具有显著的优越性.通过对多种工况进行数值模拟,分析了材料的各向异性性质、纤维方向、孔径比、加载脉冲周期等参数对孔边动应力的影响,得到了一些规律性的结果.并与现有实验结果进行对比,验证了该方法的有效性. 相似文献
10.
11.
12.
研究了计入Peierls-Nabarro (P-N)力和固体黏性效应的一维金属杆在简谐外力扰动下的动力响应,其位移波的运动规律是Sine-Gordon (SG) 型方程. 采用集结坐标 (collective coordinate)将方程的解设为未扰系统呼吸子解的形式,研究扰动作用下,组成呼吸子的扭结-反扭结波的中心的分离. 通过用集结坐标表示系统的哈密顿量,从而将SG型方程转化为常微分方程组. 分析了未扰系统的异宿轨道,并将之用于Melnikov方法对系统进行分析,给出横截异宿点出现的必要条件,从而预测混沌运动的发生. 相似文献
13.
本文应用格子Boltzmann 方法(LBM)并结合Oldroyd-B 模型,讨论了不可压缩的 Navier-Stokes 方程和平流扩散本构方程的解耦及各自求解方法,以及两类问题的边界处理格式,实现了黏弹性流体在二维1:3 扩展流道以及3:1 收缩流道中的流动的数值模拟.获得了不同雷诺数Re 和维森伯格数Wi 以及黏度vs 下流动的流线分布,计算给出了漩涡的涡心位置和大小,并分析了参数Re、Wi 和vs 对流动特点的影响.模拟结果表明本文所采用模型和边界处理方法具有良好的精度和稳定性. 相似文献
14.
基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中. 相似文献
15.
基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中. 相似文献
16.
用格子Boltzmann方法研究Burgers方程 总被引:15,自引:0,他引:15
提出了用于Burgers方程的格子Boltzmann模型,应用Chapman-Enskog展开和多重尺度技术,通过选择平衡态分布函数的高阶矩,得出了几种精度的Burgers方程,模型中的参数通过分析耗散性质和色散给出。 相似文献
17.
流体的流动可以看成是分子以上水平的粒子基本运动组合而成,任何一个粒子系统的Hamiltonian都是由动能和势能这两部分所组成.借助于Hamiltonian建立了微观粒子和宏观流体之间的能量守恒准则,发展了一个适合于热流场数值模拟的格子Boltzmann模型.从该模型可以还原出宏观的流体力学方程,所得动量方程的黏性输运项除了具有Navier-Stokes黏性力的特征外还与非定常的、非线性的动量通量和非定常的内能相关.用该模型对Benard热对流进行了数值模拟,很好地再现了Benard cell,并且克服了热格子Boltzmann模型数值稳定性差的不足. 相似文献
18.
三维非均匀介质中弹性波传播的数值模拟 总被引:5,自引:1,他引:4
提出了一种三维非均匀介质中弹性波传播数值模拟的方法,文中称为三维格子法。该算法是二维格子法(一种二维非均匀介质中P-SV波传播的数值模拟算法)向三维非均匀介质情况的推广。在空间离散上该文方法与有限元方法类似,容许根据连续体的形状和介质分界面任意剖面网格,且自然满足自由表面边界条件。不同于常规有限差分法在各个节点上满足动力学微分方程,该算法通过满足各节点周围格子的整体平衡(积分平衡方程)来对问题进行求解,三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当。算例表明,该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射和绕射。 相似文献
19.
在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程.采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间.对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效.分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度.当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑.压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显. 相似文献