A semi-analytical particle filter for identification of nonlinear oscillators

Particle filters find important applications in the problems of state and parameter estimations of dynamical systems of engineering interest. Since a typical filtering algorithm involves Monte Carlo simulations of the process equations, sample variance of the estimator is inversely proportional to the number of particles. The sample variance may be reduced if one uses a Rao-Blackwell marginalization of states and performs analytical computations as much as possible. In this work, we propose a semi-analytical particle filter, requiring no Rao-Blackwell marginalization, for state and parameter estimations of nonlinear dynamical systems with additively Gaussian process/observation noises. Through local linearizations of the nonlinear drift fields in the process/observation equations via explicit Ito-Taylor expansions, the given nonlinear system is transformed into an ensemble of locally linearized systems. Using the most recent observation, conditionally Gaussian posterior density functions of the linearized systems are analytically obtained through the Kalman filter. This information is further exploited within the particle filter algorithm for obtaining samples from the optimal posterior density of the states. The potential of the method in state/parameter estimations is demonstrated through numerical illustrations for a few nonlinear oscillators. The proposed filter is found to yield estimates with reduced sample variance and improved accuracy vis-a-vis results from a form of sequential importance sampling filter.

Nonlinear secondary whirl of an overhung rotor

Gravity critical speeds of rotors have hitherto been studied using linear analysis, and ascribed to rotor stiffness asymmetry. Here, we study an idealized asymmetric nonlinear overhung rotor model of Crandall and Brosens, spinning close to its gravity critical speed.Nonlinearities arise from finite displacements, and the rotor's staticlateral deflection under gravity is taken as small. Assuming small asymmetry and damping, slow modulations of whirl amplitudes are studied using the method of multiple scales. Inertia asymmetry appears only at second order. More interestingly, even without stiffness asymmetry, the gravity-induced resonance survives through geometric nonlinearities. The gravity resonant forcing does not influence the resonant mode at leading order, unlike the typical resonant oscillations. Nevertheless,the usual phenomena of resonances, namely saddle-node bifurcations, jump phenomena and hysteresis, are all observed. An unanticipated periodic solution branch is found. In the three-dimensional space oftwo modal coefficients and a detuning parameter, the full set of periodic solutions is found to be an imperfect version of three mutually intersecting curves: a straight line,a parabola and an ellipse.

Distance dependence of fluorescence resonance energy transfer

Deviations from the usual R (-6) dependence of the rate of fluorescence resonance energy transfer (FRET) on the distance between the donor and the acceptor have been a common scenario in the recent times. In this paper, we present a critical analysis of the distance dependence of FRET, and try to illustrate the non R (-6) type behaviour of the rate for the case of transfer from a localized electronic excitation on the donor, a dye molecule to three different energy acceptors with delocalized electronic excitations namely, graphene,two-dimensional semiconducting sheet and the case of such a semiconducting sheet rolled to obtain a nanotube. We use simple analytic models to understand the distance dependence in each case.

Universalization of access to modern energy services in Indian households-Economic and policy analysis

Provision of modern energy services for cooking (with gaseous fuels)and lighting (with electricity) is an essential component of any policy aiming to address health, education or welfare issues; yet it gets little attention from policy-makers. Secure, adequate, low-cost energy of quality and convenience is core to the delivery of these services. The present study analyses the energy consumption pattern of Indian domestic sector and examines the urban-rural divide and income energy linkage. A comprehensive analysis is done to estimate the cost for providing modern energy services to everyone by 2030. A public-private partnership-driven business model, with entrepreneurship at the core, is developed with institutional, financing and pricing mechanisms for diffusion of energy services. This approach, termed as EMPOWERS (entrepreneurship model for provision of wholesome energy-related basic services), if adopted, can facilitate large-scale dissemination of energy-efficient and renewable technologies like small-scale biogas/biofuel plants, and distributed power generation technologies to provide clean, safe, reliable and sustainable energy to rural households and urban poor. It is expected to integrate the processes of market transformation and entrepreneurship development involving government, NGOs, financial institutions and community groups as stakeholders. (C) 2009 Elsevier Ltd. All rights reserved.

R1SVM: A randomised nonlinear approach to large-scale anomaly detection

The problem of unsupervised anomaly detection arises in a wide variety of practical applications. While one-class support vector machines have demonstrated their effectiveness as an anomaly detection technique, their ability to model large datasets is limited due to their memory and time complexity for training. To address this issue for supervised learning of kernel machines, there has been growing interest in random projection methods as an alternative to the computationally expensive problems of kernel matrix construction and sup-port vector optimisation. In this paper we leverage the theory of nonlinear random projections and propose the Randomised One-class SVM (R1SVM), which is an efficient and scalable anomaly detection technique that can be trained on large-scale datasets. Our empirical analysis on several real-life and synthetic datasets shows that our randomised 1SVM algorithm achieves comparable or better accuracy to deep auto encoder and traditional kernelised approaches for anomaly detection, while being approximately 100 times faster in training and testing.

Suppression of the zero-order term in off-axis digital holography through nonlinear filtering

We present experimental validation of a new reconstruction method for off-axis digital holographic microscopy (DHM). This method effectively suppresses the object autocorrelation,namely, the zero-order term,from holographic data,thereby improving the reconstruction bandwidth of complex wavefronts. The algorithm is based on nonlinear filtering and can be applied to standard DHM setups with realistic recording conditions.We study the robustness of the technique under different experimental configurations,and quantitatively demonstrate its enhancement capabilities on phase signals.

Numerical Study of Heat Transfer From Pin-Fin Heat Sink Using Steady and Pulsated Impinging Jets

This paper investigates numerically the heat transfer characteristics of confined slot jet impingement on a pin-fin heat sink. A variety of pin-fin heat sinks is investigated, and the resulting enhancement of heat transfer studied. The distribution of heat transfer coefficient on the top surface of the base plate and that along the fin height are examined. Both steady and pulsated jets are studied. It is observed that for a steady jet impingement on a pin-fin heat sink, the effective heat transfer coefficient increases with fin height, leading to a corresponding decrease in base plate temperature for the same heat flux. In the case of pulsated jets, the influence of pulse frequency and the Reynolds number is examined, and their effect on the effective heat transfer coefficient is studied.

An Experimental Study on Energy Absorption Behavior of Polyurethane Foams

This article is concerned with a study on the energy absorption behavior of polyurethane (PU) foams such as flexible high resilience (HR), flexible viscoelastic (VE) and semi-rigid (SR) foams as a function of the overall foam density. Foam samples were prepared in the form of cubes by mixing appropriate polyol and isocyanate compounds produced by Huntsman International India Pvt. Ltd. in varying proportions leading to a range of densities for each type of foam. The cubical samples were tested under compressive load in a standard UTM. Based on the measured load-displacement behaviors, variations of peak load and energy-absorption attributes with respect to density are plotted for each type of foam and the possible existence of an optimum foam density is shown.

Analyzing changes in macro-level drivers of country-specific emissions : The role of consumption, technology and trade

Industrial ecology is an important field of sustainability science. It can be applied to study environmental problems in a policy relevant manner. Industrial ecology uses ecosystem analogy; it aims at closing the loop of materials and substances and at the same time reducing resource consumption and environmental emissions. Emissions from human activities are related to human interference in material cycles. Carbon (C), nitrogen (N) and phosphorus (P) are essential elements for all living organisms, but in excess have negative environmental impacts, such as climate change (CO2, CH4 N2O), acidification (NOx) and eutrophication (N, P). Several indirect macro-level drivers affect emissions change. Population and affluence (GDP/capita) often act as upward drivers for emissions. Technology, as emissions per service used, and consumption, as economic intensity of use, may act as drivers resulting in a reduction in emissions. In addition, the development of country-specific emissions is affected by international trade. The aim of this study was to analyse changes in emissions as affected by macro-level drivers in different European case studies. ImPACT decomposition analysis (IPAT identity) was applied as a method in papers I III. The macro-level perspective was applied to evaluate CO2 emission reduction targets (paper II) and the sharing of greenhouse gas emission reduction targets (paper IV) in the European Union (EU27) up to the year 2020. Data for the study were mainly gathered from official statistics. In all cases, the results were discussed from an environmental policy perspective. The development of nitrogen oxide (NOx) emissions was analysed in the Finnish energy sector during a long time period, 1950 2003 (paper I). Finnish emissions of NOx began to decrease in the 1980s as the progress in technology in terms of NOx/energy curbed the impact of the growth in affluence and population. Carbon dioxide (CO2) emissions related to energy use during 1993 2004 (paper II) were analysed by country and region within the European Union. Considering energy-based CO2 emissions in the European Union, dematerialization and decarbonisation did occur, but not sufficiently to offset population growth and the rapidly increasing affluence during 1993 2004. The development of nitrogen and phosphorus load from aquaculture in relation to salmonid consumption in Finland during 1980 2007 was examined, including international trade in the analysis (paper III). A regional environmental issue, eutrophication of the Baltic Sea, and a marginal, yet locally important source of nutrients was used as a case. Nutrient emissions from Finnish aquaculture decreased from the 1990s onwards: although population, affluence and salmonid consumption steadily increased, aquaculture technology improved and the relative share of imported salmonids increased. According to the sustainability challenge in industrial ecology, the environmental impact of the growing population size and affluence should be compensated by improvements in technology (emissions/service used) and with dematerialisation. In the studied cases, the emission intensity of energy production could be lowered for NOx by cleaning the exhaust gases. Reorganization of the structure of energy production as well as technological innovations will be essential in lowering the emissions of both CO2 and NOx. Regarding the intensity of energy use, making the combustion of fuels more efficient and reducing energy use are essential. In reducing nutrient emissions from Finnish aquaculture to the Baltic Sea (paper III) through technology, limits of biological and physical properties of cultured fish, among others, will eventually be faced. Regarding consumption, salmonids are preferred to many other protein sources. Regarding trade, increasing the proportion of imports will outsource the impacts. Besides improving technology and dematerialization, other viewpoints may also be needed. Reducing the total amount of nutrients cycling in energy systems and eventually contributing to NOx emissions needs to be emphasized. Considering aquaculture emissions, nutrient cycles can be partly closed through using local fish as feed replacing imported feed. In particular, the reduction of CO2 emissions in the future is a very challenging task when considering the necessary rates of dematerialisation and decarbonisation (paper II). Climate change mitigation may have to focus on other greenhouse gases than CO2 and on the potential role of biomass as a carbon sink, among others. The global population is growing and scaling up the environmental impact. Population issues and growing affluence must be considered when discussing emission reductions. Climate policy has only very recently had an influence on emissions, and strong actions are now called for climate change mitigation. Environmental policies in general must cover all the regions related to production and impacts in order to avoid outsourcing of emissions and leakage effects. The macro-level drivers affecting changes in emissions can be identified with the ImPACT framework. Statistics for generally known macro-indicators are currently relatively well available for different countries, and the method is transparent. In the papers included in this study, a similar method was successfully applied in different types of case studies. Using transparent macro-level figures and a simple top-down approach are also appropriate in evaluating and setting international emission reduction targets, as demonstrated in papers II and IV. The projected rates of population and affluence growth are especially worth consideration in setting targets. However, sensitivities in calculations must be carefully acknowledged. In the basic form of the ImPACT model, the economic intensity of consumption and emission intensity of use are included. In seeking to examine consumption but also international trade in more detail, imports were included in paper III. This example demonstrates well how outsourcing of production influences domestic emissions. Country-specific production-based emissions have often been used in similar decomposition analyses. Nevertheless, trade-related issues must not be ignored.

The tired foot. The effects of physical activity, energy expenditure and sleep on the diabetic foot population

Background The most common pathway to development of diabetes foot ulcers is repetitive daily activity stress on the plantar surface of the neuropathic foot. Studies suggest an association between different diabetic foot complications and physical activity. However, to the best of the authors knowledge the steps/day and sleep patterns of people with diabetic foot ulcers has yet to be investigated. This observational study aims to investigate the physical activity and sleep patterns of three groups of adults with type 2 diabetes and different foot complications Methods Participants with type 2 diabetes were recruited into three groups: 1. those with no reported foot complications (DNIL), 2. those with diagnosis of neuropathy (DPN) and 3. those with a neuropathic ulcer (DFU). Exclusion criteria included peripheral arterial disease and mobility aid use. Participants wore a SenseWear Pro 3 Armband continuously for 7 days and completed an Epworth Sleepiness Scale. The Armband is a validated automated measure of activity (walking steps, average Metabolic Equivalent Task (MET), physical activity (>3 METs) duration), energy expenditure(kJ) (total and physical activity (>3 METs)) and sleep (duration). Data on age, sex, BMI, diabetes duration and HbA1c were also collected. Results Sixty-Six (14 DNIL, 22 DPN and 30 DFU's participants were recruited; 71% males, mean age 61(±12) years, diabetes duration 13(±9) years, HbA1c 8.3(±2.8), BMI 32.6(±5.9), average METs 1.2(0.2). Significant differences were reported in mean(SD) steps/day (5,859(±2,381) in DNIL; 5,007(±3,349) in DPN and 3,271(±2,417) in DFU's and daily energy expenditure (10,868(±1,307)kJ in DNIL; 11,060(±1,916)kJ in DPN and 13,006(± 3,559) in DFU's(p <0.05). No significant differences were reported for average METs, physical activity duration or energy expenditure, sleep time or Epworth score (p>0.1). Conclusions Preliminary findings suggest people with diabetes are sedentary. Results indicate that patients with a diabetic foot ulcer work significantly less than those with neuropathy or nil complications and use significantly more energy to do so. Sleep Parameters showed no differences. Recruitment is still on going.

Highly efficient spin-conversion effect leading to energy up-converted electroluminescence in singlet fission photovoltaics

Free charge generation in donor-acceptor (D-A) based organic photovoltaic diodes (OPV) progresses through formation of charge-transfer (CT) and charge-separated (CS) states and excitation decay to the triplet level is considered as a terminal loss. On the other hand a direct excitation decay to the triplet state is beneficial for multiexciton harvesting in singlet fission photovoltaics (SF-PV) and the formation of CT-state is considered as a limiting factor for multiple triplet harvesting. These two extremes when present in a D-A system are expected to provide important insights into the mechanism of free charge generation and spin-character of bimolecular recombination in OPVs. Herein, we present the complete cycle of events linked to spin conversion in the model OPV system of rubrene/C60. By tracking the spectral evolution of photocurrent generation at short-circuit and close to open-circuit conditions we are able to capture spectral changes to photocurrent that reveal the triplet character of CT-state. Furthermore, we unveil an energy up-conversion effect that sets in as a consequence of triplet population build-up where triplet-triplet annihilation (TTA) process effectively regenerates the singlet excitation. This detailed balance is shown to enable a rare event of photon emission just above the open-circuit voltage (VOC) in OPVs.

A Search for Energy Minimized Sequences of Proteins

In this paper, we present numerical evidence that supports the notion of minimization in the sequence space of proteins for a target conformation. We use the conformations of the real proteins in the Protein Data Bank (PDB) and present computationally efficient methods to identify the sequences with minimum energy. We use edge-weighted connectivity graph for ranking the residue sites with reduced amino acid alphabet and then use continuous optimization to obtain the energy-minimizing sequences. Our methods enable the computation of a lower bound as well as a tight upper bound for the energy of a given conformation. We validate our results by using three different inter-residue energy matrices for five proteins from protein data bank (PDB), and by comparing our energy-minimizing sequences with 80 million diverse sequences that are generated based on different considerations in each case. When we submitted some of our chosen energy-minimizing sequences to Basic Local Alignment Search Tool (BLAST), we obtained some sequences from non-redundant protein sequence database that are similar to ours with an E-value of the order of 10(-7). In summary, we conclude that proteins show a trend towards minimizing energy in the sequence space but do not seem to adopt the global energy-minimizing sequence. The reason for this could be either that the existing energy matrices are not able to accurately represent the inter-residue interactions in the context of the protein environment or that Nature does not push the optimization in the sequence space, once it is able to perform the function.

Dynamics of loop formation in a semiflexible polymer

The dynamics of loop formation by linear polymer chains has been a topic of several theoretical and experimental studies. Formation of loops and their opening are key processes in many important biological processes. Loop formation in flexible chains has been extensively studied by many groups. However, in the more realistic case of semiflexible polymers, not much results are available. In a recent study [K. P. Santo and K. L. Sebastian, Phys. Rev. E 73, 031923 (2006)], we investigated opening dynamics of semiflexible loops in the short chain limit and presented results for opening rates as a function of the length of the chain. We presented an approximate model for a semiflexible polymer in the rod limit based on a semiclassical expansion of the bending energy of the chain. The model provided an easy way to describe the dynamics. In this paper, using this model, we investigate the reverse process, i.e., the loop formation dynamics of a semiflexible polymer chain by describing the process as a diffusion-controlled reaction. We make use of the ``closure approximation'' of Wilemski and Fixman [G. Wilemski and M. Fixman, J. Chem. Phys. 60, 878 (1974)], in which a sink function is used to represent the reaction. We perform a detailed multidimensional analysis of the problem and calculate closing times for a semiflexible chain. We show that for short chains, the loop formation time tau decreases with the contour length of the polymer. But for longer chains, it increases with length obeying a power law and so it has a minimum at an intermediate length. In terms of dimensionless variables, the closing time is found to be given by tau similar to L-n exp(const/L), where n=4.5-6. The minimum loop formation time occurs at a length L-m of about 2.2-2.4. These are, indeed, the results that are physically expected, but a multidimensional analysis leading to these results does not seem to exist in the literature so far.

Third post-Newtonian angular momentum flux and the secular evolution of orbital elements for inspiralling compact binaries in quasi-elliptical orbits

The angular-momentum flux from an inspiralling binary system of compact objects moving in quasi-elliptical orbits is computed at the third post-Newtonian (3PN) order using the multipolar post-Minkowskian wave generation formalism. The 3PN angular-momentum flux involves the instantaneous, tail, and tail-of-tails contributions as for the 3PN energy flux, and in addition a contribution due to nonlinear memory. We average the angular-momentum flux over the binary's orbit using the 3PN quasi-Keplerian representation of elliptical orbits. The averaged angular-momentum flux provides the final input needed for gravitational-wave phasing of binaries moving in quasi-elliptical orbits. We obtain the evolution of orbital elements under 3PN gravitational radiation reaction in the quasi-elliptic case. For small eccentricities, we give simpler limiting expressions relevant for phasing up to order e(2). This work is important for the construction of templates for quasi-eccentric binaries, and for the comparison of post-Newtonian results with the numerical relativity simulations of the plunge and merger of eccentric binaries.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.