一种SPH应力修正算法及自由表面流中的应用

提出了一种SPH应力修正算法,即模型中的拉应力和压应力分别采用不同的插值核函数和状态方程来处理,改善应力稳定性问题。介绍了一种改进的Quintic核函数,用于改善模型中压应力的稳定性。通过增加钟型核函数的光滑长度,改善模型中拉应力的稳定性。采用该应力修正算法模拟了无重力条件下方形液滴的震荡变形过程,对比分析了不同算法的模拟结果。此外,为进一步验证算法的适用性,模拟了溃坝算例。研究表明,改进的Quintic型核函数明显改善了粒子聚集现象,该SPH应力修正方法可以使液滴具有更均匀的粒子分布以及更光滑的自由表面,有效改善了SPH方法中的压应力不稳定作用以及自由表面流的模拟精度。     

液滴撞击具有浸润性分布的倾斜固壁铺展行为的格子Boltzmann研究

采用基于单组分多相伪势模型的格子Boltzmann方法,模拟了三维液滴撞击左右两侧浸润性不同的倾斜固壁的铺展过程,获得了液滴在壁面两侧的铺展因子、相对铺展宽度、相对高度和液滴运动速度随时间的变化情况,研究了壁面浸润性分布和壁面倾斜角度对液滴铺展过程的影响．结果表明,液滴在倾斜壁面的铺展过程受到重力和表面力的综合作用,重力影响液滴的铺展和沿壁面向下的滑动,壁面浸润性分布影响液滴向壁面亲水侧横向移动．     

磁场对双液滴运动及加速聚结的影响分析

建立了油液中单个液滴在磁场作用下的力学模型,对双液滴在磁场作用下的移动聚结进行了理论分析;研究了油液中水平双液滴的间距随着磁感应强度和容器倾斜角度的变化规律。结果表明:随着磁感应强度的增大,双液滴的移动距离均增大并且双液滴中心间距为1cm,磁感应强度为400mT～500mT时,双液滴的移动距离最快;随着倾斜角度的增大,双液滴的移动距离呈现波浪变化趋势并有最佳的倾斜角度;容器倾斜一定角度时,液滴在重力分力作用下运动并切割磁感线产生洛伦兹力和电场力,加大了液滴的碰撞几率,从而促进了液滴的聚结。     

不同重力下90°弯管内气液两相流流型及流动特性研究

研究了不同重力条件下90°弯管内气液两相流流型分布形态及流动特性. 通过建立90°弯管内气液两相流流动的三维数学物理模型,采用VOF 方法,对10-6g₀, 10-4g₀, 10-2g₀, 1g₀ (g₀= 9.8m/s2) 重力下的90°弯管内气液两相流流型分布特征、截面空隙率、滑速比及气相尾部最大斜向角进行了比较分析. 研究结果表明:所建立的模型能够正确模拟不同重力条件下90°弯管内气液两相流流型和截面空隙率,并得到气液两相弯管二次流与单相二次流的不同特性. 随着重力水平的提高,90°弯管对气相流型的影响作用减弱,气相整体向弯管内侧积聚靠拢,弯管对尾部的斜向作用减弱.      

不同重力条件下管内冷凝现象研究进展

本文对不同重力环境中水平管内冷凝气液两相流动与传热现象的研究现状进行了全面评介, 重点关注于航天应用中的小管径、低流量和以氨为工质的情形,以及部分重力条件下冷凝现象中的重力效应. 冷凝现象中管壁四周液膜的存在,导致水平管内冷凝气液两相流与沸腾或绝热气液两相流在构型上存在着明显的不同, 表面张力的作用增大,向波状和分层流转换所对应的临界Bond数也随之增大, 进而使小管径、低流量水平管内冷凝气液两相流动与传热现象中的重力效应减弱,甚至可以忽略不计. 实验发现该情形中摩擦压降要小于基于常规气液两相流实验数据的经验关联式的预测结果, 而更接近光滑环状流情形; 对冷凝两相传热系数的预测,只有采用基于空隙率与界面剪切率相互耦合的模型, 才能够给出较好的结果.      

可侵蚀地表上方含尘大气运动的数值模型

在稀相气固两相流的双流体模型框架下，导出含尘大气运动的基本方程，其中的源项考虑相间双向耦合作用；通过对基本方程进行无量纲化，求得控制气相和弥散相动力学行为的相似准则，包括弗洛德数、罗斯贝数、颗粒雷诺数、颗粒惯性参数、颗粒剪切横移参数、颗粒旋转横移参数和沙尘质量载荷率等. 作为一个模型问题，研究可侵蚀地表上方充分发展湍流大气边界层流动对土壤颗粒的气动卷扬作用，其中计及气动阻力、萨弗曼升力和重力等因素的影响. 为了克服由于沙尘跃移运动轨迹交叉导致流动参数多值性所造成的困难，引进拉格朗日坐标下弥散相连续方程. 对两种不同风速和4种不同粒径条件下的沙尘运动进行数值模拟，给出相应情况中沙尘运动特性和密度分布剖面并讨论风速和粒径等参数的影响，还细致分析了含尘大气边界层中沙尘与气流之间的能量传递过程. 所得的方程、准则和方法可以为土壤风蚀和沙尘暴等自然灾害的预报提供理论基础.     

分子动力学模拟金属玻璃板拉伸时剪切带的形成和演化过程

金属玻璃在低温高应力条件下容易生成剪切带而导致结构的破坏,大大限制了它的推广应用。本文采用分子动力学模拟研究了三种Cu64Zr36（不带缺口、一侧带缺口、两侧带缺口）金属玻璃板试样在拉伸过程中剪切带的形成和演化过程及其力学性能。结果表明：不带缺口金属玻璃板试样在低温高应力的拉伸过程中会自发出现局部剪切转变区,发生剪切局部化,继续拉伸会在与加载轴大约成45°方向上形成剪切带。剪切带的形成与剪切转变区的分布和局部化有关,带缺口比不带缺口的试样会更早出现应变局部化,即在较低的拉伸应变下便形成剪切带,其拉伸强度也相应较低。相同条件下,一侧带缺口与两侧带缺口的试样在拉伸强度上几乎相同,但两侧带缺口试样的应变局部化程度稍低,主要是两侧缺口处均出现了剪切转变区,导致其分布和局部化不够集中,这也是形成主剪切带和次剪切带的主要原因。以上结果为进一步从原子尺度理解金属玻璃剪切带的形成和演化特征提供了重要的信息。     

分子动力学模拟金属玻璃板拉伸时剪切带的形成和演化过程

金属玻璃在低温高应力条件下容易生成剪切带而导致结构的破坏,大大限制了它的推广应用。本文采用分子动力学模拟研究了三种Cu64Zr36（不带缺口、一侧带缺口、两侧带缺口）金属玻璃板试样在拉伸过程中剪切带的形成和演化过程及其力学性能。结果表明：不带缺口金属玻璃板试样在低温高应力的拉伸过程中会自发出现局部剪切转变区,发生剪切局部化,继续拉伸会在与加载轴大约成45°方向上形成剪切带。剪切带的形成与剪切转变区的分布和局部化有关,带缺口比不带缺口的试样会更早出现应变局部化,即在较低的拉伸应变下便形成剪切带,其拉伸强度也相应较低。相同条件下,一侧带缺口与两侧带缺口的试样在拉伸强度上几乎相同,但两侧带缺口试样的应变局部化程度稍低,主要是两侧缺口处均出现了剪切转变区,导致其分布和局部化不够集中,这也是形成主剪切带和次剪切带的主要原因。以上结果为进一步从原子尺度理解金属玻璃剪切带的形成和演化特征提供了重要的信息。     

垂直上升气液环状流流动参数的数值计算

提出了一个新的气核-液膜耦合模型来求解垂直上升气液环状流在充分发展段的流动参数.本模型考虑了液膜、气核以及它们之间的相互影响和作用.模型中基本的气核区域和液膜区域的质量和动量方程由Fluent6.3.26进行求解,而液滴方程以及相界面上的夹带和沉积作用通过用户自定义接口函数UDF(User Defined Functi...     

气液二相流研究概述

气液二相流研究液体和气体(或蒸气)两相介质共存条件下的流动特性。二相体系可以是液体中含有气体或(蒸气)泡或者气体,(或蒸气)中含有液体微滴,也可能因其中气泡或液滴的聚并,两相间形成更复杂的分布状态。气液二相流是自然界、日常生活中常见的现象,在许      

A New Mathematical Model for Force Gravity Drainage in Fractured Porous Media

In force gas/oil gravity drainage process in fractured porous media, gas is flowing in both matrix and fractures leading to produce a finite gas pressure gradient. Consequently, viscous force plays an important role for displacing matrix oil toward fractures in addition to gravity force that is required to be modeled appropriately. A new analytical model for estimation of steady state oil saturation distribution with assumption of fixed gas pressure gradient throughout the matrix is presented. Moreover, based on some results of this analytical model a different numerical formulation is developed to predict the performance of oil production process. Comparison of the results obtained from this numerical model with the results of a conventional simulator demonstrates that the newly developed model can be applied with satisfactory accuracy. Numerical simulations show that the viscous displacement in fractured porous media can reduce the capillary threshold height, and thus it suggests the force gravity drainage as a favorable production mechanism when the matrix length is close to the threshold height.     

Optical probe visualization of air–water flow structure through a sudden area contraction

 Single-fiber optical probes were used to investigate the time-averaged structure of gas–liquid horizontal flow through a sharp-edged sudden area contraction. The probes allowed to measure the local void fraction distribution over several cross sections of a pipe having an inner diameter of 0.08 m upstream and 0.06 m downstream of the sudden contraction. The water mass flow rate was 3 kg/s, while the gas fraction of the volume flow ranged from 0.2 to 0.8. The local void fraction was plotted as a function of its two spatial coordinates, so that a representation of the time-averaged gas distribution over the cross section could be obtained. The contraction was shown to considerably alter the distribution of the phases, so that the correlations for straight pipes appear no longer suitable. Received: 27 August 2001 / Accepted: 19 November 2001     

Size prediction of κ-carrageenan droplets formed in co-flowing immiscible liquid

The formation of κ-carrageenan droplets in channel emulsification was experimentally investigated. The dispersed phase was vertically injected into co-flowing immiscible palm oil in the direction of gravity. This study focused on predicting κ-carrageenan drop size using force balance analysis. The force balance model considers the interfacial tension to be the solitary attaching force, while a combination of the drag force from the co-flowing palm oil and the body force of the extruding κ-carrageenan liquid act as the detaching forces. The conventional model gave poor predictions for droplet size, with an average relative deviation of 23%. This large deviation could be attributed to necking phenomena and an underestimation of the drag force generated on the shear-thinning κ-carrageenan solution. By incorporating correction factors, the average relative deviation of the force balance model dropped to 4%.     

Study of two-phase flows in reduced gravity using ground based experiments

Experimental studies have been carried out to support the development of a framework of the two-fluid model along with an interfacial area transport equation applicable to reduced gravity two-phase flows. The experimental study simulates the reduced gravity condition in ground based facilities by using two immiscible liquids of similar density namely, water as the continuous phase and Therminol 59^® as the dispersed phase. We have acquired a total of eleven data sets in the bubbly flow and bubbly to slug flow transition regimes. These flow conditions have area-averaged void (volume) fractions ranging from 3 to 30% and channel Reynolds number for the continuous phase between 2,900 and 8,800. Flow visualization has been performed and a flow regime map developed which is compared with relevant bubbly to slug flow regime transition criteria. The comparison shows that the transition boundary is well predicted by the criterion based on critical void fraction. The value of the critical void fraction at transition was experimentally determined to be approximately 25%. In addition, important two-phase flow local parameters, including the void fraction, interfacial area concentration, droplet number frequency and droplet velocity, have been acquired at two axial locations using state-of-the-art multi-sensor conductivity probe. The radial profiles and axial development of the two-phase flow parameters show that the coalescence mechanism is enhanced by either increasing the continuous or dispersed phase Reynolds number. Evidence of turbulence induced particle interaction mechanism is highlighted. The data presented in this paper clearly show the marked differences in terms of bubble (droplet) size, phase distribution and phase interaction in two-phase flow between normal and reduced gravity conditions.     

Simultaneous entrainment of oil and water droplets in high Reynolds number gas turbulence in horizontal pipe flow

We develop a 1D cross sectional concentration profile model for oil and water droplets that coexist in the turbulent gas phase (of Re ∼ 10⁶) in near horizontal stratified pipe flows. Entrainment of the oil and water mixture from a liquid film near the bottom of the pipe into the gas is modeled based on earlier single-fluid entrainment correlations. A Gamma distribution for the droplet sizes based on the breakup of liquid filaments, is adopted. An explicit algebraic–exponential formula for the total concentration profile for either phase can then be derived.     

A droplet transport model for channel and pipe flow based on particle kinetic theory and a stress–ω turbulence model

A model for droplet transport in turbulent gas–liquid flow in horizontal channels or pipes is developed. We invoke the kinetic theory of Reeks and co-workers for inertial particles, and the Reynolds stress–ω model of Wilcox for the gas turbulence. We introduce a suitable turbulence boundary condition at the gas–liquid interface. The full kinetic model is compared to experimental data using high density SF₆ gas and Exxsol oil, and the model gives a good prediction of the droplet distribution.     

关于颗粒悬浮机理和悬浮动的讨论

从分析气体分子的悬浮和静水中Ｂｒｏｗｎ微粒的悬浮之机理出发,论述了重力场中粒子（分子、微粒等）的悬浮不一定需要其它外力,粒子本身的任何形式的无规则运动,达到一定强度后都能使粒子弥散悬浮．河流中的泥沙颗粒和气（水）力输送管道中的颗粒的悬浮也主要靠颗粒物的无规则运动．作用于颗粒的升力和其它力可改变颗粒悬浮沿高度的分布,但仅用这些力（若无任何无规则运动）无法解释颗粒的弥散悬浮状态．讨论了颗粒对流动阻力的双重作用：支持颗粒悬浮的湍流脉动因引入颗粒而削弱,这是颗粒的减阻作用;颗粒增阻的一个主要机制是,流体给予颗粒的水平动量在颗粒一壁面碰撞中不断地损失．用悬浮动概念解释颗粒引起的增阻是不正确的．     

The Local Analysis of Changing Force Balances in Immiscible Incompressible Two-Phase Flow

The balance of viscous, capillary and gravity forces strongly affects two-phase flow through porous media and can therefore influence the choice of appropriate methods for numerical simulation and upscaling. A strict separation of the effects of these various forces is not possible due to the nature of the nonlinear coupling between the various terms in the transport equations. However, approximate prediction of this force balance is often made by calculation of dimensionless quantities such as capillary and gravity numbers. We present an improved method for the numerical analysis of simulations which recognises the changing balance of forces – in both space and time – in a given domain. The classical two-phase transport equations for immiscible incompressible flow are expressed in two forms: (i) the convection–diffusion-gravity (CDG) formulation where convection and diffusion represent viscous and capillary effects, respectively, (ii) the oil pressure formulation where the viscous effects are attributed to the product of mobility difference and the oil pressure gradient. Each formulation provides a different perspective on the balance of forces although the two forms are equivalent. By discretising the different formulations, the effect of each force on the rate of change of water saturation can be calculated for each cell, and this can be analysed visually using a ternary force diagram. The methods have been applied to several simple models, and the results are presented here. When model parameters are varied to determine sensitivity of the estimators for the balance of forces, the CDG formulation agrees qualitatively with what is expected from physical intuition. However, the oil pressure formulation is dominated by the steady-state solution and cannot be used accurately. In addition to providing a physical method of visualising the relative magnitudes of the viscous, gravity and capillary forces, the local force balance may be used to guide our choice of upscaling method.     

Analysis of reflective and non-reflective boundary conditions in LBM simulations of flow past an obstacle in a channel

This paper analyses the effects of considering non-reflective (characteristic boundary conditions) versus reflective boundary conditions (BCs) on the flow past a square object. We observe a clear dependence of the force exerted over the obstacle on the choice of BCs. Recirculation lengths, lift and drag coefficients, phase diagrams and streamlines for several angles of incidence are compared for a range of low Reynolds numbers (50–150) and two different values of the ratio of the object cross section to the channel width, 1/8 and 1/16. We remark distinct effects depending on the combination of BCs used at the inlet and at the outlet.     

Numerical simulations of drop size evolution in a horizontal pipeline

A population balance model using a standard method of moments ($S-\gamma$) in an Eulerian–Eulerian framework has been used for oil and brine two-phase flow simulations in pipelines. Results have been compared to both numerical and experimental data from the literature. The effects of the forces constituting the momentum transfer term at the interphase between droplets and the continuous phase (drag, lift, turbulent dispersion and virtual mass), turbulence modelling, break-up and coalescence parameters are analysed; they are shown to be important for droplet mean diameter evolution. It has been demonstrated that a correct combination of models and parameters improves (47% for the best case) simulated results when compared to experimental data. Interactions between the different components of the whole model are discussed and their corresponding effects on the droplet diameter predictions are explained. In particular, the addition of the lift force tends to push the droplet toward the walls of the computational domain where turbulence and shear stress are the strongest, therefore leading to an increased break-up rate. Based on the findings of this study, recommendations for further population balance-based modelling with a standard method of moments are provided.