惯性颗粒所见被动标量统计特性的直接数值模拟研究

本文采用直接数值模拟方法,对惯性颗粒所见均匀各向同性湍流中具有平均标量梯度的被动标量场统计特性进行了研究。模拟结果表明:与颗粒所见流体速度的自相关特性不同,颗粒所见标量的自相关特性随颗粒惯性的增加而单调减少;颗粒所见标量脉动能随颗粒惯性的增加先减少再增大,在St≈1.0的临界颗粒附近达到最小值,而颗粒所见标量耗散率随颗粒惯性的变化行为与颗粒所见标量脉动能的变化行为相反;数值模拟的结果进一步揭示,在St≈1.0时颗粒所见流体标量脉动能和耗散率的极值是因为St≈1.0的临界颗粒聚集于低涡量、高应变区域和标量场在高应变区域形成强耗散的片状结构所致。     

非等温气固两相各向同性湍流的直接数值模拟

本文对非等温气固两相各向同性湍流中颗粒及颗粒所见流体温度的统计行为进行了直接数值模拟研究,系统地讨论了颗粒惯性的影响。对于Tp/Tk<1的颗粒,其温度脉动随惯性增大而减弱;对于Tp/Tk>1的颗粒,其趋势相反。颗粒及颗粒所见流体温度的自相关都随颗粒惯性的增加而减弱,而且随颗粒惯性的增加,颗粒所见流体温度的自相关比颗粒自身温度的自关联下降得快。     

单相泵驱流体回路不同拓扑结构的散热性能分析

基于单相泵驱流体回路,研究了系统在不同拓扑结构下的控温性能。采用集总参数法建立了流量计算及分配模型、流体与设备换热模型、混合点温度计算模型以及热流计算模型,并运用Openmodelica仿真平台,开展了不同拓扑结构对设备温度、系统流量和系统出口工质温度的影响分析,探究了不同载荷占空下设备散热效果、系统峰值损耗和内能变化等的动态特性。结果表明,一方面,并联拓扑体系设备散热及均温效果较好;另一方面,热惯性对系统内能变化影响较大,在系统热惯性一定的条件下,载荷占空比越大工质冷却响应时间越大,系统能量损失越小。     

高雷诺数下非混相Rayleigh-Taylor不稳定性的格子Boltzmann方法模拟

本文采用相场格子Boltzmann方法研究了竖直微通道内中等Atwoods数流体的单模Rayleigh-Taylor不稳定性问题,系统分析了雷诺数对相界面动力学行为以及扰动在各发展阶段演化规律的影响.数值结果表明高雷诺数条件下,不稳定性界面扰动的增长经历了四个不同的发展阶段,包括线性增长阶段、饱和速度阶段、重加速阶段及混沌混合阶段.在线性增长阶段,我们计算获得的气泡与尖钉振幅符合线性稳定性理论,并且线性增长率随着雷诺数的增加而增大.在第二个阶段,我们观察到气泡与尖钉将以恒定的速度增长,获得的尖钉饱和速度略高于Goncharov经典势能模型的解析解[Phys.Rev.Lett.200288134502],这归因于系统中产生了多个尺度的旋涡,而涡之间的相互作用促进了尖钉的增长.随着横向速度和纵向速度的差异扩大,气泡和尖钉界面演化诱导产生的Kelvin–Helmholtz不稳定性逐渐增强,从而流体混合区域出现许多不同层次的涡结构,加速了气泡与尖钉振幅的演化速度,并在演化后期阶段,导致界面发生多层次卷起、剧烈变形、混沌破裂等行为,最终形成了非常复杂的拓扑结构.此外,我们还统计了演化后期气泡与尖钉的无量纲加速度,发现气泡和尖钉的振幅在后期呈现二次增长规律,其增长率系数分别为0.045与0.233.而在低雷诺条件下,重流体在不稳定性后期以尖钉的形式向下运动而轻流体以气泡的形式向上升起.在整个演化过程中,界面变得足够光滑,气泡与尖钉在后期的演化速度接近于常数,未观察到后期的重加速与混沌混合阶段.     

微通道内椭球颗粒惯性聚焦行为数值研究

基于浸没边界-格子Boltzmann方法，对方形截面微通道内椭球颗粒的惯性迁移与旋转动力学行为进行数值研究，发现微通道内椭球颗粒的惯性迁移存在两种主要的运动状态：①翻转状态，即椭球颗粒前进过程中长轴始终在中心对称平面内；②滚动状态，即椭球颗粒前进过程中长轴始终垂直于中心对称平面.研究表明：在低Re数（Re=10）下颗粒以两种状态随流体迁移至平衡位置；在较大Re数（50≤Re≤200）下最终均以翻转状态随流体迁移，随Re数增加，平衡位置先逼近壁面后远离壁面.通过对不同运动状态下椭球颗粒周围的微观流场进行分析，提示该微观流动在颗粒惯性聚焦行为特征中有重要影响，并从流体和颗粒的惯性角度对颗粒不同运动状态的转换机理给出解释.     

耦合界面力的两相流相场格子Boltzmann模型

提出了一种改进的基于相场理论的两相流格子Boltzmann模型.通过引入一种新的更加简化的外力项分布函数,使得此模型克服了前人工作中界面力尺度与理论分析不一致的问题,并且通过Chapman-Enskog多尺度分析表明,所提出的模型能够准确恢复到追踪界面的Cahn-Hilliard方程和不可压的Navier-Stokes方程,并且宏观速度的计算更为简化.利用所提模型对几个经典两相流问题,包括静态液滴测试、液滴合并问题、亚稳态分解以及瑞利-泰勒不稳定性进行了数值模拟,发现本模型可以获得量级为10^-9极小的虚假速度,并且这些算例获取的数值解与解析解或已有的文献结果相吻合,从而验证了模型的准确性和可行性.最后,利用所发展的两相流格子Boltzmann模型研究了随机扰动的瑞利-泰勒不稳定性问题,并着重分析了雷诺数对流体相界面的影响.发现对于高雷诺数情形,在演化前期,流体界面出现一排“蘑菇”形状,而在演化后期,流体界面呈现十分复杂的混沌拓扑结构.不同于高雷诺数情形,低雷诺数时流体界面变得相对光滑,在演化后期未观察到混沌拓扑结构.     

纳米尺度下流体剪切流动的分子动力学模拟

利用分子动力学方法模拟流体在两无限大平板间剪切流动过程。研究通道内加入不同体积分数纳米颗粒、体积分数相同纳米颗粒数目不同以及剪切速度对流体密度、速度以及界面滑移的影响。结果表明:近壁区流体数密度呈衰减振荡分布,由近壁区到主流区振幅逐渐减小,颗粒和流体的整体数密度在中心主流区呈抛物线分布。流体剪切应变率随颗粒体积分数的增加逐渐减小,同时剪切应变率和滑移速度之间呈近似线性分布。体积分数相同颗粒数目不同,颗粒在运动过程中呈线,性排列时,剪切应变率最大。随剪切速度的增加,流体滑移速度和滑移长度随之增大,但滑移长度增加量相对较小。     

惯性颗粒所见被动标量统计特性的DNS与模型研究

本文采用DNS方法,对惯性颗粒所见各向同性湍流中具有平均标量梯度的被动标量场统计特性进行了研究。结果表明：惯性对颗粒温度脉动强度,两相温度关联,自相关特性以及颗粒热流与两相交叉热流的统计特性具有明显的影响。在PDF方法的框架下,系统地推导了非等温气固两相流的PDF方程,且基于朗之万随机体系对方程进行了封闭,并利用前面的...     

各向同性湍流中颗粒引起湍流变动的直接模拟

本文通过直接数值模拟对均匀各向同性湍流中颗粒对湍流的变动作用进行了研究.颗粒相的体积分数很小而质量载荷足够大,以至于颗粒之间的相互作用可以忽略不计,而重点考虑颗粒与湍流间能量的交换。颗粒对湍流的反向作用使得湍动能的耗散率增强,以至于湍动能的衰减速率增大.湍动能的衰减速率随颗粒惯性的增大而增大。三维湍动能谱显示,颗粒对湍动能的影响在不同的尺度上是不均匀的。在低波数段,流体带动颗粒,而高波数段则相反.     

激光脉冲能量对激光诱导铝合金等离子体物理特性的影响

为研究激光脉冲能量对激光诱导等离子体辐射特性和膨胀过程的影响,采用ICCD相机对不同激光脉冲能量激发的铝合金等离子体进行快速成像,并利用Boltzmann斜线法和Stark展宽法分析等离子体的电子温度和电子数密度随激光脉冲能量的演化规律.实验结果表明,激光诱导等离子体呈现明显的分层结构,等离子体的激发阈值约为3mJ,等离子体不同区域的面积随激光脉冲能量变化呈现不同的特征.当激光脉冲能量低于10mJ时,等离子体的分层结构不显著.激光脉冲能量从10mJ增加到100mJ过程中,等离子体电子温度从4 980K升高到7 221K,等离子体的电子数密度在1017 cm-3量级并随激光能量增加而增大且趋于饱和.     

Lagrangian view of time irreversibility of fluid turbulence

A turbulent flow is maintained by an external supply of kinetic energy, which is eventually dissipated into heat at steep velocity gradients. The scale at which energy is supplied greatly differs from the scale at which energy is dissipated, the more so as the turbulent intensity(the Reynolds number) is larger. The resulting energy flux over the range of scales, intermediate between energy injection and dissipation, acts as a source of time irreversibility. As it is now possible to follow accurately fluid particles in a turbulent flow field, both from laboratory experiments and from numerical simulations, a natural question arises: how do we detect time irreversibility from these Lagrangian data? Here we discuss recent results concerning this problem. For Lagrangian statistics involving more than one fluid particle, the distance between fluid particles introduces an intrinsic length scale into the problem. The evolution of quantities dependent on the relative motion between these fluid particles, including the kinetic energy in the relative motion, or the configuration of an initially isotropic structure can be related to the equal-time correlation functions of the velocity field, and is therefore sensitive to the energy flux through scales, hence to the irreversibility of the flow. In contrast, for singleparticle Lagrangian statistics, the most often studied velocity structure functions cannot distinguish the "arrow of time". Recent observations from experimental and numerical simulation data, however, show that the change of kinetic energy following the particle motion, is sensitive to time-reversal. We end the survey with a brief discussion of the implication of this line of work.     

Statistics of particle pair relative velocity in the homogeneous shear flow

Small scale clustering of inertial particles and relative velocity of particle pairs have been fully characterized for statistically steady homogeneous isotropic flows. Depending on the particle Stokes relaxation time, the spatial distribution of the disperse phase results in a multi-scale manifold characterized by local particle concentration and voids and, because of finite inertia, the two nearby particles have high probability to exhibit large relative velocities. Both effects might explain the speed-up of particle collision rate in turbulent flows. Recently it has been shown that the large scale geometry of the flow plays a crucial role in organizing small scale particle clusters. For instance, a mean shear preferentially orients particle patterns. In this case, depending on the Stokes time, anisotropic clustering may occur even in the inertial range of scales where the turbulent fluctuations which drive the particles have already recovered isotropy. Here we consider the statistics of particle pair relative velocity in the homogeneous shear flow, the prototypical flow which manifests anisotropic clustering at small scales. We show that the mean shear, by imprinting anisotropy on the large scale velocity fluctuations, dramatically affects the particle relative velocity distribution even in the range of small scales where the anisotropic mechanisms of turbulent kinetic energy production are sub-dominant with respect to the inertial energy transfer which drives the carrier fluid velocity towards isotropy. We find that the particles’ populations which manifest strong anisotropy in their relative velocities are the same which exhibit small scale clustering. In contrast to any Kolmogorov-like picture of turbulent transport these phenomena may persist even below the smallest dissipative scales where the residual level of anisotropy may eventually blow-up. The observed anisotropy of particle relative velocity and spatial configuration is suggested to influence the directionality of the collision probability, as inferred on the basis of the so-called “ghost collision” model.     

Large eddy simulation of orientation and rotation of ellipsoidal particles in isotropic turbulent flows

The rotational motion and orientational distribution of ellipsoidal particles in turbulent flows are of significance in environmental and engineering applications. Whereas the translational motion of an ellipsoidal particle is controlled by the turbulent motions at large scales, its rotational motion is determined by the fluid velocity gradient tensor at small scales, which raises a challenge when predicting the rotational dispersion of ellipsoidal particles using large eddy simulation (LES) method due to the lack of subgrid scale (SGS) fluid motions. We report the effects of the SGS fluid motions on the orientational and rotational statistics, such as the alignment between the long axis of ellipsoidal particles and the vorticity, the mean rotational energy at various aspect ratios against those obtained with direct numerical simulation (DNS) and filtered DNS. The performances of a stochastic differential equation (SDE) model for the SGS velocity gradient seen by the particles and the approximate deconvolution method (ADM) for LES are investigated. It is found that the missing SGS fluid motions in LES flow fields have significant effects on the rotational statistics of ellipsoidal particles. Alignment between the particles and the vorticity is weakened; and the rotational energy of the particles is reduced in LES. The SGS-SDE model leads to a large error in predicting the alignment between the particles and the vorticity and over-predicts the rotational energy of rod-like particles. The ADM significantly improves the rotational energy prediction of particles in LES.     

Velocity distributions of granular gases with drag and with long-range interactions

We study velocity statistics of electrostatically driven granular gases. For two different experiments, (i) nonmagnetic particles in a viscous fluid and (ii) magnetic particles in air, the velocity distribution is non-Maxwellian, and its high-energy tail is exponential, P(upsilon) approximately exp(-/upsilon/). This behavior is consistent with the kinetic theory of driven dissipative particles. For particles immersed in a fluid, viscous damping is responsible for the exponential tail, while for magnetic particles, long-range interactions cause the exponential tail. We conclude that velocity statistics of dissipative gases are sensitive to the fluid environment and to the form of the particle interaction.     

Settling statistics of hard sphere particles

Direct imaging of settling, non-Brownian, hard sphere, particles allows measurement of particle occupancy statistics as a function of time and sampling volume dimension. Initially random relative particle number fluctuations, (2)>/ = 1, become suppressed, anisotropic, and dependent. Fitting to a simple Gaussian pair correlation model suggests a minute long ranged correlation leads to strong if not complete suppression of number fluctuations. Calflisch and Luke predict a divergence in velocity fluctuations with increasing sample volume size based on random (Poisson) statistics. Our results suggest this is not a valid assumption for settling particles.     

Velocity statistics in excited granular media

We present an experimental study of velocity statistics for a partial layer of inelastic colliding beads driven by a vertically oscillating boundary. Over a wide range of parameters (accelerations 3-8 times the gravitational acceleration), the probability distribution P(v) deviates measurably from a Gaussian for the two horizontal velocity components. It can be described by P(v) approximately exp(-mid R:v/v(c)mid R:(1.5)), in agreement with a recent theory. The characteristic velocity v(c) is proportional to the peak velocity of the boundary. The granular temperature, defined as the mean square particle velocity, varies with particle density and exhibits a maximum at intermediate densities. On the other hand, for free cooling in the absence of excitation, we find an exponential velocity distribution. Finally, we examine the sharing of energy between particles of different mass. The more massive particles are found to have greater kinetic energy. (c) 1999 American Institute of Physics.     

Energy fluctuation and correlation in Tsallis statistics

We investigate the fluctuation of the energy in the framework of Tsallis statistics and find the correlation plays an important role in energy fluctuations. In Tsallis statistics, the correlation is induced by the nonextensivity of Tsallis entropy and exists between particles even if the particles are dynamically independent. By taking the generalized ideal gas as an example, we get that when the particle number N is large enough, the relative fluctuation of the energy is proportional to 1/N instead of in Boltzmann statistics. Thus, the relative energy fluctuation is much smaller in Tsallis statistics than that in Boltzmann statistics. Besides, we demonstrate that the introduction of correlation between particle energies leads to smaller energy fluctuations in Tsallis statistics.     

Finite-size particles in turbulent channel flow: quadrant analysis and acceleration statistics

The interaction between finite-size particles and turbulent channel flow in the absence of gravity is studied by direct numerical simulations (DNS). The study is motivated by DNS observations of a turbulent channel flow with high-density, pointwise particles, that cluster in regions of high streamwise root-mean-square (RMS) acceleration close to the wall, contrary to what is observed in homogeneous isotropic turbulence. The aim of the present study is to explore if this is still the case when size effects are taken into account in the DNS. Based on the analysis of the velocity and acceleration statistics, the present DNS shows that, close to the wall, particles with ρ_p/ρ_f ranging from 2 to 4 are surrounded by regions with low streamwise RMS velocity but high streamwise RMS acceleration. According to the normalised particle acceleration probability density functions (PDFs), size effects become important in the near-wall region. As particle inertia increases, the normalised PDFs of particle acceleration tend to a Gaussian distribution. The tails of the normalised PDFs of the fluid conditioned by the presence of particles are higher than that of the unconditioned fluid close to the wall.     

固液两相流中微对流强化的机理分析与数值模拟

本文对分散型固液两相混合物层流管流中非均匀剪切力场导致的微对流现象以及由此引起的导热系数增强效应作了机理和影响因素分析,认为除速度分布、颗粒浓度和粒径以外,还存在颗粒形状、粒径分布,壁面热流方向、颗粒表面性质及其与载流介质间的相容性等多个影响因素.模拟计算表明,由微对流导致的对流换热强化与流体在管壁面上的表观导热系数强化具有相同的数量。     

等单元长度多间隙加速结构的束流动力学特性

等单元长度多间隙加速结构是一种非同步加速结构,当粒子在加速结构中的速度变化很明显时,粒子在每个间隙的相位是不相同的,薄透镜近似下的束流纵向运动方程没有考虑粒子在加速结构中的速度变化,这在单腔的能量增益相对于粒子能量很小的情况下是合理的,但是当粒子在加速结构中的速度有明显变化时,这种处理方式是不够的。本文从带电粒子在电磁场中的运动方程出发,通过目前普遍采用的理论建立了束流在这种加速结构中的纵向运动方程,通过数值方法得到了粒子在这种结构的相运动,计算了不同初始能量的粒子在不同的电场梯度及单元数的加速腔中运动的能量变化情况。结果显示:粒子相运动与不考虑粒子速度变化的薄透镜近似下的相运动轨迹有明显差异,当粒子速度变化超过一定值后,薄透镜近似下束流纵向运动方程在计算束流能量增益时与同步相位之间满足的余弦关系不再成立,并给出了能量增益变化曲线相对余弦曲线的偏离情况与电场梯度、间隙数、初始能量等参数的关系。