Conservation Properties of the Trapezoidal Rule in Linear Time Domain Analysis of Acoustics and Structures

The trapezoidal rule, which is a special case of the Newmark family of algorithms, is one of the most widely used methods for transient hyperbolic problems. In this work, we show that this rule conserves linear and angular momenta and energy in the case of undamped linear elastodynamics problems, and an ``energy-like measure'' in the case of undamped acoustic problems. These conservation properties, thus, provide a rational basis for using this algorithm. In linear elastodynamics problems, variants of the trapezoidal rule that incorporate ``high-frequency'' dissipation are often used, since the higher frequencies, which are not approximated properly by the standard displacement-based approach, often result in unphysical behavior. Instead of modifying the trapezoidal algorithm, we propose using a hybrid finite element framework for constructing the stiffness matrix. Hybrid finite elements, which are based on a two-field variational formulation involving displacement and stresses, are known to approximate the eigenvalues much more accurately than the standard displacement-based approach, thereby either bypassing or reducing the need for high-frequency dissipation. We show this by means of several examples, where we compare the numerical solutions obtained using the displacement-based and hybrid approaches against analytical solutions.

Conformational analysis of hydroxyproline and its protected derivatives by 1H NMR

From a computer simulation of the 270 MHz 1H NMR spectra of hydroxyproline (Hyp) and its protected derivatives, precise values of ring vicinal coupling constants were obtained. These couplings were related to ring torsional angles, using a Karplus type analysis. From the NMR analysis it was observed that the pyrrolidine ring possesses a unique and highly homogeneous conformation (Cγ-exo form). Temperature dependence studies on protected dipeptides suggest that the pyrrolidine ring conformation is independent of backbone conformation. An unusual X-Hyp, β-turn was observed for Boc-Aib-Hyp-NHMe.

Genetic analysis of foot-and-mouth disease virus serotype A of Indian origin and detection of positive selection and recombination in leader protease- and capsid-coding regions

The leader protease (L-pro) and capsid-coding sequences (P1) constitute approximately 3 kb of the foot-and-mouth disease virus (FMDV). We studied the phylogenetic relationship of 46 FMDV serotype A isolates of Indian origin collected during the period 1968-2005 and also eight vaccine strains using the neighbour-joining tree and Bayesian tree methods. The viruses were categorized under three major groups - Asian, Euro-South American and European. The Indian isolates formed a distinct genetic group among the Asian isolates. The Indian isolates were further classified into different genetic subgroups (<5% divergence). Post-1995 isolates were divided into two subgroups while a few isolates which originated in the year 2005 from Andhra Pradesh formed a separate group. These isolates were closely related to the isolates of the 1970s. The FMDV isolates seem to undergo reverse mutation or onvergent evolution wherein sequences identical to the ancestors are present in the isolates in circulation. The eight vaccine strains included in the study were not related to each other and belonged to different genetic groups. Recombination was detected in the L-pro region in one isolate (A IND 20/82) and in the VP1 coding 1D region in another isolate (A RAJ 21/96). Positive selection was identified at aa positions 23 in the L-pro (P<0.05; 0.046*) and at aa 171 in the capsid protein VP1 (P<0.01; 0.003**).

Singular perturbation analysis of spherical and cylindrical N waves

Using a singular perturbation analysis the nonplanar Burgers' equation is solved to yield the shock wave-displacement due to diffusion for spherical and cylindrical N waves, thus supplementing the earlier results of Lighthill for the plane N waves. Physics of Fluids is copyrighted by The American Institute of Physics.

Analysis of microstructure and texture evolution in pure magnesium during symmetric and asymmetric rolling

Asymmetric rolling of commercially pure magnesium was carried out at three different temperatures: room temperature, 200 degrees C and 350 degrees C. Systematic analysis of microstructures, grain size distributions, texture and misorientation distributions were performed using electron backscattered diffraction in a field emission gun scanning electron microscope. The results were compared with conventional (symmetric) rolling carried out under the same conditions of temperature and strain rate. Simulations of deformation texture evolution were performed using the viscoplastic self-consistent polycrystal plasticity model. The main trends of texture evolution are faithfully reproduced by the simulations for the tests at room temperature. The deviations that appear for the textures obtained at high temperature can be explained by the occurrence of dynamic recrystallization. Finally, the mechanisms of texture evolution in magnesium during asymmetric and symmetric rolling are explained with the help of ideal orientations, grain velocity fields and divergence maps displayed in orientation space.

Free vibration and wave propagation analysis of uniform and tapered rotating beams using spectrally formulated finite elements

This paper presents a formulation of an approximate spectral element for uniform and tapered rotating Euler-Bernoulli beams. The formulation takes into account the varying centrifugal force, mass and bending stiffness. The dynamic stiffness matrix is constructed using the weak form of the governing differential equation in the frequency domain, where two different interpolating functions for the transverse displacement are used for the element formulation. Both free vibration and wave propagation analysis is performed using the formulated elements. The studies show that the formulated element predicts results, that compare well with the solution available in the literature, at a fraction of the computational effort. In addition, for wave propagation analysis, the element shows superior convergence. (C) 2007 Elsevier Ltd. All rights reserved.

Analysis of dynamics and mechanism of ligand binding to Artocarpus integrifolia agglutinin. A 13C and 19F NMR study

Binding of 13C-labeled N-acetylgalactosamine (13C-GalNAc) and N-trifluoroacetylgalactosamine (19F-GalNAc) to Artocarpus integrifolia agglutinin has been studied using 13C and 19F nuclear magnetic resonance spectroscopy, respectively. Binding of these saccharides resulted in broadening of the resonances, and no change in chemical shift was observed, suggesting that the alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc experience a magnetically equivalent environment in the lectin combining site. The alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc were found to be in slow exchange between free and protein bound states. Binding of 13C-GalNAc was studied as a function of temperature. From the temperature dependence of the line broadening, the thermodynamic and kinetic parameters were evaluated. The association rate constants obtained for the alpha-anomers of 13C-GalNAc and 19F-GalNAc (k+1 = 1.01 x 10(5) M-1.s-1 and 0.698 x 10(5) M-1.s-1, respectively) are in close agreement with those obtained for the corresponding beta-anomers (k+1 = 0.95 x 10(5) M-1.s-1 and 0.65 x 10(5) M-1.s-1, respectively), suggesting that the two anomers bind to the lectin by a similar mechanism. In addition these values are several orders of magnitude slower than those obtained for diffusion controlled processes. The dissociation rate constants obtained are 49.9, 56.9, 42, and 43 s-1, respectively, for the alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc. A two-step mechanism has been proposed for the interaction of 13C-GalNAc and 19F-GalNAc with A. integrifolia lectin in view of the slow association rates and high activation entropies. The thermodynamic parameters obtained for the association and dissociation reactions suggest that the binding process is entropically favored and that there is a small enthalpic contribution.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Analysis of gradually and spatially varied flow in sand-bed channels

Seepage through a sand bed affects the channel hydrodynamics, which in turn alters channel stability. Thus, the effect of seepage on its hydrodynamic parameters needs to be ascertained. The present work analyses spatially varied flow of a sand-bed channel subjected to seepage in the downward direction through a sand bed. Numerically calculated flow profiles affected by seepage have been verified using experimental observations. The present work also analyses the friction slope, velocity and bed shear stress variations along the channel for both seepage and no-seepage conditions. It was found that the downward seepage-induced channel flow has larger friction slope and bed shear stress than that of no-seepage.

Analysis of heat and mass transfer during bulk hydraircooling of spherical food products

Hydraircooling is a technique used for precooling food products. In this technique chilled water is sprayed over the food products while cold unsaturated air is blown over them. Hydraircooling combines the advantages of both air- and hydrocooling. The present study is concerned with the analysis of bulk hydraircooling as it occurs in a package filled with several layers of spherical food products with chilled water sprayed from the top and cold unsaturated air blown from the bottom. A mathematical model is developed to describe the hydrodynamics and simultaneous heat and mass transfer occurring inside the package. The non-dimensional governing equations are solved using the finite difference numerical methods. The results are presented in the form of time-temperature charts. A correlation is obtained to calculate the process time in terms of the process parameters.

Numerical Electromagnetic Analysis of Current and Voltage Distribution along Down Conductors and Towers Hit by Direct Lightning

In order to answer the practically important question of whether the down conductors of lightning protection systems to tall towers and buildings can be electrically isolated from the structure itself, this work is conducted. As a first step in this regard, it is presumed that the down conductor placed on metallic tower will be a pessimistic representation of the actual problem. This opinion was based on the fact that the proximity of heavy metallic structure will have a large damping effect. The post-stroke current distributions along the down conductors and towers, which can be quite different from that in the lightning channel, govern the post-stroke near field and the resulting gradient in the soil. Also, for a reliable estimation of the actual stroke current from the measured down conductor currents, it is essential to know the current distribution characteristics along the down conductors. In view of these, the present work attempts to deduce the post-stroke current and voltage distribution along typical down conductors and towers. A solution of the governing field equations on an electromagnetic model of the system is sought for the investigation. Simulation results providing the spatio-temporal distribution of the post-stroke current and voltage has provided very interesting results. It is concluded that it is almost impossible to achieve electrical isolation between the structure and the down conductor. Furthermore, there will be significant induction into the steel matrix of the supporting structure.

Theoretical and Computational Analysis of Static and Dynamic Anomalies in Water-DMSO Binary Mixture at Low DMSO Concentrations

Experiments have repeatedly observed both thermodynamic and dynamic anomalies in aqueous binary mixtures, surprisingly at low solute concentration. Examples of such binary mixtures include water-DMSO, water-ethanol, water-tertiary butyl alcohol (TBA), and water-dioxane, to name a few. The anomalies have often been attributed to the onset of a structural transition, whose nature, however, has been left rather unclear. Here we study the origin of such anomalies using large scale computer simulations and theoretical analysis in water-DMSO binary mixture. At very low DMSO concentration (below 10%), small aggregates of DMSO are solvated by water through the formation of DMSO-(H2O)(2) moieties. As the concentration is increased beyond 10-12% of DMSO, spanning clusters comprising the same moieties appear in the system. Those clusters are formed and stabilized not only through H-bonding but also through the association of CH3 groups of DMSO. We attribute the experimentally observed anomalies to a continuum percolation-like transition at DMSO concentration X-DMSO approximate to 12-15%. The largest cluster size of CH3-CH3 aggregation clearly indicates the formation of such percolating clusters. As a result, a significant slowing down is observed in the decay of associated rotational auto time correlation functions (of the S = O bond vector of DMSO and O-H bond vector of water). Markedly unusual behavior in the mean square fluctuation of total dipole moment again suggests a structural transition around the same concentration range. Furthermore, we map our findings to an interacting lattice model which substantiates the continuum percolation model as the reason for low concentration anomalies in binary mixtures where the solutes involved have both hydrophilic and hydrophobic moieties.

Analysis of geometric and strain effects in homo-Diels-Alder reactions

Molecular mechanics calculations have been carried out to quantify the key geometric and strain effects which are likely to control the homo-Diels-Alder reactivity of 1,4-dienes. The criteria considered include C1..C5 and C2..C4 distances in the diene, twist angle of the two pi units, and the magnitude of strain increase as a result of cycloaddition. By first considering these factors in a number of non-conjugated dienes with known reactivity, the ranges of values within which the reaction is favoured are proposed. Calculations are also reported on several substrates which have not been investigated so far. Promising systems for experimental study are suggested which, in addition to being intrinsically interesting, would place the present proposals on a firm basis.

Plane wave analysis of conical and exponential pipes with incompressible mean flow

The general equation for one-dimensional wave propagation at low flow Mach numbers (M less-than-or-equals, slant0·2) is derived and is solved analytically for conical and exponential shapes. The transfer matrices are derived and shown to be self-consistent. Comparison is also made with the relevant data available in the literature. The transmission loss behaviour of conical and exponential pipes, and mufflers involving these shapes, are studied. Analytical expressions of the same are given for the case of a stationary medium. The mufflers involving conical and exponential pipes are shown to be inferior to simple expansion chambers (of similar dimensions) at higher frequencies from the point of view of noise abatement, as was observed earlier experimentally.