孔板空调风口送风射流的数值模拟

介绍Ｎ点风口模型用于数值模拟室内空气流动时描述孔板类送风口的入流边界条件．然后采用该风口模型对不同的孔板风口出流条件算例进行数值计算，并就轴心速度衰减、射流扩展角以及断面流速分布等射流特性与实验数据进行了对比．比较结果表明，Ｎ点风口模型用于描述数值模拟室内空气流动的孔板类风口入流边界条件，可以获得工程上足够满意的结果．     

方型散流器空调室内空气流动的数值模拟

用N点风口动量模型和一个新零方程湍流模型对某办公室方型散流器空调的室内温度场和速度场进行了模拟,并和实验数据进行对比.结果表明,计算所得速度和温度分布与实测值吻合得很好,所用的风口模型和湍流模型能快速地将方型散流器空调通风的温度和速度场合理地模拟出来,可用于指导和优化同类空调通风气流组织设计.     

置换通风的数值模拟

首先介绍了用于数值模拟置换通风室内空气分布的N点风口模型、新零方程湍流模型，以简化和加快实际工程中对置换通风的模拟计算。然后利用这些模型对某办公室置换通风的室内温度场和速度场进行了模拟，并和实验数据进行对比。结果表明，计算所得速度和温度分布与实测值吻合得很好，所用的风口模型和湍流模型能将置换通风的温度分层特性和速度场合理地模拟出来，可用于指导和优化置换通风气流组织设计。     

多层建筑屋顶风场的数值模拟

基于计算流体动力学软件Fluent17.2,以浙江温州大学某栋带女儿墙的多层试验教学楼为研究对象,采用RNG k-ε湍流模型对其进行数值模拟.通过对不同风向角下的数值模拟结果进行对比分析,探究屋顶风场与未受扰来流风场的区别,屋顶不同高度风场的差异,以及不同风向角来流情况下近屋面风场的分布特点、变化规律.结果表明:在离屋...     

液氮过冷流动沸腾数值模拟中的双流体模型

液氮过冷流动沸腾广泛见于各种低温换热设备中.采用双流体模型分析液氮过冷流动沸腾,需要为模型提供适当的封闭方程,用于描述汽液两相间质量、动量以及能量的传输过程,封闭方程的合理性直接决定了双流体模型的准确性.将液氮流动沸腾通道划分为近壁区和主流区,分别介绍液氮核态沸腾壁面上的传热传质机理模型,以及两相流程内汽液之间的相互作用的相间传输模型,建立液氮流动沸腾过程适用的双流体模型,并分析了对模型预测能力具有显著影响的因素,指出存在的问题和解决方案.      

井内流动与传热的三维耦合数值模型

钻井过程中,井内的压力、流场及温度分布是钻井作业必需的重要参数。本文提出了一个描述井内流动和传热的三维数学模型,用SIMPLER方法获得了该数学模型的数值解。该数值模型考虑了温度、压力对钻井液的流变性、密度以及各种介质的热物理参数的影响。现场试验数据与数值模拟结果的对比表明：本文提出的数学模型及采用的数值计算方法是可行的。     

基于CFD-DEM的高炉风口回旋区形状和传热特性数值模拟研究

采用计算流体力学-离散单元法(CFD-DEM)耦合方法,对高炉风口回旋区进行了数值模拟研究。首先通过与实验结果对比,验证了CFD-DEM模型的正确性;然后考察了不同气速对风口回旋区形状和传热特性及颗粒接触的影响。数值模拟结果表明:风口回旋区的大小和形状均受气速影响较大,在较大进气速度下,颗粒受到的曳力大于颗粒间的摩擦阻力并破坏颗粒间的桥力,形成较大尺寸的回旋区;且颗粒之间接触力较小,形成较大的空隙结构,更有利于热气体向周围扩散以强化传热。目前考察的三种气速结果表明:当气速为11m/s时,热量向下方传递速度最快;当气速为13m/s时,热量向上方传递速度最快;而当气速为15m/s时,热量向右方传递速度最快;此外,气速越大流态化越明显,颗粒间接触越少,接触力也越小。     

分离流动数值模拟的几个问题

介绍了分离流动现象.给出了三维定常粘性流动分离识别准则,研究了分离流动的性态.给出了应用于分离流动的六种模型及有关的控制方程.讨论了求解NS方程的数值方法和数值模拟的前景.特别研究了应用于分离流动的Eulcr方程的数值模拟问题.      

土体颗粒物流动物质点法模拟的弹塑性和非牛顿流体本构模型比较研究

土体颗粒物流动是一种典型的大变形破坏,具有非牛顿流体的流动特征。准确模拟土体颗粒物的流动及冲击过程,对滑坡和泥石流等地质灾害的防治具有重要意义。物质点法是一种无网格粒子类方法,已在各类大变形问题中得到了广泛应用。以往土体颗粒物流动的模拟,通常采用弹塑性本构模型,但缺乏对非牛顿本构模型的模拟分析。本文引入非牛顿本构模型的模拟分析,旨在为土体颗粒物流动模拟提供一种新的方法与思路。非牛顿本构模型的模拟分析是将非牛顿广义Cross模型引入三维物质点法,通过人工阻尼力模拟颗粒间的摩擦力,对土体颗粒物的坍塌、沿斜面滑动以及冲击障碍物等问题进行了动态模拟,研究了其运动全过程,并与弹塑性本构模型的模拟结果进行了对比验证。结果表明,基于非牛顿流体本构模型的物质点法可以较好地模拟土体颗粒物加速、减速到再次稳定的流动全过程及其对障碍物的冲击效应。     

超空泡形态及其流动特性的数值模拟

对圆盘空化器分别采用CFD、"1/3法则"和空泡截面独立扩展原理三种不同方法数值模拟了水下航行体定常自然超空泡外形及其流动特性。应用CFD方法基于粘性多相流的空泡捕捉法,采用六面体网格,选择Singhal空化模型和SST湍流模式,数值求解均质超空泡流场RANS方程。研究表明:泡形态时变特性是一种行之有效的工程估算方法,应用空泡截面独立扩展原理其计算结果与CFD方法吻合较好,说明了CFD方法用于超空泡流动仿真计算的可行性和独立性原理快速估算超空泡形态的准确性;同时超空泡外形主要与头部空化器有关,空泡长度会由于航行体本身存在变长;可优先选择空泡截面独立扩展原理对水下航行体超空泡外形进行快速估算。     

空气炮碰撞实验台发射过程的数值模拟

利用计算流体动力学方法对单级空气炮碰撞实验台的发射过程进行了数值仿真,所得碰撞车发射速度与实验结果吻合。在此基础上,对碰撞车发射过程的流场变化、车体前后压力以及储气罐压力变化进行分析。结果表明,泄漏气体先于碰撞车充满整个发射空间,形成初始流场,使得车前压力出现正负交替现象,但其数值较小,对车体加速过程的影响可以忽略。当碰撞车进入泄压段后,受冲击射流作用,碰撞车仍处于加速状态,且速度增量约为2 m/s。

     

带自由表面水流与离散颗粒耦合模型研究

基于考察泥沙运动的细观行为特征,采用离散单元法(DEM)模拟泥沙颗粒运动,结合带自由表面的水动力学计算模型,建立了CFD-DEM耦合数值模型。计算程序开发基于Fortran语言来实现。耦合模型中实现了硬球模型和软球模型两种颗粒碰撞模型,应用范围较广。作为自由表面水流与泥沙颗粒流数值模型的初步研究,在模型建立的基础上,对模型做了基本的验证。分别通过单颗粒静水沉降和混合颗粒群分选两个计算工况,验证了模型的正确性及模拟精度。该耦合模型可进一步丰富带自由表面水流条件下泥沙运动的研究手段。     

A correction for Reynolds-averaged-Navier–Stokes turbulence model under the effect of shock waves in hypersonic flows

Studies on the unphysical increase of turbulent quantities for RANS simulation induced by shock waves in hypersonic flows are carried out. Numerical experiments on the hypersonic flow over a blunt body reveal that the phenomenon of unphysical increase of turbulent quantities across the detached shock wave is induced by the strain-rate-based production terms of the k-$\omega$ and k-$\omega$ SST turbulence models, which leads to the over-prediction of aerothermal prediction. While this phenomenon does not occur for Spalart–Allmaras (S–A) turbulence model because of its vorticity-based production term. In order to eliminate this unphysical phenomenon, and to maintain the accuracy of the original models for boundary layer and separation flows, a new correction method for the k-$\omega$ and k-$\omega$ SST models is proposed: by comparing the orders of magnitude between the strain-rate-based and vorticity-based production terms, the vorticity-based production term is used near the shock waves, while the original strain-rate-based production term is still used in other regions. Finally, the correction method is applied to turbulence and transition flows over blunt bodies, and the numerical results show that the correction method effectively eliminates the unphysical increase of turbulent quantities across shock waves and improves the accuracy of aerothermal and transition onset location prediction.     

不同风场下TNT炸药爆炸烟云的扩散模型及特性

为获取不同风场下TNT爆炸烟云扩散时空分布规律与高度变化模型,本文理论描述了爆炸烟云扩散过程与机理,开展了不同水平风速下烟云扩散的计算流体力学（computational fluid dynamics, CFD）仿真和外场时空分布实验,建立了不同水平风速下烟云高度随时间变化模型及烟云最终高度计算模型,分析了烟云扩散过程中形态、温度、密度、速度变化规律。研究结果显示：CFD方法仿真烟云分布结果与实验结果基本一致,大气稳定且无风条件下烟云高度随时间呈指数0.5的幂函数关系,最终高度与爆炸当量可拟合为指数0.47的幂函数模型;水平风会加快烟云与空气混合的速度,导致幂函数模型中指数参数随风速变大而呈线性减小规律,风速越大烟云上升速度衰减越快、上升时间越短、最终高度越低。

     

Multi‐resolution analysis for high accuracy and efficiency of Euler computation

A multi‐resolution analysis (MRA) is proposed for efficient flow computation with preserving the high‐order numerical accuracy of a conventional solver. In the MRA process, the smoothness of a flow pattern is assessed by the difference between original flow property values, and the values approximated by high‐order interpolating polynomial in decomposition. Insignificant data in smooth region are discarded, and flux computation is performed only where crucial features of a solution exist. The reduction of expensive flow computation improves the overall computational efficiency. In order to maintain the high‐order accuracy, modified thresholding procedure restricts the additional error introduced by the thresholding below the order of accuracy of a conventional solver. The practical applicability of the MRA method is validated in various continuous and discontinuous flow problems. The MRA stably computes the Euler equations for continuous and discontinuous flow problems and maintains the accuracy of a conventional solver. Overall, it substantially enhances the computational efficiency of the conventional third‐order accurate solver. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical investigation of a water flow in a rib-roughened channel by using Reynolds stress model

In this paper, water flow in a rib-roughened channel is investigated numerically by using Reynolds stress turbulence models (RSM) on a three-dimensional (3-D) domain. Computational results for mean streamwise velocity component and turbulent kinetic energy show good agreements with available experimental data. Five rib pitch-to-height ratios (p/h) of 1, 5, 10, 15 and 20 are analysed for six different Reynolds numbers (Re) of 3000, 7000, 12,000, 20,000, 40,000 and 65,000. Velocity vectors, streamlines and Reynolds stresses are showed for these ratios and Re numbers. Streamlines revealed that Reynolds numbers do not affect flowfield but play an important role in the Reynolds stresses.     

NUMERICAL SIMULATION OF SWIRLING FLOWS IN OXIDATION REACTORS FOR TiO2 MANUFACTURE

The oxidation reactor plays a key role in producing rutile titanium dioxide （TiO2） from vapor-phase titanium tetrachloride （TiCl4） by employing a swirling flow operation for enhanced gas mixing. This work aims to understand the effect of reactor configuration on the 3-D swirling flow field using computational fluid dynamics （CFD） simulation. Considering the anisotropic turbulence involved, the Reynolds stress model is applied to describe the complex swirling flow together with the cross-flow mixing of gases. The results show significant effect of the flow angle between the wall jet of air stream （representing TiCl4 in practice） and the axial direction on the initial flow field of cross-flow mixing, where 60° gives smooth profiles of axial velocity development while 90° may provide the fastest mixing between the jet and the axial bulk flow. The pipe shape for the reaction and developing zone, i.e., straight, expanding and shrinking, shows slight influence on the hydrodynamics.     

Characteristic analysis of mechanical thermal coupling model for bearing rotor system of high-speed train

Applied Mathematics and Mechanics - Based on Newton’s second law and the thermal network method, a mechanical thermal coupling model of the bearing rotor system of high-speed trains is...