Molecular Dynamics Simulation of Microstructure of Nanocrystalline Copper

The microstructure of computer generated nanocrystalline coppers is simulated by using molecular dynamics with the Finnis-Sinclair potential, analysed by means of radial distribution functions, coordination number, atomic energy and local crystalline order. The influence of the grain size on the nanocrystalline structure is studied. The results reveal that as the grain size is reduced, the grain boundary shows no significant structural difference, but the grain interior becomes more disordered, and their structural difference diminishes gradually; however, the density and the atomic average energy of the grain boundary present different tendencies from those of the grain interior.

Brownian dynamics simulation of the crystallization dynamics of charged colloidal particles

Crystal formation process of charged colloidal particles is investigated using Brownian dynamics (BD) simulations. The particles are assumed to interact with the pair-additive repulsive Yukawa potential. The time evolution of crystallization process and the crystal structure during the simulation are characterized by means of the radial distribution functions (RDF) and mean square displacement (MSD). The simulations show that when the interaction is featured with long-range, particles can spontaneously assemble into body-centered-cubic (BCC) arrays at relatively low particle number density. When the interaction is short-ranged, with increasing the number density particles become trapped into a stagnant disordered configuration before the crystallization could be actualized. The simulations further show that as long as the trapped configurations are bypassed, the face-centered-cubic (FCC) structures can be achieved and are actually more stable than BCC structures. (C) 2010 Elsevier Inc. All rights reserved.

Molecular dynamics study on the phase diagrams of linear and branched chain molecules

The particle transfer molecular dynamics is used to study the phase equilibria of linear and branched chain molecules. The scaling of the critical temperature versus chain length is obtained and the critical densities are found to decrease with increasing chain length, which are in agreement with the results of experiment and theory. The phase diagrams of the linear and the branched chain molecules nearly overlap with each other. Moreover, the radial distribution functions of linear and branched chain molecules in gas phase are very similar, but in the liquid phase, they are different for different kinds of chains.

Spatial distribution of soil heavy metal pollution estimated by different interpolation methods: Accuracy and uncertainty analysis

Mapping the spatial distribution of contaminants in soils is the basis of pollution evaluation and risk control. Interpolation methods are extensively applied in the mapping processes to estimate the heavy metal concentrations at unsampled sites. The performances of interpolation methods (inverse distance weighting, local polynomial, ordinary kriging and radial basis functions) were assessed and compared using the root mean square error for cross validation. The results indicated that all interpolation methods provided a high prediction accuracy of the mean concentration of soil heavy metals. However, the classic method based on percentages of polluted samples, gave a pollution area 23.54-41.92% larger than that estimated by interpolation methods. The difference in contaminated area estimation among the four methods reached 6.14%. According to the interpolation results, the spatial uncertainty of polluted areas was mainly located in three types of region: (a) the local maxima concentration region surrounded by low concentration (clean) sites, (b) the local minima concentration region surrounded with highly polluted samples; and (c) the boundaries of the contaminated areas. (C) 2010 Elsevier Ltd. All rights reserved.

Extended bidirectional reflectance distribution function for polarized light scattering from subsurface defects under a smooth surface

By introducing the scattering probability of a subsurface defect (SSD) and statistical distribution functions of SSD radius, refractive index, and position, we derive an extended bidirectional reflectance distribution function (BRDF) from the Jones scattering matrix. This function is applicable to the calculation for comparison with measurement of polarized light-scattering resulting from a SSD. A numerical calculation of the extended BRDF for the case of p-polarized incident light was performed by means of the Monte Carlo method. Our numerical results indicate that the extended BRDF strongly depends on the light incidence angle, the light scattering angle, and the out-of-plane azimuth angle. We observe a 180 degrees symmetry with respect to the azimuth angle. We further investigate the influence of the SSD density, the substrate refractive index, and the statistical distributions of the SSD radius and refractive index on the extended BRDF. For transparent substrates, we also find the dependence of the extended BRDF on the SSD positions. (c) 2006 Optical Society of America.

Lifetime distribution of spontaneous emission from line antennas in two-dimensional quasi-periodic photonic crystals

We investigate the lifetime distribution functions of spontaneous emission from line antennas embedded in finite-size two-dimensional 12-fold quasi-periodic photonic crystals. Our calculations indicate that two-dimensional quasi-periodic crystals lead to the coexistence of both accelerated and inhibited decay processes. The decay behaviors of line antennas are drastically changed as the locations of the antennas are varied from the center to the edge in quasi-periodic photonic crystals and the location of transition frequency is varied.

Spontaneous emission lifetime distribution of infinite line antennas in two-dimensional photonic crystals with finite size

The authors calculate the lifetime distribution functions of spontaneous emission from infinite line antennas embedded in two-dimensional disordered photonic crystals with finite size. The calculations indicate the coexistence of both accelerated and inhibited decay processes in disordered photonic crystals with finite size. The decay behavior of the spontaneous emission from infinite line antennas changes significantly by varying factors such as the line antennas' positions in the disordered photonic crystal, the shape of the crystal, the filling fraction, and the dielectric constant. Moreover, the authors analyze the effect of the degree of disorder on spontaneous emission. (c) 2007 American Institute of Physics.

Quark distributions in nuclei and nuclear effects on nucleon structure functions

Extended quark distribution functions are presented obtained by fitting a large amount of experimental data of the l-A DIS process on the basis of an improved nuclear density model. The experimental data of l-A DIS processes with A >= 3 in the region 0.0010 <= x <= 0.9500 axe quite satisfactorily described by using the extended formulae. Our knowledge of the influence of nuclear matter on the quark distributions is deepened.

铜-铝扩散焊及拉伸的分子动力学模拟

采用分子动力学方法模拟了铜-铝扩散焊过程,分析了理想平面铜-铝试件(001)晶面间扩散焊的过渡层厚度,并利用径向分布、键对分析方法分析了在不同的降温速率下过渡层的结构变化.降温速率大时,过渡层保持原有无序结构,降温速率小时,过渡层从无序结构向面心立方结构转变.还对扩散焊后的铜-铝试件进行了拉伸模拟,并与尺寸大小相近的单晶铜和单晶铝的拉伸模拟结果进行比较.结果发现焊接后的强度比单晶铝和单晶铜的强度都要小,最大应变值也小.

Microstructure and nanomechanical properties of the wing membrane of dragonfly

Detailed investigations on the microstructure and the mechanical properties of the wing membrane of the dragonfly were carried out. It was found that in the direction of the thickness the membrane was divided into three layers rather than as traditionally considered as a single entity, and on the surfaces the membrane displayed a random distribution rough microstructure that was composed of numerous nanometer scale columns coated by the cuticle wax secreted. The characteristics of the surfaces were accurately measured and a statistical radial distribution function of the columns was presented to describe the structural properties of the surfaces. Based on the surface microstructure, the mechanical properties of the membranes taken separately from the wings of living and dead dragonflies were investigated by the nanoindentation technique. The Young's moduli obtained here are approximately two times greater than the previous result, and the reasons that yield the difference are discussed. (C) 2007 Elsevier B.V. All rights reserved.

Study on the ions' behavior in an electron cyclotron resonance plasma

The energy, velocity, angle distribution of ions in magnetoactive electron cyclotron resonance plasma have been studied with a two-dimension hybrid mode. The dependence of these distribution functions versus position and pressure are discussed. Our simulation results are in good agreement with many experimental measurements. (C) 1997 American Institute of Physics.

双晶晶界在高温下的热稳定性的分子动力学研究

晶界结构在高温下的热稳定性问题是一个长期争论而又未能解决的问题，其争论的焦点是：在远低于熔点的温度下，晶界结构是否发生了可观察到的无序化，即是否存在一个远低于熔点的结构转化温度。为了能澄清这一争论，本文系统地研究了晶界结构的热稳定性。为了消除相互作用势的影响和系统误差，本文首先采用Morse势和经验多体势分别对铝、铜单晶的熔化过程进行了分子动力学模拟。在平衡态下，通过计算表征结构无序化的静态结构因子、径向分布函数和单晶原子位形图，获得了铝、铜单晶的熔点，结果表明：多体势计算的铝和铜的单晶熔点更接近实验值。因此，采用经验多体势应用分子动力学方法分别模拟了铝、铜Σ3、Σ5、Σ9、Σ11、Σ19、Σ33六种对称倾侧双晶晶界晶界结构由有序向无序转化的过程，计算了平衡态下的表征结构无序化的静态结构因子、径向分布函数和晶界原子位形图并将多体势获得的铝、铜单晶熔点作为晶界结构转化温度的约化熔点，获得了铝、铜Σ3、Σ5、Σ9、Σ11、Σ19、Σ33六种对称倾侧双晶晶界结构的转化温度和熔点，结果表明：1．Σ5、Σ9、Σ11、Σ19、Σ33五种对称倾侧双晶晶界均在远低于单晶熔点温度时，晶界结构发生了可观察到的无序化，而且双晶晶界结构的转变温度相差不大，双晶晶界熔点也低于单晶熔点。2．Σ3晶界在温度远低于熔点时，其晶界结构没有发生可观察到的无序化；Σ3晶界的转化温度与单晶熔点接近。所以，可以认为Σ3晶界不存在转化温度。这是由于Σ3晶界为共格孪晶，具有较低的能量。综上所述，除Σ3共格孪晶外，在远低于熔点温度下，晶界结构发生了可观察到的无序化，即：存在一个远低于熔点的转化温度，此时其静态结构因子约为0.5左右；晶界结构的熔点均低于单晶熔点，此时其静态结构因子约为0.15左右。从全文模拟结果可以看出，静态结构因子、径向分布函数、晶界原子位形图三种方法在确定晶界的结构转化温度和熔点时，静态结构因子是最有效、最准确的定量方法。

Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle- resolved simulations are briefly discussed.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.