Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNPO4 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity theta and the wavelength lambda of a plane wave; we show that PD leads to a periodic, spatial modulation of theta and a temporally periodic modulation of lambda; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNPO4 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNPO4 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Optical diagnostics of fuel-air mixing and vortex formation in a cavity combustor

The current work reports optical diagnostic measurements of fuel-air mixing and vortex structure in a single cavity trapped vortex combustor (TVC). Specifically, the mixture fraction using acetone PLIF technique in the non-reacting flow, and PIV measurements in the reacting flow are reported for the first time in trapped vortex combustors. The fuel-air momentum flux ratio, where the air momentum corresponds to that entering the cavity through a specially-incorporated flow guide vane, is used to characterize the mixing. The acetone PLIF experiments show that at high momentum flux ratios, the fuel-air mixing in the cavity is very minimal and is enhanced as the momentum flux ratio reduces, due to a favourable vortex formation in the cavity. Stoichiometric mixture fraction surfaces show that the mixing causes the reaction surfaces to shift from non-premixed to partially-premixed stratified mixtures. PIV measurements conducted in the non-reacting flow in the cavity further reinforce this observation. The scalar dissipation rates of mixture fraction were compared with the contours of RMS of fluctuating velocity and showed very good agreement. The regions of maximum mixing are observed to be along the fuel air interface. Reacting flow Ply measurements which differ substantially from the non-reacting cases primarily because of the heat release from combustion and the resulting gas expansion show that the vortex is displaced from the centre of the cavity towards the guide vane. Overall, the measurements show interesting features of the flow including the presence of the dual cavity structure and lead to a clear understanding of the underlying physics of the cavity flow highlighting the importance of the fuel-air momentum ratio parameter. (C) 2014 Elsevier Inc. All rights reserved.

Limiting two-Higgs-doublet models

We update the constraints on two-Higgs-doublet models (2HDMs) focusing on the parameter space relevant to explain the present muon g - 2 anomaly, Delta alpha(mu), in four different types of models, type I, II, ``lepton specific'' (or X) and ``flipped'' (or Y). We show that the strong constraints provided by the electroweak precision data on the mass of the pseudoscalar Higgs, whose contribution may account for Delta alpha(mu), are evaded in regions where the charged scalar is degenerate with the heavy neutral one and the mixing angles alpha and beta satisfy the Standard Model limit beta - alpha approximate to pi/2. We combine theoretical constraints from vacuum stability and perturbativity with direct and indirect bounds arising from collider and B physics. Possible future constraints from the electron g - 2 are also considered. If the 126 GeV resonance discovered at the LHC is interpreted as the light CP-even Higgs boson of the 2HDM, we find that only models of type X can satisfy all the considered theoretical and experimental constraints.

Response of shock wave deformation in AA5086 aluminum alloy

The AA5086 aluminum alloy sheets with different starting textures were subjected to shock wave deformation with an input impulse of similar to 0.2 Ns. Microstructural examination indicate no significant change in grain size; however, the evolution of substructure manifesting intra-granular misorientation was evident. The improvement in hardness indicates the absence of recovery and strain hardening during shock deformation. Shock deformed samples show characteristic texture evolution with high Brass {110}< 112 > component. The study demonstrates the viability of high velocity forming of AA5086 aluminum alloy sheet using shock wave. (C) 2014 Elsevier B.V. All rights reserved.

Wave propagation analysis in adhesively bonded composite joints using the wavelet spectral finite element method

A wavelet spectral finite element (WSFE) model is developed for studying transient dynamics and wave propagation in adhesively bonded composite joints. The adherands are formulated as shear deformable beams using the first order shear deformation theory (FSDT) to obtain accurate results for high frequency wave propagation. Equations of motion governing wave motion in the bonded beams are derived using Hamilton's principle. The adhesive layer is modeled as a line of continuously distributed tension/compression and shear springs. Daubechies compactly supported wavelet scaling functions are used to transform the governing partial differential equations from time domain to frequency domain. The dynamic stiffness matrix is derived under the spectral finite element framework relating the nodal forces and displacements in the transformed frequency domain. Time domain results for wave propagation in a lap joint are validated with conventional finite element simulations using Abaqus. Frequency domain spectrum and dispersion relation results are presented and discussed. The developed WSFE model yields efficient and accurate analysis of wave propagation in adhesively-bonded composite joints. (C) 2014 Elsevier Ltd. All rights reserved.

Diurnal and seasonal variation of mixing ratio and delta C-13 of air CO2 observed at an urban station Bangalore, India

We present here observations on diurnal and seasonal variation of mixing ratio and delta C-13 of air CO2, from an urban station-Bangalore (BLR), India, monitored between October 2008 and December 2011. On a diurnal scale, higher mixing ratio with depleted delta C-13 of air CO2 was found for the samples collected during early morning compared to the samples collected during late afternoon. On a seasonal scale, mixing ratio was found to be higher for dry summer months (April-May) and lower for southwest monsoon months (June-July). The maximum enrichment in delta C-13 of air CO2 (-8.04 +/- 0.02aEuro degrees) was seen in October, then delta C-13 started depleting and maximum depletion (-9.31 +/- 0.07aEuro degrees) was observed during dry summer months. Immediately after that an increasing trend in delta C-13 was monitored coincidental with the advancement of southwest monsoon months and maximum enrichment was seen again in October. Although a similar pattern in seasonal variation was observed for the three consecutive years, the dry summer months of 2011 captured distinctly lower amplitude in both the mixing ratio and delta C-13 of air CO2 compared to the dry summer months of 2009 and 2010. This was explained with reduced biomass burning and increased productivity associated with prominent La Nina condition. While compared with the observations from the nearest coastal and open ocean stations-Cabo de Rama (CRI) and Seychelles (SEY), BLR being located within an urban region captured higher amplitude of seasonal variation. The average delta C-13 value of the end member source CO2 was identified based on both diurnal and seasonal scale variation. The delta C-13 value of source CO2 (-24.9 +/- 3aEuro degrees) determined based on diurnal variation was found to differ drastically from the source value (-14.6 +/- 0.7aEuro degrees) identified based on seasonal scale variation. The source CO2 identified based on diurnal variation incorporated both early morning and late afternoon sample; whereas, the source CO2 identified based on seasonal variation included only afternoon samples. Thus, it is evident from the study that sampling timing is one of the important factors while characterizing the composition of end member source CO2 for a particular station. The difference in delta C-13 value of source CO2 obtained based on both diurnal and seasonal variation might be due to possible contribution from cement industry along with fossil fuel / biomass burning as predominant sources for the station along with differential meteorological conditions prevailed.

Remotely triggered micro-shock wave responsive drug delivery system for resolving diabetic wound infection and controlling blood sugar levels

A novel, micro-shock wave responsive spermidine and dextran sulfate microparticle was developed. Almost 90% of the drug release was observed when the particles were exposed to micro-shock waves 5 times. Micro-shock waves served two purposes; of releasing the antibiotic from the system and perhaps disrupting the S. aureus biofilm in the skin infection model. A combination of shock waves with ciprofloxacin loaded microparticles could completely cure the S. aureus infection lesion in a diabetic mouse model. As a proof of concept insulin release was triggered using micro-shock waves in diabetic mice to reduce the blood glucose level. Insulin release could be triggered for at least 3 days by exposing subcutaneously injected insulin loaded particles.

Wave propagation of phonon and phason displacement modes in quasicrystals: Determination of wave parameters

The paper presents the study of wave propagation in quasicrystals. Our interest is in the computation of the wavenumber (k(n)) and group speed (c(g)) of the phonon and phason displacement modes of one, two, and three dimensional quasicrystals. These wave parameter expressions are derived and computed using the elasto-hydrodynamic equations for quasicrystals. For the computation of the wavenumber and group speeds, we use Fourier transform approximation of the phonon and the phason displacement modes. The characteristic equations obtained are a polynomial equation of the wavenumber (k(n)), with frequency as a parameter. The corresponding group speeds (c(g)) for different frequencies are then computed from the wavenumber k(n). The variation of wavenumber and group speeds with frequency is plotted for the 1-D quasicrystal, 2-D decagonal Al-Ni-Co quasicrystals, and 3-D icosahedral Al-Pd-Mn and Zn-Mg-Sc quasicrystals. From the wavenumber and group speeds plots, we obtain the cut-off frequencies for different spatial wavenumber eta(m). The results show that for 1-D, 2-D, and 3-D quasicrystals, the phonon displacement modes are non-dispersive for low values of eta(m) and becomes dispersive for increasing values of eta(m). The cut-off frequencies are not observed for very low values of eta(m), whereas the cut-off frequency starts to appear with increasing eta(m). The group speeds of the phason displacement modes are orders of magnitude lower than that of the phonon displacement modes, showing that the phason modes do not propagate, and they are essentially the diffusive modes. The group speeds of the phason modes are also not influenced by eta(m). The group speeds for the 2-D quasicrystal at 35 kHz is also simulated numerically using Galerkin spectral finite element methods in frequency domain and is compared with the results obtained using wave propagation analysis. The effect of the phonon and phason elastic constants on the group speeds is studied using 3-D icosahedral Al-Pd-Mn and Zn-Mg-Sc quasicrystals. It is also shown that the phason elastic constants and the coupling coefficient do not affect the group speeds of the phonon displacement modes. (C) 2015 AIP Publishing LLC.

Optimal Pulse Width Modulation of Voltage-Source Inverter fed Motor Drives with Relaxation of Quarter Wave Symmetry Condition

Optimal switching angles for minimization of total harmonic distortion of line current (I-THD) in a voltage source inverter are determined traditionally by imposing half-wave symmetry (HWS) and quarter-wave symmetry (QWS) conditions on the pulse width modulated waveform. This paper investigates optimal switching angles with QWS relaxed. Relaxing QWS expands the solution space and presents the possibility of improved solutions. The optimal solutions without QWS are shown here to outperform the optimal solutions with QWS over a range of modulation index (M) between 0.82 and 0.94 for a switching frequency to fundamental frequency ratio of 5. Theoretical and experimental results are presented on a 2.3kW induction motor drive.

L-p estimates for the wave equation associated to the Grushin operator

We prove that the solution of the wave equation associated to the Grushin operator G = -Delta -vertical bar x vertical bar(2)partial derivative(2)(t) is bounded on L-P (Rn+1), with 1 < p < infinity, when vertical bar 1/p - 1/2 vertical bar < 1/n+2.

A spectral element for wave propagation in honeycomb sandwich construction considering core flexibility

Spectral elements are found to be extremely resourceful to study the wave propagation characteristics of structures at high frequencies. Most of the aerospace structures use honeycomb sandwich constructions. The existing spectral elements use single layer theories for a sandwich construction wherein the two face sheets vibrate together and this model is sufficient for low frequency excitations. At high frequencies, the two face sheets vibrate independently. The Extended Higher order SAndwich Plate theory (EHSaPT) is suitable for representing the independent motion of the face sheets. A 1D spectral element based on EHSaPT is developed in this work. The wave number and the wave speed characteristics are obtained using the developed spectral element. It is shown that the developed spectral element is capable of representing independent wave motions of the face sheets. The propagation speeds of a high frequency modulated pulse in the face sheets and the core of a honeycomb sandwich are demonstrated. Responses of a typical honeycomb sandwich beam to high frequency shock loads are obtained using the developed spectral element and the response match very well with the finite element results. It is shown that the developed spectral element is able to represent the flexibility of the core resulting into independent wave motions in the face sheets, for which a finite element method needs huge degrees of freedom. (C) 2015 Elsevier Ltd. All rights reserved.

Behaviour of turbulent Prandtl/Schmidt number in compressible mixing layer

The behaviour of turbulent Prandtl/Schmidt number is explored through the model-free simulation results. It has been observed that compressibility affects the Reynolds scalar flux vectors. Reduced peak values are also observed for compressible convective Mach number mixing layer as compared with the incompressible convective Mach number counterpart, indicating a reduction in the mixing of enthalpy and species. Pr-t and Sc-t variations also indicate a reduction in mixing. It is observed that unlike the incompressible case, it is difficult to assign a constant value to these numbers due to their continuous variation in space. Modelling of Pr-t and Sc-t would be necessary to cater for this continuous spatial variation. However, the turbulent Lewis number is evaluated to be near unity for the compressible case, making it necessary to model only one of the Pr-t and Sc-t..

Reynolds averaged Navier-Stokes simulations of compressible mixing layers of similar and dissimilar gases: Performance of k- turbulence model

The issue of growth rate reduction of high speed mixing layer with convective Mach number is examined for similar and dissimilar gases using Reynolds averaged Navier-Stokes (RANS) methodology with k- turbulence model. It is observed that the growth rate predicted using RANS simulations closely matches with that predicted using model free simulations. Velocity profiles do not depend on the modelled value of Pr-t and Sc-t; while the temperature and species mass fraction distributions depend heavily on them. Although basic k- turbulence model could not capture the reduced growth rate for the mixing layer formed between similar gases, it predicts very well the reduced growth rate for the mixing layer for the dissimilar gases. It appears that density ratio changes caused by temperature changes for the dissimilar gases have profound effect on the growth rate reduction.

Modeling Ultrasonic NDE and Guided Wave based Structural Health Monitoring

Structural Health Monitoring (SHM) systems require integration of non-destructive technologies into structural design and operational processes. Modeling and simulation of complex NDE inspection processes are important aspects in the development and deployment of SHM technologies. Ray tracing techniques are vital simulation tools to visualize the wave path inside a material. These techniques also help in optimizing the location of transducers and their orientation with respect to the zone of interrogation. It helps in increasing the chances of detection and identification of a flaw in that zone. While current state-of-the-art techniques such as ray tracing based on geometric principle help in such visualization, other information such as signal losses due to spherical or cylindrical shape of wave front are rarely taken into consideration. The problem becomes a little more complicated in the case of dispersive guided wave propagation and near-field defect scattering. We review the existing models and tools to perform ultrasonic NDE simulation in structural components. As an initial step, we develop a ray-tracing approach, where phase and spectral information are preserved. This enables one to study wave scattering beyond simple time of flight calculation of rays. Challenges in terms of theory and modelling of defects of various kinds are discussed. Various additional considerations such as signal decay and physics of scattering are reviewed and challenges involved in realistic computational implementation are discussed. Potential application of this approach to SHM system design is highlighted and by applying this to complex structural components such as airframe structures, SHM is demonstrated to provide additional value in terms of lighter weight and/or longevity enhancement resulting from an extension of the damage tolerance design principle not compromising safety and reliability.