Simulation Study of Aggregate Morphologies Formed by ABC Linear Triblock Copolymers in a Selective Solvent Through the Self-Consistent Field Theory

The effect of the hydrophobic properties of blocks B and C on the aggregate morphologies formed by ABC linear triblock copolymers in selective solvent was studied through the self-consistent field theory. Five typical micelles, such as core-shell-corona, hamburger-like, segmented-wormlike, were obtained by changing the hydrophobic properties of blocks B and C. The simulation results indicate that the shape and size of micelle are basically controlled by the hydrophobic degree of the middle block B, whereas the type of micelle is mainly determined by the hydrophobic degree of the end block C.

Monte Carlo simulation of the compatibility of graft copolymer compatibilized two incompatible homopolymer blends: Effect of graft structure

Compatibility of graft copolymer compatibilized two incompatible homopolymer A and B blends was simulated by using Monte Carlo method in a two-dimensional lattice model. The copolymers with various graft structures were introduced in order to study the effect of graft structure on the compatibility. Simulation results showed that incorporation of both A-g-B (A was backbone) and B-g-A (B was backbone) copolymers could much improve the compatibility of the blends. However, A-g-B copolymer was more effective to compatibilize the blend if homopolymer A formed dispersed phase. Furthermore, simulation results indicated that A-g-B copolymers tended to locate at the interface and anchor two immiscible components when the side chain is relatively long. However, most of A-g-B copolymers were likely to be dispersed into the dispersed homopolymer A phase domains if the side chains were relatively short. On the other hand, B-g-A copolymers tended to be dispersed into the matrix formed by homopolymer B. Moreover, it was found that more and more B-g-A copolymers were likely to form thin layers at the phase interface with decreasing the length of side chain.

Investigation of sol-gel transition in Pluronic F127/D2O solutions using a combination of small-angle neutron scattering and Monte Carlo simulation

Physical gelation in the concentrated Pluronic F127/D2O solution has been studied by a combination of small-angle neutron scattering (SANS) and Monte Carlo simulation. A 15% F127/D2O solution exhibits a sol-gel transition at low temperature and a gel-sol transition at the higher temperature, as evidenced by SANS and Monte Carlo simulation studies. Our SANS and simulation results also suggest that the sol-gel transition is dominated by the formation of a percolated polymer network, while the gel-sol transition is determined by the loss of bound solvent. Furthermore, different diffusion behaviors of different bound solvents and free solvent are observed. We expect that this approach can be further extended to study phase behaviors of other systems with similar sol-gel phase diagrams.

A Monte Carlo simulation for the micellization of ABC 3-miktoarm star terpolymers in a selective solvent

We have used Monte Carlo simulation to study the micellization of ABC 3-miktoarm star terpolymers in a selective solvent (good to A segment, bad to B and C segments). The simulation results reveal that the self-assembled morphology is determined by the block length, molecular architecture, terpolymer concentration and insolubility of insoluble block in the solvent. In dilute solution, symmetric terpolymers (N-B = N-C = 30) tend to aggregate into a novel wormlike pearl-necklace structure linked by an alternating arrangement of B and C spheres, whereas the asymmetric terpolymers (NB = 10, NC = 50) are likely to aggregate into spherical or cylindrical micelles (formed by C blocks) connected with some small B spheres, when the concentration of terpolymer is relatively low (chain number is 100). However, when the concentration of terpolymer is relatively high (chain number is 250), the symmetric terpolymers tend to aggregate into a netlike structure linked by an alternation of B and C spheres, whereas the asymmetric terpolymers are likely to aggregate into wormlike micelles (formed by C blocks) connected with some of small spheres (formed by B blocks). Moreover, when the insolubility of insoluble block in the solvent is weak, the insoluble blocks aggregate into some incompact micelles.

Monte Carlo simulation of the aggregation of rod-flexible triblock copolymers in a thin film

The aggregation of rod-flexible ABA and BAB triblock (A was rod block and repulsive with block B) copolymers in a thin film was studied as a function of varying the rigidity (eta) and the length of the rod block by Monte Carlo simulation. The rigidity of block A was defined as eta = R-c/R-max in this study. R-c, was the end-to-end distance below which the conformation of the block was not allowed, whereas R-max, was the longest end-to-end distance that the block could be. If eta = 0 the block was flexible, whereas if eta = 1 the block was a straight rod. The simulation results showed that the ABA triblock copolymer film were likely to form lamella structure with increasing the rigidity (eta) of block A. The lamellas were parallel each other and perpendicular to the film surface. However, the aggregation of BAB triblock copolymers tended to change from lamella to cylinder structure with increasing the rigidity (eta) of block A. Typical lamella and cylinder co-exist structure was obtained at eta = 0.504 for the BAB copolymer film. On the other hand, the simulation results indicated that the film changed from disorder to order, then to disorder structure with increasing the relative length of B block for both ABA and BAB copolymer films.

Monte Carlo simulation of miscibility of polymer blends with repulsive interactions: effect of chain structure

The effects of the chain structure and the intramolecular interaction energy of an A/B copolymer on the miscibility of the binary blends of the copolymer and homopolymer C have been studied by means of a Monte Carlo simulation. In the system, the interactions between segments A, B and C are more repulsive than those between themselves. In order to study the effect of the chain structure of the A/B copolymer on the miscibility, the alternating, random and block copolymers were introduced in the simulations, respectively. The simulation results show that the miscibility of the binary blends strongly depends on the intramolecular interaction energy ((ε) over bar (AB)) between segments A and B within the A/B copolymers. The higher the repulsive interaction energy, the more miscible the A/B copolymer and homopolymer C are. For the diblock copolymer/homopolymer blends, they tend to form micro phase domains. However, the phase domains become so small that the blend can be considered as a homogeneous phase for the alternating copolymer/ homopolymer blends. Furthermore, the investigation of the average end-to-end distance ((h) over bar) in different systems indicates that the copolymer chains tend to coil with the decrease Of (ε) over bar (AB) whereas the (h) over bar of the homopolymer chains depends on the chain structure of the copolymers.

Computer simulation of miscibility and self-assembly structure for polymer-containing systems with special interactions

The miscibility and structure of A-B copolymer/C homopolymer blends with special interactions were studied by a Monte Carlo simulation in two dimensions. The interaction between segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. In order to study the effect of copolymer chain structure on the morphology and structure of A-B copolymer/C homopolymer blends, the alternating, random and block A-B copolymers were introduced into the blends, respectively. The simulation results indicated that the miscibility of A-B block copolymer/C homopolymer blends depended on the chain structure of the A-B copolymer. Compared with alternating or random copolymer, the block copolymer, especially the diblock copolymer, could lead to a poor miscibility of A-B copolymer/C homopolymer blends. Moreover, for diblock A-B copolymer/C homopolymer blends, obvious self-organized core-shell structure was observed in the segment B composition region from 20% to 60%. However, if diblock copolymer composition in the blends is less than 40%, obvious self-organized core-shell structure could be formed in the B-segment component region from 10 to 90%. Furthermore, computer statistical analysis for the simulation results showed that the core sizes tended to increase continuously and their distribution became wider with decreasing B-segment component.

Monte Carlo simulation of the grafting of maleic anhydride onto polypropylene at higher temperature

The graft of maleic anhydride (MAH) onto isotactic polypropylene (iPP) initiated by dicumyl peroxide (DCP) at 190 degreesC was studied by means of the Monte Carlo method. The ceiling temperature theory, i.e., no possibility for the homopolymerization of MA-H to occur at higher temperatures, was used in this study. The simulation results show that most MAH monomers were grafted onto the radical chain ends arising from beta scission at a lower MAH concentration, whereas the amount of MAH monomers attached to the tertiary carbons was much larger than that grafted onto the radical chain ends at a higher MAH concentration for various DCP concentrations. This conclusion gives a good interpretation for the disagreement on the grafting sites along a PP chain. Moreover, it was found that the grafting degree increased considerably up to a peak value; thereafter, it decreased continuously with increasing MA-H concentration. The peak shifted in the lower MAH concentration direction and became lower and lower with increasing DCP concentration. When the DCP concentration was below 0.1 wt %, the peak was hardly observed. Those results are in good agreement with the experiments.

Monte Carlo simulation of phase behavior of polymer blends with special interactions

Full Paper: The phase, behavior of A-B-random copolymer/C-homopolymer, blends with special interaction was studied by a. Monte, Carlo simulation in two dimensions. The interaction between I segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. The simulation results showed that the blend became two large co-continuous phase domains at lower segment-B component compositions, indicating that the blend showed spinodal decomposition. With an increase of the segment-B component, the miscibility between the copolymer,and the polymer was gradually improved up to being miscible. In addition, it was found that segment B tended to move to the surface of the copolymer phase in the case of a lower component of segment B. On the other hand, if was observed that the average, end-to-end distances ((h) over bar) for both copolymer and polymer changed slowly with increasing segment-B component of the copolymer up to 40%, thereafter they increased considerably with increasing segment B component. Moreover, it was found that the (h) over bar of the copolymer was obviously shorter than that of the homopolymer for the segment-B composition, region from 0% to 80%. Finally, a, phase diagram showing I phase and - II phase regions under the condition of constant-temperature is presented.

Descriptive Simulation: Combining Symbolic and Numerical Methods in the Analysis of Chemical Reaction Mechanisms

The Kineticist's Workbench is a computer program currently under development whose purpose is to help chemists understand, analyze, and simplify complex chemical reaction mechanisms. This paper discusses one module of the program that numerically simulates mechanisms and constructs qualitative descriptions of the simulation results. These descriptions are given in terms that are meaningful to the working chemist (e.g., steady states, stable oscillations, and so on); and the descriptions (as well as the data structures used to construct them) are accessible as input to other programs.

Integration of micro- and macroscopic models for pedestrian evacuation simulation

Simulation of pedestrian evacuations of smart buildings in emergency is a powerful tool for building analysis, dynamic evacuation planning and real-time response to the evolving state of evacuations. Macroscopic pedestrian models are low-complexity models that are and well suited to algorithmic analysis and planning, but are quite abstract. Microscopic simulation models allow for a high level of simulation detail but can be computationally intensive. By combining micro- and macro- models we can use each to overcome the shortcomings of the other and enable new capability and applications for pedestrian evacuation simulation that would not be possible with either alone. We develop the EvacSim multi-agent pedestrian simulator and procedurally generate macroscopic flow graph models of building space, integrating micro- and macroscopic approaches to simulation of the same emergency space. By “coupling” flow graph parameters to microscopic simulation results, the graph model captures some of the higher detail and fidelity of the complex microscopic simulation model. The coupled flow graph is used for analysis and prediction of the movement of pedestrians in the microscopic simulation, and investigate the performance of dynamic evacuation planning in simulated emergencies using a variety of strategies for allocation of macroscopic evacuation routes to microscopic pedestrian agents. The predictive capability of the coupled flow graph is exploited for the decomposition of microscopic simulation space into multiple future states in a scalable manner. By simulating multiple future states of the emergency in short time frames, this enables sensing strategy based on simulation scenario pattern matching which we show to achieve fast scenario matching, enabling rich, real-time feedback in emergencies in buildings with meagre sensing capabilities.

The numerical simulation of the noncharring pyrolysis process and fire development within a compartment

A pyrolysis model for noncharring solid fuels is presented in this paper. Model predictions are compared with experimental data for the mass loss rates of polymethylmethacrylate (PMMA) and very good agreement is achieved. Using a three-dimensional CFD environment, the pyrolysis model is then coupled with a gas-phase combustion model and a thermal radiation model to simulate fire development within a small compartment. 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 post-flashover linear correlation with constant of proportionality 0.38 kg/sm5=2. The simulation results also suggest that the model is capable of predicting the onset of "flashover" and "post-flashover" type behaviour within the fire compartment.

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.

Coupled FDTD-FVM simulation of microwave heating of polymer materials for micro-systems packaging applications

Heating in an idealised polymer load in a novel open-ended variable frequency microwave oven is numerically simulated using a couple solver approach. The frequency-agile microwave oven bonding system (FAMOBS)is developed to meet rapid polymer curing requirements in microelectronics and optoelectronics manufacturing. The heating of and idealised polymer load has been investigated through numerical modelling. Assessment of the system comprises of simulation of electromagnetic fields and of temperature distribution within the load. Initial simulation results are presented and contrasted with experimental analysis of field distribution

Numerical simulation of encapsulant curing within a variable frequency microwave processing system

Curing of encapsulant material in a simplified microelectronics package using an open oven Variable Frequency Microwave (VFM) system is numerically simulated using a coupled solver approach. A numerical framework capable of simulating electromagnetic field distribution within the oven system, plus heat transfer, cure rate, degree of cure and thermally induced stresses within the encapsulant material is presented. The discrete physical processes have been integrated into a fully coupled solution, enabling usefully accurate results to be generated. Numerical results showing the heating and curing of the encapsulant material have been obtained and are presented in this contribution. The requirement to capture inter-process coupling and the variation in dielectric and thermophysical material properties is discussed and illustrated with simulation results.