fluent油水两相流动数值模拟

fluent油水两相流动数值模拟
fluent油水两相流动数值模拟

Fluent油水两相流弯管流动模拟

一、实例概述

选取某输油管道工程管径600mm的90°水平弯管道,弯径比为3,并在弯管前后各取5m直管段进行建模,其几何模型如图所示。为精确比较流体流经弯管过程中的流场变化,截取了图所示的5个截面进行辅助分析。弯管进出口的压差为800Pa,油流含水率为20%。

二、模型建立

1.启动GAMBIT,选择圆面生成面板的Plane为ZX,输入半径Radius为0.3,生成圆面,

如图所示。

2.选择圆面,保持Move被选中,在Global下的x栏输入1.8,完成该面的移动操作。

3.选取面,Angle栏输入-90,Axis选择为(0,0,0)→(0,0,1),生成弯管主体,如图所

示。

4.在Create Real Cylinder面板的Height栏输入5,在Radius1栏输入0.3,选择Axis

Location 为Positive X,生成沿x方向的5m直管段,如图所示。

5.同方法,改变Axis Location为Positive Y生成沿y方向的5m直管段,如图所示。

6.将直管段移动至正确位置,执行Volume面板中的Move/Copy命令,选中沿y轴的直管

段,在x栏输入1.8,即向x轴正向平移1.8。然后选中沿x轴的直管段,在x栏输入-5,在y栏输入-1.8,最后的模型如图所示。

7.将3个体合并成一个,弹出Unite Real Volumes面板,选中生成的3个体,视图窗口

如图所示。

三、网格划分

1.打开Create Boundary Layer面板,在Edges黄色输入栏中选取线3。选中1:1的边界

层生成方式,并设置第一个点距壁面距离为0.001m,递增比例因子为1.2,边界层为4层。绘制完边界层网格,如图所示。

2.打开Mesh Faces面板,运用Quad单元与Pave方法对该圆面进行划分,在Interval size

栏输入0.05,生成的面网格如图所示。

3.运用Hex/Wedge单元与Cooper方法对该圆面进行划分,在Interval size栏输入0.1,

生成的体网格如图所示。

4.选择面4作为进口,定义其为PRESSURE_INLET,命名为in;选择面9作为出口,定义

其为PRESSURE_OUTLET,命名为out,如图所示。

5.输出三维模型网格文件。

四、求解计算

1)启动FLUENT三维单精度计算器,读入划分好的网格文件,网格读入之后,检查网格的

信息,待最后一行出现“Done”语句。

2)从Scale Grid对话框中可以看出计算区域的尺寸范围,如图所示。

3)弹出的Solver对话框中保持默认设置,如图所示。

4)弹出的Viscous Model对话框,在Model选项中选择k-epsilon[2 eqn],在k-epsilon

Model选项下选择Realizable,保留其他默认设置。

5)如图所示,在弹出的Multiphase Model对话框中选择Mixture,保持默认参数设置。

6)在弹出的Fluent Database Materials对话框中,将数据库中的材料——水和油都拷贝

到当前工程中,完成材料物性的定义。

7)在Type选项中选择primary-phase,在列表中选择fuel-oil-liquid,在Name文本框

中输入oil。在Type选项中选择secondary-phase,在列表中选择water-liquid,并在Name文本框中输入water。

8)如图所示,在Set Units对话框中,选择Quantities中的pressure。在弹出的Define

Unit对话框中的Unit栏输入名称mpa,在Factor栏输入1000000,定义压强的新单位MPa。

9)打开Operating Conditions对话框,选中Gravity,指定重力方向为z轴,在Z[m/s2]

右侧输入-9.8。

10)定义边界条件。

11)在Solution Controls对话框中保持默认设置。

12)对流场进行初始化。

13)在Residual Monitors对话框中选中Plot,打开残差曲线图。保持FLUENT默认各参数

的收敛精度要求为0.001。

14)将文件保存。

15)设置迭代步为1000步,开始解算。

16)迭代完毕后,保存结果文件。

17)截取有助分析的几个截面,在Points下方的三点坐标中分别输入(-2.5,0,0)、

(-2.5,0,0.3)和(-2.5,0.3,0),在New Surface Name中输入x=-2.5,生成第一个截面,如图所示。

18)选取Contours of为Pressure…和Static Pressure,弹出压强云图。改变Contours of

为Velocity…和Velocity Magnitude,则显示为速度云图,如图所示。

压强云图

速度云图19)安全退出FLUENT。

midas数值模拟软件应用

某露天煤矿4-4剖面边坡稳定性分析与沿走向开采 的数值模拟 1概况 以实测4-4剖面为分析对象(如图1),根据钻孔资料确定上覆岩层属性,建立数值模拟分析模型,模型走向长300m、倾向234.17 m、高度为117.975m,模拟计算时需要考虑排土场附加荷载的影响。排土场高15.414m,其坡角35°,距离露天坡肩距离30m。具体各层参数如表1. 图1 实测4-4剖面分布图 表1岩体力学参数表 岩性 密度/ 103kg/m3 内摩擦角/° 凝聚力 /kPa 泊松比 弹性模量 /MPa 抗压强度/ MPa 表土 1.58 24 14 0.23 31.5 砂岩 2.537 33 111 0.25 5000 2.43 泥岩 2.314 34 52 0.35 1250 1.09 煤 1.45 32.7 201 0.30 1200

2二维数值模型 排土场高15m,其坡角35°,距离露天坡肩距离30m。二维模型共有1580个节点,1239个单元(如图2)。破坏判据采用莫尔-库仑准则。 2.1 二维网格划分 图2 4-4剖面二维数值模型 2.2 二维模型稳定性分析 2.2.1 稳定系数:1.3875 2.2.2 位移及应力云图如图2.2.2(a)、(b) 图2.2.2(a)4-4剖面Z方向位移变化色谱图

图2.2.2(b)4-4剖面Z方向应力变化色谱图 3三维模型 三维模型共有24692个节点,29736个单元(如图3)。破坏判据采用莫尔-库仑准则。模型参数取表1。沿走向开挖10步,前3步20m,中间4步10m,后3步20m,共开挖160m。 图3 4-4剖面三维数值模型 3.1第一步开挖 3.1.1位移云图

数值模拟软件大全

数值模拟软件大全 GEO-SLOPE Offical WebSite: www. geo-slope. com SLOPE/W: 专业的边坡稳定性分析软件, 全球岩土工程界首 选的稳定性分析软件 SEEP/W: 专业的地下渗流分析软件, 第一款全面处理非饱和土体渗流问题的商业化软件 SIGMA/W: 专业的岩土工程应力应变分析软件, 完全基于土(岩)体本构关系建立的专业有限元软件 QUAKE/W: 专业的地震应力应变分析软件, 线性、非线性土体的水平向与竖向耦合动态响应分析软件 TEMP/W: 专业的温度场改变分析软件, 首款最具权威、涵盖范围广泛的地热分析软件 CTRAN/W: 专业的污染物扩散过程分析软件, 超值实用、最具性价比的地下水环境土工软件 AIR/W:专业的空气流动分析软件, 首款处理地下水-空气-热相互作用的专业岩土软件 VADOSE/W: 专业的模拟环境变化、蒸发、地表水、渗流及地下水对某个区或对象的影响分析软件, 设计理论相当完善和全面的环境土工设计软件 Seep3D(三维渗流分析软件)是GeoStudio2007专门针对工程结构中的真实三维渗流问题, 而开发的一个专业软件, Seep3D软件将强大的交互式三维设计引入饱和、非饱和地下水的建模中, 使用户可以迅速分析各种各样的地下水渗流问题. 特点:GeoStudio其实就是从鼎鼎大名的GEO-SLOPE发展起来的, 以边坡分析出名, 扩展到整个岩土工程范围, 基于. NET平台开发的新一代岩土工程仿真分析软件, 尤其是VADOSE/W模块是极具前瞻性的, 环境岩土工程分析的利器. 遗憾的是其模块几乎都只提供平面分析功能. Rocscience Offical WebSite: www. rocscience. com Rocscience 软件的二维和三维分析主要应用在岩土工程和 采矿领域, 该软件使岩土工程师可以对岩质和土质的地表 和地下结构进行快速、准确地分析, 提高了工程的安全性并 减少设计成本. Rocscience 软件对于岩土工程分 析和设计都很方便, 可以帮助工程师们得到快速、正确的解答. Rocscience 软件对于用户最新的项目都有高效的解算结果, 软件操作界面是基于WINDOWS 系统的交互式界面. Rocscience 软件自带了基于CAD 的绘图操作界面, 可以随意输入多种格式的数据进行建模, 用户可以快速定义模型的材料属性、边界条件等, 进行计算得到自己期望的结果. Rocscience 软件包括以下十三种专业分析模块: Slide 二维边坡稳定分析模块

T型微通道内两相流动数值模拟和流场

价值工程 —————————————————————— —基金项目:西安市科技计划项目(CXY1134WL09) 。作者简介:王琳琳(1981-),女,陕西西安人,西安文理学院数学 与计算机工程学院,讲师,西安交通大学能动学院博 士生,研究方向为微通道内的两相流动。 0引言 微通道的尺寸非常小,其通道的宽度一般在之间,流量小[1],借助微通道可以进行两相流体的混合、纳米粒子合成、 蛋白质结晶等。在化工方面,要求能够控制微通道内化学物质输运的时间和物质空间的分布[2,3] 。近年来,研究者对不同结构微通道内流动的控制产生了极大的兴趣,成为 一个重要的研究方向[4]。 雷诺数是惯性力和黏性力之比,微通道内雷诺数小,两相流动受到黏性力的影响,在通道壁面约束下,表面张力和挤压力对离散相的形成起到重要作用。微通道的制作工艺精度较高,监测通道内流动的设备需要极其微小,这些都使得采用实验研究微流动的难度和费用较大,而数值模拟能够克服这些缺点。两相流动问题常见的数值模拟方法有:标记网格方法(MAC ),水平集方法(level set method ),相场方法(phase field method ),VOF 法,格子布尔兹曼方法(Lattice Boltzmann method )等,在这些的数值方法中,相场方法利用自由能量描述两相流体的界面,模拟中采用非结构化网格时容易实施,对流场的计算中不用重新初始化,物质的质量损失较小,控制方程中的变量具有一定物理意义,并能够模拟能量耗散的流动[5]。本文采用相场方法,数值模拟工程中常见的错流接触T 型微通道内离散相的形成过程,研究微通道内压强和流场的变化特点。 1控制方程 连续性方程和动量方程为: 塄· v 軆=0坠v 軆坠t +(v 軆·塄)v 軆=塄·[-p ρI+μρ(塄v 軆+(塄v 軆)T ]+1ρ F 軋σ+g 軆軋軋軋軋軋軋軋軋軋軋軋 , 其中v 軆是速度向量,p 是压强,ρ是密度,μ是动力粘性系数,F 軋σ是表面张力,σ是表面张力系数。由相场理论知,两流体间的相互作用可用自由能量密度f mix (准,塄准)=1λ 塄准2+f 0(准 )来表示,式中的第一项1λ塄准2 是两相流体界面内的能量密度,第二项f 0(准 )=λ4∈ 2(准2-1)2 是各个流体的块能量密度,λ是混合能量密度的参数,∈表示两流体界面的厚度,准是相场变量,微通道内离散相和连续 相对应的准值分别是-1和1, 准在-1和1之间变化对应的区域就是两相界面,自由能量密度反映了两相流体间的相互作用。对自由能量密度在计算区域内积分,得到自由能F ,即F=乙v f mix dv ,F 关于相场变量的变化率是化学势G , 即G=坠F 坠准 ,由自由能的定义可得到,G=f ′ 0(准)-λ塄2准。Van der Waals 假定流场中自由能最小处就是平衡的两相界面,因此两相界面满足方程坠F 坠准 =0。通过计算可得到平衡 的一维两相界面的表达式是准(x )=tanh (x 2姨∈ )。平衡 两相界面单位长度的自由能理解成表面张力系数[6],即σ= λ +∞ -∞乙 1(d 准)2 +f 0(准)dx ,结合前面的定义和公式,得到表面张力系数、两相界面厚度和混合能量密度的关系式 σ=22姨3λ∈ 。 表面张力可用化学势表示:F σ 姨=G 塄准。上面方程组结合Cahn-Hilliard 对流方程坠准坠t +v 軆·塄准=塄·(γ塄G )就是计 算微通道内两相流动的控制方程。 2T 型微通道内离散相的形成研究的物理模型是T 型微通道,通道宽度D=111μm ,离散相通道和主通道垂直,两通道长度分别是3D 和45D ,离散相和连续相分别从垂直方向和水平方向同时注入通道,在一定条件下形成间距固定的离散相。微通道内可忽略重力作用,为简化计算,设置两相流体密度相同,表面张 T 型微通道内两相流动数值模拟和流场分析 Numerical Simulation of Two-phase Flow and Flow Field Analysis in a T-junction Micro-channel 王琳琳①②WANG Lin-lin ;胡洪萍①HU Hong-ping (①西安文理学院,西安710065;②西安交通大学,西安710049) (①Xi'an University o f Arts and Science ,Xi'an 710065,China ;②Xi'an Jiaotong University ,Xi'an 710049,China ) 摘要:借助相场方法数值模拟T 型微通道内两相流动,通过改变毛细数大小,得到三种形成机理下的离散相。随着毛细数增大, 离散相形成过程对微通道内压强和速度的影响减弱。 Abstract:The two-phase flow was simulated in a T-junction micro-channel by using the phase field method,and three type droplets were obtained with different capillary number.We found that the influence of droplet formation on pressure and velocity became weak as the capillary number increases. 关键词:相场方法;数值模拟;微通道;毛细数Key words:phase field method ;numerical simulation ;micro-channel ;capillary number 中图分类号:TQ021.1文献标识码:A 文章编号:1006-4311(2012)31-0180-02 ·180·

fluent模拟基本步骤及注意事项

二维模拟: 一、模拟类型: 1、 大区域空间速度场模拟 计算区域大小设置:迎风面是建筑长度的3倍,背风面是建筑长度的12倍,两侧面是建筑宽度的3倍,高度是建筑高度的4倍。 根据相似理论:l C -几何比例尺 速度比例尺:2 10l C C =υ 风量比例尺:2520l l Q C C C C =?=υ 热量比例尺: 250l T Q C C C Cq =?=? 2、 建筑户型温度场、速度场模拟 二、基本操作步骤及注意事项: A gambit 建模 1、 建模: 方法一:直接在GAMBIT 建模; 方法二:CAD 导入gambit ; 1) 在CAD 中用PL 线将户型的基本构造画出来,创建为面域; 2) 输入命令acisoutver ,把‘70’修改为‘30’。 3) “文件”——“输出”——sat 文件 4) 在gambit 中导入Acis 文件 注意:在用PL 线构画户型时,在进口和出口边界(窗户、内户门),要各边界端点连续画线。 2、 划分网格: Interval Size :50 3、 设置边界条件 内部开口边界(门)设置为internal ,房间相邻墙壁设置为Wall 4、 保存文件,并输出mesh 文件 B 导入fluent 计算: 1、 导入mesh 文件 2、 检查网格 3、 设置单位 gambit 里可以缩小建筑比例建模,在fluent 中设置单位恢复原模型。 4、 选择计算模型 5、 设置材料类型 6、 设置边界条件 7、 设置模拟控制条件 8、 边界初始化

9、设置监视窗口 10、设置迭代次数进行计算 11、结果显示 12、保存文件 三、需解决问题: 1、湍流强度等计算; 2、层流湍流界定问题; 3、壁面湿度设置问题; 四、待提高部分: 1、户型流场模拟时,墙壁考虑采用双钱; 2、南京理工校区原始模型(不简化)模拟; 3、三维模型模拟; 五、

地下水数值模拟任务、步骤及常用软件

地下水数值模拟任务、步骤及常用软件1地下水模拟任务 大多数地下水模拟主要用于预测,其模拟任务主要有4种: 1)水流模拟 主要模拟地下水的流向及地下水水头与时间的关系。 2)地下水运移模拟 主要模拟地下水、热和溶质组分的运移速率。这种模拟要特别考虑到“优先流”。所谓“优先流”就是局部具有高和连通性的渗透性,使得水、热、溶质组分在该处的运移速率快于周围地区,即水、热、溶质组分优先在该处流动。 3)反应模拟 模拟水中、气-水界面、水-岩界面所发生的物理、化学、生物反应。 4)反应运移模拟 模拟地下水运移过程中所发生的各种反应,如溶解与沉淀、吸附与解吸、氧化与还原、配合、中和、生物降解等。这种模拟将地球化学模拟(包括动力学模拟)和溶质运移模拟(包括非饱和介质二维、三维流)有机结合,是地下水模拟的发展趋势。要成功地进行这种模拟,还需要研究许多水-岩相互作用的化学机制和动力学模型。 2模拟步骤 对于某一模拟目标而言,模拟一般分为以下步骤: 1)建立概念模型 根据详细的地形地貌、地质、水文地质、构造地质、水文地球化学、岩石矿物、水文、气象、工农业利用情况等,确定所模拟的区域大小,含水层层数,维数(一维、二维、三维),水流状态(稳定流和非稳定流、饱和流和非饱和流),介质状况(均质和非均质、各向同性和各向异性、孔隙、裂隙和双重介质、

流体的密度差),边界条件和初始条件等。必要时需进行一系列的室内试验与野外试验,以获取有关参数,如渗透系数、弥散系数、分配系数、反应速率常数等。 2)选择数学模型 根据概念模型进行选择。如一维、二维、三维数学模型,水流模型,溶质运移模型,反应模型,水动力-水质耦合模型,水动力-反应耦合模型,水动力-弥散-反应耦合模型。 3)将数学模型进行数值化 绝大部分数学模型是无法用解析法求解的。数值化就是将数学模型转化为可解的数值模型。常用数值化有有限单元法和有限差分法。 4)模型校正 将模拟结果与实测结果比较,进行参数调整,使模拟结果在给定的误差范围内与实测结果吻合。调参过程是一个复杂而辛苦的工作,所调整的参数必须符合模拟区的具体情况。所幸的是,最近国外已花费巨力开发研究了自动调参程序(如PEST),大大提高了模拟者的工作效率。 5)校正灵敏度分析 校正后的模型受参数值的时空分布、边界条件、水流状态等不确定度的影响。 灵敏度分析就是为了确定不确定度对校正模型的影响程度。 6)模型验证 模型验证是在模型校正的基础上,进一步调整参数,使模拟结果与第二次实测结果吻合,以进一步提高模型的置信度。 7)预测 用校正的参数值进行预测,预测时需估算未来的水流状态。

实体入水FLUENT模拟过程_包括划分网格方法

实体入水模拟过程 3.2.1利用GAMBIT建立计算模型 1)启动GAMBIT,打开对话框如图3.2.1选择工作目录为D:\GAMBIT working。 图 3.2.1 2)首先建立等边三角形,单击Geometry Vertex Create Real Vertex,在Create Real Vertex面板的x、y、z坐标输入(0,0,0),单击Apply按钮生成第一个点,按同样的方法建立点(0.4,0,0)。然后单击Geometry Edge Create Straight Edge,在Create Straight Edge面板中选择点1与点2,连接这两点省成线段。如图3.2.2 图3.2.2 3)单击Edge面板中的Move/Copy Edges按钮,打开如图3.2.3的面板,选择线段1,单击copy按钮,并选择Operation为Rotate,在Angle栏输入60,其他保持默认,单击Apply 按钮。即旋转复制生成第二条线段。

图3.2.3 4)剩下的一条线段只需连接右侧两点即可,如图3.2.4所示。 图3.2.4 5)创建三角形面。单击Geometry Face Create Face from Wireframe,在Create Face from Wireframe面板中利用鼠标左键框选等边三角形的三条边,然后单击Apply按钮创建面。 6)由于三角形面域的位置不对,所以还要对其位置进行调整。首先需将其旋转210度。单击Face面板中的Move/Copy Faces按钮,在Move/Copy Faces面板中,选择面1(face.1),单击Move并选择Operation为Rotate,在Angle栏输入210,其他保持默认,单击Apply 按钮。其次,需要将三角形平移,在Move/Copy Edges面板中选择面1(face.1),单击Move 并选择Operation为Translate,在x与y栏分别输入3和8.4,单击Apply按钮完成平移操作,此时的视图窗口如图3.2.5所示。

fluent油水两相流动数值模拟

Fluent油水两相流弯管流动模拟 一、实例概述 选取某输油管道工程管径600mm的90°水平弯管道,弯径比为3,并在弯管前后各取5m直管段进行建模,其几何模型如图所示。为精确比较流体流经弯管过程中的流场变化,截取了图所示的5个截面进行辅助分析。弯管进出口的压差为800Pa,油流含水率为20%。 二、模型建立 1.启动GAMBIT,选择圆面生成面板的Plane为ZX,输入半径Radius为0.3,生成圆面, 如图所示。

2.选择圆面,保持Move被选中,在Global下的x栏输入1.8,完成该面的移动操作。 3.选取面,Angle栏输入-90,Axis选择为(0,0,0)→(0,0,1),生成弯管主体,如图所 示。

4.在Create Real Cylinder面板的Height栏输入5,在Radius1栏输入0.3,选择Axis Location 为Positive X,生成沿x方向的5m直管段,如图所示。 5.同方法,改变Axis Location为Positive Y生成沿y方向的5m直管段,如图所示。

6.将直管段移动至正确位置,执行Volume面板中的Move/Copy命令,选中沿y轴的直管 段,在x栏输入1.8,即向x轴正向平移1.8。然后选中沿x轴的直管段,在x栏输入-5,在y栏输入-1.8,最后的模型如图所示。 7.将3个体合并成一个,弹出Unite Real Volumes面板,选中生成的3个体,视图窗口 如图所示。

三、网格划分 1.打开Create Boundary Layer面板,在Edges黄色输入栏中选取线3。选中1:1的边界 层生成方式,并设置第一个点距壁面距离为0.001m,递增比例因子为1.2,边界层为4层。绘制完边界层网格,如图所示。 2.打开Mesh Faces面板,运用Quad单元与Pave方法对该圆面进行划分,在Interval size 栏输入0.05,生成的面网格如图所示。

自吸泵气液两相流数值模拟分析

2009年9月 农业机械学报 第40卷第9期 自吸泵气液两相流数值模拟分析3 刘建瑞 苏起钦 (江苏大学流体机械工程技术研究中心,镇江212013) 【摘要】 采用Mixture 多相流模型、Realizable 湍流模型与SIMPL EC 算法,应用CFD 软件Fluent 对内混式自吸泵自吸过程的气液两相流进行了数值模拟。通过分析不同含气率条件下流场的压力分布、速度分布、气相分布,探讨了气液两相介质在泵内的运动情况,一定程度上揭示了内混式自吸泵自吸过程的内部流场变化规律,为自吸泵的设计提供更多的参考依据。 关键词:自吸泵 气液两相流 数值模拟 自吸性能中图分类号:TH317 文献标识码:A Numerical Simulation on G as 2liquid Two 2phase Flow in Self 2priming Pump Liu Jianrui Su Qiqin (Technology and Research Center of Fluid M achinery Engineering ,Jiangsu U niversity ,Zhenjiang 212013,China ) Abstract 32D simulation was performed for the gas 2liquid two 2phase turbulent flow in self 2priming pump by using Fluent software with Mixture model ,SIMPL EC algorithm and Realizable turbulence model.The gas 2liquid two 2phase flow in self 2priming pump was investigated in the pressure in the pump ,the velocity in the pump ,the distribution of gas and liquid phase in the pump.To some extent ,the results reveal the self 2suction process of the two 2phase flow in self 2priming pump ,and provide references for self 2priming pump design. K ey w ords Self 2priming pump ,G as 2liquid two 2phase flow ,Numerical simulation ,Self 2priming capability 收稿日期:2008210229 修回日期:20092022193国家“863”高技术研究发展计划资助项目(2006AA100211)和江苏省科技支撑计划项目(BE2008381)作者简介:刘建瑞,教授,博士生导师,主要从事流体机械及工程研究,E 2mail :ljrwjj @https://www.360docs.net/doc/5914006134.html, 引言 自吸离心泵自20世纪60年代开始研究以来发 展迅速,产品已成系列。然而目前自吸泵的理论还不完善。以往对自吸性能的研究大多采用实验方法,但由于自吸泵自吸过程是气液两相流输送过程,其内部流场较为复杂,而且影响自吸泵自吸性能的因素较多。采用实验方法势必要耗费大量的人力、物力,而且也难以为自吸泵的设计提供准确而充分的理论依据,因此需要对自吸泵自吸过程的内部流场进行深入研究[1]。近年来,随着计算流体力学和相应计算软件的发展,对自吸泵的三维数值模拟已成为可能[2]。本文借助Fluent6.2软件平台提供的Mixture 多相流模型对内混式自吸泵自吸过程的内 部流场进行三维数值模拟。分析泵进口不同含气率条件下的模拟结果,初步揭示内混式自吸泵自吸过程中气液两相流速度场分布、压力场分布、含气率分布规律,以期为自吸泵的优化设计提供更充分的理论依据。 1 数学模型的建立 数值模拟所选内混式自吸泵的结构,如图1所示。111 基本假设 (1)假定整个流场相对运动定常,绝对运动有势,且液相为不可压缩流体、气相为不可压缩理想气体。 (2)流场中气泡的直径较小,可以忽略气泡对流

Fluent数值模拟步骤

Fluent数值模拟的主要步骤 使用Gambit划分网格的工作: 首先建立几何模型,再进行网格划分,最后定义边界条件。 Gambit中采用的单位是mm,Fluent默认的长度是m。 Fluent数值模拟的主要步骤: (1)根据具体问题选择2D或3D求解器进行数值模拟; (2)导入网格(File-Read-Case),然后选择由Gambit导出的msh文件。 (3)检查网格(Grid-Check),如果网格最小体积为负值,就要重新进行网格划分。(4)选择计算模型(Define-Models-Solver)。(6) (5)确定流体的物理性质(Define-Materials)。 (6)定义操作环境(Define-Operating Conditions)。 (7)指定边界条件(Define-Boundary Conditions )。 (8)求解方法的设置及其控制(Solve-Control-Solution)。 (9)流场初始化(Solve-Initialize)。 (10)打开残插图(Solve-Monitors-Residual)可动态显示残差,然后保存当前的Case和Data文件(File-Writer-Case&Data)。 (11)迭代求解(Solve-Iterate)。 (12)检查结果。 (13)保存结果(File-Writer-Case&Data),后处理等。 在运行Fluent软件包时,会经常遇到以下形式的文件: .jou文件:日志文档,可以编辑运行。 .dbs文件:Gambit工作文件,若想修改网格,可以打开这个文件进行再编辑。 .msh文件:Gambit输出的网格文件。 .cas文件:是.msh文件经过Fluent处理后得到的文件。 .dat文件:Fluent计算数据结果的数据文件。 三维定常速度场的计算实例操作步骤 对于三维管道的速度场的数值模拟,首先利用Gambit画出计算区域,并且对边界条件进行相应的指定,然后导出Mesh文件。接着,将Mesh文件导入到Fluent求解器中,再经过一些设置就得到形影的Case文件,再利用Fluent求解器进行求解。最后,可以将Fluent 求解的结果导入到Tecplot中,并对感兴趣的结果进行进一步的处理。

Eclipse数值模拟软件问答

Eclipse数值模拟软件问答(初级) 1. ECLIPSE输出结果文件是哪些? .GRID或.FGRID: 网格文件 .EGRID: 网格文件,与GRID格式不同,文件要小的多。(用关键字GRIDFILE来控制输出类型).INIT或.FINIT: 属性文件。(用关键字INIT来控制输出) .PRT: 报告输出。文件很大,模型处理及计算结果详细报告。(RPTGRID,RPTPROP,RPTSOL,RPTSCHED控制输出) .LOG: 后台作业时的输出报告,文件比PRT要小很多。可用于错误检查。 .DBG: Debug文件,一般不用。可用于检查ECLIPSE如何处理输入参数。 .SAVE: 用于快速重启。(用关键字SAVE来控制输出) .RFT:RFT计算结果。(用关键字WRFTPLT来控制输出) .FLUX: 流动边界。(用关键字DUMPFLUX来控制输出) .Snnnn或.UNSMRY: 图形文件输出(在SUMMARY部分定义) .Xnnnn或.UNRST: 重启文件输出(用RPTRST,RPTSOL或RPTSCHED来控制输出) 2. ECLIPSE输出文件都有什么格式? 格式化输出:可读文件,文件大。(用关键字FMTOUT来控制) 非格式化输出:不可读文件,文件小。 多输出文件:每一时间步一个输出文件。 单文件输出:所有时间步输出到一个文件。(用关键字UNIFOUT来控制) ECLIPSE缺省输出:非格式化,多文件输出。 3. ECLIPSE数据文件分几部分,各部分定义什么数据类型? ECLIPSE数据类型分八部分,各部分内的关键字除几个个别的外不能混用。 RUNSPEC: 定义模型维数以及模型基本类型,包括模型网格维数,最大井数,井组数,流体类型,输出类型控制等。 GRID: 定义模型网格和属性,包括顶部深度,厚度,孔隙度,渗透率,净毛比,一般由前处理软件Flogrid或Petrel输出。 EDIT: 编辑孔隙体积,传导率。 PROPS: 流体PVT及岩石数据,包括油、气体积系数,粘度随压力变化,水的体积系数,粘度;油,气,水地面密度等。岩石数据是相渗曲线和毛管压力。 REGIONS: 分区数据,包括流体分区,岩石分区,储量区,平衡区等。 SOLUTION: 平衡区数据,包括油水界面、油气界面、参考压力、参考深度、水体参数。 SUMMARY: 计算结果输出,包括油田,井组,单井的油、气、水产量,压力输出,网格的压力,

槽道颗粒两相流的数值模拟简介

槽道颗粒两相流的数值模拟简介 使用直接数值模拟和拉格朗日粒子追踪法研究小固体颗粒对水平槽道近壁区域流体的影响,对直接数值模拟的颗粒槽道流动的运算结果,采用条件抽样的方法分析流体相干结构对颗粒的影响。把流体对颗粒的作用简化为点力模型,同时假定碰撞为完全弹性碰撞。我们主要研究自由流体的平均速度,均方根速度,涡旋角速度均方根;以及对于加入颗粒的槽道流,我们打算先加入单行的颗粒程序,研究流体对颗粒的平均速度,均方根速度,数目随高度,浓度随高度的变化。再对两者进行耦合,研究两者之间的相互作用。 关键词:直接数值模拟(DNS),拉格朗日粒子追踪法,点力模型,相干结构 含有颗粒的流体常常人们的生活中以及一些工业的加工生产中。流态化技术在化工生产中也是一项极其重要的技术,但是颗粒两相流往往呈现了多变性和复杂性;以气体固体液态化为例,当颗粒的通量,流体的速度不同时,其会呈现不同的液态形式;早在19世纪就有了关于明渠中关于泥沙的转移和沉积的两相流的研究,在60年代后,出现了对该问题的一些有关两相流的基本方程的研究及相关的著作;在过去数年间,人们从不同的动机出发,慢慢的建立了对于两相流的三大模型:连续介质模型,离散颗粒模型,流体拟颗粒模型(从刻划单颗粒尺度的运动规律入手 ,不仅仅是把宏观的离散的颗粒当成离散相处理 ,同时还把宏观的连续的流体以拟“颗粒”性质的流体微团来处理 ,来模拟非平衡态的系统,这类模型都是对流体和颗粒的运动采用拉格朗日坐标)物理概念不一样的颗粒流体两相流模型方程组的形式各异 ,即使是同一类模型 ,其方程组形式及各项参数也往往存在着差异。 在过去的几年间,有关颗粒方面的研究也有很多;尽管两相流数值模型取得了较大的进展,但是仍有许多问题亟待解决。跟我们研究比较接近的是Chris D.Dritselis 和Nicholas S.Vlachos在2008年发表的颗粒槽道流在近壁区域相干结构的研究;研究的主要内容是相干结构由于颗粒而引发的修正,其研究的方法是:直接数值模拟(DNS是通过直接求解控制流动的 N-S 方程,利用所得到的数值解来研究湍流的方法。近三十年的研究认为 N-S 方程完全可以描述包括湍流在内的复杂流动。直接求解 N-S 方程可以得到整个湍流场的全部信息,无

数值模拟报告

第一部分:数值模拟技术研究文献综述 浅析数值模拟技术 1.引言 近年来,随着我国大规模地进行“西部大开发”和“南水北调”等巨型工程,越来越多的岩土工程难题摆在我们面前,单纯依靠经验、解析法显然已不能有效指导工程问题的解决,迫切需要更强有力的分析手段来进行这些问题的研究和分析。自. Clough 上世纪60年代末首次将有限元引入某土石坝的稳定性分析以来,数值模拟技术在岩土工程领域取得了巨大的进步,并成功解决了许多重大工程问题。特别是个人电脑的普及及计算性能的不断提高,使得分析人员在室内进行岩土工程数值模拟成为可能。在这样的背景下,数值模拟特别是三维数值模拟技术逐渐成为当前中国岩土工程研究和设计的主流方法之一,也使得岩土工程数值模拟技术成为当今高校和科研院所岩土工程专业学生学习的一个热点。 采用大型通用软件对岩土工程进行数值模拟计算,在目前已成为项目科研、工程设计、风险评估等岩土类项目的必须,学习和掌握Ansys、FLAC3D、UDEC等数值计算软件已成为学校、科研院所对工程从业人员的基本要求。 数值模拟方法主要有限元法、边界元法、加权余量法、半解析元法、刚体元法、非连续变形分析法、离散元法、无界元法和流形元法等,各种方法都有其对应的软件。 2.数值模拟的发展趋势 可以说, 继理论分析和科学试验之后, 数值模拟已成为科学技术发展的主要手段之一。随着软件技术和计算机技术的发展, 目前国际上数值模拟软件发展呈现出以下一些趋势: (1). 由二维扩展为三维。早期计算机的能力十分有限,受计算费用和计算机储存能力的限制,数值模拟程序大多是一维或二维的,只能计算垂直碰撞或球形爆炸等特定问题。随着第三代、第四代计算机的出现, 才开始研制和发展更多的三维计算程序。现在,计算程序一般都由二维扩展到了三维,如LS-DYNA2D 和LS - DYNA3D、AUTODYN2D 和 AUTO-DYN3D。 (2).从单纯的结构力学计算发展到求解许多物理场问题。数值模拟分析方法最早是从结构化矩阵分析发展而来,逐步推广到板、壳和实体等连续体固体力学分析,实践证明这是一种非常有效的数值模拟方法。近年来数值模拟方法已发展到流体力学、温度场、电传导、磁场、渗流等求解计算,最近又发展到求解几个交叉学科的问题。例如内爆炸时,空气冲击波使墙、板、柱产生变形,而墙、板、柱的变形又反过来影响到空气冲击波的传播,这就需要用固体力学和流体动力学的数值模拟结果交叉迭代求解。 (3).由求解线性问题进展到分析非线性问题。随着科学技术的发展,线性理论已经远远不能满足设计的要求。诸如岩石、土壤、混凝土等,仅靠线性计算理论就不足以解决遇到的问题,只有采用非线性数值算法才能解决。众所周知,非线性的数值模拟是很复杂的,它涉及到很多专门的数学问题和运算技巧,很难为一般工程技术人员所掌握。为此,近年来国外一些公司花费了大量的人力和资金,开发了诸如LS- DYNA3D、ABAQUS和AU-TODYN等专长求解非线性问题的有限元分析软件,并广泛应用于工程实践。这些软件的共同特点是具有高效的非线性

FLUENT-6-计算模拟过程方法及步骤

FLUENT 12 模拟步骤 Problem Setup 读入网格:file read case 选择网格文件(后缀为。Mesh) 1 General 1)Mesh(网格) > Check(点击查看网格的大致情况,如有无负体积等) Maximum volume (m3)(最大体积,不能为负) Minimum volume (m3)(最小体积,不能为负) Total volume (m3)(总体体积,不能为负) > Report Quality(点击报告网格质量) Maximum cell squish(最大单元压扁,如果该值等于1,表示得到了很坏的单元) Maximum cell skewness(最大单元扭曲,该值在0到1之间,0表示最好,1表示最坏) Maximum aspect ratio(最大长宽比,1表示最好) > Scale(点击缩放网格尺寸,FLUENT默认的单位是米) Mesh Was Create In(点选mm →点击Scale按钮且只能点击一次) View Length Unit In(点选mm →直接点击Close按钮不能再点击Scale按钮) > Display(点击显示网格设定)→弹出Mesh Colors窗口 Options(选Edges和Faces) Edge Type(点选All) Surface(点选曲面) →点击Display按钮 点击Colors按钮→弹出Mesh Display窗口 Options(点选Color by ID) →点击Close按钮→再点击Display按钮 2)Solver(求解器) > Pressure-Based(压力基,压力可变,用于低速不可压缩流动) > Density-Based(密度基,密度可变,用于高速可压缩流动) 3)Velocity Formulation(速度格式) > Absolute(绝对速度) > Relative(相对速度) 4)Time(时间) > Steady(稳态) > Transient(瞬态) 5)Units(点击设置变量单位) 点击按钮→弹出Set Units窗口→在Quantities项里点选pressure →在Units项里点选atm →点击

多孔介质-Fluent模拟

7.19多孔介质边界条件 多孔介质模型适用的范围非常广泛,包括填充床,过滤纸,多孔板,流量分配器,还有管群,管束系统。当使用这个模型的时候,多孔介质将运用于网格区域,流场中的压降将由输入的条件有关,见Section 7.19.2.同样也可以计算热传导,基于介质和流场热量守恒的假设,见Section 7.19.3. 通过一个薄膜后的已知速度/压力降低特性可以简化为一维多孔介质模型,简称为“多孔跳跃”。多孔跳跃模型被运用于一个面区域而不是网格区域,而且也可以代替完全多孔介质模型在任何可能的时候,因为它更加稳定而且能够很好地收敛。见Section 7.22. 7.19.1 多孔介质模型的限制和假设 多孔介质模型就是在定义为多孔介质的区域结合了一个根据经验假设为主的流动阻力。本质上,多孔介质模型仅仅是在动量方程上叠加了一个动量源项。这种情况下,以下模型方面的假设和限制就可以很容易得到: ?因为没有表示多孔介质区域的实际存在的体,所以fluent默认是计算基于连续性方程的虚假速度。做为一个做精确的选项,你可以适用fluent 中的真是速度,见section7.19.7。 ?多孔介质对湍流流场的影响,是近似的。 ?当在移动坐标系中使用多孔介质模型的时候,fluent既有相对坐标系也可以使用绝对坐标系,当激活相对速度阻力方程。这将得到更精确的源项。 相关信息见section7.19.5和7.19.6。 ?当需要定义比热容的时候,必须是常数。 7.19.2 多孔介质模型动量方程 多孔介质模型的动量方程是在标准动量方程的后面加上动量方程源项。源项包含两个部分:粘性损失项(达西公式项,方程7.19-1右边第一项),和惯性损失项(方程7.19-1右边第二项) (7.19-1)

ECLIPSE 油藏数值模拟软件 使用经验

ECLIPSE 油藏数值模拟软件使用经验 第一:ECLIPSE 介绍 ECLIPSE100,ECLIPSE300和FrontSim是主模拟器。ECLISPE100是对黑油模型进行计算,ECLISPE300是对组分模型和热采模拟进行计算,FrontSim是流线法模拟器。前处理模块有Flogrid,PVTi,SCAL,Schedule,VFPi等。Flogrid用于为数值模拟建立模拟模型,包括油田构造模型和属性模型;PVTi用于为模拟准备流体的PVT参数,对于黑油模型,主要是流体的属性随地层压力的变化关系表,对于组分模型是状态方程;SCAL为模型准备岩石的相渗曲线和毛管压力输入参数;Schedule处理油田的生产数据,输出ECLIPSE需要的数据格式(关键字);VFPi是生成井的垂直管流曲线表,用于模拟井筒管流。ECLIPSE OFFICE和FLOVIZ 是后处理模块,进行计算曲线和三维场数据显示和分析,ECLIPSE OFFICE同时也是ECLIPSE的集成平台。 对于初学者,不但要学主模型,也需要学前后处理。对于ECLISPE的初学者,应该先从ECLISPE OFFICE学起,把ECLISPE OFFICE的安装练习做完。然后再去学Flogrid,Schedule 和SCAL。PVTi主要用于组分模型,做黑油模型可以不用。 第二:做油藏数值模拟都需要准备什么参数 在照着软件提供的安装例子做练习时经常遇到的问题是:虽然一步一步按照手册的说明做,但做的时候不明白每一步在做什么,为什么要这么做。这时候的重点在于你要知道你一开始做的工作都是为数值模拟计算提供满足软件格式要求的基础参数。有了这些基础参数你才能开始进行模拟计算。这些基础参数包括以下几个部分: 1。模拟工作的基本信息:设定是进行黑油模拟,还是热采或组分模拟;模拟采用的单位制(米制或英制);模拟模型大小(你的模型在X,Y,Z三方向的网格数);模拟模型网格类型(角点网格,矩形网格,径向网格或非结构性网格);模拟油藏的流体信息(是油,气,水三相还是油水或气水两相,还可以是油或气或水单相,有没有溶解气和挥发油等);模拟油田投入开发的时间;模拟有没有应用到一些特殊功能(局部网格加密,三次采油,端点标定,多段井等);模拟计算的解法(全隐式,隐压显饱或自适应)。 2。油藏模型:模型在X,Y,Z三方向的网格尺寸大小,每个网格的顶面深度,厚度,孔隙度,渗透率,净厚度(或净毛比)。网格是死网格还是活网格。断层走向和断层传导率。 3。流体PVT属性:油,气,水的地面密度或重度;油,气的地层体积系数,粘度随压力变化表;溶解油气比随压力的变化表;水的粘度,体积系数,压缩系数;岩石压缩系数。如果是组分模型,需要提供状态方程。 4。岩石属性:相对渗透率曲线和毛管压力曲线。如果是油,气,水三相,需要提供油水,油气相对渗透率曲线和毛管压力曲线(软件会自动计算三相流动时的相对渗透率曲线);如果是油,水两相或气,水两相,只需要提供油水或气水两相相对渗透率曲线和毛管压力曲线。5。油藏分区参数:如果所模拟的油田横向或纵向流体属性,岩性变化比较大,或者存在不同的油水界面,这时需要对模型进行PVT分区(不同区域用不同的PVT流体参数表),岩石分区(不同区域用不同的相对渗透率曲线和毛管压力曲线)或者平衡分区(不同平衡区用不同的油水界面)。另外如果想掌握油藏不同断块的储量或采收率,可以对模型进行储量分区(不同储量区可以输出不同的储量,产量,采收率,剩余储量等)。 6。初始化计算参数:油藏模型初始化即计算油藏模型初始饱和度,压力和油气比的分布,从而得到油藏模型的初始储量。这部分需要输入模型参考深度,参考深度处对应的初始压力,油水界面以及气水界面;油气比或饱和压力随深度的变化;如果是组分模型,需要输入组分随深度的变化。

某空调间fluent模拟

空调间内空气流场的模拟 模拟课题名称:空调间内空气流场的模拟 模拟的步骤: 1:gambit建立空间模型及建立网格 1.1打开gambit,选择文件的存储路径。 为了寻找的方便,我们将文件尊放在d盘fluent文件下 点击run,进入gambit界面; 1.2建立空调间的空间模型 点击体控制面板中的体控制按钮直接生成六面体,选择direction为+x+y+z方向输入房间尺寸,生成房间如图

同样的方法生成空调的六面体,并将六面体通过移动的选项将六面体移动到正确的位置上如图: 因为整个的计算区域为空调间减去空调的空间所以选择选项,将空调的空间去除,经过渲染之后得到下图: 这样计算区域的几何模型就建立起来了。 1.3网格的划分

这里我们直接对体进行网格划分点击,进入网格划分界面,选择该体,采用六面体的网格单元和submap网格划分方法进行划分如图 Interval size输入0.05,建立网格如图 如图,我们可以看到总共建立了359840个网格。接下来需要对网格单元记性检 查,经过检查发现网格符合要求。 1.4定义边界类型

本文定义房间的上下左右先后分别为壁面边界,空调的下边界为速度进口,空调的上边界为速度出口空调的其他面为壁面边界。建立的过程如图 建立完成后得到:

1.5保存文件后安全退出 执行File Export Mash,输出三维的数据文件,不要在2维的地方打勾,为 了以后查找的方便我们将文件的名字改为kongtiaojian,最后执行Exit安全退出。这里注意不要点右上角的,要选择exit键才能安全退出,否则生成的是 碎片文件。 2 fluent模拟 2.1文件的读取和检查 打开fluent三维单精度计算器如图,执行Flie Read Case,读入我们已经划分好的网格文件,当出现 表示网格文件读入完毕。执行Grid Check 检查网格,表示网格检查完毕。执行Grid Scale我们可

相关文档
最新文档