气温和颗粒密度对声场中颗粒动力学影响的数值模拟

综合考虑黏性夹带力、Basset力、虚拟质量力和压力梯度力,建立颗粒在声场中的动力学模型,利用变步长四阶RungeKutta算法和二阶隐式Adams插值算法对颗粒的受力和运动进行数值模拟。将模拟和实验得到的颗粒运动特性进行对比,验证数值模拟的正确性。在此基础上,研究气温和颗粒密度对颗粒动力学的影响规律。结果表明,黏性夹带力对颗粒运动起主导作用;气温升高,压力梯度力与黏性夹带力之间的相位差减小,Basset力、虚拟质量力与黏性夹带力之间的相位差增大。研究还发现,气温较低时,颗粒密度对颗粒运动有重要影响,夹带系数随着密度的增加而迅速下降;气温较高时,颗粒密度对颗粒运动的影响较小,颗粒位移振幅和夹带系数相对低温时明显增加。      

Photophoresis—Light induced motion of particles suspended in gas

When irradiated sideways, by visible light, a particle can perform different kinds of motion, (e.g. in direction of irradiation, opposite to irradiation, vertical movement, helicoidally, etc.). This phenomenon is called photophoresis. Photophoresis is based on momentum transfer between the aerosol particle and surrounding gas molecules. Photophoresis strongly depends on the pressure of the surrounding gas. Particles mostly influenced by photophoresis are those of μm size. Two main types of forces describe photophoretic motion: ΔT force: The thermal accommodation coefficient α is constant over the particle surface. As a result of the thermal accommodation, gas molecules on the warm side of a particle leave the surface faster than gas molecules on the cold side. This leads to ΔT force on the particle towards the colder side. Typical motion of the particle will be either away from light irradiation (positive photophoresis), or in direction of light irradiation (negative photophoresis). In the case of negative photophoresis, the back side of the particle, due the nature of light absorption, will be heated more than front side of the particle. Δα force: If the particle is at a constant temperature, which is different from the temperature of the surrounding gas, and the thermal accommodation coefficient α varies over the particle surface, the net momentum between gas molecules and particle will be transferred. In this case, the result will be body fixed Δα force. Depending on the particle surface properties, Δα force can direct the particle in any possible photophoresis could also play important role in planet formation and astrophysics.     

自由分子区内纳米颗粒的热泳力计算

基于非平衡态分子动力学模拟方法,研究了自由分子区内纳米颗粒的热泳特性.理论研究表明,纳米颗粒与周围气体分子之间的非刚体碰撞效应会明显地改变其热泳特性,经典的Waldmann热泳理论并不适用,但尚未有定量的直接验证.模拟计算结果表明:对于纳米颗粒而言,当气-固相互作用势能较弱或气体温度较高时,气体分子与纳米颗粒之间的非刚体碰撞效应可以忽略,Waldmann热泳理论与分子动力学模拟结果吻合较好;当气-固相互作用势能较强或气体温度较低时,非刚体碰撞效应较为明显,Waldmann热泳理论与模拟结果存在较大误差.基于分子动力学模拟结果,对纳米颗粒的等效粒径进行了修正,并考虑了气体分子与纳米颗粒之间的非刚体碰撞效应,理论计算结果与分子动力学模拟结果吻合较好.     

Mathematical combustion model of nanoaluminum-air suspension

This paper presents a combustion model of a nano-aluminum-air (nAl-air) suspension. The special feature of the model is performing a local mathematical model of the oxidant diffusion through an aluminum oxide layer on the particle surface taking into account the aluminum-oxidant reaction to simulate the combustion of nano-size aluminum (nAl) particles. The oxidation rate of the aluminum particles and the associated with this process the rate of heat release are determined from the solution of the local combustion problems for the entire set of nAl particles in the suspension. To obtain the suspension state parameters we solve the equation system, which includes the energy conservation equations for the gas and particles, the mass-conservation equation for the gas-dispersed mixture and the motion equations for the gas and particles controlling for the particle velocity lag. The model considers gas expansion and thus gas and particle motion. The developed model does not require setting the ignition temperature of nAl particles. The study provides the calculated propagation rate of the combustion front in the nAl-air suspension depending on the nAl mass concentration and on the initial temperature of the suspension.     

Path integral approach to quantum diffusion in metals

A general formalism for a quantum particle moving in an environment based on the path integral formulation of quantum mechanics is presented. The environment can be the valence electron gas of a metal or the phonons of the host in which the particle moves, or both. The advantage of this approach is that it gives a direct space-time picture of the motion. Integrating out the coordinates of the environment an effective action for the particle is obtained. It consist of two parts. An adiabatic term, that is local in time (and as such describable as a potential), describing the influence of the fully relaxed environment on the particle, and a term that accounts for the dynamical part of the response of the environment due to the motion of the particle. As an example I have considered a particle moving in a tight binding band with an electron gas environment. Three energy scales play a role here. They are the bare hopping amplitude, the temperature and the inverse tunneling time (i.e. the time it takes the particle to tunnel from one site to the next). Note that the Fermi energy drops out of the problem. I find, as did Kondo in his recent theory, a decreasing diffusion constant as a function of temperature.     

Drift Velocity Induced by Collisions

Stochastic motion of a hard point separating two semi-infinite subvolumes of a hard point gas in R ¹ is studied. The partitionning particle is identical to the particles of the gas and can be looked upon as a tagged particle playing the role of a (microscopic) piston. At the initial moment it is at rest having to the right and to the left of it gases in thermodynamic equilibrium. Its further motion is entirely induced by collisions. The stochastic motion of the piston is determined rigorously. The form of the stationary velocity distribution is calculated. It turns out that at equal initial pressures the piston acquires asymptotically a drift velocity oriented towards the higher temperature region. There is no drift if the temperatures and densities combine to produce on both sides equal particle fluxes. Although the qualitative agreement with Boltzmann's theory is found, the Boltzmann equation does not predict correctly the thermodynamic conditions under which the drift vanishes.     

A kinetic theory of diffusely reflecting Brownian particles

An analysis is made of the motion of a spherical Brownian particle whose surface can diffusely reflect the molecules of an equilibrium host gas. The analysis is based on Newton's second law and a limiting form of Markov's method. It is shown, both for specular and diffuse reflections, that equipartition of energy is a consequence of the dynamics and randomness of the motion. In addition, it is demonstrated that the diffusion coefficient can depend on the temperature of the particle. The entire analysis is restricted to the case for which the Knudsen number of the particle is large compared to unity.Slinn's work was supported in part by Battelle Memorial Institute and in part by the U.S. Atomic Energy Commission contract AT(45–1)-1830. Shen's work was supported in part by the U.S. Air Force Office of Scientific Research contract 49(638)–1346. Mazo's work was supported in part by the National Science Foundation, NSF Grant GP–8497.     

三维混合沙输运数值模拟

采用计算流体力学和颗粒离散元耦合的方法模拟了三维混合沙输运过程。采用体平均的Navier-stokes方程来描述气相运动,考虑了气相和颗粒相的相互作用。颗粒运动通过求解牛顿运动方程来求解,采用硬球模型描述颗粒和颗粒及颗粒和壁面的碰撞。本模型中,颗粒运动是三维的而气相运动是二维的。计算结果表明:总输沙率沿高度方向在大于2cm以上按照指数衰减,在2 cm以下则偏大;各粒径颗粒具有不同的输沙率分布,粗粒径颗粒按指数规律衰减,其它粒径颗粒输沙率随高度先指数增加后减少;各粒径颗粒平均水平速度随高度对数函数增加且同高度时随粒径增大而减小,1 cm高度以下则相反;沙粒平均粒径沿高度线性递减,2 cm以下粒径偏大。     

Particle motion in a photon gas: friction matters

The motion of a particle in the Tolman metric generated by a photon gas source is discussed. Both the case of geodesic motion and motion with nonzero friction, due to photon scattering effects, are analyzed. In the Minkowski limit, the particle moves along a straight line segment with a decelerated motion, reaching the endpoint at zero speed. The curved case shows a qualitatively different behavior; the geodesic motion consists of periodic orbits, confined within a specific radial interval. Under the effect of frictional drag, this radial interval closes up in time and in all our numerical simulations the particle ends up in the singularity at the center.     

一种新的光子多普勒速度频谱分析方法

光子多普勒速度计可给出飞层表面某一速度带内颗粒群速度随时间演化的频谱数据, 在冲击动力学实验尤其是微喷射及其混合研究中得到广泛应用. 本文提出一种新的光子多普勒频谱数据分析方法, 可推断出混合区厚度变化和前端等效颗粒尺度. 利用该方法, 对一些典型状态下喷射混合速度频谱开展分析, 获得了不同冲击压力、气体条件下颗粒度数据, 证实了气体环境下喷射颗粒的气动破碎现象, 以及破碎后尺度与初始条件的依赖性, 为喷射混合物理规律研究提供了重要依据.     

Effect of the particle temperature on lift force of nanoparticle in a shear rarefied flow

The nanoparticles suspended in a shear flow are subjected to a shear lift force, which is of great importance for the nanoparticle transport. In previous theoretical analysis on the shear lift, it is usually assumed that the particle temperature is equal to the temperature of the surrounding gas media. However, in some particular applications, the particle temperature can significantly differ from the gas temperature. In the present study, the effect of particle temperature on the shear lift of nanoparticles is investigated and the corresponding formulas of shear lift force are derived based on the gas kinetic theory. For extremely small nanoparticles(with radius R 2 nm) or large nanoparticles(R 20 nm), the influence of the particle temperature can be neglected. For the intermediate particle size, the relative error induced by the equal gas–particle temperature can be significant. Our findings can bring an insight into accurate evaluation of the nanoparticle transport properties.     

A new thermal conductivity model for nanofluids

In a quiescent suspension, nanoparticles move randomly and thereby carry relatively large volumes of surrounding liquid with them. This micro-scale interaction may occur between hot and cold regions, resulting in a lower local temperature gradient for a given heat flux compared with the pure liquid case. Thus, as a result of Brownian motion, the effective thermal conductivity, k_eff, which is composed of the particles conventional static part and the Brownian motion part, increases to result in a lower temperature gradient for a given heat flux. To capture these transport phenomena, a new thermal conductivity model for nanofluids has been developed, which takes the effects of particle size, particle volume fraction and temperature dependence as well as properties of base liquid and particle phase into consideration by considering surrounding liquid traveling with randomly moving nanoparticles.The strong dependence of the effective thermal conductivity on temperature and material properties of both particle and carrier fluid was attributed to the long impact range of the interparticle potential, which influences the particle motion. In the new model, the impact of Brownian motion is more effective at higher temperatures, as also observed experimentally. Specifically, the new model was tested with simple thermal conduction cases, and demonstrated that for a given heat flux, the temperature gradient changes significantly due to a variable thermal conductivity which mainly depends on particle volume fraction, particle size, particle material and temperature. To improve the accuracy and versatility of the k_effmodel, more experimental data sets are needed.     

Brownian motion and correlation in particle image velocimetry

In particle image velocimetry applications involving either low velocities or small seed particles, Brownian motion can be significant. This paper addresses the effects of Brownian motion. First, general equations describing cross-correlation particle image velocimetry are derived that include Brownian motion. When light-sheet illumination particle image velocimetry (PIV) is used Brownian motion diminishes the signal strength. A parameter describing this effect is introduced, and a weighting function describing the contribution to the measured velocity as a function of position is derived. The latter is unaffected by Brownian motion. Microscopic PIV Brownian motion also diminishes the signal strength. The weighting function for microscopic PIV is found to depend on Brownian motion, thus affecting an important experimental parameter, the depth of correlation. For both light-sheet illumination and microscopic PIV, a major consequence of Brownian motion is the spreading of the correlation signal peak. Because the magnitude of the spreading is dependent on temperature, PIV can, in principle, be used to simultaneously measure velocity and temperature. The location of the signal peak provides the velocity data, while the spreading of the peak yields temperature.     

固体粒子在直流式除尘器中湍流脉动的数值模拟

本文对直流除尘器涡室内固体粒子的湍流脉动现象进行了数值分析,通过气体速度随机脉动谱,把气相湍流运动对固体粒子运动的影响引入粒子的运动平衡方程中,用拉格朗日法模拟了粒子的轨迹及其扩散运动,应用四阶龙科库塔方法求解粒子的运动方程。计算结果表明对粒子的数值模拟可以较好地预测除尘器的性能,如除尘器的切割粒径。     

循环流化床内颗粒团流动的研究

循环流化床颗粒相流动具有多尺度效应:单颗粒运动的微尺度、颗粒团运动的介尺度和固相运动的宏尺度。颗粒相流动参数受单颗粒运动和颗粒聚团运动的制约,同时影响气相流动。基于气体分子运动论和颗粒动理学,建立相平均稠密气固两相流流动模型。介尺度模型考虑颗粒团与单颗粒之间、颗粒团与气相之间的动量和能量的传递和耗散。模拟计算颗粒容积份额、速度等参数与实测值相吻合。     

Geschwindigkeitsverteilung,Teilchendichte, transversale Driftgeschwindigkeit und Energiedichte des Neutralgases in der Freifallsäule unter Berücksichtigung der „Neutralgasaufzehrung”︁ durch die Ionisationsvorgänge

A neutral gas rarefaction caused by ionization processes occurs in the plasma of the low pressure gas discharges. The velocity distributions, the particle density, the transversal drift velocity and the energy density of the neutral gas are calculated both for the plane and for the cylindrical positive column under free fall conditions. The neutral gas rarefaction is taken into account. It is shown, that the velocity distributions is non-Maxwellian and anisotropic. The pressure tensor is anisotropic, too. Particle density and energy density of neutral gas decrease with increasing electron density and electron temperature relatively homogeneously over the cross section of the column. Only, if the degree of ionization is high, these densities are much smaller near the axis than at the wall. Decreasing neutral gas temperature causes a similar change in the particle density profile as increasing electron density and electron temperature do. The transverse neutral gas pressure decreases from the axis to the wall in all cases. In the steady-state column an upper limit exists for the transverse particle current density of the neutral gas and of the ion gas. This limit depends on the gas temperature, the filling pressure and the atomic mass of the filling gas. In the appendix the Boltzmann equation is given in a form, which is suitable to investigate cylindrical problems not only for simple examples.     

Measurement of particle velocity and characterization of deposition in aluminum alloy kinetic spraying process

Particle velocity is a very important parameter in kinetic spraying (or cold gas dynamic spraying). It is difficult to measure the velocity of a particle with supersonic speed at low temperature (lower than 500 °C). Thus, in many investigations only estimated values are used for evaluating coating processes. In this paper, the modeling of particle acceleration was reviewed, and the measurement of in-flight particle velocity in a kinetic spraying process was performed. Particle velocity and flux distributions from different process gas temperatures and pressures were investigated. The influences of process gas temperature and pressure on particle velocity were discussed. Characteristic of Al-Si feedstock deposition onto a mild steel substrate was described by comparing coatings structures with the in-flight particle conditions. The deposition behavior showed two critical particle velocities for Al-Si powder deposition onto a substrate and for particle-particle bonding.     

Nonlinear resonance in an acoustic system with a coagulating gas suspension

The occurrence of resonance in an acoustic system as a result of a change of the dispersion of a gas suspension filling a resonator with fixed external excitation parameters was studied. The dynamics of the medium is described by a system of equations of motion of a multi-velocity multi temperature continuum taking into account pulse and energy exchange between the carrier phase and the disperse fractions. Coagulation of different particle fractions is simulated by Smolukhovskii’s Lagrangian model, which takes into account pair collisions. Resonance occurs as a result of the resonance frequency of the system approaching a fixed external excitation frequency as a result of a change in the dispersion of the gas suspension during coagulation initiated by the addition of a small amount of the large particle fraction to the finely dispersed gas suspension. As a result of calculations, an estimate for the time of the change in dispersion of the system and the generation of resonance oscillations were obtained.     

3-D numerical analysis of EHD turbulent flow and mono-disperse charged particle transport and collection in a wire-plate ESP

The present study attempts to develop a detailed numerical approach and a simulation procedure to predict the motion of gas, ions and particles inside a simple parallel plate channel containing a single corona wire. A hybrid Finite Element (FEM)-Flux Corrected Transport (FCT)-Finite Volume (FVM) method is used: the FEM–FCT numerical algorithm is applied for modeling the steady-state corona discharge, while the turbulent gas flow and the particle motion under electrostatic forces are modeled using the commercial CFD code FLUENT. Calculations for the gas flow are carried out by solving the Reynolds-averaged Navier–Stokes equations and turbulence is modeled using the k–? turbulence model. An additional source term is added to the gas flow equation to include the effect of the electric field, obtained by solving a coupled system of the electric field and charge transport equations using User-Defined Functions (UDFs). The particle phase is simulated based on the Lagrangian approach, where a large number of particles is traced with their motion affected by the gas flow and electrostatic forces using the Discrete Phase Model (DPM) in FLUENT. The developed model is useful to gain insight into the particle collection phenomena that take place inside an ESP.