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     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

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语言进行程序设计，定量模拟狭义相结论的运动时钟变慢，并对处于两种不同时空观的时钟运行进行对比，使得时间膨胀效应变得直观。     

Numerical Analysis and Simulation Analysis for Space‐Time Data

Spatio temporal dynamics of the positive column of a dc neon glow discharge is studied and investigated experimentally and theoretically. Spatio temporal analysis by means of biorthogonal decomposition method (BOD) gives insights into the mechanism of irregularity and can be employed for characterization of spatio‐ temporal complexity. In the weak nonlinear region, the wave dynamics is approximated by an amplitude equation of the Ginzburg‐Landau equation (CGLE) with complex coefficients and an additional integral term based on a fluid model. In the present work we deal with irregular spatio‐temporal data. A comparison between the numerical analysis of the experimental data and simulation results are studied. A good agreement between the dynamical behaviour for experimental space‐time data and theoretical simulation space‐time results was obtained. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Meshless simulation of equilibrium swelling/deswelling of pH‐sensitive hydrogels

Hydrogels have been widely used in microelectromechanical systems (MEMS) and Bio‐MEMS devices. In this article, the equilibrium swelling/deswelling of the pH‐stimulus cylindrical hydrogel in the microchannel is studied and simulated by the meshless method. The multi‐field coupling model, called multi‐effect‐coupling pH‐stimulus (MECpH) model, is presented and used to describe the chemical field, electric field, and the mechanical field involved in the problem. The partial differential equations (PDEs) describing these three fields are either nonlinear or coupled together. This multi‐field coupling and high nonlinear characteristics produce difficulties for the conventional numerical methods (e.g., the finite element method or the finite difference method), so an alternative—meshless method is developed to discretize the PDEs, and the efficient iteration technique is adopted to solve the nonlinear problem. The computational results for the swelling/deswelling diameter of the hydrogel under the different pH values are firstly compared with experimental results, and they have a good agreement. The influences of other parameters on the mechanical properties of the hydrogel are also investigated in detail. It is shown that the multi‐field coupling model and the developed meshless method are efficient, stable, and accurate for simulation of the properties of the stimuli‐sensitive hydrogel. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 326–337, 2006