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

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

分形结构对随机取向烟尘团簇粒子光散射特性的影响

利用蒙特卡罗方法对不同分形维数和分形前向因子的随机取向烟尘团簇粒子的分形结构进行了仿真,采用离散偶极子近似（DDA）方法对随机取向烟尘团簇粒子的缪勒矩阵元进行了数值计算,并与球形粒子模型进行了比较,深入探讨了烟尘团簇粒子的分形维数和分形前因子对其散射特性的影响。研究表明,等效球形粒子的光散射特性与随机取向烟尘团簇粒子的光散射特性存在很大差别,并且此差别随着团簇粒子的分形维数以及分形前向因子的增大而减小;分形维数对表征团簇粒子散射特性的缪勒矩阵元的影响在一定散射角范围内均比较明显,分形前向因子对团簇粒子的缪勒矩阵元角分布的影响与分形维数的影响类似,不过其影响相对分形维数较弱。     

火灾烟颗粒分形模型和球形模型光散射的比较研究

对烟颗粒的光散射进行模拟计算是研究火灾烟颗粒光散射特性的重要手段,目前对于火灾烟颗粒光散射的数值计算多采用球形或椭球模型.实际上,火灾烟颗粒的形貌与球形和椭球均存在着显著差异.扫描电子显微镜图像表明,烟颗粒具有近似分形的结构.本文利用离散偶极近似方法计算了随机取向的火灾烟颗粒分形凝团以及同体积的球形颗粒的光散射Muller矩阵,并对两者的归一化Muller矩阵元素随散射角的分布进行了比较.研究表明:火灾烟颗粒分形模型和球形模型的归一化矩阵元素F₁₁(θ)/     

随机分布烟尘团簇粒子辐射特性研究

基于分形理论,采用蒙特卡罗方法对随机分布的烟尘团簇粒子结构进行了仿真模拟,利用离散偶极子近似(discrete dipole approximation, DDA)方法研究了随机分布的烟尘团簇粒子的辐射特性,分析讨论了分形维数、原始微粒粒径和数量以及复折射率对随机分布烟尘团簇粒子辐射特性的影响.研究表明,在给定分形维数的情况下,烟尘团簇粒子的辐射特性取决于原始微粒粒径、数量及复折射率;原始微粒较小的团簇粒子,当分形维数较小时,吸收截面变化不明显,但当分形维数大于2时,吸收截面骤然增大,然而,对于具有比较大的原始微粒粒径、数量及复折射率的烟尘团簇粒子,吸收截面随着分形维数的增大而单调递减;随着分形维数的增大,团簇粒子的散射截面、消光截面及单次散射反照率均单调递增;从整体上来讲,团簇粒子的辐射特性与等效球形粒子的辐射特性存在着比较大的差别,并且这种差别随着分形维数的增大而减小.该工作对研究气溶胶粒子的辐射及气候效应具有重要的科学价值. 关键词： 烟尘团簇粒子 辐射特性 离散偶极子近似方法     

熔石英亚表面球形杂质对入射光场的调制作用

建立激光辐照下熔石英体内球形杂质引起元件损伤的模型。基于Mie理论计算杂质的散射系数,研究杂质颗粒附近的调制光场分布,分析光强增强因子(LIEF)与杂质折射率及半径的关系。结果显示,杂质颗粒半径小于40 nm时,所有类型杂质对光场的调制作用可以忽略;颗粒半径大于40 nm时,对于折射率小于熔石英的非耗散杂质,LIEF随颗粒半径的增大和折射率的减小而增大,且后向散射远强于前向散射;对于折射率大于熔石英的非耗散杂质,随着折射率和半径的增大,LIEF整体趋势增大,局部达到102量级。对于耗散杂质,前向散射强度随颗粒半径的增大先增大后减小,当颗粒半径大于170 nm时,后向散射强于前向散射。     

熔石英亚表面球形杂质对入射光场的调制作用EI北大核心CSCD

建立激光辐照下熔石英体内球形杂质引起元件损伤的模型。基于Mie理论计算杂质的散射系数,研究杂质颗粒附近的调制光场分布,分析光强增强因子(LIEF)与杂质折射率及半径的关系。结果显示,杂质颗粒半径小于40 nm时,所有类型杂质对光场的调制作用可以忽略;颗粒半径大于40 nm时,对于折射率小于熔石英的非耗散杂质,LIEF随颗粒半径的增大和折射率的减小而增大,且后向散射远强于前向散射;对于折射率大于熔石英的非耗散杂质,随着折射率和半径的增大,LIEF整体趋势增大,局部达到102量级。对于耗散杂质,前向散射强度随颗粒半径的增大先增大后减小,当颗粒半径大于170 nm时,后向散射强于前向散射。     

双分散颗粒体系在临界堵塞态的结构特征

堵塞行为是颗粒体系中一种常见的现象,其力学性质与堆积结构的关联非常复杂.本文采用离散元法研究了由两种不同半径颗粒组成的二维双分散无摩擦球形颗粒体系在临界堵塞态所呈现的结构特征,讨论了大小颗粒粒径比与大颗粒百分比对临界堵塞态的影响.数值模拟结果表明,当粒径比小于1.4时,临界平均接触数与大颗粒百分比关系不大,当粒径比大于1.4时随着大颗粒百分比的增大临界平均接触数先减小再增大.而临界体积分数在粒径比小于1.8时随着大颗粒百分比的增加先减小后增大,大于1.8时又基本不随大颗粒百分比而变化.大颗粒百分比在接近0或1时,系统近似为单分散体系,临界平均接触数与体积分数基本不随半径比的增大而变化;在接近0.5时,临界平均接触数随着半径比的增大逐渐减小,而临界体积分数则是先减小后增大.文中对大-小颗粒这一接触类型的百分比也进行了探讨,其值随着大颗粒百分比的增大呈二次函数的变化趋势,粒径比对这一变化趋势只有较小的影响.     

非球形气溶胶粒子短波红外散射特性研究

利用T矩阵方法,以及基于扩散限制凝聚理论的广义多粒子米散射方法,研究了多种气溶胶粒子在1.6和2.0μm波段处,非球形单粒子和团簇粒子的光散射辐射特性,并分析了粒子有效半径、复折射指数、粒子形状、相对湿度等因素对非球形粒子散射特性的影响.分析表明,除了粒子有效半径和形状会在不同程度上引起粒子散射特性变化,相对湿度对其影响也比较大,球形粒子与非球形粒子在不同相对湿度下后向散射相对差异均在18%以上;当粒子体积较小时,水溶性气溶胶的后向散射强度随相对湿度的增加而增强,而当粒子体积较大时,则随相对湿度的增加而减弱;在体积相同的条件下,体积较小的团簇粒子的不对称因子比非球形单粒子平均偏大0.023,而体积较大的团簇粒子,却比非球形单粒子不对称因子平均偏小0.055;单粒子或等体积的团簇粒子,其不同波段之间单次散射反照率差异较大,最大可达0.226.该工作对研究气溶胶多次散射对CO₂浓度卫星反演精度影响具有重要的科学意义.     

梯度分布球颗粒复合介质的介电特性和光吸收特性

研究球形颗粒复合介质沿颗粒半径方向有梯度分布的有效介电响应和光吸收性质,利用微分有效偶极距近似方法计算有梯度分布球颗粒的等效介电常数.数值计算结果表明,对于一般的球形粒子,介电常数随半径按幂指数规律变化时,系统的有效介电常数随指数n的增大而减小;而对金属球形粒子,则存在等离子体共振吸收峰,且随参数mw和mv的增大,其位置将发生红移.这为复合材料的制备提供了一种新的理论依据.     

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

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

Numerical Simulation of Random Close Packings in Particle Deformation from Spheres to Cubes

Variation of packing density in particle deforming from spheres to cubes is studied. A new model is presented to describe particle deformation between different particle shapes. Deformation is simulated by relative motion of component spheres in the sphere assembly model of a particle. Random close packings of particles in deformation form spheres to cubes are simulated with an improved relaxation algorithm. Packings in both 2D and 3D cases are simulated. With the simulations, we find that the packing density increases while the particle sphericity decreases in the deformation. Spheres and cubes give the minimum (0.6404) and maximum (0.7755) of packing density in the deformation respectively. In each deforming step, packings starting from a random configuration and from the final packing of last deforming step are both simulated. The packing density in the latter case is larger than the former in two dimensions, but is smaller in three dimensions. The deformation model can be applied to other particle shapes as well.     

Collective motion of polar active particles on a sphere

Collective motion of active particles with polar alignment is investigated on a sphere. We discussed the factors that affect particle swarm motion and define an order parameter that can show the degree of particle swarm motion. In the model, we added a polar alignment strength, along with Gaussian curvature, affecting particles swarm motion. We find that when the force exceeds a certain limit, the order parameter will decrease with the increase of the force. Combined with our definition of order parameter and observation of the model, the reason is that particles begin to move side by side under the influence of polar forces. In addition, the effects of velocity, rotational diffusion coefficient, and packing fraction on particle swarm motion are discussed. It is found that the rotational diffusion coefficient and the packing fraction have a great influence on the clustering motion of particles, while the velocity has little influence on the clustering motion of particles.     

Photon sphere and phase transition of d-dimensional (d ≥ 5) charged Gauss-Bonnet AdS black holes

Motivated by recent work, nonmonotonic behaviors of photon sphere radius can be used to reflect black hole phase transition for Reissner-Nordström-AdS (RN-AdS) black holes, we study the case of five-dimensional charged Gauss-Bonnet-AdS black holes in the reduced parameter space. We find that the nonmonotonic behaviors of photon sphere radius still exist. Using the coexistence line calculated from P-V plane, we capture the photon sphere radius of saturated small and large black holes (the boundary of the coexistence phase), then illustrate the reduced coexistence region. The results show that, reduced coexistence region decreases with charge Q but increases with Gauss-Bonnet coefficient α. When the charge vanishes, reduced coexistence region does not vary with Gauss-Bonnet coefficient α any more. In this case, the Gauss-Bonnet coefficient α plays the same role as the charge of five-dimensional RN-AdS black holes. Also, the situation of higher dimension is studied in the end.     

Backscattering and negative polarization of agglomerate particles

We used the discrete dipole approximation to study the backscattering of agglomerate particles consisting of oblong monomers. We varied the aspect ratio of the monomers from approximately 1 (sphere) to 4, while we kept the total particle volume equivalent to that of an x = 10 sphere for m = 1.59 + i0 and 1.50 + i0 and considered two values of agglomerate packing density: rho = 0.25 and rho = 0.1. We found that these particles do not display a prominent brightness opposition effect but do produce significant negative polarization over a range of near-backscattering angles. Increasing the monomers' aspect ratio can make the negative polarization much more prominent. We have noted also that decreasing m and p can reduce the amplitude of the negative polarization for these particles.     

Application of Fractal Geometry in Quantitative Characterization of Aerosol Morphology

The aerosol morphologies of six types of particles plus one other from a different source were characterized quantitatively in terms of fractal dimensions based on boundary, projected area and surface roughness. Except for the particle with a smooth boundary, multifractal features were observed on the projected boundary of aerosols. It was found that the fractal dimensions of the aerosol boundaries in the high-resolution region ranged from 1.00 to 1.13, whereas those in the low-resolution region were distributed from 1.04 to 1.91. The compactness of a particle could be examined in terms of the fractal dimension of its projected area such that a solid particle gave a value of 2.00, a hollow particle with pores inside gave a value of 1.93 and a chain agglomerate gave an even lower value of 1.81. Aerosol topological features retrieved from the gray levels of a microimage were used to illustrate the difference in surface roughness or distinguish particle deposition on the substrate from that of agglomeration. A value of 2.49 was found for a hollow sphere in contrast to a value of 2.26 for a round pollen particle. Finally, a coordinate system was established by employing the computed fractal dimensions as axes to accommodate the particles characterized in this work and the Euclidean distance of a point from the origin was shown to be a potential composite index for aerosol morphology.     

Fractal substructure of a nanopowder

The structural evolution of a nanopowder by repeated dispersion and settling can lead to characteristic fractal substructures. This is shown by numerical simulations of a two-dimensional model agglomerate of adhesive rigid particles. The agglomerate is cut into fragments of a characteristic size l, which then are settling under gravity. Repeating this procedure converges to a loosely packed structure, the properties of which are investigated: (a) The final packing density is independent of the initialization, (b) the short-range correlation function is independent of the fragment size, (c) the structure is fractal up to the fragmentation scale l with a fractal dimension close to 1.7, and (d) the relaxation time increases linearly with l.     

The Fractal Characteristic of Particles in Granular Material Flows and Its Effect on Effective Thermal Conductivity

According to the fact that many pulverized particles possess fractal characteristic, a fractal model for studying fine particles in granular material flows is first proposed. An expression of particles' fractal distribution is derived to describe the relationship between the particle fractal dimensions and particle velocity distribution function. In accordance with this model, the theoretical particle effective thermal conductivity is derived. The analytical results show that for the small Biot-Fourier number, the effective thermal conductivity increases with the square root of the granular temperature. For very large Biot-Fourier number, the effective thermal conductivity linearly increases with the granular temperature. Numerically calculated results show that the thermal conductivity increases with the particle size fractal dimensions and decreases with the particle surface fractal dimensions.     

活性Janus粒子体系中的定向运动涌现行为

本文使用过阻尼郎之万方程动力学模拟方法模拟了二维体系中活性Janus粒子的群体动力学行为. 体系随着粒子密度的增加，存在一个活性粒子的定向运动涌现行为. 此现象在无各向异性相互作用或者无活性驱动的体系中无法观察到. 另外涌现行为发生的临界转变密度会随着活性大小和各向异性相互作用强度的增大而减小，且通过调节活性大小和各向异性相互作用强度大小可以控制体系不同状态的转变. 最后借助动力学理论分析，重复出了模拟结果的基本特征. 结果论证了单个粒子层次上的各向异性相互作用足以导致活性体系进入稳定的群体定向运动.