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.     

Fractional Derivatives Applied to Phase-Space Reconstructions

The concept and application of phase-space reconstructions are reviewed. Fractional derivatives are then proposed for the purpose of reconstructing dynamics from a single observed time history. A procedure is presented in which the fractional derivatives of time series data are obtained in the frequency domain. The method is applied to the Lorenz system. The ability of the method to unfold the data is assessed by the method of global false nearest neighbors. The reconstructed data is used to compute recurrences and correlation dimensions. The reconstruction is compared to the commonly used method of delays in order to assess the choice of reconstruction parameters, and also the quality of results.     

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.     

Robust handling of non‐linear constraints for GA optimization of aerodynamic shapes

A new approach to the robust handling of non‐linear constraints for GAs (genetic algorithms) optimization is proposed. A specific feature of the approach consists of the change in the conventional search strategy by employing search paths which pass through both feasible and infeasible points (contrary to the traditional approach where only feasible points may be included in a path). The method (driven by full Navier–Stokes computations) was applied to the problem of multiobjective optimization of aerodynamic shapes subject to various geometrical and aerodynamic constraints. The results demonstrated that the method retains high robustness of conventional GAs while keeping CFD computational volume to an acceptable level, which allowed the algorithm to be used in a demanding engineering environment. 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.     

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     

Dark cure studies of cationic photopolymerizations of epoxides: Characterization of kinetic rate constants at high conversions

The effective propagation rate constant (k_p; averaged over all the propagating active centers) was characterized for solvent‐free cationic photopolymerizations of phenyl glycidyl ether over the entire range of conversions, including the high conversion regime in which mass transfer limitations become important. The profile for the k_p as a function of conversion was found to exhibit a constant plateau value at low to intermediate conversions, followed by a monotonic increase above a threshold value of conversion. To explain this trend, it is proposed that at high conversion the diffusional mobility of the photoinitiator counterion is reduced whereas the mobility of the cationic active center remains high because of reactive diffusion. Therefore, with increasing conversion, the average distance between the active centers and counterions may increase, resulting in an increase in the propagation rate constant. The profiles for the k_p values were investigated as a function of the temperature, photoinitiator anion, and photoinitiator concentration. As the photoinitiator concentration was increased, the plateau value of the effective propagation rate constant decreased whereas the threshold conversion increased. All of the experimental trends are consistent with the proposed increase in ion separation at high conversions. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4409–4416, 2004     

Photopolymerization of 1,6‐hexanedioldiacrylate initiated by three‐component systems based on N‐arylphthalimides

Three‐component photoinitiators comprised of an N‐arylphthalimide, a diarylketone, and a tertiary amine were investigated for their initiation efficiency of acrylate polymerization. The use of an electron‐deficient N‐arylphthalimide resulted in a greater acrylate polymerization rate than an electron‐rich N‐arylphthalimide. Triplet energies of each N‐arylphthalimide, determined from their phosphorescence spectra, and the respective rate constants for triplet quenching by the N‐arylphthalimide derivatives (acquired via laser flash photolysis) indicated that an electron–proton transfer from an intermediate radical species to the N‐arylphthalimide (not energy transfer from triplet sensitization) is responsible for generating the initiating radicals under the conditions and species concentrations used for polymerization. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4009–4015, 2004     

Preparation of macromonomers by copolymerization of methyl acrylate dimer involving β fragmentation

The unsaturated dimer of methyl acrylate [CH₂C(CO₂CH₃)CH₂CH₂CO₂CH₃, or MAD] was copolymerized with various monomers to prepare copolymers bearing the ω-unsaturated end group [CH₂C(CO₂CH₃)CH₂ ] arising from β fragmentation of the MAD propagating radical. Copolymerizations of MAD with cyclohexyl and n-butyl acrylate resulted in copolymers with ω-unsaturated end groups, and increasing the temperature up to 180 °C resulted in an increase in the rate of β fragmentation of MAD radicals relative to propagation. Only a small amount of unsaturated end groups was introduced by copolymerization with ethyl methacrylate (EMA), and the EMA content in the copolymer increased with temperature. These findings could be explained by the reversible addition of the poly(EMA) radical to MAD. The copolymerization with ethyl α-ethyl acrylate (EEA) did yield a copolymer containing unsaturated end groups with MAD units as part of the main chain, although the steric hindrance of the ethyl group suppressed homopropagation and crosspropagation of EEA, resulting in low polymerization rates. Therefore, the copolymerization of MAD with acrylic esters at high temperatures was noted as a convenient route for obtaining acrylate–MAD copolymers bearing unsaturated end groups at the ω end (macromonomer). © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 597–607, 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