Molecular dynamics simulation of observed c(4x4) and c(4x3) C60 alignments on the Si(100) reconstructed surface.

We have implemented a large-scale classical molecular dynamics simulation at constant temperature to provide a theoretical insight into the results of a recently performed experiment on the monolayer and multi-layer formations of molecular films on the Si(100) reconstructed dimerized surface. Our simulation has successfully reproduced all of the morphologies observed on the monolayer film by this experiment. We have obtained the formation of both c(4 4) and c(4 3) structures of the molecules and have also obtained phase transitions of the former into the latter.

Molecular dynamics simulation of adsorption of Ag particles on a graphite substrate

We have performed for the first time a molecular dynamics simulation of the adsorption of gas-phase Ag particles on a graphite substrate to provide an insight into the results of a comprehensive STM-based experiment on this system. Both pair-wise and many-body interatomic potentials have been employed, and a Morse-type Ag–C potential was specifically constructed to describe the interactions at the interface. Our simulation has successfully reproduced a significant portion of the experimental findings. We have also observed the intercalation of silver in graphite.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a twodimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a two dimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Numerical modelling of adsorption of metallic particles on graphite substrate via molecular dynamics simulation

A computer-based numerical modelling of the adsorption process of gas phase metallic particles on the surface of a graphite substrate has been performed via the application of molecular dynamics simulation method. The simulation relates to an extensive STM-based experiment performed in this field, and reproduces part of the experimental results. Both two-body and many-body inter-atomic potentials have been employed. A Morse-type potential describing the metal-carbon interactions at the interface was specifically formulated for this modelling. Intercalation of silver in graphite has been observed as well as the correct alignments of monomers, dimers and two-dimensional islands on the surface. PACS numbers: 02.60.Cb, 07.05.Tp, 68.55.-a, 81.05.Tp

Multiscale numerical modelling of crack propagation in two-dimensional metal plate

A novel multiscale model of brittle crack propagation in an Ag plate with macroscopic dimensions has been developed. The model represents crack propagation as stochastic drift-diffusion motion of the crack tip atom through the material, and couples the dynamics across three different length scales. It integrates the nanomechanics of bond rupture at the crack tip with the displacement and stress field equations of continuum based fracture theories. The finite element method is employed to obtain the continuum based displacement and stress fields over the macroscopic plate, and these are then used to drive the crack tip forward at the atomic level using the molecular dynamics simulation method based on many-body interatomic potentials. The linkage from the nanoscopic scale back to the macroscopic scale is established via the Ito stochastic calculus, the stochastic differential equation of which advances the tip to a new position on the macroscopic scale using the crack velocity and diffusion constant obtained on the nanoscale. Well known crack characteristics, such as the roughening transitions of the crack surfaces, crack velocity oscillations, as well as the macroscopic crack trajectories, are obtained.

Modelling the adsorption and imaging of C60 molecules on a graphite substrate

The adsorption of a C60 monolayer on a graphite substrate was modelled via molecular dynamics simulation covering a significant period of 160 picoseconds. The final configuration of C60s agrees closely with that observed in a scanning tunnelling microscopy (STM) experiment. Clusters of adsorbed molecules were then selected and their STM-like images were computed via the Keldysh Green function method.

Dynamics of an axisymmetric electromagnetic ‘crucible’ melting

Different industrial induction melting processes involve free surface and melt-solid interface of the liquid metal subject to dynamic change during the technological operation. Simulation of the liquid metal dynamics requires to solve the non-linear, coupled hydrodynamic-electromagnetic-heat transfer problem accounting for the time development of the liquid metal free boundary with a suitable turbulent viscosity model. The present paper describes a numerical solution method applicable for various axisymmetric induction melting processes, such as, crucible with free top surface, levitation, semi-levitation, cold crucible and similar melting techniques. The presented results in the cases of semi-levitation and crucible with free top surface meltings demonstrate oscillating transient behaviour of the free metal surface indicating the presence of gravity-inertial-electromagnetic waves which are coupled to the internal fluid flow generated by both the rotational and potential parts of the electromagnetic force.

The numerical simulation of enclosure fires using a Cfd fire field model coupled with a pyrolysis based solid fuel combustion submodel—a first approximation

In this paper, we present some early work concerned with the development of a simple solid fuel combustion model incorporated within a Computational Fluid Dynamics (CFD) framework. The model is intended for use in engineering applications of fire field modeling and represents an extension of this technique to situations involving the combustion of solid cellulosic fuels. A simple solid fuel combustion model consisting of a thermal pyrolysis model, a six flux radiation model and an eddy-dissipation model for gaseous combustion have been developed and implemented within the CFD code CFDS-FLOW3D. The model is briefly described and demonstrated through two applications involving fire spread in a compartment with a plywood lined ceiling. The two scenarios considered involve a fire in an open and closed compartment. The model is shown to be able to qualitatively predict behaviors similar to "flashover"—in the case of the open room—and "backdraft"— in the case of the initially closed room.

Surface physics at the nano-scale via scanning probe microscopy and molecular dynamics simulations

Probe-based scanning microscopes, such as the STM and the AFM, are used to obtain the topographical and electronic structure maps of material surfaces, and to modify their morphologies on nanoscopic scales. They have generated new areas of research in condensed matter physics and materials science. We will review some examples from the fields of experimental nano-mechanics, nano-electronics and nano-magnetism. These now form the basis of the emerging field of Nano-technology. A parallel development has been brought about in the field of Computational Nano-science, using quantum-mechanical techniques and computer-based numerical modelling, such as the Molecular Dynamics (MD) simulation method. We will report on the simulation of nucleation and growth of nano-phase films on supporting substrates. Furthermore, a theoretical modelling of the formation of STM images of metallic clusters on metallic substrates will also be discussed within the non-equilibrium Keldysh Green function method to study the effects of coherent tunnelling through different atomic orbitals in a tip-sample geometry.

Liquid metal turbulent flow dynamics in a cylindrical container with free surface: experiment and numerical analysis

The liquid metal flow in induction crucible models is known to be unstable, turbulent and difficult to predict in the regime of medium frequencies when the electromagnetic skin-layer is of considerable extent. We present long term turbulent flow measurements by a permanent magnet incorporated potential difference velocity probe in a cylindrical container filled with eutectic melt In-Ga-Sn. The parallel numerical simulation of the long time scale development of the turbulent average flow is presented. The numerical flow model uses an implicit pseudo-spectral code and k-w turbulence model, which was recently developed for the transitional flow modelling. The results compare reasonably to the experiment and demonstrate the time development of the turbulent flow field and the turbulence energy.

On the coupling of Navier–Stokes and linearised Euler equations for aeroacoustic simulation

Aerodynamic generation of sound is governed by the Navier–Stokes equations while acoustic propagation in a non-uniform medium is effectively described by the linearised Euler equations. Different numerical schemes are required for the efficient solution of these two sets of equations, and therefore, coupling techniques become an essential issue. Two types of one-way coupling between the flow solver and the acoustic solver are discussed: (a) for aerodynamic sound generated at solid surfaces, and (b) in the free stream. Test results indicate how the coupling achieves the necessary accuracy so that Computational Fluid Dynamics codes can be used in aeroacoustic simulations.

The numerical simulation of noncharring thermal degradation and its application to the prediction of compartment fire development

In this paper we present some work concerned with the development and testing of a simple solid fuel combustion model incorporated within a Computational Fluid Dynamics (CFD) framework. The model is intended for use in engineering applications of fire field modeling and represents an extension of this technique to situations involving the combustion of solid fuels. The CFD model is coupled with a simple thermal pyrolysis model for combustible solid noncharring fuels, a six-flux radiation model and an eddy-dissipation model for gaseous combustion. The model is then used to simulate a series of small-scale room fire experiments in which the target solid fuel is polymethylmethacrylate. The numerical predictions produced by this coupled model are found to be in very good agreement with experimental data. Furthermore, numerical predictions of the relationship between the air entrained into the fire compartment and the ventilation factor produce a characteristic linear correlation with constant of proportionality 0.38 kg/sm5/12. The simulation results also suggest that the model is capable of predicting the onset of "flashover" type behavior within the fire compartment.

Liquid metal turbulent flow dynamics in a cylindrical container with free surface: experiment and numerical analysis

The liquid metal flow in inducation crucible models is known to be higly unstable and turbutlen in the regim e of medium frequecies when the elctronmagnetic skin-layer is of considerable extent. We present long term turbulent flow measurements by a permanent magnet incorporated potential difference veolocity probe in a cylindirical container filled with eutecti mlt In-Ga-SN. The parallel numerical simulation of the long time scale development of the turbulen average flow is presented. The numerical lfow model uses a pseud-spectral code and k-w turbulence model, which was recently developed for the transitional flow modelling. The result compare reasonably to the experiment and demonstrate the time development of the turbulent flow field.

Utilising computational fluid dynamics (CFD) for the modelling of granular material in large-scale engineering processes

In this paper, the framework is described for the modelling of granular material by employing Computational Fluid Dynamics (CFD). This is achieved through the use and implementation in the continuum theory of constitutive relations, which are derived in a granular dynamics framework and parametrise particle interactions that occur at the micro-scale level. The simulation of a process often met in bulk solids handling industrial plants involving granular matter, (i.e. filling of a flat-bottomed bin with a binary material mixture through pneumatic conveying-emptying of the bin in core flow mode-pneumatic conveying of the material coming out of a the bin) is presented. The results of the presented simulation demonstrate the capability of the numerical model to represent successfully key granular processes (i.e. segregation/degradation), the prediction of which is of great importance in the process engineering industry.