液—液旋流分离管中强旋湍流的Reynold应力输运方程数值模拟

应用湍流Ｒｅｙｎｏｌｄ应力输运方程模型（ＤＳＭ）对液－液旋流分离管中的强旋湍流进行了数值模拟,并与ＬＤＶ测量结果进行了比较,结果表明：ＤＳＭ模型不仅较合理地预报出 切向速度的Ｒａｎｋｉｎｅ涡结构及其位置,而且提示了液－液旋流管中切向速度所特有的双峰分布现象;对轴向速度预报给出了近壁下行流,近轴上行流以及介于两者之间的零速区等;湍流动能在肇流管中的继上游近壁在、中部近轴大和下游近似均布,旋流管中的静     

超燃冲压发动机燃烧室内湍流燃烧流场的数值模拟研究

应用经显式可压缩修正的SST湍流模型,对超燃冲压发动机燃烧室流场进行了三维数值模拟,并与实验结果进行了对比.通过对燃烧室流场数值模拟结果的分析给出了湍流燃烧流场的特征.为了定量分析超声速流动条件下湍流燃烧的作用特点,基于湍流燃烧理论与湍流燃烧的数值模拟结果,通过对控制无量纲参数在流场中变化的研究发现：超声速燃烧发生的区域为充分发展的湍流区;在喷口附近湍流和燃烧的作用最强烈;超声速燃烧流动中,湍流引起的火焰非定常效应和局部火焰熄灭现象均可忽略,定常火焰面的近似假设是成立的.     

阶梯溢流坝水流数值模拟及特性分析

针对Mixture模型,考察了Realizable k-ε模型,SSTk-ω模型,v2-f模型,LES模型等4种湍流模型对阶梯溢流坝水流数值模拟的适用性.通过对非掺气区域坝面平均流速、旋度大小、边界层厚度的计算,并与实验结果进行分析比较,结果表明：Realizable k-ε模型考虑了流体微团的旋转效应,提高了对较大曲率流动、旋流流动和涡旋运动的计算精度,对阶梯溢流坝有很好的适用性.以Realizable k-ε模型计算结果为基础,进一步分析了边界层内的速度分布、阶梯面压力分布等坝面水流特性.     

百叶窗式翅片散热器的数值模拟和实验研究

先对百叶窗式翅片散热器进行了理论计算;然后对其物理模型进行合理的简化处理,利用ANSYS软件,采用SIMPLE算法和标准k-ε湍流模型,通过求解三维N-S方程和能量方程模拟了散热器的散热过程;最后进行了试验研究。试验测试和理论计算、数值模拟三者数据基本吻合,说明数值模拟方法能真实反映百叶窗式翅片散热器的传热特性,为散热器的设计和改进提供理论依据。     

熔液法晶体生长流场的数值模拟

本文给出了对二种熔液法生长晶体的液相流场进行的计算机模拟的结果,对ＣＺ法晶体生长的一常流场数值计算,采用了二阶迎风ＱＵＩＣＫ差分格式的有限控制体积法,通过对１２种不同参数的流场情况所作的计算表明,二阶ＱＵＩＣＫ格式具有更高的计算精度和数值稳定性。     

适用于超声速湍流扩散燃烧流动的火焰面模型

为了建立适用于超声速流动的湍流扩散燃烧模型,本文首先分析了火焰面模型应用于超声速流动的物理基础,然后数值模拟了轴对称超声速射流形成的氢气/空气扩散燃烧流场,利用实验数据校正了火焰面模型中重要物理量标量耗散率的模型系数.计算结果与实验数据的对比表明,本文修正后的火焰面模型对超声速湍流扩散燃烧流动的模拟能力是令人满意的.基于火焰面模型理论以及数值模拟结果,本文首次研究了湍流脉动对平均状态方程以及化学反应源项的影响机理,得到以下结论：组分浓度和温度脉动相关项对平均状态方程影响很小;温度脉动会降低水的生成速率,但其影响较小;组分浓度脉动在接近于氧化剂一侧区域增加水的生成速率,而在另外的大部分区域会降低水的生成速率;组分浓度与温度的脉动相关作用会很大程度上降低水的生成速率.     

基于欧拉模型的淹没射流冲刷数值模拟

采用Eulerian双流体模型对二维淹没射流作用下的无粘性沙床冲刷进行数值模拟,水沙两相间以及沙粒之间的相互作用通过动量交换和相间作用力体现.分别针对Aderibigbe等人的轴对称实验和本文设计的平面二维实验进行计算,动态平衡状态下冲坑剖面形态的计算结果与实验结果吻合较好.分析结果表明：不同冲刷强度下射流水体速度场具有不同的特征;沙粒受水体作用力、重力和沙粒间的摩擦力的综合作用,在沙粒的密度、直径、孔隙度等参数均相同的条件下,水体作用力决定沙粒的运动方式,从而影响冲坑的最终形态;孔隙度是影响冲坑形态的重要因素,在本文的研究范围内,随着孔隙度的增加,孔隙水的流速增大,冲坑深度增加.     

矩形液池中热毛细对流起振过程的数值模拟

本文数值研究微重力条件下浅液池中大Pr数流体(Pr=105.6)热毛细对流的转捩过程.当液池两侧的温度差超过临界温度差时,在液池中会出现热毛细振荡对流.振荡对流的涨落值远小于时间平均值,涨落值表现为从冷端向热端传播的热流体波.起振时,涨落值的胞元分布在液层的中心区域;随着外加温度差的增加,涨落值的胞元逐渐充满液层.本文的结果显示了有限延伸矩形容器液层中热毛细对流及其起振过程与其它类型热毛细对流不同的特征.     

颗粒惯性对气固两相各向同性湍流中温度统计特性的影响

颗粒运动轨迹上流体的温度统计特性对于理解非等温/反应气粒两相湍流的机理, 特别是对于检验非等温气粒两相湍流Lagrangian模型是十分重要的. 对带有平均标量梯度的气固两相各向同性湍流中颗粒及颗粒所见流体温度的统计行为进行了直接数值模拟研究, 讨论了颗粒惯性对于颗粒温度以及颗粒所见流体温度的Lagrangian统计特性的影响. 结果显示, 对于τ_p/τ_k<1的颗粒, 颗粒所见流体温度的脉动强度随τ_p/τ_k的增大而减小; 而对于τ_p/τ_k>1的颗粒, 其趋势相反. 小颗粒(τ_p/τ_k<5)温度的Lagrangian自相关系数R_p^T也随颗粒惯性(τ_p/τ_k)的增大而减小, 对于大颗粒这一趋势也相反. 颗粒运动轨迹上流体温度的自相关系数 都随颗粒惯性的增加而减弱, 而且随颗粒惯性的增加, 颗粒运动轨迹上流体温度的自相关比颗粒温度的自关联下降得快. 平均温度梯度的存在使得在沿平均温度梯度的方向上颗粒速度和温度有很强的关联性. 当τ_p/τ_k<1时, 其关联系数随颗粒惯性的增加而增大; 当τ_p/τ_k>1时, 这一系数的值与颗粒惯性无关.     

玻壳压制成型中残余应力的数学建模与模拟方法

玻壳成型过程中产生的残余应力对制品质量的影响很大．在分析玻壳成型特点的基础上，建立了玻璃压制成型过程中残余应力计算的数学模型，其中材料假定为热流变简单黏弹性材料，忽略了成型中的流动应力，讨论了模型问题的平衡推论及相容性方程，并详细分析了成型中不同阶段的边界条件．模型的数值求解采用了薄层理论，并通过空间上的分层离散和时间上的有限差分来进行．所提出的模型及求解方法可以很容易地扩展到通常的玻璃压制成型过程，也可用于分析许多与玻璃压制成型相关的问题，具有较广泛的参考价值．     

玻壳压制成型中残余应力的数学建模与模拟方法

玻壳成型过程中产生的残余应力对制品质量的影响很大。在分析玻壳成型特点的基础上，建立了玻璃压制成型过程中残余应力计算的数学模型。其中材料假定为热流变简单粘弹性材料，忽略了成型中的流动应力，讨论了模型问题的平衡推论及兼容性方程，并详细分析了成型中不同阶段的边界条件。模型的数值求解采用了薄层理论，并通过空间方向上的有限元和时间上的有限差分来进行。所提出的模型及求解方法可以很容易地扩展到通常的玻璃压制成型过程，也可用于分析许多与玻璃压制成型相关的问题，具有较广泛的参考价值。     

导波在全长粘结锚杆中传播的有限差分数值模拟

基于导波理论,把锚固锚杆当波导质来考虑,利用有限差分数值模拟软件FLAC3D4.0的动力学模块建立了三维锚杆锚固模型,模拟了导波在全长粘结锚杆中的传播特征,同时结合模型试验验证了数值模拟的正确性.在前人的研究基础上,深入和全面地研究了导波在锚杆锚固结构中传播的波速、波系和衰减特征,发现了一些新的规律和现,并对前人没有发现的一些现了探讨.研究结果表明：导波在锚杆锚固结构中传播的时,当体波衰减掉,还存在面波—stonely波,而stonely波在锚杆轴是不衰减的.研究结果可为锚杆检测及锚杆检测最优激发波的选取提供一定的鉴意义.     

三维可压大气中重力波波包非线性传播的数值模拟

采用全隐Euler格式(FICE)对重力波波包在风场中的三维非线性传播进行了数值模拟, 给出了重力波波包三维非线性传播的全过程, 分析了重力波的传播特性及风场对重力波传播的影响. 结果表明: 二维全隐Euler格式能够推广到三维情况下使用, 并能保持计算的长期稳定; 波振幅的增长比线性条件下的指数增长要快, 非线性传播使水平扰动速度大幅度增长, 会导致局地背景风场的增强; 波包非线性传播的路径、能量传输速度不同于线性重力波理论给出的结果, 非线性效应导致了重力波传播特性的改变；风场会改变重力波传播的路径和速度.     

三维Hall MHD数值模拟中初始载流子对场向电流分布的影响

用Hall MHD数值模拟的方法研究了Hall等离子体中不同初始粒子载流情形下磁场拓扑形态结构的改变以及场向电流与Alfven波的产生．在考虑了初始离子载流子的影响后,模拟结果中磁场的拓扑形态结构更加复杂．模拟结果中除了传统的By四极结构以外,还出现了一个与传统By四极结构相反的反四极结构,这种结构的出现使Hall MHD理论能解释完全电子载流情形下不能解释的观测现象,作为事例给出了Cluster卫星观测事例．同时还得出以下几个非常有意义的结果：1）受Hall效应影响的区域（空间变化尺度小于或相当离子回旋半径的区域）电子与离子分离．在非Hall效应影响的大部分区域,受初始离子＋y方向运动的影响整个磁结构向＋y方向偏移;而在受Hall影响的较小区域,受电子运动影响磁力线向-y方向弯曲．随之,By产生;2）由于By的出现,场向电流（FACs）产生．与完全电子载流的情形相比,结果中场向电流分布的中心随离子载流比例的增加向＋y方向偏移,场向电流主要分布在y〉0的区域;3）模拟结果中Ae≈0．76,Ai≈1．36,Ae×Ai≈1．03,完全符合Hall等离子体中的瓦伦关系,证实了Alfven波的存在．