热对流条件下颗粒沉降的格子Boltzmann方法模拟

采用格子Boltzmann方法模拟了在热对流条件下的颗粒沉降问题, 在研究单颗粒在等温流体、热流体和冷流体中运动的基础上, 进一步模拟了两个不同温度的颗粒在流体中的沉降.结果表明:两等温颗粒的沉降方式与雷诺数Re以及格拉晓夫数Gr密切相关, 而两不同温度的颗粒与两等温颗粒的沉降规律有显著不同.无论初始位置如何, 冷颗粒最终总位于热颗粒下方运动, Re较大时, 发生连续的拖曳、接触现象, 而Re较小时, 冷颗粒会以较大的沉降速度远离热颗粒. 关键词： 格子Boltzmann方法 颗粒沉降 热对流     

Structural characterization of pseudo-binary semiconducting alloys using molecular simulations

Monte Carlo simulations with the Keating model have been performed to predict the lattice constant and bond length variations with composition for pseudo-binary semiconducting alloys. In general, it is observed that the deviations of the lattice constants from Vegard's law predictions are larger as the lattice mismatch between the constituent binaries increases. Further, it is noted that these alloys have partial virtual crystal model characteristics and tend to be more towards the flexible (floppy) crystal limit as compared to the rigid crystal limit. The topological rigidity parameters are bond-type dependent. The angular deviations from perfect tetrahedral structure are also measured.     

Entropic lattice Boltzmann method for turbulent flow simulations: Boundary conditions

We introduce a new concept of boundary conditions for realization of the lattice Boltzmann simulations of turbulent flows. The key innovation is the use of a universal distribution function for particles, analogous to the Tamm–Mott-Smith solution for the shock wave in the classical Boltzmann kinetic equation. Turbulent channel flow simulations demonstrate that the new boundary enables accurate results even with severely under-resolved grids. Generalization to complex boundary is illustrated with an example of turbulent flow past a circular cylinder.     

Non-local order in Mott insulators,duality and Wilson loops

It is shown that the Mott insulating and superfluid phases of bosons in an optical lattice may be distinguished by a non-local ‘parity order parameter’ which is directly accessible via single site resolution imaging. In one dimension, the lattice Bose model is dual to a classical interface roughening problem. We use known exact results from the latter to prove that the parity order parameter exhibits long range order in the Mott insulating phase, consistent with recent experiments by Endres et al. [M. Endres, M. Cheneau, T. Fukuhara, C. Weitenberg, P. Schauß, C. Gross, L. Mazza, M.C. Bañuls, L. Pollet, I. Bloch, et al., Science 334 (2011) 200]. In two spatial dimensions, the parity order parameter can be expressed in terms of an equal time Wilson loop of a non-trivial U(1)$U\left(1\right)$ gauge theory in 2+1$2+1$ dimensions which exhibits a transition between a Coulomb and a confining phase. The negative logarithm of the parity order parameter obeys a perimeter law in the Mott insulator and is enhanced by a logarithmic factor in the superfluid.     

Etude de la restauration de l'aluminium irradié aux electrons (Etape III)

The nature of the electron bombarded aluminium electrical recovery process has been studied by combining quenches with irradiations. Measured Frenkel pair formation rates per incident electrons as well as annealing kinetics emphasize the role of impurity interstitial trapping.     

Harmonics in the stress response of a WCA-Lennard-Jones liquid under oscillatory elongational flow

The WCA fluid shows a nonlinear behaviour toward a periodic elongational flow through the presence of higher harmonics in stress responses. At the same time the liquid remains purely viscous and isotropic. The absence of numerical artefacts is checked and the physical origin of the harmonics is discussed.     

Lattice-dynamical functions of solids at very high pressures

Abstract

Vibrational spectra of face-centered cubical (f.c.c.) solids are investigated at the pressure region, where the quantum-statistical model for the electron energy is applicable. Phonon dispersion curves are obtained for different values of the specific volume. The integration of the phonon spectra over the Brillouin zone yields the zero-point vibrational energy, the Debye temperature, the Gruneisen coefficient, mean-square amplitude of zero-point and thermal vibration, the melting temperature. All quantities under consideration are calculated from the universal functions of the reduced volume, using the scaling relations for given atomic number 2.     

Dynamics of surface dipole and tripole solitons in nonlocal nonlinear media with optical lattice field

We find the existence conditions for stationary dipole and tripole surface solitons formed at the interface of a nonlocal nonlinear medium and a lattice with linearly modulated frequency. We investigate how the degree of nonlocality, the depth, and the modulation frequency of the optical lattice field affect on the existence of the surface solitons and their dynamics. The relationship between the power and the model parameters is identified. The stability of the surface dipole and tripole solitons is numerically investigated.     

Structure of precipitates and orientation relationships in Al,Ge, Mo-doped higher manganese silicide crystals

Abstract

The structure of Al, Ge, Mo-doped Higher Manganese Silicide (HMS) crystals with the general formulas Mn(Si_0.99Ge_0.01)_1.75, Mn(Si_0.995Ge_0.005)_1.75 and (Mn_0.98Mo_0.02)[(Si_0.98Ge_0.02)_1.75]_0.99Al_0.01 was investigated by scanning and transmission electron microscopy, electron diffraction and X-ray energy dispersive spectrometry in a wide scale range from a few mm to several Å. Several secondary phases were identified in the Mn₄Si₇ matrix: Ge_1?xSi_x (0.1 < x < 0.9) solid solution precipitates with Ge concentration ranging from 5 at. % up to 93 at.%, MoSi₂ platelets, MnSi and Mn₅Si₃ precipitates. Their morphology, structure and crystallographic relationships with the HMS matrix were determined. Mostly local strains in the matrix and precipitates due to lattice misfits at interfaces derived from crystallographic relationships were found two orders of magnitude higher than deformation induced by thermal expansion mismatch. Only a few exceptions of specific relationships were found when the lattice misfit and thermal mismatch have close values. The largest misfit of about 22% was observed between MnSi and Mn₄Si₇ what led to big and numerous cracks in crystals. Therefore, doping can improve the material performance (1) by preventing the formation of MnSi precipitates with metallic properties and (2) by reduction of cracking and crack propagation because of larger MnSi /Mn₄Si₇ lattice misfit compared to Ge_1?xSi_x /Mn₄Si₇ or MoSi₂/Mn₄Si₇ misfits.     

Uniform random number generation

In typical stochastic simulations, randomness is produced by generating a sequence of independent uniform variates (usually real-valued between 0 and 1, or integer-valued in some interval) and transforming them in an appropriate way. In this paper, we examine practical ways of generating (deterministic approximations to) such uniform variates on a computer. We compare them in terms of ease of implementation, efficiency, theoretical support, and statistical robustness. We look in particular at several classes of generators, such as linear congruential, multiple recursive, digital multistep, Tausworthe, lagged-Fibonacci, generalized feedback shift register, matrix, linear congruential over fields of formal series, and combined generators, and show how all of them can be analyzed in terms of their lattice structure. We also mention other classes of generators, like non-linear generators, discuss other kinds of theoretical and empirical statistical tests, and give a bibliographic survey of recent papers on the subject.