在水平圆管两相流数值模拟中的应用

在水平圆管两相流数值模拟中的应用
在水平圆管两相流数值模拟中的应用

气液两相流

气液两相流流型识别理论的研究进展 摘要:介绍了气液两相流的识别理论,探讨了气液两相流流型的划分方法。叙述了两相流流型软测量方法,并重点介绍了图像处理识别、在线流型技术识别、神经网络、基于压差波动理论、混沌理论等识别流型的新方法。 关键词:气液两相流;流型识别 0 引言 相的概念通常是指某一系统中具有相同成分及相同物理、化学性质的均匀物质成分,各相之间有明显可分的界面。从宏观的角度出发,可以把自然界的物质分为三种,即:气相、液相和固相。单相物质的流动称为单相流,如气体流或液体流。所谓两相流(Two-Phase Flow)或多相流(Multiphase Flow)是指同时存在两种或多种不同相的物质的流动。 近年来随着国内外石油和天然气工业的发展,迫切需要开发出精度较高的油气水三相流量在线测量仪,以便掌握各个油井的生产动态。然而,多年来尽管在这方面进行了大量的研究工作,取得了一些进展,但是仍然没有彻底清晰地认识和了解油气水三相混合物的流动型态。在现今的多相流检测技术领域中,流型的识别问题变得越来越重要。 1 两相流流型 由于存在一个形状和分布在时间和空间里是随机可变的相界面,而相间实际上又存在一个不可忽略的相对速度,致使流经管道的分相流量比和分相所占的管截面比并不相等。这就导致了两相流动结构多种多样,流型十分复杂。流型是影响两相流压力损失和传热特性的重要因素。两相流各种参数的准确测量也往往依赖于对流型的了解。因此为了对两相流的特征参数进行测量,必须了解它们的流型。 1.1垂直上升管中气液两相流流型 (1)、泡状流(Bubbly Flow):气泡以不同尺寸的小气泡形式随机离散分布在流动的液体中。显然,此时气体为离散相,而液体为连续相。随着气速的增加,气泡尺寸会不断增大。 (2)、段塞流(Slug Flow):在气泡流动中当气泡的浓度增高时,气泡聚合为直径接近于管内径的塞状或炮弹状气泡,气泡前端部分呈现为抛物线形状。在这些塞状气泡之间可带有小气泡的液团。当气泡快速上升时,液体在气泡与管内壁间的间隙中流动。 (3)、混状流(Churn Flow):当气泡速度进一步增大时,段塞流中的气泡速度也随之增加并产生破裂、碰撞、聚合和变形,与液体混合成为一种不稳定的上下翻滚的湍动混合物。此时气液两相界为离散相。 (4)、环状流(Annular Flow):液流沿着管道的内壁形成一层液体薄膜,而气流则在管道中央流动。这样,气液两相都变成了连续相。不过,在这种情况下,管道中央的气体通常还夹带着一些液滴一起流动。 (5)、液丝环状流(Wispy-Annular Flow):当气液两相流为环状流时,继续增加液相流量,管壁的液膜将加厚且含有小气泡,中心的液滴浓度增加,被中心

《冶金工程概论》课程大纲

东北大学本科课程教学大纲 课程名称:冶金工程概论 开课单位:材料与冶金学院 制订时间:2004年3月 修订时间:2013年3月

《冶金工程概论》课程教学大纲一.课程基本信息

二.内容结构 基于《冶金工程概论》课程性质,依照东北大学冶金人才培养目标,设计《冶金工程概论》课程

内容,共6章、24学时(其中2学时“职业发展规划”内容,此处未列入),具体分配如下:第一章走进冶金行业(4学时),介绍冶金行业的特点及培养冶金人才知识结构,介绍钢铁生产的现状、最新前沿研究及热点问题,介绍冶金史、历史重要人物及事件。本章的主要内容结构为: 1.1 冶金专业的选择与设置 1.1.1我们为什么选择冶金专业 1.1.2 为何设置冶金专业 1.1.3 合格的冶金工程师是什么样 1.2 怎样走进冶金领域 1.2.1 我们怎样走进冶金领域 1.2.2 在冶金领域我们应该做什么 1.3 学习冶金的任务及目的 1.4 冶金史 1.4.1 冶金工艺的发展历史:过去、现在、将来 1.4.2 我国古代和当代钢铁冶金的地位 第二章钢铁冶金概述(5学时),介绍钢铁冶炼的基本原理、工艺流程、主要设备及新一代钢铁冶金流程。本章的主要内容结构为: 2.1 钢铁冶金流程概述 2.1.1 高炉炼铁-转炉炼钢流程 2.1.2 废钢电炉炼钢流程 2.1.3 非高炉-电炉炼钢流程 2.2 高炉炼铁 2.2.1 高炉炼铁基本任务 2.2.2 高炉炼铁系统 2.2.3 高炉内的主要物理化学过程 2.3 铁水预处理 2.3.1 铁水预处理基本任务 2.3.2 铁水预处理设备及处理剂 2.4 转炉炼钢 2.4.1 转炉炼钢基本任务 2.4.2 转炉内的主要物理化学过程 2.5 电炉炼钢

颗粒流介绍

年后 颗粒流(PFC)简介 2009-10-07 11:14:48| 分类:岩土工程| 标签:|字号大中小订阅 注:今天偶然间见到颗粒流的概念,以前一直不了解,今天查了查,贴在这里,以备以 后可以温故知新。 本文内容源自浙江大学罗永先生的博士论文,使得吾辈能花较少的时间看到广博的知 识,在此特向其表示感谢! 岩土工程数值计算总体上可以分为两大类:一类是基于连续介质力学理论的方法,如有限元法(FEM)和快速拉格朗日法 (FLAC(1tasea,2002))等;另一类是不连续介质力学的方法,如离散元法UDEC(1tasca,2000)、3DEC(Itasea,1998)、PFC(Itasea,2002)和块体理论DDA(石根华,1988)等。离散元方法按其用途又可以分为宏观离散元方法和细观离散元方法,前者主要针对解决规模相对较大的不连续面,如断层节理结构与基础之间的结合面等引起的问题(UDEC,3DEC),后者则着重于数目众多具有不连续特性的接触面或点,如破碎岩体中的破裂面、砂土中的接触面(点)和材料中颗粒之间的接触面(点)等。PFC(Particle Flow Code)是在著名学者Peter Cundall 主持下采用细观离散元理论(又称为粒子流理论)开发的一种数值计算平台,可以广泛地应用于研究细观结构控制问题。目前,PFC在世界上的应用并不广泛,成果报道也主要集中在PFC国际会议论文集中。

颗粒流PFC2D (Particle Flow Code in 2 Dimensions)平台数值模拟单元有两种:颗粒圆筒和颗粒(disc or particle),主要用于平面应力和平面应变的特殊情况;颗粒流PFC3D(Particle Flow Code in 3 Dimensions)的数值模拟单元是三维球体颗粒(granular),主要用于三维受力分析。 Cundall(2002)博士认为PFC在描述岩土体介质特殊特性方面有着其他常用数值方法不可比拟的优势,主要表现在如下方面: (l)能自动模拟介质基本特性随应力环境的变化; (2)能实现岩土体对历史应力一应变记忆特性的模拟(屈服面变化 Kaiser效等); (3)反映剪胀及其对历史应力等的依赖性; (4)自动反映介质的连续非线行应力一应变关系屈服强度和此后的应 变软化或硬化过程; (5)能描述循环加载条件下的滞后效应; (6)描述中间应力增大时介质特性的脆性一塑性转化;

气液两相流

热物理量测试技术1 概述 两相流广泛应用于热能动力工程、核能工程、低温工程以及航天领域等许多领域。所谓两相流,广义上讲是指一种物质或两种物质在不同状态下的流动,其中气体和液体一起流动称为气液两相流。对于两相流中的气液混合物,它们可以是同一种物质,即汽—液(如水和水蒸气),也可以是两种不同的物质,即气—液(如水和空气混合物)。气液两相流是一个相当复杂的问题,。在单相流中,经过一段距离之后,就会建立一个稳定的速度场。但对于两相流,例如蒸汽和水,则很难建立一个稳定的流动,因为在管道流动中有压降产生,由于此压降作用会产生液体的蒸发,所以在研究气液两相流时必须考虑两相间的传热与传质问题。 两相流学科还处于半经验半理论阶段,对于两相流的流动和传热规律进行研究时,除了依靠各种数学物理模型外,还要依靠实验,这就需要两者相结合从而更好地进行研究。 2 两相流压降测量[1] 压降,即两相流通过系统时产生的压力变化,是两相流体流动过程中的一个重要参数。保持两相流体流动所需的动力以及动力系统的容量和功率就取决于压降的大小。一般说来,两相流体流动时产生的压降一般由三部分组成,即摩擦阻力压降、重位压降、加速压降,管道系统出现阀门、孔板等管件时,还需测量局部压降。目前,常用差压计或传感器来测量两相流压降。 2.1 利用差压计测量压降 应用差压计测量气液两相流压降的测量原理图如图1所示。所测压降为下部抽头的压力与上部抽头压力之差。在差压计的Z1截面上可列出压力平衡式如下: (2.1)式中,为取压管中的流体密度;为差压计的流体密度。 由(2.1)可得: (2.2)由上式可知,要算出压降的值,必须知道取压管中的流体密度和差压计读数。 当管中流体不流动时:

数值模拟收敛性问题

如何加快数模计算以及如何解决数模计算的收敛性问题? 收敛性问题可以说是所有从事数模工作的人员都会面临的问题,本文将以ECLIPSE软件为例从两方面介绍收敛性问题。第一方面介绍数值模拟计算与收敛有关的一些概念。第二部分介绍如何通过修改模型数据来加速计算,解决收敛性问题。 一、数模计算的收敛性:在了解收敛性之前,应该首先了解几个基本概念: 1、报告步:一个数模作业包括多个报告步,报告步是用户设置要求多长时间输出运行报告,比如可以每个月,每季度或每年输出运行报告,运行报告包括产量报告和动态场(重启)报告。在ECLIPSE软件中,报告步是通过DATES和TSTEP关键字来设置的。 2、时间步:一个报告步包括多个时间步,时间步是软件自动设置(VIP需要用户设置)即通过多个时间步的计算来达到下一个报告步,以ECLIPSE为例,假如报告步为一个月,在缺省条件下,ECLISPE第一个时间步取一天,然后以三倍增加,即第二个时间步取三天,然后取九天,下一个时间步是17天来达到30天的报告步,然后会以每30天的时间步来计算。时间步可以通过TUNING关键字来修改。 3、非线形迭代:一个时间步包括多次非线形迭代。在缺省情况下,ECLIPSE如果通过12次的非线形迭代没有收敛,ECLIPSE将对时间步减小10倍。比如下一个时间步应该是30天,如果通过12次的迭代计算不能达到收敛,ECLIPSE将把时间步缩短为3天。下一个时间步将以1.25倍增长,即3.75天,4.68天,。。。。如果在计算过程中经常发生时间步的截断,计算将很慢。 4、线形迭代:一个非线形迭代包括多次线形迭代。线形迭代是解矩阵。 在ECLIPSE输出报告PRT文件中可以找到时间步,迭代次数的信息, STEP 10 TIME= 100.00 DAYS ( +10.0 DAYS REPT 5 ITS) (1-FEB-2008) “STEP 10” : 说明这是第10个时间步。 “TIME= 100.00 DAYS”:说明现在模拟到第100天。 “+10.0 DAYS”:说明这个时间步是10天。 “REPT" :说明为什么选10天做为时间步,REPT是指由于到了下一个报告步。 “5 ITS": 说明此10天时间步需要5次非线形迭代。 ”(1-FEB-2008)“:现在的模拟时间。 模拟计算的时间取决于时间步的大小,如果模型没有发生时间步的截断而且能保持长的时间步,那表明该模型没有收敛性问题,反之如果经常发生时间步截断,那模型计算将很慢,收敛性差。时间步的大小主要取决于非线形迭代次数。如果模型只用一次非线形迭代计算就可以收敛,那表明模型很容易收敛,如果需要2到3次,模型较易收敛,如果需要4到9次,那模型不易收敛,大于10次的化模型可能有问题,如果大于12次,时间步将截

各种钢管重量计算公式

各种钢管(材)重量计算通用公式大全 (圆管)钢管的重量=0.25×π×(外径平方-内径平方)×L×钢铁比重其中:π= 3.14 ;L=钢管长度钢铁比重取7.8 所以,钢管的重量=0.25×3.14×(外径平方-内径平方)×L×7.8 * 如果尺寸单位取米(M),则计算的重量结果为公斤(Kg) 钢的密度为: 7.85g/cm3 (圆管)钢管每米的理论重量(钢的密度为7.85kg/dm3)计算公式: W=0.02466(D-S)S 式中:W--钢管每米理论重量,kg/m; D--钢管的公称外径,mm; S--钢管的公称壁厚,mm。 钢材理论重量计算 钢材理论重量计算的计量单位为公斤( kg )。其基本公式为: W(重量,kg )=F(断面积mm2)×L(长度,m)×ρ(密度,g/cm3)×1/1000 各种钢材理论重量计算公式如下: 圆钢盘条(kg/m) W= 0.006165 ×d×d d = 直径mm 直径100 mm 的圆钢,求每m 重量。每m 重量= 0.006165 ×1002=61.65kg 螺纹钢(kg/m) W= 0.00617 ×d×d d= 断面直径mm 断面直径为12 mm 的螺纹钢,求每m 重量。每m 重量=0.00617 ×12 2=0.89kg 方钢(kg/m) W= 0.00785 ×a ×a a= 边宽mm

边宽20 mm 的方钢,求每m 重量。每m 重量= 0.00785 ×202=3.14kg 扁钢(kg/m) W= 0.00785 ×b ×d b= 边宽mm d= 厚mm 边宽40 mm ,厚5mm 的扁钢,求每m 重量。每m 重量= 0.00785 ×40 ×5= 1.57kg 六角钢(kg/m) W= 0.006798 ×s×s s= 对边距离mm 对边距离50 mm 的六角钢,求每m 重量。每m 重量= 0.006798 ×502=17kg 八角钢(kg/m) W= 0.0065 ×s ×s s= 对边距离mm 对边距离80 mm 的八角钢,求每m 重量。每m 重量= 0.0065 ×802=41.62kg 等边角钢(kg/m) = 0.00785 ×[d (2b – d )+0.215 (R2 – 2r 2 )] b= 边宽 d= 边厚 R= 内弧半径 r= 端弧半径 求20 mm ×4mm 等边角钢的每m 重量。从冶金产品目录中查出4mm ×20 mm 等边角钢的R 为3.5 ,r 为1.2 ,则每m 重量= 0.00785 ×[4 ×(2 ×20 –4 )+0.215 ×(3.52 – 2 ×1.2 2 )]=1.15kg 不等边角钢(kg/m) W= 0.00785 ×[d (B+b – d )+0.215 (R2 – 2 r 2 )] B= 长边宽 b= 短边宽 d= 边厚 R= 内弧半径 r= 端弧半径 求30 mm ×20mm ×4mm 不等边角钢的每m 重量。从冶金产品目录中查出30 ×20 ×4 不等边角钢的R 为3.5 ,r 为1.2 ,则每m 重量= 0.00785 ×[4 ×(30+20 – 4 )+0.215 ×(3.52 – 2 ×1.2 2 )]=1.46kg

浅析气液两相流及其应用

浅析气液两相流及其应用 浅析气液两相流及其应用 摘要:气液两相流存在于石油、天然气、动力、化工、水利、航天、环境保护等工业中,其研究已成为国内外学者广泛关注前沿学科。本文概要性的描述了气液两相流的应用背景、流动型式,并介绍了气液两相流参数检测的手段和两相流计算的基本方法。 关键词:气液两相流流动型式参数检测计算方法 1.气液两相流的应用背景 近些年来,石油、天然气、动力、化工、水利、航天、环境保护等工业的迅速发展促进了气液两相流的研究和应用。在实际应用中可以将凝析天然气简化的看作气相为甲烷,液相为水的气液两相流[3]。为了在实现天然气井口对凝析天然气气、液两相流量的实时在线测量,需要对其进行相应研究。再如,火力发电厂中锅炉的汽水分离、蒸发管中的汽水混合物的流动都属于气液两相流问题[1]。 2.气液两相流的流动型式 气液两相流中气液两相的分界面多变,其流动结构受各相的物理特性、各相流量、压力、受热、管道布置等影响。在不同的流型下,两相流的流体力学特性不同,因此为了研究两相流的运动规律,必须研究其运动型式。 在水平管道中,气液两相流常见流动形态如图1所示。 图1 水平管道中气液两相流流型 水平管中,气泡流的特征为液相中带有散布的细小气泡,由于受到重力的影响,气泡多位于管子上部。随着泡状流中的气相流量的增加,气泡聚结成为气塞,气塞一般较长,且多沿管子上部流动。当气、液两相流速均较小,会受到重力分离效应产生分层流,而当分层流动中气相速度较大时,气液的交界面将产生扰动波形成波状流。若气相速度再增大,则气液分界面由于剧烈波动将有一部分与管道顶部接触,分隔气相成为气弹,从而形成弹状流,大气弹则将在管道上部高速运动。

冶金论文

选课课号:(2012-2013-1)-BG11191-320401-1课程类别:必修课 《冶金工程概论》课程考核 (课程论文) 题目:论我国冶金工业的发展现状及存在的问题 学生姓名:谢安静 学号:2011443343 授课教师:张明远 班级:酒店管理11-1 教师评语: 成绩: 重庆科技学院冶金与材料工程学院 2012年11月中国重庆

论我国冶金工业的发展现状及存在的问题 谢安静酒店管理11-1 2011443343 摘要:我国冶金行业在高速发展的同时,仍然同国际先进水平存在着很大差距。面对差距,我们需要认清现状,我们应充分调动和利用各种积极因素,着力解决当前冶金行业运行和发展中的突出问题,确保其平稳较快发展。 关键词:现状、问题、发展

1.冶金工业的发展现状 1.1钢铁生产工艺流程逐步优化 20世纪90年代以来,世界钢铁工业在激烈的国际市场竞争中,由20世纪80年代以前的以扩大规模、增加产量为主转向降低消耗、降低成本、提高质量、增加品种和保护环境。博士论文,高速钢轧辊。钢铁工业技术进步的主流是缩短生产流程,减少工序,提高质量,降低消耗,提高效率。技术进步中有两大主要趋向:一是寻找可以替代传统工艺的新工艺流程的研究开发;二是现有工艺和技术装备的完善化。两大技术进步趋向互相竞争、相互渗透,促使钢铁工业不断提高钢材质量、减少消耗、降低成本、减轻对环境的污染,进一步走向集约化。 1.2新兴技术的不断发展 传统的钢铁生产工艺流程是一种“冷态”下间歇式生产的工艺流程。日本在20世纪60年代建设的10多个大型钢铁厂都是采用这种工艺流程。20世纪80年代以后,世界钢铁业已逐步将上述传统的钢铁生产工艺流程改造成为现代化“热态”连续生产工艺流程。这种工艺流程具有高效、连续、紧凑、智能等特点。20世纪80年代末期,德国、法国、日本、意大利、美国等钢铁工业发达国家开发成功接近最终钢材产品形状的连铸、连轧技术,如带钢、型钢的连铸连轧等。由于该技术具有工艺流程紧凑、生产周期短、物料消耗少、生产效率高等一系列优点,在近十多年来得到了快速发展。自从1989年世界第一条薄板坯连铸连轧生产线在美国纽柯公司克劳福兹维尔厂投产以来, 经过10多年发展,到2002年底,世界上已有38个薄板坯连铸连轧生产厂共56条生产线,总生产能力已超过5 500万吨。我国现已有5个钢铁企业建成8条薄板坯连铸连轧生产线,到目前为止又有5个钢铁企业正在建设厚板坯连铸连轧生产线,不久的将来总生产能力将达2000万吨,预计届时将占全世界同类生产线能力的1/4以上。博士论文,高速钢轧辊。2001年我国连铸比达到89.71%,已经超过了2000年的世界平均水平。2003年达到了96.96%,目前,全国重点大中型企业中,连铸比达到99%以上的企业已达41家。 带钢连铸连轧技术是世界主要钢铁生产国家正在积极开发应用的一项重大钢铁生产前沿技术,它将是21世纪钢铁生产技术的一个主要发展方向。 1.3钢铁产量不断增长

实验三气液两相流实验

气液两相流实验 实验1 垂直上升管中气液两相流特性实验 一、实验目的 1. 在大型电站锅炉中垂直布置的锅炉水冷壁管被广泛应用,本实验将模拟其两相流现象和水动力特性; 2. 通过观察垂直上升管中气液两相流的流型,进一步加深了解垂直上升管中气液两相流型的特点; 3. 对垂直上升管中气液两相流的压力降有比较直观的认识,并掌握垂直上升管中气液两相流的压力降的计算方法。 二、实验仪器 仪器名称型号参数范围 水泵FS40 11m3/h 气泵HG-1100 180m3/h 电磁流量计DXLD-25 0.53-21m3/h 转子气体流量计LZB-4 0-400L/h 转子气体流量计LZB-15 0-4m3/h 转子气体流量计LZB-25 0-50m3/h 三、实验原理图 1 水箱 2 空气压缩机 3 磁力泵 4 转子流量计 5电磁流量计 6 气液混合器 7 减压阀 8 调节阀 9截止阀 10球阀 11 水集箱 12 针阀 13 过滤器 四、实验任务 1.观察垂直上升管中气液两相流的流型:

(1)打开系统电源,使气体、液体流量计预热2分钟;然后打开采集程序,记下采集程序上显示的气路 和水路温度(根据此温度查出水和空气的密度); (2)打开磁力泵,将主路的调节阀开度调小和旁路的调节阀开度调大,同时将垂直上升管实验段水路的球 阀开启,使水缓慢地流过实验段,直到取压管内大体上充满水为止; (3)关闭磁力泵和水路的球阀,打开空气压缩机和气路的球阀,将50-500L/min转子流量计一路的针阀开 启,调节针阀开度,使转子气体流量计所测得的体积流量保持在300L/min;打开磁力泵,调节主路和旁路的调节阀开度,将主路阀门开度达到最小,旁路阀门开度达到最大。按下表调节气量和水量,观察并记录垂直上升管中气液两相流的流型的变化; 水流量(L/min)0.6-1.6(96L/h,LZB-10)0.7-1.4 空气(L/min)160-220(13.2m3/h,LZB-25) 18-36 流型环状流块状流 水流量(L/min)0.65-1.5 3.6-5.6 空气(L/min)5-158.6-15.6 流型弹状流泡状流 (4)实验完毕时,先关闭磁力泵,然后关闭实验段水路的球阀,再关闭气路的球阀,最后关闭空气压缩机同时关闭气路的针阀。

我国国内连铸结晶器数值模拟近况

我国国内连铸结晶器数值模拟近况 摘要:本文就近年来我国在连铸铸坯在结晶器内热力生长行为方面的数值模拟研究做了一个简要的概述,着重叙述了几种常用的方法及各自的特点。 关键词:连铸,结晶器,数值模拟 1 前言 鉴于钢铁生产的特殊原因和现实检测手段的限制,我们只能通过相似原理的水模拟和数值模拟来研究现场的工艺过程。由于数学计算方法和相应大型计算软件的发展,数值模拟在冶金工业的应用得到了空前发展。数值模拟相对于其他方法具有用时短、节省人力、物力,尤其是可大幅度的改变各种参数的范围,对研究工艺过程,改进工艺方案,甚至是开发新工艺,有着得天独厚的条件。本文着重总结了近年来在我国人们对钢水在结晶器中的凝固行为过程的温度场和应力场方面的数值模拟。随着人们对结晶器中钢水凝固行为的研究,相关的数值模拟越来越全面,越来越细致,可谓是硕果累累。 结晶器是连铸机中的核心技术,被称为连铸机的心脏。结晶器内钢水的、传热和应力变形行为对连铸工艺以及铸坯质量有着重要影响,结晶器内初生坯壳的温度场、应力场及变形分布影响着连铸工艺的顺行和铸坯质量的提高。同时,研究结晶器内钢水的流场、温度场、应力场及变形对结晶器的设计和生产也有着重要的指导意义。 2 结晶器数值模拟的现状 王刚恩等[1]通过有限元方法模拟研究了角部气隙的产生对钢水凝固行为和出生坯壳力学行 为的影响。王哲等[2]建立的有限元数学模型,采用热通量系数法反映实际的坯壳角部凝固特征现象,并在方坯验证的基础上对H连铸坯的凝固行为进行了数值模拟。康丽等[3]建立了三维稳态模型,模拟出了整个结晶器内的热和力学状态,研究了坯壳产生裂纹的可能性和结晶器锥度对气隙产生的影响。张炯明等[4]建立数学模型来研究结晶器的合理锥度,建立了热力直接耦合模型对气隙的生成、大小做了动态分析。荆德君等[5]建立了完全热力耦合的二维热-热弹塑性有限元模型,为优化结晶器锥度,高拉速曲面结晶器提供理论依据和技术基础。张家泉等[6]模拟了板坯结晶内钢水的凝固温度分布和结晶器变形。陈志平等[7]利用MSC.Marc 接触理论建立了连铸结晶器内热和应力状态的有限元耦合分析模型。崔立新等[8]再考虑结晶器水槽结构尺寸和分布的基础上建立了热和应力的耦合模型。兰岳光等[9]通过调整气隙厚度模拟结晶器角部锥度变化, 使用有限单元法计算出不同锥度下方坯坯壳内的温度及应力场, 明确了结晶器角部锥度对方坯温度及应力分布的影响。谢严敏等[10]模拟了不同结晶器铜板厚度的应力,结晶器壁越薄,变形越严重,应力越大。得出了合适的结晶器厚度为10mm。张立夫等[11]建立二维完全热力耦合的有限元分析模型,研究划分沟槽对坯壳形成过程的影响。 3 几种计算模型的分析 由于连铸结晶器的传热特点,在传热过程当中铸坯的横向传热远远大于拉坯方向的传热,故而在研究结晶器坯壳凝固行为的数值模型常用二维切片法进行研究。 (1) 计算模型只计算结晶器区域内的钢水和铸坯坯壳,不包括结晶器器壁。对于铸坯坯壳与保护渣和结晶器器壁的传热边界条件,采用平均热流量法,常用的平均热流量有Savage 提出的见公式(1):

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位移云图

气液两相流传热实验

气液两相流传热实验 一、实验目的 1、通过测定换热器冷、热流体的流量,测定换热器的进、出口温度,熟悉换热器性能的测试方法; 2、了解套管换热器的结构特点及性能。 3、通过测定参数计算换热器流体的热量;计算换热器的传热系数;并整理成准数关联式形式。 二、基本原理 1、概述 本换热器性能测试实验装置,主要对应用较广的套管式换热器进行其性能的测试。其中,对套管式换热器可以进行顺流和逆流两种方式的性能测试。 换热器性能实验的内容主要为测定换热器的总传热系数,对数传热温差和热平衡误差等,并对实验数据进行整理,分析流体无相变时的对流传热系数与Dittus-Boelter 关联式。 2、实验装置参数 本实验所用的热水加热采用电加热方式,采用热水加热常温空气。冷—热流体的进出口温度采用pt100加智能多路液晶巡检仪表进行测量显示,实验台参数如下: (1)电加热管总功率:3KW (2)冷热流体风机:允许工作温度:<80℃,额定流量:76 m 3/h 电机电压:220V 电机功率:750W (3)孔板流量计: 流量:8-30m 3/h 允许工作温度:0-80℃ 3、对流传热系数α的测定: 根据传热总方程,用实验测定。 m Q S t α= ? 式中:α-管内流体对流传热系数,W/(m 2·℃); Q -传热速率W ;

S -管内换热面积, m 2 ; ?t m -对数平均温度差,℃。 本实验中,具体的计算过程如下: ,,56()m h p h Q q c t t =-,热水的物性数据取定性温度 56 2 t t +下的数值,计算质量流量, /m c V t q q kg s ρ=。 换热面积2 o S d l m π=,此处管内径0.016m ,壁厚0.0015m ,管长1.3m 。 {}()2121/ln /T T T T t m ???-?=? 851t T T -=? 762t T T -=? t 5,t 6为热流体进出口温度, T 7,T 8为冷流体进出口温度。 依次可以求出对流传热系数α。 4.关联式的整理 求出努塞尔准数i u d N αλ=,普朗特准数p r c P μλ=,雷诺准数i e ud R ρμ =。准数关联式为0.3m u e r N AR P =,其中的A ,m 可以用图解法得到。具体做法是:先把关联式变换成 0.3m u e r N AR P =,然后两边取对数0.3 ln ln ln u e r N m R A P =+,求出m 和A 。 三、实验装置与流程 1.实验装置流程: 本实验装置采用空气和用阀门换向进行并逆流实验;工作流程如图2-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 二维边坡稳定分析模块

颗粒流(PFC)简介

颗粒流(PFC)简介 2009-10-07 11:14:48| 分类:岩土工程| 标签:|字号大中小订阅 注:今天偶然间见到颗粒流的概念,以前一直不了解,今天查了查,贴在这里,以备以后可以温故知新。 本文内容源自浙江大学罗永先生的博士论文,使得吾辈能花较少的时间看到广博的知识,在此特向其表示感谢! 岩土工程数值计算总体上可以分为两大类:一类是基于连续介质 力学理论的方法,如有限元法(FEM)和快速拉格朗日法(FLAC(1tasea,2002))等;另一类是不连续介质力学的方法,如离散元法UDEC(1tasca,2000)、3DEC(Itasea,1998)、PFC(Itasea,2002)和块体理论DDA(石 根华,1988)等。离散元方法按其用途又可以分为宏观离散元方法和 细观离散元方法,前者主要针对解决规模相对较大的不连续面,如断层节理结构与基础之间的结合面等引起的问题(UDEC,3DEC),后者则着重于数目众多具有不连续特性的接触面或点,如破碎岩体中的破裂面、砂土中的接触面(点)和材料中颗粒之间的接触面(点)等。 PFC(Particle Flow Code)是在著名学者Peter Cundall主持下采用细观 离散元理论(又称为粒子流理论)开发的一种数值计算平台,可以广泛地应用于研究细观结构控制问题。目前,PFC在世界上的应用并不广泛,成果报道也主要集中在PFC国际会议论文集中。

颗粒流PFCZD (Particle Flow Code in 2 Dimensions)平台数值模拟单元有两种:颗粒圆筒和颗粒(disc or particle),主要用于平面应力和平面应变的特殊情况;颗粒流PFC3D(Particle Flow Code in 3 Dimensions)的数值模拟单元是三维球体颗粒(granular),主要用于三维受力分析。 Cundall(2002)博士认为PFC在描述岩土体介质特殊特性方面有着其他常用数值方法不可比拟的优势,主要表现在如下方面: (l)能自动模拟介质基本特性随应力环境的变化; (2)能实现岩土体对历史应力一应变记忆特性的模拟(屈服面变化Kaiser效等); (3)反映剪胀及其对历史应力等的依赖性; (4)自动反映介质的连续非线行应力一应变关系屈服强度和此后的应变软化或硬化过程; (5)能描述循环加载条件下的滞后效应; (6)描述中间应力增大时介质特性的脆性一塑性转化; (7)能考虑增量刚度对中间应力和应力历史的依赖性; (8)能反映应力一应变路径引起的刚度和强度的各向异性问题; (9)描述了强度包线的非线性特征;

气液两相流流型实验报告

气液两相流流型实验报告 实验名称:气液两相流流型 实验目的: 1. 熟悉台架,掌握流量测量仪表的使用; 2. 掌握常见两相流流型的划分方法及相关规律,观察水平管中不同流型的特点; 3. 根据各工况点实验数据绘制两相流流型图,并与典型流型图做比较。 实验任务: 实验测量数据: ,,,. (1) 测取不同情况下气相,液相流量;记录P P t t w 气减室 (2) 判别流型 要求: (1) 实验数据汇总表; (2) 绘制αβ -曲线 (3) 根据实验数据用Weisman图判别流型 实验原理 1、水平管道中气液两相流流型的划分及各流型特征 在水平管道中的气液两相流,由于重力影响使流型结构呈现不对称性,因而水平管中的流型特征变得较为复杂。Oshinowo流型划分原理使流型变得相对简单,根据Oshinowo的划分原则,一般把水平管道中的流型划分为六种,泡状流、塞状流、层状流、波状流、弹状流、环状流。 (1)泡状流 在泡状流中,气相是以分离的气泡散布在连续的液相内,气泡趋向于沿管道上半部流动,这种流型在含气率低时出现。 (2)塞状流 在塞状流中,小气泡结合大气泡,如栓塞状,分布在连续的液相内,大气泡也是趋向于沿管道上部流动,并且在大气泡之间还存在一些小气泡。 (3)层状流 在层状流中,两个相的波动被一层较光滑的分界面隔开,由于重力和密度不同,气相在上部液相在下部分开流动。层状流只有在气相和液相的速度都很低时才出现。 (4)波状流

当气流速度增大时,在气、液分界面上掀起了扰动的波浪,分界面由于受到沿流动方向的波浪作用而变得波动不止。 (5)弹状流 当气体流速更高时,分界面处的波浪被激起与管道上部管壁接触,并形成以高速沿管道向前推进的弹状块。 (6)环状流 当气体流速进一步增高时,就形成气核和环绕管周的一层液膜,液膜不一定连续均匀的环绕整个管周,管子的下部液膜较厚,在气芯中也夹带有液滴。 表1水平绝热管中的流型变化 A表示环状流(annular);B表示气泡(bubble); BTS表示中空气弹(blow through slug);D表示液滴(droplet); F表示液膜(film);IW表示平缓波(inertial wave); LRW表示大翻卷波(large roll wave);PB表示气栓加气泡(plug&bubble);PF 表示气栓加泡沫(plug&froth);R表示涟漪波(ripple); RW表示翻卷波(roll wave);S表示气弹(slug);ST表示层状流(stratified)。

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·

气液两相流

---------------------------------------------------------------最新资料推荐------------------------------------------------------ 气液两相流 气液两相流流型识别理论的研究进展摘要:介绍了气液两相流的识别理论,探讨了气液两相流流型的划分方法。 叙述了两相流流型软测量方法,并重点介绍了图像处理识别、在线流型技术识别、神经网络、基于压差波动理论、混沌理论等识别流型的新方法。 关键词:气液两相流;流型识别0 引言相的概念通常是指某一系统中具有相同成分及相同物理、化学性质的均匀物质成分,各相之间有明显可分的界面。 从宏观的角度出发,可以把自然界的物质分为三种,即:气相、液相和固相。 单相物质的流动称为单相流,如气体流或液体流。 所谓两相流(Two-Phase Flow)或多相流(Multiphase Flow)是指同时存在两种或多种不同相的物质的流动。 近年来随着国内外石油和天然气工业的发展,迫切需要开发出精度较高的油气水三相流量在线测量仪,以便掌握各个油井的生产动态。 然而,多年来尽管在这方面进行了大量的研究工作,取得了一些进展,但是仍然没有彻底清晰地认识和了解油气水三相混合物的流动型态。 在现今的多相流检测技术领域中,流型的识别问题变得越来越重 1/ 10

要。 1 两相流流型由于存在一个形状和分布在时间和空间里是随机可变的相界面,而相间实际上又存在一个不可忽略的相对速度,致使流经管道的分相流量比和分相所占的管截面比并不相等。 这就导致了两相流动结构多种多样,流型十分复杂。 流型是影响两相流压力损失和传热特性的重要因素。 两相流各种参数的准确测量也往往依赖于对流型的了解。 因此为了对两相流的特征参数进行测量,必须了解它们的流型。 1.1 垂直上升管中气液两相流流型(1)、泡状流(Bubbly Flow):气泡以不同尺寸的小气泡形式随机离散分布在流动的液体中。 显然,此时气体为离散相,而液体为连续相。 随着气速的增加,气泡尺寸会不断增大。 (2)、段塞流(Slug Flow):在气泡流动中当气泡的浓度增高时,气泡聚合为直径接近于管内径的塞状或炮弹状气泡,气泡前端部分呈现为抛物线形状。 在这些塞状气泡之间可带有小气泡的液团。 当气泡快速上升时,液体在气泡与管内壁间的间隙中流动。 (3)、混状流(Churn Flow):当气泡速度进一步增大时,段塞流中的气泡速度也随之增加并产生破裂、碰撞、聚合和变形,与液体混合成为一种不稳定的上下翻滚的湍动混合物。 此时气液两相界为离散相。

相关文档
最新文档