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.     

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     

Hollow latex particles as submicrometer reactors for polymerization in confined geometries

A method was developed for free‐radical polymerization in the confines of a hollow latex particle. Hollow particles were prepared via the dynamic swelling method from polystyrene seed and divinylbenzene and had hollows of 500–1000 nm. So that these hollow poly(divinylbenzene) particles could function as submicrometer reactors, the particles were filled with a monomer (N‐isopropylacrylamide) via the dispersion of the dried particles in the molten monomer. The monomer that was not contained in the hollows was removed by washing and gentle abrasion. Free‐radical polymerization was then initiated by γ radiolysis in the solid state. Transmission electron microscopy showed that poly(N‐isopropylacrylamide) formed in the hollow interior of the particles, which functioned as submicrometer reactors. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5706–5713, 2004     

Plasticized microporous poly(vinylidene fluoride) separators for lithium‐ion batteries. II. Poly(vinylidene fluoride) dense membrane swelling behavior in a liquid electrolyte—characterization of the swelling kinetics

Some microporous poly(vinylidene fluoride) (PVdF) separators for lithium‐ion batteries, used in liquid organic electrolytes based on a mixture of carbonate solvents and lithium salt LiPF₆, were characterized by the study of the swelling phenomena on dense PVdF membranes. Various aspects of the kinetics of the carbonate solvents and the solvent mixture sorption in dense PVdF slabs were studied at different temperatures. Non‐Fickian behavior, characterized by S‐shaped sorption curves, was highlighted, and a salt effect, which resulted in two‐stage sorption, was studied. Diffusion coefficients and activation energies were calculated for the Fickian portions of the sorption curves, that is, at short times and low swelling ratios. A strong influence of the different interaction parameters was shown for the swelling kinetics. This study proved that the swelling of microporous PVdF membranes could be considered instantaneous. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 544–552, 2004     

Rheooptical investigation of the crystallization of poly(ethylene terephthalate) under tensile strain

The crystallization of poly(ethylene terephthalate) under uniaxial tensile strain at different extension rates was investigated with optical polarimetry in a temperature range between the glass-transition temperature and the quiescent crystallization temperature. The evolution of the optical properties of the polymer, including the turbidity, birefringence, and dichroism, were monitored simultaneously with the mechanical parameters. To complete the semicrystalline microstructure characterization of the polymer under strain, an online wide-angle X-ray diffraction (WAXD) technique was used in separate experiments, which were performed under the same thermomechanical conditions. For real-time measurements, a high-energy synchrotron radiation source was used. The optical properties provided information about both the crystalline and amorphous phases, whereas the WAXD patterns essentially gave information about the crystalline phase. The two experimental techniques were then used in a complementary way to characterize the semicrystalline microstructure. Significant deviations from the stress-optical rule were found. This was attributed to both transient effects and the appearance of crystallites, which consisted of highly oriented molecular segments that could contribute to the optical anisotropy but not necessarily to the stress. The behavior of the optical dichroism was found to be qualitatively different from that of the birefringence. The latter monotonically increased with the strain, whereas the former first increased with the strain, passed through a maximum, and then decreased to a steady-state value. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1915–1927, 2004     

Plasticized microporous poly(vinylidene fluoride) separators for lithium-ion batteries. III. Gel properties and irreversible modifications of poly(vinylidene fluoride) membranes under swelling in liquid electrolytes

Microporous poly(vinylidene fluoride) (PVdF) separators for lithium-ion batteries, used in liquid organic electrolytes, have been characterized with respect to the swelling phenomena on dense PVdF membranes (obtained through hot pressing). In the first and second parts of this study, we have described the swelling equilibria and swelling kinetics of dense PVdF. Here the thermal properties of PVdF gels and their irreversible modifications induced by swelling are characterized. Particular attention is paid to crystallinity modifications, polymer plasticization, and membrane degradation. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2308–2317, 2004