均匀各向同性湍流中液滴团聚的实验研究

湍流与惯性颗粒相互作用在云物理、沙尘暴、气力输送和喷雾燃烧等地球物理现象和工业应用中广泛存在。由于湍流掺混效应和颗粒惯性的联合作用,颗粒在湍流中呈现非均匀的空间分布,这种现象被称为颗粒倾向聚集现象~([1][2]);具体表现为在颗粒空间分布上形成局部浓度高于平均值的团聚区和低于平均值的空隙区。此现象对众多工程应用均造成显著影响;例如,在聚集区中油滴的蒸发速率会降低,从而导致局部燃油富集区的出现,由此造成颗粒物排放的增加和燃烧效率的降低。因此,了解液滴聚集的机理对于优化众多的工程实际应用至关重要。颗粒-湍流相互作用是一个多尺度、非稳态的复杂物理过程,研究颗粒-湍流之间的关联耦合效应以及颗粒团聚的物理机制尤为困难。均匀各向同性湍流(Homogeneous isotropic turbulence,HIT)是一种对复杂湍流的较为理想的简化~([3]),由于没有平均剪切流动贡献,特别适合于研究多相流动。如图1所示,是中科院力学所的均匀各向同性湍流实验装置,其中HIT由八个扬声器的薄膜振动产生的空气射流形成;HIT区域在图中框架中心约5cm3的范围内。通过调节各个扬声器的幅值可以调节射流强度,从而改变HIT的泰勒雷诺数(Re_λ)。液滴由位于框架上方的喷雾产生,喷头位于测量区域(框架中心)正上方1m;通过调节喷雾的空气流量可以调节液滴粒径大小,调节喷雾的液体流量可以调节测量区域液滴的体积分数。在对流场进行标定后,本实验先采用高分辨率(4400×4400 pixel~2)CCD相机测量无液滴的各向同性湍流二维流场,采样频率为1Hz,测量区域为53 mm (x_1)×53 mm(x_2),其中x_1为垂直重力方向,x_2为平行重力方向;相应的速度分别为μ_1和μ_2,如图1所示;采用基于多层迭代变形窗口的互相关算法计算气相速度场,最后一层的查询窗口尺寸为32×32pixel~2,重叠率为75%;向量场空间分辨率为0.38 mm/vector;向量间距约为0.1mm或0.5η。如表1所示为两个泰勒雷诺数情况下各向同性湍流的基本参数。其后,用同样的实验配置测量液滴在HIT中的运动。在本实验中测试了三种不同的空气流量,并保持液体流量不变,液滴的体积分数约为1×10~(-5),可以认为液滴对流场基本没有影响。为了消除传统PTV算法计算液滴速度时出现的像素锁定现象,本实验中液滴速度采用PTV互相关算法获得液滴的离散速度场。作为对照实验,也测量了没有湍流(即未开扬声器)时,相同喷雾条件(液体和空气流量)下测量区域内液滴的运动和分布状态。最后用粒度仪测量了待测区域内液滴的粒径分布,并计算得到液滴的斯托克斯数St_p,发现液滴粒径表现出明显的多分散性。通过对比存在和不存在各向同性湍流时液滴的运动和空间分布,发现:在湍流作用下,液滴出现了明显的团聚作用,且团聚强度随Re_λ的增加而增加。采用Vorono?分析提取出液滴空间分布中的团聚区域和空隙区域,如图2所示。进一步的分析发现,团聚区中的液滴的沉降速度(μ_2)明显高于空隙区的液滴,后者更接近于流场速度;团聚区中的液滴沉降速度随浓度增加而增加,而空隙区中的沉降速度与浓度只呈弱相关。产生这种现象的原因是惯性液滴的局部扫掠现象,即在湍流的作用下,液滴倾向于朝重力方向运动,继而形成沿重力方向的团聚。另外,对液滴团聚区进行了几何特性的分析,发现液滴团聚区存在自相似特性。对单个团聚区中的液滴空间坐标进行奇异值分解,从而得到单个团聚区域的主轴长度、长细比、倾角等;发现团聚尺度越大,团聚区的长细比越大,且越倾向于重力方向。基于Vorono?分析,将液滴的速度信息映射到相应的Vorono?多边形内,可以对液滴的速度场进行两点相关分析,同时对湍流场也进行此分析,发现液滴的横向速度分量μ_1的相干结构的横向和纵向尺度同步减少,其尺度比例与湍流的μ_1分量基本一致,说明液滴μ_1分量遵循湍流的尺度相关性;而对于液滴μ_2分量的相干结构,相比流场,横向、纵向尺度均有增加,且横向尺度增加更多,这有可能是重力对液滴团聚结构的产生了显著的影响

求解运动方程的一种不等时间步长的显式数值积分方法

在波动有限元模拟中, 若采用传统的显式数值积分方法求解运动方程, 计算时间步长需采用计算区内满足稳定条件要求的最小时间步长. 然而, 对于大部分计算区域, 这一时间步长过小, 是不必要的. 本文提出了一种不等时间步长的显式数值积分方法, 其基本思想是不同的计算区域采用满足各自稳定条件的计算时间步长. 最后, 本文通过数值试验检验了这一方法的可行性及其对数值计算精度的影响

A varying time-step explicit numerical integration algorithm for solving motion equation

If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simu- lation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stabil- ity criterion in computational region. However, the excessivel..

Rayleigh波入射情形下的断层效应

本文运用数值模拟方法对Rayleigh波入射情形下的断层效应进行分析,结果表明:断层对Rayleigh面波具有隔震作用,随断层深度的增加其隔震效果更加明显,且水平向的隔震效果强于竖向.在入射波一侧,由于断层的作用使接触带附近的反应增强,随着断层破碎带宽度的增加,竖向反应变化特征明显不同于水平向反应变化特征

Research and Development on the Microbial Furfural-Tolerance Mechanism

由木糖脱水生成的糠醛会严重抑制微生物的生长与后续发酵过程。高温液态水和稀酸处理等是去除生物质中半纤维素较好的方法,却不可避免地会导致大量糠醛产生。本文简述糠醛的产生和降解途径及糠醛对微生物危害的研究进展,对糠醛耐受过程中涉及到的细胞形态分析与氧化胁迫、蛋白响应、基因响应和氧化还原平衡等进行了概述

中国球状绿藻新纪录属—索囊藻属(Choricystis)

在西藏那曲的错鄂和湖北武汉的紫阳湖采集到两株淡水绿色海绵。形态观察显示,两株海绵呈现绿色是由于其内共生有球状绿藻。基于形态观察和分子系统发育分析确定为该球状绿藻为索囊藻属的寄生索囊藻(Choricystis parasitica(K.Brandt)Pr?schold&Darienko)。结果表明:(1)两株藻均为纺锤形或倒卵形,长1.2—2.8μm,宽0.8—1.5μm;具一个侧位片状色素体;(2)依据18S r DNA和rbc L cp DNA的进化分析显示这两株藻均位于胶球藻科(Coccomyxaceae),索囊藻属(Choricystis);(3)基于ITS2的二级结构分析显示这两株藻与寄生索囊藻具有相同的二级结构,表明这两株内共生的球状绿藻均为寄生索囊藻。寄生索囊藻为中国新纪录种,所在属索囊藻属为中国新纪录属。本研究以形态性状和分子性状描述了采自中国的寄生索囊藻标本,为今后研究提供依据和参考

激光毛化形貌对高速轮轨冰润滑黏着系数的影响

轮轨黏着系数对于列车运行的安全性与稳定性具有决定性的影响。在冬季,钢轮表面的薄冰会大幅降低轮轨黏着系数。根据水的相图和压力融化原理,将轮轨接触区域划分为钢-冰固体接触区和钢-水混合润滑接触区,基于统一雷诺方程模型,研究了温度在-21℃至0℃,运动速度为100～500 km/h范围内,高速轮轨黏着系数随速度的变化规律,并提出采用激光毛化技术在车轮表面形成规则形貌(纵纹、横纹、菱形),将有助于提高冰水混合润滑状态下的轮轨黏着系数,为高速列车黏着控制技术提供理论基础

中国球状绿藻新纪录属—索囊藻属(Choricystis)

在西藏那曲的错鄂和湖北武汉的紫阳湖采集到两株淡水绿色海绵。形态观察显示,两株海绵呈现绿色是由于其内共生有球状绿藻。基于形态观察和分子系统发育分析确定为该球状绿藻为索囊藻属的寄生索囊藻(Choricystis parasitica(K.Brandt)Pr?schold&Darienko)。结果表明:(1)两株藻均为纺锤形或倒卵形,长1.2—2.8μm,宽0.8—1.5μm;具一个侧位片状色素体;(2)依据18S r DNA和rbc L cp DNA的进化分析显示这两株藻均位于胶球藻科(Coccomyxaceae),索囊藻属(Choricystis);(3)基于ITS2的二级结构分析显示这两株藻与寄生索囊藻具有相同的二级结构,表明这两株内共生的球状绿藻均为寄生索囊藻。寄生索囊藻为中国新纪录种,所在属索囊藻属为中国新纪录属。本研究以形态性状和分子性状描述了采自中国的寄生索囊藻标本,为今后研究提供依据和参考