两种颗粒湍流扩散模型数值模拟气液两相流泄漏扩散的比较

采用相间耦合的欧拉-拉格朗日方法,模拟了装有液化气（丙烷）的容器出现小孔或裂缝时,发生泄漏后的气液两相扩散过程。分别应用随机轨道模型及颗粒群模型来考察湍流对液滴扩散的影响,并与确定性轨道进行了比较。给出了数学物理模型,计算结果与实验数据进行了对比。结果表明：采用随机轨道模型能较好地描述液滴的湍流扩散,适用于有液相蒸发的两相流扩散问题。     

关于二相流、多相流、多流体模型和非牛顿流等概念的探讨

本文分析了单相流、二相流和多相流等概念上的差异,也分析了单流体模型、双流体模型和多流体模型等概念上的差异,指出前面三种概念是按流动介质的客观物理构成划分的,而后者是按主观采用的研究方法划分的．目前这些概念在使用中存在一些混乱,如二相流与多相流,多相流与多流体模型等．本文还研究了扩散模型、非牛顿流模型和颗粒流模型等,指出前两种模型在分类上属于单流体模型,分析了非牛顿流模型、扩散模型和双（多）流体模型的特点和应用范围,最后,以泥石流为例讨论了以上概念的应用．      

二维数值模拟爆轰中化学反应率的应用

续定常爆轰数值模拟中化学反应率与人为粘性的相关性,本文对体积起爆函数进行了一维、二维系统考查,并与Cochran反应率做了二维对比计算给出数值结果。数值模拟爆轰的复杂相互作用,应用体积起爆函数为好。     

气液两相流泄漏扩散的数值模拟

采用Eulerian-Lagrangian方法数值模拟过热液体丙烷从喷嘴释放后射流的发展过程。其中,对丙烷蒸汽与空气的气相混合物采用欧拉法求解,利用拉格朗日法跟踪丙烷液滴的运动轨迹,并加入颗粒随机轨道模型来考虑液滴的湍流扩散效应,数值计算所得的物理参量与实验测量结果相吻合。在此基础上,本文将此数值模型和方法应用于开放空间、常温高压下储存液化气(丙烷)的容器发生小孔泄漏时,两相射流在大气中的扩散过程,计算结果符合实际情况。文中还讨论了影响两相流扩散过程的主要因素,为预测可能发生的火灾爆炸事故提供理论依据。     

交联葡聚糖凝胶中反常扩散的Hausdorff分形导数扩散成像模型

本文通过考虑复杂介质的分形特征,并结合交联葡聚糖凝胶的磁共振信号,给出了基于Hausdorff分形导数的磁共振信号衰减模型.该模型包含三个参数,其中分形导数的阶数可以表征介质的分形特性,扩散系数可以表征介质的复杂程度.通过交联葡聚糖凝胶的磁共振实验验证了模型的有效性.结果表明,模型参数的大小与介质内部结构和外部结构有关.与传统的模型相比,基于分形导数的模型参数更符合介质的物理特性,可以准确地描述水分子在复杂介质中的反常扩散,并用于区分不同介质结构特征.     

固体炸药爆轰与惰性介质相互作用的一种扩散界面模型

提出一种保持热力学一致性的扩散界面模型,用来数值模拟固体炸药爆轰与惰性介质的相互作用问题。基于混合网格内各组分物质间可以达到力学平衡状态而不能达到热学平衡状态的假设,由混合网格能量守恒以及压力相等条件,推导出每种组分物质的体积分数演化方程。由此获得的扩散界面模型包括组分物质的质量守恒方程、混合物质的动量及总能量守恒方程,同时包括组分物质的体积分数演化方程和混合物质的压力演化方程。该扩散界面模型的主要特点是考虑了化学反应以及热学非平衡的影响。提出的扩散界面模型在物质界面附近不会出现物理量的非物理振荡现象、适用于任意表达形式的物质状态方程以及任意数目的惰性介质。

     

边界元方法模拟二维势流问题时角点处处理方法的研究

研究了用三节点等参曲线单元时边界元方法模拟二维势流问题中角点处法向量不连续的处理方法.该方法的核心思想是在任意一点处速度矢量是唯一确定的且在角点处速度矢量的方向是已知的.该方法的特点是不改变所求解的线性方程组的大小,不降低计算精度.通过模拟二维平行管道中直匀流绕圆柱的流动验证了该方法的可行性.     

求解二维污染扩散方程的时空任意阶的三点紧致显格式

通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间二阶导数任意精度的三点紧致的表达式,并在半高散方程中通过二维扩散方程本身把时间导数转换为空间导数,从而推导出了时空任意阶的三点紧致显格式.数值实验表明,本文格式的精度很高,而且具有使用简单,易于编程的优点,对求解二维污染扩散方程具有很好的应用前景.     

基于WENO-THINC/WLIC模型的水气二相流数值模拟

水气二相流与诸多领域的实际工程问题密切相关.对二相流运动进行高精度的数值模拟是计算流体力学研究的难点和热点.针对开敞水域的自由表面流运动问题,将水和空气均视为不可压缩流体,采用五阶加权基本无震荡(weighted essentially non-oscillatory,WENO)格式求解描述流体运动的纳维斯托克斯(Na...     

基于五阶WENO格式的燃气在药床中流动过程二维两相流研究

为研究内弹道初始阶段中心点火管燃气在膛内药床中的流动特性和传播规律,设计了可视化点传火实验平台,并进行了膛内假药床的点传火实验。基于加权本质无震荡(weighted essentially non-oscillatory, WENO)格式,构造了膛内轴对称二维内弹道两相流模型,对膛内燃气在假药床中的流动过程进行数值模拟。计算结果与可视化实验结果符合较好,全局压力平均误差为5.35%。表明数值计算准确地描述了燃气流动特性,完整地呈现了点火管燃气在假药床中的发展过程。在点火初始阶段,膛内压力径向效应明显,气相沿径向传播较快,药床药粒基本不会发生运动;随着燃气逐渐在膛内传播,膛内压力呈现径向一致、轴向梯度分布的特征,在压力梯度作用下,气相轴向速度开始占据主导,径向速度在膛底和中部区域减小为零,而固相速度随气相速度变化而变化;气相在到达弹底前,由于固相颗粒的壅塞,会提前出现速度反向波动现象。     

Molecular transfer in turbulent flows

For modeling the molecular transfer of a passive scalar in a known turbulent field, the equations for the average scalar value and the correlation function for the scalar field are written in a form which makes it possible to examine the effect of molecular transfer on turbulent transfer and scalar dissipation. For the closure of the equation for the correlation function, the Prandtl hypothesis is used. The statistical reliability of this closure is demonstrated. The system proposed makes it possible to predict the dynamics of a decaying uniform scalar field and to explain why the effect of the real value of the molecular-transfer coefficient on the decaying scalar field is weak. Specific features of the transport process in a plane layer with prescribed scalar values on the layer boundaries are considered.     

Mapping closure approximation for reactive scalars in random flows

The Mapping Closure Approximation (MCA) approach is developed to describe the statistics of both conserved and reactive scalars in random flows. The statistics include Probability Density Function (PDF), Conditional Dissipation Rate (CDR) and Conditional Laplacian (CL). The statistical quantities are calculated using the MCA and compared with the results of the Direct Numerical Simulation (DNS). The results obtained from the MCA are in agreement with those from the DNS. It is shown that the MCA approach can predict the statistics of reactive scalars in random flows. The project supported by the National Committee of Science and Technology, China, under the Special Funds for Major Basic Research Project (G2000077305 and G1999032801), and the National Natural Science Foundation of China (10325211)     

Direct numerical simulation of wall turbulent flows with microbubbles

The marker‐density‐function (MDF) method has been developed to conduct direct numerical simulation (DNS) for bubbly flows. The method is applied to turbulent bubbly channel flows to elucidate the interaction between bubbles and wall turbulence. The simulation is designed to clarify the structure of the turbulent boundary layer containing microbubbles and the mechanism of frictional drag reduction. It is deduced from the numerical tests that the interaction between bubbles and wall turbulence depends on the Weber and Froude numbers. The reduction of the frictional resistance on the wall is attained and its mechanism is explained from the modulation of the three‐dimensional structure of the turbulent flow. Copyright © 2001 John Wiley & Sons, Ltd.     

The distribution of solid particles suspended in a turbulent flow: A stochastic approach

Basic fluid mechanics and stochastic theories are applied to show that the concentration distribution of suspended solid particles in a direction normal to the mean streamlines of a two-dimensional turbulent flow is greatly influenced by the lift force exerted on them in the vicinity of the wall. Analytic solution shows that, when the direction of the mean flow is horizontal, the probability density functionp (y, t) for random displacements of the particles will have a maximum value at a point from the wall where the perpendicular component of the lift force precisely balances particle gravity. Interpretation of experimental observations is presented using this theory.     

Differential and algebraic models for the second moments of particle velocity and temperature fluctuations in turbulent flows

Differential and algebraic models are constructed for the dispersed-phase turbulent stresses and heat fluxes and for the mixed moments of the velocity and temperature fluctuations in the continuous and dispersed phases. The models are based on a kinetic equation for the joint probability density of the particle velocity and temperature in an anisotropic turbulent flow. The results are compared with the available direct numerical simulation (DNS) data.     

Assessment of the validity of a log-law for wall-bounded turbulent bubbly flows

There has been considerable discussion in recent years concerning whether a log-law exists for wall-bounded, turbulent bubbly flows. Previous studies have argued for the existence of such a log-law, with a modified von Kármán constant, and this is used in various modelling studies. We provide a critique of this idea, and present several theoretical reasons why a log-law need not be expected in general for wall-bounded, turbulent bubbly flows. We then demonstrate using recent data from interface-resolving Direct Numerical Simulations that when the bubbles make a significant contribution to the channel flow dynamics, the mean flow profile of the fluid can deviate significantly from the log-law behaviour that approximately holds for the single-phase case. The departures are not surprising and the basic reason for them is simple, namely that for bubbly flows, the mean flow is affected by a number of additional dynamical parameters, such as the void fraction, that do not play a role for the single-phase case. As a result, the inner/outer asymptotic regimes that form the basis of the derivation of the log-law for single-phase flow do not exist in general for bubbly turbulent flows. Nevertheless, we do find that for some cases, the bubbles do not cause significant departures from the unladen log-law behaviour. Moreover, we show that if departures occur these cannot be understood simply in terms of the averaged void fraction, but that more subtle effects such as the bubble Reynolds number and the competition between the wall-induced turbulence and the bubble-induced turbulence must play a role.     

On basic equations of turbulent swirling gas-solid flows and their application in cyclones

The basic equations of turbulent gas-solid flows are derived by using the pseudo-fluid model of particle phase with a refined two-phase turbulence model. These equations are then applied to swirling gas-particle flows for analyzing the collection efficiency in cyclone separators.     

Schlieren measurement of turbulent structure in a diffusion flame

The regular and random mixing structures in a turbulent diffusion flame were investigated using the quantitative, dynamic crossed-beam schlieren method. Evidence was found close to the nozzle relating to the vortexlike structure of eddies surrounding the central fuel jet flow. The observations also make possible resolution of turbulent intensity, scales, convection, and spectra within the diffusion flame without the use of seeding or intrusion of measuring probes. It is found that length scales and other turbulence parameters in the diffusion flame progressively revert to values similar to those expected and observed in scalar passive mixing as the combustion reaction intensity reduces with axial distance from the nozzle system.