Self-regularized pseudo time-marching schemes for structural system identification with static measurements

We explore the application of pseudo time marching schemes, involving either deterministic integration or stochastic filtering, to solve the inverse problem of parameter identification of large dimensional structural systems from partial and noisy measurements of strictly static response. Solutions of such non-linear inverse problems could provide useful local stiffness variations and do not have to confront modeling uncertainties in damping, an important, yet inadequately understood, aspect in dynamic system identification problems. The usual method of least-square solution is through a regularized Gauss-Newton method (GNM) whose results are known to be sensitively dependent on the regularization parameter and data noise intensity. Finite time, recursive integration of the pseudo-dynamical GNM (PD-GNM) update equation addresses the major numerical difficulty associated with the near-zero singular values of the linearized operator and gives results that are not sensitive to the time step of integration. Therefore, we also propose a pseudo-dynamic stochastic filtering approach for the same problem using a parsimonious representation of states and specifically solve the linearized filtering equations through apseudo-dynamic ensemble Kalman filter (PD-EnKF). For multiple sets ofmeasurements involving various load cases, we expedite the speed of the PD-EnKF by proposing an inner iteration within every time step. Results using the pseudo-dynamic strategy obtained through PD-EnKF and recursive integration are compared with those from the conventional GNM, which prove that the PD-EnKF is the best performer showing little sensitivity to process noise covariance and yielding reconstructions with less artifacts even when the ensemble size is small. Copyright (C) 2009 John Wiley & Sons, Ltd.

Estimation of signals in colored non gaussian noise based on Gaussian Mixture models

Non-Gaussianity of signals/noise often results in significant performance degradation for systems, which are designed using the Gaussian assumption. So non-Gaussian signals/noise require a different modelling and processing approach. In this paper, we discuss a new Bayesian estimation technique for non-Gaussian signals corrupted by colored non Gaussian noise. The method is based on using zero mean finite Gaussian Mixture Models (GMMs) for signal and noise. The estimation is done using an adaptive non-causal nonlinear filtering technique. The method involves deriving an estimator in terms of the GMM parameters, which are in turn estimated using the EM algorithm. The proposed filter is of finite length and offers computational feasibility. The simulations show that the proposed method gives a significant improvement compared to the linear filter for a wide variety of noise conditions, including impulsive noise. We also claim that the estimation of signal using the correlation with past and future samples leads to reduced mean squared error as compared to signal estimation based on past samples only.

A self-consistent formulation for analysis and generation of non-uniform DNA structure

A fully self-consistent formulation is described here for the analysis and generation of base-pairs in non-uniform DNA structures, in terms of various local parameters. It is shown that the internal "wedge parameters" are mathematically related to the parameters describing the base-pair orientation with respect to an external helix axis. Hence any one set of three translation and three rotation parameters are necessary and sufficient to completely describe the relative orientation of the base-pairs comprising a step (or doublet). A general procedure is outlined for obtaining an average or global helix axis from the local helix axes for each step. A graphical representation of the local helix axes in the form of a polar plot is also shown and its application for estimating the curvature of oligonucleotide structures is illustrated, with examples of both A and B type structures.

Molecular structure-activity relationship study of some non-steroidal antiinflammatory agents using electrostatic potential mapping

A series of 6,11-dihydro-11-oxodibenz[b,e]oxepin-2-acetic acids (DOAA) which are known to be anti-inflammatory agents were studied. The geometries of some of the molecules obtained from X-ray crystallography were used in the calculations as such while the geometries of their derivatives were obtained by local, partial geometry optimization around the Sites of substitution employing the AMI method, keeping the remaining parts of the geometries the same as those in the parent molecules. Molecular electrostatic potential (MEP) mapping was performed for the molecules using optimized hybridization displacement charges (HDC) combined with Lowdin charges, as this charge distribution has been shown earlier to yield near ab initio quality results. A good correlation has been found between the MEP values near the oxygen atoms of the hydroxyl groups of the carboxy groups of the molecules and their anti-inflammatory activities. The result is broadly in agreement with the model proposed earlier by other authors regarding the structure-activity relationship for other similar molecules.

Molecular mechanisms of hypertension-induced heart failure : Experimental studies with special emphasis on local renin-angiotensin system, cardiac metabolism and levosimendan

Hypertension is a major risk factor for stroke, ischaemic heart disease, and the development of heart failure. Hypertension-induced heart failure is usually preceded by the development of left ventricular hypertrophy (LVH), which represents an adaptive and compensatory response to the increased cardiac workload. Biomechanical stress and neurohumoral activation are the most important triggers of pathologic hypertrophy and the transition of cardiac hypertrophy to heart failure. Non-clinical and clinical studies have also revealed derangements of energy metabolism in hypertensive heart failure. The goal of this study was to investigate in experimental models the molecular mechanisms and signalling pathways involved in hypertension-induced heart failure with special emphasis on local renin-angiotensin-aldosterone system (RAAS), cardiac metabolism, and calcium sensitizers, a novel class of inotropic agents used currently in the treatment of acute decompensated heart failure. Two different animal models of hypertensive heart failure were used in the present study, i.e. hypertensive and salt-sensitive Dahl/Rapp rats on a high salt diet (a salt-sensitive model of hypertensive heart failure) and double transgenic rats (dTGR) harboring human renin and human angiotensinogen genes (a transgenic model of hypertensive heart failure with increased local RAAS activity). The influence of angiotensin II (Ang II) on cardiac substrate utilization and cardiac metabolomic profile was investigated by using gas chromatography coupled to time-of-flight mass spectrometry to detect 247 intermediary metabolites. It was found that Ang II could alter cardiac metabolomics both in normotensive and hypertensive rats in an Ang II receptor type 1 (AT1)-dependent manner. A distinct substrate use from fatty acid oxidation towards glycolysis was found in dTGR. Altered cardiac substrate utilization in dTGR was associated with mitochondrial dysfunction. Cardiac expression of the redox-sensitive metabolic sensor sirtuin1 (SIRT1) was increased in dTGR. Resveratrol supplementation prevented cardiovascular mortality and ameliorated Ang II-induced cardiac remodeling in dTGR via blood pressure-dependent pathways and mechanisms linked to increased mitochondrial biogenesis. Resveratrol dose-dependently increased SIRT1 activity in vitro. Oral levosimendan treatment was also found to improve survival and systolic function in dTGR via blood pressure-independent mechanisms, and ameliorate Ang II-induced coronary and cardiomyocyte damage. Finally, using Dahl/Rapp rats it was demonstrated that oral levosimendan as well as the AT1 receptor antagonist valsartan improved survival and prevented cardiac remodeling. The beneficial effects of levosimendan were associated with improved diastolic function without significantly improved systolic changes. These positive effects were potentiated when the drug combination was administered. In conclusion, the present study points to an important role for local RAAS in the pathophysiology of hypertension-induced heart failure as well as its involvement as a regulator of cardiac substrate utilization and mitochondrial function. Our findings suggest a therapeutic role for natural polyphenol resveratrol and calcium sensitizer, levosimendan, and the novel drug combination of valsartan and levosimendan, in prevention of hypertension-induced heart failure. The present study also provides a better understanding of the pathophysiology of hypertension-induced heart failure, and may help identify potential targets for novel therapeutic interventions.

Inhomogeneity in sintered ZnO based non-linear resistors

Investigations have been made on the inhomogeneous characteristics of sintered ZnO based non-linear resistors caused by localized particle growth on the surface and by internal flaws. The presence of Sb2O3 was responsible for the observed particle growth. A part of the Bi2O3 on the surface was found to be in the reduced state. Two kinds of failure mode, cracking and puncturing, were observed when the samples were subjected to high-energy pulses. The puncture mode is caused by local melting around the regions of high current density, whereas the cracking mode results from thermally induced stresses.

Non-Newtonian boundary layers on a moving cylinder

The analysis of steady laminar forced convection boundary layer of power-law non-Newtonian fluids on a continuously moving cylinder with the surface maintained at a uniform temperature or uniform heat flux is presented. Of interest were the effects of power-law viscosity index, transverse curvature, generalized Prandtl number and streamwise coordinate on the local Nusselt number as well as on the velocity and temperature profiles. The two thermal boundary conditions yield quite similar results. Comparison of the calculated results with available series expansion solutions for a Newtonian fluid shows a very good performance of the present numerical procedure.

Hopping and non-hopping localized ionic motion: β-AgI and β-Ag3SI

In β-AgI and β-Ag3SI the ionic conductivity has been measured at frequencies from 1kHz to 2.6 GHz and from 10 MHz to 10 THz, respectively. In both phases we observe a conductivity increase of some orders of magnitude, due to localized types of motion of the silver ions. In β-AgI the increase is found at about 1 MHz and reflects cooperative back-and-forth hopping processes between adjacent tetrahedral sites. In β-Ag3SI the phenomenon occurs at microwave frequencies. Here it is caused by a non-hopping, non-periodic localized Ag+-motion within shallow potentials.

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.

Autoignition in a non-premixed medium: DNS studies on the effects of three-dimensional turbulence

Direct numerical simulation (DNS) results of autoignition in anon-premixed medium under an isotropic, homogeneous, and decaying turbulence are presented. The initial mixture consists of segregated fuel parcels randomly distributed within warm air, and the entire medium is subjected to a three-dimensional turbulence. Chemical kinetics is modeled by a four-step reduced reaction mechanism for autoignition of n-heptane/air mixture. Thus, this work overcomes the principal limitations of a previous contribution of the authors on two-dimensional DNS of autoignition with a one-step reaction model. Specific attention is focused on the differences in the effects of two- and three-dimensional turbulence on autoignition characteristics. The three-dimensional results show that ignition spots are most likely to originate at locations jointly corresponding to the most reactive mixture fraction and low scalar dissipation rate. Further, these ignition spots are found to originate at locations corresponding to the core of local vortical structures, and after ignition, the burning gases move toward the vortex periphery Such a movement is explained as caused by the cyclostrophic imbalance developed when the local gas density is variable. These results lead to the conclusion that the local ignition-zone structure does not conform to the classical stretched flamelet description. Parametric studies show that the ignition delay time decreases with an increase in turbulence intensity. Hence, these three-dimensional simulation results resolve the discrepancy between trends in experimental data and predictions from DNSs of two-dimensional turbulence. This qualitative difference between DNS results from three- and two-dimensional simulations is discussed and attributed to the effect of vortex stretching that is present in the former, but not in the latter.

Two stage iterative wiener filtering for speech enhancement

We formulate a two-stage Iterative Wiener filtering (IWF) approach to speech enhancement, bettering the performance of constrained IWF, reported in literature. The codebook constrained IWF (CCIWF) has been shown to be effective in achieving convergence of IWF in the presence of both stationary and non-stationary noise. To this, we include a second stage of unconstrained IWF and show that the speech enhancement performance can be improved in terms of average segmental SNR (SSNR), Itakura-Saito (IS) distance and Linear Prediction Coefficients (LPC) parameter coincidence. We also explore the tradeoff between the number of CCIWF iterations and the second stage IWF iterations.

Linear Filtering in MDCT domain

In this paper, expressions for convolution multiplication properties of MDCT are derived starting from the equivalent DFT representations. Using these expressions, methods for implementing linear filtering through block convolution in the MDCT domain are presented. The implementation is exact for symmetric filters and approximate for non-symmetric filters in the case of rectangular window based MDCT. For a general MDCT window function, the filtering is done on the windowed segments and hence the convolution is approximate for symmetric as well as non-symmetric filters. This approximation error is shown to be perceptually insignificant for symmetric impulse response filters. Moreover, the inherent $50 \%$ overlap between adjacent frames used in MDCT computation does reduce this approximation error similar to smoothing of other block processing errors. The presented techniques are useful for compressed domain processing of audio signals.

Fully Comprehensive Geometrically Non-Linear Dynamic Analysis of Multi-Body Beam Systems with Elastic Couplings

This paper is concerned with the dynamic analysis of flexible,non-linear multi-body beam systems. The focus is on problems where the strains within each elastic body (beam) remain small. Based on geometrically non-linear elasticity theory, the non-linear 3-D beam problem splits into either a linear or non-linear 2-D analysis of the beam cross-section and a non-linear 1-D analysis along the beam reference line. The splitting of the three-dimensional beam problem into two- and one-dimensional parts, called dimensional reduction,results in a tremendous savings of computational effort relative to the cost of three-dimensional finite element analysis,the only alternative for realistic beams. The analysis of beam-like structures made of laminated composite materials requires a much more complicated methodology. Hence, the analysis procedure based on Variational Asymptotic Method (VAM), a tool to carry out the dimensional reduction, is used here.The analysis methodology can be viewed as a 3-step procedure. First, the sectional properties of beams made of composite materials are determined either based on an asymptotic procedure that involves a 2-D finite element nonlinear analysis of the beam cross-section to capture trapeze effect or using strip-like beam analysis, starting from Classical Laminated Shell Theory (CLST). Second, the dynamic response of non-linear, flexible multi-body beam systems is simulated within the framework of energy-preserving and energy-decaying time integration schemes that provide unconditional stability for non-linear beam systems. Finally,local 3-D responses in the beams are recovered, based on the 1-D responses predicted in the second step. Numerical examples are presented and results from this analysis are compared with those available in the literature.

A DCT Based Filtering of Biomedical Images

Image filtering techniques have numerous potential applications in biomedical imaging and image processing. The design of filters largely depends on the a-priori knowledge about the type of noise corrupting the image and image features. This makes the standard filters to be application and image specific. The most popular filters such as average, Gaussian and Wiener reduce noisy artifacts by smoothing. However, this operation normally results in smoothing of the edges as well. On the other hand, sharpening filters enhance the high frequency details making the image non-smooth. An integrated general approach to design filters based on discrete cosine transform (DCT) is proposed in this study for optimal medical image filtering. This algorithm exploits the better energy compaction property of DCT and re-arrange these coefficients in a wavelet manner to get the better energy clustering at desired spatial locations. This algorithm performs optimal smoothing of the noisy image by preserving high and low frequency features. Evaluation results show that the proposed filter is robust under various noise distributions.

Cold model study of local liquid holdup in the dropping zone using x-ray

It is important to know and to quantify the liquid holdups both dynamic and static at local levels as it will lead to understand various blast furnace phenomena properly such as slag/metal.gas.solid reactions, gas flow behaviour and interfacial area between the gas/solid/liquid. In the present study, considering the importance of local liquid holdup and non-availability of holdup data in these systems, an attempt has been made to quantify the local holdups in the dropping and around raceway zones in a cold model study using a non-wetting packing for liquid. In order to quantify the liquid holdups at microscopic level, a previously developed technique, X-ray radiography, has been used. It is observed that the liquid flows in preferred paths or channels which carry droplets/rivulets. It has been found that local holdup in some regions of the packed bed is much higher than average at a particular flow rate and this can have important consequences for the correct modelling of such systems.