一维双流体模型的压力修正算法

针对一维双流体模型,通过推导双流体模型的压力修正方程,及相含率修正方程,将原本应用于单相不可压流动的压力修正系列算法推广至双流体模型求解,提出双流体模型的压力修正算法.在离散过程中运用高阶有界格式,在保证二阶以上精度的基础上克服了由相含率分布的阶跃所造成的数值结果非物理震荡.与公开发表算例进行对比,验证求解的可靠性.     

T-T吸收剂脱硫模型的改进及与双流体模型的耦合

原有T-T吸收剂固硫模型以时间作为反应状态参数,无法用于欧拉法脱硫模拟.为此,基于钙转化率的概念,对已有的T-T吸收剂脱硫模型进行了改进.改进的模型形式简单,反应速率与钙转化率的平方呈双线性分布,模型中快慢反应的转折点和反应终止点均随温度升高而增大,表征了反应活性的提高.进一步推导了吸收剂钙转化率的输运方程并对两相质量守恒方程和SO2组分方程进行了修改,将脱硫模型与两相流模型进行了耦合,为循环流化床脱硫数值模拟打下了基础.     

双流体等离子体模型的动力学可容变分

动力学可容变分方法是一种广义哈密顿系统中的李扰动变换方法,能自动保证卡西米尔函数在相应阶数上的守恒性质.通过动力学可容方法得到了双流体在欧拉描述中的一组约束变分,而后利用这组变分对双流体哈密顿量取极值得到了平衡方程.     

三维煤粉燃烧全双流体模型的数值模拟

1引言目前煤粉燃烧综合模型最为广泛的是以轨道模型为基础，美国杨伯翰大学的先进燃烧工程中心（CERC）自1980年起研制二维煤粉燃烧程序PCGC－2，从1990年起研制三维煤粉燃烧的PCGC－3程序[1]。与轨道模型发展的同时，美国Rabcock＆Wilcox公司Fiveland[2]等人研制了FURMO程序，用无滑移模型对560MW侧墙喷燃煤粉炉进行了三维全模拟。其特点是首次用全欧拉的处理方法计算三维煤粉燃烧过程，其不足之处是不考虑气粒两相间的速度滑移和温度滑移．总的看来，用轨道模型模拟煤粉燃烧，易于考虑颗粒反应经历，也可给出两相之间的速度及温…     

炉内两相流动和煤粉燃烧的双流体-轨道模型

本文首次用双流体一轨道（连续介质－轨道ＣＴ）模型对大尺寸四角喷燃炉内气粒两相流动及煤粉燃烧进行了模拟。该模型基于欧拉气相方程组、欧拉颗粒连续方程和动量方程以及拉氏颗粒能量和质量变化的方程,对各子模型用ｋ－ε－ｋｐ两相湍流模型,ＥＢＵ－Ａｒｒｈｅｎｉｕｓ湍流燃烧模型,煤粉颗粒的水分蒸发,热解挥发和焦炭燃烧的扩散－动力模型,ＤＯ（离散坐标）辐射模型。采用了将坐标扭转一定角度的方法减小入口射流和网格斜交造成的伪扩散。编制了ＬＥＡＧＡＰ－ＦＵＲＮＡＣＥ－３程序,分别对冷态模型炉内两相流动和大尺寸炉内三维两相流动和煤粉燃烧进行了模拟,并与颗粒轨道（ＳＴ）模型的模拟结果进行了对照。采用ＰＤＰＡ对冷态模型炉内气粒两相流场进行了测量。冷态两相流动的模拟与实验结果的对比表明ＣＴ模型的模拟结果和实验符合较好,ＳＴ模型所得颗粒浓度分布和实验山入较大。热态模拟的结果给出了两相速度,气相温度,组分浓度和壁面热流分布。模拟结果定性合理。模拟结果显示在出口处由于气流旋转,有一局部高温区存在。     

分体空调室内机的双流体数值模拟

本文采用Navier-Stokes方程和标准k-ε两方程,对分体空调室内机(包括换热器和横流风机)整体系统,以水蒸汽和干空气两种气体为介质,进行了等温流动和实际工况流动的双流体三维数值模拟和研究.两种气体在叶轮的进口面上不均匀分布,叶轮轴向的压强和速度分布存在强度和位置的差异.采用密度干扰方法对水蒸汽的凝结参数近似处理,给出了水蒸汽体积分数和湍流强度在换热器域的分布、混合物温度分布等内流特性.与实验结果相比,双流体模拟的工作点外特性匹配较好.同时,对蜗舌间隙参数变化对性能的影响做了说明.     

一维随机成核生长模型的边界效应

通过计算机模拟，研究周期边界与固定边界对一维随机成生长模型的聚集生长图形形状及相应的分形维数的影响。     

快速成型中有效保留模型特征的自适应分层方法

目前,虽然各种自适应分层方法层出不穷,但针对模型特征偏移和丢失的自适应分层方法还相对较少。针对此点,提出一种以体积误差为基础的自适应分层方法。首先,在成型方向上将单个切片层对三角面片的切割分为未跨越、跨越到相邻、跨越到不相邻三角面片3种情况,并建立了数学模型。分析了模型特征的保留特性,并结合3种情况依照体积误差公式给出了具体的算法流程。最后,为验证方法的有效性,对一零件应用3种不同方法进行分层,得到相应的模拟成型件,并与CAD原件进行对比。结果表明,体积自适应算法较其余2种方法与CAD模型轮廓最为接近,在采样点处的偏离距离分别为0.09 mm和0.10 mm,形状精度最高。     

中温循环流化床烟气脱硫基础-曳力模型与流动的数值研究

为了准确模拟中温烟气脱硫的循环流化床(CFB)反应器内的脱硫过程,本文首先采用欧拉双流体方法研究了床内的稠密气固两相流动,讨论了表征气固两相相互作用的关键参数-曳力模型,并对适于均匀流动的经典曳力模型进行了颗粒团聚所致的非均匀效应的修正.计算的床内存料量和总压降与实测值之间的误差均小于5%.本研究为中温脱硫过程的优化和降低能耗奠定了流体力学基础.     

一维Al/SrF2低温超导体-电介质光子晶体能带结构

基于二流体电子模型和平面波展开法,计算了一维Al/SrF2超导体-电介质光子晶体的能带结构.结果表明:随着超导层厚度的增加,第一光子带隙中心频率和截止频率均发生蓝移,且第一带隙宽度逐渐增加到一个峰值后又逐渐变窄.更重要的是,在低于临界温度的超低温环境中,温度的微小变化,对该类光子晶体的带隙宽度、中心频率以及截止频率均有...     

双流体颗粒-壁面碰撞模型用于旋流流动

对旋流两相流动,用考虑壁面粗糙度的双流体颗粒-壁面碰撞模型,结合二阶矩两相湍流模型和不同颗粒相边界条件,进行了数值模拟。结果表明,由于考虑了各方向雷诺应力之间的相互转化,雷诺应力从平均运动中得到能量,以及壁面对颗粒运动的衰减作用等因素,该颗粒壁面碰撞模型给出的模拟结果与实验符合较好,而广泛应用的零梯度颗粒壁面边界条件给出的模拟结果最差。     

用直接模拟法分析量子阱半导体激光器瞬态响应

根据光增益与载流子密度的对数关系,通过适应于多量子阱激光器的速率方程的直接模拟分析,得到了注入电流、阱数和腔长对多量子阱激光器的激射阈值、开关延迟时间、弛豫振荡频率和光输出等参量之间的依赖关系.运用相图确立了在瞬态过程中,载流子数密度和光子数密度之间的转化过程.从而为改善量子阱激光器的高频调制特性以及优化设计器件结构参数提供了理论依据.     

用直接模拟法分析量子阱半导体激光器瞬态响应

根据光增益与载流子密度的对数关系 ,通过适应于多量子阱激光器的速率方程的直接模拟分析 ,得到了注入电流、阱数和腔长对多量子阱激光器的激射阈值、开关延迟时间、弛豫振荡频率和光输出等参量之间的依赖关系 .运用相图确立了在瞬态过程中 ,载流子数密度和光子数密度之间的转化过程 .从而为改善量子阱激光器的高频调制特性以及优化设计器件结构参数提供了理论依据     

空化紊流流动的数值计算模型及其验证

本文介绍了一种空化紊流流动的数值计算模型,并计算了射流放水阀内部的空化流动现象。该模型基于均相平衡模型和液相与汽相传输方程,基本方程采用N-S方程,空化模型采用Kunz等提出的汽液质量转换模型。紊流封闭采用标准的κ-ε素流模型;稳定流动计算采用扩展的SIMPLE压力修正方法,非稳定计算采用PISO算法。为评价该数值模型,计算了绕射流放水阀的空化流动,并与实验结果进行了对比,计算结果与实验结果取得较好的一致,说明该流动模型和计算方法是可靠的。     

水平管段塞流持液率波动规律研究

以空气和水为工质,采用双平行电导探针,在水平多相流环道(d=50.00mm)中研究了段塞流持液率与压力波动的关系以及折算气速、折算液速变化对持液率的影响,并将持液率的试验值与预测模型进行了比较.结果表明:持液率波动能更真实地反映段塞流动特性,可以用来确定液塞频率;增加折算气速时,液塞持液率减小;增加折算液速时,持液率增大.     

两介质段瞬变流的修正特征线法

本文针对变波速可压缩两介质段非定常流的数值模拟,提出了修正特征线法。为使计算网格保持矩形,将一个大于1的修正因子α引入动量方程,从而在一定的范围内对特征线斜率进行调整,并据此确定时间和空间步长。其中因不满足稳定性条件而造成的误差由修正因子对动量方程进行修正。这种方法适用于一般瞬变流和振荡流。为解决介质界面连续移动的问题,本文提出了可自动调节空间步长的浮动网格,这使计算单元介质得以保持单一。通过在压力项中加入人工粘性改善了计算的稳定性。文章最后给出了计算实例。     

着火油罐燃烧过程预测的通用模型

本文根据油罐火灾的燃烧特点,从基本的传热规律和燃烧机理出发,以油品表面热反馈的能量平衡为前提,并利用化工热力学方法计算了火焰的平均温度,通过迭代方法,并结合相关的实验数据和结果,提出了油罐火灾燃烧过程的通用模型,该模型可适用于多种碳氢类燃料。本文以汽油为例,通过计算可以得到着火油罐的多项燃烧特性,包括燃烧速度、平均火焰温度等,是获取油罐周围热辐射规律的基础和前提,计算结果可为火灾现场战术的决策以及预案的制订提供依据。     

气液段塞流气弹区相界面结构理论研究

根据气液段塞流气弹区相界面结构特征将气弹分为气弹头、气弹体、水跃面和气弹尾四部分,并根据各自的流动和界面结构特征分别进行模化,建立了描述不同倾角的圆管内气液段塞流气弹区相界面结构的一维理论模型.该模型的计算结果表明气液混合Froude数、管道倾角和气弹长度显著影响气液段塞流气弹区相界面结构,计算与实验结果吻合良好.     

光电经纬仪动力学模型分析

从动量定理和动量矩定理两方面分析光电经纬仪动力学特性，建立光电经纬仪的非线性动力学模型。结合某型号光电经纬仪的具体参数对模型进行了分析。