FLUENT培训教材02求解器基础

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中的求解器、算法和离散方法

FLUENT中的求解器、算法和离散方法 作为一个非科班出身的CFD工程师,一开始常常被CFD软件里各种概念搞的晕头转向。最近终于静下心来看了看CFD理论的书,理清了一些概念。就此写一遍博文,顺便整理一下所学内容。 I 求解器: FLUENT中求解器的选择在如下图所示界面中设置: FLUENT中的求解器主要是按照是否联立求解各控制方程来区分的,详见下图:

II 算法: 算法是求解时的策略,即按照什么样的方式和步骤进行求解。FLUENT中算法的选择在如下图所示的界面中设置:

这里简单介绍一下SIMPLE、SIMPLEC、PISO等算法的基本思想和适用范围。 SIMPLE算法:基本思想如前面讲求解器的那张图中解释分离式求解器的例子所示的一样,这里再贴一遍: 1.假设初始压力场分布。 2.利用压力场求解动量方程,得到速度场。 3.利用速度场求解连续性方程,使压力场得到修正。 4.根据需要,求解湍流方程及其他方程 5.判断但前计算是否收敛。若不收敛,返回第二步。 简单说来,SIMPLE算法就是分两步走:第一步预测,第二步修正,即预测-修正。 SIMPLC算法:是对SIMPLE算法的一种改进,其计算步骤与SIMPLE算法相同,只是压力修正项中的一些系数不同,可以加快迭代过程的收敛。 PISO算法:比SIMPLE算法增加了一个修正步,即分三步:第一步预测,第二步修正得到一个修正的场分布,第三步在第二步基础上在进行一侧修正。即预测-修正-修正。PISO算法在求解瞬态问题时有明显优势。对于稳态问题可能SIMPLE 或SIMPLEC更合适。 如果你实在不知道该如何选择,就保持FLUENT的默认选项好了。因为默认选项可以很好解决70%以上的问题,而且对于大部分出了问题的计算来说,也很少是因为算法选择不恰当所致。 III 离散方法: 离散方法是指按照什么样的方式将控制方程在网格节点离散,即将偏微分格式的控制方程转化为各节点上的代数方程组。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技术基础与应用实例 4.2.2 fluent数值模拟步骤简介 主要步骤: 1、根据实际问题选择2D或3Dfluent求解器从而进行数值模拟。 2、导入网格(File→Read→Case,然后选择有gambit导出的.msh文件) 3、检查网格(Grid→Check)。如果网格最小体积为负值,就要重新 进行网格划分。 4、选择计算模型。 5、确定流体物理性质(Define→Material)。 6、定义操作环境(Define→operating condition) 7、制定边界条件(Define→Boundary Conditions) 8、求解方法的设置及其控制。 9、流场初始化(Solve→Initialize) 10、迭代求解(Solve→Iterate) 11、检查结果。 12、保存结果,后处理等。 具体操作步骤: 1、fluent2d或3d求解器的选择。 2、网格的相关操作 (1)、读入网格文件 (2)、检查网格文件 文件读入后,一定要对网格进行检查。上述的操作可以得到网格信息,从中看出几何区域的大小。另外从minimum volume 可以知道最小网格的体积,若是它的值大于零,网格可以用于计算,否则就要重新划 分网格。 (3)、设置计算区域 在gambit中画出的图形是没有单位的,它是一个纯数量的模型。故 在进行实际计算的时候,要根据实际将模型放大或缩小。方法是改变fluent总求解器的单位。 (4)、显示网格。 Display→Grid 3、选择计算模型

(1)、基本求解器的定义 Define→Models→Solver Fluent中提供了三种求解方法: ·非耦合求解 segregated ·耦合隐式求解 coupled implicit ·耦合显示求解 coupled explicit 非耦合求解方法主要用于不可压缩流体或者压缩性不强的流体。 耦合求解方法用在高速可压缩流体 fluent默认设置是非耦合求解方法,但对于高速可压缩流动,有强的体积力(浮力或离心力)的流动,求解问题时网格要比较密集,建 议采用耦合隐式求解方法。耦合能量和动量方程,可以较快的得到收敛值。耦合隐式求解的短板:运行所需要的存比较大。若果必须要耦合求解而机器存不够用,可以考虑采用耦合显示求解方法。盖求解方法也耦合了动量,能量和组分方程,但是存却比隐式求解方法要小。 需要指出的是,非耦合求解器的一些模型在耦合求解器里并不一定都有。耦合求解器里没有的模型包括:多相流模型、混合分数/PDF燃烧模型、预混燃烧模型。污染物生成模型、相变模型、Rosseland辐射模型、确定质量流率的周期性流动模型和周期性换热模型。 %%%有点重复,但是可以看看加深理解 Fluent提供三种不同的求解方法;分离解、隐式耦合解、显示耦合解。分理解和耦合解的主要区别在于:连续方程、动量方程、能量方程和 组分方程解的步骤不同。 分离解按照顺序解,耦合解是同时解。两种解法都是最后解附加的标量方程。隐式解和显示解的区别在于线性耦合方程的方式不同。 Fluent默认使用分离求解器,但是对于高速可压流动,强体积力导致 的强烈耦合流动(流体流动耦合流体换热耦合流体的混合,三者相互耦合的过程—文档整理者注)(浮力或者旋转力),或者在非常精细的网格上的流动,需要考虑隐式解。这一解法耦合了流动和能量方程, 收敛很快。%%% (2)、其他求解器的选择 在实际问题中,除了要计算流场,有时还要计算温度场或者浓度场等,因此还需要其他的模型。主要的模型有: Multiphase(多相流动)viscous(层流或湍流)energy(是否考虑传热)species(反应及其传热相关) (3)操作环境的设置 Define→operation→condition

利用FLUENT 3D求解器求解

利用FLUENT 3D求解器求解 一、在FLUENT中读入网格文件,检查网格并定义长度单位 1、启动FLUENT,进入3D模式 操作:开始→程序→FLUENT→3d→Run,进入FLUENT。 2、读入网格文件 操作:File→Read→Case,选择在Gambit中绘制的网格文件.msh文件,点击OK完成数据读入。 3、调整网格尺寸比例 操作:Grid→Scale 打开“Scale Grid”对话框 (1)在Units Conversion 下的Grid Was Created In 右侧列表中选择合适的单位如:cm (在gambit中一般是以m为单位,要转化成fluent对应的单位cm); (2)点击Change length Units: 此时左侧的Scale Factors下的X,Y,Z项都变为0.01。 (3)点击下边的Scale按钮:此时,Domain Extents下的单位由m变成cm;并给出区域的范围; (4)点击Close关闭对话框。 4、检查网格 操作:Grid→Check Fluent会对网格进行各种检查并在信息反馈窗口显示检查过程和结果,其中要注意保持最小体积为正值。 5、显示网格 操作:Display→Grid 打开网格显示对话框后,点击Display。 注意:用鼠标右键点击边界线,则在信息反馈窗口内将显示此边界的类型等信息。也可用此方法检查任何内部节点和网格线的信息。 二、创建计算模型 1、设置求解器 操作:Define→Models→Solver (1)在Solver项选择Segregated; (2)在Formulation项选择Implicit; (3)在Space项选择3D; (4)在Time项选择Unsteady; (5)Velocity Formulation,Unsteady Formulation保持默认值; (6)点击OK。 2、定义多相流模型 操作:Define→Models→Multiphase (1)在Model项选择Volume of Fluid; (2)在Number of Phase下选2; (3)在VOF Scheme项选择Geo-Reconstruct,Courant Number保持默认值; (fluent6.3.26里边VOF Scheme选expicity。) 在Body Force Formulation项选择Implicit Body Force; (4)点击OK。 3 、设置标准的k﹣ε湍流模型 操作:Define→Models→Viscous

FLUENT中求解方法的区别与选择

FLUENT中求解方法的区别与选择 2012-02-11 12:08:37| 分类:数值模拟|字号订阅 1.非耦合求解( Segregated );2.耦合隐式求解( Coupled Implicit );3.耦合显式求 解( Coupled Explicit ) 非耦合求解方法主要用于不可压缩或压缩性不强的流体流动。耦合求解则可以用在高速可压缩流动。FLUENT默认设置是非耦合求解,但对于高速可压流动,有强的体积力(浮力或离心力)的流动,求解问题时网格要比较密,建议采用耦合隐式求解方法,可以耦合求解能量和动量方程,能比较快地得到收敛解。缺点是需要的内存比较大(是非耦合求解迭代时间的1.5-2倍)。如果必须要耦合求解,但是你的机器内存不够,这时候可以考虑用耦合显式解法器求解问题。该解法器也耦合了动量,能量及组分方程,但内存却比隐式求解方法小。缺点是收敛时间比较长。 这里需要指出的是非耦合求解的一些模型在耦合求解解法器里并不都有。耦合解法器没有的模型包括:多相流模型,混合分数/PDF燃烧模型,预混燃烧模型,污染物生成模型,相变模型,Rosseland辐射模型,确定质量流率的周期性流动模型及周期性换热模型等。 隐式( Implicit ):对于给定变量,单元内的未知值用邻近单元的已知和未知值计算得出。因此,每一个未知值会在不止一个方程中出现,这些方程必须同时解来给出未知量。 显式( Explicit ):对于给定变量,每一个单元内的未知量用只包含已知量的关系式计算得到。因此未知量只在一个方程中出现,而且每一个单元内的未知量的方程只需解一次就可以给出未知量的值。 一阶迎风格式( First Order Upwind ):当需要一阶精度时,我们假定描述单元内变量平均值的单元中心变量就是整个单元内各个变量的值,而且单元表面的量等于单元内的量。因此,当选择一阶迎风格式时,表面值被设定等于迎风单元的单元中心值。 二阶迎风格式( Second Order Upwind ):当需要二阶精度时,使用多维线性重建方法来计算单元表面处的值。在这种方法中,通过单元中心解在单元中心处的泰勒展开来实现单元表面的二阶精度值。因此,当使用二阶迎风格式时,用下面的方程来计算表面值; QUICK格式:对于四边形和六面体网格,我们可以确定它们唯一的上游和下游表面以及单元。FLUENT还提供了计算对流变量在表面处高阶值的QUICK格式。QUICK类型的格式是通过变量的二阶迎风与中心插值加上适当的权因子得到的; 亚松驰( Under-Relaxation ):由于FLUENT所解方程组的非线性,我们有必要控制的变化。一般用亚松驰方法来实现控制,该方法在每一部迭代中减少了f的变化量。亚松驰最简单的形式为:单元内变量f等于原来的值f_old加上亚松驰因子a与f变化的积. SIMPLE:SIMPLE算法使用压力和速度之间的相互校正关系来强制质量守恒并获取压力场。

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可能会得到更好的结果,尤其是可以应用增加的亚松驰迭代时,具体介绍如下: 对于相对简单的问题(如:没有附加模型激活的层流流动),其收敛性已经被压力速

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 (流体常识及软件计算参数设置)

luent中一些问题----(目录) 1 如何入门 2 CFD计算中涉及到的流体及流动的基本概念和术语 2.1 理想流体(Ideal Fluid)和粘性流体(Viscous Fluid) 2.2 牛顿流体(Newtonian Fluid)和非牛顿流体(non-Newtonian Fluid) 2.3 可压缩流体(Compressible Fluid)和不可压缩流体(Incompressible Fluid) 2.4 层流(Laminar Flow)和湍流(Turbulent Flow) 2.5 定常流动(Steady Flow)和非定常流动(Unsteady Flow) 2.6 亚音速流动(Subsonic)与超音速流动(Supersonic) 2.7 热传导(Heat Transfer)及扩散(Diffusion) 3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不 同? 3.1 离散化的目的 3.2 计算区域的离散及通常使用的网格 3.3 控制方程的离散及其方法 3.4 各种离散化方法的区别 4 常见离散格式的性能的对比(稳定性、精度和经济性) 5 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么? 6 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难? 6.1 可压缩Euler及Navier-Stokes方程数值解 6.2 不可压缩Navier-Stokes方程求解 7 什么叫边界条件?有何物理意义?它与初始条件有什么关系? 8 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别? 9 在网格生成技术中,什么叫贴体坐标系?什么叫网格独立解? 10 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量?及其在做网格时大致注意到哪些细节? 11 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢? 12 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理? b、计算域内的内部边界如何处理(2D)? 13 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些? 14 20 何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT是怎样使用区域的? 15 21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收

fluent计算讨论

建议大家一起讨论一下湍流边界条件该如何设置 本人也是刚刚接触Fluent系列软件不久,在应用它来求解CFD问题时遇到了不少问题,也得到了很多宝贵经验,其中湍流边界条件的设置就是一个很棘手的问题。 最近对该问题总结经验如下: 在入口、出口或远场边界流入流域的流动,FLUENT需要指定输运标量的值。本节描述了对于特定模型需要哪些量,并且该如何指定它们。也为确定流入边界值最为合适的方法提供了指导方针。 使用轮廓指定湍流参量 在入口处要准确的描述边界层和完全发展的湍流流动,你应该通过实验数据和经验公式创建边界轮廓文件来完美的设定湍流量。如果你有轮廓的分析描述而不是数据点,你也可以用这个分析描述来创建边界轮廓文件,或者创建用户自定义函数来提供入口边界的信息。一旦你创建了轮廓函数,你就可以使用如下的方法: λSpalart-Allmaras模型:在湍流指定方法下拉菜单中指定湍流粘性比,并在在湍流粘性比之后的下拉菜单中选择适当的轮廓名。通过将m_t/m和密度与分子粘性的适当结合,F LUENT为修改后的湍流粘性计算边界值。 λk-e模型:在湍流指定方法下拉菜单中选择K和Epsilon并在湍动能(Turb. Kinetic E nergy)和湍流扩散速度(Turb. Dissipation Rate)之后的下拉菜单中选择适当的轮廓名。λ雷诺应力模型:在湍流指定方法下拉菜单中选择K和Epsilon并在湍动能(Turb. Kin etic Energy)和湍流扩散速度(Turb. Dissipation Rate)之后的下拉菜单中选择适当的轮廓名。在湍流指定方法下拉菜单中选择雷诺应力部分,并在每一个单独的雷诺应力部分之后的下拉菜单中选择适当的轮廓名。 湍流量的统一说明 在某些情况下流动流入开始时,将边界处的所有湍流量指定为统一值是适当的。比如说,在进入管道的流体,远场边界,甚至完全发展的管流中,湍流量的精确轮廓是未知的。 在大多数湍流流动中,湍流的更高层次产生于边界层而不是流动边界进入流域的地方,因此这就导致了计算结果对流入边界值相对来说不敏感。然而必须注意的是要保证边界值不是非物理边界。非物理边界会导致你的解不准确或者不收敛。对于外部流来说这一特点尤其突出,如果自由流的有效粘性系数具有非物理性的大值,边界层就会找不到了。 你可以在使用轮廓指定湍流量一节中描述的湍流指定方法,来输入同一数值取代轮廓。你也可以选择用更为方便的量来指定湍流量,如湍流强度,湍流粘性比,水力直径以及湍流特征尺度,下面将会对这些内容作一详细叙述。 湍流强度I定义为相对于平均速度u_avg的脉动速度u^'的均方根。 小于或等于1%的湍流强度通常被认为低强度湍流,大于10%被认为是高强度湍流。从外界,测量数据的入口边界,你可以很好的估计湍流强度。例如:如果你模拟风洞试验,自由流的湍流强度通常可以从风洞指标中得到。在现代低湍流风洞中自由流湍流强度通常低到0.0

FLUENT学习经验总结(狠珍贵,学长传授)

1对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢? 答:学习任何一个软件,对于每一个人来说,都存在入门的时期。认真勤学是必须的,什么是最好的学习方法,我也不能妄加定论,在此,我愿意将我三年前入门FLUENT心得介绍一下,希望能给学习FLUENT的新手一点帮助。 由于当时我需要学习FLUENT来做毕业设计,老师给了我一本书,韩占忠的《FLUENT流体工程仿真计算实例与应用》,当然,学这本书之前必须要有两个条件,第一,具有流体力学的基础,第二,有FLUENT 安装软件可以应用。然后就照着书上二维的计算例子,一个例子,一个步骤地去学习,然后学习三维,再针对具体你所遇到的项目进行针对性的计算。不能急于求成,从前处理器GAMBIT,到通过FLUENT进行仿真,再到后处理,如TECPLOT,进行循序渐进的学习,坚持,效果是非常显著的。如果身边有懂得FLUENT的老师,那么遇到问题向老师请教是最有效的方法,碰到不懂的问题也可以上网或者查找相关书籍来得到答案。另外我还有本《计算流体动力学分析》王福军的,两者结合起来学习效果更好。 2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。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): 在流体的运动过程中,由于压力、温度等因素的改变,流体质点的体积(或密度,因质点的质量一定),或多或少有所改变。流体质点的体积或密度在受到一定压力差或温度差的条件下可以改变的这个性质称为压缩性。真实流体都是可以压缩的。它的压缩程度依赖于流体的性质及外界的条件。例如水在100个大气压下,容积缩小0.5%,温度从20°变化到100°,容积降低4%。因此在一般情况下液体可以近似地看成不可压的。但是在某些特殊问题中,例如水中爆炸或水击等问题,则必须把液体看作是可压缩的。气体的压缩性比液体大得多,所以在一般情形下应该当作可压缩流体处理。但是如果压力差较小,运动速度较小,并且没有很大的温度差,则实际上气体所产生的体积变化也不大。此时,也可以近似地将气体视为不可压缩的。 在可压缩流体的连续方程中含密度,因而可把密度视为连续方程中的独立变量进行求解,再根据气体的状态方程求出压力。不可压流体的压力场是通过连续方程间接规定的。由于没有直接求解压力的方程,不可压流体的流动方程的求解具有其特殊的困难。 D. 层流(Laminar Flow)和湍流(Turbulent Flow):

FLUENT知识点解读(良心出品必属精品)

一、基本设置 1.Double Precision的选择 启动设置如图,这里着重说说Double Precision(双精度)复选框,对于大多数情况,单精度求解器已能很好的满足精度要求,且计算量小,这里我们选择单精度。然而对于以下一些特定的问题,使用双精度求解器可能更有利[1]。 a.几何特征包含某些极端的尺度(如非常长且窄的管道),单精度求解器可能不能足够精确地表达各尺度方向的节点信息。 b.如果几何模型包含多个通过小直径管道相互连接的体,而某一个区域的压力特别大(因为用户只能设定一个总体的参考压力位置),此时,双精度求解器可能更能体现压差带来的流动(如渐缩渐扩管的无粘与可压缩流动模拟)。 c.对于某些高导热系数比或高宽纵比的网格,使用单精度求解器可能会遇到收敛性不佳或精确度不足不足的问题,此时,使用双精度求解器可能会有所帮助。 [1] 李鹏飞,徐敏义,王飞飞.精通CFD工程仿真与案例实战:FLUENT GAMBIT ICEM CFD Tecplot[M]. 北京,人民邮电出版社,2011:114-116

2.网格光顺化 用光滑和交换的方式改善网格:通过Mesh下的Smooth/Swap来实现,可用来提高网格质量,一般用于三角形或四边形网格,不过质量提高的效果一般般,影响较小,网格质量的提高主要还是在网格生成软件里面实现,所以这里不再用光滑和交换的方式改善网格,其原理可参考《FLUENT全攻略》(已下载)。 3.Pressure-based与Density-based 求解器设置如图。下面说一说Pressure-based和Density-based 的区别: Pressure-Based Solver是Fluent的优势,它是基于压力法的求解器,使用的是压力修正算法,求解的控制方程是标量形式的,擅长求解不可压缩流动,对于可压流动也可以求解;Fluent 6.3以前的版本求解器,只有Segregated Solver和Coupled Solver,其实也是Pressure-Based Solver的两种处理方法;

Fluent求解参数设置

求解参数设置(Solution Methods/Solution Controls): 在设置完计算模型和边界条件后,即可开始求解计算了,因为常会出现求解不收敛或者收敛速度很慢的情况,所以就要根据具体的模型制定具体的求解策略,主要通过修改求解参数来完成。在求解参数中主要设置求解的控制方程、选择压力速度耦合方法、松弛因子、离散格式等。 在VOF模型中,PISO比较适合于不复杂的流体,SIMPLE和SIMPLEC适合于可压缩的流体或者处于封闭域中的流体。 ? 求解的控制方程: 在求解参数设置中,可以选择所需要求解的 控制方程。可选择的方程包括Flow(流动方程)、 Turbulence(湍流方程)、Energy(能量方程)、 V olume Fraction(体积分数方程)等。在求解过程 中,有时为了得到收敛的解,先关闭一些方程, 等一些简单的方程收敛后,再开启复杂的方程 一起计算。 ? 选择压力速度耦合方法: 在基于压力求解器中,FLUENT提供了压力 速度耦合的4种方法,即SIMPLE、 SIMPLEC(SIMPLE.Consistent)、PISO以及 Coupled。定常状态计算一般使用SIMPLE或者 SIMPLEC方法,对于过渡计算推荐使用PISO 方法。PISO方法还可以用于高度倾斜网格的定 常状态计算和过渡计算。需要注意的是压力速 度耦合只用于分离求解器,在耦合求解器中不 可以使用。 在FLUENT中,可以使用标准SIMPLE算法和SIMPLEC算法,默认是SIMPLE算法,但对于许多问题如果使用SIMPLEC可能会得到更好的结果,尤其是可以应用增加的亚松弛迭代时。 对于相对简单的问题(如没有附加模型激活的层流流动),其收敛性可以被压力速度耦合所限制,用户通常可以使用SIMPLEC算法很快得到收敛解。在SIMPLEC算法中,压力校正亚松弛因子通常设为1.0,它有助于收敛,但是,在有些问题中,将压力校正松弛因子增加到1.0可能会导致流动不稳定,对于这种情况,则需要使用更为保守的亚松弛或者使用SIMPLE 算法。对于包含湍流或附加物理模型的复杂流动,只要用压力速度耦合做限制,SIMPLEC 就会提高收敛性,它通常是一种限制收敛性的附加模拟参数,在这种情况下,SIMPLE和SIMPLEC会给出相似的收敛速度。 对于所有的过渡流动计算,推荐使用PISO算法邻近校正。它允许用户使用大的时间步,而且对于动量和压力都可以使用亚松弛因子1.0。对于定常状态问题,具有邻近校正的PISO 并不会比具有较好的亚松弛因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。 当使用PISO邻近校正时,对所有方程都推荐使用亚松弛因子为1.0或者接近1.0。如果只对高度扭曲的网格使用PISO倾斜校正,则要设定动量和压力的亚松弛因子之和为1.0(例如,压力亚松弛因子0.3,动量亚松弛因子0.7)。

fluent经验

Fluent 问题集锦 问题1: 如果体网格做好后,感觉质量不好,然后将体网格删除,在其面上重新作网格,结果发现网格都脱离面,不再附体了,比其先前的网格质量更差了. 原因: 删除体网格时,也许连同较低层次的网格都删除了.上面的脱离面可能是需要的体的面. 解决方法: 重新生成了面,在重新划分网格 问题2: 在gambit下做一虚的曲面的网格,结果面上的网格线脱离曲面,由此产生的体网格出现负体积. 原因: 估计是曲面扭曲太严重造成的 解决方法: 可以试试分区域划分体网格,先将曲面分成几个小面,生成各自的面网格,再划体网格。 问题3: 当好网格文件的时候,并检查了网格质量满足要求,但输出*.msh时报错误. 原因: 应该不是网格数量和尺寸.可能是在定义边界条件或continuum type时出了问题. 解决方法: 先把边界条件删除重新导出看行不行.其二如果有两个几何信息重合在一起, 也可能出现上诉情况,将几何信息合并掉. 问题4: 当把两个面(其中一个实际是由若干小面组成,将若干小面定义为了group了)拼接在一起,也就是说两者之间有流体通过,两个面各属不同的体,网格导入到fluent时,使用interface时出现网格check的错误,将interface的边界条件删除,就不会发生网格检查的错误.如何将两个面的网格相连. 原因: interface后的两个体的交接面,fluent以将其作为内部流体处理(非重叠部分默认为wall,合并后网格会在某些地方发生畸变,导致合并失败.也可能准备合并的两个面几何位臵有误差,应该准确的在同一几何位臵(合并的面大小相等时),在合并之前要合理分块。 解决方法: 为了避免网格发生畸变(可能一个面上的网格跑到另外的面上了),可以一面网格粗,一面网格细,避免; 再者就是通过将一个面的网格直接映射到另一面上的,两个面默认为interior.也可以将网格拼接一起. 上述语言有些模糊不清,仅供参考,并希望高手批评指正,^_^

求解器的使用

求解器的使用 FLUENT提供了三种不同的求解器 Segregated,coupled implicit,coupled explicit(显式格式主要用于激波等波动解的捕捉问题) 传统上,分离解法(segregated)主要用于不可压缩以及适度压缩的流动中。相反,耦合算法是为高速可压流体设计的。 默认情况下,fluent使用分离求解器。对于高速可压流体,与很强的体积力高度耦合的流动,或者是在非常精确的网格上求解流动情况,可以考虑使用耦合隐式算法代替。 对于需要使用耦合隐式算法(coupled implicit)的case,如果电脑没有足够的内存,可以使用分离解法(segregated)或者耦合显式算法(coupled explicit)代替。(显示算法节约内存,但是需要更多的计算步数达到收敛。) 选择离散格式 1.一阶迎风格式v.s. 二阶迎风格式 当流动与网格匹配(校准)时,一阶迎风格式是可以接受的。对于三角形和四面体网格,由于流动不会与网格匹配,通常使用二阶离散格式会得到更准确的结果。对于四边形/六面体网格,使用二阶离散格式会取得更好的结果,尤其是复杂的流动情况。对于大多数情况,可以在计算初始,使用二阶的离散格式。然而在一些情况下,可以开始使用一阶的离散格式然后在一些计算之后转变为二阶格式。例如,如果正在运行一个高马赫数的流动计算,这个的初始解与期望的解相差很大,最终,如果二阶离散格式很难收敛,应该尝试使用一阶离散格式。 2.Quick格式v.s. Upwind(Quick格式适用于网络结构,流动方向与网格一致,对于非结 构网格推荐使用二阶迎风) 对于在四边形或者六面体网格中的旋转或者回旋流,Quick离散格相比于二阶离散格式可以提供更准确的结果。对于存在震动的可压缩流动(网格为四边形,六面体或者混合网格),推荐对所有的变量使用Quick离散格式,包括密度。 3.中心差分格式v.s. 迎风格式 当使用LES湍流模型时,是可以使用中心差分格式的,并且只有当网格间距足够好,以至于局部的Peclet数的大小小于1时才可以使用。 4.power法则(power law)v.s. 迎风 power法则是可以使用的,但是总体上产生与一阶格式相同的准确度 选择压力离散格式 当使用分离式求解器时(segregated),可以使用多种压力离散格式。在大多数情况下,(默认的)标准的格式是可以接受的,但是一些类型的模型会在其他格式中取得更好的结果:对于包含大的体积力的问题,推荐使用体积力重量(body-force-weighted)格式。 对于包含高的漩涡数,高瑞利数的自然对流,高速的旋转流动,包含多孔介质的流动,在高度弯曲区域中的流动,使用PRESTO!格式 对于可压流,应使用二阶格式 当其他格式并不合适时,使用二阶格式以提高准确性。 选择密度离散格式(求解一个单相可压流动) 如果计算一个含有震动的可压流动时,一阶迎风格式可以平滑震动;对于这样的流动,应该使用二阶迎风或者Quick格式

相关文档
最新文档