fluent基础(学习入门篇~).doc

上传人:小** 文档编号:630103 上传时间:2019-04-22 格式:DOC 页数:14 大小:220.50KB
返回 下载 相关 举报
fluent基础(学习入门篇~).doc_第1页
第1页 / 共14页
fluent基础(学习入门篇~).doc_第2页
第2页 / 共14页
点击查看更多>>
资源描述

《fluent基础(学习入门篇~).doc》由会员分享,可在线阅读,更多相关《fluent基础(学习入门篇~).doc(14页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。

1、-_1 单精度与双精度的区别大多数情况下,单精度解算器高效准确,但是对于某些问题使 用双精度解算器更合适。下面举几个例子: 如果几何图形长度尺度相差太多(比如细长管道) ,描述节点 坐标时单精度网格计算就不合适了;如果几何图形是由很多层小直 径管道包围而成(比如:汽车的集管)平均压力不大,但是局部区 域压力却可能相当大(因为你只能设定一个全局参考压力位置) , 此时采用双精度解算器来计算压差就很有必要了。对于包括很大热 传导比率和(或)高比率网格的成对问题,如果使用单精度解算器 便无法有效实现边界信息的传递,从而导致收敛性和(或)精度下 降2 分离解与耦合解的区别选择解的格式 FLUENT 提

2、供三种不同的解格式:分离解;隐式耦合解;显式耦 合解。三种解法都可以在很大流动范围内提供准确的结果, 但是它 们也各有优缺点。 分离解和耦合解方法的区别在于, 连续性方程、动量方程、 能量方程以及组分方程的解的步骤不同, 分离解是按 顺序解,耦合解是同时解。两种解法都是最后解附加的标量方程 (比如:湍流或辐射)。隐式解法和显式解法的区别在于线化耦合方 程的方式不同。分离解以前用于 FLUENT 4 和 FLUENT/UNS,耦合显 式解以前用于 RAMPANT。分离解以前是用于不可压流和一般可压流 的。而耦合方法最初是用来解高速可压流的。现在,两种方法都适 用于很大范围的流动(从不可压到高速可

3、压),但是计算高速可压流 时耦合格式比分离格式更合适。FLUENT 默认使用分离解算器,但是 对于高速可压流(如上所述) ,强体积力导致的强烈耦合流动(比 如浮力或者旋转力), 或者在非常精细的网格上的流动, 你需要考 虑隐式解法。这一解法耦合了流动和能量方程, 常常很快便可以收 敛。 耦合隐式解所需要内存大约是分离解的 1.5 到 2 倍,选择时 可以通过这一性能来权衡利弊。在需要隐式耦合解的时候,如果计 算机的内存不够就可以采用分离解或者耦合显式解。 耦合显式解虽 然也耦合了流动和能量方程,但是它还是比耦合隐式解需要的内存 少,但是它的收敛性相应的也就差一些。 注意:注意: 分离解中提供的

4、几个物理模型, 在耦合解中是没有的:多项流模型; 混合组分/PDF 燃烧模型/预混合燃烧模型 /Pollutant formation models/相变模型/Rosseland 辐射模型/指 定质量流周期流动模型/周期性热传导模型。分离求解器是默认的-_3 欧拉方程用于解决无粘流动,4 操纵、相对、绝对压力对于不同的流动状态, 操作压力的指定以不同的方式影响你的 计算。在低马赫数流动中压力计算的数值截断的影响。在低马赫数 可压流动中, 全部的压降和绝对静压相比很小, 因此数值截断会 对其有很大的影响。比方说吧,考虑 M 1 的可压流动。压力变 化 D p 与动压头(1/2) c p M2 有

5、关,其中 p 是静压,c 是指定 的比热比。这就给出了 D p/p 和 M2 的关系式,以至于 M 0 时 D p/p0。因此,除非给予足够的注意,否则低马赫数流 动计算结果往往很容易会受到截断误差截断误差的影响。 操作压力、标准压力和绝对压力 FLUENT 通过从绝对压力中减去操作压力(一般说来大的压力粗 略的等于流动中绝对压力的平均值)来避免截断误差(见在低马赫 数流动中压力计算的数值截断的影响一节)产生的问题,并使用得 到的压力来计算,这个压力称作标准压力。下面是操作压力,标准 压力和绝对压力之间的关系式。绝对压力是操作压力和标准压力之 和: 你所指定的所有压力以及 FLUENT 所报告

6、和计算的压力都是 标准压力。 操作压力的意义 操作压力对于不可压理想气体流动来说是十分重要的, 因为它 直接决定了不可压理想气体定律所计算出来的密度, 不可压理想气 体定律计算密度的关系式为: r = (p_op/R T)。 因此,你必须保 证适当的设定操作压力。 操作压力在低马赫数可压流动中具有十分重要的意义, 因为它 在避免截断误差问题中扮演了重要的角色,如操作压力,标准压力 和绝对压力一节所述。同样地,你必须保证适当地设定操作压力。 对于高马赫数可压流动,操作压力的意义就不是很明显了。在 这种情况下,压力的变化比低马赫数可压流动中压力的变化大得多,因此截断误差不会产生什么实际的问题, 因

7、此也就不真正需要使 用标准压力。事实上,在这种计算中使用绝对压力通常会更方便。 因为 FLUENT 总是使用标准压力,所以你可以简单的设定操作压力 为零,而使标准压力和绝对压力相等。 如果密度假定为常数, 或者密度是从温度的轮廓函数中推导出 来, 那么根本就不使用操作压力。需要注意的是:默认的操作压力 为 101325 Pa.。 如何设定操作压力-_选择合适的操作压力的判据是基于流动马赫数的区域以及确定 密度的关系式。 例如: 如果你在不可压流动的计算中使用理想气 体定律(如自然对流问题) ,你应该使用平均流动压力的典型值。 下表是设定操作压力的推荐方法。请记住默认的操作压力为 101325

8、Pa。 你需要在操作压力面板中设定操作压力。菜单: Define/Operating Conditions.。 对于不包括任何压力边界的不可压流动,FLUENT 会在每次迭代 之后调节标准压力场以避免它浮动。这一操作是通过在(或接近) 参考压力位置的单元中使用的压力实现的。在完全的压力场中减去 单元内的压力值, 从而保证参考压力位置的标准压力总为零。 如 果包含了压力条件,就不需要调节了,参考压力位置也忽略了。参 考压力位置默认为单元的中心或者接近点(0,0,0)。有时候你可能 想要移动参考压力位置,也许要将它定位于绝对压力已知的点处 (比如:如果你想将计算结果和实验数据比较) 。要改变位置,

9、请 在操作压力面板中输入参考压力位置的新的坐标值(X,Y,Z) 。 菜5 全隐格式的优点是它关于时间步无条件稳定.使用分离或者耦合隐式解算器时。显式 时间步形式只能用在耦合显式解算器中。不可压流动。显式时间步 不能用于计算时间精度不可压流动(即:除了理想气体的气体定律)。在每一个时间步内,不可压解必须迭代直至收敛。 。显式形式主 要用于解决捕捉诸如激波之类的运动波的过渡问题。6 湍流模型选择湍流模型的要考虑的因素:流体是否可压、建立特殊的可 行的问题、精度的要求、计算机的能力、时间的限制。The Spalart-Allmaras 模型对于解决动力漩涡粘性,Spalart-Allmaras 模型

10、是相对简单的 方程。它包含了一组新的方程,在这些方程里不必要去计算和剪应 力层厚度相关的长度尺度。Spalart-Allmaras 模型是设计用于航空 领域的,主要是墙壁束缚流动,而且已经显示出和好的效果。在透 平机械中的应用也愈加广泛。 在原始形式中 Spalart-Allmaras 模型对于低雷诺数模型是十 分有效的,要求边界层中粘性影响的区域被适当的解决。在 FLUENT 中,Spalart-Allmaras 模型用在网格划分的不是很好时。这将是最 好的选择,当精确的计算在湍流中并不是十分需要时。再有,在模-_型中近壁的变量梯度比在 k-e 模型和 k- 模型中的要小的多。这 也许可以使

11、模型对于数值的误差变得不敏感。 需要注意的是 Spalart-Allmaras 模型是一种新出现的模型,现在 不能断定它适用于所有的复杂的工程流体。例如,不能依靠它去预 测均匀衰退,各向同性湍流。还有要注意的是,单方程的模型经常因 为对长度的不敏感而受到批评,例如当流动墙壁束缚变为自由剪切 流。k-e 模型标准标准 k-ek-e 模型模型最简单的完整湍流模型是两个方程的模型, 要解两个变量, 速度和长度尺度。 在 FLUENT 中,标准 k-e 模型自从被 Launder and Spalding 提出之后,就变成工程流场计算中主要的工具了。适 用范围广、经济、合理的精度,这就是为什么它在工业

12、流场和热交 换模拟中有如此广泛的应用了。它是个半经验的公式,是从实验现 象中总结出来的。由于人们已经知道了 k-e 模型适用的范围,因此 人们对它加以改造,出现了 RNG k-e 模型和带旋流修正 k-e 模型RNGRNG k-ek-e 模型模型RNG k-e 模型来源于严格的统计技术。它和标准 k-e 模型很相 似,但是有以下改进:RNG 模型在 e 方程中加了一个条件,有效 的改善了精度。 考虑到了湍流漩涡,提高了在这方面的精度。 RNG 理论为湍流 Prandtl 数提供了一个解析公式,然而标准 k-e 模型使用的是用户提供的常数。 然而标准 k-e 模型是一种高雷诺数的模型, RNG

13、理论提供 了一个考虑低雷诺数流动粘性的解析公式。这些公式的效用依靠正 确的对待近壁区域这些特点使得 RNG k-e 模型比标准 k-e 模型在更 广泛的流动中有更高的可信度和精度。带旋流修正的带旋流修正的 k-ek-e 模型模型带旋流修正的 k-e 模型是近期才出现的,比起标准 k-e 模型 来有两个主要的不同点。 带旋流修正的 k-e 模型为湍流粘性增加了一个公式。 为耗散率增加了新的传输方程, 这个方程来源于一个为层 流速度波动而作的精确方程术语“realizable” ,意味着模型要确保 在雷诺压力中要有数学约束,湍流的连续性。带旋流修正的 k-e 模 型直接的好处是对于平板和圆柱射流的

14、发散比率的更精确的预测。-_而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次 流有很好的表现。 带旋流修正的 k-e 模型和 RNG k-e 模型都显现出比标准 k-e 模型在强流线弯曲、漩涡和旋转有更好的表现。由于带旋流修正的 k-e 模型是新出现的模型,所以现在还没有确凿的证据表明它比 RNG k-e 模型有更好的表现。但是最初的研究表明带旋流修正的 k- e 模型在所有 k-e 模型中流动分离和复杂二次流有很好的作用。带 旋流修正的 k-e 模型的一个不足是在主要计算旋转和静态流动区域 时不能提供自然的湍流粘度。这是因为带旋流修正的 k-e 模型在定 义湍流粘度时考虑了平均旋度

15、的影响。这种额外的旋转影响已经在 单一旋转参考系中得到证实,而且表现要好于标准 k-e 模型。由于 这些修改,把它应用于多重参考系统中需要注意。k- 模型标准标准 k-k- 模型模型标准 k- 模型是基于 Wilcox k- 模型,它是为考虑低雷诺数、 可压缩性和剪切流传播而修改的。Wilcox k- 模型预测了自由剪 切流传播速率,像尾流、混合流动、平板绕流、圆柱绕流和放射状 喷射,因而可以应用于墙壁束缚流动和自由剪切流动。标准 k-e 模 型的一个变形是 SST k- 模型,它在 FLUENT 中也是可用的。剪切压力传输(剪切压力传输(SSTSST) k-k- 模型模型SST k- 模型由

16、 Menter 发展,以便使得在广泛的领域中可以 独立于 k-e 模型,使得在近壁自由流中 k- 模型有广泛的应用范围 和精度。为了达到此目的,k-e 模型变成了 k- 公式。SSTk- 模 型和标准 k- 模型相似,但有以下改进: SST k- 模型和 k-e 模型的变形增长于混合功能和双模型加 在一起。混合功能是为近壁 区域设计的,这个区域对标准 k- 模型有效,还有自由表面, 这对 k-e 模型的变形有效。 SST k- 模型合并了来源于 方程中的交叉扩散。 湍流粘度考虑到了湍流剪应力的传波。 模型常量不同 这些改进使得 SST k- 模型比标准 k- 模型在在广泛的流动 领域中有更高的

17、精度和可信度。雷诺压力模型(RSM)在 FLUENT 中 RSM 是最精细制作的模型。放弃等方性边界速度假-_设,RSM 使得雷诺平均 N-S 方程封闭,解决了关于方程中的雷诺压 力,还有耗散速率。这意味这在二维流动中加入了四个方程,而在 三维流动中加入了七个方程。由于 RSM 比单方程和双方程模型更加 严格的考虑了流线型弯曲、漩涡、旋转和张力快速变化, 它对于复 杂流动有更高的精度预测的潜力。 但是这种预测仅仅限于与雷诺压 力有关的方程。压力张力和耗散速率被认为是使 RSM 模型预测精度 降低的主要因素。RSM 模型并不总是因为比简单模型好而花费更多 的计算机资源。 但是要考虑雷诺压力的各向

18、异性时,必须用 RSM 模 型。例如飓风流动、燃烧室高速旋转流、管道中二次流。k-e 模型间的区别标准、 RNG 和带旋流修正 k-e 模型这三种模型有相似的形式, 有 k 方程和 e 方程,它们主要的不同点是: 计算湍流粘性的方法 湍流 Prandtl 数由 k 和 e 方程的湍流扩散决定 在 e 方程中湍流的产生和消失 每个模型计算湍流粘性的方法和模型的常数不一样。但从本质 上它们在其它方面是一样的。标准 k-e 模型是个半经验公式,主要 是基于湍流动能和扩散率。k 方程是个精确方程,e 方程是个由经验 公式导出的方程。k-e 模型假定流场完全是湍流,分之之间的粘性 可以忽略。标准 k-e

19、 模型因而只对完全是湍流的流场有效。 RNG 模型相比于标准 k-e 模型对瞬变流和流线弯曲的影响能作 出更好的反应, 这也可以解释 RNG 模型在某类流动中有很好的表现。这个模型(realize)对于和广泛的的流动有效, 包括旋转均 匀剪切流, 自由流中包括喷射和混合流,管道和边界流,还有分离 流。由于这些原因,这种模型比标准 k-e 模型要好。尤其需要注意 的是这种模型可以解决圆柱射流。 比如, 它预测了轴对称射流的 传播速率, 和平板射流一样。k- 模型间的区别俩种模型有相似的形式,有方程 k 和 。SST 和标准模型的不 同之处是 从边界层内部的标准 k- 模型到边界层外部的高雷诺数的

20、 ke 模型的逐渐转变 考虑到湍流剪应力的影响修改了湍流粘性公式标准 k- 模型 是一种经验模型,是基于湍流能量方程和扩散速率方程。由于 k- 模型已经修改多年,k 方程和 方程都增加了项,这样增加了模型 的精度。FLUENT 还提供了 SST 模型。它更适合对流减压区的计算。-_另外它还考虑了正交发散项从而使方程在近壁面和远壁面都适合7 非定常流动如果流动参数随时间而变,则这种流动成为非定常流动否则成 为定常流动(v p )8 分离和耦合离式求解器在计算不可压流体时, 不考虑压力做功和动能的影 响。在处理可压缩流动或使用任何耦合式求解器时,总是考虑压力 做功和动能。在使用分离式求解器时,FL

21、UENT 的缺省能量方程不包 含这一项(粘性耗散项) (因为一般灰忽略粘性产生的热量) 。当 Brinkman 数,Br,接近或大于 1 时,流体粘性生成的热量不可忽 视:当用户的问题需要考虑粘性耗散项,并且使用的是分离式求解 器时,用户应该在 Viscous Model panel 面板中使用 Viscous Heating 选项激活此项。对于一般的可压缩流动,需要注意的是, 如果用户已经定义了可压缩流动,但使用了分离式求解器,FLUENT 不会自动激活粘性耗散选项。对于任一种耦合式求解器,在求解能 量方程时,粘性耗散项总是被考虑进去。 考虑组分扩散项:当使用分离式求解器时,FLUET 允许

22、用户对 固体介质设定各向异性的导热率。如果使用分离式求解器,用户可 以禁止掉入口的能量扩散来达到此目的。如果使用耦合式求解器, 那么,入口区的扩散传热不能禁止掉。如果模拟的是粘性流动,并 且希望在能量方程中包含粘性生成热,请在 Viscous Model panel. 面板中激活 Viscous Heating 选项 (可选的, 且仅适用于离散求 解器, ) 。在使用离散求解器时,缺省情况下,FLUENT 在能量方程中 忽略了粘性生成热(各个耦合求解器总是包含有粘性生成热) 。对 于流体剪切应力较大(例如流体润滑问题) ,和或高速可压缩流 动(请参阅方程 11.2-9) ,用户应该考虑粘性耗散

23、。 传热计算的求解过程中,能量方程的亚松弛当用户使用分离式 求解器时,FLUENT 使用亚松弛参数(系数,在 Solution Controls panel 面板中由用户定义,请参阅 22.9 节)来对能量方程进行亚 松弛处理。如果使用分离式求解器时收敛困难, 那么, 用户可以 考虑在 Species Model panel.面板中禁止 Diffusion Energy Source 选项。需要注意的是,对于任何耦合求解器,总是包含有组 分扩散的影响。 对于传热计算,更有效的计算策略是先计算绝热流动,然后再 考虑能量方程的计算。依据流动与传热是否耦合,求解过程稍有不 同。若用户的问题是非耦合的

24、流动-传热过程(没有依赖于温度的介 质属性或浮力) ,那么,用户可以先求解绝热流动 (屏蔽掉能量-_方程) 以得到收敛的流场, 然后再单独求解能量输运方程。 由于耦合求解器同时求解流动与能量方程,所以,上述的能量 方程单独求解过程仅对分离式求解器有效。流动与传热的耦合计算 对于流动与传热耦合问题(例如,模型中包含有依赖于温度的介质 属性或浮力) ,在计算能量方程之前,用户可以首先求解流动方程。 获得收敛的流场计算结果之后,用户可以再选择能量方程,然后同 时求解流动与传热方程,最终获得问题的完整解。在 FLUENT 的非 耦合求解器中,入口的物质净输送量由对流量和扩散量组成,对耦 合解算器,只包

25、括对流部分。对流部分由你指定的物质浓度确定。 扩散部分依赖于计算得到的物质浓度场。因此,扩散部分(从而使 净输送量)不预先指定。9 多相流模型VOF 模型所谓 VOF 模型是一种在固定的欧拉网格下的表面跟踪方法。当 需要得到一种或多种互不相融流体间的交界面时,可以采用这种模 型。在 VOF 模型中, 不同的流体组分共用着一套动量方程,计算 时在全流场的每个计算单元内,都记录下各流体组分所占有的体积 率。VOF 模型的应用例子包括分层流,自由面流动,灌注,晃动, 液体中大气泡的流动,水坝决堤时的水流,对喷射衰竭(jet breakup)(表面张力)的预测,以及求得任意液气分界面的稳态或 瞬时分界

26、面混合物模型混和物模型可用于两相流或多相流(流体或颗粒) 。因为在欧拉 模型中, 各相被处理为互相贯通的连续体,混和物模型求解的是混 合物的动量方程, 并通过相对速度来描述离散相。混合物模型的应 用包括低负载的粒子负载流,气泡流,沉降,以及旋风分离器。混 合物模型也可用于没有离散相相对速度的均匀多相流。欧拉模型欧拉模型是 Fluent 中最复杂的多相流模型。它建立了一套包 含有 n 个的动量方程和连续方程来求解每一相。 压力项和各界面 交换系数是耦合在一起的。耦合的方式则依赖于所含相的情况,颗 粒流(流固)的处理与非颗粒流(流-流)是不同的。对于颗粒流, 可应用分子运动理论来求得流动特性。不同

27、相之间的动量交换也依 赖于混合物的类别。通过 FLUENT 的客户自定义函数(user- defined functions) ,你可以自己定义动量交换的计算方式。欧拉-_模型的应用包括气泡柱,上浮,颗粒悬浮,以及流化床。 对于体积率小于 10的气泡、液滴和粒子负载流动,采用离散 相模型。 对于离散相混合物或者单独的离散相体积率超出 10的气泡、 液滴和粒子负载流动,采用混合物模型。或者欧拉模型。 对于活塞流,采用 VOF 模型。 对于分层/自由面流动,采用 VOF 模型。 对于气动输运, 如果是均匀流动, 则采用混合物模型; 如果 是粒子流见,则采用欧拉模型。 对于流化床,采用欧拉模型模拟粒

28、子流。 对于泥浆流和水力输运,采用混合物模型或欧拉模型 对于沉降,采用欧拉模型。 对于更加一般的,同时包含若干种多相流模式的情况,应根据 最感兴趣的流动特征,选择合适的流动模型。此时由于模型只是对 部分流动特征做了较好模拟,其精度必然低于只包含单个模式的流 动。 对于分层流和活塞流,最直接的就是选择 VOF 模型,离散相模 型只适用于低体积率的情况。但同时,也只有离散相模型才允许你 指定颗粒的分布或者在多相流模型中同时加入燃烧模型。9通用多相流模型VOF 模型适合于分层的或自由表面流,而 mixture 和 Eulerian 模型适合于流动中有相混合或分离,或者分散相的 volume frac

29、tion 超过 10%的情形。 (流动中分散相的 volume fraction 小 于或等于 10%时可使用第 19 章讨论过的离散相模型) 。为了在 mixture 模型和 Eulerian 模型之间作出选择,除了 Section18.4 中详细的指导外, 你还应考虑以下几点: 如果分散相有着宽广的分布,mixture 模型是最可取的。如 果分散相只集中在区域的一部分,你应当使用 Eulerian 模型。 如果应用于你的系统的相间曳力规律是可利用的(either within FLUENT or through a user-defined function) ,Eulerian 模型通常

30、比 mixture 模型能给出更精确的结果。如果 相间的曳力规律不知道或者它们应用于你的系统是有疑问的, mixture 模型可能是更好的选择。 如果你想解一个需要计算付出较少的简单的问题,mixture 模型可能是更好的选择,因为它比 Eulerian 模型要少解一部分方 程。如果精度比计算付出更重要,Eulerian 模型是更好的选择。但 是请记住,复杂的 Eulerian 模型比 mixture 模型的计算稳定性要-_差。10 松弛因子为了改进耦合计算的收敛性,用户可减小松弛系数。 由于流体力学中要求解非线性的方程,在求解过程中,控制变量 的变化是很必要的,这就通过松弛因子来实现的.它控

31、制变量在每次 迭代中的变化.也就是说,变量的新值为原值加上变化量乘以松弛因 子. 如:A1=A0+B*DETA A1 新值 A0 原值 B 松弛因子 DETA 变化量 松弛因子可控制收敛的速度和改善收敛的状况! 为 1,相当于不用松弛因子 大于 1,为超松弛因子,加快收敛速 度 小于 1,欠松弛因子,改善收敛的条件 一般来讲,大家都是在收敛 不好的时候, 采用一个较小的欠松弛因子。 Fluent 里面用的是欠松弛,主要防止两次迭代值相差太大引起 发散。 松弛因子的值在 01 之间,越小表示两次迭代值之间变化 越小,也就越稳 定,但收敛也就越慢。a a 亚孙持因子亚孙持因子 1、亚松驰(Unde

32、r Relaxation):所谓亚松驰就是将本层次计 算结果与上 一层次结果的差值作适当缩减, 以避免由于差值过大 而引起非线性迭代过程的发 散。用通用变量来写出时,为松驰因子 (Relaxation Factors) 。 数值传热学 -2142、FLUENT 中的亚松驰:由于 FLUENT 所解方程组的非线性, 我们有必要控制的变化。一般用亚松驰方法来实现控制,该方法在 每一部迭代中减少了的变 化量。亚松驰最简单的形式为:单元内变 量 等于原来的值 加上亚松驰因子 a 与 变化的积分离解算器使用 亚松驰来控制每一步迭代中的计算变量的更新。这 就意味着使用分 离解算器解的方程,包括耦合解算器所

33、解的非耦合方程(湍流和 其 他标量)都会有一个相关的亚松驰因子。注注:在 FLUENT 中,所有变量的默认亚松驰因子都是对大多数 问题的最优 值。这个值适合于很多问题,但是对于一些特殊的非线 性问题(如:某些湍流或 者高 Rayleigh 数自然对流问题) ,在计 算开始时要慎重减小亚松驰因子。 使用默认的亚松驰因子开始计算 是很好的习惯。如果经过 4 到 5 步的迭代 残差仍然增长,你就需 要减小亚松驰因子。 有时候,如果发现残差开始增加,你可以改变 亚松驰因子重新计算。在亚松 驰因子过大时通常会出现这种情况。 最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文 件,并对解的算法做几

34、步迭代以调节到新的参数。最典型 的情况是, 亚松驰因子的增加会使残差有少量的增加,但是随着解的进行残差 的 增加又消失了。 如果残差变化有几个量级你就需要考虑停止计-_算并回到最后保存 的较好的数据文件。 注意:注意:粘性和密度的亚松驰是在每一次迭代之间的。而且,如 果直接解焓方 程而不是温度方程(即:对 PDF 计算) ,基于焓的温 度的更新是要进行亚松驰 的。要查看默认的亚松弛因子的值,你可 以在解控制面板点击默认按钮。 对于大多数流动,不需要修改默认 亚松弛因子。但是,如果出现不稳定或者 发散你就需要减小默认的 亚松弛因子了,其中压力、动量、k 和 e 的亚松弛因子 默认值分 别为 0.

35、2,0.5,0.5 和 0.5。 对于 SIMPLEC 格式一般不需要减小 压力的亚松弛因子。 在密度和温度强烈 耦合的问题中,如相当高 的 Rayleigh 数的自然或混合对流流动,应该对温度和/ 或密度 (所用的亚松弛因子小于 1.0)进行亚松弛。相反,当温度和动量 方程没 有耦合或者耦合较弱时,流动密度是常数, 温度的亚松弛 因子可以设为 1.0。 对于其它的标量方程,如漩涡,组分, PDF 变量,对于某些问题默认的亚松弛可能过大,尤其是对于初始计算。 你可 以将松弛因子设为 0.8 以使得收敛更容易。 11 库朗数Courant number 实际上是指时间步长和空间步长的相对关系,系

36、统 自动减小 courant 数,这种情况一般出现在存在尖锐外形的计算域, 当局部的流速过大或者压差过大时出错,把局部的网格加密再试一 下。 在 FLUENT 中,用 courant number 来调节计算的稳定性与收敛性。 一般来说,随着 courant number 的从小到大的变化,收敛速度逐渐 加快,但是稳定性逐渐降低。所以具体的问题,在计算的过程中, 最好是把 courant number 从小开始设置,看看迭代残差的收敛情况, 如果收敛速度较慢而且比较稳定的话,可以适当的增加 courant number 的大小,根据自己具体的问题,找出一个比较合适的 courant number,让收敛速度能够足够的快,而且能够保持它的稳定性。 FLUENT 计算开始迭代最好使用较小的库朗数,否则容易导致迭代 发散,修改办法 slovecontrolssolution,修改 courant Number 默 认值为 1,开始没有经验的改小点,比如 0.01,然后逐渐加大。接 9通用多相流-_-_-_

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 教案示例

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知得利文库网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号-8 |  经营许可证:黑B2-20190332号 |   黑公网安备:91230400333293403D

© 2020-2023 www.deliwenku.com 得利文库. All Rights Reserved 黑龙江转换宝科技有限公司 

黑龙江省互联网违法和不良信息举报
举报电话:0468-3380021 邮箱:hgswwxb@163.com