颗粒介质的结构及热力学

颗粒介质具有远程无序和近程有序的结构, 是产生动力学不均匀性(dynamical heterogeneity) 和复杂不可逆过程的根源. 本文分析了颗粒介质的结构特征、变形和能量耗散之间的内在关联, 讨论了颗粒介质的弹性, 提出了流变应变增量、耦合应变增量和弹性应变增量的应变增量分解方式. 沿用非平衡热力学框架, 引入表征运动无序的动理学颗粒温度T^k和表征弹性应力涨落的构型温度T^c, 作为非平衡态变量, 建立了双颗粒温度热力学(two-granular-temperature thermodynamics, TGT理论), 注重分析了不可逆过程中的热力学力和流, 并与著名的砂土内变量热力学进行了对比.     

准静态颗粒介质的弹性势能弛豫分析

颗粒体系是典型的多体相互作用体系, 具有多重的能量亚稳态. 对于准静态颗粒体系, 引入构型颗粒温度T_c描述弹性势能涨落. 本文认为平衡的体系具有一定的构型颗粒温度T_a, 其量值反映了其结构特征. 当外界扰动激发的构型颗粒温度超出T_a时, 产生不可逆过程. 通过对应力松弛过程的分析, 发现(T_c-T_a)激发了弹性弛豫, 且(T_c-T_a)越大则松弛过程中应力变化越大, 最终构型颗粒温度T_c→T_a时,宏观应力松弛结束,体系达到新的能量亚稳态.     

锆基非晶合金的动态弛豫和应力松弛

非晶合金动态弛豫、变形和微观结构非均匀性之间的关联是非晶态物理领域重要研究内容之一.本文选取玻璃形成能力优良,热稳定性好的Zr₄₈(Cu_5/6Ag_1/6)₄₄Al₈块体非晶合金作为研究载体,借助于动态力学分析及应力松弛实验探究了温度和结构弛豫对非晶合金动态弛豫和变形行为的影响.研究结果表明非晶合金动态弛豫行为对温度极其敏感,随温度升高表现为等构型、老化、玻璃转变和晶化4个阶段.基于玻璃转变温度之下等温测试,结构弛豫行为导致非平衡高能量状态非晶合金向低能量状态迁移,内耗随时间演化规律可用Kohlrausch-Williams-Watts (KWW)扩展指数方程描述.此外,基于KWW方程和激活能谱等方法分析了模型合金应力松弛过程中非均匀结构激活,该过程涉及弹性变形向非弹性变形的转变.由于非晶合金存在微观尺度的结构非均匀性与能量起伏,变形单元的激活并非为单一特征时间,而是服从一定分布.通过考虑对数时间尺度和变形单元激活的特征时间分布分别为对称Gauss分布和非对称Gumbel分布,可...     

晶界滞弹性弛豫理论的现代进展

上世纪中期,人们通过扭摆试验测量内耗,发现晶界滞弹性弛豫峰．后来尽管很多学者提出了各种理论模型,但晶界滞弹性弛豫的微观机理仍然是不清楚的．最近,Xu根据弹性应力引起晶界溶质偏聚或贫化实验结果,提出了晶界滞弹性弛豫的微观机理是晶界吸收或发射空位,建立了晶界滞弹性弛豫的平衡方程和动力学方程,解析地表述了晶界滞弹性弛豫过程,并成功地阐明了普遍存在于金属中的中温脆性峰温度移动现象．本文将综述晶界滞弹性弛豫理论的这一现代发展．     

准静态颗粒流流动规律的热力学分析

本文分析了颗粒流的介观结构及其特征,提出了颗粒流的双颗粒温度概念Tkin和Tconf,表征颗粒无序运动和构型无序演化的程度;进而作为非平衡变量,与经典非平衡热力学(classical irreversible thermodynamics,CIT)变量共同构成颗粒流的热力学状态变量集,确定了颗粒流的能量转换规律和熵产生率等,发展了颗粒流双颗粒温度(two granular temperate,TGT)模型.以体积恒定的简单剪切准静态颗粒流为例,结合离散元模拟(discrete element method,DEM),确定了双颗粒温度模型所需的材料参数,分析了颗粒流发展段的规律和稳恒段的有效摩擦系数.     

力和扩散机理下外延形貌的演化分析

本文提出了一个新的基于扩散界面的相场模型来描述外延生长中岛的形核、生长及熟化过程. 该模型同时考虑了弹性场、表面能、沉积、扩散、解吸和能量势垒等热力学及动力学过程对表面纳米形貌的影响. 采用经典的BCF模型来描述生长中的扩散形核过程, 而采用一个新的包含弹性应变能的自由能函数, 通过变分得到一个描述多层岛生长的相场方程, 该方法可以有效地描述外延生长中复杂的外延形貌. 采用有限差分格式对非线性耦合方程组进行求解. 数值结果显示, 该模型可以真实地再现外延生长中多层岛结构(即山丘状形貌)的演化过程, 模拟结果与已有实验结果一致. 同时模拟了生长过程中随外延形貌演化而形成的复杂生长应力, 研究表明, 在生长过程中, 岛中存在着复杂的应力分布, 且在岛边界处应力达到局部最大, 这与实验结果定性一致. 此外, 本文的重要发现是, 外延生长中的应力演化明显地影响原子的扩散过程, 当应力存在时, 外延结构变化较无弹性场时变快. 该项研究对理解外延生长中各物理机理的协同作用有重要的指导意义.     

颗粒物质的多尺度结构及其研究框架

颗粒物质力学的研究刚刚起步,目前很多方面还不完善.文章作者认为,该学科是应用性很强的基础学科,因此提出了从实际应用中去发现科学问题,进而促进颗粒物质力学基础理论发展的研究思路,并以土力学为例做了说明.土是岩石经风化作用,由重力、流水和风力等搬运和沉积而成的产物,是密集颗粒物质体系之一.土粒是构成土体骨架、传递载荷的基本单元,与颗粒间复杂分布的孔隙水、气体共同决定了土体的非线件本构关系、剪胀(剪缩)和应力路径相关等复杂特征.以研究密集分布颗粒体系的颗粒物质力学在近20年内得到充分发展,它侧重现象的机理分析和实验的精细检测,为土力学的基础研究提供了重要启示.基于文章作者多年土力学和颗粒物质力学的研究经验,提出了土体具有多尺度结构的观点:除微观的单颗粒尺度和宏观土体尺度外,细观尺度的力链是颗粒接触力传递的路径,足存在于土体内的相对稳定的结构体.建立了初步的理论研究框架,提出了力链网络的复杂动力学响应决定土体复杂本构关系的基本没想.下一阶段将从理论分析、物理试验和基于自主开发的颗粒离散元模拟3个方面进行研究,逐渐充实土力学多尺度理论体系,以期取得突破.     

激光等离子体去除微纳颗粒的热力学研究

微杂质污染一直是影响精密器件制造质量和使用寿命的关键因素之一.对于微纳米杂质颗粒用传统的清洗方式(超声清洗等)难以去除,而激光等离子体冲击波具有高压特性,可以实现纳米量级杂质颗粒的去除,具有很大的应用潜力.本文主要研究了激光等离子体去除微纳米颗粒过程中的热力学效应:实验研究了激光等离子体在不同脉冲数下对Si基底上Al颗粒去除后的颗粒形貌变化,发现大颗粒会发生破碎而转变成小颗粒,一些颗粒达到熔点后发生相变形成光滑球体,这源于等离子体的热力学效应共同作用的结果.为了研究微粒物态转化过程,基于冲击波传播理论研究,得到冲击波压强与温度特性的演化规律;同时,利用有限元模拟方式研究激光等离子冲击波压强和温度对微粒作用规律,得到了颗粒内随时间变化的应力分布和温度分布,并在此基础上得到等离子体对颗粒的热力学作用机制.     

稠密气固流动二阶矩摩擦应力模型和模拟

基于统一二阶矩两相湍流应力模型和稠密颗粒动力学理论,建立了欧拉-欧拉双流体二阶矩颗粒摩擦应力模型。模型充分反映各向异性的气固两相相间雷诺应力相互作用,并引入有效颗粒弹性恢复系数,考虑了因颗粒表面不光滑产生的摩擦力对湍流流动结构和颗粒弥散特性的影响,对于下降管的计算结果与实验吻合较好。因颗粒摩擦产生能量耗散降低了颗粒温度和导致颗粒雷诺应力再分布。在入口和出口区域因颗粒碰撞频率较高而产生的能量耗散对流动结构影响明显。     

激光二极管阵列端面抽运复合棒状激光器热效应的有限元法分析

建立了激光二极管阵列(LDA)端面抽运棒状激光介质的数值模型.考虑到介质与空气的对流换热和介质的热力学参数随温度的变化,根据经典热传导方程和热弹性方程,运用有限元法得出了复合棒状介质和未复合棒状介质内瞬态温度、热应力和应变的时空分布,分析了温度、热应力和应变随抽运功率、换热系数和时间的变化规律.结果表明,复合棒的最高温度、最大张应力和最大轴向应变的位置与未复合棒不同,并且数值分别为未复合棒的73%,60%和33%.由此可知,利用复合棒可极大地减小热效应的影响.理论分析结果可为LDA抽运固体激光器的结构优化设计和实验研究提供理论参考.     

A new statistical mechanics approach to dense granular media

A new statistical mechanics approach to dense granular media is presented. The thermodynamic formalism is set out directly in terms of elastic potential energy, such that the configurational temperature (the intensive property which defines the steady state) relates to a quadratic function of the stresses (rather than other linear functions used in recent developments). Dense granular media are considered as canonical ensembles of noninteracting clusters, which can be identified with repeatable equilibrium configurations. Then, particles can be located in a new proposed phase space (conceived to separate the elastic potential energy levels). Although the importance of this paper lies in the method itself, it has been illustratively applied to the simple case of two-dimensional (2D) dense granular media (an arrangement of frictionless monodisperse elastic disks under isotropic horizontal stress compression). In this case, the temperature is directly replaced by the squared external pressure, and the packing ratio of the most probable microstate is close to the reported value of random close packing. Moreover, some interesting general conclusions arise.     

垂直载荷下湿玻璃珠样品的声速和衰减测量及非线性声学响应

为了分析垂直载荷下颗粒物质的声速、声衰减系数、谐波非线性等特性,本工作采用飞行时间法测量了不同含水量下声速随压强的变化规律,并利用傅里叶变换法分析了干、湿玻璃珠样品的声衰减和非线性声学特性。结果表明:干、湿玻璃珠样品中的声速、声衰减系数以及谐波非线性均随压强呈幂律变化;湿颗粒样品中随着液体含量增多,声速逐渐增加,超声波透过湿颗粒样品时的能量耗散和非线性逐渐减小。分析原因表明,压强和孔隙流体改变了颗粒之间的接触分布,使得颗粒体系的声速、声衰减以及谐波非线性等特性都随之发生变化。     

Some properties of an ultimately nonideal metastable supercooled plasma

The formation of a metastable plasma state from an initial strongly nonideal state and certain properties of this metastable supercooled state are considered by directly modeling the dynamics of a system of many Coulomb particles. It is shown that the relaxation of the average kinetic energy of the particles can be characterized by a universal dimensionless function and, in particular, that an ultimate value degree of plasma nonidealness occurs which can be reached in the metastable state when no external action is present. The calculated binary correlation functions are in good agreement with the results of the Debye model, even outside its range of applicability. The time dependence of the total dipole moment of the system of particles is investigated. It is shown that there are oscillations of the total dipole moment. These collective oscillations take place at a frequency which is somewhat below the Langmuir frequency, the oscillations of the free and bound electrons occurring in antiphase. A hypothesis is advanced to the effect that recombination relaxation is frozen on account of the interaction of quasibound electrons with the Langmuir oscillations of the free electrons. Institute of General Physics, Academy of Sciences of the Russian Federation. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 10, pp. 3–16, October, 1996.     

Possible test of the thermodynamic approach to granular media

We study the steady state distribution of the energy of the Sherrington-Kirkpatrick model driven by a tapping mechanism which mimics the mechanically driven dynamics of granular media. The dynamics consists of two phases: a zero temperature relaxation phase which leads the system to a metastable state, then a tapping which excites the system thus reactivating the relaxational dynamics. Numerically, we investigate whether the distribution of the energies of the blocked states obtained agrees with a simple canonical form of the Edwards measure. It is found that this canonical measure is in good agreement with the dynamically measured energy distribution. A possible experimental test of the Edwards measure based on the study here is proposed.     

Limiting nonideal metastable supercooled plasma

Direct modeling of the dynamics of a system of many Coulomb particles is applied to analyze the formation stage of a metastable plasma state from an initial, highly nonideal state, and also to consider some properties of this metastable supercooled state. It is shown that relaxation of the average particle kinetic energy may be characterized by a universal dimensionless function and in particular, there is a limiting degree of plasma nonideality which may be achieved in the metastable state, in the absence of any external influence. The calculated pair correlation functions agree with the results of the Debye model, even outside its limits of validity. The time dependence of the total dipole moment of the particle system is investigated. It is shown that oscillations of the total dipole moment are observed. These collective oscillations take place at a frequency slightly below the Langmuir frequency and the oscillations of free and bound electrons are in antiphase. The hypothesis is put forward that recombination relaxation is frozen as a result of interaction between quasibound electrons and Langmuir oscillations of free electrons. Zh. Tekh. Fiz. 67, 42–52 (August 1997)     

Spatial Density Distributions and Correlations in a Quasi-one-Dimensional Polydisperse Granular Gas

By Monte Carlo simulations, the effect of the dispersion of particle size distribution on the spatial density distributions and correlations of a quasi one-dimensional polydisperse granular gas with fractal size distribution is investigated in the same inelasticity. The dispersive degree of the particle size distribution can be measured by a fractal dimension dr, and the smooth particles are constrained to move along a circle of length L, colliding inelastically with each other and thermalized by a viscosity heat bath. When the typical relaxation time τ of the driving Brownian process is longer than the mean collision time To, the system can reach a nonequilibrium steady state. The average energy of the system decays exponentially with time towards a stable asymptotic value, and the energy relaxation time τB to the steady state becomes shorter with increasing values of df. In the steady state, the spatial density distribution becomes more clusterized as df increases, which can be quantitatively characterized by statistical entropy of the system. Furthermore, the spatial correlation functions of density and velocities are found to be a power-law form for small separation distance of particles, and both of the correlations become stronger with the increase of df. Also, tile density clusterization is explained from the correlations.     

Self-Similar Asymptotics for the Boltzmann Equation with Inelastic and Elastic Interactions

We consider some questions related to the self-similar asymptotics in the kinetic theory of both elastic and inelastic particles. In the second case we have in mind granular materials, when the model of hard spheres with inelastic collisions is replaced by a Maxwell model, characterized by a collision frequency independent of the relative speed of the colliding particles. We first discuss how to define the n-dimensional (n = 1,2,...) inelastic Maxwell model and its connection with the more basic Boltzmann equation for inelastic hard spheres. Then we consider both elastic and inelastic Maxwell models from a unified viewpoint. We prove the existence of (positive in the inelastic case) self-similar solutions with finite energy and investigate their role in large time asymptotics. It is proved that a recent conjecture by Ernst and Brito devoted to high energy tails for inelastic Maxwell particles is true for a certain class of initial data which includes Maxwellians. We also prove that the self-similar asymptotics for high energies is typical for some classes of solutions of the classical (elastic) Boltzmann equation for Maxwell molecules. New classes of (not necessarily positive) finite-energy eternal solutions of this equation are also studied.     

The dynamics of multiparticle collisions in motion of a granular material

We consider the complex problem of how to calculate particle motions taking into account multiparticle collisions. Multiparticle contacts occur when a particle collides with neighbouring particles, so that those contacts have a direct influence on each other. We will focus on the molecular dynamics method. Particularly, we will analyse what happens in cohesive materials during multiparticle contacts. We investigated the expression of repulsive force formulated under fractional calculus which is able to control dynamically the transfer and dissipation of energy in granular media. Such approach allows to perform simulations of arbitrary multiparticle collisions and also granular cohesion dynamics.     

Radiative and radiationless processes of Er3+ in MnF2

We have designed a method for the measurement of radiative and radiationless transition rates of fluorescent systems with several metastable levels. This method applies to systems with adjacent metastable levels connected by multiphonon relaxation processes. The measurements required are those of the integrated fluorescence intensities and decay times at different temperatures; no absorption or excitation measurements are required. We have applied this method to the study of the fluorescence of Er³⁺ in MnF₂. The knowledge of both radiative and non- radiative decays has allowed us to reconstruct the kinetics of the excitation and de-excitation processes and to find the dependence on temperature and energy gap of multiphonon relaxation processes in this system.