发布时间 : 星期一 文章CFXANSYS使用高手流固耦合至理名言更新完毕开始阅读e9f1d0d428ea81c758f57873
可以用简单的方法测试出来： 1)定义一个简单的边界条件，计算一下，看看流场是不是连续。事实上，如果网格不连续的话，在CFX-Solver刚开始求解时，会提示有isolated fluid region。 2)如果网格是不连续的，CFX-Pre会自动把界面处的网格放在Domain Default边界条件里面。如果网格是连续的，界面处的网格会被默认为是连通的，被忽略掉。因此，这要选中Domain Default，看看里面有没有界面网格就行了。 CGNS网格里的界面会被忽略，而Patran网格的界面会被放在Domain Default里，当成是壁面条件。 另外我试用了用ICEM CFD把CGNS structured mesh转换成unstructured mesh。ICEM CFD虽然能完成转换，但是不能更新part名，还需要看看有什么其他解决方法。 Technorati : CFX, CGNS, 网格
http://blog.singinfo.com/?p=102
Tue, 21 Aug 2007 05:17:00 +0000
Workbench环境下CFX-Post背景色的设定
CFX 11.0提供的report功能可以快速的在CFX-Post里后处理结果并且发布为html格式。用户可以方便的更新模型的方向、显示方式等选项。Report里的设定可以重复使用在相同的模型上，而不必每次重复操作。 用户可以更改report开头的ANSYS图标(logo)，但是每个图片右上角的ANSYS图表是去不掉的(如果有谁知道怎么去掉这个图标，欢迎email给我)。 问题是这样的，默认的图片背景色是Workbench默认的蓝白过渡背景色，通常情况下，我们可以通过CFX-Post -> Edit -> Options来修改(如下图所示)，但是实际上用户看不到改变背景色的选项(另一幅截图)。 当CFX单独运行时，选项如第一幅图；当CFX在Workbench里运行时，选项如第二幅图。 当在Workbench里运行时，需要通过修改Workbench的背景色来设定CFX背景色。[Project] tab -> Tools -> Options -> Common Settings -> Graphics Style。 一点题外话，CFX-Post的report功能还是新功能，还不太稳定。比如在comment里插入的图片在发布后不能正常显示。这个问题可能在以后的版本中得到解决。 Technorati : CFX, Workbench, 背景色 http://blog.singinfo.com/?p=96
Tue, 14 Aug 2007 03:21:24 +0000
Flomerics发布EFD 8
新闻连接在这里。 EFD是Engineering Fluid Dynamics的缩写，基于和通用CFD同样的数学理论。同时和CAD软件有更好的协作性，使用也相对容易。 Technorati : EFD, Flomerics
http://blog.singinfo.com/?p=95
Sun, 12 Aug 2007 12:36:51 +0000
CFX：最大堆栈长度
在求解很大的CFX模型时，可能会遇到下面的错误信息#001100279(其实在正常求解时，求解器也会给出内存分布的信息)。 这个信息其实并不是因为系统物理内存不够。而是由于实型堆栈(real stack)容量不足。在CFX里，堆栈的最大容量是2^31-1。任何堆栈大于这个长度，系统都会给出这个错误信息，即使系统内存足够。这个时候解决问题的唯一办法是使用并行计算。 怎么解读这个错误信息？ 错误消息的第一行是：实型堆栈长度为2481110.5千字(word)，平均每个节点需要318.47字，每个单元需要330.24字，总共需要的内存是9691838.0千字节，平均每个节点需要1273.89字节。 这个模型使用的是六面体网格，因此单元和节点的比例为1:1。模型有大概8百万节点。需要9.6GB内存来
存储实型变量，2.7GB内存来存储整型变量，总共需要12GB多内存。 在CFX里，任何一个堆栈的长度都不能超过2147483647。 +——————————————————————–+ | Memory Allocated for Run (Actual usage may be less) | +——————————————————————–+ Data Type Kwords Words/Node Words/Elem Kbytes Bytes/Node Real 2481110.5 318.47 330.24 9691838.0 1273.89 Integer 712176.8 91.41 94.79 2781940.8 365.66 Character 2480.0 .32 .33 2421.9 .32 Logical 40.0 .01 .01 156.2 .02 Double 1208.0 .16 .16 9437.5 1.24 +——————————————————————–+ | ERROR #001100279 has occurred [...] http://blog.singinfo.com/?p=90
Sun, 05 Aug 2007 06:33:08 +0000
CFX：浮力和对称边界条件
前面有两篇文章(文章1、文章2)提到大气边界层的定义。流场的顶面(天空)可以定义成对称边界条件。 在有些情况下，可能需要考虑浮力项，比如在大气边界层内有一个热源。在Domain上加上浮力项后，求解器可能会给出下面的错误信息： ERROR #001100279 has occurred in subroutine ErrAction. Message:When buoyancy is active, all symmetry plane normals must be perpendicular to the gravity vector. Writing crash recovery file ERROR #001100279 has occurred in subroutine ErrAction. Message: Stopped in routine ASS_GLOB_SYM 起初我以为这是由于软件问题，在把expert parameter里的vector parallel tolerance设为90或大于90后，可以让求解器继续运行。 这个错误信息说的是，如果激活浮力项，所有对称面的法线方向必须要和重力矢量垂直。也就是说对称面必须和重力方向平行。因为有浮力项时，在重力方向会有对流速度。如果对称面和重力不平行的话，在对称面上会有一个垂直于该面的速度分量。而对称面的物理意义是速度场只和该面相平行。这两点之间是矛盾的，因此对称面必须和重力方向相平行。 因此要解决上面的问题，需要定义一个合适的边界条件，而不是修改expert parameter。 Technorati : CFX, 对称边界条件, 浮力 http://blog.singinfo.com/?p=91
Sat, 04 Aug 2007 07:22:00 +0000
本博客主要介绍ANSYS，CFX，Fluent的最新更新，使用技巧等。同时也有CFD和FEA基本理论，及其他CAE软件的介绍。
Airpak收敛问题的解决
由于Airpak使用Fluent求解器，默认情况下对流项使用一阶迎风(first order upwind)离散格式。通常情况下，收敛性还是很好的。最近帮忙解决了几个收敛问题，现在把解决方法列在下面，本帖将持续更新。 1. slack value。在使用Airpak 3.0做non-conformal mesh的时候，需要填写slack值。slack值可以让assembly的范围比Airpak对象稍大。有时候，assembly和room之间会有一个很小的缝隙，使网格质量很差。 由于slack值是一个相对值，当更新了模型的尺寸后，assembly的大小也会随之变化，这样也容易导致小的缝隙。 2. 多面体(ploygon)高度的差别。比如有两个多面体，起始高度都为z = 0 m。如果第一个多面体的高度为10 m，第二个多面体的高度为10.01 m。这两者之间会有0.01 m的缝隙，网格质量会很低。 要检查网格质量，可以检查face alignment。把网格质量范围设成0.0到0.15之间，检查有没有质量低于0.15的网格。如果有的话，可以点中相应的柱状图，检查低质量网格所在的位置。另外align face功能可以帮助精确对齐表面。
另外使用较小的动量和压力松驰因子也有助于加快收敛。 http://blog.singinfo.com/?p=92
Fri, 03 Aug 2007 07:41:00 +0000
LES/DES里可以用对称或周期边界条件吗？
很多时候，我们可以使用对称性边界条件或者周期性边界条件来简化问题、节省计算时间。通常情况下，对称性或周期性是指几何形状和流场两者都具有对称性/周期性。当然，对称性也适用于几何形状不对称，而流场对称的情形(速度和边界平行，标量在垂直于边界的方向上梯度为0)。 在做LES/DES时，通常不能使用对称或周期性边界条件。因为平均速度场也许是对称或周期性，但是湍流是没有规律、紊乱的，它们不可能是对称或周期性的。 如果用LES/DES模拟透平机械内部或者汽车外部流动，需要模拟所有的叶片以及整个汽车，而不可以用对称或周期性条件来进行简化。 但是当边界距离重要区域比较远的时候，对称性/周期性边界条件是可以使用的。比如做汽车外部流场分析时，需要建立一个相当大的外部流畅，除了进口、入口和地面外，两个侧面和顶面可以使用对称条件。虽然这种情况下对称边界条件并不成立，但是相对而言，对称条件比较接近实际的远场流场。 同样，有的时候位移周期性边界条件也是可以接受的。比如，用DNS/LES模拟管道内流场时，带有压力差的位移性边界条件是很常用的。 Technorati : DES, LES, 周期性, 对称, 边界条件 http://blog.singinfo.com/?p=93
Thu, 02 Aug 2007 09:47:00 +0000
Flowizard 3.0.4发布
Fluent的产品发布好像都是无声无息的，只有经常查看用户服务中心才知道。最新版本的Flowizard是3.0.4，发布于7月。 New Features in FloWizard V3.0.4 CAD Connections FloWizard V3 provides CAD connections for Pro/ENGINEER, SolidWorks, Autodesk Inventor, and UGS’ NX users. These CAD connections let you send your geometry from your CAD software to FloWizard with the click of a button. This reduces the time it takes to generate CFD models from your CAD models. One license [...] http://blog.singinfo.com/?p=94
Wed, 01 Aug 2007 13:32:00 +0000
我的email地址
昨天夜里把我的email地址写在一片帖子里，不到10个小时就有垃圾邮件。 终于找了一个网站，把email地址做成图片。 刚刚看了一些CFD的幻灯片，CFD的应用还是很浩渺的！有问题大家一起讨论。 Technorati : email http://blog.singinfo.com/?p=87
Sun, 29 Jul 2007 06:53:07 +0000
CFX：求解时实时显示阻力系数和升力系数
问题是怎么用CEL计算阻力系数(drag coefficient, cd)和升力系数(lift coefficient, cl)，并且在求解时实时显示cd和cl。 碰巧遇到一个例子，现在把CEL贴在下面。应该把这个CEL稍微修改一下，可以轻松应用到其他模型上。 在CFX-Pre里定义一个新的monitor point的方法是点击”Create Output Files and Monitor Points”图标，然后到Monitor tab下面把CEL表达式定义成一个monitor point。 LIBRARY:CEL:EXPRESSIONS:DragForc
e = 2*(force_x()@DLRF6*cos(alpha)+force_z()@DLRF6*sin(alpha))DynamicPres = 0.5*ReferenceDensity*Vinlet^2LiftForce = max(2*(force_z()@DLRF6*cos(alpha)-force_x()@DLRF6*sin(alpha)), 1[N])MachNumber = 0.75PitchingMoment = torque_y()@DLRF6*2PitchingMomentatQuaterChord=-1*(-torque_y()@DLRF6*2/LiftForce-QuaterChordLocation)*LiftForce QuaterChordLocation = 0.1579 [m]ReferenceArea = 0.1453 [m^2]ReferenceChord = 0.1412 [m]ReferenceDensity = StaticPres/(287 [J K^-1 kg^-1]*Temp)ReferenceViscosity = ReferenceDensity*ReferenceChord*Vinlet/ReynoldsNumberReynoldsNumber = 3000000StaticPres = 119750 [Pa]Temp = 293.0 [K]Timestep = 1e-5 [s] +TimeStepRamp(atstep)*2e-4 [s]ToTemp = [...] http://blog.singinfo.com/?p=89
Sun, 24 Jun 2007 15:58:00 +0000
ANSYS：部分浸没固体模态分析
模态分析是常见的动态分析形式，通常我们都忽略周围流体对固体模态的影响。在一些结构中，固体会部分浸没在油中。液体可能会影响固体的频率。ANSYS Verification Manual里的VM177介绍了固体全部浸没在液体中时的模态分析方法。可以使用频率扫描或者模态分析的方法，找到共振频率。仿照这个例子，我做了一个圆柱形固体部分浸没在水中的模态分析的例子。APDL帖在下面。 具体的理论可以在流体分析部分的声学模块里找到。模块分析时，要使用unsymmetric方法来提取模。 我对结构分析不是很了解，不知道这样分析是否合理； 这只是一个简单的练习，网格质量和求解控制不是很好。欢迎指教。 我的朋友帮我比较了一个有水和没有水时候模态分析结果。其实对这个例子，流体对结果的影响很小。 提示：把文件保存为demo.inp，运行会比较方便。 finish/clear,start*ABBR,INPUT,/input,demo,inp/prep7et,1,fluid30 ! interface fluidet,2,solid45 ! solidet,3,fluid30,,1 ! non-interface fluidmp,ex,1,1mp,dens,1,1000mp,sonc,1,1482mp,ex,2,2e11mp,dens,2,7850mp,nuxy,2,0.3r,1,1CYL4, , ,0.01CYL4, , ,0.05AOVLAP,allVEXT,all, , ,0,0,0.05,,,,SMRT,4type,3MAT,1REAL,1VSWEEP,2type,2MAT,2VSWEEP,1EXTOPT,ESIZE,10,0VOFFST,2,0.05CSYS,1NSEL,S,LOC,X,0.01!NPLOTESLN,sESEL,R,TYPE,,3!EPLOTEMODIF,all,TYPE,1,allselCSYS,1NSEL,S,LOC,X,0.01NSEL,R,LOC,Z,0,0.05ESEL,S,TYPE,,1SF,ALL,FSIallselNSEL,S,LOC,X,0.05D,ALL,PRES,0.0NSEL,S,LOC,Z,0NSEL,R,LOC,x,0,0.01D,all, , , , , ,ALL, , , , ,Allselfinish/solutionantype,modalmodopt,unsym,10,-1mxpand,10solvefinish Technorati : ANSYS, 声学, 模态分析 http://blog.singinfo.com/?p=88
Sat, 23 Jun 2007 15:48:00 +0000
Workbench里怎么划分non-manifold geometry网格？
在Workbench里，manifold geometry是指对三维体而言，每个边只和两个面相连接。这是几何形状里最常见的一种拓朴结构。但是有些情况下，可能会生成non-manifold geometry，就是说有的边会和超过两个面相连接。 通常的基本形状如长方体、圆柱体都是manifold geometry。但是如下图所示的复杂情况，会有non-manifold geometry。左图最上面一条边和四个面相连接。要解决这个问题，可以把体切成两份，这样就把几何形状转换成manifold geometry。 对于二维情况而言，manifold geometry是指一个点只和两条边相连接；而non-manifold geometry是指一个点和超过两条边相连接。 在CFX-Mesh的帮助里明确指出CFX-Mesh是不支持non-manifold geometry。要划分non-manifold geometry生成CFD网格，除了上面提到的切割方法外，还可以试试用patch indepedent方法来划