当前位置:文档之家› 计算流体力学ICEM CFD 网格生成基础教程0204192301

计算流体力学ICEM CFD 网格生成基础教程0204192301

计算流体力学ICEM CFD 网格生成基础教程0204192301
计算流体力学ICEM CFD 网格生成基础教程0204192301

第一章介绍

ICEM CFD 工程

Tutorials目录中每个工程是一个次级子目录。每个工程的目录下有下列子目录:import, parts, domains, mesh, 和transfer。他们分别代表:

? import/: 要导入到ICEMCFD中的集合模型交换文件,比如igs,STL等;

? parts/: CAD模型

? domains/: 非结构六面体网格文件(hex.unstruct), 结构六面体网格分区文件(domain.n), 非结构四面体网格文件(cut_domain.1)

? mesh/: 边界条件文件(family_boco, boco),结构网格的拓扑定义文件(family_topo, topo_mulcad_out), 和Tetin几何文件(tetin1).

? transfer/: 求解器输入文件(star.elem), 用于Mom3d.的分析数据

mesh目录中Tetin文件代表将要划分网格的几何体。包含B-spline曲面定义和曲线信息,以及分组定义

Replay 文件是六面体网格划分的分块的脚本

鼠标和键盘操作

鼠标或键盘操作功能

鼠标左键点击和拖动旋转模型

鼠标中键点击和拖动平移模型

鼠标右键点击和上下拖动缩放模型

鼠标右键点击和左右拖动绕屏幕Z轴旋转模型

F9 按住F9,然后点击任意鼠标键进行操作的时候进行模型运动

F10 按F10 紧急图象Reset

第二章ICEM CFD Mesh Editor界面

The Mesh Editor, 创建修改网格的集成环境,包含三个窗口

? The ICEM CFD 主窗口

? 显示窗口

? The ICEM CFD 消息窗口

主窗口

主窗口中除了图形显示区域,外,还有6个radio按钮:File, Geometry, Meshing, Edit Mesh and Output. The File Menu

The File menu 包含

? Open, Save, Save as, Close, Quit, Project dir, Tetin file,

Domain file, B.C file, Import geo, Export geo, Options, Utilities,

Scripting, Annotations, Import mesh, DDN part.

The Geometry Menu

The Geometry menu 模型修补和编辑,边界条件的设置,调用ICEM CFD DDN。它包含

? DDN tools, Bound conds, Repair, Utilities, Global setup.

模型编辑模式,由一排彩色radiobuttons控制

? Surface, Curve, Point, Material, Density, Loop.

它们又分别具有次级菜单

? Create, Delete, Modify, Mesh params, Change family,

Copy/move.

The Meshing Menu

The meshing menu 包含ICEM CFD各种网格划分模块。只有用户拥有楼阁模块的license,他才可以使用这个模块。目前ICEMCFD提供下列模块:

? Hexa, Tetra, Global, Prism, Quad, AutoHexa, Mulcad/Padamm,

P-Cube

按不同的按钮会调用不同的模块。

The Edit Mesh Menu

The Edit Mesh menu包含必要的网格编辑功能, 粗化、平滑合并等。具体的操作有:

? Copy/move, Smooth, Refine, Coarsen, Merge, Extrude,

Diagnostics, Uncouple, Bandwidth, Change type, Change

family, Utilities, Edit nodes, Edit elements, Edit edges, Edit

Blocks, Edit subfaces, Repair

The Output Menu

The Output menu针对不同求解器进行边界条件的设置。此外用户可以调用Mom3d 和Visual3 进行网格自适

应和可视化。可使用的功能有:

? Select solver, Bound conds, Solver params, Solver input,

Run solver, Visual3, Mom3d, RAMM-ICE

Note: The ICEM CFD 后处理模块Visual3, 网格优化模块Mom3d, 以及与RAMM-ICE的内燃机网格接口必须有另外的license

The Utilities Cluster

这些功能都位于主窗口的右上角:

? Help: connects the user to the ICEM CFD on-line help

? Orient: 控制显示坐标

? View: 定制显示属性

? Undo: 操作反悔

? Redo: 撤销反悔

? Print: 打印

? Shell: 调用一个X-Term命令窗口

显示窗口

The display window, 处于屏幕右边, 允许控制按照family, geometric entity, element

type and user-defined subsets来显示.

Important: Since some functions are performed only on the entities

shown, the Display window is a very important feature to use when

isolating the particular entities to be modified.

Families

Color-coded and customizable, display by family, as defined in the

mesher interface, is available in this window.

Geometry

When this button is toggled on, the user has control over the

display of geometric entities. The Opts button beside for each

entity type allows control over the display characteristics.

Note: Based on the domain file chosen, this toggle button will

either be labeled “unstructured” or “structured,” or be absent

altogether if no domain file is specified.

? Structured or Unstructured: Aside from display of the mesh

based on element type, for even greater power, ICEM CFD

allows the user to create custom subsets; the careful display of

which allows a great deal of control in the mesh operations.

? Display Subset: This toggle button allows you to define subset

of your model based on a set of numerous attributes, including

element type, screen selection, coordinate position, etc.

第三章ICEM CFD Tetra四面体网格划分器

这是Tetra生成的棱拄和四面体混合网格,包含55万四面

体网格和12层33万棱拄网格

介绍

Tetra采用8叉树算法来对体积进行四面体填充并生成表面网格。用户必须事先规定一些关键的点和曲线。Tetra 具有强大的网格平滑算法,以及局部适应性加密和粗化算法。

对于复杂模型,ICEM CFD Tetra具有如下优点:

? 基于8叉树算法的网格生成

? 快速模型set-up

? 网格与表面拓扑独立

? 无需表面的三角形划分

? 可以直接从CAD模型和STL数据进行网格生成

? 对CAD surfaces and/or STL Representation定义网格尺寸

? 控制体积内部的网格尺寸

? 四面体的节点和曲线与事先的规定匹配

? 采用Natural size 单独的决定几何特征上的四面体网格尺寸

? 体积网格和表面网格的平滑、节点合并和边交换

? 四面体网格能够合并到混合网格中,并实施平滑操作

? 单独区域的粗化

? Enforcement of mesh periodicity, both rotational and translational

? 表面网格编辑和诊断工具

? 局部细化和粗化

? 为多种材料提供一个统一的网格

? 快速算法: 1500 cells/second

ICEM CFD Tetra的输入

ICEM CFD Tetra的输入方法有:

? Sets of B-Spline curves and trimmed B-Spline surfaces with

prescribed points

? Triangular surface meshes as geometry definition

B-Spline曲线和曲面当B-spline曲线和曲面输入后,采用线段和三角形对曲面和曲线进

行近似,并在预先规定的点上设置顶点。

B-Spline曲线允许Tetra 处理表面上的间断。如果在表面的边

界上没有定义曲线,Tetra划分的三角形会自由的越过间断。

而预先规定的点会使得它认出曲线上尖锐的拐弯。Tetra

中有工具来自动的在尖锐的特征上获取点和曲线。

三角形表面网格对于三角形表面网格,关键点和曲线能够自动的被识别。虽然Tetra

生成的网格上的节点不完全和原始的网格重合,但是它会符合模

型的形状。这对于从别的网格数据或者立体扫描数据中导入几何

模型是很有用处的。

ICEM CFD Tetra中的智能几何

ICEM CFD的CAD接口,保留有CAD几何模型的参数化描述,网格可以在修改过的几何模型上重新生成。在CAD 中选中被导入的模型带有附加的信息,它们与主几何模型一起存储,几何模型的参数改变后,用户要重新生成网格只需简单的File > Save,就可以立即进行非结构四面体网格重新计算。

八叉树算法

Tetra网格生成是基于如下的空间划分算法:这种算法需要的区域保证必要的网格密度,但是为了快速计算尽量采用大的单元。

1. 在几何模型的曲线和表面上规定网格尺寸

2. 构造一个初始单元来包围整个几何模型

3. 单元被不断细分来达到最大网格尺寸(每个维的尺寸按照1/2分割,对于三维就是1/8)

4. 均一化网格来消除悬挂网格现象

5.构造出最初的最大尺寸单元网格来包围整个模型

6.节点调整以匹配几何模型形状

7.剔除材料外的单元

8.进一步细分单元以满足规定的网格尺寸要求

9. 通过节点的合并、移动、交换和删除进行网格平滑

Family设置, Material Points, 预先规定Points

在ICEMCFD中可以对几何模型中的各种几何元素进行分组,形成不通的families。这样允许用户对不通的families 进行不同的参数设置。除了对不通families设置边界条件,用户还可以对每个family定义控制网格尺寸的参数:maximum length, initial height, and height ratio。此外还可以对单独的曲线和表面定义网格尺寸参数。这样通过预先规定的网格尺寸参数,用户可以控制四面体节点和关键区域的边的位置。

如上面的八叉树算法已经描述的,当网格划分器调整网格节点以匹配几何模型外形时,它首先匹配预先规定点和曲线。

在单元剔除步骤,Tetra要求在每个需要进行分析的不同材料内定义一个material point。对于CFD,material point可以来定义流体区域;对于结构有限元分析,可以来定义固体区域;或者对于耦合传热问题,用来定义两种区域。

ICEM CFD Tetra重要特性

Natural Size如果对一个family定义的最大四面体尺寸比某个成员几何元素大,用户必须使用

natural size限制。The user can specify a natural size

that is proportional to the reference tetra size; it should be assigned

a value that is slightly smaller than the smallest gap in the model,

so that the mesher will further subdivide the tetrahedra to match

this geometric feature.

The natural size limit 将是整个网格中的最小尺寸。Only if the user defines a

smaller max size on a geometry entity will Tetra continue to subdivide until meeting

the maximum size request. The effect of the natural size limit is a

geometry-based adaptation of the mesh.

Tetrahedral Mesh Smoother the tetrahedral smoother计算每个单元的质量。然后根据用户指定的质量要求进行平滑。

Smoother修改所有不符合标准的单元——节点被移动、合并,边被交换,或者单元

被删除。如此重复直到满足要求或者达到一定的迭代次数。

有时候用户不愿对某个family进行网格平滑,ICEM CFD Tetra 提供fix_families

功能来指定某个family保持不变。此外,另外一个相似的功能fix_prisms,用来排除

棱拄单元的平滑。

Tetrahedral Mesh Coarsener Coarsening和smoother的过程相似。此外也可以采用fix_families来排除某个family的

coarsen操作。

如果coarsen操作中,size_check选项打开,得到的网格在所有几何元素上会满足预设的

尺寸标准。而且Tetra包含了完全的投影projection和smoothing工具tools,以及单元创

建、删除、分裂、交换和uniform enforcement of orientation.

Triangular

Surface Mesh

Smoother

In the interest of minimizing grid points, the coarsener reduces the

number of triangles in a mesh by merging triangles. This operation

is based on the maximum deviation of the resultant triangle center

from the surface, the aspect ratio of the merged triangle, and the

maximum size of the merged triangle.

Triangular

Surface Editing

Tools

ICEM CFD Tetra 提供交互式网格编辑功能,可以实现节点在几何模型表面的移动、

合并或者删除。三角形单元可以被单独的分割和标记以不同的名称,并进行质量检

查以进行评价。

表面网格诊断工具允许用户快速填充表面网格中的“洞”,查找重叠的三角形和non-

manifold顶点以及detection of single/multiple edge and duplicate elements.

Prism Layer Offset对于支持非结构四面体和棱拄混合网格的求解器来说,ICEM CFD Tetra 的一个独

特功能Prisms 能够在壁面附近产生偏移的棱拄层网格。

Mesh Periodicity采用周期性网格ICEM CFD Tetra 网格能够很好的适应叶轮机问题的分析。

第三章四面体网格划分示范1:三维管道回合

figure 3.1.1

本示范介绍的操作有:

Geometry

? 设置网格参数Mesh params >Model 和Mesh params > Families.

Meshing menu

? 自动网格生成和平滑Tetra > Batch from CAD

Edit mesh menu

? 检查网格质量Diagnostics.

打开工程

打开File > Open,选择Select project,选择3DpipeJunct project,再选择tetin1。如果出现一个提示窗口要求加载一个domain,则选择Cancel 来开始一个新的模型。

生成四面体网格

定义family

1 如果families还未定义,那么首先需要对几何模型的不同区域分类以创建family。最初整个几何模型被设置为

一个family,名称为GEOM。我们将按照figure 3.1.1所示来修改family设置。表面family有INL, OUT, CYL1, CYL2, and SYM。

2 选择Geometry > Surface> Change family.来修改family设置

3 用鼠标左键选择某个表面,用鼠标中键确认选择,然后会出现一个family设置窗口。或者为所选择的面选择一个已经存在的family,或者在New family中输入一个新的名称来为这个(些)面创建一个family

4 一旦所有的surface被设置重新命名,采用相似的操作继续对Curves, Points, 和Materials 进行设置。

5 在Geometry > Curve> Change family几何模型中所有的曲线将被设置为CUR。

用户可以按住鼠标左键拖动进行box-select。

6 所有的点放在PTS family中。注意将光标放在主窗口中,按"v"键,则可以选择

所有点元素。

定义Material点

1 模型的material将被设置为LIVE family。它处于圆柱的内部。

2 选择Geometry > Material > Create > Screen midpoint.

3 在Select family 窗口中, 输入LIVE 作为New family 名称,按Accept 接受

4 用鼠标左键在圆柱上选择两个相对的顶点,然后按鼠标中键确认。LIVE 应该

在几何模型的内部。旋转模型以确认它确实在体积内部。

5 要删除空的families,在Display菜单中选择Families > Edit,点击Delete

unused。这样列表中的无用family会自动清除。按Accept关闭这个窗口。

设置网格参数

1 要确认Surfaces 和Curves 都已经被toggled on, 在display窗口的分别按Surfaces 和Curves的Opts > Tetra sizes,四面体形状的图标就会显示在几何模型上。

2 这时用户可以用Mesh params > Model 来修改整个模型的网格尺寸(右图)。将Reference size设置为5,

size

不同Reference size的效果

3

注意改变Reference size 并按Apply 会改变几何模型上的四面体图标的尺寸。

Surface>Mesh params > Families。用户就可以从Select

family窗口中选择family(如图)。本示范中选择CYL2 然后

按Accept。CYL2 对应的表面就会以加亮显示。

5 在Surface family mesh params 窗口中(如图),对Size

指定为0.5,其他参数为0,因为Tetra 生成网格时不用这些

参数,只有Hexa用。选择Show size params using

refsize 会使得按下Apply后参数调整可见。选择Dismiss 关

闭窗口。

6 当对网格尺寸设置满意后,按File > Save 来存储。

Running the Tetra Mesh

1. 选择Meshing > Tetra > from CAD/Surface Mesh

来生成四面体网格。Tetra parameters 窗口将会

出现(如右图);

2. 按照图中设置后,按Mesh进行网格划分。

Message Window将会报告网格划分过程信息。

译者注:这个步骤和PDF文档中所述不相同,原文中的菜单和Tetra parameters窗口与4.2版不太一样。原文和原图如下:

Choose Meshing > Tetra > Batch from CAD to generate the tetrahedral mesh. The Tetra > Interactive from CAD mode indicates that Tetra will not automatically save the volume mesh (Write volume) or surface mesh (Write surface). The Tetra parameters window will appear (refer to Figure 3.1.4):

最后显示的网格为:

Saving the Mesh

用File > Save >Domain file 存盘

第五章Sphere Cube

一个半球体,内含一个立方体空腔。图中以标明表面的family名称。

1.打开工程

在Select project 窗口,选择SphereCube工程后,再打开tetin1文件。

3. 生成四面体网格

Family设置

1. 与前面的示范相同,最初所有的表面属于一个family——GEOM。用户需要按照上图对表面重新创建family

(SPHERE, CUBE, SYM) 。

2. 在Geometry > Surface> Change family菜单中修改family设置。

3. 用鼠标左键选择需要的表面,然后中键确认选择。在Select family 窗口出现后,或者选择已经存在的

family,或者创建新的family。

4. 本示范中,对于所有的点和曲线,设为一个family, GEOM。

定义Material点

同上个示范一样,material点被设置为单独的family,LIVE。采用Create > Screen midpoint功能设置material点。用鼠标左键在半球上选择两个点,按中键确认,LIVE就会出现在模型上。

Reassigning Mesh Parameters

1. 用户必须在网格生成之前设置网格参数。在Geometry 分类中,

按Mesh Params > Model 来进入Model Mesh Params 窗口。

2. 注意: 为了方便,建议打开Reference size旁边的display选项来显示

标记为ref 的四面体图标。从Display菜单中,打开Surfaces > Opts

> Tetra sizes和Curves > Opts > Tetra sizes,以显示代表网格尺寸

的图标。

3. 回到Model mesh params 窗口,输入如图参数;

4. 在Mesh params > Families中可以对不通的family设置不同的网格

参数。对于SYM,按照下图设置。

5.存盘

6.参考第四章选择Meshing > Tetra > from

CAD/Surface Mesh 来生成四面体网格。在Tetra

parameters 窗口做设置后,按Mesh进行网格划分。

Message Window将会报告网格划分过程信息。

译者注:原PDF文档所述与此不同。ICEMCFD4.2版的

Tetra parameters 窗口中无法关闭Smoother。因此在这

里设置Smoother不要做很好的平滑。以便后面进行平滑

操作的对比更明显。

7. 网格划分好后存盘。

中国石油大学流体力学实验报告

中国石油大学(流体力学)实验报告 实验日期:2012-2-15 成绩: 班级:学号:姓名:教师: 同组者: 实验一、流体静力学实验 一、实验目的 1.掌握用液式测压计测量流体静压强的技能; 2.验证不可压缩流体静力学基本方程,加深对位置水头、压力水头和测压管水头的理解; 3. 观察真空度(负压)的产生过程,进一步加深对真空度的理解; 4.测定油的相对密度; 5.通过对诸多流体静力学现象的实验分析,进一步提高解决实际问题的能力。 二、实验装置 1、在图1-1-1下方的横线上正确填写实验装置各部分的名称 本实验的装置如图所示。 1.测压管; 2.带标尺的测压管; 3.连通管; 4.通气阀; 5.加压打气球; 6.真空测压管; 7 截止阀.;8. U形测压管;9.油柱;

10.水柱;11.减压放气阀 图1-1-1流体静力学实验装置图 2、说明 1.所有测管液面标高均以测压管2标尺零读数为基准; 2.仪器铭牌所注B ?、C ?、D ?系测点B 、C 、D 标高;若同时取标尺零点作为静力学基本方程的基准,则B ?、C ?、D ?亦为B z 、C z 、D z ; 3.本仪器中所有阀门旋柄均以顺管轴线为开。 三、实验原理在横线上正确写出以下公式 1.在重力作用下不可压缩流体静力学基本方程 形式之一: const p =+ γ z (1-1-1a ) 形式之二: h p p γ+=0(1-1b ) 式中z ——被测点在基准面以上的位置高度; p ——被测点的静水压强,用相对压强表示,以下同; 0p ——水箱中液面的表面压强; γ——液体重度; h ——被测点的液体深度。 2. 油密度测量原理 当U 型管中水面与油水界面齐平(图1-1-2),取其顶面为等压面,有 01w 1o p h H γγ==(1-1-2) 另当U 型管中水面和油面齐平(图1-1-3),取其油水界面为等压面,则有 02w o p H H γγ+= 即 02w 2o w p h H H γγγ=-=-(1-1-3)

高等流体力学重点

1.流体的连续介质模型:研究流体的宏观运动,在远远大于分子运动尺度的范围里考察流体运动,而不考虑个别分子的行为,因此我们可以把流体视为连续介质。 它有如下性质: (1)流体是连续分布的物质,它可以无限分割为具有均布质量的宏观微元体。 (2)不发生化学反应和离解等非平衡热力学过程的运动流体中,微元体内流体状态服 从热力学关系 (3)除了特殊面外,流体的力学和热力学状态参数在时空中是连续分布的,并且通常 认为是无限可微的 2.应力:有限体的微元面积上单位面积的表面力称为表面力的局部强度,又称为应力,定义如下:=n T A F A δδδlim 0→ 3.流体的界面性质:微元界面两侧的流体的速度和温度相等,应力向量的大小相等.方向相反或应力分量相等。 4.流体具有易流行和压缩性。 5.应力张量具有对称性。 6.欧拉描述法:在任意指定的时间逐点描绘当地的运动特征量(如速度、加速度)及其它的物理量的分布(如压力、密度等)。 7.拉格朗日描述法:从某个时刻开始跟踪质点的位置、速度、加速度和物理参数的变化,这种方法是离散质点的运动描述法称为拉格朗日描述法。 8.流线:速度场的向量线,该曲线上的任意一点的切向量与当地的的速度向量重合。 迹线:流体质点点的运动迹象。 差别:迹线是同一质点在不同时刻的位移曲线。 流线是同一时刻、不同质点连接起来的速度场向量线。 流线微分方程:ω dz v dy u dx == 迹线微分方程:t x U i i ??= 9.质点加速度:质点速度向量随时间的变化率。 U U t U a )(??+??= 质点加速度=速度的局部导数+速度的迁移导数。 物理量的质点导数=物理量的局部导数+物理量的对流导数。

流体力学实验报告

流体力学 实验指导书与报告 静力学实验 雷诺实验 中国矿业大学能源与动力实验中心

学生实验守则 一、学生进入实验室必须遵守实验室规章制度,遵守课堂纪律,衣着整洁,保持安静,不得迟到早退,严禁喧哗、吸烟、吃零食和随地吐痰。如有违犯,指导教师有权停止基实验。 二、实验课前,要认真阅读教材,作好实验预习,根据不同科目要求写出预习报告,明确实验目的、要求和注意事项。 三、实验课上必须专心听讲,服从指导教师的安排和指导,遵守操作规程,认真操作,正确读数,不得草率敷衍,拼凑数据。 四、预习报告和实验报告必须独自完成,不得互相抄袭。 五、因故缺课的学生,可向指导教师申请一次补做机会,不补做的,该试验以零分计算,作为总成绩的一部分,累计三次者,该课实验以不及格论处,不能参加该门课程的考试。 六、在使用大型精密仪器设备前,必须接受技术培训,经考核合格后方可使用,使用中要严格遵守操作规程,并详细填写使用记录。 七、爱护仪器设备,不准动用与本实验无关的仪器设备。要节约水、电、试剂药品、元器件、材料等。如发生仪器、设备损坏要及时向指导教师报告,属责任事故的,应按有关文件规定赔偿。 八、注意实验安全,遵守安全规定,防止人身和仪器设备事故发生。一旦发生事故,要立即向指导教师报告,采取正确的应急措施,防止事故扩大,保护人身安全和财产安全。重大事故要同时保护好现场,迅速向有关部门报告,事故后尽快写出书面报告交上级有关部门,不得隐瞒事实真相。 九、试验完毕要做好整理工作,将试剂、药品、工具、材料及公用仪器等放回原处。洗刷器皿,清扫试验场地,切断电源、气源、水源,经指导教师检查合格后方可离开。 十、各类实验室可根据自身特点,制定出切实可行的实验守则,报经系(院)主管领导同意后执行,并送实验室管理科备案。 1984年5月制定 2014年4月再修订 中国矿业大学能源与动力实验中心

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、观察液体流动时的层流与紊流现象。区分两种不同流态的特征,搞清两种流态产生的条件。分析圆管流态转化的规律,加深对雷诺数的理解。 2、测定颜色水在管中的不同状态下的雷诺数及沿程水头损失。绘制沿程水头损失与断面平均流速的关系曲线,验证不同流态下沿程水头损失的规律就是不同的。进一步掌握层流、紊流两种流态的运动学特性与动力学特性。 3、通过对颜色水在管中的不同状态的分析,加深对管流不同流态的了解。学习古典流体力学中应用无量纲参数进行实验研究的方法,并了解其实用意义。 (二)、实验原理 1、液体在运动时,存在着两种根本不同的流动状态。当液体流速较小时,惯性力较小,粘滞力对质点起控制作用,使各流层的液体质点互不混杂,液流呈层流运动。当液体流速逐渐增大,质点惯性力也逐渐增大,粘滞力对质点的控制逐渐减弱,当流速达到一定程度时,各流层

的液体形成涡体并能脱离原流层,液流质点即互相混杂,液流呈紊流运动。这种从层流到紊流的运动状态,反应了液流内部结构从量变到质变的一个变化过程。 液体运动的层流与紊流两种型态,首先由英国物理学家雷诺进行了定性与定量的证实,并根据研究结果,提出液流型态可用下列无量纲数来判断: Re=Vd/ν Re 称为雷诺数。液流型态开始变化时的雷诺数叫做临界雷诺数。 在雷诺实验装置中,通过有色液体的质点运动,可以将两种流态的根本区别清晰地反映出来。在层流中,有色液体与水互不混惨,呈直线运动状态,在紊流中,有大小不等的涡体振荡于各流层之间,有色液体与水混掺。 2、在如图所示的实验设备图中,取1-1,1-2两断面,由恒定总流的能量方程知: f 2222221111h g 2V a p z g 2V a p z ++γ+=+γ+ 因为管径不变V 1=V 2 ∴=γ +-γ+=)p z ()p z (h 2211f △h 所以,压差计两测压管水面高差△h 即为1-1与1-2两断面间的沿程水头损失,用重量法或体积浊测出流量,并由实测的流量值求得断面平均流速A Q V =,作为lgh f 与lgv 关系曲线,如下图所示,曲线上EC 段与BD 段均可用直线关系式表示,由斜截式方程得: lgh f =lgk+mlgv lgh f =lgkv m h f =kv m m 为直线的斜率

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 生成网格的过程

计算流体力学教案

计算流体力学教案 Teaching plan of computational fluid mechanics

计算流体力学教案 前言:本文档根据题材书写内容要求展开,具有实践指导意义,适用于组织或个人。便于学习和使用,本文档下载后内容可按需编辑修改及打印。 一、流体地基本特征 1.物质地三态 在地球上,物质存在地主要形式有:固体、液体和气体。 流体和固体地区别:从力学分析地意义上看,在于它们对外力抵抗地能力不同。 固体:既能承受压力,也能承受拉力与抵抗拉伸变形。 流体:只能承受压力,一般不能承受拉力与抵抗拉伸变形。 液体和气体地区别:气体易于压缩;而液体难于压缩; 液体有一定地体积,存在一个自由液面;气体能充满任意形状地容器,无一定地体积,不存在自由液面。 液体和气体地共同点:两者均具有易流动性,即在任何 微小切应力作用下都会发生变形或流动,故二者统称为流体。 2.流体地连续介质模型

微观:流体是由大量做无规则运动地分子组成地,分子之间存在空隙,但在标准状况下,1cm3液体中含有3.3×1022个左右地分子,相邻分子间地距离约为3.1×10-8cm。1cm3气体中含有2.7×1019个左右地分子,相邻分子间地距离约为3.2×10-7cm。 宏观:考虑宏观特性,在流动空间和时间上所采用地一切特征尺度和特征时间都比分子距离和分子碰撞时间大得多。 (1)概念 连续介质(continuum/continuous medium):质点连续充满所占空间地流体或固体。 连续介质模型(continuum continuous medium model):把流体视为没有间隙地充满它所占据地整个空间地一种连续介质,且其所有地物理量都是空间坐标和时间地连续函数地一种假设模型:u =u(t,x,y,z)。 (2)优点 排除了分子运动地复杂性。物理量作为时空连续函数,则可以利用连续函数这一数学工具来研究问题。 3.流体地分类

第二章计算流体力学的基本知识

第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中,所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式,最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂,在考虑粘性作用时更是如此,如果不靠计算机,就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30~40 年代,对于复杂而又特别重要的流体力学问题,曾组织过人力用几个月甚至几年的时间做数值计算,比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展,计算机的不断进步,以及流体力学各种计算方法的发明,使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性,这又促进了流体力学计算方法的发展,并形成了"计算流体力学" 。 从20 世纪60 年代起,在飞行器和其他涉及流体运动的课题中,经常采用电子计算机做数值模拟,这可以和物理实验相辅相成。数值模拟和实验模拟相互配合,使科学技术的研究和工程设计的速度加快,并节省开支。数值计算方法最近发展很快,其重要性与日俱增。 自然界存在着大量复杂的流动现象,随着人类认识的深入,人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前,流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题,无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现,从而催生了计算流体力

流体力学-伯努利方程实验报告

中国石油大学(华东)工程流体力学实验报告 实验日期:2014.12.11成绩: 班级:石工12-09学号:12021409姓名:陈相君教师:李成华 同组者:魏晓彤,刘海飞 实验二、能量方程(伯诺利方程)实验 一、实验目的 1.验证实际流体稳定流的能量方程; 2.通过对诸多动水水力现象的实验分析,理解能量转换特性; 3.掌握流速、流量、压强等水力要素的实验量测技能。 二、实验装置 本实验的装置如图2-1所示。 图2-1 自循环伯诺利方程实验装置 1.自循环供水器; 2.实验台; 3.可控硅无极调速器;4溢流板;5.稳水孔板; 6.恒压水箱; 7.测压机;8滑动测量尺;9.测压管;10.试验管道; 11.测压点;12皮托管;13.试验流量调节阀 说明 本仪器测压管有两种: (1)皮托管测压管(表2-1中标﹡的测压管),用以测读皮托管探头对准点的总水头; (2)普通测压管(表2-1未标﹡者),用以定量量测测压管水头。 实验流量用阀13调节,流量由调节阀13测量。

三、实验原理 在实验管路中沿管内水流方向取n 个过水断面。可以列出进口断面(1)至另一断面(i )的能量方程式(i =2,3,…,n ) i w i i i i h g v p z g p z -++ + =+ + 1222 2 111 1αγυαγ 取12n 1a a a ==???==,选好基准面,从已设置的各断面的测压管中读出 z+p/r 值,测 出透过管路的流量,即可计算出断面平均流速,从而即可得到各断面测压管水头和总水头。 四、实验要求 1.记录有关常数实验装置编号 No._4____ 均匀段1d = 1.40-210m ?;缩管段2d =1.01-210m ?;扩管段3d =2.00-2 10m ?; 水箱液面高程0?= 47.6-2 10m ?;上管道轴线高程z ?=19 -2 10m ? (基准面选在标尺的零点上) 2.量测(p z γ + )并记入表2-2。 注:i i i p h z γ =+ 为测压管水头,单位:-2 10m ,i 为测点编号。 3.计算流速水头和总水头。

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技术太重要了。

室外风环境模拟计算报告123

新项目 室外风环境模拟计算报告 计算软件:风模拟分析软件PKPM-CFD 开发单位:中国建筑科学研究院 建研科技股份 合作单位:Software Cradle Co., Ltd. 韵能建筑科技 应用版本:Ver1.00 2015.10.19

室外风环境模拟分析报告 项目名称:新项目 项目地址: 建设单位: 设计单位: 参与单位: 规标准参考依据: 1、《绿色建筑评价标准》(GB/T 50378-2014) 2、《民用建筑设计通则》(GB 50352-2005) 3、《绿色建筑评价技术细则》

一、项目概述 1.1计算模型概况 1.2建筑物概况 图 1 建筑群平面图,红线建筑为目标建筑

二、指标要求 针对室外风环境评价依据为《绿色建筑评价标准》(GB/T 50378-2014)中有关室外风环境的条目要求。 2.1规的评价要求 《绿色建筑评价标准》(GB/T 50378-2014)中有关室外风环境的具体要求如下: 4.2.6 场地风环境有利于室外行走、活动舒适和建筑的自然通风。评分规则如下: 1 冬季典型风速和风向条件下,建筑物周围人行区风速低于5m/s,且室外风速放大系数小于2,得2分;除迎风第一排建筑外,建筑迎风面与背风面表面风压差不超过5Pa,再得1分。 2 过渡季、夏季典型风速和风向条件下,场地人活动区不出现涡旋或无风区,得2分;50%以上可开启外窗室外表面的风压差大于0.5Pa,得1分。 2.2模拟条件设置要求 1、室外风环境模拟的边界条件和基本设置需满足以下规定: 1)计算区域:建筑覆盖区域小于整个计算域面积3%;以目标建筑为中心,半径5H围为水平计算域。建筑上方计算区域要大于3H;H为建筑主体高度; 2)网格划分:建筑的每一边人行高度区1.5m或2m高度应划分10个网格或以上; 3)湍流模型选择:标准k-ε模型。高精度要求时采用Durbin模型或MMK模型。

计算流体力学过渡到编程的傻瓜入门教程

借宝地写几个小短文,介绍CFD的一些实际的入门知识。主要是因为这里支持Latex,写起来比较方便。 CFD,计算流体力学,是一个挺难的学科,涉及流体力学、数值分析和计算机算法,还有计算机图形学的一些知识。尤其是有关偏微分方程数值分析的东西,不是那么容易入门。大多数图书,片中数学原理而不重实际动手,因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力,看了还不知道怎么实现。本人当年虽说是学航天工程的,但是那时本科教育已经退步,基础的流体力学课被砍得只剩下一维气体动力学了,因此自学CFD的时候也是头晕眼花。不知道怎么实现,也很难找到教学代码——那时候网络还不发达,只在教研室的故纸堆里搜罗到一些完全没有注释,编程风格也不好的冗长代码,硬着头皮分析。后来网上淘到一些代码研读,结合书籍论文才慢慢入门。可以说中间没有老师教,后来赌博士为了混学分上过CFD专门课程,不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛,不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统,而是希望能结合实际,阐明一些最基本的概念和手段,其中一些复杂的道理只是点到为止。目前也没有具体的计划,想到哪里写到哪里,因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西,看过之后,会知道一个CFD代码怎么炼成的(这“炼”字好像很流行啊)。欢迎大家提出意见,这样我尽可能的可以追加一些修改和解释。

言归正传,第一部分,我打算介绍一个最基本的算例,一维激波管问题。说白了就是一根两端封闭的管子,中间有个隔板,隔板左边和右边的气体状态(密度、速度、压力)不一样,突然把隔板抽去,管子内面的气体怎么运动。这是个一维问题,被称作黎曼间断问题,好像是黎曼最初研究双曲微分方程的时候提出的一个问题,用一维无粘可压缩Euler方程就可以描述了。 这里 这个方程就是描述的气体密度、动量和能量随时间的变化()与它们各自的流量(密度流量,动量流量,能量流量 )随空间变化()的关系。 在CFD中通常把这个方程写成矢量形式 这里 进一步可以写成散度形式

计算流体力学螺旋管分析报告

重庆大学《计算流体力学与计算传热学基础》上机实验水平螺旋管内的对流换热过程 学生:刘伟文 学号:20123000 指导教师:李隆键 专业:热能与动力工程 重庆大学动力工程学院 二O一五年六月

一、前言 螺旋管在热力、化工、石油及核工业等领域得到了广泛应用,螺旋管换热器也具有结构简单、传热系数高等优点。它的传热系数比直管高,在相同空间里可得到更大的传热面积,布置更长的管道,减少了焊缝,提高了安全性。尽管螺旋管的流体阻力增大,压降增大,但是其传热效率的提高导致能量的节约要高于因阻力增大而消耗的能量。因此,螺旋管在许多行业得到普遍应用而倍受青睐。在工程应用中,由于工艺要求,往往需将流体加热至规定的温度范围,传热是其中的基本单元操作,所以有必要对螺旋管的传热与流动特性进行研究。从理论知识我们知道由于向心力的作用,流体从管中心部分由螺旋管内侧流向外侧壁面,因而造成了螺旋管内侧的低压区。在压差作用下,流体从外侧沿着圆管的上部和下部壁面流回内侧。这种流动是与管的轴向垂直的,也就是与流体的主体流动相垂直,称为二次流。流体的这种二次流与轴向主流复合成螺旋式的前进运动。这样,对于流体的传热传质,不仅可依靠流体的径向扩散,还有径向二次流的作用,相当于边界层进行了破坏,增强了流体传质。 二、GAMBIT建模

1、先建立一个半径为6的圆面。 2、将该圆面向X轴正方向移动120。 3、用圆面sweep形成螺旋柱体。(绕Y轴正方向)

4、重复以上操作,得到如图所示几何体弯管。 5、设置边界层。

并应用至每个截面:

6、设置圆面的网格,选择pave方式,interval size 选择0.6,这样边界层网格与圆面中心网格过渡较平缓。 7、依次建立体网格。 8、检查网格质量。 最差网格为0.41,满足要求。 8、输出网格。

流体力学报告

流体力学报告 每一门力学学科的建立,都需要建立模型,也就是把实际的问题抽象化,而抽象过程就是把现实中对所研究问题不重要的因素忽略掉,也就是模型假设,从而建立于这个问题相适应的模型进行研究,如果有意义有价值,也就慢慢深入研究,从而形成一门学科,它们都是随社会的发展而发展形成的.比如现如今最前沿的力学学科"纳米力学"就是如此。我们土木工程常说的三大力学有:1.理论力学---分析力学,振动力学,水力学或称为流体力学(这些研究对材料都不太侧重 )2.材料力学---弹性力学,塑性力学(都是又材料特性而分的) 3.结构力学:就是分析复杂的结构的情形。在此我重点叙述我对流体力学这门课学科的学习和认知。 一·流体的基本信息解释: 流体,是与固体相对应的一种物体形态,是液体和气体的总称. 由大量的、不断地作热运动而且无固定平衡位置的分子构成的,它的基本特征是没有一定的形状并且具有流动性。流体都有一定的可压缩性,液体可压缩性很小,而气体的可压缩性较大,在流体的形状改变时,流体各层之间也存在一定的运动阻力(即粘滞性)。当流体的粘滞性和可压缩性很小时,可近似看作是理想流体,它是人们为研究流体的运动和状态而引入的一个理想模型。是液压传动和气压传动的介质。大气和水是最常见的两种流体,大气包围着整个地球,地球表面的70%是水面。大气运动、海水运动(包括波浪、潮汐、中尺度涡旋、环流等)乃至地球深处熔浆的流动都是流体的研究内容。

二·流体力学的阐述: 流体力学是连续介质力学的一门分支,是研究流体(包含气体,液体以及等离子态)现象以及相关力学行为的科学。可以按照研究对象的运动方式分为流体静力学和流体动力学,还可按流动物质的种类分为水力学,空气动力学等等。对流体力学学科的形成作出第一个贡献的是古希腊的阿基米德,他建立了包括物理浮力定律和浮体稳定性在内的液体平衡理论,奠定了流体静力学的基础,特别是从20世纪以来,流体力学已发展成为基础科学体系的一部分,同时又在工业、农业、交通运输、天文学、地学、生物学、医学等方面得到广泛应用。流体力学既包含自然科学的基础理论,又涉及工程技术科学方面的应用。此外,如从流体作用力的角度,则可分为流体静力学、流体运动学和流体动力学;从对不同"力学模型"的研究来分,则有理想流体动力学、粘性流体动力学、不可压缩流体动力学、可压缩流体动力学和非牛顿流体力学等。 三·对流体的研究假设: 连续体假设 物质都由分子构成,尽管分子都是离散分布的,做无规则的热运动.但理论和实验都表明,在很小的范围内,做热运动的流体分子微团的统计平均值是稳定的.因此可以近似的认为流体是由连续物质构成,其中的温度,密度,压力等物理量都是连续分布的标量场。 质量守恒 质量守恒目的是建立描述流体运动的方程组。欧拉法描述为:流进

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”的实体 操作中,左键选择,中键确认,右键完成并退出——类似的操作方法很多地方用到,要多练习,今后

高等流体力学试题

1.简述流体力学有哪些研究方法和优缺点? 实验方法就是运用模型实验理论设计试验装置和流程,直接观察流动现象,测量流体的流动参数并加以分析和处理,然后从中得到流动规律。实验研究方法的优点:能够直接解决工程实际中较为复杂的流动问题,能够根据观察到的流动现象,发现新问题和新的原理,所得的结果可以作为检验其他方法的正确性和准确性。实验研究方法的缺点主要是对于不同的流动需要进行不同的实验,实验结果的普遍性稍差。 理论方法就是根据流动的物理模型和物理定律建立描写流体运动规律的封闭方程组以及相应初始条件和边界条件,运 用数学方法准确或近似地求解流场,揭示流动规律。理论方法的优点是:所得到的流动方程的解是精确解,可以明确地给出各个流动参数之间的函数关系。解析方法的缺点是:数学上的困难比较大,只能对少数比较简单的流动给出解析解,所能得到的解析解的数目是非常有限的。 数值方法要将流场按照一定的规则离散成若干个计算点,即网格节点;然后,将流动方程转化为关于各个节点上流动 参数的代数方程;最后,求解出各个节点上的流动参数。数值方法的优点是:可以求解解析方法无能为力的复杂流动。数值方法的缺点是:对于复杂而又缺乏完整数学模型的流动仍然无能为力,其结果仍然需要与实验研究结果进行对比和验证。 2.写出静止流体中的应力张量,解释其中非0项的意义. 无粘流体或静止流场中,由于不存在切向应力,即p ij =0(i ≠j ),此时有 P =00000 0xx yy zz p p p ??????????=000000p p p -????-????-??=-p 00000011????1?????? = -p I 式中I 为单位张量,p 为流体静压力。 流体力学中,常将应力张量表示为 p =-+P I T (2-9) 式中p 为静压力或平均压力,由于其作用方向与应力定义的方向相反,所以取负值;T 称为偏应力张量,即 T =xx xy xz yx yy yz zx zy zz τττττττττ?????????? (2-10) 偏应力张量的分量与应力张量各分量的关系为:i =j 时,p ij 为法向应力,τii = p ij - p ;当i ≠j 时p ij 为粘性剪切应力,τij =p ij 。τii =0的流体称为非弹性流体或纯粘流体,τii ≠0的流体称为粘弹性流体。 3.分析可压缩(不可压缩)流体和可压缩(不可压缩)流动的关系. 当气体速度流动较小(马赫数小于0.3)时,其密度变化不大,或者说对气流速度的变化不十分敏感,气体的压缩性没有表现出来。因此,在处理工程实际问题时,可以把低速气流看成是不可压缩流动,把气体可以看作是不可压缩流体。而当气体以较大的速度流动时,其密度要发生明显的变化,则此时气体的流动必须看成是可压缩流动。 流场任一点处的流速v 与该点(当地)气体的声速c 的比值,叫做该点处气流的马赫数,用符号Ma 表示: Ma /v c v == (4-20) 当气流速度小于当地声速时,即Ma<1时,这种气流叫做亚声速气流;当气流速度大于当地声速时,即Ma>l 时,这种气流称为超声速气流;当气流速度等于当地声速时,即Ma=l 时,这种气流称为声速气流。以后将会看到,超声速气流和亚声速气流所遵循的规律有着本质的不同。 马赫数与气流的压缩性有着直接的联系。由式(4-11)可得 所以有 222Ma d ρv dv dv ρc v v =-=-。 (4-21) 当Ma≤0.3时,dρ/ρ≤0.09dv /v 。由此可见,当速度变化一倍时,气体的密度仅仅改变9%以下,一般可以不考虑密度的变化,即认为气流是不可压缩的。反之,当Ma>0.3时,气流必须看成是可压缩的。 4.试解释为什么有时候飞机飞过我们头顶之后才能听见飞机的声音. 5.试分析绝能等熵条件下截面积变化对气流参数(v ,p ,ρ,T )的影响.

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导入信息可以看到,完全没有问题。

计算流体力学结课报告

计算流体力学结课报告200Km/h列车fluent仿真计算 学部:化、环、生学部 学院:化工机械与安全学院 学号:31507095 班级:化1512班 学生姓名:孙金

引言 数值仿真就是对所建立的数值模型进行数值实验和求解的过程。而计算流体力学CFD (Computational Fluid Dynamics)就是在工程仿真实验领域中应用最广泛的一门学科。任何流体运动的规律都是以质量守恒定律、动量守恒定律和能量守恒定律为基础的。这些基本定律可由数学方程组来描述,如欧拉方程、N-S方程。采用数值计算方法,通过计算机求解这些控制流体流动的数学方程,进而研究流体的运动规律这就是CFD研究问题的方法。在实际计算流体力学方面,采用通用的CFD软件来完成工程上的一些流体力学问题,有极为广泛的应用前景。近年来,随着计算机技术以及相关技术的发展,CFD技术已经在工程领域内取得重大的进步,特别是在高速列车的外型设计方面起了很大作用。随着国家经济的发展,国家运输业也有了很大的发展,特别是列车经过几次提速后,高速列车在国家运输行业中所占比例不断提高。高速列车的特点是庞大、细长、在地面轨道上运行,其空气动力学问题非常复杂。空气在列车表面形成空气流场,空气阻力急剧增加,作用在列车的阻力大部分来自压强阻力,而一部分来自表面磨擦阻力,这就使能耗过大,同时列车可能出现较大的空气升力,导致列车产生“飘”的现象,激发列车脱轨事故的发生,因此研究高速列车气动力性能非常重要。用CFD仿真可以详细了解高速列车的空气动力特性,从而设计出阻力小、噪音低等各方面性能完善的高质量列车。本文采用CFD学科中的常用商业软件Fluent仿真一个时速200km/h的二维流线型车头的外流场,对其空气动力性能进行分析,从而得到不同车辆形状其周围流场的不同,进而分析哪种车型更适合。

高等流体力学

高等流体力学 第一章 流体力学的基本概念 连续介质:流体是由一个紧挨着一个的连续的质点所组成的,没有任何空隙的连续体,即所 谓的连续介质。 流体质点:是指微小体积内所有流体分子的总和。 欧拉法质点加速度:时变加速度与位变加速度和 z u u y u u x u u t u dt du a x z x y x x x x x ??+??+??+??== 质点的随体导数:质点携带的物理量随时间的变化率称为质点的随体导数,用dt d 表示。在欧拉法描述中的任意物理量Q 的质点随体导数表述如下: x k k Q u t Q dt dQ ??+??= 式中Q 可以是标量、矢量、张量。质点的随体导数公式对任意物理量都成立,故将质点的 随体导数的运算符号表示如下: x k k u t dt d ??+??= 其中 t ?? 称为局部随体导数,x k k u ??称为对流随体导数,即在欧拉法描述的流动中,物理 量的质点随体导数等于局部随体导数与对流随体导数之和。 体积分的随体导数:质点携带的物理量随时间的变化率称为质点的随体导数。则在由流体质点组成的流动体积V 中标量函数Φ(x, t )随时间的变化率就是体积分的随导函数。 由两部分组成①函数Φ 对时间的偏导数沿体积V 的积分,是由标量场的非恒定性引起的。②函数Φ通过表面S 的通量。由体积V 的改变引起的。 ()dV divv dt d dV v div t dS u dV t dV dt d v v n s v v ?? ? ???Φ+Φ=??????Φ+?Φ?=Φ+?Φ?=Φ??????????????()dV adivv dt da dV av div t a dS au dV t a adV dt d v v n s v v ?? ????+=??????+??=+??=?????????????? 变形率张量: 11ε 12ε13ε D ij = 21ε 22ε 23ε 31ε 32ε 33ε

相关主题
文本预览
相关文档 最新文档