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

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

液滴在气体介质中剪切破碎的数值模拟研究

液滴变形和破碎是燃料抛撒问题的重要过程.本文将VOF方法和标准k-ε湍流模型组合,建立了计算液滴在气流中变形破碎过程的数值方法.数值模拟了相关的实验,计算得到的液滴破碎过程与实验结果符合较好.在此基础上,分析了几个关键参数(Weber数、Ohnesorge数、液气密度比)对液滴破碎过程的影响.计算结果表明,Weber数...     

云状空化非定常脱落机理的数值与实验研究

结合数值计算和实验技术研究了云状空化的非定常脱落机理. 实验采用高速录像技术观察了绕Clark-y型水翼云状空化形态随时间的变化. 数值计算采用了汽-液两相的当地均相流模型,湍流封闭采用了一种修正的RNG k-\varepsilon方程,利用商业软件的二次开发技术,引入了一种与空化区域水汽相密度相关的系数,对湍流模型进行了修正. 计算和实验结果均表明:云状空化尾部存在着准周期性的涡状空化团的脱落;局部压强的增加是引起空穴断裂进而脱落的直接原因;压强升高是由于近壁处的反向射流引起的;空穴尾部的大规模的旋涡运动引起的近壁处的逆压梯度是引起反向射流的主要原因.      

不同界面计算方法对液滴在电场作用下变形数值模拟精度的影响

液滴在电场作用下的变形是电流体动力学的基础课题之一，表面张力的计算精度对液滴变形量的模拟结果有重要影响。本文以开源计算流体动力学平台OpenFOAM的VOF模型为框架，研究了MULES和isoAdvector两类界面更新算法与相分数梯度和RDF函数两类曲率算法对电场作用下液滴变形模拟精度的影响。研究表明，isoAdvector算法相比MULES算法对网格密度的要求更低，但其耦合相分数梯度算法计算表面张力的误差较高。isoAdvector算法耦合RDF函数算法计算误差较低，并且在使用轴对称网格时，只有该算法能够同时处理液滴平行于电场和垂直于电场方向的变形，得到的数值结果与解析解吻合较好。     

气固两相流场的湍流颗粒浓度理论模型

本文进行了气固两相流动颗粒湍流扩散现象的理论分析，提出了颗粒湍流扩散系数和气流弥散效应二个颗粒湍流模化新概念，在此基础上建立了气固两相流场湍流颗粒浓度模型。理论模型包括离心力和其它外加力场作用下颗粒运动和浓度分布的计算方法。运用湍流颗粒浓度模型，对直管气固两相流动、受限射流气固两相流动和９０°弯管气固两相流动等三种流动做了数值模拟，计算获得颗粒速度、颗粒浓度等主要流动参数。讨论了湍流颗粒浓度模型的适用性。     

SPH方法在聚能装药射流三维数值模拟中的应用

采用有限元软件LS-DYNA中的SPH算法实现聚能装药射流形成过程的三维数值模拟.将SPH方法和FEM方法与理论计算结果相比较.研究结果表明:SPH数值模拟计算过程十分稳定,避免了有限元数值模拟过程中的网格扭曲、缠绕和物质穿透等问题,而且计算精度比有限元方法更高;采用SPH方法计算的射流速度、射流长度与有限元计算结果和...     

气相中双液滴正碰过程数值模拟研究

液滴碰撞现象普遍存在于动力装置燃烧室喷嘴的下游区域,影响燃料的雾化性能。为了揭示相同直径的双液滴中心碰撞机理,求解了轴对称坐标系下的N-S方程,采用VOF(Volume of Fluid)方法捕捉液滴碰撞过程中气液自由表面的演化规律。利用Qian等提供的实验结果对计算模型进行数值校验,验证了模型的准确性。在此基础上,研究了环境压强对液滴碰撞反弹后不同结果(分离和融合)的影响,分析了环境压强和Weber数对液滴碰撞分离的影响。结果表明,液滴在碰撞反弹后的状态(分离或融合)是由液滴间气膜压强与环境气动阻力共同作用的结果,环境压强对液滴碰撞分离过程基本没有影响;Weber数越大,碰撞过程中变形的幅度越大。     

用格子Boltzmann方法模拟液滴撞击固壁动力学行为

首次用格子Boltzmann方法中的伪势模型对液滴撞击固壁的动力学行为进行了数值模拟.详细研究了液滴在壁面上的流动状态以及各种因素对撞击过程的影响.通过数值模拟得到:壁面的可润湿性越小,液滴越容易发生反弹,液滴的回缩速度越快;液滴的撞击速度越大,所得到的相对直径越大,回缩速度越快;液滴的粘性越小,所得到的相对直径越大;液滴的表面张力越大,液滴越容易发生反弹现象.另外,液滴的最大相对直径与We数满足一定的线性关系,这些结果与前人的理论预测和实验结果完全吻合.     

同轴旋转可压缩流动中液体射流稳定性

基于线性稳定性理论,建立了描述同轴旋转可压缩流动中超空化条件下液体射流稳定性的数学模型,并对数学模型及其求解方法进行了验证;在此基础上,对模型中考虑的射流及气体可压缩性、气体同轴旋转以及超空化等因素对射流稳定性的影响进行了分析. 分析结果表明,模型中考虑射流及气体的可压缩性后,与不考虑可压缩性相比,计算得到的射流稳定性明显变差,最小液滴直径减小,分裂液滴直径变化范围变宽,且小液滴数量增多. 气体的同轴旋转在轴对称与非轴对称扰动下对射流稳定性的影响完全相反;轴对称扰动时,气体旋转使射流稳定性增强,而非轴对称扰动时则正好相反;气体旋转有可能导致影响射流稳定性的扰动模式发生根本性变化. 超空化使射流稳定性变差;超空化程度较弱时,超空化使分裂液滴最小直径减小,分裂液滴直径变化范围增大;而超空化达到一定程度后,进一步提高超空化程度,分裂液滴最小直径几乎保持不变.     

用火焰面模型模拟甲烷/空气湍流射流扩散火焰

以层流对撞扩散火焰为基础，利用层流火焰面模型(laminar flamelet model)的方法生成层流火焰面数据库，分别采用预先设定的几率密度函数(propabality density function, PDF)模型和混合物分数-湍流频率的联合几率密度函数输运模型，将火焰面方法应用于甲烷/空气湍流射流扩散火焰结构的模拟计算中.两个模型的计算结果和实验结果进行了比较和分析.     

A PDF propagation approach to model turbulent dispersion in swirling flows

A computationally efficient approach that solves for the spatial covariance matrix along the dense particle ensemble-averaged trajectory has been successfully applied to describe turbulent dispersion in swirling flows. The procedure to solve for the spatial covariance matrix is based on turbulence isotropy assumption, and it is analogous to Taylor's approach for turbulent dispersion. Unlike stochastic dispersion models, this approach does not involve computing a large number of individual particle trajectories in order to adequately represent the particle phase; a few representative particle ensembles are sufficient to describe turbulent dispersion. The particle Lagrangian properties required in this method are based on a previous study (Shirolkar and McQuay, 1998). The fluid phase information available from practical turbulence models is sufficient to estimate the time and length scales in the model. In this study, two different turbulence models are used to solve for the fluid phase – the standard k–ε model, and a multiple-time-scale (MTS) model. The models developed here are evaluated with the experiments of Sommerfeld and Qiu (1991). A direct comparison between the dispersion model developed in this study and a stochastic dispersion model based on the eddy lifetime concept is also provided. Estimates for the Reynolds stresses required in the stochastic model are obtained from a set of second-order algebraic relations. The results presented in the study demonstrate the computational efficiency of the present dispersion modeling approach. The results also show that the MTS model provides improved single-phase results in comparison to the k–ε model. The particle statistics, which are computed based on the fundamentals of the present approach, compare favorably with the experimental data. Furthermore, these statistics closely compare to those obtained using a stochastic dispersion model. Finally, the results indicate that the particle predictions are relatively unaffected by whether the Reynolds stresses are based on algebraic relations or on the turbulence isotropy assumption.     

A Lagrangian–Eulerian model of particle dispersion in a turbulent plane mixing layer

A Lagrangian–Eulerian model for the dispersion of solid particles in a two‐dimensional, incompressible, turbulent flow is reported and validated. Prediction of the continuous phase is done by solving an Eulerian model using a control‐volume finite element method (CVFEM). A Lagrangian model is also applied, using a Runge–Kutta method to obtain the particle trajectories. The effect of fluid turbulence upon particle dispersion is taken into consideration through a simple stochastic approach. Validation tests are performed by comparing predictions for both phases in a particle‐laden, plane mixing layer airflow with corresponding measurements formerly reported by other authors. Even though some limitations are detected in the calculation of particle dispersion, on the whole the validation results are rather successful. Copyright © 2002 John Wiley & Sons, Ltd.     

A Novel Framework for Integration of Random-Walk Particle-Tracking Simulation in Subsurface Multi-phase Flow Modeling

Coarse-scale models are generally preferred in the numerical simulation of multi-phase flow due to computational constraints. However, capturing the effects of fine-scale heterogeneity on flow and isolating the impacts of numerical (artificial) dispersion, which increases with scale, are not trivial. In this paper, a particle-tracking method is devised and integrated in a scale-up workflow to estimate the conditional probability distributions of multi-phase flow functions, which can be considered as inputs in coarse-scale simulations with existing commercial packages. First, a novel particle-tracking method is developed to solve the saturation transport equation. The transport calculation is coupled with a velocity update, following the implicit pressure, explicit saturation framework, to solve the governing equations of two-phase immiscible flow. Each phase particle is advanced in a deterministic convection step according to the phase velocity, as well as in a stochastic dispersion step based on the random Brownian motion. A kernel-based formulation is proposed for computation of fluid saturation in accordance with the phase particle distribution. A novel aspect is that this method employs the kernel approach to construct saturation from phase particle distribution, which is an important improvement to the conventional box method that necessitates a large number of particles per grid cell for consistent saturation interpolation. The model is validated against various analytical solutions. Finally, the validated model is integrated in a statistical scale-up procedure to calibrate effective, or “pseudo,” multi-phase flow functions (e.g., relative permeability functions). The proposed scale-up framework does not impose any length scale requirement regarding the distribution of sub-grid heterogeneities.     

Large eddy simulation of polydispersed inertial particles using two-way coupled PDF-PBE

This paper presents large eddy simulation (LES) of polydispersed turbulent recirculating flows using a two-way coupled probability density function of the population balance equation (PDF-PBE). A stochastic Monte Carlo method is adopted to solve the PDF-PBE on an ensemble of notional Lagrangian particles and the method of Stokes binning, that was recently developed by Salehi et al. (2017) is employed to explicitly treat effects of inertia. The PDF-PBE is applied to the experiment of Boŕee et al. (2001) which studied dispersion of polysized inertial particles in a bluff body configuration. The particle mass loading is 22% where the dispersed elements affect the carrier phase velocity. The simulations are performed using both the standard Smagorinsky and the Wall-Adapting Local Eddy-viscosity (WALE) subgird turbulent models. It is found that the subgrid model has a significant impact on the results, particularly on the carrier velocity. The WALE model shows a better agreement with the measurements. Different boundary conditions are tested for injection of notional particles. It is demonstrated both dispersed and carrier phase velocities are initially sensitive to the particle boundary condition whereas the difference between tested conditions becomes marginal further downstream. The best results are obtained for the particle boundary condition that accounts for effects of inertia at the inlet. Finally, the sensitivity of the PDF-PBE simulations to the number of Stokes bins is studied. It is found that eight Stokes bins are enough to accurately model the polysized particles dispersion in a bluff body configuration. The results indicate that particle dispersion is notably sensitive to the number of Stokes bins but the particle velocity predictions are much less sensitive with small variation in velocity results over the range of Stokes bins that were tested.     

A Fickian diffusion model for the spreading of liquid plumes infiltrating in heterogeneous media

Infiltration of water and non-aqueous phase liquids (NAPLs) in the vadose zone gives rise to complex two- and three-phase immiscible displacement processes. Physical and numerical experiments have shown that ever-present small-scale heterogeneities will cause a lateral broadening of the descending liquid plumes. This behavior of liquid plumes infiltrating in the vadose zone may be similar to the familiar transversal dispersion of solute plumes in single-phase flow. Noting this analogy we introduce a mathematical model for ‘phase dispersion’ in multiphase flow as a Fickian diffusion process. It is shown that the driving force for phase dispersion is the gradient of relative permeability, and that addition of a phase-dispersive term to the governing equations for multiphase flow is equivalent to an effective capillary pressure which is proportional to the logarithm of the relative permeability of the infiltrating liquid phase. The relationship between heterogeneity-induced phase dispersion and capillary and numerical dispersion effects is established. High-resolution numerical simulation experiments in heterogeneous media show that plume spreading tends to be diffusive, supporting the proposed convection-dispersion model. Finite difference discretization of the phase-dispersive flux is discussed, and an illustrative application to NAPL infiltration from a localized source is presented. It is found that a small amount of phase dispersion can completely alter the behavior of an infiltrating NAPL plume, and that neglect of phase-dispersive processes may lead to unrealistic predictions of NAPL behavior in the vadose zone.     

梯度介质半空间Rayleigh面波传播性能研究

对剪切弹性模量沿深度以指数函数变化的非均质半空间，本文用摄动法得到了Rayleigh面波的波函数解答及相速度方程。以不同金属与陶瓷复合而成的几种梯度材料为例，用数值方法求解了相速度方程，给出了相应的波的弥散曲线，结果表明，梯度介质半空间自由表面附近的Rayleigh波通常有两种不同的弥散形式，即正常弥散和非正常弥散。     

薄膜涂层材料中的广义瑞利波

基于弹性动力学的线性理论，建立了涂层材料中广义瑞利波传播的理论分析模型，并 且由波动方程和边界条件推导了波的频散方程．分析了慢层和快层对相速度频散的影响，给 出了不同层厚-波长比和不同涂层-基体密度比情况下广义瑞利波相速度的理论解．算例分 析分别比较了慢层和快层结构中波的相速度、群速度，以及随深度衰减的位移与应力振 幅．另外，相速度曲线和位移振幅曲线与文献中给出的结果吻合，验证了理论模型和分析过 程的正确性．     

Lagrangian particle dispersion in turbulent flow over a wall mounted obstacle

Large-eddy simulations (LES) of particle-laden turbulent flows are presented in order to investigate the effects of particle response time on the dispersion patterns of a space developing flow with an obstruction, where solid particles are injected inside the wake of an obstacle [Vincont, J.Y., Simoens, S., Ayrault M., Wallace, J.M., 2000. Passive scalar dispersion in a turbulent boundary layer from a line source at the wall and downstream of an obstacle. J. Fluid Mech. 424, 127–167]. The numerical method is based on a fully explicit fractional step approach and finite-differences on Cartesian grids, using the immersed boundary method (IBM) to represent the existence of solid obstacles. Two different turbulence models have been tested, the classical Smagorinsky turbulence model and the filtered structure function model. The dispersed phase was modelled either by an Eulerian approach or a Lagrangian particle tracking scheme of solid particles with Stokes numbers in the range St = 0–25, assuming one-way coupling between the two phases. A very good agreement was observed between the Lagrangian and Eulerian approaches. The effect of particle size was found to significantly differentiate the dispersion pattern for the inhomogeneous flow over the obstacle. Although in homogeneous flows like particle-laden turbulent channels near-wall particle clustering increases monotonically with particle size, for the examined flow over an obstacle, preferential concentration effects were stronger only for an intermediate range of Stokes numbers.     

Modeling droplet dispersion and interphase turbulent kinetic energy transfer using a new dual-timescale Langevin model

Dispersion of spray droplets and the modulation of turbulence in the ambient gas by the dispersing droplets are two coupled phenomena that are closely linked to the evolution of global spray characteristics, such as the spreading rate of the spray and the spray cone angle. Direct numerical simulations (DNS) of turbulent gas flows laden with sub-Kolmogorov size particles, in the absence of gravity, report that dispersion statistics and turbulent kinetic energy (TKE) evolve on different timescales. Furthermore, each timescale behaves differently with Stokes number, a non-dimensional flow parameter (defined in this context as the ratio of the particle response time to the Kolmogorov timescale of turbulence) that characterizes how quickly a particle responds to turbulent fluctuations in the carrier or gas phase. A new dual-timescale Langevin model (DLM) composed of two coupled Langevin equations for the fluctuating velocities, one for each phase, is proposed. This model possesses a unique feature that the implied TKE and velocity autocorrelation in each phase evolve on different timescales. Consequently, this model has the capability of simultaneously predicting the disparate Stokes number trends in the evolution of dispersion statistics, such as velocity autocorrelations, and TKE in each phase. Predictions of dispersion statistics and TKE from the new model show good agreement with published DNS of non-evaporating and evaporating droplet-laden turbulent flow.