搜档网
当前位置:搜档网 › CFD 计算对计算网格有特殊的要求

CFD 计算对计算网格有特殊的要求

CFD 计算对计算网格有特殊的要求
CFD 计算对计算网格有特殊的要求

CFD 计算对计算网格有特殊的要求,一是考虑到近壁粘性效应采用较密的贴体网格,二是网格的疏密程度与流场参数的变化梯度大体一致。对于面网格,可以设置平行于给定边的边界层网格,可以指定第二层与第一层的间距比,及总的层数。对于体网格,也可以设置垂直于壁面方向的边界层,从而可以划分出高质量的贴体网格。而其它通用的CAE 前处理器主要是根据结构强度分析的需要而设计的,在结构分析中不存在边界层问题,因而采用这种工具生成的网格难以满足CFD 计算要求,而Gambit 软件解决了这个特殊要求。

如果先在一条边上画密网格再在之上画边界层,边界层与网格能很好的对应起来

如果直接在一条边上画边界层,则边界层横向之间的距离很宽

怎么设置边界层横向之间的距离,即不用先画网格也能画出横向距离很密的边界层来?

在划分边界层网格之前,用粘性网格间距计算器,计算出想要的y+值对应的第一层网格高度;

第一层高度出来之后,关于网格的纵横向网格间距之比,也就是边界层第一层网格高度与横向间距之比,大概在1/sqrt(Re),最为适宜;

先在你要划边界层网格的边上划分线网格,然后再划分边界层。

gambit本人也用了一段时间,六面体网格四面体网格我都画过,但是最头疼的还是三维边界层网格的生成。

用gambit自带的边界层网格生成功能画出来的边界层网格经常达不到好的效果,或者对于复杂的外形根本就无法生成边界层网格。为此我就采用手动设置边界层,但是比较费时间,效果还一般。

不知道大家是不是也遇到相似的问题,或者有更好的方法,请指点一下,先谢谢了!

22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?

1、亚松驰(Under Relaxation):所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减,以避免由于差值过大而引起非线性迭代过程的发散。用通用变量来写出时,为松驰因子(Relaxati on Factors)。《数值传热学-214》

2、FLUENT中的亚松驰:由于FLUENT所解方程组的非线性,我们有必要控制的变化。一般用亚松驰方法来实现控制,该方法在每一部迭代中减少了的变化量。亚松驰最简单的形式为:单元内变量等于原来的值加上亚松驰因子a与变化的积, 分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程,包括耦合解算器所解的非耦合方程(湍流和其他标量)都会有一个相关的亚松驰因子。在FLUENT中,所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题,但是对

于一些特殊的非线性问题(如:某些湍流或者高Rayleigh数自然对流问题),在计算开始时要慎重减小亚松驰因子。使用默认的亚松驰因子开始计算是很好的习惯。如果经过4到5步的迭代残差仍然增长,你就需要减小亚松驰因子。有时候,如果发现残差开始增加,你可以改变亚松驰因子重新计算。在亚松驰因子过大时通常会出现这种情况。最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文件,并对解的算法做几步迭代以调节到新的参数。最典型的情况是,亚松驰因子的增加会使残差有少量的增加,但是随着解的进行残差的增加又消失了。如果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好的数据文件。注意:粘性和密度的亚松驰是在每一次迭代之间的。而且,如果直接解焓方程而不是温度方程(即:对P DF计算),基于焓的温度的更新是要进行亚松驰的。要查看默认的亚松弛因子的值,你可以在解控制面板点击默认按钮。对于大多数流动,不需要修改默认亚松弛因子。但是,如果出现不稳定或者发散你就需要减小默认的亚松弛因子了,其中压力、动量、k和e的亚松弛因子默认值分别为0.2,0.5,0.5和0.5。对于SIMPLEC格式一般不需要减小压力的亚松弛因子。在密度和温度强烈耦合的问题中,如相当高的Rayleigh数的自然或混合对流流动,应该对温度和/或密度(所用的亚松弛因子小于1.0)进行亚松弛。相反,当温度和动量方程没有耦合或者耦合较弱时,流动密度是常数,温度的亚松弛因子可以设为1.0。对于其它的标量方程,如漩涡,组分,PDF变量,对于某些

问题默认的亚松弛可能过大,尤其是对于初始计算。你可以将松弛因子设为0.8以使得收敛更容易。

SIMPLE与SIMPLEC比较

在FLUENT中,可以使用标准SIMPLE算法和SIMPLEC(SIMP LE-Consistent)算法,默认是SIMPLE算法,但是对于许多问题如果使用SIMPLEC可能会得到更好的结果,尤其是可以应用增加的亚松驰迭代时,具体介绍如下:

对于相对简单的问题(如:没有附加模型激活的层流流动),其收敛性已经被压力速度耦合所限制,你通常可以用SIMPLEC算法很快得到收敛解。在SIMPLEC中,压力校正亚松驰因子通常设为1.0,它有助于收敛。但是,在有些问题中,将压力校正松弛因子增加到1.0可能会导致不稳定。对于所有的过渡流动计算,强烈推荐使用PISO 算法邻近校正。它允许你使用大的时间步,而且对于动量和压力都可以使用亚松驰因子1.0。对于定常状态问题,具有邻近校正的PISO

并不会比具有较好的亚松驰因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。当你使用PISO邻近校正时,对所有方程都推荐使用亚松驰因子为1. 0或者接近1.0。如果你只对高度扭曲的网格使用PISO倾斜校正,请设定动量和压力的亚松驰因子之和为1.0比如:压力亚松驰因子0.3,动量亚松驰因子0.7)。如果你同时使用PISO的两种校正方法,推荐参阅PISO邻近校正中所用的方法。

实体、实面与虚体、虚面的区别

在建模中,经常会遇到实...与虚...,而且虚体的计算域好像也可以进行计算并得到所需的结果,对二者的根本区别及在功能上的不同对于求解是没有任何区别的,只要你能在虚体或者实体上划分你需要的网格Gambit的实体和虚体在生成网格和计算的时候对于结果没有任何影响,实体和虚体的主要区别有以下几点:

1.实体可以进行布尔运算但是虚体不能,虽然不能进行布尔运算,但是虚体存在merge,split等功能;

2.实体运算在很多cad软件里面都有,但是虚体是gambit的一大特色,有了虚体以后,Gambit的建模和网格生成的灵活性增加了很多。

3.在网格生成的过程中,如果有几个相对比较平坦的面,你可以把它们通过merge合成一个,这样,作网格的时候,可以节省步骤,对于曲率比较大的面,可能生成的网格质量不好,这时候,你可以采取用split的方式把它划分成几个小面以提高网格质量。

对于虚体生成的计算网格,和实体生成的计算网格,在计算的时候没有区别,关键是看网格生成的质量如何,与实体虚体无关。经常在作复杂模型计算的时候,大部分都是用的虚体,特别是从其他的建模软件里面导进来的复杂模型,基本上不能够生成实体。至于计算的效果如何,与Fluent的设置和网格的质量有关,与模型无关。

21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在F LUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,

并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?

可以采用残差控制面板来显示;或者采用通过某面的流量控制;如监控出口上流量的变化;采用某点或者面上受力的监视;涡街中计算达到收敛时,绕流体的面上受的升力为周期交变,而阻力为平缓的直线。

怎样判断计算结果是否收敛?

1、观察点处的值不再随计算步骤的增加而变化;

2、各个参数的残差随计算步数的增加而降低,最后趋于平缓;

3、要满足质量守恒(计算中不牵涉到能量)或者是质量与能量守恒(计算中牵涉到能量)。

特别要指出的是,即使前两个判据都已经满足了,也并不表示已经得到合理的收敛解了,因为,如果松弛因子设置得太紧,各参数在每步计算的变化都不是太大,也会使前两个判据得到满足。此时就要再看第三个判据了。

还需要说明的就是,一般我们都希望在收敛的情况下,残差越小越好,但是残差曲线是全场求平均的结果,有时其大小并不一定代表计算结果的好坏,有时即使计算的残差很大,但结果也许是好的,关键是要看计算结果是否符合物理事实,即残差的大小与模拟的物理现象本身的复杂性有关,必须从实际物理现象上看计算结果。比如说一个全机模型,在大攻角情况下,解震荡得非常厉害,而且残差的量级也

总下不去,但这仍然是正确的,为什么呢,因为大攻角下实际流动情形就是这样的,不断有涡的周期性脱落,流场本身就是非定常的,所以解也是波动的,处理的时候取平均就可以呢:)

26 什么叫问题的初始化?在FLUENT中初始化的方法对计算结果有什么样的影响?初始化中的“patch”怎么理解?

问题的初始化就是在做计算时,给流场一个初始值,包括压力、速度、温度和湍流系数等。理论上,给的初始场对最终结果不会产生影响,因为随着跌倒步数的增加,计算得到的流场会向真实的流场无限逼近,但是,由于Fluent等计算软件存在像离散格式精度(会产生离散误差)和截断误差等问题的限制,如果初始场给的过于偏离实际物理场,就会出现计算很难收敛,甚至是刚开始计算就发散的问题。因此,在初始化时,初值还是应该给的尽量符合实际物理现象。这就要求我们对要计算的物理场,有一个比较清楚的理解。

初始化中的patch就是对初始化的一种补充,比如当遇到多相流问题时,需要对各相的参数进行更细的限制,以最大限度接近现实物理场。这些就可以通过patch来实现,patch可以对流场分区进行初始化,还可以通过编写简单的函数来对特定区域初始化。

以上是我对这个问题的理解,欢迎大家批评指正。

94 把带网格的几个volume,copy到另一处,但原来split的界面,现在都变成了wall,怎么才能把wall变成内部流体呢?

===〉直接边界面定义为interior即可

102在FLUENT模拟以后用display下的操作都无法显示,不过刚开始用的是好的,然后就不行了,为什么?

===> DirectX 控制面板中的“加速”功能禁用即可

128 在gambit中对一体积成功的进行了体网格,网格进行了exa mine mesh,也没有什么问题,可当要进行边界类型(boundary type)的设定时,却发现type 只有node,element_side两项,没有什么w all,pressure_outlet等。为何无法定义边界?

===> 因为没有选择求解器为fluent 5/6

43 FLUENT中常用的文件格式类型:dbs,msh,cas,dat,trn,j ou,profile等有什么用处?

在Gambit目录中,有三个文件,分别是default_id.dbs,jou,trn 文件,对Gambit运行save,将会在工作目录下保存这三个文件:def ault_id.dbs,default_id.jou,default_id.trn。

jou文件是gambit命令记录文件,可以通过运行jou文件来批处理gambit命令;

dbs文件是gambit默认的储存几何体和网格数据的文件;

trn文件是记录gambit命令显示窗(transcript)信息的文件;

msh文件可以在gambit划分网格和设置好边界条件之后export中选择msh文件输出格式,该文件可以被fluent求解器读取。

Case文件包括网格,边界条件,解的参数,用户界面和图形环境。

Data文件包含每个网格单元的流动值以及收敛的历史纪录(残差值)。Fluent自动保存文件类型,默认为date和case文件Profile文件边界轮廓用于指定求解域的边界区域的流动条件。例如,它们可以用于指定入口平面的速度场。

读入轮廓文件,点击菜单File/Read/Profile...弹出选择文件对话框,你就可以读入边界轮廓文件了。

写入轮廓文件,你也可以在指定边界或者表面的条件上创建轮廓文件。例如:你可以在一个算例的出口条件中创建一个轮廓文件,然后在其它算例中读入该轮廓文件,并使用出口轮廓作为新算例的入口轮廓。要写一个轮廓文件,你需要使用Write Profile面板(Figure 1),菜单:File/Write/Profile...

第23题:

Let's take care of the warning "turbulent viscosity limited to visc osity ratio****" which is not physical. This problem is mainly due to one of the following:

1)Poor mesh quality(i.e.,skewness > 0.85 for Quad/Hex, or skew ness > 0.9 for Tri/Tetra elements). {what values do you have?}

2)Use of improper turbulent boudary conditions.

3)Not supplying good initial values for turbulent quantities.

出现这个警告,一般来讲,最可能的就是网格质量的问题,尤其是Y+值的问题;在划分网格的时候要注意,第一层网格高度非常重要,可以使用NASA的Viscous Grid Space Calculator来计算第一层网格高度;如果这方面已经注意了,那就可能是边界条件中有关湍流量的设置问题,关于这个,本版中已经有专门的帖子进行了讨论,Fl uent培训的教程中也有讲到,请大家参考。

第24题:

这个问题的意思是出现了回流,这个问题相对于湍流粘性比的警告要宽松一些,有些case可能只在计算的开始阶段出现这个警告,随着迭代的计算,可能会消失,如果计算一段时间之后,警告消失了,那么对计算结果是没有什么影响的,如果这个警告一直存在,可能需要作以下处理:

1.如果是模拟外部绕流,出现这个警告的原因可能是边界条件取得距离物体不够远,如果边界条件取的足够远,该处可能在计算的过程中的确存在回流现象;对于可压缩流动,边界最好取在10倍的物体特征长度之处;对于不可压缩流动,边界最好取在4倍的物体特征长度之处。

2.如果出现了这个警告,不论对于外部绕流还是内部流动,可以使用pressure-outlet边界条件代替outflow边界条件改善这个问题。

第30题:

这个问题我也一直在想,看到simwe上有人这么回答的:(顺便把相类似的问题的解答也放在一起,方便大家一起解决这类的问题。)

一. 残差波动的主要原因:1、高精度格式;2、网格太粗;3、网格质量差;4、流场本身边界复杂,流动复杂;5、模型的不恰当使用。

二. 问:在进行稳态计算时候,开始残差线是一直下降的,可是到后来各种残差线都显示为波形波动,是不是不收敛阿?

答:有些复杂或流动环境恶劣情形下确实很难收敛。计算的精度(2阶),网格太疏,网格质量太差,等都会使残差波动。经常遇到,一开始下降,然后出现波动,可以降低松弛系数,我的问题就能收敛,但如果网格质量不好,是很难的。通常,计算非结构网格,如果问题比较复杂,会出现这种情况,建议作网格时多下些功夫。理论上说,残差的震荡是数值迭代在计算域内传递遭遇障碍物反射形成周期震荡导致的结果,与网格亚尺度雷诺数有关。例如,通常压力边界是主要的反射源,换成OUTFLOW 边界会好些。这主要根据经验判断。所以我说网格和边界条件是主要因素。

三. 1、网格问题:比如流场内部存在尖点等突变,导致网格在局部质量存在问题,影响收敛。

2、可以调整一下courant number,courant number实际上是指时

间步长和空间步长的相对关系,系统自动减小courant数,这种情况一般出现在存在尖锐外形的计算域,当局部的流速过大或者压差过大时出错,把局部的网格加密再试一下。

在fluent中,用courant number来调节计算的稳定性与收敛性。一般来说,随着courant number的从小到大的变化,收敛速度逐渐加快,但是稳定性逐渐降低。所以具体的问题,在计算的过程中,最好是把courant number从小开始设置,看看迭代残差的收敛情况,如果收敛速度较慢而且比较稳定的话,可以适当的增加courant number的大小,根据自己具体的问题,找出一个比较合适的courant number,让收敛速度能够足够的快,而且能够保持它的稳定性。

四. 另外,Fluent官方论坛上有一个管理员这么解释的:

1. Residual pattern in an unsteady simulation will be oscillatory. At each time-step, Fluent will try to converge residual values belo w the default limits set in the Residual Monitors panel.

2. In a steady simulation, ideally the residuals should not oscilla te.

Oscillating residuals may be due to wrong case setup. Following co mments should help in reducing convergence trouble.

a. Make sure that grid is of good quality. Please refer the sectio n 6.2.2 of Fluent 6.3 Users Guide for more details on the grid qua

lity. https://www.sodocs.net/doc/f916629400.html,/fluent/doc/ori/html/ug/node155.htm

b. Recheck all the boundary condition.

c. If you are using turbulence model, please check the turbulence p arameters on all inlet / outlet boundaries. Please refer to the follow ing white paper for more details regarding the different turbulent bo undary conditions:

https://www.sodocs.net/doc/f916629400.html,/software/university/whitepapers/turbulent.pdf

d. Try with better initial guess.

3. The detail about Grid Adaption is available in Chapter 26 of Fluent Users Guide:

https://www.sodocs.net/doc/f916629400.html,/fluent6326/doc/ori/html/ug/node1071.htm 本人觉得可以重点参考第四个回答。另外,如果出现连续方程残差很高收敛慢的情况,首先应该检查的是网格质量;由于现在大量使用分块网格,这时要看看两相邻块处的网格大小是不是相差较大,也就是看看有没有出现cell jump的情况,相邻网格的大小最好不要超过2倍的关系,这时出现高连续方程残差的一个主要原因,这需要在划分网格时做好规划。

第31题:

假扩散(false diffusion)的含义:

基本含义:由于对流—扩散方程中一阶导数项的离散格式的截断误差小于二阶而引起较大数值计算误差的现象。有的文献中将人工粘性(artificial viscosity)或数值粘性(numerical viscosity)视为它的同义词。

拓宽含义:现在通常把以下三种原因引起的数值计算误差都归在假扩散的名称下

1.非稳态项或对流项采用一阶截差的格式;

2.流动方向与网格线呈倾斜交叉(多维问题);

3.建立差分格式时没有考虑到非常数的源项的影响。

克服或减轻假扩散的格式或方法,

为克服或减轻数值计算中的假扩散(包括流向扩散及交叉扩散)误差,应当:

1. 采用截差阶数较高的格式;

2. 减轻流线与网格线之间的倾斜交叉现象或在构造格式时考虑到来流方向的影响。

3. 至于非常数源项的问题,目前文献中,还没有为克服这种影响而专门构造的格式,但是高阶格式显然对减轻其影响是有利的。

第34题:

GAUGE PRESSURE 就是静压。

GAUGE total PRESSURE 是总压。

这里需要强调一下Gauge为名义值,

什么意思呢?如果,INITIAL Gauge PRESSURE =0

那么GAUGE PRESSURE 就是实际的静压Pinf。

GAUGE total PRESSURE 是实际的总压Pt。

如果INITIAL Gauge PRESSURE 不等于零

GAUGE PRESSURE =Pinf - INITIAL Gauge PRESSURE GAUGE total PRESSURE =Pt - INITIAL Gauge PRESSURE 第132题:

大概有这样一个估计:“1k网格=1M内存。”

对于一台有1G内存的计算机,你能接受的计算网格数最好少于1 00万,当然这只是一个粗略的说法,影响计算速度的因素还有Fluen t计算的设置等。

第11题:

关于贴体网格的定义,随便找本CFD的书上面基本上都有介绍,在此就不多言。在此谈谈网格独立性的问题。

数值计算的与实验值之间的误差来源只要有这几个:物理模型近似误差(无粘或有粘,定常与非定常,二维或三维等等),差分方程的截断误差及求解区域的离散误差(这两种误差通常统称为离散误差),迭代误差(离散后的代数方程组的求解方法以及迭代次数所产生的误差),舍入误差(计算机只能用有限位存储计算物理量所产生的误差)等等。在通常的计算中,离散误差随网格变细而减小,但由于网格变细时,离散点数增多,舍入误差也随之加大。

由此可见,网格数量并不是越多越好的。

再说说网格无关性的问题,由上面的介绍,我们知道网格数太密或者太疏都可能产生误差过大的计算结果,网格数在一定的范围内的结果才与实验值比较接近,这样在划分网格时就要求我们首先依据已有的经验大致划分一个网格进行计算,将计算结果(当然这个计算结果必须是收敛的)与实验值进行比较(如果没有实验值,则不需要比较,后面的比较与此类型相同),再酌情加密或减少网格,再进行计算,再与实验值进行比较,并与前一次计算结果比较,如果两次的计算结果相差较小(例如在2%),说明这一范围的网格的计算结果是可信的,说明计算结果是网格无关的。再加密网格已经没有什么意义

(除非你要求的计算精度较高)。但是,如果你用粗网格也能得到相差很小的计算结果,从计算效率上讲,你就可以完全使用粗网格去完成你的计算。加密或者减少网格数量,你可以以一倍的量级进行。

第154题:

在网上看到很多关于Y+的估算公式,Fluent公司的培训教程上也有相关的估算方法;在使用这些公式在给定的Y+要求下反过来计算第一层网格高度时,如果估算的不准,划分出来的网格得到的Y+值可能就相差的有点儿大。

关于Y+的计算公式,在这里也就不多啰嗦了,请参考Fluent的培训教程,关于湍流模型的章节。

在这里,建议大家不要自己动手估算,最好是使用NASA Viscou s Grid Space Calculator计算,网址为:

https://www.sodocs.net/doc/f916629400.html,/APPS/YPlus/

另外,网上还有另一个计算器:

网址为:https://www.sodocs.net/doc/f916629400.html,/cfd-shop/goodies/yplus.htm

通过比较这两个计算器的源代码,个人感觉NASA的计算器适用范围更广,考虑的因素更多,因此,建议大家优先使用NASA的计算器。

第49题:

从字面的意思很好理解axisymmetric和axisymmetric swirl的差别:

axisymmetric:是轴对称的意思,也就是关于一个坐标轴对称,2 D的axisymmetric问题仍为2D问题。

而axisymmetric swirl:是轴对称旋转的意思,就是一个区域关于一条坐标轴回转所产生的区域,这产生的将是一个回转体,是3D的问题。在Fluent中使用这个,是将一个3D的问题简化为2D问题,减少计算量,需要注意的是,在Fluent中,回转轴必须是x轴。

第13题:

具体内容可以参考Gambit Documentation中的Quality Type Defi nitions章节。

判断网格质量的方面有:

Area单元面积,适用于2D单元,较为基本的单元质量特征。

Aspect Ratio长宽比,不同的网格单元有不同的计算方法,等于1是最好的单元,如正三角形,正四边形,正四面体,正六面体等;一般情况下不要超过5:1.

Diagonal Ratio对角线之比,仅适用于四边形和六面体单元,默认是大于或等于1的,该值越高,说明单元越不规则,最好等于1,也就是正四边形或正六面体。

Edge Ratio长边与最短边长度之比,大于或等于1,最好等于1,解释同上。

EquiAngle Skew通过单元夹角计算的歪斜度,在0到1之间,0为质量最好,1为质量最差。最好是要控制在0到0.4之间。

EquiSize Skew通过单元大小计算的歪斜度,在0到1之间,0为质量最好,1为质量最差。2D质量好的单元该值最好在0.1以内,3 D单元在0.4以内。

MidAngle Skew通过单元边中点连线夹角计算的歪斜度,仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。

Size Change相邻单元大小之比,仅适用于3D单元,最好控制在2以内。

Stretch伸展度。通过单元的对角线长度与边长计算出来的,仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。

Taper锥度。仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。

Volume单元体积,仅适用于3D单元,划分网格时应避免出现负体积。

Warpage翘曲。仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。

以上只是针对Gambit帮助文件的简单归纳,不同的软件有不同的评价单元质量的指标,使用时最好仔细阅读帮助文件。

另外,在Fluent中的窗口键入:grid quality 然后回车,Fluent能检查网格的质量,主要有以下三个指标:

1.Maxium cell squish: 如果该值等于1,表示得到了很坏的单元;

2.Maxium cell skewness: 该值在0到1之间,0表示最好,1表示最坏;

3.Maxium 'aspect-ratio': 1表示最好。

以上的一些只是简略提要,具体的请参考相关资料。

第159题:在fluent中如何设置工作目录?在Gambit中如何设置工作目录?

找到桌面上的Fluent或者Gambit图标,右键图标,“属性”->"起始位置”...

将起始位置设置为你想要的文件夹目录就可以了。

这种设置对于Fluent有效,但有时对Gambit无效,不知道是什么原因。

第16题:在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?

有限元网格划分的基本原则

有限元网格划分的基本原则 划分网格是建立有限元模型的一个重要环节,它要求考虑的问题较多,需要的工作量较大,所划分的网格形式对计算精度和计算规模将产生直接影响。为建立正确、合理的有限元模型,这里介绍划分网格时应考虑的一些基本原则。 1 网格数量 网格数量的多少将影响计算结果的精度和计算规模的大小。一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因数综合考虑。图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。所以应注意增加网格的经济性。实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。 图1 位移精度和计算时间随网格数量的变化 在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格,如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。 2 网格疏密 网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。这样,整个结构便表现出疏密不同的网格划分形式。图2是中心带圆孔方板的四分之一模型,其网格反映了疏密不同的划分原则。小圆孔附近存在应力集中,采用了比较密的网格。板的四周应力梯度较小,网格分得较稀。其中图b中网格疏密相差更大,它比图a中的网格少48个,但计算出的孔缘最大应力相差1%,而计算时间却减小了36%。由此可见,采用疏密不同的网格划分,既可以保持相当的计算精度,又可使网格数量减

ICEM-CFD-网格划分入门

WorkBe nch ICEM CFD网格划分入门111A nsysWB!集成了一个非常重要的工具:ICEM CFD 它是一个建模、划分网格的集成工具,功能非常强大。我也只是蜻蜓点水的用了几次,感觉确实非常棒,以前遇到复杂的模型,用过几个划分网格的工具。但这是我觉得最方便和最具效率的。 网格划分很大程度上影响着后续的仿真分析一一相信各位都有所体会。而ICEMCFD特别长于划分六面体网格,相信无论是结构或流体(当然铁别是流体),都会得益于它的威力。 ICEMCFD建模的能力不敢恭维,但划分网格确实有其独到之处。教程开始前,作一个简单的原理介绍,方面没有使用过ICEM CFD的朋友理解主要的任务:111如下图: 1:白色的物体是我们需要划分网格的,但是它非常不规则。 2:这时候你一定想:怎么这个不规则呢,要是它是一个方方正正的形状多好(例如红色的那个形状) 01 111于是有了这样一种思想: 1:对于异型,我们用一种规则形状去描述它。 2:或者说:如果目标形状非常复杂,我们就用很多规则的,简单的形状单元合成在一起,去描述它 之后,将网格划分的设置,做到规则形状上。 最后,这些规则,通过最初的“描述”关系,自动的“映射”到原先的复杂形状上一一问题就得到了解决!!! ICEM CFD E是使用了这种思想 如下是一个三通管,在ProE里做得

在ProE里面直接启动WB 进入WE后,选择如下图: 111如下: 1:代表工作空间里的实体 2:代表某实体的子实体,可以控制它们的开关状态3:控制显示的地方

下面需要创建一个Body实体 这个实体代表了真实的物体。这个真实的物体的外形由我们导入的外形来定义。 ――我们导入的外形并不是真实的实体。这个概念要清楚。 但是今后基本上不会对这个真实的实体作什么操作。这种处理方式主要是为工作空间内有多个物体的时候准备的。 1:点击“创建Body” 2、3:点选这两个点 4:于是创建出一个叫“ Body”的实体 操作中,左键选择,中键确认,右键完成并退出一一类似的操作方法很多地方用到,要多练习,今后就不特别说明了

网格划分模版

生成的网格所能达到的基本指标 1概述 1.1控制网格质量的必要性 在CFD计算中数值误差,也即数值解与微分方程精确解之间的偏差,主要是由截断误差及网格划分不够细密所造成的。而当离散格式的截断误差确定以后,网格的疏密及其分布特性就成了决定离散误差的关键因素。一般在CFD计算中,第一步就是生成计算网格,流场的主要信息都存储在计算网格的节点或者界面上,网格生成质量的高低直接影响着数值分析结果的精度与稳定性。特别是近壁处及通量梯度较大的区域的网格分布最为关键。粗糙的网格会导致数值模拟精度的降低,甚至不能得到收敛解;而过细的网格一方面会耗费过多的计算资源,另一方面也可能导致离散误差的增加,选择适宜的精密网格对于提高计算精度非常关键。因此生成高质量的、适宜的精密网格是获得高精度数值模拟结果的必要条件,在进行CFD计算中必须控制网格的数量及质量。 1.2对计算网格的基本要求 网格分为结构化和非结构化两大类,由于结构化网格在计算精度、计算时间等方面存在相对优势,目前在CFD计算中广泛采用的仍是结构型网格。因此为确保计算结果的正确性及模拟的精度,本课题组要求尽量使用结构化网格,除非在极个别的情况下(如几何结构过于复杂,很难生成结构化网格)才允许使用非结构化网格。 对生成的六面体结构化网格的质量有以下几方面的要求: 首先计算网格中不允许存在负体积,这是保障计算网格正确性的基本要求。 网格单元的总体分布应尽量与主流方向保持一致。 有叶片的区域,应采用绕叶片的O型网格来处理边界层内的流动,另外,O型网格对网格加密很有利。 在所有计算区域的边界处的计算网格线应最大程度的与边界正交,角度最小应大于45°。 计算单元的纵横比不能过大,一般应控制在[1,100]之间,不应高于100。(Aspect Ratio,[1,∞],越接近于1表明网格质量越高)

CFD网格及其生成方法概述

CFD网格及其生成方法概述 作者:王福军 网格是CFD模型的几何表达形式,也是模拟与分析的载体。网格质量对CFD计算精度和计算效率有重要影响。对于复杂的CFD问题,网格生成极为耗时,且极易出错,生成网格所需时间常常大于实际CFD计算的时间。因此,有必要对网格生成方式给以足够的关注。 1 网格类型 网格(grid)分为结构网格和非结构网格两大类。结构网格即网格中节点排列有序、邻点间的关系明确,如图1所示。对一于复杂的儿何区域,结构网格是分块构造的,这就形成了块结构网格(block-structured grids)。图2是块结构网格实例。 图1 结构网格实例 图2 块结构网格实例 与结构网格不同,在非结构网格(unstructured grid)中,节点的位置无法用一个固定的法则予以有序地命名。图3是非结构网格示例。这种网格虽然生成过程比较复杂,但却有着极好的适应性,尤其对具有复杂边界的流场计算问题特别有效。非结构网格一般通过专门的

程序或软件来生成。 图3 非结构网格实例 2 网格单元的分类 单元(cell)是构成网格的基本元素。在结构网格中,常用的ZD网格单元是四边形单元,3D网格单元是六面体单元。而在非结构网格中,常用的2D网格单元还有三角形单元,3D 网格单元还有四面体单元和五面体单元,其中五面体单元还可分为棱锥形(或楔形)和金字塔形单元等。图4和图5分别示出了常用的2D和3D网格单元。 图4 常用的2D网格单元 图5 常用的3D网格单元

3 单连域与多连域网格 网格区域(cell zone)分为单连域和多连域两类。所谓单连域是指求解区域边界线内不包含有非求解区域的情形。单连域内的任何封闭曲线都能连续地收缩至点而不越过其边界。如果在求解区域内包含有非求解区域,则称该求解区域为多连域。所有的绕流流动,都属于典型的多连域问题,如机翼的绕流,水轮机或水泵内单个叶片或一组叶片的绕流等。图2及图3均是多连域的例子。 对于绕流问题的多连域内的网格,有O型和C型两种。O型网格像一个变形的圆,一圈一圈地包围着翼型,最外层网格线上可以取来流的条件,如图6所示。C型网格则像一个变形的C字,围在翼型的外面,如图7所示。这两种网格部属于结构网格。 图6 O型网格 图7 C型网格 4 生成网格的过程

CFD网格的通俗介绍

CFD网格的分类,如果按照构成形式分,可以分为结构化和非结构化 结构化:只能有六面体一种网格单元,六面体顾名思义,也就是有六个面,但这里要区分一下六面体和长方体。长方体(也就是所有边都是两两正交的六面体)是最理想完美的六面体网格。但如果边边不是正交,一般就说网格单元有扭曲(skewed).但绝大多数情况下,是不可能得到完全没有扭曲的六面体网格的。一般用skewness来评估网格的质量,sknewness=V/(a*b*c).这里V是网格的体积,a,b,c是六面体长,宽和斜边。sknewness越接近1,网格质量就越好。很明显对于长方体,sknewness=1.那些扭曲很厉害的网格,sknewness很小。一般说如果所有网格sknewness>0.1也就可以了。结构化网格是有分区的。简单说就是每一个六面体单元是有它的坐标的,这些坐标用,分区号码(B),I,J,K四个数字代表的。区和区之间有数据交换。比如一个单元,它的属性是B=1, I=2,J=3,K=4。其实整个结构化单元的概念就是CFD计算从物理空间到计算空间mapping的概念。I,J,K可以认为是空间x,y,z 在结构化网格结构中的变量。三维网技术论坛! p9 T0 u2 z+ @, i6 c 非机构化:可以是多种形状,四面体(也就三角的形状),六面体,棱形。对任何网格,都是希望网格单元越规则越好,比如六面体希望是长方形,对于四面体,高质量的四面体网格就是正四面体。sknewness的概念这里同样适用,sknewness越小,网格形状相比正方形或者正四面体就越扭曲。越接近1就越好。 很明显非结构化网格也可以是六面体,但非结构化六面体网格没有什么B,IJK的概念,他们就是充满整个空间。 对于复杂形状,结构化网格比较难以生成。主要是生成时候要建立拓扑,拓扑是个外来词,英语是topology,所以不要试图从字面上来理解它的意思。其实拓扑就是指一种有点和线组成的结构。工人建房子,需要先搭房粱,立房柱子,然后再砌砖头。拓扑其实就是房子的结构。这么理解拓扑比较容易些,以后认识多了,就能彻底通了。 生成结构化网格的软件gridgen,icem等等都是需要你去建立拓扑,也就是结构,然后软件好根据你的机构来建立网格,或者砌砖头,呵呵。 非结构化网格的生成相对简单,四面体网格基本就是简单的填充。非结构化六面体网格生成还有些复杂的。但仍然比结构化的建立拓扑简单多。比如gambit的非结构化六面体网格是建立在从一个面到另外一个面扫描(sweep)的基础上的。Numeca公司的hexpress的非结构化六面体网格是用的一种吸附的方法。反正你还是要花点功夫。另外一点就是,结构化网格可以直接应与于各种非结构化网格的CFD软件,比如你在gridgen里面生成了一个结构化网格,用fluent读入就可以了。fluent是非结构化网格CFD软件,它会忽略那些结构化网格的结构信息(也就是B,I,J,K),当成简单的非结构网格读入。非结构化六面体网格就不能用在结构化网格的CFD求解器了 结构化网格仍然是CFD工程师的首选。非结构化六面体网格也还凑合,四面体网格我就不喜欢了。数量多,计算慢,后处理难看。简单说,如果非结构化即快又好,结构化网格早就被淘汰了。 总结一下, 结构化六面体:建立拓扑(所有软件gridgen,icem什么的都是一种拓扑概念,界面不一样罢了),生成网格 非结构化六面体:学习软件,gambit用扫描方法,hexpress用吸附方法,按照步骤就行了。非结构化四面体:简单,看两页教程,搞定,就是简单填充,没什么技术含量!三维网技术论坛) z9 F9 N* `; y" I' l 其他非结构化网格,棱形等等:学习软件,按照步骤,很容易。 不管用什么网格软件,我们最好有比较扎实的CAD(pro/e, solidworks, UG什么的)基础。熟练的CAD技术太重要了。

高性能计算、分布式计算、网格计算、云计算概念与区别

高性能计算、分布式计算、网格计算、云计算--概念和区别 《程序员》2009-02 P34 “见证高性能计算21年” 高性能计算(High Performance Computing)HPC是计算机科学的一个分支,研究并行算法和开发相关软件,致力于开发高性能计算机(High Performance Computer)。 分布式计算是利用互联网上的计算机的中央处理器的闲置处理能力来解决大型计算问题的一种计算科学。 网格计算也是一种分布式计算。网格计算的思路是聚合分布资源,支持虚拟组织,提供高层次的服务,例如分布协同科学研究等。网格计算更多地面向科研应用,商业模型不清晰。网格计算则是聚合分散的资源,支持大型集中式应用(一个大的应用分到多处执行)。 云计算(Cloud Computing)是分布式处理(Distributed Computing)、并行处理(Parallel Computing)和网格计算(Grid Computing)的发展,或者说是这些计算机科学概念的商业实现。云计算的资源相对集中,主要以数据中心的形式提供底层资源的使用,并不强调虚拟组织(VO)的概念。云计算从诞生开始就是针对企业商业应用,商业模型比较清晰。云计算是以相对集中的资源,运行分散的应用(大量分散的应用在若干大的中心执行);

目录 高性能计算、分布式计算、网格计算、云计算--概念和区别 (1) 高性能计算 (3) 百科名片 (3) 概念 (3) 服务领域 (3) 网格 (5) 百科名片 (5) 网格的产生 (5) 网格技术的特征及其体系结构 (5) 高性能计算机的发展与应用 (17) 我国高性能计算机应用前景及发展中的问题 (17) 高性能计算机与大众生活息息相关 (17) 高性能计算机发展任重道远 (18) 分布式计算、网格计算和云计算 (21) 分布式计算 (21) 网格计算 (21) 云计算 (22) 网格计算和云计算的概念和区别 (24) 目标不同 (24) 分配资源方式的不同 (25) 殊途同归 (26) 钱德沛教授:云计算和网格计算差别何在? (27) 云计算与网格计算的概念 (27) 网格计算的特点是什么呢? (27) 云计算与网格计算区别何在 (28)

网格划分的几种基本处理方法

网格划分的几种基本处理方法 贴体坐标法: 贴体坐标是利用曲线坐标,并使其坐标线与燃烧室外形或复杂计算区域边界重合,这样所有边界点能够用网格点来表示,不需要任何插值。一旦贴体坐标生成通过变换,偏微分方程求解可以不在任意形状的物理平面上,而在矩形或矩形的组合(空间问题求解域为长方体或它们的组合)转换平面上进行。这样计算与燃烧室外形无关,也与在物理平面上网格间隔无关。 而是把边界条件复杂的问题转换成一个边界条件简单的问题;这样不仅可避免因燃烧室外形与坐标网格线不一致带来计算误差,而且还可节省计算时间和内存,使流场计算较准确,同时方便求解,较好地解决了复杂形状流动区域的计算,在工程上比较广泛应用。 区域法: 虽然贴体坐标系可以使坐标线与燃烧室外形相重合,从而解决复杂流动区域计算问题。但有时实际流场是一个复杂的多通道区域,很难用一种网格来模拟,生成单域贴体网格,即使生成了也不能保证网格质量,影响流场数值求解的效果。因此,目前常采用区域法或分区网格,其基本思想是,根据外形特点把复杂的物理域或复杂拓扑结构的网格,分成若干个区域,分别对每个子区域生成拓扑结构简单的网格。由这些子区域组合而成的网格,或结构块网格。对区域进行分区时,若相邻两个子域分离边界是协调对接,称为对接网格;若相邻两子域有相互重叠部分,则此分区网格称为重叠网格。根据实际数值模拟计算的需要,把整个区域(燃烧室)分成几个不同的子区域,并分别生成网格。这样不仅可提高计算精度,而且还可节省计算机内存,提高收敛精度。但是计算时,必须考虑各区域连接边界处耦合以及变量信息及时、准确地传递问题。处理各个区域连接有多种方法,其中一个办法是在求解各变量时各区域可以单独求解若干次而对压力校正方程.设压力校正值在最初迭代时为零,为了保证流量连续各个区域应同时求解,然后对各个速度和压力进行校正。或者采用在两个区域交界处有一个重叠区,两个区域都对重叠区进行计算,重叠区一边区域内的值,要供重叠区另一边区域求解时用。或通过在重叠内建立两个区域坐标对应关系,实现数据在重叠区内及时传递。如果两个区采用网格疏密分布不相同,要求重叠区二边流量相等。区域法能合理解决网格生成问题,已被大量用来计算复杂形状区域流动。 区域分解法: 对于复杂几何形状的实际燃烧装置,为了保证数值求解流场质量,目前常采用区域分解法。该法基本要点是:根据燃烧室形状特点和流场计算需要,把计算区域分成一个主区域和若干个子区域,对各个区域(块)分别建立网格,并对各个区域分别进行数值求解。区域分解原则是尽量使每个子区域边界简便以便于网格建立,各个子区域大小也尽可能相同,使计算负载平衡有利于平行计算。各区域的网格间距数学模型以及计算方法都可以不同,通常在变量变化梯度大的区域,可以布置较细网格,并采用高阶紊流模型和描述复杂反应的紊流燃烧模型,以便更合理模拟实际流场。对于变量变化不太大区域,可采用较疏的网格和较简单的数学模型,这样可节省计算时间。各子区域的解在相邻子区域边界处通过耦合条件来实现光滑,相邻子区域连接重叠网格或对接网格来实现,在各子区域交界处通过插值法提供各子域求解变量的信息传递,满足各子域流场计算要求通量和动量守恒条件以便实现在交界面处各子域流场解的匹配和耦合,从而取得全流场解。 非结构网格法: 上述各方法所生成的网格均属于结构化网格,其共同特点是网格中各节点排列有序,每个节点与邻点之间关系是固定的,在计算区域内网格线和平面保持连续。特别是其中分区结构网格生成方法已积累了较多经验,计算技术也较成熟,目前被广泛用来构造复杂外形区域

网格和单元的基本概念

网格和单元的基本概念 前记:首先说明,和一般的有限元或者计算力学的教材不一样,本人也不打算去抄袭别人的著作,下面的连载是一个阶段的学习或者专业感悟集大成,可以说深入浅出,也可以说浅薄之极——如果你认为浅薄,很好,说明我理解透了,也祝贺你理解透了!好了,废话少说,书归正传。 无论是CSD(计算结构力学)、CTD(计算热力学)还是CFD(计算流体动力学)——我们统一称之为工程物理数值计算技术。支撑这个体系的4大要素就是:材料本构、网格、边界和荷载(荷载问题可以理解为数学物理方程的初值问题),当然,如果把求解技术也看作一个要素,则也可以称之为5大要素。网格是一门复杂的边缘学科,是几何拓补学和力学的杂交问题,也是支撑数值计算的前提保证。本番连载不做任何网格理论的探讨(网格理论是纯粹的数学理论),仅限于尽量简单化的应用技术揭秘。 网格出现的思想源于离散化求解思想,离散化把连续求解域离散为若干有限的子区域,分别求解各个子区域的物理变量,各个子区域相邻连续与协调,从而达到整个变量场的协调与连续。离散网格仅仅是物理量的一个“表征符号”,网格是有形的,但被离散对象既可以是有形的(各类固体),也可以是无形的(热传导、气体),最关键的核心在于网格背后隐藏的数学物理列式,因此,简单点说,看得见的网格离散是形式,而看不见的物理量离散才是本质核心。 对计算结构力学问题,网格剖分主要包含几个内容:杆系单元剖分(梁、杆、索、弹簧等)、二维板壳剖分(曲面或者平面单元)、三维实体剖分(非结构化全六面体网格、四面体网格、金字塔网格、结构化六面体网格、混合网格等),计算热力学和计算流体动力学的网格绝大部分是三维问题。对于CAE工程师而言,任何复杂问题域最终均直接表现为网格的堆砌,工程师的任务等同于上帝造人的过程,网格是一个机体,承载着灵魂(材料本构、网格、边界和荷载),求解技术则是一个思维过程。 网格基本要素是由最基本的节点(node)、单元线(edge)、单元面(face)、单元体(body)构成,实质上,线、面、体只不过是为了让网格看起来更加直观,在分析求解过程中,线、面、体本质上并没有起多大的作用,数值离散的落脚点在节点(node)上,所有的物理变量均转化为节点变量实现连续和传递。在所有的CAE环境下,网格的基本要素均可以直接构成,但对于复杂问题而言,这是一个在操作上很难实现的事情,因此,基于几何要素的网格划分技术成为现代网格剖分应用的支点,和网格基本要素完全相同,对应的几何要素分别称之为点(point)、线(curve)、面(surface)和实体(solid)。 数值离散求解器是不能识别几何元素的,要对其添加“饲料”,工程师必须对几何元素进行“精加工”,因此,从这个意义上来说,网格剖分的本质就是把几何要素转换为若干离散的元素组,这些元素组堆砌成形态上近似逼近原有几何域的简单网格集合体。因此,这里说明了一个网格“加工”质量的基本判别标准——和几何元素的拟合逼近程度,理论上,越逼近几何元素的网格质量越好,当然,几何逼近只是一个基本的判别标准,网格质量判别有一系列复杂的标准,后文详细阐述。 本篇将专门解释几个基本概念:点网格;一维线网格;二维三角形面网格、二维四边形面网格;三维四面体网格(tetrahedra)、三维金字塔单元(pyramid)、五面体单元(prism)、三维六面体单元(hexahedra);结构化网格(structural grid)、非结构化网格(nonstructural grid)、混合网格(blend grid)。需要专门

ICEMCFD网格划分入门基础

WorkBench ICEM CFD 网格划分入门 111AnsysWB里集成了一个非常重要的工具:ICEM CFD。 它是一个建模、划分网格的集成工具,功能非常强大。我也只是蜻蜓点水的用了几次,感觉确实非常棒,以前遇到复杂的模型,用过几个划分网格的工具。但这是我觉得最方便和最具效率的。 网格划分很大程度上影响着后续的仿真分析——相信各位都有所体会。而ICEM CFD特别长于划分六面体网格,相信无论是结构或流体(当然铁别是流体),都会得益于它的威力。 ICEM CFD建模的能力不敢恭维,但划分网格确实有其独到之处。教程开始前,作一个简单的原理介绍,方面没有使用过ICEM CFD的朋友理解主要的任务: 111如下图: 1:白色的物体是我们需要划分网格的,但是它非常不规则。 2:这时候你一定想:怎么这个不规则呢,要是它是一个方方正正的形状多好(例如红色的那个形状)01 111于是有了这样一种思想: 1:对于异型,我们用一种规则形状去描述它。 2:或者说:如果目标形状非常复杂,我们就用很多规则的,简单的形状单元合成在一起,去描述它。 之后,将网格划分的设置,做到规则形状上。 最后,这些规则,通过最初的“描述”关系,自动的“映射”到原先的复杂形状上——问题就得到了解决!!! ICEM CFD正是使用了这种思想。 如下是一个三通管,在ProE里做得

02 在ProE里面直接启动WB 进入WB后,选择如下图: 03 111如下: 1:代表工作空间里的实体 2:代表某实体的子实体,可以控制它们的开关状态3:控制显示的地方

04 下面需要创建一个Body实体 这个实体代表了真实的物体。这个真实的物体的外形由我们导入的外形来定义。 ——我们导入的外形并不是真实的实体。这个概念要清楚。 但是今后基本上不会对这个真实的实体作什么操作。这种处理方式主要是为工作空间内有多个物体的时候准备的。 05 1:点击“创建Body” 2、3:点选这两个点 4:于是创建出一个叫“Body”的实体 操作中,左键选择,中键确认,右键完成并退出——类似的操作方法很多地方用到,要多练习,今后

ICEM_CFD混合网格

ICEM CFD中合并多个网格 对于结构十分复杂的几何模型,若能够将几何体分割成多个部分由多人分别进行网格划分,生成网格后能够对网格进行组装,这恐怕是很多人梦寐以求的功能了。其实很多前处理软件都具有此功能。今天要说的是如何在ICEM CFD中实现此功能。 为了简单起见,这里用一个非常简单的模型进行演示。当然复杂的模型的处理方式也是相同的。我们要处理的几何模型如图1所示。一个L型整体块被切割成3份。分别导出为3个不同的几何文件。按图中标示的顺序分别导出为1.x_t,2.x_t,3.x_t,当然其他的格式也无妨。但是最好是在同一个体上进行切割,否则网格组装的过程中会存在定位的问题。同一个体上切割的几何则不会存在几何坐标定位的问题。 图1 原始几何图2 几何1生成的网格图3 保存网格 1、将几何1.x_t导入到ICEM CFD中进行网格划分。注意千万保证单位的一致,切记。 这里是一个长方体,网格划分方法就不多说了。预览网格如图2所示。选择菜单File > Mesh > Load From Blocking生成网格。 2、保存网格。 选择File > Mesh >Save Mesh As…,我们这里保存已生成的网格为1.uns,后面组装的时候要用到此文件。 3、按照相同的步骤对模型2与模型3进行网格文件,同时保存网格文件为2.uns与3.uns。

图4 模型2的网格图5 模型3的网格 4、网格组装 先导入1.uns,点击菜单File > Mesh >Open Mesh…,选择第2步保存的网格文件1.uns,导入模型1的网格。 以同样的菜单,选择2.uns,会弹出对话框如图6所示。注意此时选择Merge,否则如果选择Replace的话,则只会导入模型2的网格,将模型1的网格替换掉,这不是我们想要的。接下来我们以相同的步骤导入3.uns,同样选择Merge。导入后网格如图7所示。 图6 对话框图7 全部倒入后的模型 5、导出网格 以常规方式导出网格。我们这里测试将网格导入至少fluent中。从图8导入信息可以看到,完全没有问题。

云计算和网格计算有什么本质区别

云计算和网格计算有什么本质区别 https://www.sodocs.net/doc/f916629400.html,/z/q157731426.htm?w=%CD%F8%B8%F1%BC%C6%CB%E3%BC%BC%CA %F5&spi=1&sr=1&w8=%E7%BD%91%E6%A0%BC%E8%AE%A1%E7%AE%97%E6%8A%80%E6%9 C%AF&qf=10&rn=360 [标签:云计算,本质区别,区别] 我对云了解的比较深入,对网格 计算不太了解,但是初步观察发现相 似之处很多,求解两者本质区别 限量版回答:4 人气:108 解决 时间:2009-10-03 20:35 满意答案 耐心看吧 您可能非常关注云计算和网格计 算的比较。本文介绍了云计算服务类 型,云计算和网格计算的相似与不同。 同时本文探讨了云计算优于网格计算 的地方,两者面临的共同问题以及一 些安全方面的问题。本文以Amazon Web Services 为例。 实现云计算需要三个部分:瘦客 户机(或者能够在胖瘦之间切换的客 户机)、网格计算和效用计算。网格 计算将独立的计算机连接成一个大的 基础设施,充分利用闲置的资源。效用计算就是支付在共享服务器上使用的服务,就好象支付公共事业一样(比如电力、天然气等)。 通过网格计算,可以把计算资源作为能够开启关闭的公用事业来提供。云计算更进一步,可以随需提供计算资源。这样在使用公用定价时就可以避免过度供给。在满足数百万用户的需求时也消除了过度供给的需要。 基础设施即服务 消费者通过Internet 可以从完善的计算机基础设施获得服务。这类服务称为基础设施即服务(Infrastructure as a Service,IaaS)。基于Internet 的服务(如存储和数据库)是IaaS 的一部分。Internet 上其他类型的服务包括平台即服务(Platform as a Service,PaaS)和软件即服务(Software as a Service,SaaS)。PaaS 提供了用户可以访问的完整或部分的应用程序开发,SaaS 则提供了完整的可直接使用的应用程序,比如通过Internet 管理企业资源。 作为Infrastructure as a Service (IaaS) 在实际应用中的一个例子,The New York Times 使用成百上千台Amazon EC2 实例在36 小时内处理TB 级的文档数据。如果没有EC2,The New York Times 处理这些数据将要花费数天或者数月的时间。 IaaS 分为两种用法:公共的和私有的。Amazon EC2 在基础设施云中使用公共服务器池。更加私有化的服务会使用企业内部数据中心的一组公用或私有服务器池。如果在企业数据中心环境中开发软件,那么这两种类型都能使用,而且使用EC2 临时扩展资源的成本也很低—比方说测试。结合使用两者可以更快地开发应用程序和服务,缩短开发和测试周期。 Amazon Web 服务

网格技术与网格计算

引言:网格作为一种能带来巨大处理、存储能力和其他IT资源的新型网络,可以应付临时之用。网格计算通过共享网络将不同地点的大量计算机相联,从而形成虚拟的超级计算机,将各处计算机的多余处理器能力合在一起,可为研究和其他数据集中应用提供巨大的处理能力。有了网格计算,那些没有能力购买价值数百万美元的超级计算机的机构,也能利用其巨大的计算能力。 一. 网格技术 1.计算的“乌托邦”? Gartner公司的Rob Batchelder认为,网格的构想一直是计算领域的“乌托邦”,在科技应用上虽有巨大前景,但最大的缺陷是缺乏明显的商业应用。自20世纪90年代在欧美出现以来,网格主要被用于帮助分散的大学研究人员分析粒子加速器和巨型望远镜的数据。但在过去的两年中,网格的概念和GlobusToolkit已在研究和教育领域得到广泛应用,数十项全球性的大项目采用这些技术,以挑战科学计算中的海量计算问题。 目前网格技术虽主要为学术机构所控制,但企业也在陆续跟进。事实上,全球网格论坛(GlobalGridForum)的主要赞助企业就包括Unilever——一家以经销肥皂、冰淇淋著称的企业。与许多正在研究和评估网格技术的企业一样,Unilever自己对于如何利用此技术仍秘而不宣。而Johnson&Johnson与Merck等制药公司、BMW与波音等制造企业却已利用这一技术的处理能力和存储空间进行仿真试验,例如药品能否保护细胞免受病毒侵袭?飞机机翼是否会在暴风雨中折断? 基因研究是网格技术的自然应用,这一领域所需的投资很难由一家企业来承担,生物科技企业可用网格技术来分析基因数据;医生可以用网格技术制作出病人器官的三维模型,作为诊断疾病的辅助手段;网格可以处理来自商店现金记录或金融市场的数据流。其他行业,如航空、保险、运输和国防,也会从中受益。如此看来,网格计算并非是可望不可及的乌托邦,其商业应用的广阔前景就在眼前。

将ICEM CFD网格导入Fluent的方法

将ICEM CFD网格导入Fluent的方法 将ICEM CFD网格导入Fluent的方法及生成CFD边界条件的2种方法 ICEM CFD软件具有强大的流体网格生成能力,这是一个业内普遍的认识。但是,ICEM CFD创建网格的基本理念与大多数软件不一样,比如说ICEM CFD中基本上没有几何体的概念,又比如说网格不直接依赖于几何模型而是更依赖于Block,这和Patran、Ansys等常用软件的几何模型+网格模型、网格依赖于几何模型等的普遍概念完全不同。因此,对于初学ICEM CFD的人来说,这个软件入门很难,这也是作者Dengguide这半年学习ICEM CFD的深切体会。在这篇短文中,作者主要想介绍一下比较困扰新手的2个问题:(1)怎么将ICEM CFD网格导入Fluent中?(2)怎么在ICEM CFD中生成出入口、耦合面等边界条件? 1.将ICEM CFD网格导入Fluent ICEM CFD默认输出网格文件的后缀是.uds,不是Fluent能够读入的.msh格式。设置方法如下:(1)点击output->点击select solver,在Output Solver选项中选择Fluent_V6,Common Structural Solver通常默认为ANSYS,不用理它;(2)点击Write input,就是第四个图标,会弹出“Save current project first”对话框,选Yes保存uns文件,然后弹出“Fluent V6”对话框,其他选项都可不改动,可在Output file这一项中改一个你想要的文件名称,最后点Done选项就好了,输出的是mesh文件,可以直接读入Fluent中。 2.在ICEM CFD中生成CFD边界条件 迄今为止,作者总共摸索出了2种在ICEM CFD中生成CFD边界条件生成CFD边界条件的方法。在介绍之前,先说明2个基本提示:(1)ICEM输出的mesh文件只是包含了节点坐标、节点之间的连接信息以组成网格和表面(faces)、所有表面的区域类型(zone type)和数量,并不包含几何信息;(2)mesh文件导入后Fluent后,Fluent自动将体单元组的内部单元和表面分成2组,体单元1个组,面单元1个组。下面继续介绍三种方法。 方法一 导入几何模型后,在右上角树状图中右键点击Parts,选择Creat Part项,选择对应于边界条件的面,定义一个Part,有n个边界条件面,就定义n-1个Part,生成网格后各个part 就对应于边界网格面,最后一个边界条件面Fluent会自动生成。这种方法要求在划分网格之前下定义后边界对应面Part,要是忘了怎么办呢?没关系还有第二种方法。 方法二 导入几何模型后,定义网格参数,划分网格,然后再来定义边界面网格。根据提示一,mesh文件值输出网格信息,没有几何信息。在划分好网格后,我们可以手动选择对应的边界网格,并将其定义为Part,这样输出导入Fluent后就会保留在ICEM CFD中定义的Part。定义网格Part的方法与方法一相同,需要提示大家的是“Select geometry”工具默认显示几何模型选择界面,需要先点击最后一个图框代表的“Toggle mesh tools”,转化到“Select mesh elements”界面才能选择网格,通过圆圈或者方框选项手动选择边界网格,建立对应的Part 即可。这种方法的缺点是选起来比较多比较累。 那有没有更好的方法呢?也许有,不过作者还没有找到,如果找到了再告诉大家,或者大家找到了,再告诉我吧!另外关于耦合面,需要定义不同块之后才能生成,具体实现方法以后再整理吧!

网格划分的原则

划分网格是建立有限元模型的一个重要环节,它要求考虑的问题较多,需要的工作量较大,所划分的网格形式对计算精度和计算规模将产生直接影响。为建立正确、合理的有限元模型,这里介绍划分网格时应考虑的一些基本原则。 1网格数量 网格数量的多少将影响计算结果的精度和计算规模的大小。一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因数综合考虑。 图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。所以应注意增加网格的经济性。实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。 图1位移精度和计算时间随网格数量的变化 在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格,如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。 2网格疏密 网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。这样,整个结构便表现出疏密不同的网格划分形式。 图2是中心带圆孔方板的四分之一模型,其网格反映了疏密不同的划分原则。小圆孔附近存在应力集中,采用了比较密的网格。板的四周应力梯度较小,网格分得较稀。其中图b 中网格疏密相差更大,它比图a中的网格少48个,但计算出的孔缘最大应力相差1%,而计算时间却减小了36%。由此可见,采用疏密不同的网格划分,既可以保持相当的计算精度,又可使网格数量减小。因此,网格数量应增加到结构的关键部位,在次要部位增加网格是不必要的,也是不经济的。

浅谈网格计算相关技术与应用(一)

浅谈网格计算相关技术与应用(一) 论文关键词]数据库浮点运算虚拟化资源共享论文摘要]论述网格计算的发展概况,在科学领域的应用范围,网格服务的特点以及在未来网络下场中的发展潜力。 一、网格计算的由来与发展 网格计算是伴随着互联网技术而迅速发展起来的,是将地理上分布的计算资源(包括数据库、贵重仪器等各种资源)充分运用起来,协同解决复杂的大规模问题,特别是解决仅靠本地资源无法解决的复杂问题,是专门针对复杂科学计算的新型计算模式。这种计算模式是利用互联网把分散在不同地理位置的电脑组织成一个“虚拟的超级计算机”,其中每一台参与计算的计算机就是一个“节点”,而整个计算机是由成千上万个“节点”组成的“一张网格”,所以这种计算方式叫网格计算。这样组织起来的“虚拟的超级计算机”有两个优势,一个是数据处理能力超强,另一个是能充分利用网上的闲置处理能力。简单地讲,网格是把整个网络整合成一台巨大的超级计算机,实现计算资源、存储资源、数据资源、信息资源、知识资源、专家资源的全面共享。 近几年,随着计算机计算能力的迅速增长,互联网络的普及和高速网络成本的大幅降低以及传统计算方式和计算机的使用方式的改变,网格计算已经逐渐成为超级计算发展的一个重要趋势。网格计算是一个崭新而重要的研究领域,它以大粒度资源共享、高性能计算和创新性应用为主要特征,必将成为21世纪经济发展的重要推动力。 二十世纪九十年代以来,世界各个国家,尤其是发达国家,建立了很多超级计算应用中心和工程研究中心,美国还制定了新一轮规划的先进计算框架(ACIP),发展面向21世纪的先进计算技术。我国在科技部的领导和主持下,经过专家组及相关单位的努力,作为我国高性能计算和信息服务战略性基础设施的国家高性能计算环境发展很快。在已经建成的5个国家级高性能计算中心的基础上,又于中南、西北等地建立了新的国家高性能计算中心,科技部加强了网络节点的建设,形成了以科学院为主体的计算网格。教育部也启动了网格计算工程,第一批12个网点正在建设之中,国家基金委也列出专项基金资助网格计算。 网格是借鉴电力网的概念出来的,网格的最终目的是希望用户在使用网格的计算能力时,就如同现在使用电力一样方便简单。 二、网格计算的应用 (一)分布式超级计算 网格计算可以把分布式的超级计算机集中起来,协同解决复杂大规模的问题。是大量的闲置计算机资源得到有效的组织,提高了资源的利用效率,节省了大量的重复投资,使用户的需求能够得到及时满足。 (二)高吞吐率计算机 网络技术能够十分有效地提高计算的吞吐率,它利用CPU周期窃取技术,将大量闲置计算机的计算资源集中起来,提供给对时间不太敏感的问题,作为计算资源的重要来源。 (三)数据密集型计算 数据密集型计算的问题求解通常同时产生很大的通讯和计算需求,需要网格能力才可以解决。网格已经在药物分子设计、计算力学、计算材料、电子学、生物学、核物理反映、航空航天等众多领域得到广泛应用。 (四)给予更广泛信息共享的人与人交互 网格的出现更急突破了人与人之间地理界线的限制,使得科技工作者之间的交流更加的方便,从某种程度上说,可以实现人与人之间的智慧共享。 (五)更广泛的资源贸易 随着大型机性能的提高和微机的更加普及,其资源的闲置问题越来越突出,网络技术可以有效地组织这些闲置资源,使得有大量的计算需求用户能够获得这些资源,而资源提供者的应

浅析网格计算技术的应用与发展

浅析网格计算技术的应用与发展 郑琦1卢德利2 (1:吉林建筑工程学院,长春130021;2:吉林医药设计院有限公司,长春130021) 摘要:网格计算是伴随着互联网而迅速发展起来的、专门针对复杂科学计算的新型计算模式。我们可以将服务器、存储系统和网络联合在一起,数据文件、应用程序和系统看起来就像是一个巨大的虚拟计算系统,为用户提供功能强大的多系统资源来处理特定的任务。本文简述了网格计算的概念、关键组件,介绍了网格计算在若干领域的应用,概述了网格计算现状和发展前景。 关键字:网格计算;分布式系统;资源共享 中图分类号:文献标识符:文章编号: Analyze Briefly The Application and Development of Grid Computing ZhengQi1 Lu De-li2 (1:JiLin Institute of Architecture and Civil Engineering ,Changchun 130021; 2:JiLin Pharmaceutical Design Institute Co.Ltd,Changchun 130021) Abstract:Grid computing is a new computation pattern that aims at the complex science computation specifically, accompanied by the rapid development of the Internet. We can be able to servers, memory systems and networks together. It can provide users with a powerful multi-system resources to handle specific tasks, like a huge virtual computing systems. This paper has summarized the concept and the key components of grid computing.Introduced a number of grid computing in the field of application. Grid computing provides an overview of the current situation and development prospects. Key Words: Grid Computing; Distributional System; Resource Sharing 0 引言 随着超级计算机的不断发展,它已经成为复杂科学计算领域的主宰。但以超级计算机为中心的计算模式存在明显的不足,而且目前正在经受挑战。超级计算机虽然具有非常强大处理能力,但它造价极高,通常只有一些国家级的如航天、国防等部门才有能力配置这样的设备。而随着人们日常工作遇到的商业计算、科学计算越来越复杂,人们越来越需要数据处理能力更强大的计算机,而超级计算机的价格使得它进入普通人的工作领域成为一种奢望。于是人们在寻找造价低廉而数据处理作者简介:郑琦(1978~),女,吉林省长春市人,助教,硕士

相关主题