Structure of solid monolayers and multilayers of n -hexane on graphite

We present all-atom molecular dynamics simulations ofn-hexane on the basal plane of graphite at monolayer and multilayer coverages. In keeping with experimental data, we find the presence of ordered adsorbed layers both at single monolayer coverage and when the adsorbed layer coexists with excess liquid adsorbate. Using a simulation method that does not impose any particular periodicity on the adsorbed layer, we quantitatively compare our results to the results of neutron diffraction experiments and find a structural transition from a uniaxially incommensurate lattice to a fully commensurate structure on increasing the coverage from a monolayer to a multilayer. The zig-zag backbone planes of all the alkane molecules lie parallel to the graphite surface at the multilayer coverage, while a few molecules are observed to attain the perpendicular orientation at monolayer coverage. Dedicated to Professor C N R Rao on his 70th birthday     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

Theoretical and numerical analysis on multispectral bioluminescence tomography

Recently, molecular imaging has been rapidly developed to studyphysiological and pathological processes in vivo at the cellularand molecular levels. Among molecular imaging modalities, opticalimaging has attracted a major attention for its unique advantages.In this paper, we establish a mathematical framework for multispectralbioluminescence tomography (BLT) that allows simultaneous studiesof multiple optical reporters. We show solution existence, uniquenessand continuous dependence on data as well as the limiting behaviourswhen the regularization parameter approaches zero or when thepenalty parameter approaches infinity. Then, we propose twonumerical schemes for multispectral BLT and derive error estimatesfor the corresponding solutions.     

A RAPID APPROXIMATION FOR PREDICTING THE HYDROCARBON DISCOVERY RATE: PART 1. ASSESSING THE ACCURACY

The hydrocarbon discovery prediction problem is important to firms having to make decisions about the deployment of scarce exploration resources. Traditional methods for estimating the discovery rate rely on the completion of time consuming simulation experiments. A rapid approximation that does not require the completion of simulation exists and has been shown to have some promise as a prediction tool. This paper investigates the accuracy of the approximation method under a wide variety of distributional and drilling efficiency assumptions. The results indicate that the approximation produces predictions close to those of simulation under most of the tested conditions. This suggests that resource exploration firms could conveniently use the method for a wide variety of planning purposes without incurring the same costs in time and personnel required for simulation.     

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

Measurements and simulations of particle dispersion in a turbulent flow

A particle imaging technique has been used to collect droplet displacement statistics in a round turbulent jet of air. Droplets are injected on the jet axis, and a laser sheet and position-sensitive photomultiplier tube are used to track their radial displacement and time-of-flight. Dispersion statistics can be computed which are Lagrangian or Eulerian in nature. The experiments have been simulated numerically using a second-order closure scheme for the jet and a stochastic simulation for the particle trajectories. Results are presented for non-vaporizing droplets of sizes from 35 to 160 μm. The simulations have underscored the importance of initial conditions and early droplet displacement history on the droplet trajectory for droplets with large inertia relative to the turbulence. Estimates of initial conditions have been made and their effect on dispersion is quantified.     

狭义相对论效应──运动时钟变慢的计算机模拟演示实验

利用精炼的C语言进行程序设计，定量模拟狭义相结论的运动时钟变慢，并对处于两种不同时空观的时钟运行进行对比，使得时间膨胀效应变得直观。     

Stationary states of crystal growth in three dimensions

A new Markov process describing crystal growth in three dimensions is introduced. States of the process are configurations of the crystal surface, which has a terrace-edge-kink structure. The states are continuous along edges but discrete across edges, in accordance with the very different rates for the two types of captures of particles. Stationary distributions, describing steady crystal growth, are found in general. To our knowledge, these are the first examples of stationary distributions for layered crystal growth in three dimensions. The steady growth rate and other quantities are obtained explicitly for two interacting edges. For many interacting edges, growth behavior is determined (a) in various asymptotic regimes including thermodynamic limits, (b) via simulations, and (c) using series (cluster) expansions in the slope of the surface, the first three coefficients being computed. The theoretical growth rates show a marked dependence on surface dimensions. This may contribute to the size dependence and dispersion in the observed growth rate of small crystals.     

Short Communication: Application of WENO (Weighted Essentially Non-oscillatory) Approach to Navier–Stokes Computations

The direct implementation of the essentially non-oscillatory schemes for flow simulation over complex geometries sometimes results in insufficiently robust numerical algorithms. In order to overcome this difficulty, it is suggested to use the weighted essentially non-oscillatory approach for multidimensional Navier–Stokes computations. The results indicate a significant improvement in accuracy and robustness, especially for low Mach and high supersonic flows.     

Solid deformation induced by the adsorption of methane and methanol under sub- and supercritical conditions

Carbonaceous materials with some degree of flexibility in their physical structure can expand or contract under the influence of the forces exerted by adsorbed molecules. To gain insight into how adsorption of non-polar and polar fluids could deform a carbon solid, we present GCMC simulations of sub- and supercritical adsorption of methane and methanol in slit-shaped pores whose walls are made of graphene layers. Our extensive simulation study shows that there is a strong correlation between solvation pressure and solid deformation, and that the expansion or contraction of the pore strongly depends on adsorbate loading, temperature and pore size.