文档库 最新最全的文档下载
当前位置:文档库 › 两个Fluent实验

两个Fluent实验

两个Fluent实验
两个Fluent实验

一、实验目的

1、计算圆柱绕流问题,比较不同雷诺数下的异同,并截图查看相关涡街。

2、模拟不同雷诺数下的方腔环流,并比较不同的网格设置对方腔环流的影响。

二、实验原理

1、Re=20的圆柱绕流问题。

Setup and solution of a classical problems of fluid mechanics such as a steady flow past circular cylinder. (Re =20) We know D = 2 m. To obtain Re = 20, we can arbitrarily set ρ, V and μ. For our case, let's set ρ = 1 kg/m 3 , V = 1 m/s and μ = 0.1 k g/ms.

我们取大圆半径为10,小圆半径为2,矩形区域为10*40。在Gambit软件中建模,并导出为msh文件。将msh 导入fluent软件中进行求解,相应的参数如上所述,作出相应结果的图。

2、Re=150的圆柱绕流问题。

setup and solution of a classical problems of fluid mechanics such as a steady flow past circular cylinder. (Re =150) We know D = 2 m. To obtain Re = 150, we can arbitrarily set ρ, V and μ. For our case, let's set ρ = 75 kg/m 3 , V = 1 m/s and μ = 1 kg/ms.

与上一个问题的建模相同,在fluent软件中需要调整部分参数。

3、方腔环流问题。

使用fluent计算之前的方腔环流问题,并与Fortran程序计算结果进行对比。

(1)均匀网格计算Re=100,1000,5000的情况

(2)加密网格(边界加密)计算Re=100,1000,5000的情况

4、注:本实验部分采用Fluent 6.3,部分采用Fluent 15 完成。

三、实验结果

1、Re=20的圆柱绕流问题

我们在Gambit中按照如下标准划分网格:

划分的结果如下图:

定义边界条件:A: Velocity Inlet;BC & DE: Periodic;FGH: Pressure-Outlet (or Outflow);IJ: Wall 输出msh文件,用fluent进行求解,参数定义如下:

◆Solver: default setting (steady)

◆Viscous: default setting (laminar)

◆Energy: default setting (energy equation not

activated)

◆Material: set the density and viscosity for air ◆Operation Condition: reference pressure = 1 atm ◆Boundary Conditions: velocity = 1m/s

◆Discretization for convective terms: Second Order

Upwind

Initilization

◆Monitor Residuals & Force

◆Save the case and data files ◆Setting Reference Values for pressure coefficient

然后我们将求解的结果图放在下面:

Plot Velocity Vectors

Display Pressure Coefficient on the cylinder

Plot Stream Function

Plot Vorticity Magnitude

2、Re=150的圆柱绕流问题

与Re=20的问题类似,网格划分不变,将fluent中的一些参数做一下修改,便可以进行重新求解。

定义好各种参数后,我们可以求解,并像上面一样作出Contours of Cp以及Contours of Velocity Magnitude。

我们对比Re=20和Re=150的结果,可以看到,后者流场已经出现涡街的雏形,事实上当雷诺数更大如Re=2000时,涡街就已经十分明显。

3、方腔环流

我们取10*10边长的正方形网格,通过改变上边界板的速度从而改变雷诺数。在Gambit中进行网格划分,方腔网格分别采用均匀网格和double side 的设置,从而加密边界网格,如下图所示。

均匀网格的情况,Re=100,作出流函数和速度矢量图,这里精度取为10?3。

Re=1000情形

接下来我们加密网格,考察Re=1000和Re=5000的情形,精度为10?6。Re=1000

Re=5000

比较可知,网格加密,精度提高,会使结果更加精确,Re=1000下的二次涡变得明显,Re=5000在高精度下二次涡十分明显。由此,也可以佐证Fortran程序及Matlab作图的结果。

四、实验总结

1、圆柱绕流实验需要注意网格的划分,要在圆柱周围加密网格,圆柱的上下部分也要加密网格,并且要在尾流部分按照由密到疏的程度划分网格,这样才能看到比较明显的涡街雏形。

2、Fluent 实验中,精度的选择有时会影响图像的观察。精度越高,需要的计算步更多,花费的时间更长,但同时计算的结果也更加精确,图像也更明显。

FLUENT中文全教程1-250

FLUENT 教程 赵玉新 I、目录 第一章、开始 第二章、操作界面 第三章、文件的读写 第四章、单位系统 第五章、读入和操作网格 第六章、边界条件 第七章、物理特性 第八章、基本物理模型 第九章、湍流模型 第十章、辐射模型 第十一章、化学输运与反应流 第十二章、污染形成模型 第十三章、相变模拟 第十四章、多相流模型 第十五章、动坐标系下的流动 第十六章、解算器的使用 第十七章、网格适应 第十八章、数据显示与报告界面的产生 第十九章、图形与可视化 第二十章、Alphanumeric Reporting 第二十一章、流场函数定义 第二十二章、并行处理 第二十三章、自定义函数 第二十四章、参考向导 第二十五章、索引(Bibliography) 第二十六章、命令索引 II、如何使用该教程 概述 本教程主要介绍了FLUENT 的使用,其中附带了相关的算例,从而能够使每一位使用 者在学习的同时积累相关的经验。本教程大致分以下四个部分:第一部分包括介绍信息、用户界面信息、文件输入输出、单位系统、网格、边界条件以及物理特性。第二和第三部分包含物理模型,解以及网格适应的信息。第四部分包括界面的生成、后处理、图形报告、并行处理、自定义函数以及FLUENT 所使用的流场函数与变量的定义。 下面是各章的简略概括 第一部分: z开始使用:本章描述了FLUENT 的计算能力以及它与其它程序的接口。介绍了如何对具体的应用选择适当的解形式,并且概述了问题解决的大致步骤。在本章中,我们给出

了一个可以在你自己计算机上运行的简单的算例。 z使用界面:本章描述了用户界面、文本界面以及在线帮助的使用方法。同时也提供了远程处理与批处理的一些方法。(请参考关于特定的文本界面命令的在线帮助) z读写文件:本章描述了FLUENT 可以读写的文件以及硬拷贝文件。 z单位系统:本章描述了如何使用FLUENT 所提供的标准与自定义单位系统。 z读和操纵网格:本章描述了各种各样的计算网格来源,并解释了如何获取关于网格的诊断信息,以及通过尺度化(scale)、分区(partition)等方法对网格的修改。本章还描述了非一致(nonconformal)网格的使用. z边界条件:本章描述了FLUENT 所提供的各种类型边界条件,如何使用它们,如何定义它们and how to define boundary profiles and volumetric sources. z物理特性:本章描述了如何定义流体的物理特性与方程。FLUENT 采用这些信息来处理你的输入信息。 第二部分: z基本物理模型:本章描述了FLUENT 计算流体流动和热传导所使用的物理模型(包括自然对流、周期流、热传导、swirling、旋转流、可压流、无粘流以及时间相关流)。以及在使用这些模型时你需要输入的数据,本章也包含了自定义标量的信息。 z湍流模型:本章描述了FLUENT 的湍流模型以及使用条件。 z辐射模型:本章描述了FLUENT 的热辐射模型以及使用条件。 z化学组分输运和反应流:本章描述了化学组分输运和反应流的模型及其使用方法。本章详细的叙述了prePDF 的使用方法。 z污染形成模型:本章描述了NOx 和烟尘的形成的模型,以及这些模型的使用方法。 第三部分: z相变模拟:本章描述了FLUENT 的相变模型及其使用方法。 z离散相变模型:本章描述了FLUENT 的离散相变模型及其使用方法。 z多相流模型:本章描述了FLUENT 的多相流模型及其使用方法。 z Flows in Moving Zones(移动坐标系下的流动):本章描述了FLUENT 中单一旋转坐标系,多重移动坐标系,以及滑动网格的使用方法。 z Solver 的使用:本章描述了如何使用FLUENT 的解法器(solver)。 z网格适应:本章描述了explains the solution-adaptive mesh refinement feature in FLUENT and how to use it 第四部分: z显示和报告数据界面的创建:本章描述了explains how to create surfaces in the domain on which you can examine FLUENT solution data z图形和可视化:本章描述了检验FLUENT 解的图形工具 z Alphanumeric Reporting:本章描述了如何获取流动、力、表面积分以及其它解的数据。 z流场函数的定义:本章描述了如何定义FLUENT 面板内出现的变量选择下拉菜单中的流动变量,并且告诉我们如何创建自己的自定义流场函数。 z并行处理:本章描述了FLUENT 的并行处理特点以及使用方法 z自定义函数:本章描述了如何通过用户定义边界条件,物理性质函数来形成自己的FLUENT 软件。 如何使用该手册 z根据你对CFD 以及FLUENT 公司的熟悉,你可以通过各种途径使用该手册 对于初学者,建议如下:

《FLUENT中文手册(简化版)》

FLUENT中文手册(简化版) 本手册介绍FLUENT的使用方法,并附带了相关的算例。下面是本教程各部分各章节的简略概括。 第一部分: ?开始使用:描述了FLUENT的计算能力以及它与其它程序的接口。介绍了如何对具体的应用选择适当的解形式,并且概述了问题解决的大致步骤。在本章中给出了一个简单的算例。 ?使用界面:描述用户界面、文本界面以及在线帮助的使用方法,还有远程处理与批处理的一些方法。?读写文件:描述了FLUENT可以读写的文件以及硬拷贝文件。 ?单位系统:描述了如何使用FLUENT所提供的标准与自定义单位系统。 ?使用网格:描述了各种计算网格来源,并解释了如何获取关于网格的诊断信息,以及通过尺度化(scale)、分区(partition)等方法对网格的修改。还描述了非一致(nonconformal)网格的使用. ?边界条件:描述了FLUENT所提供的各种类型边界条件和源项,如何使用它们,如何定义它们等 ?物理特性:描述了如何定义流体的物理特性与方程。FLUENT采用这些信息来处理你的输入信息。 第二部分: ?基本物理模型:描述了计算流动和传热所用的物理模型(包括自然对流、周期流、热传导、swirling、旋转流、可压流、无粘流以及时间相关流)及其使用方法,还有自定义标量的信息。 ?湍流模型:描述了FLUENT的湍流模型以及使用条件。 ?辐射模型:描述了FLUENT的热辐射模型以及使用条件。 ?化学组分输运和反应流:描述了化学组分输运和反应流的模型及其使用方法,并详细叙述了prePDF 的使用方法。 ?污染形成模型:描述了NOx和烟尘的形成的模型,以及这些模型的使用方法。 第三部分: ?相变模拟:描述了FLUENT的相变模型及其使用方法。 ?离散相变模型:描述了FLUENT的离散相变模型及其使用方法。 ?多相流模型:描述了FLUENT的多相流模型及其使用方法。 ?移动坐标系下的流动:描述单一旋转坐标系、多重移动坐标系、以及滑动网格的使用方法。 ?解法器(solver)的使用:描述了如何使用FLUENT的解法器。 ?网格适应:描述了如何优化网格以适应计算需求。 第四部分: ?显示和报告数据界面的创建:本章描述了explains how to create surfaces in the domain on which you can examine FLUENT solution data ?图形和可视化:本章描述了检验FLUENT解的图形工具 ?Alphanumeric Reporting:本章描述了如何获取流动、力、表面积分以及其它解的数据。 ?流场函数的定义:本章描述了如何定义FLUENT面板内出现的变量选择下拉菜单中的流动变量,并且告诉我们如何创建自己的自定义流场函数。 ?并行处理:本章描述了FLUENT的并行处理特点以及使用方法 ?自定义函数:本章描述了如何通过用户定义边界条件,物理性质函数来形成自己的FLUENT软件。 如何使用该手册 对于初学者,建议从阅读“开始”这一章起步。 对于有经验的使用者,有三种不同的方法供你使用该手册:按照特定程序的步骤从按程序顺序排列的目录列表和主题列表中查找相关资料;从命令索引查找特定的面板和文本命令的使用方法;从分类索引查找特定类别信息(在线帮助中没有此类索引,只能在印刷手册中找到它)。 什么时候使用Support Engineer:Support Engineer能帮你计划CFD模拟工程并解决在使用FLUENT 中所遇到的困难。在遇到困难时我们建议你使用Support Engineer。但是在使用之前有以下几个注意事项:●仔细阅读手册中关于你使用并产生问题的命令的信息 ●回忆导致你产生问题的每一步 ●如果可能的话,请记下所出现的错误信息 ●对于特别困难的问题,保存FLUENT出现问题时的日志以及手稿。在解决问题时,它是最好的资源。

Fluent经典问题及解答

Fluent经典问题及解答 1 对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢?(#61) 2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。(13楼) 3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不同?(#80) 4 常见离散格式的性能的对比(稳定性、精度和经济性)(#62) 5 在利用有限体积法建立离散方程时,必须遵守哪几个基本原则?(#81) 6 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么?(#130) 7 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难?(#55) 8 什么叫边界条件?有何物理意义?它与初始条件有什么关系?(#56) 9 在一个物理问题的多个边界上,如何协调各边界上的不同边界条件?在边界条件的组合问题上,有什么原则? 10 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别?(#143) 11 在网格生成技术中,什么叫贴体坐标系?什么叫网格独立解?(#35) 12 在GAMBIT的foreground和background中,真实体和虚实体、实操作和虚操作四个之间是什么关系? 13 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量?及其在做网格时大致注意到哪些细节?(#38) 14 画网格时,网格类型和网格方法如何配合使用?各种方法有什么样的应用范围及做网格时需注意的问题?(#169) 15 对于自己的模型,大多数人有这样的想法:我的模型如何来画网格?用什么样的方法最简单?这样做网格到底对不对?(#154) 16 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?(#40) 17 依据实体在GAMBIT建模之前简化时,必须遵循哪几个原则?(#170) 18 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理?b、计算域内的内部边界如何处理(2D)?(#128) 19 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些?(#127) 20 何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT是怎样使用区域的?(#41) 21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?(9楼) 22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?(7楼)

Fluent多相流模型选择

FLUENT多相流模型 分类 1、气液或液液流动 气泡流动:连续流体中存在离散的气泡或液泡 液滴流动:连续相为气相,其它相为液滴 栓塞(泡状)流动:在连续流体中存在尺寸较大的气泡 分层自由流动:由明显的分界面隔开的非混合流体流动。 2、气固两相流动 粒子负载流动:连续气体流动中有离散的固体粒子 气力输运:流动模式依赖,如固体载荷、雷诺数和例子属性等。最典型的模式有沙子的流动,泥浆流,填充床以及各相同性流 流化床:有一个盛有粒子的竖直圆筒构成,气体从一个分散器进入筒内,从床底不断冲入的气体使得颗粒得以悬浮。 3、液固两相流动 泥浆流:流体中的大量颗粒流动。颗粒的stokes数通常小于1。大于1是成为流化了的液固流动。 水力运输:在连续流体中密布着固体颗粒 沉降运动:在有一定高度的盛有液体的容器内,初始时刻均匀散布着颗粒物质,随后,流体会出现分层。 4、三相流 以上各种情况的组合 多相流动系统的实例 气泡流:抽吸、通风、空气泵、气穴、蒸发、浮选、洗刷。 液滴流:抽吸、喷雾、燃烧室、低温泵、干燥机、蒸发、气冷、洗刷。 栓塞流:管道或容器中有大尺度气泡的流动 分层流:分离器中的晃动、核反应装置沸腾和冷凝 粒子负载流:旋风分离器、空气分类器、洗尘器、环境尘埃流动 气力输运:水泥、谷粒和金属粉末的输运 流化床:流化床反应器、循环流化床 泥浆流:泥浆输运、矿物处理 水力输运:矿物处理、生物医学、物理化学中的流体系统 沉降流动:矿物处理。 多相流模型的选择原则 1、基本原则

1)对于体积分数小于10%的气泡、液滴和粒子负载流动,采用离散相 模型。 2)对于离散相混合物或者单独的离散相体积率超出10%的气泡、液滴 和粒子负载流动,采用混合模型或欧拉模型。 3)对于栓塞流、泡状流,采用VOF模型 4)对于分层/自由面流动,采用VOF模型 5)对于气动输运,均匀流动采用混合模型,粒子流采用欧拉模型。 6)对于流化床,采用欧拉模型 7)泥浆和水力输运,采用混合模型或欧拉模型。 8)沉降采用欧拉模型 9)对于更一般的,同时包含多种多相流模式的情况,应根据最感兴趣 的流动特种,选择合适的流动模型。此时由于模型只是对部分流动特 征采用了较好的模拟,其精度必然低于只包含单个模式的流动。 2、混合模型和欧拉模型的选择原则 VOF模型适合于分层的或自由表面流,而混合模型和欧拉模型适合于流动中有相混合或分离,或者分散相的体积分数超过10%的情况(小于10%可使用离散相模型)。 1)如果分散相有宽广的分布(如颗粒的尺寸分布很宽),最好采用混 合模型,反之使用欧拉模型。 2)如果相间曳力规律一直,欧拉模型通常比混合模型更精确;若相间 曳力规律不明确,最好选用混合模型。 3)如果希望减小计算了,最好选用混合模型,它比欧拉模型少解一部 分方程;如果要求精度而不在意计算量,欧拉模型可能是更好的选择。 但是要注意,复杂的欧拉模型比混合模型的稳定性差,可能会遇到收 敛困难。

fluent经典问题整理

网格质量与那些因素有关? 网格质量本身与具体问题的具体几何特性、流动特性及流场求解算法有关。因此,网格质量最终要由计算结果来评判,但是误差分析以及经验表明,CFD计算对计算网格有一些一般性的要求,例如光滑性、正交性、网格单元的正则性以及在流动变化剧烈的区域分布足够多的网格点等。对于复杂几何外形的网格生成,这些要求往往并不可能同时完全满足。例如,给定边界网格点分布,采用Laplace 方程生成的网格是最光滑的,但是最光滑的网格不一定满足物面边界正交性条件,其网格点分布也很有可能不能捕捉流动特征,因此,最光滑的网格不一定是最好的网格。对计算网格的一个最基本的要求当然是所有网格点的Jacobian必须为正值,即网格体积必须为正,其他一些最常用的网格质量度量参数包括扭角(skew angle)、纵横比(aspect ratio、Laplacian)、以及弧长(arc length)等。通过计算、检查这些参数,可以定性的甚至从某种程度上定量的对网格质量进行评判。Parmley等给出了更多的基于网格元素和网格节点的网格质量度量参数。有限元素法关于插值逼近误差估计的理论,实际上也对网格单元的品质给出了基本的规定:即每个单元的内切球半径与外切球半径之,应该是一个适当的,与网格疏密无关的常数。 实体与虚体的区别 在建模中,经常会遇到实体、实面与虚体、虚面,虚体的计算域也可以进行计算并得到所需的结果。那么它们的区别是什么呢? 对于求解是没有任何区别的,只要你能在虚体或者实体上划分你需要的网格。关键是看你网格生成的质量如何,与实体虚体无关。 gambit的实体和虚体在生成网格和计算的时候对于结果没有任何影响,实体和虚体的主要区别有以下几点: 1.实体可以进行布尔运算但是虚体不能,虽然不能进行布尔运算,但是虚体存在merge,split 等功能。 2.实体运算在很多cad软件里面都有,但是虚体是gambit的一大特色,有了虚体以后,gambit 的建模和网格生成的灵活性增加了很多。 3.在网格生成的过程中,如果有几个相对比较平坦的面,你可以把它们通过merge合成一个,这样,作网格的时候,可以节省步骤,对于曲率比较大的面,可能生成的网格质量不好,这时候,你可以采取用split的方式把它划分成几个小面以提高网格质量。 在Fluent中进行非稳态(unsteady)计算时如何设置步长?

Fluent多相流模型选择与设定

1.多相流动模式 我们可以根据下面的原则对多相流分成四类: ?气-液或者液-液两相流: o 气泡流动:连续流体中的气泡或者液泡。 o 液滴流动:连续气体中的离散流体液滴。 o 活塞流动: 在连续流体中的大的气泡 o 分层自由面流动:由明显的分界面隔开的非混合流体流动。 ?气-固两相流: o 充满粒子的流动:连续气体流动中有离散的固体粒子。 o 气动输运:流动模式依赖诸如固体载荷、雷诺数和粒子属性等因素。最典型的模式有沙子的流动,泥浆流,填充床,以及各向同性流。 o 流化床:由一个盛有粒子的竖直圆筒构成,气体从一个分散器导入筒内。从 床底不断充入的气体使得颗粒得以悬浮。改变气体的流量,就会有气泡不断 的出现并穿过整个容器,从而使得颗粒在床内得到充分混合。 ?液-固两相流 o 泥浆流:流体中的颗粒输运。液-固两相流的基本特征不同于液体中固体颗 粒的流动。在泥浆流中,Stokes 数通常小于1。当Stokes数大于1 时,流动成为流化(fluidization)了的液-固流动。 o 水力运输: 在连续流体中密布着固体颗粒 o 沉降运动: 在有一定高度的成有液体的容器内,初始时刻均匀散布着颗粒物 质。随后,流体将会分层,在容器底部因为颗粒的不断沉降并堆积形成了淤 积层,在顶部出现了澄清层,里面没有颗粒物质,在中间则是沉降层,那里 的粒子仍然在沉降。在澄清层和沉降层中间,是一个清晰可辨的交界面。 ?三相流(上面各种情况的组合) 各流动模式对应的例子如下: ?气泡流例子:抽吸,通风,空气泵,气穴,蒸发,浮选,洗刷 ?液滴流例子:抽吸,喷雾,燃烧室,低温泵,干燥机,蒸发,气冷,刷洗?活塞流例子:管道或容器内有大尺度气泡的流动 ?分层自由面流动例子:分离器中的晃动,核反应装置中的沸腾和冷凝 ?粒子负载流动例子:旋风分离器,空气分类器,洗尘器,环境尘埃流动 ?风力输运例子:水泥、谷粒和金属粉末的输运

Fluent动网格----layering个一个简单实例(作者Snow)

Fluent动网格----layering个一个简单实例我这几天看了点动网格技术方面的东西,在学习过程中发现这方面的例子很少,自己也走了一些弯路。现在还好,弄明白了一些,能够应付现在我的工作。为了让更多学习者快速了解动网格,我打算尽量把我学习心得在这里和大家分享,这里给出一个layering的一个简单例子。 1.Gambit画网格 本例很简单,在Gambit里画一个10*10的矩形,网格间隔为1,也就是有100个网格,具体见下图。都学动网格的人了,不至于这个不会做! 这里需要注意一个问题:设置边界条件的时候,一定要把要移动的边单独设定,本例中一右边界作为移动的边,设成wall就可以,这里再后面需要制定。 2.编写UDF #include "udf.h" #include "unsteady.h" #include "stdio.h" #include "stdlib.h" /************************************************************/ real current_time = 0.0 ; Domain * domain ; Thread * thread ; real NV_VEC( origin ),NV_VEC( force ),NV_VEC( moment ) ; /************************************************************/ DEFINE_CG_MOTION(throttle,dt,vel,omega,time,dtime) { current_time = CURRENT_TIME ; vel[0] = 30; Message("time=%f omega=%f\n",current_time) ; }

fluent多相流算例

Tutorial:Dam-Break Simulation Using FLUENT’s Volume of Fluid Model Purpose This tutorial examines the dam-break problem using the Volume of Fluid(VOF)multiphase model. This tutorial demonstrates how to do the following: ?Set up a dam-break problem. ?Choose the time step by estimating the maximum possible velocity of the interface and the grid cell dimension. ?Solve the problem using the VOF model. ?Manipulate the solution parameters. Prerequisites This tutorial assumes that you are familiar with the FLUENT interface and that you have a good understanding of basic setup and solution procedures.In this tutorial,you will use VOF multiphase model,so you should have some experience with it.This tutorial will not cover the mechanics of using this model;instead,it will focus on the application of this model to solve a dam-break problem. If you have not used this model before,it would be helpful to?rst refer to the FLUENT6.3 User’s Guide and the FLUENT6.2Tutorial Guide. Problem Description The initial setup of the dam-break problem is shown in Figure1. In this problem,a rectangular column of water,in hydrostatic equilibrium,is con?ned between two walls.Gravity is acting downwards with a magnitude of-9.81m/s2.At the beginning of the calculation,the right wall is removed and the water is allowed to?ow out to the horizontal wall.

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的优势在于收敛速度快,但是低速流场计算,其大多数的着眼点在于对流场结构的探索,所以计算得到的

FLUENT中两相流、多相流中模型的的选择问题

两相流:通常把含有大量固体或液体颗粒的气体或液体流动称为两相流;其中含有多种尺寸组颗粒群为一个“相”,气体或液体为另一“相”,由此就有气—液,气—固,液—固等两相流之分。 两相流的研究:对两相流的研究有两种不同的观点:一是把流体作为连续介质,而把颗粒群作为离散体系;而另一是除了把流体作为连续介质外,还把颗粒群当作拟连续介质或拟流体。 引入两种坐标系:即拉格朗日坐标和欧拉坐标,以变形前的初始坐标为自变量称为拉格朗日Langrangian 坐标或物质坐标;以变形后瞬时坐标为自变量称为欧拉Eulerian 坐标或空间坐标。 一.离散相模型 FLUENT在求解连续相的输运方程的同时,在拉格朗日坐标下模拟流场中离散相的第二相;← 离散相模型解决的问题:煤粉燃烧、颗粒分离、喷雾干燥、液体燃料的燃烧等;← 应用范←围:FLUENT中的离散相模型假定第二相体积分数一般说来要小于10-12%(但颗粒质量承载率可以大于10-12%,即可模拟离散相质量流率等/大于连续相的流动);不适用于模拟在连续相中无限期悬浮的颗粒流问题,包括:搅拌釜、流化床等; 颗粒-颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑;← 湍流中颗粒处理的两种模型:Stochastic Tracking,应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响;Cloud Tracking,运用统计方法来跟踪颗粒围绕某一平均轨道的湍流扩散。通过计算颗粒的系统平均运动方程得到颗粒的某个“平均轨道”← 二.多相流模型 FLUENT中提供的模型: VOF模型(Volume of Fluid Model)← 混合模型(Mixture Model)← 欧拉模型(Eulerian Model)← 1.VOF模型(Volume of Fluid Model) VOF模型用来处理没有相互穿插的多相流问题,在处理两相流中,假设计算的每个控制容积中第一相的体积含量为α1,如果α1=0,表示该控制容积中不含第一相,如果α1=1,则表示该控制容积中只含有第一相,如果0←<α1<1,表示该控制容积中有两相交界面; VOF方法是用体积率函数表示流体自由面的位置和流体所占的体积,其方法占内存小,是一种简单而有效的方法。← 2.混合模型(Mixture Model) 用混合特性参数描述的两相流场的场方程组称为混合模型;← 考虑了界面传递特性以及两相间的扩散作用和脉动作用;使用了滑移速度的概念,允许相以不同的速度运动;← 用于模拟各相有不同速度的多相流;也用于模拟有强烈耦合的各向同性多相流和各相以相同速度运动的多相流;← 缺点:界面特性包括不全,扩散和脉动特性难于处理。← 3.欧拉模型(Eulerian Model) 欧拉模型指的是欧拉—欧拉模型;← 把颗粒和气体看成两种流体,空间各点都有这两种流体各自不同的速度、温度和密度,这些流体其存在在同一空间并相互渗透,但各有不同的体积分数,相互间有滑移;←

fluent全攻略(探索阶段)

GAMBIT使用说明 GAMBIT是使用FLUENT进行计算的第一个步骤。在GAMBIT 中我们将完成对计算模型的基本定义和初始化,并输出初始化结果供FLUENT的计算需要。以下是使用GAMBIT的基本步骤。 1.1定义模型的基本几何形状 如左图所示的按钮就是用于构造模型的基本几何形状的。当按下这个按钮时,将出现 如下5个按钮,它们分别是用以定义点、线、面、体的几何形状的。 值得注意的是我们定义这些基本的几何元素的一般是依照以下的顺序: 点——线(两点确定一线)——面(3线以上确定一面)——体(3面以上确定体)对各种几何元素的操作基本方式是:首先选中所要进行的操作,再定义完成操作所要的其他元素,作后点“APPL Y”按钮完成操作。以下不一一重复。 下面我们分别介绍各个几何元素的确定方法: 1.1.1点的操作 对点的操作在按下点操作按钮后进行(其他几何元素的操作也是这样)。点有以下几种主要操作 定义点的位置按钮,按下后出现下面对话框 Coordinate Sys.:用以选择已有坐标系中进行当前操 作的坐标系 Type:可以选择3种相对坐标系为当前坐标系:笛卡 儿坐标、柱坐标、球坐标。 以下通过在Global 中直接输入点的x、y、z值定义点, 注意这里的坐标值是绝对坐标值,而Local中输入的是相 对坐标值,一般我们使用绝对坐标值。 Label:为所定义的点命名。 在完成以上定义后就可以通过进行这个点 的定义,同时屏幕左半部的绘图区中将出现被定义的点。 用关闭此对话框。 查看所有点的几何参数按钮(在以后的操作中也可以查看其他元素的几何参数) 在Vertices栏中选择被查询的点,有两种选择方式(其他几 何元素的选择与此类似): ①按住shift键的同时用鼠标左键取点

Fluent经典问题及答疑2

Fluent经典问题及答疑2 51 对于出口有回流的问题,在出口应该选用什么样的边界条件(压力出口边界条件、质量出口边界条件等)计算效果会更好?(#42) 52 对于不同求解器,离散格式的选择应注意哪些细节?实际计算中一阶迎风差分与二阶迎风差分有什么异同?(#69) 53 对于FLUENT的耦合解算器,对时间步进格式的主要控制是Courant数(CFL),那么Courant 数对计算结果有何影响?(#43) 54 在分离求解器中,FLUENT提供了压力速度耦和的三种方法:SIMPLE,SIMPLEC及PISO,它们的应用有什么不同?(#44) 55 对于大多数情况,在选择选择压力插值格式时,标准格式已经足够了,但是对于特定的某些模型使用其它格式有什么特别的要求? (#60) 56 计算流体力学中在设定初始条件和边界条件的时候总是要先选择一组湍流参数,并给出其初值。如何选择并给出这些初值呢?有什么经验公式或者别的好的办法吗?(#73) 57 讨论在数值模拟过程中采用四面体网格计算效果好,还是采用六面体网格更妙呢?(#70) 58 如何将自己用C语言编辑的程序导入到FLUENT中?在利用UDF编写程序时需注意哪些问题?(#157) 59 在UDF中compiled型的执行方式和interpreted型的执行方式有什么不同?(#72) 60 在用gambit的时候,导入pro/e的stp文件后,在消去最短边的时候,有些最短边不能消去,其是空间线段,用面merge的方法和连接点的方法都不行,请问该怎么消去这类短边?(#144) 61 FLUENT help和GAMBIT help能教会我们(特别是刚入门的新手)学习什么基本知识?(#126) 62 FLUENT如何做汽车外流场计算的模拟?并且怎么可以得到汽车的阻力系数和升力系数?(#170) 63 FLUENT模拟飞行器外部流场,最高MA多少时就不准确了?MA达到一定的程度做模拟需注意哪些问题?(#125) 64 在用gambit建模,保存成*.msh文件时总是出现No entity的错误:Continuum Entity fluid does not contain any valid entity and is not written! Boundary Entity wall does not contain any validentity and is not written! 不知道是什么问题?产生的原因是什么?如何解决?(#150) 65 在做燃烧模拟的时候,入口燃料温度定义为蒸发/离解开始时的温度(也就是,为离散相材料指定的蒸发温度“Vaporization Temperature”),这是指水分蒸发温度吗?一般是多少?(#196) 66 在计算煤粉燃烧时遇到这样的问题: Warning: volatile + combustible fraction for lignite is greater than 1.0shell conduction zones 如何解决? 67 FLUENT控制方程是无因次的还是有因次的?如果是无因次的,怎么无因次的? 68 做飞机设计时,经常计算一些翼型,可是经常出现计算出来的阻力是负值,出现负值究竟是什么原因,是网格的问题还是计算参数设置的问题?(#71) 69 FLUENT中的Turbulent intensify是如何定义的?该值应该是小于等于100%,可是我的计算中该值达到400%,不知为何? 70 边界条件中湍流强度怎么设置:入口边界条件中的湍流强度和出口边界条件中的回流湍流强度怎么设置?是取默认值10%吗?(#135) 71 关于Injection中的Total Flow rate:injection 选surface,此时选了好几个面(面积不一定完全相同,但颗粒的入口速度相同),那Total Flow Rate 是指几个面的总流量还是某一个面的啊?只能处理完全相同的面吗?(#160) 72 FLUENT中能不能做插值:在ansys中的模型节点坐标和FLUENT中模型的节点坐标不一致,能

FLUENT菜鸟入门-不可不知的50个经典问题

Fluent必知的一些基本概念! 连续性方程不收敛是怎么回事? 在计算过程中其它指数都收敛了,就continuity不收敛是怎么回事 这和Fluent程序的求解方法SIMPLE有关。SIMPLE根据连续方程推导出压力修正方法求解压力。由于连续方程中流场耦合项被过渡简化,使得压力修正方程不能准确反映流场的变化,从而导致该方程收敛缓慢。 你可以试验SIMPLEC方法,应该会收敛快些。 湍流与黏性有什么关系? 湍流和粘性都是客观存在的流动性质。 湍流的形成需要一定的条件,粘性是一切流动都具有的。 流体流动方程本身就是具非线性的。 NS方程中的粘性项就是非线性项,当然无粘的欧拉方程也是非线性的。 粘性是分子无规则运动引起的,湍流相对于层流的特性是由涡体混掺运动引起的。 湍流粘性是基于湍流体的parcel湍流混掺是类比于层流体中的分子无规则运动,只是分子无规则运动遥远弱些吧了。不过,这只是类比于,要注意他们可是具有不同的属性。粘性是耗散的根源,实际流体总是有耗散的。 而粘性是制约湍流的。 LANDAU说,粘性的存在制约了湍流的自由度。 湍流粘性系数和层流的是不一样的,层流的粘性系数基本可认为是常数,可湍流中层流底层中粘性系数很小,远小于层流时的粘性系数;而在过渡区,与之相当,在一个数量级;在充分发展的湍流区,又远大于层流时的粘性系数.这是鮑辛内斯克1987年提出的。 1 FLUENT的初始化面板中有一项是设置从哪个地方开始计算(compute from),选择从不同的边界开始计算有很大的区别吗?该怎样根据具体问题选择从哪里计算呢?比如有两个速度入口A和B,还有压力出口等等,是选速度入口还是压力出口?如果选速度入口,有两个,该选哪个呀?有没有什么原则标准之类的东西? 一般是选取ALL ZONE,即所有区域的平均处理,通常也可选择有代表性的进口(如多个进口时)进行初始化。对于一般流动问题,初始值的设定并不重要,因为计算容易收敛。但当几何条件复杂,而且流动速度高变化快(如音速流动),初始条件要仔细选择。如果不收敛,还应试验不同的初始条件,甚至逐次改变边界条件最后达到所要求的条件。 2 要判断自己模拟的结果是否是正确的,似乎解的收敛性要比那些初始条件和边界条件更重要,可以这样理解吗?也就是说,对于一个具体的问题,初始条件和边界条件的设定并不是唯一的,为了使解收敛,需要不断调整初始条件和边界条件直到解收敛为止,是吗?如果解收敛了,是不是就可以基本确定模拟的结果是正确的呢? 对于一个具体的问题,边界条件的设定当然是唯一的,只不过初始化时可以选择不同的初始条件(指定常流),为了使解的收敛比较好,我一般是逐渐的调节边界条件到额定值("额定值"是指你题目中要求的入口或出口条件,例如计算一个管内流动,要求入口压力和温度为10MPa和3000K,那么我开始叠代时选择入口压力和温度为1MPa和500K(假设,这看你自己问题了),等流场计算的初具规模、收敛的较好了,再逐渐调高压力和温度,经过好几次调节后最终到达额定值10MPa和3000K,这样比一开始就设为10MPa和3000K收敛的要好些)这样每次叠代可以比较容易收敛,每次调节后不用再初始化即自动调用上次的解为这次的初始解,然后继续叠代。即使解收敛了,这并

Ansys Workbench Fluid Flow(FLUENT)经典问题

1 对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢? 学习任何一个软件,对于每一个人来说,都存在入门的时期。认真勤学是必须的,什么是最好的学习方法,我也不能妄加定论,在此,我愿意将我三年前入门FLUENT心得介绍一下,希望能给学习FLUENT的新手一点帮助。 由于当时我需要学习FLUENT来做毕业设计,老师给了我一本书,韩占忠的《FLUENT流体工程仿真计算实例与应用》,当然,学这本书之前必须要有两个条件,第一,具有流体力学的基础,第二,有FLUENT 安装软件可以应用。然后就照着书上二维的计算例子,一个例子,一个步骤地去学习,然后学习三维,再针对具体你所遇到的项目进行针对性的计算。不能急于求成,从前处理器GAMBIT,到通过FLUENT进行仿真,再到后处理,如TECPLOT,进行循序渐进的学习,坚持,效果是非常显著的。如果身边有懂得FLUENT 的老师,那么遇到问题向老师请教是最有效的方法,碰到不懂的问题也可以上网或者查找相关书籍来得到答案。另外我还有本《计算流体动力学分析》王福军的,两者结合起来学习效果更好。 2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。 https://www.wendangku.net/doc/1615016269.html,/dvbbs/viewFile.asp?BoardID=61&ID=1411 A.理想流体(Ideal Fluid)和粘性流体(Viscous Fluid): 流体在静止时虽不能承受切应力,但在运动时,对相邻的两层流体间的相对运动,即相对滑动速度却是有抵抗的,这种抵抗力称为粘性应力。流体所具备的这种抵抗两层流体相对滑动速度,或普遍说来抵抗变形的性质称为粘性。粘性的大小依赖于流体的性质,并显著地随温度变化。实验表明,粘性应力的大小与粘性及相对速度成正比。当流体的粘性较小(实际上最重要的流体如空气、水等的粘性都是很小的),运动的相对速度也不大时,所产生的粘性应力比起其他类型的力如惯性力可忽略不计。此时我们可以近似地把流体看成无粘性的,这样的流体称为理想流体。十分明显,理想流体对于切向变形没有任何抗拒能力。这样对于粘性而言,我们可以将流体分为理想流体和粘性流体两大类。应该强调指出,真正的理想流体在客观实际中是不存在的,它只是实际流体在某些条件下的一种近似模型。 B.牛顿流体(Newtonian Fluid)和非牛顿流体(non-Newtonian Fluid): 日常生活和工程实践中最常遇到的流体其切应力与剪切变形速率符合下式的线性关系,称为牛顿流体。而切应力与变形速率不成线性关系者称为非牛顿流体。图2-1(a)中绘出了切应力与变形速率的关系曲线。其中符合上式的线性关系者为牛顿流体。其他为非牛顿流体,非牛顿流体中又因其切应力与变形速率关系特点分为膨胀性流体(Dilalant),拟塑性流体(Pseudoplastic),具有屈服应力的理想宾厄流体(Ideal Bingham Fluid)和塑性流体(Plastic Fluid)等。通常油脂、油漆、牛奶、牙膏、血液、泥浆等均为非牛顿流体。非牛顿流体的研究在化纤、塑料、石油、化工、食品及很多轻工业中有着广泛的应用。图2-1(b)还显示出对于有些非牛顿流体,其粘滞特性具有时间效应,即剪切应力不仅与变形速率有关而且与作用时间有关。当变形速率保持常量,切应力随时间增大,这种非牛顿流体称为震凝性流体(Rheopectic Fluid)。当变形速率保持常量而切应力随时间减小的非牛顿流体则称为触变性流体(Thixotropic Fluid)。 C.可压缩流体(Compressible Fluid)和不可压缩流体(Incompressible Fluid):

相关文档