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.     

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.     

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     

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     

Completely biodegradable composites of poly(propylene carbonate) and short,lignocellulose fiber Hildegardia populifolia

The composites of biodegradable poly(propylene carbonate) (PPC) reinforced with short Hildegardia populifolia natural fiber were prepared by melt mixing followed by compression molding. The mechanical properties, thermal properties, and morphologies of the composites were studied via static and dynamic mechanical measurements, thermogravimetric analysis, and scanning electron microscopy (SEM) techniques, respectively. Static tensile tests showed that the stiffness and tensile strength of the composites increased with an increasing fiber content. However, the elongation at break and the energy to break decreased dramatically with the addition of short fiber. The relationship between the experimental results and the compatibility or interaction between the PPC matrix and fiber was correlated. SEM observations indicated good interfacial contact between the short fiber and PPC matrix. Thermogravimetric analysis revealed that the introduction of short Hildegardia populifolia fiber led to a slightly improved thermooxidative stability of PPC. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 666–675, 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