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.     

Stability and slightly supercritical oscillatory regimes of natural convection in a 8:1 cavity: solution of the benchmark problem by a global Galerkin method

The global Galerkin method is applied to the benchmark problem that considers an oscillatory regime of convection of air in a tall two‐dimensional rectangular cavity. The three most unstable modes of the linearized system of the Boussinesq equations are studied. The converged values of the critical Rayleigh numbers together with the corresponding oscillation frequencies are calculated for each mode. The oscillatory flow regimes corresponding to each of the three modes are approximated asymptotically. No direct time integration is applied. Good agreement with the previously published results obtained by solution of the time‐dependent Boussinesq equations is reported. Copyright © 2004 John Wiley & Sons, Ltd.     

High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Standardization of BSE rapid test performances and experiences gathered during the implementation of large-scale testing

The accumulation of pathological prion protein is used as a diagnostic marker for transmissible spongiform encephalopathies. According to European Union (EU) regulations cattle older than 30 months of age (Germany, France, Italy, and Spain by national law >24 months) and slaughtered for human consumption must be tested by using rapid tests for bovine spongiform encephalopathy (BSE). Likewise fallen stock and clinically affected animals must be tested. This article gives a short overview of the incidence of BSE in Europe. The diagnostic hierarchy, i.e., the officially approved methodology for the confirmation of suspect rapid test cases, and the organization of the numerous laboratories involved in this large-scale testing for BSE are described. Special emphasis is given to necessary quality control measures currently in place for BSE rapid testing laboratories and to measures intended to assure a consistent performance of the commercially available rapid test kits.Presented at BERM-9—Ninth International Symposium on Biological and Environmental Reference Materials, 15–19 June 2003, Berlin, Germany.     

Polyelectrolyte complexes. VI. Polycation structure,polyanion molar mass,and polyion concentration effects on complex nanoparticles based on poly(sodium 2‐acrylamido‐2‐methylpropanesulfonate)

The formation and characterization of some interpolyelectrolyte complex (IPEC) nanoparticles based on poly(sodium 2‐acrylamido‐2‐methylpropanesulfonate) (NaPAMPS), as a function of the polycation structure, polyanion molar mass, and polyion concentration, were followed in this work. Poly(diallyldimethylammonium chloride) and two polycations (PCs) containing (N,N‐dimethyl‐2‐hydroxypropyleneammonium chloride) units in the backbone (PCA₅ and PCA₅D₁) were used as starting polyions. The complex stoichiometry, (n^?/n⁺)_iso, was pointed out by optical density at 500 nm (OD₅₀₀), polyelectrolyte titration, and dynamic light scattering. IPEC nanoparticle sizes were influenced by the polycation structure and polyanion molar mass only before the complex stoichiometry, which was higher for the more hydrophilic polycations (PCA₅ and PCA₅D₁) and for a higher NaPAMPS molar mass, and were almost independent of these factors after that, at a flow rate of the added polyion of about 0.28 mL × (mL PC)^?1 × h^?1. The IPEC nanoparticle sizes remained almost constant for more than 2 weeks, both before and after the complex stoichiometry, at low concentrations of polyions. NIPECs as stable colloidal dispersions with positive charges in excess were prepared at a ratio between charges (n^?/n⁺) of 0.7, and their storage colloidal stability, as a function of the polycation structure and polyion concentration (from 0.8 to ca. 7.8 mmol/L), was demonstrated. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 2495–2505, 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     

Relationship between the dynamic softening transition and wet sliding friction of elastomer compounds

For properly chosen elastomer compounds, thermorheological characterization is combined with an examination of the variation of the wet sliding friction with temperature. A conceptual argument leads to the assumption that the wet sliding friction should maximize at the energy dissipation peak associated with the dynamic softening transition at a characteristic frequency determined by the sliding speed and the effective smallest surface asperity scale. The dynamic softening transition is characterized with the peak in tan δ/G′ⁿ, where tan δ is the loss tangent, G′ is the elastic modulus, and n is a constant between 0 and 1. The William–Landel–Ferry transform is uncritically applied for extrapolating the position of the peak in tan δ/G′ⁿ at high frequencies. Even based on the criterion of tan δ, the results obtained on a concrete surface indicate that the effective smallest asperity scale is of order of 100 μm. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2467–2478, 2004     

Continuation of travelling‐wave solutions of the Navier–Stokes equations

An efficient way of obtaining travelling waves in a periodic fluid system is described and tested. We search for steady states in a reference frame travelling at the wave phase velocity using a first‐order pseudospectral semi‐implicit time scheme adapted to carry out the Newton's iterations. The method is compared to a standard Newton–Raphson solver and is shown to be highly efficient in performing this task, even when high‐resolution grids are used. This method is well suited to three‐dimensional calculations in cylindrical or spherical geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

Quality assurance for the analytical data of micro elements in food

The micro element content of food is an important quality index due to the action of these elements on human health. In this article, we discuss how to ensure the reliability of analytical data on micro elements in order to truly represent the condition of food. Sampling, treatment of the analytical sample, selection of the analytical method, standard solution, and certified reference material, blank test, calibration of the instrument and equipment, application of the quality control chart, assessment of the final analytical result, and quality assurance system are briefly described. Received: 5 July 2001 Accepted: 19 November 2001     

Order conditions for linearly implicit fractional step Runge Kutta methods

^{In this paper, we study the consistency of a variant of fractionalstep Runge–Kutta methods. These methods are designed tointegrate efficiently semi-linear multidimensional parabolicproblems by means of linearly implicit time integration processes.Such time discretization procedures are also related to a splittingof the space differential operator (or the spatial discretizationof it) as a sum of ‘simpler’ linear differentialoperators and a nonlinear term.}