Kinetic energy conserving integrators for Gaussian thermostatted SLLOD

A new integration scheme is developed for nonequilibrium molecular dynamics simulations where the temperature is constrained by a Gaussian thermostat. The utility of the scheme is demonstrated by its application to the SLLOD algorithm which is the standard nonequilibrium molecular dynamics algorithm for studying shear flow. Unlike conventional integrators, the new integrators are constructed using operator-splitting techniques to ensure stability and that little or no drift in the kinetic energy occurs. Moreover, they require minimum computer memory and are straightforward to program. Numerical experiments show that the efficiency and stability of the new integrators compare favorably with conventional integrators such as the Runge-Kutta and Gear predictor-corrector methods. (C) 1999 American Institute of Physics. [S0021-9606(99)50125-6].

Analysis of the free vibration of rectangular plates with central cut-outs using the discrete Ritz method

We present an analysis of the free vibration of plates with internal discontinuities due to central cut-outs. A numerical formulation for a basic L-shaped element which is divided into appropriate sub-domains that are dependent upon the location of the cut-out is used as the basic building element. Trial functions formed to satisfy certain boundary conditions are employed to define the transverse deflection of each sub-domain. Mathematical treatments in terms of the continuities in displacement, slope, moment, and higher derivatives between the adjacent sub-domains are enforced at the interconnecting edges. The energy functional results, from the proper assembly of the coupled strain and kinetic energy contributions of each sub-domain, are minimized via the Ritz procedure to extract the vibration frequencies and. mode shapes of the plates. The procedures are demonstrated by considering plates with central cut-outs that are subjected to two types of boundary conditions. (C) 2003 Elsevier Ltd. All rights reserved.

Prediction of protein B-factor profiles

The polypeptide backbones and side chains of proteins are constantly moving due to thermal motion and the kinetic energy of the atoms. The B-factors of protein crystal structures reflect the fluctuation of atoms about their average positions and provide important information about protein dynamics. Computational approaches to predict thermal motion are useful for analyzing the dynamic properties of proteins with unknown structures. In this article, we utilize a novel support vector regression (SVR) approach to predict the B-factor distribution (B-factor profile) of a protein from its sequence. We explore schemes for encoding sequences and various settings for the parameters used in SVR. Based on a large dataset of high-resolution proteins, our method predicts the B-factor distribution with a Pearson correlation coefficient (CC) of 0.53. In addition, our method predicts the B-factor profile with a CC of at least 0.56 for more than half of the proteins. Our method also performs well for classifying residues (rigid vs. flexible). For almost all predicted B-factor thresholds, prediction accuracies (percent of correctly predicted residues) are greater than 70%. These results exceed the best results of other sequence-based prediction methods. (C) 2005 Wiley-Liss, Inc.

Mechanisms influencing levitation and the scaling laws in nanopores: Oscillator model theory

We provide here a detailed theoretical explanation of the floating molecule or levitation effect, for molecules diffusing through nanopores, using the oscillator model theory (Phys. Rev. Lett. 2003, 91, 126102) recently developed in this laboratory. It is shown that on reduction of pore size the effect occurs due to decrease in frequency of wall collision of diffusing particles at a critical pore size. This effect is, however, absent at high temperatures where the ratio of kinetic energy to the solid-fluid interaction strength is sufficiently large. It is shown that the transport diffusivities scale with this ratio. Scaling of transport diffusivities with respect to mass is also observed, even in the presence of interactions.

Direct numerical simulation of turbulent flow and mass transfer around a rotating circular cylinder

Turbulent flow around a rotating circular cylinder has numerous applications including wall shear stress and mass-transfer measurement related to the corrosion studies. It is also of interest in the context of flow over convex surfaces where standard turbulence models perform poorly. The main purpose of this paper is to elucidate the basic turbulence mechanism around a rotating cylinder at low Reynolds numbers to provide a better understanding of flow fundamentals. Direct numerical simulation (DNS) has been performed in a reference frame rotating at constant angular velocity with the cylinder. The governing equations are discretized by using a finite-volume method. As for fully developed channel, pipe, and boundary layer flows, a laminar sublayer, buffer layer, and logarithmic outer region were observed. The level of mean velocity is lower in the buffer and outer regions but the logarithmic region still has a slope equal to the inverse of the von Karman constant. Instantaneous flow visualization revealed that the turbulence length scale typically decreases as the Reynolds number increases. Wavelet analysis provided some insight into the dependence of structural characteristics on wave number. The budget of the turbulent kinetic energy was computed and found to be similar to that in plane channel flow as well as in pipe and zero pressure gradient boundary layer flows. Coriolis effects show as an equivalent production for the azimuthal and radial velocity fluctuations leading to their ratio being lowered relative to similar nonrotating boundary layer flows.

Explicit derivation of the relativistic mass-energy relation for internal kinetic and potential energies of a composite system

A straightforward derivation of relativistic expressions for the mechanical momentum, kinetic and total energies, and mass-energy equivalence (including potential energy) which does not require any knowledge of the energy-momentum relation for electromagnetic waves or consideration of elastic collisions, but is directly based on Newton's second law and Lorentz's transformations, is presented in this paper. The existence of an invariant force is shown to be important for the validity of the relativistic mechanics.

Stoichiometric and kinetic characterisation of Nitrobacter in mixed culture by decoupling the growth and energy generation processes

The growth, maintenance and lysis processes of Nitrobacter were characterised. A Nitrobacter culture was enriched in a sequencing batch reactor (SBR). Fluorescent in situ hybridisation showed that Nitrobacter constituted 73% of the bacterial population. Batch tests were carried out to measure the oxygen uptake rate and/or nitrite consumption rate when both nitrite and CO2 were in excess, and in the absence of either of these two substrates. The results obtained, along with the SBR performance data, allowed the determination of the maintenance coefficient and in situ cell lysis rate of Nitrobacter. Nitrobacter spends a significant amount of energy for maintenance, which varies considerably with the specific growth rate. At maximum growth, Nitrobacter consume nitrite at a rate of 0.042 mgN/mgCOD(biomass)center dot h for maintenance purposes, which increases more than threefold to 0.143 mgN/mgCOD(biomass)center dot h in the absence of growth. In the SBR, where Nitrobacter grew at 40% of its maximum growth rate, a maintenance coefficient of 0.113 mgN/mgCOD center dot h was found, resulting in 42% of the total amount of nitrite being consumed for maintenance. The above three maintenance coefficient values obtained at different growth rates appear to support the maintenance model proposed in Pirt (1982). The in situ lysis rate of Nitrobacter was determined to be 0.07/day under aerobic conditions at 22 C and pH 7.3. Further, the maximum specific growth rate of Nitrobacter was estimated to be 0.02/h (0.48/day). The affinity constant of Nitrobacter with respect to nitrite was determined to be 1.50 mgNO(2)(-)-N/L, independent of the presence or absence of CO2. (c) 2006 Wiley Periodicals, Inc.

Stoichiometric and kinetic characterisation of Nitrosomonas sp in mixed culture by decoupling the growth and energy generation processes

A novel method that relies on the decoupling of the energy production and biosynthesis processes was used to characterise the maintenance, cell lysis and growth processes of Nitrosomonas sp. A Nitrosolnonas culture was enriched in a sequencing batch reactor (SBR) with ammonium as the sole energy source. Fluorescent in situ hybridization (FISH) showed that Nitrosomonas bound to the NEU probe constituted 82% of the bacterial population, while no other known ammonium or nitrite oxidizing bacteria were detected. Batch tests were carried out under conditions that both ammonium and CO, were in excess, and in the absence of one of these two substrates. The oxygen uptake rate and nitrite production rate were measured during these batch tests. The results obtained from these batch tests, along with the SBR performance data, allowed the determination of the maintenance coefficient and the in situ cell lysis rate, as well as the maximum specific growth rate of the Nitrosomonas culture. It is shown that, during normal growth, the Nitrosomonas culture spends approximately 65% of the energy generated for maintenance. The maintenance coefficient was determined to be 0.14 - 0.16 mgN mgCOD(biomass)(-1) h(-1), and was shown to be independent of the specific growth rate. The in situ lysis rate and the maximum specific growth rate of the Nitrosomonas culture were determined to be 0.26 and 1.0 day(-1) (0.043 h(-1)), respectively, under aerobic conditions at 30 degrees C and pH7. (c) 2006 Elsevier B.V. All rights reserved.

Kinetic study of the thermal degradation of high density polyethylene

The thermal degradation of high density polyethylene has been modelled by the random breakage of polymer bonds, using a set of population balance equations. A model was proposed in which the population balances were lumped into representative sizes so that the experimentally determined molecular weight distribution of the original polymer could be used as the initial condition. This model was then compared to two different cases of the unlumped population balance which assumed unimolecular initial distributions of 100 and 500 monomer units, respectively. The model that utilised the experimentally determined molecular weight distribution was found to best describe the experimental data. The model fits suggested a second mechanism in addition to random breakage at slow reaction rates. (c) 2005 Elsevier Ltd. All rights reserved.