fluent中多孔粘性阻力与内部阻力的计算方法

Experimental data that is available in the form of pressure drop against velocity through

the porous component, can be extrapolated to determine the coecients for the porous

media. To e ect a pressure drop across a porous medium of thickness, n, the coecients

of the porous media are determined in the manner described below.If the experimental data is:

Velocity Pressure Drop

(m/s) (Pa)

20.0 78.0

50.0 487.0

80.0 1432.0

110.0 2964.0

多孔介质在fluent中的操作方法 网络上传版本

如何在Fluent中实现多孔介质双能量方程(LNTE) How to use Non-equilibrium Thermal equation (LNTE) model for Porous media in Fluent Software ●请参照本人发表的文章： ●Please refer to the following papers： 1)Wang Fu–Qiang*，Shuai Yong*，Wang Zhi–Q iang，Leng Yu，Tan He–Ping. Thermal and chemical reaction performance analyses of steam methane reforming in porous media solar thermochemical reactor，International Journal of Hydrogen Energy，39(2)：718-730，2014 关键词：Porous, Solar, Hydrogen, Methane, Reforming, P1 approximation, radiative heat transfer 2)Wang Fu–Qiang*，Shuai Yong*，Tan He–Ping，Zhang Xiao-Feng，Mao Qian-Jun，Heat transfer analyses of porous media receiver with multi–dish collector by coupling MCRT and FVM method，Solar Energy，93：158–168，2013 关键词：Solar, Porous, dish concentrator, Receiver, Monte Carlo 3)Wang Fu–Qiang*，Shuai Yong*，Tan He–Ping，Yu Chun–Liang，Thermal Performance Analysis of Porous Media Receiver with Concentrated Solar Irradiation，International Journal of Heat and Mass Transfer，62：247–254，2013 关键词：Solar, Porous, dish concentrator, Receiver, Monte Carlo

基于fluent的阻力计算

基于fluent的兴波阻力计算本文主要研究内容 本文的工作主要涉及小型航行器在近水面航行时的绕流场及兴波模拟和阻力的数值模拟两个方面。在阅读大量文献资料的基础上，通过分析、比较上述领域所采用的理论和方法，针对目前需要解决的问题，选择合理的方法加以有机地综合运用。具体工作体现在以下几个方面： 1．本人利用FLUENT软件的前处理软件GAMBIT自主建立简单回转体潜器模型，利用FLUENT求解器进行计算，得出在不同潜深下潜器直线航行的绕流场、自由面形状及阻力系数的变化情况。 2．通过对比潜器在不同潜深情况下的阻力系数，论证了增加近水面小型航行器的深度可以有效降低阻力。通过对模型型线的改动，为近水面小型航行器的型线设计提供了一定的参考。通过改变附体形状和位置计算了附体对阻力的影响程度，为附体的优化设计提供了一定的依据。 计算模型

航行器粘性流场的数值计算理论 水动力计算数学模型的建立 根据流体运动时所遵循的物理定律，基于合理假设（连续介质假设）用定量的数学关系式表达其运动规律，这些表达式成为流体运动的数学模型，它们是对流体运动的一种定量模型化，称为流体运动控制方程组。根据控制方程组，结合预先给定的初始条件和边界条件，就可以求解反映流体运动的变量值，从而实现对流体运动的数值模拟预报，形成分析报告。 基于连续介质假设的流体力学中流体运动必须满足要遵循的物理定律: 1) 质量守恒定律 2）动量守恒定律 3）能量守恒定律 4）组分质量守恒方程 针对具体研究的问题，有选择的满足上述四个定律。船体的粘性不可压缩绕流运动，如果不考虑水温对水物理性质的影响，水的密度和分子粘性系数都是常数，同时没有能量的转换，就仅仅需要满足质量守恒定律、动量守恒定律。在满足这些定律下所建立的数学模型称为Navier-Stokes方程。 另外，自由液面的存在也需要建立合适的数学模型。本文是利用FLUENT 进行数值模拟，而软件里面关于自由液面模拟是用界面追踪方法的一种－流体体积法（VOF）,基于该方法所建立的数学模型称为流体体积分数方程。另外，高雷诺数下的水动力问题还需要考虑粘性不可压缩流体的湍流运动。对于湍流运动的数值模拟一直是流体力学数值计算的一个难点。直接数值模拟（DNS）目前还仅仅在院校中研究，而且也仅限于二维流体问题。大涡模拟（LES）向工程应用的过渡似乎还没有完成，并且就高雷诺数问题而言，对计算机硬件要求很苛刻。目前，从算法的可行性、硬件要求的可实现性、完成任务所消耗时间和人力等方面看，基于湍流模型的数值计算更为工程实际所接受。本章将会对各种湍流模型加以介绍。 粘性不可压缩流体流动数学模型 连续方程 任何流动问题都必须满足质量守恒方程即：连续方程。根据连续介质假设，单位时间内流体微团的质量变化等于同时间间隔内进入微团的总净质量。按照这一定律，连续方程数学表达式写为： （2.1） 以上是在笛卡尔直角坐标系下表示，上面给出的是瞬态可压流体连续方程。由于对于潜艇粘性流场介质的不可压缩，密度ρ为常数，引入散度算子，则方程（2.1）变成为： （2.2）式中：速度矢量V= { u ,v, w }。上式为粘性不可压缩流体运动的连续方程。 动量方程

fluent计算错误分析

1. FlUENT 1.1 求解方面 1.1.1 floating point error是什么意思？怎样避免它？ Floating point error已经提过很多次了并且也已经对它讨论了许多。下面是在Fluent论坛上的一些答案： 从数值计算方面看，计算机所执行的运算在计算机内是以浮点数（floating point number）来表示的。那些由于用户的非法数值计算或者所用计算机的限制所引起的错误称为floating point error。 1）非法运算：最简单的例子是使用Newton Raphson方法来求解f(x)=0的根时，如果执行第N次迭代时有，x=x(N)，f’(x(N))=0，那么根据公式x(N+1)=x(N)-f(x(N))/ f’(x(N))进行下一次迭代时就会出现被0除的错误。 2）上溢或下溢：这种错误是数据太大或太小造成的，数据太大称为上溢，太小称为下溢。这样的数据在计算机中不能被处理器的算术运算单元进行计算。 3）舍入错误：当对数据进行舍入时，一些重的数字会被丢失并且不可再恢复。例如，如果对0.1进行舍入取整，得到的值为0，如果再对它又进行计算就会导致错误。 避免方法 计算和迭代我认为设一个比较小的时间步长会比较好的。或者改成小的欠松驰因子也会比较好。从我的经验来看，我把欠松驰因子设为默认值的1/3；降低欠松驰因子或使用耦合隐式求解；改变欠松驰因子，如果是非稳态问题可能是时间步长太大；改善solver-control-limits 比例或许会有帮助；你需要降低Courant数；如果仍然有错误，不选择compute from初始化求解域，然后单击init。再选择你想从哪个面初始化并迭代，这样应该会起作用。另外一个原因可能是courant数太大，就样就是说两次迭代之间的时间步太大并且计算结果变化也较大（残差高）。 网格问题当我开始缩放网格时就会发生这个错误。在Gambit中，所有的尺寸都是以mm 为单位，在fluent按scale按钮把它转换成m，然后迭代几百次时就会发生这种错误。但是当我不把网格缩放到m时，让它和在Gambit中一样，迭代就会成功；我认为你应当检查网格，你的网格数太多了，使用较少的网格问题就会解决；网格太多，计算机资源不够用，使使比较粗的网格。 边界条件在我的分析中，我设了一个wall边界条件来代迭axis边界条件，结果fluent拒绝计算并告诉我floating point error。你的边界条件不能代表真实的物理现象；错误的边界条件定义可能会导致floating point error。例如把内边界设成interior；一次我使用对称边界条件模拟2D区间时也遇到这种问题，我把symmetry设为axe symmetric，就发生了floating point error；检查你设的湍流参数，减小湍流强度，先进行50次迭代。 多处理器问题我近来在进行多处理器模拟时也遇到相似的问题。问题的解决方法是在单个处理器上运行，这样就运算得很好。 错误迭代以错误的条件来初始化，在开始迭代时就会发生floating point error。 1.1.2 coupled和segregated求解有什么区别？ Coupled会同时求解所有的方程（质量守恒方程、动量守恒方程和能量守恒方程）而不是单个方程求解（方程互相分离）。当速度和压力高度耦合（高压和高速）时应该使用耦合求解，但这样会需要较长的计算时间。 在耦合求解中，能量方程中总是包含组分扩散（Species Diffusion Term）项。

3_流体流动时摩擦阻力系数的测定

生物系统传输过程实验报告 实验报告 课程名称：生物系统传输过程实验 指导老师：叶章颖 成绩：______ 实验名称：流体流动时摩擦阻力系数的测定 实验类型：__探究型实验__ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 一、实验目的 测定流体流动时的沿程阻力系数和局部阻力系数及不同流型下直管沿程阻力系数λ随雷诺数Re 的变化关系。掌握流量压强的测量方法。 二、实验装置 1.实验设备的基本情况 实验流程示意图见图1。 水泵2将储水槽1中的水抽出，送入实验系统，首先经玻璃转子流量计15、16测量流量，然后送入被测直管段测量流体在光滑管或粗糙管的流动阻力，或经10测量局部阻力后回到储水槽，水循环使用。被测直管段流体流动阻力△p 可根据其数值大小分别采用变送器12或空气-水倒置∪型管22来测量。 图1 流动阻力实验流程示意图 1-水箱；2-离心泵；3、4-放水阀；5、13-缓冲罐；6-局部阻力近端测压阀；7、15-局部阻力远端测压阀；8、20-粗糙管测压回水阀；9、19-光滑管测压回水阀；10-局部阻力管阀；11-U 型管进水阀；12-压力传感器；14-流量调节阀； 15、16-水转子流量计；17-光滑管阀；18-粗糙管阀； 21-倒置U 型管放空阀；22-倒置U 型管；23-水箱放水阀；24-放水阀； 专业：生物系统工程 姓名：邵建智 学号：3110100122 日期：2013.9.30 地点： 院楼D228

2.设备的主要技术数据 （1） 被测光滑直管段： 管径d —0.008m ； 管长L —1.69m ； 材料—不锈钢管 被测粗糙直管段： 管径 d —0.010m ； 管长L —1.69m ； 材料—不锈钢管 （2）被测局部阻力直管段: 管径 d —0.015m ；管长 L —1.2m ； 材料—不锈钢管 （3）压力传感器: 型号:LXWY 测量范围: 200 KPa （4）直流数字电压表: 测量范围: 0 ～ 200 KPa （5）离心泵: 型号: WB70/055 流量: 8(m 3／h) 扬程: 12(m) 电机功率: 550(W) （6）玻璃转子流量计: 型号 测量范围 精度 LZB —40 100～1000(L ／h) 1.5 LZB —10 10～100(L ／h) 2.5 3.实验设备的功能与特点 本实验装置可用于实验教学和科研。利用该实验装置，可学习和掌握光滑直管、粗糙直管的阻力系数与雷诺准数的测量方法；也可学习局部阻力的测量方法；学习几种压差测量方法；加深对流体流动阻力概念的理解。 本实验装置的特点： ⑴ 本实验装置数据稳定，重现性好，能给实验者明确的流体流动阻力概念。 ⑵ 雷诺准数的数据范围宽，可作出102～104三个数量级。能够测量出光滑管、粗糙管的阻力系数与雷诺准数的关系，同时也可以测量阀门局部阻力。 ⑶ 实验采用循环水系统，节约实验费用。 ⑷ 测量系统采用量程不同的两种流量计和压差测量仪表，测量精度较高。 ⑸ 采用压力传感器—数字表系统,测量大流量下的流体流动阻力, 实验数据稳定可靠。 三、实验原理 1. 直管摩擦系数λ与雷诺数Re 的测定 直管的摩擦阻力系数是雷诺数和相对粗糙度的函数，即)/(Re,d f ελ=，对一定的相对粗糙度而言，(Re)f =λ。 流体在一定长度等直径的水平圆管内流动时，其管路阻力引起的能量损失为： g P g P P h f f ρρ?=-=21 （1） 又因为摩擦阻力系数与阻力损失之间有如下关系（范宁公式） g u d l h g P f f 22λρ== ? （2） 整理（1）（2）两式得 2 2u P l d f ???=ρλ （3）

Fluent性能分析

Fluent性能分析 仅仅就我接触过得谈谈对fluent的认识，并说说哪些用户适合用，哪些不适合 fluent对我来说最麻烦的不在里面的设置，因为我本身解决的就是高速流动可压缩N-S方程，而且本人也是学力学的，诸如边界条件设置等概念还是非常清楚的同时我接触的流场模拟，都不会有很特别的介质，所以设置起来很简单 对我来说，颇费周折的是gambit做图和生成网格，并不是我不会，而是gambit对作图要求的条件很苛刻，也就是说，稍有不甚，就前功尽弃，当然对于计算流场很简单的用户，这不是问题。有时候好几天生成不了的图形，突然就搞定了，逐渐我也总结了一点经验，就是要注意一些小的拐角地方的图形，有时候做布尔运算在图形吻合的地方，容易产生一些小的面最终将导致无法在此生成网格， fluent里面的计算方法是有限体积法，而且我觉得它在计算过程中为了加快收敛速度，采取了交错网格，这样，计算精度就不会很高。同时由于非结构网格，肯定会导致计算精度的下降，所以我一贯来认为在fluent里面选取复杂的粘性模型和高精度的 格式没有任何意义，除非你的网格做的非常好。 而且fluent5.5以前的版本（包括5。5），其物理模型，（比如粘性流体的几个模型）都是预先设定的，所以，对于那些做探索性或者检验新方法而进行的模拟，就不适合 用。 同时gambit做网格，对于粘性流体，特别是计算湍流尺度，或者做热流计算来说其网格精度一般是不可能满足的，除非是很小的计算区域 所以，用fluent做的比较复杂一点的流场（除了经典的几个基本流场） 其计算所得热流，湍流，以及用雷诺应力模拟的粘性都不可能是准确的， 这在物理上和计算方法已经给fluent判了死刑，有时候看到很多这样讨论的文章，觉得 大家应该从物理和力学的本质上考虑问题。 但是，fluent往往能计算出量级差不多的结果，我曾经做了一个复杂的飞行器热流计算，高超音速流场，得到的壁面热流，居然在量级上是吻合的，但是，从计算热流需要的壁面网格精度来判断，gambit所做的网格比起壁面网格所满足的尺寸的要大了至少2个数量级， 我到现在还不明白fluent是怎么搞的。 综上，我觉得，如果对付老板的一些工程项目，可以用fluent对付过去但是如果真的做论文，或者需要发表文章，除非是做一些技术性工作，比如优化计算 一般用fluent是不适合的。 我感觉fluent做力的计算是很不错的，做流场结构的计算，即使得出一些涡也不是流场本身性质的反应，做低速流场计算，fluent的优势在于收敛 速度快，但是低速流场计算，其大多数的着眼点在于对流场结构的探索，所以计算得到的结果就要好好斟酌一下了，高速流场的模拟中，一般着眼点在于气动力的结果，压力分布

阻力损失的计算方法

1.5阻力损失 1.5.1两种阻力损失 直管阻力和局部阻力 化工管路主要由两部分组成：一种是直管，另一种是弯头、三通、阀门等各种管件。 直管造成的机械能损失称为直管阻力损失（或称沿程阻力损失） 管件造成的机械能损失称为局部阻力 注意 将直管阻力损失与固体表面间的摩擦损失相区别 阻力损失表现为流体势能的降低 由机械能衡算式（1-42）可知： ρρρ212211P P g z p g z p h f -=??? ? ??+-???? ??+= (1-71) 层流时直管阻力损失 流体在直管中作层流流动时，因阻力损失造成的势能差可直接由式（1-68）求出： 232d lu μ?= ? (1-72) 此式称为泊稷叶（Poiseuille)方程。层流阻力损失遂为： 232d lu h f ρμ= (1-73) 1.5.2湍流时直管阻力损失的实验研究方法 实验研究的基本步骤如下： （1）析因实验－寻找影响过程的主要因素

对所研究的过程作初步的实验和经验的归纳，尽可能的列出影响过程的主要因素。对湍流时直管阻力损失f h ，经分析和初步实验获知诸影响因素为： 流体性质：密度ρ、粘度μ； 流动的几何尺寸：管径d 、管长l 、管壁粗糙度ε（管内壁表面高低不平）： 流动条件：流速u 。 于是待求的关系式为： ) ,,,,,(ερμu l d f h f = (1-74) （2）规划实验－减少实验工作量 因次分析法的基础是：任何物理方程的等式两边或方程中的每一项均具有相同的因次，此称为因次和谐或因次的一致性。 以层流时的阻力损失计算式为例，式（1-73）可写成如下形式 ???? ????? ??=??? ? ??dup d l u h f μ322 (1-75) 式中每一项都为无因次项，称为无因次数群。 换言之，未作无因次处理前，层流时阻力的函数形式为： ) ,,,,(u l d f h f ρμ= (1-76) 作无因次处理后，可写成

通风摩擦阻力系数

中华人民共和国煤炭工业部 矿井通风巷道摩擦阻力系数（a标）表 （试行） 主编部门：沈阳煤矿设计研究院 批准部门：煤炭工业部规划设计总院 试行日期：1985年1月1日 整理： 校核： 二ΟΟ三年一月

说明 1．井巷道通风摩擦阻力系数表，是我国自行实测的矿井巷道通风阻力系数，（除锚喷支护外其它各种支护巷道系验证测定）于1983年3月由煤炭工业部设计管理局主持召开了鉴定会，本表系根据鉴定会纪要精神，进行修改后，汇编而成。 2．表中摩擦阻力系数a标是标准状态下（t=20℃，P=760mmHg，ψ=60%）空气重率r=1.2kg ?/m3时的a值。 3．巷道类别划分原则，以支护特征、巷道壁面特征、巷道装备等与摩擦阻力系数相关的影响因素分类，不以巷道使用名称和进、回风道等分类。 4．表中凡是平巷的皆包含无行人台阶的倾斜巷道，凡是斜巷皆指设有行人台阶而言，通风行人巷为不铺轨的巷道，胶带输送机巷均铺设一条单轨轨道。 5．无轨道的锚喷胶带输送机巷道的a值，未能实测，暂可参照锚喷通风行人巷（无轨道、台阶）的a值与胶带机的附加a值综合选取。即光爆凸凹度<150mm，a=（10.9～17.6）×10-4；普爆凸凹度>150mm，a=(11.6～19.9)×10-4。 6．光面爆破与壁面凸凹度划分的标准以煤炭部制订的“煤矿井巷工程光面爆破、锚杆、喷浆、喷射混凝土支护施工试行规程”为准，普通爆破系指采用光面爆破的煤矿一般常用的爆破方法。 7．巷道壁面平滑与粗糙的划分标准，以粗糙度的平均突起高度为准。混凝土井巷壁面，壁面平滑的粗糙度平均突起高度为0.00025m，壁面粗糙的粗糙度平均突起高度为0.0007m，为测量和选取方便，将壁面经过抹光或粉刷的视为壁面平滑，壁面未经过抹光或未粉刷的视为壁面粗糙。 8．系数值的来源依据，除已注明资料出处之外的实测值，均可查找本资料的附件部分，以便于选取系数值时参考现场条件。 9．本表所给出的a值，应用时需要乘以10-4，并不需再考虑装有设备、台阶和工作面采煤机的a附加值。 10．经实测、资料统计提供各类的a附加值：装有胶带输送机的巷道,a附加值(4～10)×10-4；没有行人台阶的巷道，a附加值(1～3)×10-4；巷道堵塞较严重时，a附加值(3～10)×10-4；弯曲的巷道，a附加值(2～5)×10-4；巷道断面局部变化（单、双轨）a附加值3×10-4；铺轨无道渣填充的平巷a附加值(1～3)×10-4；工作面采煤机的a附加值(6～9)×10-4. 11．1mmH2O=9.80665Pa h摩=(a×L×U/S3)×Q2 =R×Q2

FLUENT 15.0 VOF模型测试报告

ANSYS 15.0 系列测试报告 FLUENT 15.0 VOF模型 测试人：崔亮安世亚太公司 测试时间：2013.12.01

1、仿真平台 HP Z820工作站，Intel Xeon E5-2690 * 2，内存64GB，2TB SATA硬盘。安装ANSYS 15.0 Preview3版本。 2、仿真模型 对某车型上带有底部隔板的油箱，在车辆加速时油箱内燃油晃动的瞬态过程进行瞬态仿真分析，网格单元数约10万，使用FLUENT的VOF模型计算空气和燃油的两相交界面。重点考察FLUENT 15.0中VOF模型的计算效率和两相交界面捕捉精度的提升。 测试案例的几何形状 测试案例的网格模型 3、试用情况 1）.稳定性 在整个试用过程中，软件保持稳定，未出现任何不流畅、死机、系统崩溃等情况。2）.流畅度 模型拖动、旋转、缩放等操作十分流畅，模型设定及求解过程操作十分流畅。 3）.效率 该模型使用0.0005秒的时间步长进行瞬态计算，共计算了2000步，共计1.0秒时长。使用15.0 Preview3版本所用的计算时间为3693秒。之前使用13.0版本计算该模型所用计算时间为4381秒。新版本提速15.7%。 4）.硬件资源调用情况 由于该模型网格数量较少，仅使用单核进行求解计算。在整个计算过程中，单核占用率达到100%，内存占用峰值约为400 MB。之前使用13.0版本计算该模型的内存占用峰值约

为450兆。新版本对内存的峰值占用约为旧版本的90%左右。 5）.计算精度 VOF模型的计算精度体现在两相交界面捕捉的清晰程度，15.0版本的交界面捕捉清晰程度比旧版本略有提升，对于一些较小的气泡有着更好的捕捉能力。 t=0.45s时，15.0版本和13.0版本计算的两相交界面对比 t=0.45s时，15.0版本和13.0版本计算的两相交界面对比 4、总结 在ANSYS 15.0 Preview3版本的试用过程中，对FLUENT 15.0中VOF模型的计算效率提升感到满意，相比较于旧版本，约有15%的计算速度提升，这对缩短仿真分析的周期有极大帮助；还有约10%的内存峰值占用量下降，这对于合理利用现有硬件资源进行更大规模的模型计算有着重要意义。此外，新版本VOF模型的计算精度也有所提升，两相交界面捕捉更加锐利，对于一些较小的气泡，相对于旧版本有着更好的捕捉能力

FLUENT多孔介质中平面面板(plane surface)工具的使用

1、输出grid图形 2、选择surface---plane,打开plane surface面板 3、通过确定三个点来确定平面位置。单击slect point，出现提示，不点选cancel.在grid 图形的多孔介质区域任意位置右键点选3个点。 4、回到plane surface面板，勾选plane tool，则在grid图形的多孔介质区域出现一个平面。

若出现的平面与我们的预期相差比较大的话，可以单击reset points,可以获得一个特殊位置的平面。 5、打开多孔介质的控制面板，选择porou zone标签，点击update from plane tool按钮，获得方向矢量1，和方向矢量2的原始值，并与左下角的坐标系统比较，确定我们大概的旋转方向。 6、对比grid图形左下角的坐标系统，红线和红色箭头代表的是方向矢量1，绿线和绿色箭头代表的是方向矢量2 应该使红线和X正方向平行，绿线和Y正方向平行。具体的操作应该是： 一：先单击白线的蓝色箭头，固定了该方向在旋转过程中不变，可以保证在旋转的过程比较有规律，然后右键点选白线的红色箭头旋转红线的红色箭头到X的正轴； 二: 接下来应该是单击白线的红色箭头，固定该方向不变，单击白线的蓝色箭头，旋转绿线的绿色箭头指向Y的正轴。（所以多孔介质区域我们一般是设置在坐标系统里面，轴线等 与坐标系统无非直角角度关系）。把平面移动到图形外有利于旋转，比较清楚。平面

法线方向的移动是用鼠标右键单击平面阴影部分并拖动，横向移动则需按下shift并进行如上操作。 7、旋转到适当的位置后（鼠标右键拖动箭头），再次点击update from plane tool按钮，获得方向矢量1，和方向矢量2。 得到的数值很可能不是整数，这个时候我们可以把他简化为整数。例如：0.9123可以简化为1，0.01245可以简化为0，以此类推。

管道阻力的基本计算方法

管道阻力计算 空气在风管内的流动阻力有两种形式：一是由于空气本身的黏滞性以及空气与管壁间的摩擦所产生的阻力称为摩擦阻力；另一是空气流经管道中的管件时(如三通、弯头等)，流速的大小和方向发生变化，由此产生的局部涡流所引起的阻力，称为局部阻力。 一、摩擦阻力 根据流体力学原理，空气在管道内流动时，单位长度管道的摩擦阻力按下式计算： ρ λ 242 v R R s m ?= (5—3) 式中 Rm ——单位长度摩擦阻力，Pa ／m ； υ——风管内空气的平均流速，m ／s ； ρ——空气的密度，kg ／m 3 ； λ——摩擦阻力系数； Rs ——风管的水力半径，m 。 对圆形风管： 4D R s = (5—4) 式中 D ——风管直径，m 。 对矩形风管 )(2b a ab R s += (5—5) 式中 a ，b ——矩形风管的边长，m 。 因此，圆形风管的单位长度摩擦阻力 ρ λ 22 v D R m ?= (5—6) 摩擦阻力系数λ与空气在风管内的流动状态和风管内壁的粗糙度有关。计算摩擦阻力系数的公式很多，美国、日本、德国的一些暖通手册和我国通用通风管道计算表中所采用的公式如下： ) Re 51 .27.3lg( 21 λλ +-=D K (5—7) 式中 K ——风管内壁粗糙度，mm ；

Re ——雷诺数。 υvd = Re (5—8) 式中 υ——风管内空气流速，m ／s ； d ——风管内径，m ； ν——运动黏度，m 2 ／s 。 在实际应用中，为了避免烦琐的计算，可制成各种形式的计算表或线解图。图5—2是计算圆形钢板风管的线解图。它是在气体压力B ＝101．3kPa 、温度t=20℃、管壁粗糙度K ＝0．15mm 等条件下得出的。经核算，按此图查得的Rm 值与《全国通用通风管道计算表》查得的λ／d 值算出的Rm 值基本一致，其误差已可满足工程设计的需要。只要已知风量、管径、流速、单位摩擦阻力4个参数中的任意两个，即可利用该图求得其余两个参数，计算很方便。 图5—2 圆形钢板风管计算线解图 [例] 有一个10m 长薄钢板风管，已知风量L ＝2400m 3 ／h ，流速υ＝16m ／s ，管壁粗糙

基于FLUENT的水下子弹数值模拟

基于FLUENT的水下子弹数值模拟 一、问题介绍 Fluent是目前处于领先地位的CFD软件包之一，它为流体力学领域提供了强大丰富的计算平台。针对各种复杂流动的物理现象，fluent可以采用不同的离散格式和插值方法，在特定区领域内使计算速度、稳定性和精度等方面达到最佳组合。但是针对特殊的物理情景，其自身无法提供完全近似的物理模型，必须通过其自带的自定义函数UDF(user Defined function)对其进行二次开发。例如，对于水下航行的子弹，现有的模块不能模拟子弹运动过程，只能给定一个逆向的速度入口边界条件来实现，即水动子弹不动。为了更近似的模拟子弹运动过程，需要对FLUENT进行二次开发。 本文采取的方法是利用用户自定义函数UDF，对fluent进行二次改造，采用动网格技术，对子弹匀速运动进行模拟。 二、物理平台 本文用到的软件包括：前处理软件gambit,计算软件fluent12.0,UDF支持语言microsoft visual C++6.0,后处理软件tecplot. Gambit用于建立几何模型，划分网格。本算例采用二维轴对称模型，子弹以一矩形近似代替。 Fluent用于计算子弹运动流场特性，包括速度，压力，阻力等。本算例涉及到子弹运动，需要采用动网格技术，即网格在计算过程不断变化更新，因此通过UDF对fluent进行二次开发。其中UDF是fluent自带的宏命令，其运行环境是C语言，所以必须在计算机中提空支持C语言的运算环境microsoft visual C++6.0（也可以是其他更高版本）。 Tecplot为计算结果后处理软件，可以制作动画，云图等。 三、动网格简介 1.spring based smoothing 弹簧近似光滑模型 2.dynamic layering 动态分层模型 3.local remeshing 局部重构模型 三种动网格技术有自己的使用范围，其中弹簧近似光滑模型对于结构和非结构化网格都适用，对于模拟物体小范围运动或者变形有巨大优势，而对于大范围移动却无法实现，一般情况下此方法很少单独使用，通常是配合后面两种动网格技术使用。网格局部重构模型只适合四面体网格和三角形网格，可以用于模拟大范围平移和旋转，同时在模拟复杂边界移动具有很大优势。动态分层模型只是用于结构化网格，对于三角形和四面体网格则毫无用处。由于本文采用此方法，这里讲具体介绍此模型： 动态分层模型（dynamic layering），对于规则计算域，划分规则四边形网格后，采用此模型有利于实现物体的单方向平动，其优点是不需要网格重构，不会改变运动体周围的网格形态，节约了网格更新的时间，同时能保证网格一致性，提高计算效率。此方法的弱点在于无法模拟物体旋转运动，无论对于二维还是三维模型，只能解决物体平动问题。 此方法的思想是当网格边界缩短小于设定值时自动与上一个网格合并成，而当网格被拉伸小于设定的最大网格线长度时，又自动分裂成两个网格。 无论是哪种动网格模型，对内部的参数设置要求十分严格，若参数设置不当，可能出现

管道阻力损失计算

管道的阻力计算 风管内空气流动的阻力有两种，一种是由于空气本身的粘滞性及其与管壁间的摩擦而产生的沿程能量损失，称为摩擦阻力或沿程阻力；另一种是空气流经风管中的管件及设备时，由于流速的大小和方向变化以及产生涡流造成比较集中的能量损失，称为局部阻力。通常直管中以摩擦阻力为主，而弯管以局部阻力阻力为主（图6-1-1）。 图6-1-1 直管与弯管 （一）摩擦阻力 1．圆形管道摩擦阻力的计算 根据流体力学原理，空气在横断面形状不变的管道内流动时的摩擦阻力按下式计 算： （6-1-1） 对于圆形风管，摩擦阻力计算公式可改为： （6-1-2） 圆形风管单位长度的摩擦阻力（又称比摩阻）为：

（6-1-3） 以上各式中 λ——摩擦阻力系数； v——风秘内空气的平均流速，m/s； ρ——空气的密度，kg/m3； l——风管长度，m； Rs——风管的水力半径，m； f——管道中充满流体部分的横断面积，m2； P——湿周，在通风、空调系统中即为风管的周长，m； D——圆形风管直径，m。 摩擦阻力系数λ与空气在风管内的流动状态和风管管壁的粗糙度有关。在通风和空调系统中，薄钢板风管的空气流动状态大多数属于紊流光滑区到粗糙区之间的过渡区。通常，高速风管的流动状态也处于过渡区。只有流速很高、表面粗糙的砖、混凝土风管流动状态才属于粗糙区。计算过渡区摩擦阻力系数的公式很多，下面列出的公式适用范围较大，在目前得到较广泛的采用： （6-1-4） 式中 K——风管内壁粗糙度，mm； D——风管直径，mm。 进行通风管道的设计时，为了避免烦琐的计算，可根据公式（6-1-3）和（6-1-4）制成各种形式的计算表或线解图，供计算管道阻力时使用。只要已知流量、管径、流速、阻力四个参数中的任意两个，即可利用线解图求得其余的两个参数。线解图是按过渡区的λ值，在压力B0=101.3kPa、温度t0=20℃、宽气密度ρ0=1.204kg/m3、运动粘度v0=15.06×10－6m2/s、管壁粗糙度K=0.15mm、圆形风管等条件下得出的。当实际使用条件下上述条件不相符时，应进行修正。 （1）密度和粘度的修正

基于FLUENT的某微型面包车外流场数值仿真分析

基于FLUENT的某微型面包车外流场数值仿真分析 摘要：利用UG对某微型轿车进行三维实体建模，将其导入前处理软件ANSYS ICEM中，建立计算域后网格化。用CFD 软件FLUENT对汽车模型的外流场进行三维稳态流动数值模拟，得出汽车周围流场的气流速度和压力分布，并通过计算得到了该车的阻力系数，该仿真分析的数据为进行汽车气动特性分析提供基础，可进一步指导汽车的设计开发。 关键字：汽车空气动力学；计算流体动力学；FLUENT；外流场 ABSTRACT：The 3D model of a mini car is carried out by UG, and then it is introduced into the pre-processing software ANSYS ICEM to establish the computational domain grid. By using the software of CFD FLUENT is to the automobile model flow field numerical simulation of three-dimensional steady flow, flow velocity and pressure distribution of the flow field around the car, the car and the drag coefficient is obtained by calculating, the simulation data for automobile gas analysis provide the basic dynamic characteristics, design and development can further guide the car. Keywords: Automobile aerodynamics；CFD；FLUENT；Outflow field 引言 空气动力学特性是汽车的重要特性之一。汽车行驶时与空气产生复杂的相互作用，承受着强大的气动力，对汽车的行驶状态有着重大影响；汽车行驶时受到的空气阻力与汽车速度平方成正比，汽车克服空气阻力所消耗的功率和燃料与车速的三次方成正比。因此，对汽车外流场空气动力学的研究，不仅可以提高汽车动力性和安全性，还可以提高汽车的燃料经济性。 目前，汽车空气动力学的研究主要有三种方法，即风洞实验、理论分析和计算流体动力学（Computational Fluid Dynamics，CFD）分析。随着计算机技术的发展，计算流体动力学相对于实验和理论计算具有成本低、周期短等特点，因此受到越来越广泛的应用。CFD方法对于预测和改进汽车的气动性能，指导汽车产品设计具有重要意义[1-2]。故本文采用大型商业化CFD 软件的FLUENT对某微型汽车的外流场进行数值仿真分析。 1 汽车空气动力学特性与CFD 理论基础1.1 汽车空气动力学特性 在正常道路行驶过程中的汽车，通常受到两种力的作用，这两种力分别为路面与汽车之间的相互作用力和来自空气的力与力矩。其中，前者主要由汽车自身的物理属性和轮胎的滚动阻力系数等决定；另一种则是来自空气作用的力和力矩，取决于汽车的外形设计、行驶速度以及横摆角[3-4]。 气流作用于汽车上分相互垂直的三个方向的力和绕三个轴的力矩，如图1 所示；在图示坐标系中，X 方向是汽车直线行驶方向，通常的气动阻力就是指来流沿X 方向的作用力；Y 轴方向为汽车的侧向力，还有沿Z 轴方向趋于抬起汽车的升力。汽车在道路行驶过程中这三个方向的空气作用力同时存在，相互影响。除上述三个方向的力外，还有绕三个轴方向的力矩，分别为绕X 轴方向的侧倾力矩、绕Y 轴方向的俯仰力矩以及绕Z 轴方向的横摆力矩。

fluent经验总结

1什么叫松弛因子？松弛因子对计算结果有什么样的影响？它对计算的收敛情况又有什 么样的影响？ 1、亚松驰（Under Relaxation）：所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减，以避免由于差值过大而引起非线性迭代过程的发散。用通用变量来写 出时，为松驰因子（Relaxation Factors）。《数值传热学-214》 2、FLUENT中的亚松驰：由于FLUENT所解方程组的非线性，我们有必要控制的变化。一般用亚松驰方法来实现控制，该方法在每一部迭代中减少了的变化量。亚松驰最简 单的形式为：单元内变量等于原来的值加上亚松驰因子a与变化的积, 分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程，包 括耦合解算器所解的非耦合方程（湍流和其他标量）都会有一个相关的亚松驰因子。在FLUENT中，所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题，但是对于一些特殊的非线性问题（如：某些湍流或者高Rayleigh数自然对流问题），在计算开始时要慎重减小亚松驰因子。使用默认的亚松驰因子开始计算是很好的习惯。如 果经过4到5步的迭代残差仍然增长，你就需要减小亚松驰因子。有时候，如果发现残差 开始增加，你可以改变亚松驰因子重新计算。在亚松驰因子过大时通常会出现这种情况。 最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文件，并对解的算法做几 步迭代以调节到新的参数。最典型的情况是，亚松驰因子的增加会使残差有少量的增加， 但是随着解的进行残差的增加又消失了。如果残差变化有几个量级你就需要考虑停止计算 并回到最后保存的较好的数据文件。注意：粘性和密度的亚松驰是在每一次迭代之间的。 而且，如果直接解焓方程而不是温度方程（即：对PDF计算），基于焓的温度的更新是要进行亚松驰的。要查看默认的亚松弛因子的值，你可以在解控制面板点击默认按钮。对于 大多数流动，不需要修改默认亚松弛因子。但是，如果出现不稳定或者发散你就需要减小 默认的亚松弛因子了，其中压力、动量、k和e的亚松弛因子默认值分别为0.2，0.5，0.5和0.5。对于SIMPLEC格式一般不需要减小压力的亚松弛因子。在密度和温度强烈耦合 的问题中，如相当高的Rayleigh数的自然或混合对流流动，应该对温度和/或密度（所用 的亚松弛因子小于1.0）进行亚松弛。相反，当温度和动量方程没有耦合或者耦合较弱时，流动密度是常数，温度的亚松弛因子可以设为1.0。对于其它的标量方程，如漩涡，组分，PDF变量，对于某些问题默认的亚松弛可能过大，尤其是对于初始计算。你可以将松弛因子设为0.8以使得收敛更容易。 SIMPLE与SIMPLEC比较 在FLUENT中，可以使用标准SIMPLE算法和SIMPLEC（SIMPLE-Consistent）算法，默认是SIMPLE算法，但是对于许多问题如果使用SIMPLEC可能会得到更好的结果，尤其是可以应用增加的亚松驰迭代时，具体介绍如下: 对于相对简单的问题（如：没有附加模型激活的层流流动），其收敛性已经被压力速

矿井通风阻力计算方法

矿井通风阻力 第一节通风阻力产生的原因 当空气沿井巷运动时，由于风流的粘滞性和惯性以及井巷壁面等对风流的阻滞、扰动作用而形成通风阻力，它是造成风流能量损失的原因。 井巷通风阻力可分为两类：摩擦阻力(也称为沿程阻力)和局部阻力。 一、风流流态（以管道流为例） 同一流体在同一管道中流动时，不同的流速，会形成不同的流动状态。当流速较低时，流体质点互不混杂，沿着与管轴平行的方向作层状运动，称为层流(或滞流)。当流速较大时，流体质点的运动速度在大小和方向上都随时发生变化，成为互相混杂的紊乱流动，称为紊流(或湍流)。(降低风速的原因) (二)、巷道风速分布 由于空气的粘性和井巷壁面摩擦影响，井巷断面上风速分布是不均匀的。 在同一巷道断面上存在层流区和紊区，在贴近壁面处仍存在层流运动薄层，即层流区。在层流区以外，为紊流区。从巷壁向巷道轴心方向，风速逐渐增大，呈抛物线分布。 巷壁愈光滑，断面上风速分布愈均匀。 第二节摩擦阻力与局部阻力的计算 一、摩擦阻力 风流在井巷中作沿程流动时，由于流体层间的摩擦和流体与井巷壁面之间的摩擦所形成的阻力称为摩擦阻力(也叫沿程阻力)。 由流体力学可知，无论层流还是紊流，以风流压能损失（能量损失）来反映的摩擦阻力可用下式来计算： H f =λ×L/d×ρν2/2pa λ——摩擦阻力系数。 L——风道长度，m

d——圆形风管直径，非圆形管用当量直径； ρ——空气密度，kg/m3 ν2——断面平均风速，m/s； 1、层流摩擦阻力：层流摩擦阻力与巷道中的平均流速的一次方成正比。因井下多为紊流，故不详细叙述。 2、紊流摩擦阻力：对于紊流运动，井巷的摩擦阻力计算式为： H f =α×LU/S3×Q2 =R f×Q2pa R f=α×LU/S3 α——摩擦阻力系数，单位kgf·s2/m4或N·s2/m4，kgf·s2/m4=9.8N·s2/m4 L、U——巷道长度、周长，单位m； S——巷道断面积，m2 Q——风量，单位m/s R f——摩擦风阻，对于已给定的井巷，L，U，S都为已知数，故可把上式中的α，L，U，S 归结为一个参数R f，其单位为：kg/m7 或N·s2/m8 3、井巷摩擦阻力计算方法 新建矿井：查表得α→h f→R f 生产矿井：已测定的h f→R f→α，再由α→h f→R f 二、局部阻力 由于井巷断面，方向变化以及分岔或汇合等原因，使均匀流动在局部地区受到影响而破坏，从而引起风流速度场分布变化和产生涡流等，造成风流的能量损失，这种阻力称为局部阻力。由于局部阻力所产生风流速度场分布的变化比较复杂性，对局部阻力的计算一般采用经验公式。 1、几种常见的局部阻力产生的类型： （1）、突变 紊流通过突变部分时，由于惯性作用，出现主流与边壁脱离的现象，在主流与边壁之间形成涡漩区，从而增加能量损失。

矿井通风阻力计算方法

矿井通风阻力 第一节通风阻力产生的原因当空气沿井巷运动时，由于风流的粘滞性和惯性以及井巷壁面等对风流的阻滞、扰动作用而形成通风阻力，它是造成风流能量损失的原因。 井巷通风阻力可分为两类：摩擦阻力（也称为沿程阻力）和局部阻力。 一、风流流态（以管道流为例）同一流体在同一管道中流动时，不同的流速，会形成不同的流动状态。当流速较低时，流体质点互不混杂，沿着与管轴平行的方向作层状运动，称为层流（或滞流）。当流速较大时，流体质点的运动速度在大小和方向上都随时发生变化，成为互相混杂的紊乱流动，称为紊流（或湍流）。（降低风速的原因） （二）、巷道风速分布 由于空气的粘性和井巷壁面摩擦影响，井巷断面上风速分布是不均匀的。在同一巷道断面上存在层流区和紊区，在贴近壁面处仍存在层流运动薄层，即层流区。在层流区以外，为紊流区。从巷壁向巷道轴心方向，风速逐渐增大，呈抛物线分布。 巷壁愈光滑，断面上风速分布愈均匀。 第二节摩擦阻力与局部阻力的计算 一、摩擦阻力风流在井巷中作沿程流动时，由于流体层间的摩擦和流体与井巷壁面之间的摩擦所形成的阻力称为摩擦阻力（也叫沿程阻力）。 由流体力学可知，无论层流还是紊流，以风流压能损失（能量损失）来反映的摩擦阻力可用下式来计算： 2 H = λ×L/d ×ρν/2 Pa λ——摩擦阻力系数。 L ---- 风道长度，m d――圆形风管直径，非圆形管用当量直径；

空气密度，kg/m3 断面平均风速，m/s； 1、层流摩擦阻力：层流摩擦阻力与巷道中的平均流速的一次方成正比。因井下多为紊流，故不详细叙述。 2、紊流摩擦阻力：对于紊流运动，井巷的摩擦阻力计算式为： H = α ×LU∕S3×Q2 =R f ×Q2 Pa 3 R f=α× LU∕S3 α --- 摩擦阻力系数，单位kgf ?s2∕m4或N ? s7m4, kgf ?s7m4=9.8N ? s7m4 L、U――巷道长度、周长，单位m S—巷道断面积，m Q ---- 风量，单位m/s R ——摩擦风阻，对于已给定的井巷，L，U S都为已知数，故可把上式中的α, L, U, S归结为一个参数R,其单位为：kg∕m7或N ?s7m8 3、井巷摩擦阻力计算方法 新建矿井：查表得α→ h f → R f 生产矿井：已测定的h f → R f → α, 再由α→ h f → R f 二、局部阻力 由于井巷断面,方向变化以及分岔或汇合等原因, 使均匀流动在局部地区受到影响而破坏, 从而引起风流速度场分布变化和产生涡流等,造成风流的能量损失,这种阻力称为局部阻力。由于局部阻力所产生风流速度场分布的变化比较复杂性,对局部阻力的计算一般采用经验公式。 1、几种常见的局部阻力产生的类型： （1）、突变紊流通过突变部分时,由于惯性作用,出现主流与边壁脱离的现象,在主流与边壁之间形成涡漩区,从而增加能量损失。 （2）、渐变 主要是由于沿流动方向出现减速增压现象, 在边壁附近产生涡漩。因为压差