Evaporation of Lennard-Jones fluids

Evaporation and condensation at a liquid/vapor interface are ubiquitous interphase mass and energy transfer phenomena that are still not well understood. We have carried out large scale molecular dynamics simulations of Lennard-Jones (LJ) fluids composed of monomers, dimers, or trimers to investigate these processes with molecular detail. For LJ monomers in contact with a vacuum, the evaporation rate is found to be very high with significant evaporative cooling and an accompanying density gradient in the liquid domain near the liquid/vapor interface. Increasing the chain length to just dimers significantly reduces the evaporation rate. We confirm that mechanical equilibrium plays a key role in determining the evaporation rate and the density and temperature profiles across the liquid/vapor interface. The velocity distributions of evaporated molecules and the evaporation and condensation coefficients are measured and compared to the predictions of an existing model based on kinetic theory of gases. Our results indicate that for both monatomic and polyatomic molecules, the evaporation and condensation coefficients are equal when systems are not far from equilibrium and smaller than one, and decrease with increasing temperature. For the same reduced temperature T/T(c), where T(c) is the critical temperature, these two coefficients are higher for LJ dimers and trimers than for monomers, in contrast to the traditional viewpoint that they are close to unity for monatomic molecules and decrease for polyatomic molecules. Furthermore, data for the two coefficients collapse onto a master curve when plotted against a translational length ratio between the liquid and vapor phase.     

Cover Image

A hybrid computational method coupling the lattice-Boltzmann (LB) method and a Langevin-dynamics (LD) method is developed to simulate nanoscale particle and polymer (NPP) suspensions in the presence of both thermal fluctuation and long-range many-body hydrodynamic interactions (HIs). Brownian motion of the NPP is explicitly captured by a stochastic forcing term in the LD method. The LD method is two-way coupled to the nonfluctuating LB fluid through a discrete LB forcing source distribution to capture the long-range HI. To ensure intrinsically linear scalability with respect to the number of particles, a Eulerian-host algorithm for short-distance particle neighbor search and interaction is developed and embedded to LB-LD framework. The validity and accuracy of the LB-LD approach are demonstrated through several sample problems. The simulation results show good agreements with theory and experiment. The LB-LD approach can be favorably incorporated into complex multiscale computational frameworks for efficiently simulating multiscale multicomponent particulate suspension systems such as complex blood suspensions.     

A conformal decomposition finite element method for modeling stationary fluid interface problems

A method is developed for modeling fluid transport in domains that do not conform to the finite element mesh. One or more level set functions are used to describe the fluid domain. A background, non‐conformal mesh is decomposed into elements that conform to the level set interfaces. Enrichment takes place by adding nodes that lie on the interfaces. Unlike other enriched finite element methods, the proposed technique requires no changes to the underlying element assembly, element interpolation, or element quadrature. The complexity is entirely contained within the element decomposition routines. It is argued that the accuracy of the method is no less than that for eXtended Finite Element Methods (XFEM) with Heaviside enrichment. The accuracy is demonstrated using multiple numerical tests. In all cases, optimal rates of convergence are obtained for both volume and surface quantities. Jacobi preconditioning is shown to remove the ill‐conditioning that may result from the nearly degenerate conformal elements. Copyright © 2009 John Wiley & Sons, Ltd.     

Three-dimensional shear in granular flow

The evolution of granular shear flow is investigated as a function of height in a split-bottom Couette cell. Using particle tracking, magnetic-resonance imaging, and large-scale simulations, we find a transition in the nature of the shear as a characteristic height H* is exceeded. Below H* there is a central stationary core; above H* we observe the onset of additional axial shear associated with torsional failure. Radial and axial shear profiles are qualitatively different: the radial extent is wide and increases with height, while the axial width remains narrow and fixed.     

Simple Zeros of the Riemann Zeta-Function

Assuming the Riemann Hypothesis, Montgomery showed by meansof his pair correlation method that at least two-thirds of thezeros of Riemann's zeta-function are simple. Later he and Taylorimproved this to 67.25 percent and, more recently, Cheer andGoldston increased the percentage to 67.2753. Here we proveby a new method that if the Riemann and Generalized LindelöofHypotheses hold, then at least 70.3704 percent of the zerosare simple and at least 84.5679 percent are distinct. Our methoduses mean value estimates for various functions defined by Dirichletseries sampled at the zeros of the Riemann zeta-function. 1991Mathematics Subject Classification: 11M26.     

Existence and Uniqueness of Fixed Points for Markov Operators and Markov Processes

This paper concerns a Markov operator T on a space L₁, and aMarkov process P which defines a Markov operator on a spaceM of finite signed measures. For T, the paper presents necessaryand sufficient conditions for:

a the existence of invariant probabilitydensities (IPDs)

b the existence of strictly positive IPDs,and

c the existence and uniqueness of IPDs.

Similar resultson invariant probability measures for P are presented. The basicapproach is to pose a fixed-point problem as the problem ofsolving a certain linear equation in a suitable Banach space,and then obtain necessary and sufficient conditions for thisequation to have a solution. 1991 Mathematics Subject Classification:60J05, 47B65, 47N30.     

库仑加线性位基态解析解

应用新发展的单一轨迹积分方法求解库仑加线性位的基态量子波函数,得到基态能量和波函数的一般解析表达式,并讨论了解的收敛性.应用此方法讨论了重夸克偶素系统.