Detecting instabilities in flows of viscoelastic fluids

There is a growing interest in developing numerical tools to investigate the onset of physical instabilities observed in experiments involving viscoelastic flows, which is a difficult and challenging task as the simulations are very sensitive to numerical instabilities. Following a recent linear stability analysis carried out in order to better understand qualitatively the origin of numerical instabilities occurring in the simulation of flows viscoelastic fluids, the present paper considers a possible extension for more complex flows. This promising method could be applied to track instabilities in complex (i.e. essentially non‐parallel) flows. In addition, results related to transient growth mechanism indicate that it might be responsible for the development of numerical instabilities in the simulation of viscoelastic fluids. Copyright © 2003 John Wiley & Sons, Ltd.     

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 numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. 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.     

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