GeoDict软件仿真流1.docx

1、GeoDict软件仿真流1GeoDict软件仿真流程简介中国北方车辆研究所 卢进军整理2012年4月GeoDict软件仿真过程包括下述几步：1.明确仿真条件：确定仿真模型边界条件，如流速、粒子直径和质量分布等，建模之前应该明确。了解仿真介质的化学特性，例如粒子与仿真介质的粘结力等。2.构建模型，即构建过滤介质或者过滤单元的几何模型，并对模型进展几何编辑。3.计算气体流场:获得初始压力降、透析度等过滤材料技术参数。4.加粒子(应用FilterDict模块)计算粒子的运动与分布：计算粒子通过过程与分布状态。5.计算模型过滤效率与压力降。6.通过计算结果改变模型参数，改良几何模型，重新进展25步。具

2、体流程如如下图1所示:图1 仿真流程图下面分步骤进展说明：1.明确仿真条件基于我们专业的特点，在进展过滤材料仿真计算之前，应该明确一些仿真模型边界条件，如进口或出口流速、粒子直径和质量分布等仿真输入，这样在仿真过程中才能够做到有的放矢。同时针对该软件的特点，作为微观仿真过滤材料的工具，我们在仿真计算前也应该适当了解仿真介质的化学特性，例如粒子与仿真介质的粘结力、粒子在这种材料的透过率等。2.构建模型在培训中，我们主要学习了GeoDict两种根本建模方法，它是由两种模型来实现的，第一种方法是直接纤维建模，通过FiberGeo模块来实现。其建模界面如如下图2所示：FiberGeo模块是建立在设计人

3、员对纤维材料形成过程与结构构造非常了解的根底上的，主要适用于纤维材料开发人员，而我们专业定位不在这一点上，故这种建模方法我们一般不采用。第二种建模方法是我们主要采用的方法，GeoDict有一个模块是模块是ImportModule，该模块可以导入图片文件，例如.jpg文件，如如下图3所示，而这里的图片文件就可以是我们用电镜拍摄的过滤介质单元的电镜照片，例如图4所示。图3 ImportModule中导入图片图4 电镜图片示例生成模型的厚度决定于电镜照片的X数X数少，厚度薄，在生成模型的过程中要调整纤维色彩与色彩的比照度，通过调整threshold值去除照片中的杂波，通过裁剪，最终生成更准确的纤维模

4、型。通过这种途径生成的模型可以进展材料特性分析像孔径分布，流电阻率(渗透性，压降)，热、电传导率，弹性，扩散能力，过滤效率，过滤粒子储存能力，等等。导入图片生成初始的几何模型后，就要对模型进展几何处理了，这时我们就要用到ProcessGeo模块，模块进入界面如如下图5所示：图5 ProcessGeo模块开始窗口ProcessGeo模块包含如下几个选项需要我们特别注意也是我们经常用到的：1Cleanse：去除附着在纤维上的小的连接局部，通常就是我们常说的杂质或杂波，例如在由一些电镜照片生成的模型中，threshold值设置之后，照片会产生一些杂波，可以应用该选项去除掉。2Dilate：字面的意思

5、是膨胀，在这个软件中，它是在已有的纤维外表包裹上一层纤维的意思，也可以反向包围那就是Erosion了，具体的例子可以参考软件手册。3Adding Binder：字面的意思是粘结到一起，在这个软件中，这个选项的作用是连接纤维，并将连接后的纤维重新生成为一个新的模型。Embed、Repeat、Mirror和Crop用来修饰微观机构的模型，对模型进展一些几何上的处理，下面分别说明：1Embed：字面的意思是把*嵌入到*上，在这个软件中，一般的操作时在模型的两侧添加两个无纤维的计算域，使模型延长，例如在粒子通过纤维模型的计算中，空的计算域在粒子计算时可以放粒子使用，加载的方向可以通过三维设定。2Rep

6、eat：顾名思义，就是在原模型的根底上再生成一个同样的模型，二者连接到一起形成一个新的模型。如如下图6所示：图6 Repeat选项的使用3Mirror：镜像的意思，其操作是在原模型的根底上按一定方向做对称镜像，二者连接成一体形成一个新的模型，如如下图7所示：图7 Mirror选项的使用4Crop：中文中的意思为修剪、剪短。在本软件中它的作用就是减去模型的一局部，比如前文中Repeat或者Mirror生产的一局部模型，可以通过Crop减去。也可以自己定义修剪模型的尺寸、方向等。3.计算流场，获得模型的流场特性。在已经获得模型的条件下，对模型进展特性计算。 如果模型是多孔介质模型，我们需要应用Po

7、roDict模块，如计算多孔介质模型的外表面积、孔径与渗透路径，这些参数可以决定多孔介质的根本特性。仿真计算后，结果被保存为*.txt格式文件，一些计算结果如几何孔径分布和孔隙率都可以通过PoroDexcel进展数据分析和图表显示。下面依据手册中的实例介绍PoroDict模块仿真过程：1根据手册中指示步骤创建如图8所示纤维模型：图8 PoroDict模块仿真演示模型2运行PoroDict模块，选择Geometrical Pore Size Distribution选项，在Geometrical PSD Parameters Bin Size中设置尺寸值，随后编辑gdt文件名，勾选write .

8、gdt for Visualization选项，点击Run按钮。3计算结果已经保存到初始设置的工作目录中，其中有些结果是.txt格式的文件，可以通过写字板或PoroDexcel进展数据分析和图表显示。但是在处理之前，应该进展Enable Macros的操作。4打开PoroDexcel，打开计算仿真的结果文件，在底部选择GeomPSD标签，并点击 Load Results File. 选择我们计算得到的，点击CreatChart还可以创建柱状图等，结果显示如图9所示。图9 PoroDict模块结果处理5同样的步骤还可以对材料进展Porosimetry Simulation、Percolation

9、 Path、surface area、Porosity determination、Chord Length Distribution等项目的计算。3.2FlowDict是用来计算流体通过纤维模型而得到纤维模型的各种性能参数的，如在特定流体介质穿透下的压力降、流阻和渗透率。在计算之前，首先要选择选用哪一种计算模型来计算，GeoDict提供如如下图10中的6种计算模型图10 FlowDict模块初始界面下面以第一种计算模型为例简要说明FlowDict模块计算过程：1构建如图8所示的纤维模型。2选择计算模型为Stokes (EJ)，设置计算求解器选项和边界条件如如下图11所示：图11 设置求解器选

10、项3设置仿真流场的物理参数，如流体样式和流体的具体参数等，设置界面如如下图12所示：图12 设置FlowDict计算的物理参数4点击Run进展计算，在Effective coefficients选项卡中得到了本次计算的结果，如选择Flow permeability，得到的计算结果值如如下图13显示：图13 FlowDict计算结果显示5FlowDict计算结果一般保存在临时文件中，如果采用其他计算模型分别计算的话，如此计算文件如此自动编号保存，如，保存在工作目录，需要对哪一组结果进展显示或编辑，读取该文件即可。FlowDict计算结果除了可以在选项卡中数值显示外，还可以通过流线或者流场等方式进

11、展3D图形显示，在显示之前可以设定显示的一些参数，如流线的数量、颜色等。具体显示样式如如下图14所示：图14 FlowDict计算结果的两种显示上流线，下流场4.对模型进展流场下加粒子仿真，以期得到某种过滤材料对指定粒子的过滤效率等数据这是我们专业仿真的重点所在，我要用较多的篇幅说明如何应用FilterDict对模型进展仿真计算，依旧沿用图8所建立的纤维模型，知识模块调整到FilterDict，界面如如下图15所示，首先我们注意到，仿真运算的第一个选项就是选择单通过仿真与多通过仿真。图15 FilterDict模块运算界面4.1 仿真计算求解器选项设置，点击Solver Options后的Ed

12、it按钮，会看到图16所示的求解器选项设置窗口，依据仿真条件设置后求解器后，再点击Process Parameters后的Edit按钮，进入仿真过程环境参数设置窗口如图17所示，设置好仿真过程环境参数之后，点击save保存，完成设置。进入粒子分布设置，点击Particle Distribution后的Edit选项，进入粒子分布的编辑窗口，如如下图18所示，在窗口中可以对粒子进展编辑和导入，系统默认的粒子设置时参照的SAEFine粉尘即SAE细灰，如果我们要在仿真中使用270目石英粉或者ISO的细灰和粗灰的话，我们应该获得这些粉尘的相应参数粒径、某种颗粒所占百分比、密度和粘性力模型，这四个参数中

13、，前三个是可以通过实验手段获取的，我们只需委托有检验粉尘能力的单位，就可以获取，第四个参数粘聚力模型是一种理论模型，其粘聚力系数也可以通过堆积角的测试方法获得，将这些参数导入到Particle Distribution中，就可完成仿真粒子的设置。图16 求解器选项设置窗口图17 仿真过程环境参数设置窗口图18 仿真粒子分布设置窗口在仿真之前还要注意的一点就是需要对模型进展一下处理，需要在模型的前端和后端加上入流和出流区域，以便计算时粒子有进入的空间和脱离的空间。这一操作可以在前面提到的ProcessGeo模块中实现，也可以直接通过增加一定方向和数量的voxel数量得以实现。在FilterDic

14、t中添加空白voxels的提示窗如如下图19所示如果不添加貌似是无法计算的，我没试过：图19 FilterDict中添加空白voxels的提示窗模型添加空白voxel之后的效果如如下图20所示：图20 模型添加空白voxel之后的效果图一切设置就绪之后，就可以点击Run进展运算了。4.2 计算完成之后，右下侧的主显示窗口就会显示带有初始粒子的模型，粒子就是我们设置的初始粒子。如图21所示：图21 仿真完成后显示模型在显示模型的上部就是粒子显示动画控制选项区Particle Animation Control，吧粒子选项由filtered调整到all，这样就不仅可以看到留在结构模型中的粒子，还可

15、以看到通过模型的粒子。如图22所示：图22 粒子显示动画控制选项区Particle Animation Control4.3 点击Play可以动画显示粒子穿过纤维模型的动态过程。在此过程中可以调整播放速度，暂停和保存输出照片。过程如图23所示：图23 粒子穿过纤维模型动态过程4.4 例子1：粒子穿越双纤维模型构造如下列图的双纤维模型如图24所示具体建模过程可以参考GeoDict用户手册，建模完成之后，切换到FilterDict模块，依次设置Solver Options、Process Parameters和Particle Distribution等选项，其中在Process Paramete

16、rs中设置Number of batches为1，设置Particles per Batch为5000，设置速度为10，运行仿真。图24 仿真前后的双纤维模型可以通过不选择Structure选项隐藏模型在菜单中点击View-3D Objects，然后运行Play观看粒子运行动画，如图25所示：图25 模型隐藏状态下的粒子运动过程动画显示点击Particle Distributions Edit按钮，进入Edit Particle Distribution窗口，可以手动编辑不同的粒子分布，如图26所示，可以分别调整粒子的粒径与百分比分布，通过Remove Last Row按钮可以去除系统默认的粒

17、子分布设置，自己定义的粒子分布可以保存为*.pdd2文件，同样不同的粒子分布的*.pdd2文件还可以加载到Edit Particle Distribution窗口中。图26 粒子编辑串口点击Process Parameters Edit，编辑过程参数，粒子束数设为3，速度设为20m/s,温度设为30介质为空气，每束粒子为5000，设置完毕。温度会影响到流体的年度和密度两个参数，在过滤过程的计算中，预设的流体粘度和流体密度如空气和水都是有文件获取的，*.gts文件可以通过很多编辑器记事本、写字板、word等直接进展编辑，所以，在很广泛的温度X围内，流体的粘度和密度，是可以定义的。开始计算，同样，

18、计算之前要在模型的两侧添加10个单元作为入流区域和出流区域。计算完毕，几种结果显示形式如图27所示：图27 几种粒子显示状态4.5 过滤过程流场的可视化为了研究过滤过程中流场的一些特性变化，我们需要可视化模型过滤过程中的整个流场，首先我们需要先加载结果文件，在FilterDict选项中点击LoadResults，在打开的窗口中选择两个文件。在菜单栏选择File Show Flow Field.打开Load flow solution file窗口，选择在计算中生成的文件，在视图区域就会看到以绿色箭头显示的流场矢量场，把粒子显示设置为off，隐藏粒子，也可以隐藏纤维模型，显示样式如如下图28所示

19、：图28 矢量显示流场4.6 过滤过程中流场的计算在FiberGeo中构造如下列图模型：切换到FilterDict模块，选择默认的SAEfine.pdd2粉尘作为粒子，设置仿真参数：Number of Batches为3，PressureDrop Increase为200 Pa，Particles per Batch 为 5000，速度为10m/s，温度设为20。然后在FilterDict中计算，计算结果如如下图显示：运行FiterDexcel可以将计算结果进展图表显示，参照以下步骤：The FilterDexcel spreadsheet that opens is divided into

20、 sections, observable by scrolling down:a. In 1. Select Results Files, click Select FiltrationLogFile and choose the file. Click Select ParticleLogFile and choose the particlesLogFile.gfp file from thechosen project directory (or folder). The file names and the paths to them in theproject directory

21、(or folder) appear in the yellow fields. Click Read Files.b. In 2. Visualization of Filtration Efficiency, check that the Filtration Batch is set to1. Choose a Chart Title and a Spreadsheet Name Base or keep the default names(Filtration Efficiency).Tutorial Sessions- 126 - GeoDict User Guidec. In 3.

22、 Visualization of Layerwise Filtration Efficiency (FE), check that the FiltrationBatch is set to 1. Choose a Chart Title and a Spreadsheet Name Base or keep thedefault names (Layerwise FE).d. In 4. Visualization of Particle Deposition, enter 32 for the Number of Layers andset the Largest Particle Di

23、ameter to 22. Choose a Chart Title and a SpreadsheetName Base or keep the default names (Particle Deposition).e. In 5. Visualization of Pressure Drop Evolution, choose a Chart Title and aSpreadsheet Name Base or keep the default names (Pressure Drop).f. For all these sections, click Create Chart sequentially.For all of them, a chart is plotted and a results spreadsheet page is created. Switch betweenthem and the Master page at your convenience by clicking on the tabs at the bottom of theFilterDexcel spreadsheet.Filtration EfficiencyLayerwise FEParticle DepositionPressure Drop