Dissipative particle dynamics simulation on the rheological properties of heavy crude oil

The rheological properties of heavy crude oil have a significant impact on the production, refining and transportation. In this paper, dissipative particle dynamics (DPD) simulations were performed to study the effects of the addition of light crude oil and emulsification on the rheological properties of heavy crude oil. The simulation results reflected that the addition of light crude oil reduced the viscosity effectively. The shear thinning behaviour of crude oil mixtures were becoming less distinct as the increase of the mass fraction of light crude oil. According to the statistics, the shear had an influence on the aggregation and spatial orientation of asphaltene molecules. In addition, the relationship between the viscosity and the oil mass fraction was investigated in the simulations of emulsion systems. The viscosity increased with the oil mass fraction slowly in oil-in-water emulsions. When the oil mass fraction was higher than 50%, the increase became much faster since systems had been converted into water-in-oil emulsions. The equilibrated morphologies of emulsion systems were shown to illustrate the phase inversion. The surfactant-like feature of asphaltenes was also studied in the simulations.     

Effect of shear equilibrium flow in Tokamak plasma on resistive wall modes

A code named LARWM with non-ideal magnetohydrodynamic equations in cylindrical model is used to describe the instability in Tokamak plasma surrounded by a conducting wall with finite resistivity. We mainly take three factors related to the shear equilibrium plasma flow into consideration to study the stabilizing effect of the shear flow on the resistive wall modes (RWMs). The three factors are the velocity amplitude of flow, the shear rate of flow on plasma surface, and the inertial energy of equilibrium plasma flow. In addition, a local shear plasma flow is also calculated by the LARWM code. Consequently, it is found that the inertial energy of the shear equilibrium plasma flow has an important role in the stabilization of the RWMs.     

碳纳米锥力学特性的分子动力学研究

结合原子间短程作用势(Brenner势)和长程作用势(Lennard-Jones势),利用分子动力学方法对各种锥角的碳纳米锥进行拉伸和压缩实验,获得其载荷-应变关系曲线、受拉/压载荷极限、应变极限和构形演变等力学特性,并与等量原子组成的碳纳米管进行比较研究.研究结果表明,等量碳原子组成的碳纳米锥的受拉/压载荷极限随着锥角的增大先是增大后减小,受拉/压应变极限则随着锥角的增大而增大.与碳纳米锥相比,等量碳原子组成的碳纳米管的受拉/压载荷极限和应变极限显得既不突出也不逊色.在受压构形演化方面,与碳纳米管丰富的径向屈曲/扭转/侧向屈曲组合形变不同,112.88°和83.62°锥角的碳纳米锥受压沿轴向完美内陷,而60.0°和38.94°锥角的碳纳米锥受压发生侧向屈曲.     

铜原子纳米团簇热力学性质的分子动力学模拟研究

利用分子动力学模拟方法,研究了CuN(N=80,140,216,312,408,500,628和736)纳米团簇在热化和冷凝过程中结构和热力学性质的变化,模型采用的是Johnson的EAM作用势.模拟结果表明:铜团簇的熔点和凝固点随其尺寸线性增加,并逐渐向大块晶体靠拢;所有团簇的凝固点都低于熔点,出现凝固过程中的滞后现象;在熔点和凝固点附近,团簇都具有负热容特性,负热容是由相变前后团簇内部结构突变引起的.     

Calculation of the thermodynamic properties of copper by molecular dynamics simulation

It is shown that the thermodynamic characteristics of a system can be accurately described using many-body potentials of the interatomic interaction in molecular dynamics calculations. Original Russian Text ? O.V. Stepanyuk, D.B. Alekseev, A.M. Saletskii, 2009, published in Vestnik Moskovskogo Universiteta. Fizika, 2009, No. 2, pp. 115–116.     

二维颗粒流通道宽度效应的分子动力学模拟

通过用分子动力学方法对颗粒物质流的计算机模拟，研究发现增大通道宽度可以使二维颗粒流从稀疏流转变为密集流状态.通过对不同通道宽度下，固定开口为9.5d的颗粒流和 漏斗口以上9.5d×8d区域记录的模拟结果分析，发现随通道宽度增大，密度变大、温 度降低.当“颗粒温度”T较低时(T/m<0.05 J/kg),颗粒流内部接触数开始超过1.2 ，同时出现较为牢固的横向链状颗粒团簇，是造成流量突变以及密集流的原因. 关键词： 颗粒物质 颗粒流 计算机模拟     

Lattice Boltzmann simulation of the two-dimensional flow of a magnetic suspension between two parallel walls under a non-uniform magnetic field

We have developed a lattice Boltzmann method based on fluctuation hydrodynamics that is applicable to the flow problem of a particle suspension. In this method, we have introduced the viscosity-modifying method, rather than the velocity-scaling method, in which a modified viscosity is used for generating random forces in lattice Boltzmann simulations. The viscosity-modifying method is found to be applicable to the simulation of a magnetic particle suspension. We have applied this method to the two-dimensional Poiseuille flow of a magnetic suspension between two parallel walls in order to investigate the behavior of magnetic particles in a non-uniform applied magnetic field. From the results of the snapshots, the pair correlation function between the magnetic pole and the magnetic particles and the averaged local particle velocity and magnetization distributions, it was observed that the behavior of the magnetic particles changes significantly depending upon which factor dominates the phenomenon in the balance between the magnetic particle–particle interaction, the non-uniform applied magnetic field and the translational and rotational Brownian motion.     

Analytic solution to the problem of the couette flow in a plane channel with infinitely large parallel walls

An analytic (in the form of a Neuman series) solution to the problem of the Couette flow in a plane channel with infinitely large parallel walls is constructed using the kinetic approach in the isothermal approximation. For the basic equation, the Bhatnagar-Gross-Krook (BGK) model of the kinetic Boltzmann equation is used, while the boundary condition is determined by the diffuse reflection model. The mass flux through half the channel thickness in the direction parallel to the channel walls as well as the nonzero component of the viscous stress tensor are calculated taking into account the constructed distribution function. The results are compared with analogous data obtained by numerical methods.     

Shear zones and wall slip in the capillary flow of concentrated colloidal suspensions

We image the flow of a nearly random close packed, hard-sphere colloidal suspension (a "paste") in a square capillary using confocal microscopy. The flow consists of a "plug" in the center while shear occurs localized adjacent to the channel walls, reminiscent of yield-stress fluid behavior. However, the observed scaling of the velocity profiles with the flow rate strongly contrasts yield-stress fluid predictions. Instead, the velocity profiles can be captured by a theory of stress fluctuations originally developed for chute flow of dry granular media. We verified this both for smooth and rough walls.     

大尺寸贵金属纳米团簇热力学性质的分子动力学模拟研究

采用分子动力学方法和原子嵌入势模拟了大尺寸金(n=1136～1556)、银(n=1088～1724)、铜(n=1000～1600)、铂(n=1004～1800)原子纳米团簇的熔化过程,得出了相应纳米团簇的势能随温度的变化曲线以及热容量随温度的变化曲线,研究了各种原子纳米团簇熔点与其团簇尺寸的关系.模拟结果表明团簇的熔点随团簇尺寸增大而升高,并逐渐向大块晶体靠拢.所有纳米团簇在熔化过程中在熔点附近都出现负热容现象,通过对团簇熔化前后结构的比较,分析了导致这种现象的原因.     

第十五届全国原子与分子物理学术会议:大尺寸贵金属纳米团簇热力学性质的分子动力学模拟研究

采用分子动力学方法和原子嵌入模型势模拟了大尺寸金（n=1136--1556）、银（n=1088--1724）、铜（n=1000--1600）、铂（n=1004--1800）原子纳米团簇的熔化过程,得出了相应纳米团簇的势能随温度的变化曲线以及热容量随温度的变化曲线,研究了各种原子纳米团簇熔点与其团簇尺寸的关系。模拟结果表明团簇的熔点随团簇尺寸增大而升高,并逐渐向大块晶体靠拢。所有纳米团簇在熔化过程中在熔点附近都出现负热容现象,通过对团簇熔化前后结构的比较研,分析了导致这种现象的原因。     

Effect of slow dynamics on elastic properties of materials with residual and shear strains

The effects of fast and slow dynamics in acoustic resonators from rock and metal samples from the D16 microcrystalline aluminum alloy are studied using resonance ultrasonic spectroscopy. Before the experiment, residual shear strains were artificially created in the metal samples. A decrease in the elasticity modulus in the fast dynamics process has been revealed in resonators from rock and the D16 alloy with residual strains. Based on an analysis of the experimental results, the following conclusion was drawn: the experimentally observed slow dynamics effect cannot be explained by thermoelastic effects alone. The slow dynamics effect is to a great extent related to metastable states of the defect structure the latter passes through to due to the force action (dynamic or static) applied to the sample. After its removal, the defect structure slowly relaxes to its equilibrium state.     

Effect of thermal annealing on nanoimprinted Cu-Ni alloys using molecular dynamics simulation

The mechanical behaviors of nanoimprinted Cu-Ni alloys before and after annealing were studied using molecular dynamics simulation with a tight-binding potential. The results showed that when the punch is advancing, the punching force obtained from the simulation with a tight-binding potential is lower than with the Morse potential. During and after withdrawing the punch from the specimen, the adhesive phenomena are observed and the large residual stress in the Cu-Ni alloys is induced. During the annealing process, the internal energy of Cu-Ni alloys decreased with increasing the temperature and the component of Cu. In addition, comparing the maximum residual stress in the Cu-Ni alloys with and without annealing treatment, the stress is significantly released after annealing, especially in the higher component of Ni.     

Nonequilibrium molecular dynamics simulation of coupling between nanoparticles and base-fluid in a nanofluid

The intent of this study is to examine nonequilibrium heat transfer in a copper-argon nanofluid by molecular dynamics simulation. Two different methods, the physical definition method and the curve fitting method, are introduced to calculate the coupling factor between nanoparticles and base fluid. The results show that the coupling factors obtained by these two methods are consistent. The coupling factor is proportional to the volume fraction of the nanoparticle and inversely proportional to nanoparticle diameter. In the temperature range of 90-200 K, the coupling factor is not affected by temperature. The nanoparticle aggregation results in a decrease of the coupling factor.     

On the effect of an anisotropy-resolving subgrid-scale model on large eddy simulation predictions of turbulent open channel flow with wall roughness

A large eddy simulation (LES) was conducted of turbulent flow in a channel with a rough wall on one side and a free surface on the other by adopting an anisotropy-resolving subgrid-scale (SGS) model. A shear Reynolds number of Re_τ = 395 was used based on the mean friction velocity and channel height. To investigate the grid dependency of the LES results caused by the SGS model, three grid resolutions were tested under the same definition of a roughness shape by using the immersed boundary method. The results obtained were compared with direct numerical simulation data with and without the wall roughness and those without the extra anisotropic term. The primary focus was on how the present anisotropic SGS model with coarser grid resolutions can properly provide the effects of roughness on the mean velocity and turbulent stresses, leading to a considerable reduction of the computational cost of LES.     

接触角与液固界面热阻关系的分子动力学模拟

本文利用分子动力学方法模拟了液体在固体表面的 接触角及液固界面热阻, 并探讨了二者之间的关系. 通过分别改变液固结合强度和固体的原子性质来分析接触角和界面热阻的关系及变化趋势. 模拟结果显示增强液固间相互作用时, 接触角减小的同时界面热阻也随之单调减小; 而改变固体原子间结合强度和原子质量时, 接触角几乎保持不变, 但界面热阻显著改变. 固体原子间结合强度和原子质量影响界面热阻的原因是其改变了固体的振动频率分布, 导致液固原子间的振动耦合程度发生变化. 本文的结果表明界面热阻不仅与由接触角所表征的液固结合强度有关, 还与液固原子间的振动耦合程度有关. 接触角与界面热阻间不存在单值的对应关系, 不能单一地将接触角作为液固界面热阻的评价标准. 关键词： 液固界面 接触角 界面热阻 分子动力学模拟     

Coupling effect of Brownian motion and laminar shear flow on colloid coagulation:a Brownian dynamics simulation study

Simultaneous orthokinetic and perikinetic coagulations(SOPCs) are studied for small and large Peclet numbers(P e) using Brownian dynamics simulation.The results demonstrate that the contributions of the Brownian motion and the shear flow to the overall coagulation rate are basically not additive.At the early stages of coagulation with small Peclet numbers,the ratio of overall coagulation rate to the rate of pure perikinetic coagulation is proportional to P 1/2 e,while with high Peclet numbers,the ratio of overall coagulation rate to the rate of pure orthokinetic coagulation is proportional to P 1/2 e.Moreover,our results show that the aggregation rate generally changes with time for the SOPC,which is different from that for pure perikinetic and pure orthokinetic coagulations.By comparing the SOPC with pure perikinetic and pure orthokinetic coagulations,we show that the redistribution of particles due to Brownian motion can play a very important role in the SOPC.In addition,the effects of redistribution in the directions perpendicular and parallel to the shear flow direction are different.This perspective explains the behavior of coagulation due to the joint effects of the Brownian motion(perikinetic) and the fluid motion(orthokinetic).