A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

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.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Photochemical attachment of polymers on planar surfaces with a covalently anchored monolayer of a novel naphthyl ketone photochemical radical generator

A monolayer of covalently anchored, novel, binaphthyl ketone is used as a surface‐confined photochemical radical generator (PRG) for anchoring a variety of polymers to silicon surfaces. The precursor PRG is synthesized by the application of a facile and novel method for the oxidation of sterically hindered benzylic hydrocarbons to carbonyl compounds. Oxidation was carried out with a stoichiometric amount of potassium peroxydisulfate, in the presence of a catalytic amount of copper sulfate in an acetonitrile/water mixture. The PRG synthesized is characterized by ¹H NMR, UV, and Fourier transform infrared (FTIR). The covalently attached monolayers are characterized by X‐ray photoelectron spectroscopy, ellipsometry, and water contact angle measurements. The method developed is applicable to the preparation of a monolayer of a variety of polymers on a wide range of substrates carrying surface hydroxyl groups. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5413–5423, 2004     

Tetrafunctional initiators for cationic polymerization of olefins

3,3′,5,5′‐Tetrakis(2‐chloro‐2‐propyl)biphenyl (biphenyl tetracumyl chloride, BPTCC) and 1,3‐bis[3,5‐bis(2‐chloro‐2‐propyl)phenoxy]propane (diphenoxypropane tetracumyl chloride, DPPTCC) were synthesized as initiators for quasiliving cationic polymerization of isobutylene (IB). In the synthesis of BPTCC, tetrafunctionality was achieved via the coupling of dimethyl 5‐bromoisophthalate (DMBI) using nickel dibromide bis(triphenylphosphine) and zinc in the presence of a base; in the synthesis of DPPTCC, two equivalents of dimethyl 5‐hydroxyisophthalate were linked via reaction with 1,3‐dibromopropane in the presence of potassium carbonate. Both initiators were used to initiate the polymerization of IB under quasiliving cationic polymerization conditions. PIB initiated from BPTCC revealed a chain end/molecule value (as determined by ¹H‐NMR) of 3.85, verifying the nearly exclusive production of 4‐arm polyisobutylene (PIB). GPC analysis revealed a narrow peak representing the target four‐arm PIB, with a slight shoulder at high elution volumes (low molecular weights). GPC analysis of the PIB initiated by DPPTCC revealed multimodal distributions, suggesting the formation of two‐, three‐, and four‐arm star polymers during the polymerization. This behavior was attributed to Friedel–Crafts alkylation of the initiator core after the addition of one IB unit, which was activated by the electron‐donating oxytrimethyleneoxy linking moiety. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5942–5953, 2004     

New liquid‐crystalline poly(ester amide)s: The role of nitro groups in the phase behavior

New hydrogen‐bonded liquid‐crystalline poly(ester amide)s (PEA)s were obtained from 1,4‐terephthaloyl[bis‐(3‐nitro‐N‐anthranilic acid)] (5) or 1,4‐terephthaloyl[bis‐(N‐anthranilic acid)] (6), with or without nitro groups, respectively, through the separate condensation of each with hydroquinone or dihydroxynaphthalene. The dicarboxylic monomers were synthesized from 2‐aminobenzoic acid. The phase behavior of the monomers and polymers were studied with differential scanning calorimetry, polarized light microscopy, and wide‐angle X‐ray diffraction methods. Monomer 5, containing nitro groups, exhibited a smectic liquid‐crystalline phase, whereas the texture of monomer 6 without nitro groups appeared to be nematic. The PEAs containing nitro groups exhibited polymorphism (smectic and nematic), whereas those without nitro groups exhibited only one phase transition (a nematic threaded texture). The changes occurring in the phase behavior of the polymers were explained by the introduction of nitro groups. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1289–1298, 2004     

The effect of compatibilization and rheological properties of polystyrene and poly(dimethylsiloxane) on phase structure of polystyrene/poly(dimethylsiloxane) blends

The compatibilization effect of polystyrene (PS)‐poly(dimethylsiloxane) (PDMS) diblock copolymer (PS‐b‐PDMS) and the effect of rheological properties of PS and PDMS on phase structure of PS/PDMS blends were investigated using a selective extraction technique and scanning electron microscopy (SEM). The dual‐phase continuity of PS/PDMS blends takes place in a wide composition range. The formation and the onset of a cocontinuous phase structure largely depend on blend composition, viscosity ratio of the constituent components, and addition of diblock copolymers. The width of the concentration region of the cocontinuous structure is narrowed with increasing the viscosity ratio of the blends and in the presence of the small amount diblock copolymers. Quiescent annealing shifts the onset values of continuity. The experimental results are compared with the volume fraction of phase inversion calculated with various theoretical models, but none of the models can account quantitatively for the observed data. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 898–913, 2004