A family of Runge-Kutta based explicit methods for rotational dynamics

The present paper develops a family of explicit algorithms for rotational dynamics and presents their comparison with several existing methods. For rotational motion the configuration space is a non-linear manifold, not a Euclidean vector space. As a consequence the rotation vector and its time derivatives correspond to different tangent spaces of rotation manifold at different time instants. This renders the usual integration algorithms for Euclidean space inapplicable for rotation. In the present algorithms this problem is circumvented by relating the equation of motion to a particular tangent space. It has been accomplished with the help of already existing relation between rotation increments which belongs to two different tangent spaces. The suggested method could in principle make any integration algorithm on Euclidean space, applicable to rotation. However, the present paper is restricted only within explicit Runge-Kutta enabled to handle rotation. The algorithms developed here are explicit and hence computationally cheaper than implicit methods. Moreover, they appear to have much higher local accuracy and hence accurate in predicting any constants of motion for reasonably longer time. The numerical results for solutions as well as constants of motion, indicate superior performance by most of our algorithms, when compared to some of the currently known algorithms, namely ALGO-C1, STW, LIEMID[EA], MCG, SUBCYC-M.

Pseudo-time particle filtering for diffuse optical tomography

We recast the reconstruction problem of diffuse optical tomography (DOT) in a pseudo-dynamical framework and develop a method to recover the optical parameters using particle filters, i.e., stochastic filters based on Monte Carlo simulations. In particular, we have implemented two such filters, viz., the bootstrap (BS) filter and the Gaussian-sum (GS) filter and employed them to recover optical absorption coefficient distribution from both numerically simulated and experimentally generated photon fluence data. Using either indicator functions or compactly supported continuous kernels to represent the unknown property distribution within the inhomogeneous inclusions, we have drastically reduced the number of parameters to be recovered and thus brought the overall computation time to within reasonable limits. Even though the GS filter outperformed the BS filter in terms of accuracy of reconstruction, both gave fairly accurate recovery of the height, radius, and location of the inclusions. Since the present filtering algorithms do not use derivatives, we could demonstrate accurate contrast recovery even in the middle of the object where the usual deterministic algorithms perform poorly owing to the poor sensitivity of measurement of the parameters. Consistent with the fact that the DOT recovery, being ill posed, admits multiple solutions, both the filters gave solutions that were verified to be admissible by the closeness of the data computed through them to the data used in the filtering step (either numerically simulated or experimentally generated). (C) 2011 Optical Society of America

Pricing a Shared Access Link for Fair and Efficient Operation with Variable User Data Rates

The literature on pricing implicitly assumes an "infinite data" model, in which sources can sustain any data rate indefinitely. We assume a more realistic "finite data" model, in which sources occasionally run out of data; this leads to variable user data rates. Further, we assume that users have contracts with the service provider, specifying the rates at which they can inject traffic into the network. Our objective is to study how prices can be set such that a single link can be shared efficiently and fairly among users in a dynamically changing scenario where a subset of users occasionally has little data to send. User preferences are modelled by concave increasing utility functions. Further, we introduce two additional elements: a convex increasing disutility function and a convex increasing multiplicative congestion-penally function. The disutility function takes the shortfall (contracted rate minus present rate) as its argument, and essentially encourages users to send traffic at their contracted rates, while the congestion-penalty function discourages heavy users from sending excess data when the link is congested. We obtain simple necessary and sufficient conditions on prices for fair and efficient link sharing; moreover, we show that a single price for all users achieves this. We illustrate the ideas using a simple experiment.

Fluctuation quench and hydrophobic collapse of polymers and biopolymers in water-DMSO binary mixture at low DMSO concentration

Binary mixtures have strong influence on activities of polymers and biopolymers even at low cosolvent concentration. Among the several aqueous binary mixtures studied, water-DMSO especially stands out for its unusual behavior at certain specific concentrations of DMSO. In the present work, we study the effect of water-DMSO binary mixture on polymers and biopolymers by taking a simple linear hydrocarbon chain of intermediate length (n = 30) and the protein lysozyme, respectively. We find that at a mole fraction of 0.05 of DMSO (x(DMSO) = 0.05) in aqueous solution, the hydrocarbon chain adopts the collapsed conformation as the most stable and rigid state. In this case of 0.05 mole fraction of DMSO in bulk, the DMSO concentration in the first hydration layer around the polymer is found to be as large as 17%. Formation of such hydrophobic environment around the polymer is the reason for the collapsed state gaining so much stability. Interestingly, similar quench of conformational fluctuation is also observed for the protein investigated. It is observed that in the case of alkane polymer chains, long wavelength fluctuation gets easily quenched, the polymer being purely hydrophobic. However, in case of the protein, quench of fluctuation is prominent only at the hydrophobic surface, and quench of long wavelength fluctuation becomes insignificant for the full protein. As protein contains both hydrophobic and hydrophilic moieties, the extent of quench of conformational fluctuation with respect to that in pure water is almost half for the biopolymer complex (16.83%) than the same for pure hydrophobic polymer chain (32.43%).

Resource Pricing for Fair and Efficient Link Sharing with a Finite Data Model

The literature on pricing implicitly assumes an "infinite data" model, in which sources can sustain any data rate indefinitely. We assume a more realistic "finite data" model, in which sources occasionally run out of data. Further, we assume that users have contracts with the service provider, specifying the rates at which they can inject traffic into the network. Our objective is to study how prices can be set such that a single link can be shared efficiently and fairly among users in a dynamically changing scenario where a subset of users occasionally has little data to send. We obtain simple necessary and sufficient conditions on prices such that efficient and fair link sharing is possible. We illustrate the ideas using a simple example

Size and temperature dependent stability and phase transformation in single-crystal zirconium nanowire

A novel size dependent FCC (face-centered-cubic) -> HCP (hexagonally-closed-pack) phase transformation and stability of an initial FCC zirconium nanowire are studied. FCC zirconium nanowires with cross-sectional dimensions < 20 are found unstable in nature, and they undergo a FCC -> HCP phase transformation, which is driven by tensile surface stress induced high internal compressive stresses. FCC nanowire with cross-sectional dimensions > 20 , in which surface stresses are not enough to drive the phase transformation, show meta-stability. In such a case, an external kinetic energy in the form of thermal heating is required to overcome the energy barrier and achieve FCC -> HCP phase transformation. The FCC-HCP transition pathway is also studied using Nudged Elastic Band (NEB) method, to further confirm the size dependent stability/metastability of Zr nanowires. We also show size dependent critical temperature, which is required for complete phase transformation of a metastable-FCC nanowire.

Spectral Finite Element Modeling of Complex Structures for Applications in Real-time Structural Health Monitoring

In this paper we discuss the recent progresses in spectral finite element modeling of complex structures and its application in real-time structural health monitoring system based on sensor-actuator network and near real-time computation of Damage Force Indicator (DFI) vector. A waveguide network formalism is developed by mapping the original variational problem into the variational problem involving product spaces of 1D waveguides. Numerical convergence is studied using a h()-refinement scheme, where is the wavelength of interest. Computational issues towards successful implementation of this method with SHM system are discussed.

Control of Molecular Weight and Polydispersity of Hyperbranched Polymers Using a Reactive B(3) Core: A Single-Step Route to Orthogonally Functionalizable Hyperbranched Polymers

Molecular weight and polydispersity are two structural features of hyperbranched polymers that are difficult to control because of the statistical nature of the step-growth polycondensation of AB(2) type monomers; the statistical growth also causes the polydispersity index to increase with percent conversion (or molecular weight). We demonstrate that using controlled amounts of a specifically designed B(3) core, containing B-type functionality that are more reactive than those present in the AB(2) monomer, both the molecular weight and the polydispersity can be readily controlled; the PDI was shown to improve with increasing mole-fraction of the B(3) core while the polymer molecular weight showed an expected decrease. Incorporation of a ``clickable'' propargyl group in the B(3) core unit permitted the generation of a core-functionalizable hyperbranched polymer. Importantly, this clickable core, in combination with a recently developed AB(2) monomer, wherein the B-type groups are allyl ethers and A is an hydroxyl group, led to the generation of a hyperbranched polymer carrying orthogonally functionalizable core and peripheral groups, via a single-step melt polycondensation. Selective functionalization of the core and periphery using two different types of chromophores was achieved, and the occurrence of fluorescence resonance energy transfer (FRET) between the donor and acceptor chromophores was demonstrated.

Control of Molecular Weight and Polydispersity of Hyperbranched Polymers Using a Reactive B(3) Core: A Single-Step Route to Orthogonally Functionalizable Hyperbranched Polymers

Molecular weight and polydispersity are two structural features of hyperbranched polymers that are difficult to control because of the statistical nature of the step-growth polycondensation of AB(2) type monomers; the statistical growth also causes the polydispersity index to increase with percent conversion (or molecular weight). We demonstrate that using controlled amounts of a specifically designed B(3) core, containing B-type functionality that are more reactive than those present in the AB(2) monomer, both the molecular weight and the polydispersity can be readily controlled; the PDI was shown to improve with increasing mole-fraction of the B(3) core while the polymer molecular weight showed an expected decrease. Incorporation of a ``clickable'' propargyl group in the B(3) core unit permitted the generation of a core-functionalizable hyperbranched polymer. Importantly, this clickable core, in combination with a recently developed AB(2) monomer, wherein the B-type groups are allyl ethers and A is an hydroxyl group, led to the generation of a hyperbranched polymer carrying orthogonally functionalizable core and peripheral groups, via a single-step melt polycondensation. Selective functionalization of the core and periphery using two different types of chromophores was achieved, and the occurrence of fluorescence resonance energy transfer (FRET) between the donor and acceptor chromophores was demonstrated.

Microstructure and texture evolution during accumulative roll bonding of aluminium alloy AA5086

In the present investigation, a strongly bonded strip of an aluminium-magnesium based alloy AA5086 is successfully produced through accumulative roll bonding (ARB). A maximum of up to eight passes has been used for the purpose. Microstructural characterization using electron backscatter diffraction (EBSD) technique indicates the formation of submicron sized (similar to 200-300 nm) subgrains inside the layered microstructure. The material is strongly textured where individual layers possess typical FCC rolling texture components. More than three times enhancement in 0.2% proof stress (PS) has been obtained after 8 passes due to grain refinement and strain hardening. (C) 2011 Elsevier B.V. All rights reserved.

Crystallization and preliminary X-ray diffraction studies of sortase A from Streptococcus pneumoniae

Sortases are cell-membrane-anchored cysteine transpeptidases that are essential for the assembly and anchoring of cell-surface adhesins in Gram-positive bacteria. Thus, they play critical roles in virulence, infection and colonization by pathogens. Sortases have been classified into four types based on their primary sequence and the target-protein motifs that they recognize. All Gram-positive bacteria express a class A housekeeping sortase (SrtA). Sortase A from Streptococcus pneumoniae (NP_358691) has been crystallized in two crystal forms. Diamond-shaped crystals of Delta N(59)SrtA diffracted to 4.0 angstrom resolution and belonged to a tetragonal system with unit-cell parameters a = b = 122.8, c = 86.5 angstrom, alpha = beta = gamma = 90 degrees, while rod-shaped crystals of Delta N(81)SrtA diffracted to 2.91 angstrom resolution and belonged to the monoclinic space group P2(1) with unit-cell parameters a = 66.8, b = 103.47, c = 74.79 angstrom, alpha = gamma = 90, beta = 115.65 degrees. The Matthews coefficient (V(M) = 2.77 angstrom(3) Da(-1)) with similar to 56% solvent content suggested the presence of four molecules in the asymmetric unit for Delta N(81)SrtA. Also, a multi-copy search using a monomer as a probe in the molecular-replacement method resulted in the successful location of four sortase molecules in the asymmetric unit, with statistics R = 41.61, R(free) = 46.44, correlation coefficient (CC) = 64.31, CC(free) = 57.67.

Ultrasound modulated optical tomography: Young's modulus of the insonified region from measurement of natural frequency of vibration

We demonstrate a method to recover the Young's modulus (E) of a tissue-mimicking phantom from measurements of ultrasound modulated optical tomography (UMOT). The object is insonified by a dualbeam, confocal ultrasound transducer (US) oscillating at frequencies f(0) and f(0) + Delta f and the variation of modulation depth (M) in the autocorrelation of light traversed through the focal region of the US transducer against Delta f is measured. From the dominant peaks observed in the above variation, the natural frequencies of the insonified region associated with the vibration along the US transducer axis are deduced. A consequence of the above resonance is that the speckle fluctuation at the resonance frequency has a higher signal-to-noise to ratio (SNR). From these natural frequencies and the associated eigenspectrum of the oscillating object, Young's modulus (E) of the material in the focal region is recovered. The working of this method is confirmed by recovering E in the case of three tissue-mimicking phantoms of different elastic modulus values. (C) 2011 Optical Society of America

Optimised form of acceleration correction algorithm within SPH-based simulations of impact mechanics

In the context of SPH-based simulations of impact dynamics, an optimised and automated form of the acceleration correction algorithm (Shaw and Reid, 2009a) is developed so as to remove spurious high frequency oscillations in computed responses whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. A rational framework for an insightful characterisation of the erstwhile acceleration correction method is first set up. This is followed by the proposal of an optimised version of the method, wherein the strength of the correction term in the momentum balance and energy equations is optimised. For the first time, this leads to an automated procedure to arrive at the artificial viscosity term. In particular, this is achieved by taking a spatially varying response-dependent support size for the kernel function through which the correction term is computed. The optimum value of the support size is deduced by minimising the (spatially localised) total variation of the high oscillation in the acceleration term with respect to its (local) mean. The derivation of the method, its advantages over the heuristic method and issues related to its numerical implementation are discussed in detail. (C) 2011 Elsevier Ltd. All rights reserved.