冷却速率对Lennard-Jones体系凝固过程中结构与动力学性质的影响

用分子动力学模拟方法研究了五种不同冷却速率对Lennard-Jones体系凝固过程中结构与动力学性质的影响.采用两种不同的方法来确定玻璃转变温度Tg,并且对结晶温度Tc、径向分布函数g(r)、均方位移函数MSD与扩散系数D、平均配位数进行比较分析.结果表明:冷却速率影响Lennard-Jones体系凝固过程中的结构.当使用足够高的冷却速率冷却时,体系发生玻璃化转变,而且冷却速率越快,玻璃转变温度越高;当冷却速率较小时,体系形成晶体,而且冷却速率越慢,结晶温度越高,结晶程度也越高.同时发现,冷却速率对扩散系数和平均配位数也有很大影响,二者在体系发生玻璃转变时都有一个缓变的过程,表明了过冷液相区的存在.     

真空冷却过程中实验条件对真空冷却速率的影响

以熟肉为实验材料,对实验条件对真空冷却速率的影响进行了理论分析和实验研究。实验结果表明:真空室有效容积越小、真空泵抽速越高,则真空冷却时间就会越短;冷阱温度对真空冷却速率有着明显的影响;当真空室内的最终压力在0.4～0.61kPa变化时,熟肉的表面温度一直在0℃以上,其真空冷却的时间随着真空室内压力的升高而增加。而真空室内的最终压力在0.3kPa左右时,熟肉的表面温度在真空冷却过程会低于0℃。     

原子尺寸对Lennard-Jones液体玻璃转变过程中动力学不均匀性的影响

采用分子动力学模拟技术,通过对非高斯参数α_2(t)、粘度η等动力学参数的计算探究了原子尺寸对Lennard-Jones(LJ)液体在玻璃转变过程中动力学不均匀性的影响.结果表明在玻璃转变温度T_g附近,原子尺寸越小,粘度值增加越显著,且粘度随温度的变化满足Vogel–Fulcher–Tammann (VFT)方程;β弛豫阶段的非高斯参数α_2(t)与时间的关系满足幂律函数,且随着原子尺寸的减小动力学不均匀性越来越明显.     

对“不良导体导热系数测量”实验中冷却速率求法之我见

张祖寿同志在1990年第2期(《物理实验》发表的《对“不良导体导热系数的测量”实验的研究》中,提出了计算自然冷却速率的一种新方法。文中根据牛顿冷却定律推导出计算冷却速率的理论公式 dT/dτ=1/τ′(T-T_0)1n(T′-T_0)/(T_c-T_0) (1) 式中T_0为环境温度,T_c为开始计时的温度,T′为任一时刻τ′时的温度。只要测得T_c和T′即可由(1)式计算出dT/dτ|_(r=r2)·笔者对此方法持有异见。     

导热系数的测定实验中冷却速率的测量

分析了导热系数的测定实验中冷却速率的测量方法。建议在学生实验时,采用作图与近似计算两种方法求冷却速率,可以收到较好的教学效果。     

磁感应强度和冷却速率对Tb_(0.27)Dy_(0.73)Fe_(1.95)合金凝固过程中取向行为的影响

实验研究了磁感应强度和冷却速率对Tb_(0.27)Dy_(0.73)Fe_(1.95)合金凝固过程中(Tb,Dy)Fe_2相取向行为及合金磁性能的影响.结果表明,将强磁场作用于Tb_(0.27)Dy_(0.73)Fe_(1.95)合金的凝固过程可以制备出(Tb,Dy)Fe_2相沿111取向的组织,同时显著提高了合金的磁致伸缩性能;通过提高磁感应强度可以在更快的冷却速率下得到111取向的组织;在4-10 T范围内,随着冷却速率的增加,(Tb,Dy)Fe_2相沿111取向所需的磁感应强度增加,而发生(110)取向的磁感应强度减小.随着冷却速率的增加,合金的饱和磁化强度增加,而强磁场的施加对合金饱和磁化强度的变化没有明显影响.(Tb,Dy)Fe_2相的取向行为受*Tb,Dy)Fe_3相取向行为的影响,且由磁晶各向异性能与磁场作用时间共同控制.     

冷速对液态合金Ca50Zn50快速凝固过程中微观结构演变的影响

采用分子动力学方法对六种不同冷速对原子尺寸相差较大的液态合金Ca50Znso凝固过程中微观结构演变的影响进行了模拟研究,并采用双体分布函数、Honeycutt—Andersen（HA）键型指数法、原子团类型指数法（CTIM一2）、可视化等方法进行了深入分析,结果表明：系统存在一个临界冷速,介于和5×10^11K／s与1×10^11K／s之间,在临界冷速以上（如1×10^11K／s,1×10^11K／s,1×10^11K／s和5×10^11K／s）时,系统形成以1551,1541,1431键型或二十面体基本原子团（12012000）等为主体的非晶态结构;在临界冷速以下时,系统形成以1441和1661键型或bcc基本原子团（1460800）为主体（含有少量的hcp（1200066）和fcc（12000120）基本原子团）的部分晶态结构．在非晶形成的冷速范围内,其总双体分布函数的第一峰明显分裂成与近邻分别为Zn—Zn,Ca—Zn,Ca—Ca相对应的三个次峰;且随着冷速的下降,同类原子近邻的次峰峰值升高、异类原子近邻的次峰峰值下降;Zn原子容易偏聚,随着冷速降低,二十面体的数量增多,非晶态结构也越稳定．在晶态形成的冷速范围内,Zn原子己大量偏聚形成大块bcc晶态结构,Ca原子也部分形成hcp和fcc晶态结构．     

冷却速率对液态金属Cu凝固过程中微观结构演变影响的模拟研究

采用Quantum Sutton-Chen(Q-SC)多体势对液态金属Cu在四个不同冷却速率下的凝固过程进行了分子动力学模拟研究. 通过双体分布函数、键型指数、配位数、均方位移及可视化分析, 结果表明：冷却速率对液态金属Cu的微观结构演变有决定性影响. 当冷却速率为1.0×10¹⁴K/s时得到非晶态结构；当冷速分别为1.0×10¹³K/s，1.0×10¹²K/s和1.3×10¹¹K/s时，系统形成以1421键型为主体的面心立方(fcc)与六角密集(hcp)共存的混合晶体结构；且其结晶温度分别为373K，773K和873K，即冷速越慢，其结晶温度越高，结晶程度也越高；且冷速越慢，1421键型越多，混合晶体中面心立方(fcc)结构所占的比例越高. 同时发现，原子的平均配位数的变化与1551，1441，1661键型的变化密切相关, 反映出体系对称性结构的变化规律与配位数的变化有关. 在可视化分析中，进一步采用中心原子法展现出非晶态与晶体结构的2D截面，及在3D下混合晶体中两个基本原子团分别为面心立方(fcc)与六角密集(hcp)基本原子团的具体结构. 关键词： Q-SC多体势 液态金属Cu 凝固过程 分子动力学模拟     

高温高压下Ag-Mg-Zn合金中金属间化合物的微观结构与热动力学性质的第一性原理计算

采用基于密度泛函理论的第一性原理计算方法,研究了Ag-Mg-Zn合金中金属间化合物AgMg,Mg4Zn8和Ag8Mg4Zn4在高温高压下的结构稳定性、弹性性能和热动力学性质.理论计算结果与实验值和其他的理论结果符合得非常好.研究表明:金属间化合物AgMg,Mg4Zn8和Ag8Mg4Zn4在零温零压下是力学稳定的;Mg4Zn8和Ag8Mg4Zn4为延性相,而AgMg则为脆性相;在这三种金属间化合物中,Ag8Mg4Zn4的塑性最好,AgMg的塑性最差.利用准谐Debye模型,讨论了高温高压下Ag-Mg-Zn合金中金属间化合物的摩尔振动内能Uvib,m,摩尔Helmholtz振动自由能Avib,m,摩尔振动熵Svib,m,摩尔定容热容Cv,m,摩尔定压热容Cp,m,热膨胀系数α,Griüneisen 参量γ和Debye温度(O).     

Spray flow-network flow transition of binary Lennard-Jones particle system

We simulate gas-liquid flows caused by rapid depressurization using a molecular dynamics model. The model consists of two types of Lennard-Jones particles, which we call liquid particles and gas particles. These two types of particles are distinguished by their mass and strength of interaction: a liquid particle has heavier mass and stronger interaction than a gas particle. By simulations with various initial number densities of these particles, we found that there is a transition from a spray flow to a network flow with an increase of the number density of the liquid particles. At the transition point, the size of the liquid droplets follows a power-law distribution, while it follows an exponential distribution when the number density of the liquid particles is lower than the critical value. The comparison between the transition of the model and that of models of percolation is discussed. The change of the average droplet size with the initial number density of the gas particles is also presented.     

Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟

本文用分子动力学模拟研究了Ni熔体以不同冷速凝固后微观结构的演变规律, 并通过理论计算确定出了Ni熔体凝固后获得理想非晶的临界条件. 模拟结果发现冷速小于10¹¹ K/s时, Ni 熔体凝固后形成晶态组织; 冷速在10¹¹ K/s到10^14.5 K/s之间时, Ni熔体凝固后形成由晶态结构与非晶态结构组成的混合组织. 冷速小于10¹⁰ K/s, Ni 熔体凝固后形成的晶态组织具有fcc结构; 冷速在10¹⁰ K/s到10^14.5 K/s之间时, Ni熔体凝固后组织中的晶态由fcc和hcp结构层状镶嵌排列构成. 通过分析模拟结果和计算结果, 确定出了Ni熔体凝固后形成理想非晶的临界冷速为10^14.5 K/s. 并发现Ni熔体中临界晶核(冷速等于10^14.5 K/s)和亚临界晶核(冷速大于10^14.5 K/s) 均由fcc和hcp组成层状偏聚结构, 这表明Ni熔体中生长的晶体、临界晶核和晶胚的结构是相同的. 关键词： 分子动力学模拟 晶体团簇 临界冷速 结构     

A simulation study of microstructure evolution during solidification process of liquid metal Ni

A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method (EAM) model potential. The bond-type index method of Honeycutt--Andersen (HA) and a new cluster-type index method (CTIM-2) have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution: below the crystallization temperature $T_{\rm c}$, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of $2.0\times10^{13}$\,K\,$\cdot$\,s$^{-1}$ and $1.0\times10^{12}$\,K\,$\cdot$\,s$^{-1}$, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp (expressed by (12 0 0 0 6 6) in CTIM-2) and the fcc (12 0 0 0 12 0) basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc (12 0 0 0 12 0) basic clusters consisting of 1421 bond-type, and their crystallization temperatures $T_{\rm c}$ would be 1073 and 1173\,K, respectively.     

激光器中压气体风冷系统的实验研究

本文通过模拟传热模型,建立起激光器风冷系统的换热基本方程式.分析指出:仅当热阻主要集中于光泵浦——冷却质一侧时,增加循环工作物质密度才有提高换热效果、减小流动阻力、降低风机耗功的明显效果.在此基础上研制成功重复率为10c/s,热耗散功率200W的激光器中压气体封闭循环风冷系统.实验研究所用气体压力P、流速W、壳体外壁侧热阻等诸因素对冷却效能影响的结果亦与分析结论一致.     

Lennard-Jones binary fluids: A comparative study between the molecular dynamics and Monte Carlo descriptions of their structural properties

Molecular dynamics and Monte Carlo techniques are employed for the study of binary Lennard-Jones fluids. Systematic comparisons between the predictions of both techniques are discussed, with particular emphasis on the dependency of the structural properties with respect to temperature and Lennard-Jones potential parameters.     

Influence of dispersive long-range interactions on properties of vapour–liquid equilibria and interfaces of binary Lennard-Jones mixtures

The influence of dispersive long-range interactions on properties of vapour–liquid equilibria and interfaces of six binary Lennard-Jones (LJ) mixtures was studied by molecular dynamics (MD) simulations and density gradient theory (DGT). The mixtures were investigated at a constant temperature T, at which the low-boiling component, which is the same in all mixtures, is subcritical. Two different high-boiling components were considered: one is subcritical, the other is supercritical at T. Furthermore, the unlike dispersive interaction was varied such that mixtures with three different types of phase behaviour were obtained: ideal, low-boiling azeotrope, and high-boiling azeotrope. In a first series of simulations, the full LJ potential was used to describe these mixtures. To assess the influence of the long-range interactions, these results were compared with simulations carried out with the LJ truncated and shifted (LJTS) potential applying the corresponding states principle. The dispersive long-range interactions have a significant influence on the surface tension and the interfacial thickness of the studied mixtures, whereas the relative adsorption and the enrichment are hardly affected. Furthermore, the influence of the long-range interactions on Henry's law constants and the phase envelopes of the vapour–liquid equilibrium was investigated. The long-range interactions have practically no influence on the composition dependency of the investigated mixture properties.