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.     

A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

An improved numerical simulator for multiphase flow in porous media

In basin modelling the thermodynamics of a multicomponent multiphase fluid flux are computationally too expensive when derived from an equation of state and the Gibbs equality constraints. In this article we present a novel implicit molar mass formulation technique using binary mixture thermodynamics. The two proposed solution methods, with and without cross derivative terms between components, are based on a preconditioned Newton‐GMRES scheme for each time‐step with analytical computation of the derivatives. These new algorithms reduce significantly the numerical effort for the computation of the molar masses, and we illustrate the behavior of these methods with numerical computations. Copyright © 2004 John Wiley & Sons Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

Particle correlations in e + e − → W + W − events

Preliminary results are reported on the two-particle correlation function R(Q) in hadronic Z decays, fully hadronic WW decays, and mixed hadronic-leptonic WW decays using data collected by the DELPHI detector at LEP at energies between 189 and 206 GeV. Evidence for Bose-Einstein correlations was observed in all three cases. The event mixing technique was used to determine correlations between particles arisingfrom different W bosons in fully hadronic WW decays. An excess of like-sign particle pairs with low four-momentum difference in fully hadronic WW events is observed, consistent with the effect expected from correlations between identical particles from different W bosons.     

Testing underwater bottom-moored antenna arrays in the sea and in a man-made lake

Measurements of response, gain, and noise immunity are carried out for an underwater compensated additive receiving array with randomly spaced hydrophones that is moored at the bottom of a man-made lake with multimode sound propagation. The in-sea locating ability of a similar array is demonstrated with the sources of noiselike signals at frequencies of 5–100 Hz. A dedicated numerical processor is developed and tested for processing the signals received by a random underwater array.     

Photoinitiators with functional groups. VII. Covalently bonded camphorquinone—amines

Camphorquinone (CQ), a widely used photoinitiator (PI) in dental applications, was covalently bonded to aromatic amines to enhance the rate of electron and proton transfer effect due to the close vicinity of the diketone and the amine group. 10‐bromocamphorquinone and 10‐bromomethylcamphorquinone were selected as suitable precursors for esterification with the carboxyl group containing aromatic amines based on 4‐dimethylaminobenzoic acid. Properties of the new photoinitiating systems were investigated by UV spectroscopy and differential scanning photocalorimetry in lauryl acrylate. Compared to physical mixtures, in all cases similar or even better performance was obtained. Surprisingly, 10‐acetyl derivatives 7 – 9 and 18 especially, were found to be highly reactive. Compared to CQ/ethyl 4‐dimethylaminobenzoate, the rate of photopolymerization was increased by a factor of up to 2. Intramolecular reaction was confirmed by photo‐differential scanning calorimetry experiments with varying PI concentrations. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4948–4963, 2004     

Atom transfer radical polymerization of methyl methacrylate by copper(I)‐pyridine‐2‐carboximidate catalysts

Pyridine‐2‐carboximidates [methyl ( 1a ), ethyl ( 1b ), isopropyl ( 1c ), cyclopentyl ( 1d ), cyclohexyl ( 1e ), n‐octyl ( 1f ), and benzyl ( 1g )] were prepared from the reaction of 2‐cyanopyridine with the corresponding alcohols. Cyclopentyl‐substituted 1d was found to be a highly effective ligand for copper‐catalyzed atom transfer radical polymerization (ATRP) of methyl methacrylate (MMA). For example, the observed rate constant for a CuBr/ 1d catalytic system was found to be nearly twice as high as the cyclohexyl‐substituted CuBr/ 1e catalytic system [k_obs = (1.19 vs 0.56) × 10^?4 s^?1). The effects of the solvents, temperature, catalyst/initiator, and solvent/monomer ratio on the ATRP of MMA were studied systematically for the CuBr/ 1d catalytic system. The optimum condition for the ATRP of MMA was found to be a 1:2:1:400 [CuBr]_o/[ 1d ]_o/[ethyl 2‐bromoisobutyrate]_o/[MMA]_o ratio at 60 °C in veratrole solution, which yielded well‐defined poly(MMA) with a narrow molecular weight distribution of 1.14. The catalytically active copper complex 2d was isolated from the reaction of CuBr with 1d . Narrow molecular weight distributions as low as 1.06 were achieved for the CuBr/ 1d catalytic system by employing 10% of the deactivator CuBr₂. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 2747–2755, 2004     

Viscoelastic behavior of polypropylene/nitrile rubber thermoplastic elastomer blends: Application of Kerner's models for reactively compatibilized and dynamically vulcanized systems

The viscoelastic properties of binary blends of nitrile rubber (NBR) and isotactic polypropylene (PP) of different compositions have been calculated with mean‐field theories developed by Kerner. The phase morphology and geometry have been assumed, and experimental data for the component polymers over a wide temperature range have been used. Hashin's elastic–viscoelastic analogy principle is used in applying Kerner's theory of elastic systems for viscoelastic materials, namely, polymer blends. The two theoretical models used are the discrete particle model (which assumes one component as dispersed inclusions in the matrix of the other) and the polyaggregate model (in which no matrix phase but a cocontinuous structure of the two is postulated). A solution method for the coupled equations of the polyaggregate model, considering Poisson's ratio as a complex parameter, is deduced. The viscoelastic properties are determined in terms of the small‐strain dynamic storage modulus and loss tangent with a Rheovibron DDV viscoelastometer for the blends and the component polymers. Theoretical calculations are compared with the experimental small‐strain dynamic mechanical properties of the blends and their morphological characterizations. Predictions are also compared with the experimental mechanical properties of compatibilized and dynamically cured 70/30 PP/NBR blends. The results computed with the discrete particle model with PP as the matrix compare well with the experimental results for 30/70, 70/30, and 50/50 PP/NBR blends. For 70/30 and 50/50 blends, these predictions are supported by scanning electron microscopy (SEM) investigations. However, for 30/70 blends, the predictions are not in agreement with SEM results, which reveal a cocontinuous blend of the two. Predictions of the discrete particle model are poor with NBR as the matrix for all three volume fractions. A closer agreement of the predicted results for a 70/30 PP/NBR blend and the properties of a 1% maleic anhydride modified PP or 3% phenolic‐modified PP compatibilized 70/30 PP/NBR blend in the lower temperature zone has been observed. This may be explained by improved interfacial adhesion and stable phase morphology. A mixed‐cure dynamically vulcanized system gave a better agreement with the predictions with PP as the matrix than the peroxide, sulfur, and unvulcanized systems. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1417–1432, 2004     

Behavior and surface energies of polybenzoxazines formed by polymerization with argon,oxygen, and hydrogen plasmas

Polybenzoxazine (PBZZ) thin films can be fabricated by the plasma‐polymerization technique with, as the energy source, plasmas of argon, oxygen, or hydrogen atoms and ions. When benzoxazine (BZZ) films are polymerized through the use of high‐energy argon atoms, electronegative oxygen atoms, or excited hydrogen atoms, the PBZZ films that form possess different properties and morphologies in their surfaces. High‐energy argon atoms provide a thermodynamic factor to initiate the ring‐opening polymerization of BZZ and result in the polymer surface having a grid‐like structure. The ring‐opening polymerization of the BZZ film that is initiated by cationic species such as oxygen atoms in plasma, is propagated around nodule structures to form the PBZZ. The excited hydrogen atom plasma initiates both polymerization and decomposition reactions simultaneously in the BZZ film and results in the formation of a porous structure on the PBZZ surface. We evaluated the surface energies of the PBZZ films polymerized by the action of these three plasmas by measuring the contact angles of diiodomethane and water droplets. The surface roughness of the films range from 0.5 to 26 nm, depending on the type of carrier gas and the plasma‐polymerization time. By estimating changes in thickness, we found that the PBZZ film synthesized by the oxygen plasma‐polymerization process undergoes the slowest rate of etching in CF₄ plasma. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 4063–4074, 2004