Monte Carlo simulation for the modification of polymer via grafting

Monte Carlo method has been applied to investigate the kinetic of grafting reaction in free radical copolymerization. The simulation is quits in agreement with that of theoretical and experimental results. It proves that the Monte Carlo simulation is an effective method for investigating the grafting reaction of free radical copolymerization. (C) 2000 Elsevier Science Ltd. All rights reserved.

Rarefied gas dynamics: fundamentals, simulations and micro flows

Development of an efficient particle approach for micro-scale gas flow simulations

The micro-scale gas flows are usually low-speed flows and exhibit rarefied gas effects. It is challenging to simulate these flows because traditional CFD method is unable to capture the rarefied gas effects and the direct simulation Monte Carlo (DSMC) method is very inefficient for low-speed flows. In this study we combine two techniques to improve the efficiency of the DSMC method. The information preservation technique is used to reduce the statistical noise and the cell-size relaxed technique is employed to increase the effective cell size. The new cell-size relaxed IP method is found capable of simulating micro-scale gas flows as shown by the 2D lid-driven cavity flows.

Stress transfer and damage evolution simulations of fiber-reinforced polymer-matrix composites

The stress transfer from broken fibers to unbroken fibers in fiber-reinforced thermosetting polymer-matrix composites and thermoplastic polymer-matrix composites was studied using a detailed finite element model. In order to check the validity of this approach, an epoxy-matrix monolayer composite was used as thermosetting polymer-matrix composite and a polypropylene (PP)-matrix monolayer composite was used as thermoplastic polymer-matrix composite, respectively. It is found that the stress concentrations near the broken fiber element cause damage to the neighboring epoxy matrix prior to the breakage of other fibers, whereas in the case of PP-matrix composites the fibers nearest to the broken fiber break prior to the PP matrix damage, because the PP matrix around the broken fiber element yields. In order to simulate composite damage evolution, a Monte Carlo technique based on a finite element method has been developed in the paper. The finite element code coupled with statistical model of fiber strength specifically written for this problem was used to determine the stress redistribution. Five hundred samples of numerical simulation were carried out to obtain statistical deformation and failure process of composites with fixed fiber volume fraction.

Exact solutions to collisionless external gas flows

This paper presents exact density, velocity and temperature solutions for two problems of collisionless gas flows around a flat plate or a spherical object. At any point off the object, the local velocity distribution function consists of two pieces of Maxwellian distributions: one for the free stream which is characterized by free stream density, temperature and average velocity, n0, T0, U0; and the other is for the wall and it is characterized by density at wall and wall temperature, nw,Tw. Directly integrating the distribution functions leads to complex but exact flowfield solutions. To validate these solutions, we perform numerical simulations with the direct simulation Monte Carlo (DSMC) method. In general, the analytical and numerical results are virtually identical. The evaluation of these analytical solutions only requires less than one minute while the DSMC simulations require several days.

Size distributions of quantum islands on stepped substrates

The size distributions of self-assembled quantum islands on stepped substrates are studied using kinetic Monte Carlo simulations. It is found that the energy barrier E-SW between the step and the terrace region is the key factor in affecting the size distribution of islands. With small E-SW (<= 0.1 eV), lines of uniform islands can be obtained at relative low surface coverage. As the surface coverage is increased, wirelike islands can be obtained. Scaling behavior is obeyed for the size distributions of the wirelike islands. When the size distributions are separated into their width and length components, however, scaling is only observed in the length distribution of the wirelike islands. With larger E-SW, the size distribution of islands shows a clear bimodal size distribution and anomalous growth temperature dependent island size evolutions are observed. The simulation results reproduce qualitatively the phenomena observed in the cases of InAs islands grown on stepped GaAs substrates. (c) 2009 American Institute of Physics. [doi:10.1063/1.3248367]

Formation trends of ordered self-assembled nanoislands on stepped substrates

The growth of ordered self-assembled nanoislands on stepped substrates is studied systematically by kinetic Monte Carlo simulations. As the terrace width is small, the formation of nanoislands is confined in the steps and nanoislands ordered in lines or nanowires can be obtained. The Schwoebel barrier at the step edges has a great influence on the evolution of both the size and space distributions of the islands. When the terrace width is relatively large, self-ordering of nanoislands in the center regions of the terraces happens. An unexpected trend of the nanoisland self-ordering is found as the deposition thickness is larger than 0.2 ML, which can be related to the attractive migrations between nearby islands.

Experimental study of two-proton correlated emission from S-29 excited states

An experiment of a S-29 beam bombarding a Au-197 target at an energy of 49.2 MeV/u has been performed to study the two-proton correlated emission from S-29 excited states. Complete-kinematics measurements were carried out in the experiment. The relative momentum, opening angle, and relative energy of two protons, as well as the invariant mass of the final system, were deduced by relativistic-kinematics reconstruction. The Si-27-p-p coincident events were picked out under strict conditions and the phenomenon of p-p correlations was observed among these events. The mechanisms of two-proton emission were analyzed in a simple schematic model, in which the extreme decay modes like He-2 cluster emission, three-body phase-space decay, and two-body sequential emission were taken into account. Associated with the Monte Carlo simulations, the present results show that two protons emitted from the excited states between 9.6 MeV and 10.4 MeV exhibit the features of He-2 cluster decay with a branching ratio of 29(-11)(+10)%.

Investigation of two-proton emission from excited states of the odd-Z nucleus P-28 by complete-kinematics measurements

An experiment to study exotic two-proton emission from excited levels of the odd-Z nucleus P-28 was performed at the National Laboratory of Heavy Ion Research-Radioactive Ion Beam Line (HIRFL-RIBLL) facility. The projectile P-28 at the energy of 46.5 MeV/u was bombarding a Au-197 target to populate the excited states via Coulomb excitation. Complete-kinematics measurements were realized by the array of silicon strip detectors and the CsI + PIN telescope. Two-proton events were selected and the relativistic-kinematics reconstruction was carried out. The spectrum of relative momentum and opening angle between two protons was deduced from Monte Carlo simulations. Experimental results show that two-proton emission from P-28 excited states less than 17.0 MeV is mainly two-body sequential emission or three-body simultaneous decay in phase space. The present simulations cannot distinguish these two decay modes. No obvious diproton emission was found.

Design studies on the 4 pi gamma-ray calorimeter for the ETF experiment at HIRFL-CSR

A high detection efficiency calorimeter which is used to detect gamma-rays with energies from 1 MeV up to 10 MeV as well as light charged particles has been proposed. Design of the geometry, results of the crystal tests and Monte Carlo simulations are presented in this paper. The simulation results confirm that the calorimeter can obtain high detection efficiency and good energy resolution with the current designed geometry. And the calorimeter is competent for the future External Target Facility (ETIF) experiments.

Effect of the nature of annealing solvent on the morphology of diblock copolymer blend thin films

The morphologies and structures for the thin film of blend systems consisting of two asymmetric polystyrene-block-polybutadiene (SB) diblock copolymers induced by annealing in the vapor of different solvents, namely, cyclohexane, benzene, and heptane, which have different selectivity or preferential affinity for a certain block, were investigated by tapping mode atomic force microscopy (AFM) and transmission electron microscopy (TEM). The results revealed that even a slight preferential affinity of good solvent for one block would strongly alter the morphology of the blend thin film.

Microphase Separation of Mixed Polymer Brushes: Dependence of the Morphology on Grafting Density, Composition, Chain-Length Asymmetry, Solvent Quality, and Selectivity

Microphase separation of binary mixed A/B polymer brushes exposed to different solvents is studied using Single-Chain-in-Mean-Field simulations. Effects of solvent quality and selectivity, grafting density, composition, and chain-length asymmetry are systematically investigated, and diagrams of morphologies in various solvents are constructed as a function of grafting density and composition or chain-length asymmetry. The structure of the microphase segregated morphologies lacks long-range periodic order, and it is analyzed quantitatively Using Minkowski measures.

Effects of polar group saturation on physical gelation of Amphiphilic polymer solutions

Monte Carlo simulation on the basis of the comblike coarse grained nonpolar/polar (NP) model has been carried out to study the polar group saturation effect on physical gelation of amphiphilic polymer solutions. The effects of polar group saturation due to hydrogen bonding or ion bridging on the sol-gel phase diagram, microstructure of aggregates, and chain conformation of amphiphilic polymer solutions under four different solvent conditions to either the nonpolar backbone or the polar side chain in amphiphilic polymer chains have been investigated. It is found that an increase of polar group saturation results in a monotonically decreased critical concentration of gelation point, which can be qualitatively supported by the dynamic theological measurements on pectin aqueous solutions. Furthermore, various solvent conditions to either the backbone or the side chain have significant impact on both chain conformation and microstructure of aggregates. When the solvent is repulsive to the nonpolar backbone but attractive to the polar side chain, the polymer chains are collapsed, and the gelation follows the mechanism of colloidal packing; at the other solvent conditions, the gelation follows the mechanism of random aggregation.

Nano-reactors for controlling the selectivity of the free radical grafting of maleic anhydride onto polypropylene in the melt

This paper reports on a successful application of the concept of nanoreactors to effectively controlling the selectivity of the free radical grafting of maleic anhydride (MAH) onto polypropylene (PP) in the melt, an industrially relevant process. More specifically, a free radical initiator of type ROOR was first confined into (or encapsulated by) the galleries of an organically modified montmorillonite (o-MMT) whose interdistance was 2.4 nm. Primary free radicals (RO center dot) formed inside the o-MMT galleries had to diffuse out before they could react with the PP backbone. The controlled release of the primary free radicals significantly increased the grafting degree of MAH onto PP and greatly reduced the level of the chain scission of the latter. Those results were better understood by electron spin resonance studies on model systems and by Monte Carlo simulations.