文档库 最新最全的文档下载
当前位置:文档库 › 用Fluent模拟平均停留时间分布

用Fluent模拟平均停留时间分布

用Fluent模拟平均停留时间分布
用Fluent模拟平均停留时间分布

fluent 软件介绍

百科名片 Fluent是目前国际上比较流行的商用CFD软件包,在美国的市场占有率为60%,凡是和流体、热传递和化学反应等有关的工业均可使用。它具有丰富的物理模型、先进的数值方法和强大的前后处理功能,在航空航天、汽车设计、石油天然气和涡轮机设计等方面都有着广泛的应用。 简介 Fluent算例 CFD商业软件FLUENT,是通用CFD软件包,用来模拟从不可压缩到高度可压缩范围内的复杂流动。由于采用了多种求解方法和多重网格加速收敛技术,因而FLUENT能达到最佳的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网格技术及成熟的物理模型,使FLUENT在转换与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面有广泛应用。 基本特点 FLUENT软件具有以下特点: FLUENT软件采用基于完全非结构化网格的有限体积法,而且具有基于网格节点和网格单元的梯度算法; 定常/非定常流动模拟,而且新增快速非定常模拟功能; Fluent 前处理网格划分 FLUENT软件中的动/变形网格技术主要解决边界运动的问题,用户只需指定初始网格和运动壁面的边界条件,余下的网格变化完全由解算器自动生成。网格变形方式有三种:弹簧压缩式、动态铺层式以及局部网格重生式。其局部网格重生式是FLUENT所独有的,而

且用途广泛,可用于非结构网格、变形较大问题以及物体运动规律事先不知道而完全由流动所产生的力所决定的问题; FLUENT软件具有强大的网格支持能力,支持界面不连续的网格、混合网格、动/变形网格以及滑动网格等。值得强调的是,FLUENT软件还拥有多种基于解的网格的自适应、动态自适应技术以及动网格与网格动态自适应相结合的技术; FLUENT软件包含三种算法:非耦合隐式算法、耦合显式算法、耦合隐式算法,是商用软件中最多的; FLUENT软件包含丰富而先进的物理模型,使得用户能够精确地模拟无粘流、层流、湍流。湍流模型包含Spalart-Allmaras模型、k-ω模型组、k-ε模型组、雷诺应力模型(RSM)组、大涡模拟模型(LES)组以及最新的分离涡模拟(DES)和V2F模型等。另外用户还可以定制或添加自己的湍流模型; 适用于牛顿流体、非牛顿流体; 含有强制/自然/混合对流的热传导,固体/流体的热传导、辐射; 化学组份的混合/反应; 自由表面流模型,欧拉多相流模型,混合多相流模型,颗粒相模型,空穴两相流模型,湿蒸汽模型; 融化溶化/凝固;蒸发/冷凝相变模型; 离散相的拉格朗日跟踪计算; 非均质渗透性、惯性阻抗、固体热传导,多孔介质模型(考虑多孔介质压力突变); 风扇,散热器,以热交换器为对象的集中参数模型; 惯性或非惯性坐标系,复数基准坐标系及滑移网格; 动静翼相互作用模型化后的接续界面; 基于精细流场解算的预测流体噪声的声学模型; 质量、动量、热、化学组份的体积源项; 丰富的物性参数的数据库; 磁流体模块主要模拟电磁场和导电流体之间的相互作用问题; 连续纤维模块主要模拟纤维和气体流动之间的动量、质量以及热的交换问题; 高效率的并行计算功能,提供多种自动/手动分区算法;内置MPI并行机制大幅度提高并行效率。另外,FLUENT特有动态负载平衡功能,确保全局高效并行计算; FLUENT软件提供了友好的用户界面,并为用户提供了二次开发接口(UDF); FLUENT软件采用C/C++语言编写,从而大大提高了对计算机内存的利用率。 在CFD软件中,Fluent软件是目前国内外使用最多、最流行的商业软件之一。Fluent 的软件设计基于"CFD计算机软件群的概念",针对每一种流动的物理问题的特点,采用适合于它的数值解法在计算速度、稳定性和精度等各方面达到最佳。由于囊括了Fluent Dynamical International比利时PolyFlow和Fluent Dynamical International(FDI)的全部技术力量(前者是公认的在黏弹性和聚合物流动模拟方面占领先地位的公司,后者是基于有限元方法CFD软件方面领先的公司),因此Fluent具有以上软件的许优点 软件简介

停留时间分布的测量

停留时间分布的测量 1 固体粒子停留时间分布测量 固体粒子停留时间测量方法有间接法和直接法两种.间接的测量是基于总的固体相速度和相分率,而直接法大多借助示踪剂进行测量. 采用示踪剂测量时.除要求示踪剂具有与被测体系有相同的流动行为外.还要求具有可辨别的其他物理或化学性质.诸如荧光性.导电性.红外或介电性等.最常用的示踪剂是颜色示踪剂.化学示踪剂.磁性示踪剂.放射性示踪剂等. (1)颜色示踪剂方法:对于具有透明壁的设备.颜色示踪剂无疑是简单易行的方法.对于多孔粒子.只需将颜料包埋到粒子中即可.对于非多孔粒子.最简单的方法是将颜料涂在表面.但由于表面的颜料容易脱落.影响测量精度.为此.最好能将颜料均匀地混到粒子中.当然无论是多孔或非多孔体系.颜料须与液相不相溶. 这种方法的优点是形象.直接.其缺点或局限性是要求设备具有透明壁和良好的观察设备.如遥控可移动的摄像机及其良好的光学系统,如果示踪剂被非示踪剂粒子包围.则无法检测到.从面影响精度,同时能够观察到的主要是壁附近的信息.有时难以代表整体的行为,实验后的粒子重新分拣通常是困难的.所以物料难以重复使用.废料的处理也是个问题. (2)化学示踪剂方法:这种示踪剂的化学组成是不同的.所有其他的物理性质与非示踪剂粒子相同.示踪剂必须不与体系中的其他组分发生反应.取样必须有代表性. 2 液相停留时间分布测量 液相停留时间分布常用的测量方法有染料示踪剂测量法.电导示踪法.折射指数法.放射性法和热示踪剂法等. 在采用示踪剂测量方法时.一般要求示踪剂的流动性质.特别是密度.粘度.界面张力和互溶性等必须尽可能接近于被测试介质.同时为了保证相似性.对于水溶性体系建议用水溶性电解质作为示踪剂.对于有机体系则采用可溶于有机物的物质作示踪剂. 如果混合体系中含有多孔固体.或者容器内壁.或内部构件易于吸收示踪剂.那么一些示踪剂就会被吸附到这些材料上.从而增大测量难度.此时就要设法选择一种不易被吸附的材料作示踪剂. (1)染料示踪剂测量法:这种测量方法的原理很简单.首先将示踪剂在进口区注入.如果釜壁是透明的.可以目测示踪剂的运动轨迹.通常是将相应的检测器置于进口和出口处.根据两处浓度的变化曲线计算停留时间分布.这种方法的优点是测量直接.迅速.成本低,缺点是对于非常深色的体系或重油的体系.难以找到合适的染料. (2)电导示踪法:这种方法通常用于被测体系是非电导的.面引入电导介质作示踪剂.电导率随浓度的不同而变化.常用的电解质有NaCl.KCl等.该方法具有简单.迅速.易于使用等优点,但一般只用于水溶性体系. (3)折射指数法:这种方法引入了折射指数.明显不同于被测体系的介质作为示踪剂的方法.测量原理类似于前面的方法.该方法使用简单.响应迅速,但也只适用于透明体系. 3 气相停留时间分布测量 气体停留时间的测量非常类似干前面用于固体粒子和液相体系的测量方法.其中化学示踪剂法是使用最广泛的一种测量方法.常用的示踪气体是氦.

实验一多釜串联连续流动反应器中停留时间分布的测定

实验一 多釜串联连续流动反应器中停留时间分布的测定 一、实验目的 本实验通过单釜与三釜反应器中停留时间分布的测定,将数据计算结果用多釜串联模型来描述返混程度,从而认识限制返混的措施。 1、掌握停留时间分布的测定方法; 2、了解停留时间分布与多釜串联模型的关系; 3、掌握多釜串联模型参数N 的物理意义及计算方法。 二、实验原理 在连续流动的反应器内,不同停留时间的物料之间的混和称为返混。返混程度的大小,一般很难直接测定,通常是利用物料停留时间分布的测定来研究。然而在测定不同状态的反应器内停留时间分布时,可以发现,相同的停留时间分布可以有不同的返混情况,即返混与停留时间分布不存在一一对应的关系,因此不能用停留时间分布的实验测定数据直接表示返混程度,而必须借助于反应器数学模型来间接表达。 物料在反应器内的停留时间完全是一个随机过程,须用概率分布方法来定量描述。所用的概率分布函数为停留时间分布密度函数E (t)和停留时间分布函数F (t)。 停留时间分布密度函数E (t )的物理意义是:同时进入的N 个流体粒子中,停留时间介于t 到t +dt 间的流体粒子所占的分率dN/N 为E (t )dt 。 停留时间分布函数F (t )的物理意义是:流过系统的物料中停留时间小于t 的物料所占的分率。 停留时间分布的测定方法有脉冲输入法、阶跃输入法等,常用的是脉冲输入法。当系统达到稳定后,在系统的入口处瞬间注入一定量Q 的示踪物料,同时开始在出口流体中检测示踪物料的浓度变化。 由停留时间分布密度函数的物理含义,可知: E (t )dt =VC (t )/Q (1) ?∞ =0)(dt t VC Q (2) 所以 ? ? ∞ ∞ = = )() ()() ()(dt t C t C dt t VC t VC t E (3)

FLUENT算例 (5)搅拌桨底部十字挡板的流场分析

搅拌桨底部十字挡板的流场分析搅拌设备在各个行业运用的十分广泛,搅拌就是为了更够更快速更高效的将物质与介质充分混合,发生充分的反应,而搅拌中存在着许多不利于混合的情况,比如液体旋流。为了解决这个问题,之前很多人提出在罐体的侧壁上增加挡板,可以抵消大部分旋流,然后大部分都是研究侧挡板的,对于底部挡板的研究十分少,本文就在椭圆底部挡板增加十字型挡板,对罐体中进行流场分析。 1.Gambit建模 首先用Gambit建模图形如下: 图1:Gambit建立的模型 分为两个区域,里面的圆柱为动区域,外面包着的大圆柱设为静区域,静区域划分网格大,划分粗糙,内部动区域划分网格小,划分精细。边界条件主要设置了轴,搅拌桨,底部挡板,上层液面。以下就是fluent进行数值模拟。 2.fluent数值模拟 2.1导入case文件

2.2对网格进行检查 Minimum volume的数值大于0即可。 图2网格检查2.3调节比例 单位选择mm单位。 图3比例调节2.4定义求解器参数 设置如图4所示

图4设置求解器参数2.5设置能量线 图5能量线 2.6设置粘度模型,选择k-e模型 k-e模型对该模型模拟十分实用。

图6粘度模型2.7定义材料 介质选择液体水。 2.8定义操作条件

由于存在着终于,建模时的方向向上,所以在Z轴增加一个重力加速度。 图8操作条件 2.9定义边界条件 在边界设置重,动区域如图所示,将材料设成水,motion type设成moving reference frame (相对滑动),转速设为10rad/s,单位可在Define中的set unit中的angular-velocity设置。而在在轴的设置中,如上图所示,将wall motion设成moving wall,motion设成Absolute,速度设成-10,由于轴跟动区域速度是相对的,所以设成反的。

实验1连续搅拌釜式反应器停留时间分布的测定

实验一 连续搅拌釜式反应器停留时间分布的测定 一、 实验目的 (1) 加深对停留时间分布概念的理解; (2) 掌握测定液相停留时间分布的方法; (3) 了解停留时间分布曲线的应用。 (4)了解停留时间分布于多釜串联模型的关系,了解模型参数N 的物理意义及计算方法。 (5) 了解物料流速及搅拌转速对停留时间分布的影响。 二、 实验原理 (1)停留时间分布 当物料连续流经反应器时,停留时间及停留时间分布是重要概念。停留时间分布和流动模型密切相关。流动模型分平推流,全混流与非理想流动三种类型。 对于平推流,流体各质点在反应器内的停留时间均相等,对于全混流,流体各质点在反应器内的停留时间是不一的,在0~∞范围内变化。对于非理想流动,流体各质点在反应器内的停留时间分布情况介乎于以上两种理想状态之间,总之,无论流动类型如何,都存在停留时间分布与停留时间分布的定量描述问题。 (2)停留时间分布密度函数E (t ) 停留时间分布密度函数E (t )的定义: 当物料以稳定流速流入设备(但不发生化学变化)时,在时间t =0时,于瞬时间dt 进入设备的N 个流体微元中,具有停留时间为t 到(t +dt )之间的流体微元量dN 占当初流入量N 的分率为E (t )dt ,即 ()=dN E t dt N (1) E (t )定义为停留时间分布密度函数。 由于讨论的前提是稳定流动系统,因此,在不同瞬间同时进入系统的各批N 个流体微元均具有相同的停留时间分布密度,显然,流过系统的全部流体,物料停留时间分布密度为同一个E (t )所确定。根据E (t )定义,它必然具有归一化性质:

()1∞ =? E t dt (2) 不同流动类型的E (t )曲线形状如图1所示。根据E (t )曲线形状,可以定性分析物料在反应器(设备)内停留时间分布。 平推流 全混流 非理想流动 图1 各种流动的E (t )~t 关系曲线图 (3)停留时间分布密度函数E (t )的测定 停留时间分布密度函数E (t )的测定,常用的方法是脉冲法。此法采用的示踪剂,既不与被测流体发生化学反应,又不影响流体流动特性,也就是说,示踪物在反应器(设备)内的停留时间分布与被测流体的停留时间分布相同。所以,当注入一定量Q 的示踪物时,经过t →(t +dt )时间间隔流出的示踪物量占示踪物注入总量Q 的分率就是与示踪物注入同时进入系统的物料中,停留时间为t →(t +dt )的那部分流体物料占总流体的物料的分率, 即: 亦即: ()()??=V C t dt E t dt Q 或 () ()?= V C t E t Q (3) V ——流体体积流量,(ml/s) Q ——加入的示踪物总量,(mg) C (t )——示踪物的出口浓度,(mg/ml)

FLUENT算例 (9)模拟燃烧

计算流体力学作业FLUENT 模拟燃烧 问题描述:长为2m、直径为0.45m的圆筒形燃烧器结构如图1所示,燃烧筒壁上嵌有三块厚为0.0005 m,高0.05 m的薄板,以利于甲烷与空气的混合。燃烧火焰为湍流扩散火焰。在燃烧器中心有一个直径为0.01 m、长为0.01 m、壁厚为0.002 m的小喷嘴,甲烷以60 m/s的速度从小喷嘴注入燃烧器。空气从喷嘴周围以0.5 m/s的速度进入燃烧器。总当量比大约是0.76(甲烷含量超过空气约28%),甲烷气体在燃烧器中高速流动,并与低速流动的空气混合,基于甲烷喷嘴直径的雷诺数约为5.7×103。 假定燃料完全燃烧并转换为:CH4+2O2→CO2+2H2O 反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。利用FLUENT的finite-rate化学反应模型对一个圆筒形燃烧器内的甲烷和空气的混合物的流动和燃烧过程进行研究。 1、建立物理模型,选择材料属性,定义带化学组分混合与反应的湍流流动边界条件 2、使用非耦合求解器求解燃烧问题 3、对燃烧组分的比热分别为常量和变量的情况进行计算,并比较其结果 4、利用分布云图检查反应流的计算结果 5、预测热力型和快速型的NO X含量 6、使用场函数计算器进行NO含量计算 一、利用GAMBIT建立计算模型 第1步启动GAMBIT,建立基本结构 分析:圆筒燃烧器是一个轴对称的结构,可简化为二维流动,故只要建立轴对称面上的

二维结构就可以了,几何结构如图2所示。 (1)建立新文件夹 在F盘根目录下建立一个名为combustion的文件夹。 (2)启动GAMBIT (3)创建对称轴 ①创建两端点。A(0,0,0),B(2,0,0) ②将两端点连成线 (4)创建小喷嘴及空气进口边界 ①创建C、D、E、F、G点

停留时间分布与反应器答案(李绍芬)

5 停留时间分布与反应器 5.1设F(θ)及E(θ)分别为闭式流动反应器的停留时间分布函数及停留时间分布密度函数,θ为对比时间。 (1) (1) 若该反应器为活塞流反应器,试求 (a ) (a ) F(1)(b)E(1)(c)F(0.8)(d)E(0.8)(e)E(1.2) (2)若该反应器为全混流反应器,试求 (a )F(1)(b)E(1)(c)F(0.8)(d)E(0.8)(e)E(1.2) (3) 若该反应器为一个非理想流动反应器,试求 (a )F(∞)(b)F(0)(c)E(∞)(d)E(0)(e)0∞ ? ()E d θ? (f)0∞ ? ()E d θθ? 解:(1)因是活塞流反应器,故符合理想活塞流模型的停留时间分布,由(5.33-5.36)式可得: (a)F(1)=1.0(b)E(1)=∝(c)F(0.8)=0(d)E(0.8)=0(e)E(1.2)=0 (2) (2) 因是全混流反应器,故符合理想全混流模型的停留时间分布,由 (5.33-5.36)式可得: (a )F(1)=1-e -1=0.6321 (b)E(1)=e -1=0.3679 (c)F(0.8)=1- e -0.8=0.5507 (d)E(0.8)= e -0.8=0.4493 (e)=E(1.2)=0.3012 (3) (3) 因是一个非理想流动反应器,故可得: (a )F(∞)=1 (b)F(0)=0 (c)E(∞)=0 (d)E(0)>1 (e)0 ∞?()E d θ? =1 (f) 0 ∞ ? ()E d θθ? =1 5.2用阶跃法测定一闭式流动反应器的停留时间分布,得到离开反应器的示踪剂与时间的关系如下: 022 2313 ≤??=-≤≤??≥? ()t c t t t t 试求: (1) (1) 该反应器的停留时间分布函数F(θ)及分布密度函数E(θ)。 (2) (2) 数学期望θ及方差2 θσ。

化工专业实验:停留时间分布的测定

停留时间分布的测定 一、实验目的 1. 了解利用电导率测定停留时间分布的基本原理和实验方法; 2. 掌握停留时间分布的统计特征值的计算方法; 3. 掌握用理想反应器串联模型来描述实验系统的流动特性。 二、实验原理 停留时间分布测定所采用的方法主要是示踪响应法。它的基本思路是:在反应器入口以一定的方式加入示踪剂,然后通过测量反应器出口处示踪剂浓度的变化,间接地描述反应器内流体的停留时间。常用的示踪剂加入的方式有脉冲输入、阶跃输入和周期输入等。本实验选用的是脉冲输入法。 脉冲输入法是在极短的时间内,将示踪剂从系统的入口处注入主流体,在不影响主流体原有流动特性的情况下随之进入反应器。与此同时,在反应器出口检测示踪剂浓度c(t)随时间的变化。整个过程可以用图1形象地描述。 示踪剂加入示踪剂检测 Q∞Q (a) 脉冲输入法 C C t=0 t (b) 脉冲输入(c) 出口响应 图1 脉冲法测停留时间分布 脉冲输入法测得的停留时间分布代表了物料在反应器中的停留时间分布密度即E(t)。若加入示踪剂后混合流体的流率为Q,出口处示踪剂浓度为C(t),在dt时间里示踪剂的流出量为Qc(t)dt,由E(t)定义知E(t)dt是出口物料中停留时间在t与t+dt之间示踪剂所占分率,若在反应器入口加入示踪剂总量为m 对反应器出口作示踪剂的物料衡算,即 (1) 示踪剂的加入量可以用下式计算 (2) 在Q值不变的情况下,由(1)式和(2)式求出:

(3) 关于停留时间的另一个统计函数是停留时间分布函数F(t),即 (4) 用停留时间分布密度函数E(t)和停留时间分布函数F(t)来描述系统的停留时间,给出了 很好的统计分布规律。但是为了比较不同停留时间分布之间的差异,还需要引入另外两个统 计特征值,即数学期望和方差。 数学期望对停留时间分布而言就是平均停留时间_ t ,即 (5) 方差是和理想反应器模型关系密切的参数,它的定义是: (6) 若采用无因次方差2Θσ则有22 2/t t σσΘ=对活塞流反应器02=Θσ,而对全混流反应器 12=Θσ;对介于上述两种理想反应器之间的非理想反应器可以用多釜串联模型描述。多釜串 联模型中的模型参数N可以由实验数据处理得到的 2Θσ来计算。 21N σΘ= (7) 当N为整数时,代表该非理想流动反应器可以用N个等体积的全混流反应器的串联来建 立模型。当N为非整数时,可以用四舍五入的方法近似处理,也可以用不等体积的全混流反 应器串联模型。 三、实验设备的特点 1. 本实验装置数据稳定,重现性好, 使用方便,安全可靠。 2. 本装置设备紧凑,功能齐全。 四、装置、流程及试剂

大涡模拟的FLUENT算例2D

Tutorial:Modeling Aeroacoustics for a Helmholtz Resonator Using the Direct Method(CAA) Introduction The purpose of this tutorial is to provide guidelines and recommendations for the basic setup and solution procedure for a typical aeroacoustic application using computational aeroacoustic(CAA)method. In this tutorial you will learn how to: ?Model a Helmholtz resonator. ?Use the transient k-epsilon model and the large eddy simulation(LES)model for aeroacoustic application. ?Set up,run,and perform postprocessing in FLUENT. Prerequisites This tutorial assumes that you are familiar with the user interface,basic setup and solution procedures in FLUENT.This tutorial does not cover mechanics of using acoustics model,but focuses on setting up the problem for Helmholtz-Resonator and solving it.It also assumes that you have basic understanding of aeroacoustic physics. If you have not used FLUENT before,it would be helpful to?rst review FLUENT6.3User’s Guide and FLUENT6.3Tutorial Guide. Problem Description A Helmholtz resonator consists of a cavity in a rigid structure that communicates through a narrow neck or slit to the outside air.The frequency of resonance is determined by the mass of air in the neck resonating in conjunction with the compliance of the air in the cavity. The physics behind the Helmholtz resonator is similar to wind noise applications like sun roof bu?eting. We assume that out of the two cavities that are present,smaller one is the resonator.The motion of the?uid takes place because of the inlet velocity of27.78m/s(100km/h).The ?ow separates into a highly unsteady motion from the opening to the small cavity.This unsteady motion leads to a pressure?uctuations.Two monitor points(Point-1and Point-2) act as microphone points to record the generated sound.The acoustic signal is calculated within FLUENT.The?ow exits the domain through the pressure outlet.

连续均相反应器停留时间分布的测定实验报告

实验日期 2015.5.29 成绩 同组人×××(2)、×××(3)、×××(4)、×××(5)、×××(6)闽南师范大学应用化学专业实验报告题目:连续均相反应器停留时间分布的测定 12应化1 ××12060001××B1组 0 前言 实验目的:1,、了解管式反应器的特点和原理;2、掌握脉冲示踪法测定管式反应器和釜式反应器内物料停留时间分布测定和数据处理方法;3、掌握活塞管式反应器和全混流反应器内物料的停留时间分布密度函数E(t)和停留时间分布函数F(t)的特点及其数学特征;4、学会用理想反应器的串联模型来描述实验的流动性。[1] 实验原理:由于反应器内流体速度分布不均匀,或某些流体微元运动方向与主体流动方向相反,因此使反应器内流体流动产生不同程度的返混。在反应器设计、放大和操作时,往往需要知道反应器中返混程度的大小。停留时间分布能定量描述返混程度的大小,而且能够直接测定。因此停留时间分布测定技术在化学反应工程领域中有一定的地位。 停留时间分布可用分布函数F(t)和分布密度E(t)来表示,两者的关系为: 测定停留时间分布最常用的方法是阶跃示踪法和脉冲示踪法。阶跃法: 脉冲法: 式中:C(t)——示踪剂的出口浓度。Co——示踪剂的入口浓度。[2] Vs———流体的流量Qλ——示踪剂的注入量。

由此可见,若采用阶跃示踪法,则测定出口示踪物浓度变化,即可得到F(t)函数;而采用脉冲示踪法,则测定出口示踪物浓度变化,就可得到E(t)函数。 1 实验方案 1.1 实验材料 三釜串联反应器 1.2 实验流程与步骤 实验流程图: 实验步骤: (1)准备工作 ②饱和KNO3液体注入标有KNO3的储瓶内。 ②连接好入水管线,打开自来水阀门,使管路充满水。 ③查电极导线是否正确。 (2)操作 ①打开总电源开关,,开启水阀门,向朝内注水,将回流阀开到最大,启动水泵,慢慢打开进水转子阀门通过“管式/釜式”阀门转向“釜式”一侧,直管坏了,只做三釜串联反应实验,调节水流量维持在30-40L/h之间的某值,直至釜1内有少量水,并能正常流出。 ②开启电磁阀开关和电导仪总开关。 ③“示踪剂/清洗水”阀门转向示踪剂一侧。

大涡模拟的fluent算例

Introduction:This tutorial demonstrates how to model the2D turbu-lent?ow across a circular cylinder using LES(Large Eddy Simula-tion),and compute?ow-induced noise(aero-noise)using FLUENT’s acoustics model. In this tutorial you will learn how to: ?Perform2D Large Eddy Simulation(LES) ?Set parameters for an aero-noise calculation ?Save surface pressure data for an aero-noise calculation ?Calculate aero-noise quantities ?Postprocess an aero-noise solution Prerequisites:This tutorial assumes that you are familiar with the menu structure in FLUENT,and that you have solved or read Tu-torial1.Some steps in the setup and solution procedure will not be shown explicitly. Problem Description:The problem considers turbulent air?ow over a2D circular cylinder at a free stream velocity U of69.19m/s. The cylinder diameter D is1.9cm.The Reynolds number based on the?ow parameters is about90000.The computational do-main(Figure3.0.1)extends5D upstream and20D downstream of the cylinder,and5D on both sides of it.If the computational domain is not taken wide enough on the downstream side,so that no reversed?ow occurs,the accuracy of the aero-noise prediction may be a?ected.The rule of thumb is to take at least20D on the downstream side of the obstacle. c Fluent Inc.June20,20023-1

FLUENT算例 (3)三维圆管紊流流动状况的数值模拟分析

三维圆管紊流流动状况的数值模拟分析 在工程和生活中,圆管内的流动是最常见也是最简单的一种流动,圆管流动有层流和紊流两种流动状况。层流,即液体质点作有序的线状运动,彼此互不混掺的流动;紊流,即液体质点流动的轨迹极为紊乱,质点相互掺混、碰撞的流动。雷诺数是判别流体流动状态的准则数。本研究用CFD 软件来模拟研究三维圆管的紊流流动状况,主要对流速分布和压强分布作出分析。 1 物理模型 三维圆管长2000mm l =,直径100mm d =。 流体介质:水,其运动粘度系数6 2 110m /s ν-=?。 Inlet :流速入口,10.005m /s υ=,20.1m /s υ= Outlet :压强出口 Wall :光滑壁面,无滑移 2 在ICEM CFD 中建立模型 2.1 首先建立三维圆管的几何模型Geometry 2.2 做Blocking 因为截面为圆形,故需做“O ”型网格。

2.3 划分网格mesh 注意检查网格质量。 在未加密的情况下,网格质量不是很好,如下图 因管流存在边界层,故需对边界进行加密,网格质量有所提升,如下图

2.4 生成非结构化网格,输出fluent.msh等相关文件 3 数值模拟原理 紊流流动

当以水流以流速20.1m /s υ=,从Inlet 方向流入圆管,可计算出雷诺数10000υd Re ν ==,故圆管内流动为紊流。 假设水的粘性为常数(运动粘度系数62 110m /s ν-=?)、不可压流体,圆管光滑,则流动的控制方程如下: ①质量守恒方程: ()()()0u v w t x y z ρρρρ????+++=???? (0-1) ②动量守恒方程: 2()()()()()()()()()()[]u uu uv uw u u u t x y z x x y y z z u u v u w p x y z x ρρρρμμμρρρ??????????+++=++??????????'''''????+---- ???? (0-2) 2 ()()()()()()()()()()[]v vu vv vw v v v t x y z x x y y z z u v v v w p x y z y ρρρρμμμρρρ??????????+++=++??????????'''''????+- ---???? (0-3) 2 ()()()()()()()()()()[]w wu wv ww w w w t x y z x x y y z z u w v w w p x y z z ρρρρμμμρρρ??????????+++=++??????????'''''????+- ---???? (0-4) ③湍动能方程: ()()()()[())][())][())]t t k k t k k k ku kv kw k k t x y z x x y y k G z z μμρρρρμμσσμμρεσ????????+++=+++????????? ?+ ++-?? (0-5) ④湍能耗散率方程: 212()()()()[())][())][())]t t k k t k k u v w t x y z x x y y C G C z z k k εεμμρερερερεεεμμσσμεεεμρσ??????? ?+++=+++??????????+++-?? (0-6) 式中,ρ为密度,u 、ν、w 是流速矢量在x 、y 和z 方向的分量,p 为流体微元体上的压强。 方程求解:采用双精度求解器,定常流动,标准ε-k 模型,SIMPLEC 算法。 4 在FLUENT 中求解计算紊流流动 4.1 FLUENT 设置 除以下设置为紊流所必须设置的外,其余选项和层流相同,不再详述。

实验一 填充管式反应器液体停留时间分布及其流动模型参数的测定

实验一 填充管式反应器 液体停留时间分布及其流动模型参数的测定 一、实验目的 当流体(气体或液体)流经填充层进行均相反应,或者流体通过固体颗粒层(固定床)进行非均相反应或非均相催化反应时,由于各种原因造成流体质点在反应器内停留时间不一,形成不同的停留时间分布。不同的停留时间分布直接影响反应结果,如反应的最终转化不同。填充管式反应器或固定床反应器均可视为连续流动的管式反应器,其理想流动模型为活塞流模型。这类反应器的理想流动模型能够的检验,实现理想流动的边值操作条件的确定,以及非理想流动反应器的流动模型和模型参数的确定,都应先通过实验测定流体流经反应器的停留 时间分布。 停留时间分布的实验测定方法通常用两种方法:脉冲激发——响应法和阶跃激发——响应法。本实验以水为主体流体,以氯化钾饱和溶液为示踪剂,采用脉冲输入的方法测定流体流经填充层或固定床层的停留时间分布。这种方法不仅用于检验或确定填充管式均相反应器和固定床均相反应器的流动模型,也适用于填料塔等传质设备。 通过本实验掌握一种测定停留时间分布的实验技术,同时初步掌握对流体流经固体颗粒层这类是设备的流动模型检验和模型参数的实验测定方法。毫无疑问。通过实验对于数学模 型方法和流动模型等方面的有关概念,原理和方法会有更深入的理解。 二、实验原理 采用脉冲激发——响应法测定停留时间分布的实验方法,是当主流流体以恒定的体积流率流经具有一定堆积的填充层时,在反应器如口出瞬时脉冲注入一定量的示踪剂,与此同时在反应器出口处检测示踪物浓度与时间的关系曲线,即t t c -)(曲线,并可转化为停留时间 分布密度与时间的关系曲线,即t t E -)(曲线。 由停留时间分布实验曲线可以定性地诊断流体流经反应器的流动状况。停留时间分布属于随机变量的分布,概率上还可以定量地用数字特征加以描述,表征这种随机分布的数字特 征主要是数学期望和方差。 (1) 停留时间分布的数学期望,t ) 随机变量的数学期望也就是该变量的平均数。流体流经反应器停留时间分布的数学期 望也就是停留时间的平均值。停留时间分布数学期望的定义式为 ? ?∞∞ = 0)()(dt t E dt t tE t ) (1)

FLUENT算例(9)模拟燃烧(优选.)

最新文件---------------- 仅供参考--------------------已改成-----------word文本 --------------------- 方便更改 计算流体力学作业FLUENT 模拟燃烧 问题描述:长为2m、直径为0.45m的圆筒形燃烧器结构如图1所示,燃烧筒壁上嵌有三块厚为0.0005 m,高0.05 m的薄板,以利于甲烷与空气的混合。燃烧火焰为湍流扩散火焰。在燃烧器中心有一个直径为0.01 m、长为0.01 m、壁厚为0.002 m的小喷嘴,甲烷以60 m/s的速度从小喷嘴注入燃烧器。空气从喷嘴周围以0.5 m/s的速度进入燃烧器。总当量比大约是0.76(甲烷含量超过空气约28%),甲烷气体在燃烧器中高速流动,并与低速流动的空气混合,基于甲烷喷嘴直径的雷诺数约为5.7×103。 假定燃料完全燃烧并转换为:CH4+2O2→CO2+2H2O 反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。利用FLUENT的finite-rate化学反应模型对一个圆筒形燃烧器内的甲烷和空气的混合物的流动和燃烧过程进行研究。 1、建立物理模型,选择材料属性,定义带化学组分混合与反应的湍流流动边界条件 2、使用非耦合求解器求解燃烧问题 3、对燃烧组分的比热分别为常量和变量的情况进行计算,并比较其结果 4、利用分布云图检查反应流的计算结果 5、预测热力型和快速型的NO X含量 6、使用场函数计算器进行NO含量计算 一、利用GAMBIT建立计算模型 第1步启动GAMBIT,建立基本结构

分析:圆筒燃烧器是一个轴对称的结构,可简化为二维流动,故只要建立轴对称面上的二维结构就可以了,几何结构如图2所示。 (1)建立新文件夹 在F盘根目录下建立一个名为combustion的文件夹。 (2)启动GAMBIT (3)创建对称轴 ①创建两端点。A(0,0,0),B(2,0,0) ②将两端点连成线 (4)创建小喷嘴及空气进口边界

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、三维模型模拟; 五、

(完整版)如何在gambit中提高网格质量

如何在gambit中提高网格质量 经常在网上看到一些网友为gambit划分不出好的网格质量而烦恼。 要生成一套好的网格,我觉得以下几点是很必要的: 1.选择一款好的网格生成软件; 2.确保实体尽量简洁; 3.合理布置线上节点; 但是,对于一些初学者来说,gridgen等专业点的网格划分软件在短时间内是很难掌握的,所以大部分人还是喜欢用gambit。对于gambit来说,有的时候满足了条件2,3,仍然有可能生成质量很差的网格,这个时候就需要手动调整以提高网格质量了。下面我将以一个例子来详细讲解一下如何在gambit中提高网格质量。 这个是个简单的楔形体,包括附面层网格。该网格满足实体简单,节点的布置也合理,但是生成的网格质量很差,主要是在楔形体尾部附面层网格与三角形网格交接的地方。 该图为放大图,从中可以看出有一个网格基本上已经退化成一条线了,从而导致整个网格最大的倾斜率超过了0.99。

解决方法一: 由于质量差的网格集中在附面层与三角形网格过渡的地方,可以从改变附面层网格分布入手。 改变楔形体三个顶点的类型,将其改为side,从而改变附面层网格。 改变附面层网格分布后,重新生成的网格质量提高了不少。 解决方法二: 改变三角形网格分布。

选择调整面网格的节点分布。 手动调整质量差的网格的节点,使其分布合理。

通过调整后,最大倾斜率小于0.91了。该质量的网格基本上就能导入fluent计算了,通过fluent中的smooth/swap功能,还能进一步提高网格质量。 以上例子只是给网友一个在gambit中调整网格的思路,希望能解决一部分人的问题。 其实,提高网格质量最好的办法就是将坏的网格merge到好的网格中,可惜我目前还没有在gambit中发现该功能。有机会再跟大家探讨一下在tgrid中如何用merge功能提高网格质量。

实体入水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所示。

相关文档
相关文档 最新文档