二维球流运动摩擦系数影响的数值模拟

格式：pdf
大小：263.48 KB
文档页数：5

下载文档原格式

超椭圆颗粒群绕流的二维数值模拟及结果分析

第 35 卷第 9 期 2019 年 9 月
科技通报
BULLETIN OF SCIENCE AND TECHNOLOGY
Vol．35 No．9 Sep． 2019
超椭圆颗粒群绕流的二维数值模拟及结果分析
张宁，曾庆敦
（华南理工大学力学系，广州 510640）
摘要：采用有限体积法和自适应网格技术，对超椭圆颗粒群在粘性、不可压流场中的绕流过程进行了数值模拟。在 3 ≤ Ｒe ≤ 200，0. 303≤ k ≤ 0. 836 and 2 ≤ n ≤ 5 的工况下，研究了雷诺数（Ｒe）、孔隙率（ k）和形状指数（ n）对颗粒群绕流过程中的流动特性（流线分布和速度矢量）、阻力特性（阻力系数 Cd ）和传热特性（平均努塞尔数 Nu）的影响。结果表明，随着Ｒe 从 3 增大至 200，超椭圆颗粒群中的 Cd 逐渐减小，而 Nu 逐渐增大。同样，随着 k 的增大， Cd 先急剧减小，当 k 超过某一临界值时，Cd 趋于稳定，而 Nu 逐始终随 k 的增大而增大。此外，无论 n 如何变化，e 始终与 k 呈正相关，并在 n = 2 时达到最大值。关键词：自适应网格；超椭圆颗粒群；阻力系数；努塞尔数
中图分类号： O357. 1
文献标识码： A
DOI： 10. 13774 / j．cnki．kjtb．2019. 09. 002
文章编号： 1001－7119（ -dimensional Numerical Simulation of Flow Past Multi-superelliptic Particles
流体流经由多个颗粒堆叠而成的颗粒群现象，广泛存在于自然界及工业生产中，并在化工工程、能源和环保等诸多领域具有重大的研究价值［1］。了

球底筒形件表面摩擦系数对拉深成形影响的数值模拟

球底筒形件表面摩擦系数对拉深成形影响的数值模拟杨守军;陈政弘;王红彬;赵福成;王瑞平【摘要】通过ABAQUS软件对球底筒形件拉深过程进行数值模拟,得到球底筒形件模具表面摩擦系数对成形件厚度和材料流动性具有重要影响,且模具表面存在最优的摩擦系数分布.基于此,采用均匀试验,对模具表面的摩擦系数进行优化组合.模具表面摩擦系数优化组合后,成形件的最大减薄率降低了26.54％,成形件的厚度分布更加均匀,均匀度提高24.87％;模拟结果为模具表面摩擦系数的主动优化设计制造提供理论依据.【期刊名称】《锻压装备与制造技术》【年(卷),期】2014(049)006【总页数】4页(P62-65)【关键词】冲压成形;筒形件;拉深;数值模拟;摩擦系数;优化设计【作者】杨守军;陈政弘;王红彬;赵福成;王瑞平【作者单位】宁波吉利罗佑发动机零部件有限公司,浙江宁波315336;宁波吉利罗佑发动机零部件有限公司,浙江宁波315336;宁波吉利罗佑发动机零部件有限公司,浙江宁波315336;宁波吉利罗佑发动机零部件有限公司,浙江宁波315336;宁波吉利罗佑发动机零部件有限公司,浙江宁波315336;浙江吉利罗佑发动机有限公司,浙江宁波315800【正文语种】中文【中图分类】TG386.3+21 引言金属板材冲压技术被广泛应用于国民生产的各行各业，拉深是最重要的成形工序之一，被国内外学者所关注。

拉深成形件常见的失效形式有起皱和拉裂，为了降低成形件的失效率，科研工作者做出了很多尝试和努力，但结果仍不能满足现实的要求。

经研究发现，对板料变形有影响的诸多因素中，摩擦和润滑状态占有举足轻重的地位。

摩擦显著影响着成形件的厚度和应变分布。

接触表面间的有效的润滑能够提高板料的成形性能和成形件的质量，另外还能够减少模具的磨损。

已有许多学者开始研究摩擦在金属板料成形过程中的作用。

但是绝大多数数值模型中，摩擦系数均设为常数，并没有考虑模具表面与板料之间摩擦系数随拉深过程的变化。

二维流体流动行为的数值模拟与分析

二维流体流动行为的数值模拟与分析引言流体力学是研究流体（包括气体和液体）在空间中运动的力学分支。

对于复杂的流体流动问题，数值模拟成为了解决问题的重要方法之一。

本文将针对二维流体流动行为展开数值模拟与分析的研究。

数值模拟方法在二维流体流动模拟中，常用的数值模拟方法包括有限差分法、有限体积法和有限元法等。

这些方法本质上都是将流动区域分割为离散的网格，然后利用数学算法进行计算。

不同的方法有不同的适用范围和精度，根据实际问题的特点选择合适的数值模拟方法十分重要。

有限差分法有限差分法是数值模拟中最常用的方法之一。

它将流动区域离散化成网格，然后利用差分近似方法求解偏微分方程。

有限差分法的主要优点是简单易懂，计算速度较快。

但是在处理复杂边界条件和非结构化网格时存在一定困难。

有限体积法有限体积法是一种将流动区域离散化为有限个控制体积的方法。

它将流体流动问题转化为在控制体积内的质量、动量和能量守恒方程求解问题。

有限体积法适用范围广，特别适合处理复杂边界条件和非结构化网格。

有限元法有限元法是一种将流动区域划分为无数小单元的方法。

每个小单元内的流场参数通过一组基函数在整个流动区域内插值得到。

有限元法适用于处理非线性问题和接触问题，但计算量较大。

数值模拟案例管道内的流动考虑一根长度为L的水平管道内的流动问题。

假设流体为牛顿流体，流动速度不变，管道内壁无滑移条件。

我们可以采用有限差分法对该问题进行数值模拟。

模拟设置：•管道长度L = 2m•管道内径d = 0.1m•流体动力粘度ν = 1.0e-6 m2/s•入口压力P_in = 2MPa数值模拟步骤：1.确定网格划分：将管道长度L等分为n个小段，确定网格划分尺寸Δx。

2.初始化流场参数：设置初始速度和压力分布。

3.迭代求解：利用差分近似方法求解控制方程，迭代得到收敛的流场解。

4.分析结果：根据模拟结果分析流体在管道内的流动行为。

翼型的气动性能研究翼型的气动性能是航空航天领域的重要课题之一。

自由流场二维圆柱阻力系数数值模拟

基于有限体积方法fluent将计算区域分割成很多不重复的控制体积每个网格点周围有一个控制体积在控制体积上对待解的微分方程积分得出一组离散方程未知数为网格点上的因变量加上边界条件和初始条件求解方程组中的未知量针对不同的收敛准则利用fluent软件进行模拟计算的一般流程如下定义边界类型和生成网然后将gamit中的网格输出文件导入fluent求解器设置边界条件和初始条件对流动区域进行求解最后对计算结果进行后处理fluent的前处理软件主要作用是建立网格模可以建立各种几何体划分几何体为很小的网格通过设立边界条件fluen它也可以导入在cadsolidwork模软件中生成的复杂图形然后对其进行网格划分建立网格模型在二维计算中圆柱用圆来代替圆的直径为cm

黏性流体绕流圆柱时产生两种力, 一种是圆柱边壁引起的边壁剪切应力, 另一种是圆柱体前后压差引起的正应力。这两种力的合力在入流方向上的分量即常说的阻力, 在垂直于入流方向上的分量即升阻力 [ 1]。
圆柱在黏性不可压流体流场中的阻力系数随雷诺数 (R e) 的变化而变化 [ 2] 。当 Re < 1时, 阻力系数很大; 当 1 < Re < 3 @ 105 时, 阻力系数随 R e增大而近乎单调减小。
当 0 < R e < 4时, 流线基本对称, 特征是在任意流体单元, 因压力差异而产生的力和摩擦力几乎相等, 流态如图 1( a) 所示; 当 4 < R e < 40时, 在圆柱后出现层流分流, 即形成两个稳定的方向相反的驻涡, 如图 1( b) 所示; 当 Re > 40时, 圆柱后部的驻涡开始不稳定地摆动、脱落, R e达到 100 左右时这对不稳定的对称旋涡分裂, 最后形成几乎稳定的、非对称的、规则的、旋转方向相反的交替旋涡, 如图 1( c) 所示; 当雷诺数再增大时, 流动分离严重, 约从 R e = 10 000起, 层流边界层从圆柱的前部距离停滞点约 80b位置处开始分离, 流动变为湍流, 形成很宽的分离区, 如图 1( d) 所示; 当雷诺数超过 100 000后, 层流边界层变为湍流边界层, 分离点向后推移 , 距离停滞点约 120b 位置处, 分离区变小, 如图 1( e) 所示 [3] 。

半球形件模具表面摩擦特性影响的数值模拟与优化设计_符永宏-诚也相关 - 江西是关于测量摩擦系数的方法设计

向压应力超过材料临界起皱压应力时，板料就会失稳起皱。防外皱其实是在特定条件下如何施加压边力只要施加适当的压边力，就可防止外皱的产的问题，生，因此内皱就成为研究半球形零件成形的关键。半球形件成形过程复杂，成形质量受诸多因素影响，如压边力、摩擦特性、原始毛坯的尺寸与形状、板料材料性能以及凹模圆角半径等，且这些因素之间有相关性，并不是彼此独立的，这就给研究带来了一定的难度。国内外的有关学者对半球形件拉深成形进行了大量的研究，取得了一定的研究成果。相关文献指出了板料上最易发生起皱的部位为原始板料上靠近凹模口内侧的金属，称此点为起皱关键点，若起皱关键点在成形过程中不起皱，则整个毛坯就不会发生起皱；板料上最易发生破裂的部位为板料上处于平面应变状态的点，这些点成为破
不受模具约束，当板料所受的纬板料处于悬空状态，
引言
半球形件是典型的曲面零件，因此研究半球形件的拉深成形过程具有一定的普遍意义。半球形件拉深成形是拉深变形和胀形变形的复合，其成形机理决定了球形零件的成形缺陷。半球形零件拉深成形的失效形式有起皱和拉裂，其中起皱包括法兰区。由于悬空区的起皱（外皱）和悬空区的起皱（内皱）

两相流对摩擦片间流体传动特性影响的数值模拟

LI Xiaojian，ZHANG Hong （ College of Mechanical Engineering，Taiyuan University of Technology，Taiyuan Shanxi 030024，China ）
Abstract： In order to study the effect of two phases of hydro-viscous drive on the fluid transmission characteristics of the friction film，also to reduce the two phase flow and improve the ability of the oil film to transfer shear torque，the flow field of oil film based on the VOF model was simulated numerically by fluid dynamics simulation software． The production and development rules of two phase flow regions in oil film flow field were gotten． The results show that the oil groove is the main area that produces two phase flow，and the air occupies a volume ratio of 1 to the side transition position of the oil groove near the outer diameter． What’s more，the production of two phase flow regions is mainly influenced by fluid pressure，and its production starts from near the outer diameter and gradually develops towards the inner diameter． The volume of air in the oil film flow field decreases with the increase of oil supply，but increases with the increase of rotational speed． The oil film transmitted torque increases with the increase of oil supply and input speed．

流场中互相靠近的两物体的水动力特性的数值模拟研究

关键词：圆柱；水动力特性；动网格；相对运动；计算流体力学
ABSTRACT
In the flow field, when two or more objects in close proximity are in relative motion, they will interfere with each other, which is widespread in the enginering. To solve these problems, it is necessary to deeply research the characteristics of the interaction among two or more objects in the flow field. In this dissertation, the hydrodynamic characteristics of two objects approaching to each other in viscous flow were studied using the Computational Fluid Dynamics (CFD) tool based on dynamic mesh technology and some results are obtained. In the two-dimensional (2D) viscid flow, the cylinder force characteristics in the viscid flow were firstly compared with that in the inviscid flow. Then, the characteristics of forces acting on the cylinders in viscid flow under different interstice ratios, different approaching velocities and different background flows were analyzed. The numerical results show that when the background flow velocity is zero, the drag coefficient increases with the decrease of interstice ratio and the increase of approaching velocity. As the background flow velocity is not zero, the drag coefficient increases with the increase of background flow velocity. While the ratio of the background flow velocity to the approaching velocity is large, the force acting on the cylinder vibrates obviously. Further on, the force of the two objects close to each other in the three-dimensional (3D) flow field was analyzed based on the method in the 2D flow field. Compared with the 2D case, the repulsion significantly reduced when the cylinders approached to each other in the 3D case, which shows that the wall effect of cylinders approaching to each other in the 3D case is less than that in the 2D case. We also verify that the calculation method used to simulate rotation ellipsoid problem is reliable. When simulating the situation that underwater ellipsoid was close to the structures and performing rotational movement, we find that the torque exerted on the ellipsoid changes greatly when the ellipsoid rotates around itself and the structures. It should be noticed when the ellipsoid is controlled to perform the precise motion.

二维不可压缩粘性流体绕钝体流动的数值模拟

Ｕ／，其中Ｄ为圆柱直径；Ｄ￣，Ｕ为自由来流速度，过通改变来流速度来得到不同的雷诺数。３２计算域及网格划分．
数时，光滑圆柱体的周期性尾流只是雷诺数的函数。按圆柱体直径计算的雷诺数很小时，体并不脱离流
界区，此时柱体表面上的边界层为层流，而柱体后面
的涡街已完全转变为湍流，并按一定的频率发放漩涡； × ０＜ × ０称为过渡区，３１３１。此时柱体表面上的
边界层也已变为湍流，离点向后移，力显著下分阻
符合重新分离，这一类物体的绕流在较大的雷诺数范围内变化不明显。
３计算实例
３１控制方程．
通过控制流体由计算域的流速来控制雷诺数，
得到的尾迹及漩涡脱落图如图３所示。当Ｒｅ＝１，流中有一对稳定的弗普尔旋时尾
第８期
曹广龙等：二维不可压缩粘性流体绕钝体流动的数值模拟
Ｃ＝ｐ竿堕
式中：和为来流的静压和速度；底部压强。Ｐｐ为
对于不可压缩粘性流体，直角坐标系下，在其运动规律受Ｎ—Ｓ方程控制，续性方程和动量方程连
分别为：
Ｏｔｕ
＿
：
对于三维钝体绕流，压差阻力仍然是总阻力的主要部分。若对物体形状进行流线型处理；压差阻力则
定义是：
．＝ｓ
决定圆柱绕流流态的是雷诺数（ｅ的值，Ｒ）当

二维圆柱绕流数值仿真模拟报告

二维圆柱绕流数值仿真模拟报告
二维圆柱绕流数值仿真模拟是流体力学领域的一个重要研究课题。

在这个问题中，我们将从多个角度来讨论这个课题，包括数值
模拟的背景和意义、数值方法、模拟结果及其分析等方面。

首先，让我们来谈谈数值模拟的背景和意义。

在工程和科学研
究中，流体力学的数值模拟在许多领域都发挥着重要作用，例如飞
行器设计、汽车空气动力学、海洋工程等。

而二维圆柱绕流作为流
体力学中经典的基本问题，对于理解流体运动、探索流体特性具有
重要意义。

通过数值仿真模拟，我们可以更好地理解圆柱绕流的流
场特性，为相关工程和科学问题提供重要参考。

其次，数值方法是进行二维圆柱绕流数值仿真模拟的关键。

常
用的数值方法包括有限元法、有限体积法和边界元法等。

这些方法
可以通过数值计算求解流体力学方程，如Navier-Stokes方程，从
而得到圆柱绕流的流场分布、压力分布等信息。

在选择数值方法时，需要考虑计算精度、计算效率以及模拟的适用范围等因素。

接下来，我们来讨论模拟结果及其分析。

通过数值仿真模拟，
我们可以得到圆柱绕流的流速场、压力分布、升力和阻力等重要参
数。

通过分析这些参数，我们可以深入理解圆柱绕流的特性，比如卡门涡的形成与演变、升力和阻力的变化规律等。

这些分析结果对于优化工程设计、改进流体力学理论具有重要意义。

总的来说，二维圆柱绕流数值仿真模拟是一个复杂而有意义的课题，通过深入研究和全面分析，可以为相关领域的工程和科学问题提供重要的参考和指导。

希望以上讨论能够对你有所帮助。

基于统计分析和数值模拟方法的滑坡视摩擦系数影响因素及特征研究

０引言
滑坡的最大垂直运动距离Ｈ和最大水平运动距离Ｌ是滑坡运动学研究的重要参数。前者能决定滑坡的总能量，后者是滑坡灾害防灾减灾重点考虑的因素。Ｈ／Ｌ被称为视摩擦系数，其能够衡量滑坡的运动能力（Ｊｏｈｎｓｏｎｅｔａｌ．，２０１７）。不少学者研究了坡脚角度、滑面长度、偏转角等对滑坡运动的影响，詹威威等（２０１７）选取了３８处汶川地震中触发的地震滑坡通过逐步回归方法建立了滑坡－碎屑流最大水平运动距离Ｌ的最优多元回归模型；杨海龙等（２０１９）利用９０处地震滑坡，基于主要地形参数和滑坡规模建立了最大水平运动距离的预测统计模型；郑光等（２０１９）考虑滑坡碎屑流运动距离与势能的关系，得到远程运动距离Ｌ和势能ＶＨ间的方程；汪发武等（２００２）分析了土粒破碎特性在高速滑坡运动过程的作用，结合视摩擦系数模型和滑坡运动量方程及连续性方程，对滑坡的运动范围进行了预测。
ＩＮＦＬＵＥＮＣＥＦＡＣＴＯＲＳＡＮＤＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＡＰＰＡＲＥＮＴＦＲＩＣＴＩＯＮＣＯＥＦＦＩＣＩＥＮＴＯＦＬＡＮＤＳＬＩＤＥＢＡＳＥＤＯＮＳＴＡＴＩＳＴＩＣＡＬＡＮＡＬＹＳＩＳＡＮＤＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮ
ＷＡＮＧＲａｎ①②③ ＷＡＮＧＸｕｅｌｉａｎｇ①②③ ＹＵＡＮＨｏｎｇｈｕ④ ＳＵＮＪｕａｎｊｕａｎ①②③ ＬＩＵＨａｉｙａｎｇ①②③ Ｚ最大水平运动距离Ｌ的比值（Ｈ／Ｌ）能很好地衡量滑坡的运动能力。作者采用数理统计方法分析了藏东南区滑坡及地震型滑坡中影响因素（滑坡体积Ｖ、斜坡角度 θ、滑坡最大垂直运动距离Ｈ）及试验因素（Ｈ／Ｌ值）的水平特征，进一步分析了３个影响因素对该试验因素的影响特征，并采用基于连续介质力学法的数值模拟方法，模拟了滑坡运动的时空演化过程。得到以下结论：θ是Ｈ／Ｌ值的显著影响因素且主要集中在２０°～４０°。随着 θ的减小、Ｖ的增大和Ｈ的减小，Ｈ／Ｌ具有减小的趋势。藏东南滑坡Ｈ／Ｌ值普遍大于地震型滑坡Ｈ／Ｌ值，地震型滑坡的远程特征显著，地震作用会弱化 θ、Ｖ、Ｈ与Ｈ／Ｌ值的相关程度。在坡脚附近滑坡的运动速度往往达到峰值，同时坡脚对滑坡的运动具有显著的阻止效应。研究结果对于滑坡运动特征研究及防灾减灾具有一定参考意义。关键词滑坡；视摩擦系数；统计分析；连续介质力学法；数值模拟中图分类号：Ｐ６４２２２文献标识码：Ａｄｏｉ：１０．１３５４４／ｊ．ｃｎｋｉ．ｊｅｇ．２０２０－４４９

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

关键词：球床式高温气冷堆；球流运动；离散单元法；摩擦系数
中图分类号：ＴＬ３３４；ＴＬ４２４文献标志码：Ａ文章编号：１０００ — ６９３１（２０１３）０７ — １１４１一Ｏ５
ｄｏｉ：１０．７５３８／ｙｚｋ．２Ｏ１３．４７．０７．１１４ｌ
从表中可看出随着球摩擦系数的增大球流流场内速度分布更加均匀处于对称中心附近的球流流速变化不大而靠表1不同球摩擦系数下统计平均时间table1meanresidencetimeatdifferentballfrictions球摩擦系数统计平均时间单位时间标记点0标记点1标记点2标记点301534856428561213850254125718144202900453085626754518877近壁面的流速越来越快
ｓｅｅｍｓｔｏｂｅｍｏｒｅｕｎｉｆｏｒｍ．Ｔｈｅｗａｌ１ｆｒｉｃｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｈａｓ１ｉｔｔｌｅｅｆｆｅｃｔｏｎｔｈｅｈｏｒｉｚｏｎｔａｌ
动进行了数值模拟，并对从标记点产生的球流区域的均值流线、标准方差和平均停留时间等进行了比较
和分析。结果表明：球的摩擦系数对球流流场影响不大；球的摩擦系数越大，水平方向扩散越小，流动越均匀；壁面摩擦系数对水平方向扩散影响不大；壁面摩擦系数越大，流动越不均匀。
第４７卷第７５１
２０１３年７月
原
子
能
Байду номын сангаас
科
学
技
术
Ｖｏ１．４７，ＮＯ．７
ＡｔｏｍｉｃＥｎｅｒｇｙＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ
Ｊｕ１．２０１３
维球流运动摩擦系数影响的数值模拟
ＷＡＮＧＷｅｎ — ｊｕｎ，ＹＡＮＧＸｉｎｇ — ｔｕａｎ，ＪＩＡＮＧＳｈｅｎｇ — ｙａｏ，
ＴＡＮＧＺｈｅｎ — ｗｅｉ，ＬＩＵＺｈｉ — ｙｏｎｇ
（ＫｅｙＬａｂｏｒａｔｏｒｙｏｆＡｄｖａｎｃｅｄＲｅａｃｔｏｒＥｎｇｉｎｅｅｒｉｎｇａｎｄＳａｆｅｔｙｏｆＭｉｎｉｓｔｒｙｏｆＥｄｕｃａｔｉｏｎ，ＩｎｓｔｉｔｕｔｅｏｆＮｕｃｌｅａｒａｎｄＮｅｗＥｎｅｒｇｙＴｅｃｈｎｏｌｏｇｙ，ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ１０００８４，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｉｎｏｒｄｅｒｔｏｉｎｖｅｓｔｉｇａｔｅｔｈｅｐｅｂｂｌｅｆｌｏｗｄｙｎａｍｉｃｓｉｎｔｈｅｈｉｇｈ— ｔｅｍｐｅｒａｔｕｒｅｒｅａｃｔｏｒｃｏｒｅａｎｄｂａｓｅｄｏｎｔｈｅｔｗｏ — ｄｉｍｅｎｓｉｏｎａｌｅｘｐｅｒｉｍｅｎｔｓｏｆｐｅｂｂｌｅｆｌｏｗｄｙｎａｍｉｃｓ，ｄｉｓｃｒｅｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＤＥＭ）ｗａｓｕｓｅｄｔｏｓｉｍｕｌａｔｅｔｈｅｐｅｂｂｌｅｆｌｏｗｄｙｎａｍｉｃｓ．Ｔｈｅｍｅａｎｆｌｏｗｓｔｒｅａｍｌｉｎｅｓ，ｓｔａｎｄａｒｄｄｅｖｉａｔｉｏｎａｎｄｔｈｅｍｅａｎｒｅｓｉｄｅｎｃｅｔｉｍｅｏｆｔｈｅｐｅｂｂｌｅｆｌｏｗｚｏｎｅｇｅｎｅｒａｔｅｄｂｙｍａｒｋｅｒｓｗｅｒｅｃｏｍｐａｒｅｄａｎｄａｎａｌｙｚｅｄ．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｔｈａｔｂａｌｌｆｒｉｃｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｈａｓｌｉｔｔｌｅｅｆｆｅｃｔｏｎｔｈｅｐｅｂｂｌｅｆｌｏｗｆｉｅｌｄ．Ｗｉｔｈｔｈｅｐｅｂｂｌｅｆｒｉｃｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｉｎｃｒｅａｓｉｎｇ，ｔｈｅｈｏｒｉｚｏｎｔａｌｄｉｆｆｕｓｉｏｎｏｆｐｅｂｂｌｅｓｄｅｃｒｅａｓｅｓａｎｄｔｈｅｐｅｂｂｌｅｆｌｏｗ
ＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎｏｎＦｒｉｃｔｉｏｎＣｏｅｆｆｉｃｉｅｎｔＥｆｆｅｃｔ
ｏｆＰｅｂｂｌｅＦｌｏｗＤｙｎａｍｉｃｓｉｎＴｗｏ ‘ －ｄｉｍｅｎｓｉｏｎａｌＰｅｂｂｌｅ ‘ －ｂｅｄＲｅａｃｔｏｒ
王文俊，杨星团，姜胜耀，汤振伟，刘志勇
（清华大学核能与新能源技术研究院先进反应堆与安全教育部重点实验室，北京１０００８４）
二
摘要：为研究摩擦系数对球流运动的影响，参照二维堆芯球流运动实验台架，采用离散单元法对球流运