Comparisons of diffusive and viscous contributions to transport coefficients of light gases in single-walled carbon nanotubes

We examine here the relative importance of different contributions to transport of light gases in single walled carbon nanotubes, using methane and hydrogen as examples. Transport coefficients at 298 K are determined using molecular dynamics simulation with atomistic models of the nanotube wall, from which the diffusive and viscous contributions are resolved using a recent approach that provides an explicit expression for the latter. We also exploit an exact theory for the transport of Lennard-Jones fluids at low density considering diffuse reflection at the tube wall, thereby permitting the estimation of Maxwell coefficients for the wall reflection. It is found that reflection from the carbon nanotube wall is nearly specular, as a result of which slip flow dominates, and the viscous contribution is small in comparison, even for a tube as large as 8.1 nm in diameter. The reflection coefficient for hydrogen is 3-6 times as large as that for methane in tubes of 1.36 nm diameter, indicating less specular reflection for hydrogen and greater sensitivity to atomic detail of the surface. This reconciles results showing that transport coefficients for hydrogen and methane, obtained in simulation, are comparable in tubes of this size. With increase in adsorbate density, the reflection coefficient increases, suggesting that adsorbate interactions near the wall serve to roughen the local potential energy landscape perceived by fluid molecules.

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy ECG simulation when it is used in heart torso model-based body surface potential simulation studies.

Interfacial interactions and structure of polyurethane intercalated nanocomposite

Understanding the interfacial interactions between the nanofiller and polymer matrix is important to improve the design and manufacture of polymer nanocomposites. This paper reports a molecular dynamic Study on the interfacial interactions and structure of a clay-based polyurethane intercalated nanocomposite. The results show that the intercalation of surfactant (i.e. dioctadecyldlmethyl ammonium) and polyurethane (PU) into the nanoconfined gallery of clay leads to the multilayer structure for both surfactant and PU, and the absence of phase separation for PU chains. Such structural characteristics are attributed to the result of competitive interactions among the surfactant, PU and the clay surface, including van der Waals, electrostatic and hydrogen bonding.

Quantum effects on adsorption and diffusion of hydrogen and deuterium in microporous materials

Monte Carlo and molecular dynamics simulations and neutron scattering experiments are used to study the adsorption and diffusion of hydrogen and deuterium in zeolite Rho in the temperature range of 30-150 K. In the molecular simulations, quantum effects are incorporated via the Feynman-Hibbs variational approach. We suggest a new set of potential parameters for hydrogen, which can be used when Feynman-Hibbs variational approach is used for quantum corrections. The dynamic properties obtained from molecular dynamics simulations are in excellent agreement with the experimental results and show significant quantum effects on the transport at very low temperature. The molecular dynamics simulation results show that the quantum effect is very sensitive to pore dimensions and under suitable conditions can lead to a reverse kinetic molecular sieving with deuterium diffusing faster than hydrogen.

The effect of box shape on the dynamic properties of proteins simulated under periodic boundary conditions

The effect of the box shape on the dynamic behavior of proteins simulated under periodic boundary conditions is evaluated. In particular, the influence of simulation boxes defined by the near-densest lattice packing (NDLP) in conjunction with rotational constraints is compared to that of standard box types without these constraints. Three different proteins of varying size, shape, and secondary structure content were examined in the study. The statistical significance of differences in RMSD, radius of gyration, solvent-accessible surface, number of hydrogen bonds, and secondary structure content between proteins, box types, and the application or not of rotational constraints has been assessed. Furthermore, the differences in the collective modes for each protein between different boxes and the application or not of rotational constraints have been examined. In total 105 simulations were performed, and the results compared using a three-way multivariate analysis of variance (MANOVA) for properties derived from the trajectories and a three-way univariate analysis of variance (ANOVA) for collective modes. It is shown that application of roto-translational constraints does not have a statistically significant effect on the results obtained from the different simulations. However, the choice of simulation box was found to have a small (5-10%), but statistically significant effect on the behavior of two of the three proteins included in the study. (c) 2005 Wiley Periodicals, Inc.

Adsorption of Supercritical fluids on graphitised thermal carbon black: Molecular layer structure theory versus grand canonical Monte Carlo simulation

This paper presents a detailed analysis of adsorption of supercritical fluids on nonporous graphitized thermal carbon black. Two methods are employed in the analysis. One is the molecular layer structure theory (MLST), proposed recently by our group, and the other is the grand canonical Monte Carlo (GCMC) simulation. They were applied to describe the adsorption of argon, krypton, methane, ethylene, and sulfur hexafluoride on graphitized thermal carbon black. It was found that the MLST describes all the experimental data at various temperatures well. Results from GCMC simulations describe well the data at low pressure but show some deviations at higher pressures for all the adsorbates tested. The question of negative surface excess is also discussed in this paper.

A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Correct procedures for the calculation of heats of adsorption for heterogeneous adsorbents from molecular simulation

Several procedures for calculating the heat of adsorption from Monte Carlo simulations for a heterogeneous adsorbent are presented. Simulations have been performed to generate isotherms for nitrogen at 77 K and methane at 273.15 K in graphitic slit pores of various widths. The procedures were then applied to calculate the heat of adsorption of an activated carbon with an arbitrary pore size distribution. The consistency of the different procedures shows them to be correct in calculating interaction energy contributions to the heat of adsorption. The currently favored procedure for this type of calculation, from the literature, is shown to be incorrect and in serious error when calculating the heat of adsorption of activated carbon.

Adsorption in mesopores: A molecular-continuum model with application to MCM-41

The classical model of surface layering followed by capillary condensation during adsorption in mesopores, is modified here by consideration of the adsorbate solid interaction potential. The new theory accurately predicts the capillary coexistence curve as well as pore criticality, matching that predicted by density functional theory. The model also satisfactorily predicts the isotherm for nitrogen adsorption at 77.4 K on MCM-41 material of various pore sizes, synthesized and characterized in our laboratory, including the multilayer region, using only data on the variation of condensation pressures with pore diameter. The results indicate a minimum mesopore diameter for the surface layering model to hold as 14.1 Å, below which size micropore filling must occur, and a minimum pore diameter for mechanical stability of the hemispherical meniscus during desorption as 34.2 Å. For pores in-between these two sizes reversible condensation is predicted to occur, in accord with the experimental data for nitrogen adsorption on MCM-41 at 77.4 K.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Reactivity of chars and carbons: New insights through molecular modeling

A new model proposed for the gasification of chars and carbons incorporates features of the turbostratic nanoscale structure that exists in such materials. The model also considers the effect of initial surface chemistry and different reactivities perpendicular to the edges and to the faces of the underlying crystallite planes comprising the turbostratic structure. It may be more realistic than earlier models based on pore or grain structure idealizations when the carbon contains large amounts of crystallite matter. Shrinkage of the carbon particles in the chemically controlled regime is also possible due to the random complete gasification of crystallitic planes. This mechanism can explain observations in the literature of particle size reduction. Based on the model predictions, both initial surface chemistry and the number of stacked planes in the crystallites strongly influence the reactivity and particle shrinkage. Its test results agree well with literature data on the air-oxidation of Spherocarb and show that it accurately predicts the variation of particle size with conversion. Model parameters are determined entirely from rate measurements.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Hemichordates and deuterostome evolution: robust molecular phylogenetic support for a hemichordate plus echinoderm clade

Hemichordates were traditionally allied to the chordates, but recent molecular analyses have suggested that hemichordates are a sister group to the echinoderms, a relationship that has important consequences for the interpretation of the evolution of deuterostome body plans. However, the molecular phylogenetic analyses to date have not provided robust support for the hemichordate + echinoderm clade. We use a maximum likelihood framework, including the parametric bootstrap, to reanalyze DNA data from complete mitochondrial genomes and nuclear 18S rRNA. This approach provides the first statistically significant support for the hemichordate + echinoderm clade from molecular data. This grouping implies that the ancestral deuterostome had features that included an adult with a pharynx and a dorsal nerve cord and an indirectly developing dipleurula-like larva.

Disagreement between correlations of quantum mechanics and stochastic electrodynamics in the damped parametric oscillator

Intracavity and external third order correlations in the damped nondegenerate parametric oscillator are calculated for quantum mechanics and stochastic electrodynamics (SED), a semiclassical theory. The two theories yield greatly different results, with the correlations of quantum mechanics being cubic in the system's nonlinear coupling constant and those of SED being linear in the same constant. In particular, differences between the two theories are present in at least a mesoscopic regime. They also exist when realistic damping is included. Such differences illustrate distinctions between quantum mechanics and a hidden variable theory for continuous variables.