Bayesian parameter identification in dynamic state space models using modified measurement equations

When Markov chain Monte Carlo (MCMC) samplers are used in problems of system parameter identification, one would face computational difficulties in dealing with large amount of measurement data and (or) low levels of measurement noise. Such exigencies are likely to occur in problems of parameter identification in dynamical systems when amount of vibratory measurement data and number of parameters to be identified could be large. In such cases, the posterior probability density function of the system parameters tends to have regions of narrow supports and a finite length MCMC chain is unlikely to cover pertinent regions. The present study proposes strategies based on modification of measurement equations and subsequent corrections, to alleviate this difficulty. This involves artificial enhancement of measurement noise, assimilation of transformed packets of measurements, and a global iteration strategy to improve the choice of prior models. Illustrative examples cover laboratory studies on a time variant dynamical system and a bending-torsion coupled, geometrically non-linear building frame under earthquake support motions. (C) 2015 Elsevier Ltd. All rights reserved.

EXPERIENCE OF OPERATING OPEN TOP BIOMASS GASIFIER USING A LOW DENSITY FUEL-COCONUT FRONDS

Use of fuel other than woody generally has been limited to rice husk and other residues are rarely tried as a fuel in a gasification system. With the availability of woody biomass in most countries like India, alternates fuels are being explored for sustainable supply of fuel. Use of agro residues has been explored after briquetting. There are few feedstock's like coconut fronts, maize cobs, etc, that might require lesser preprocessing steps compared to briquetting. The paper presents a detailed investigation into using coconut fronds as a fuel in an open top down draft gasification system. The fuel has ash content of 7% and was dried to moisture levels of 12 %. The average bulk density was found to be 230 kg/m3 with a fuel size particle of an average size 40 mm as compared to 350 kg/m3 for a standard wood pieces. A typical dry coconut fronds weighs about 2.5kgs and on an average 6 m long and 90 % of the frond is the petiole which is generally used as a fuel. The focus was also to compare the overall process with respect to operating with a typical woody biomass like subabul whose ash content is 1 %. The open top gasification system consists of a reactor, cooling and cleaning system along with water treatment. The performance parameters studied were the gas composition, tar and particulates in the clean gas, water quality and reactor pressure drop apart from other standard data collection of fuel flow rate, etc. The average gas composition was found to be CO 15 1.0 % H-2 16 +/- 1% CH4 0.5 +/- 0.1 % CO2 12.0 +/- 1.0 % and rest N2 compared to CO 19 +/- 1.0 % H-2 17 +/- 1.0 %, CH4 1 +/- 0.2 %, CO2 12 +/- 1.0 % and rest N2. The tar and particulate content in the clean gas has been found to be about 10 and 12 mg/m3 in both cases. The presence of high ash content material increased the pressure drop with coconut frond compared to woody biomass.

A faster algorithm for testing polynomial representability of functions over finite integer rings

Given a function from Z(n) to itself one can determine its polynomial representability by using Kempner function. In this paper we present an alternative characterization of polynomial functions over Z(n) by constructing a generating set for the Z(n)-module of polynomial functions. This characterization results in an algorithm that is faster on average in deciding polynomial representability. We also extend the characterization to functions in several variables. (C) 2015 Elsevier B.V. All rights reserved.

``Conformational Simulation'' of Sulfamethizole by Molecular Complexation and Insights from Charge Density Analysis: Role of Intramolecular S center dot center dot center dot O Chalcogen Bonding

Intramolecular S center dot center dot center dot O chalcogen bonding and its potential to lock molecular conformation have been examined in the crystal forms of sulfamethizole, a sulfonamide antibiotic. Molecular complexes of sulfamethizole, including salts and cocrystal, have been synthesized, and their crystal structures were analyzed in order to examine the possible conformational preferences of the molecule in various ionic states and supramolecular environments (neutral/cocrystal, anionic salt, and cationic salt forms). The electrostatic potential mapped on Hirshfeld surfaces generated for these crystal forms provides insights into the possible binding modes of the drug in different environments. Further, the observed conformation locking feature has been rationalized in terms of the experimental charge density features of the intramolecular S center dot center dot O chalcogen bonding in sulfamethizole. The study quantitatively illustrates and rationalizes an intriguing case of a local minimum of molecular conformation being exclusively preferred over the global minimum, as it facilitates more efficient intermolecular interactions in a supramolecular environment.

Energy Sharing for Multiple Sensor Nodes With Finite Buffers

We consider the problem of finding optimal energy sharing policies that maximize the network performance of a system comprising of multiple sensor nodes and a single energy harvesting (EH) source. Sensor nodes periodically sense the random field and generate data, which is stored in the corresponding data queues. The EH source harnesses energy from ambient energy sources and the generated energy is stored in an energy buffer. Sensor nodes receive energy for data transmission from the EH source. The EH source has to efficiently share the stored energy among the nodes to minimize the long-run average delay in data transmission. We formulate the problem of energy sharing between the nodes in the framework of average cost infinite-horizon Markov decision processes (MDPs). We develop efficient energy sharing algorithms, namely Q-learning algorithm with exploration mechanisms based on the epsilon-greedy method as well as upper confidence bound (UCB). We extend these algorithms by incorporating state and action space aggregation to tackle state-action space explosion in the MDP. We also develop a cross entropy based method that incorporates policy parameterization to find near optimal energy sharing policies. Through simulations, we show that our algorithms yield energy sharing policies that outperform the heuristic greedy method.

A Reliable Residual Based A Posteriori Error Estimator for a Quadratic Finite Element Method for the Elliptic Obstacle Problem

A residual based a posteriori error estimator is derived for a quadratic finite element method (FEM) for the elliptic obstacle problem. The error estimator involves various residuals consisting of the data of the problem, discrete solution and a Lagrange multiplier related to the obstacle constraint. The choice of the discrete Lagrange multiplier yields an error estimator that is comparable with the error estimator in the case of linear FEM. Further, an a priori error estimate is derived to show that the discrete Lagrange multiplier converges at the same rate as that of the discrete solution of the obstacle problem. The numerical experiments of adaptive FEM show optimal order convergence. This demonstrates that the quadratic FEM for obstacle problem exhibits optimal performance.

Crystal landscape in the orcinol:4,4 `-bipyridine system: synthon modularity, polymorphism and transferability of multipole charge density parameters

Polymorphism in the orcinol: 4,4'-bipyridine cocrystal system is analyzed in terms of a robust convergent modular phenol...pyridine supramolecular synthon. Employing the Synthon Based Fragments Approach (SBFA) to transfer the multipole charge density parameters, it is demonstrated that the crystal landscape can be quantified in terms of intermolecular interaction energies in the five crystal forms so far isolated in this complex system. There are five crystal forms. The first has an open, divergent O-H...N based structure with alternating orcinol and bipyridine molecules. The other four polymorphs have different three-dimensional packing but all of them are similar at an interaction level, and are based on a modular O-H...N mediated supramolecular synthon that consists of two orcinol and two bipyridine molecules in a closed, convergent structure. The SBFA method, which depends on the modularity of synthons, provides good agreement between experiment and theory because it takes into account the supramolecular contribution to charge density. The existence of five crystal forms in this system shows that polymorphism in cocrystals need not be considered to be an unusual phenomenon. Studies of the crystal landscape could lead to an understanding of the kinetic pathways that control the crystallization processes, in other words the valleys in the landscape. These pathways are traditionally not considered in exercises pertaining to computational crystal structure prediction, which rather monitors the thermodynamics of the various stable forms in the system, in other words the peaks in the landscape.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Experimentally validated Raman Monte Carlo simulation for a cuboid object to obtain Raman spectroscopic signatures for hidden material

In conventional Raman spectroscopic measurements of liquids or surfaces the preferred geometry for detection of the Raman signal is the backscattering (or reflection) mode. For non-transparent layered materials, sub-surface Raman signals have been retrieved using spatially offset Raman spectroscopy (SORS), usually with light collection in the same plane as the point of excitation. However, as a result of multiple scattering in a turbid medium, Raman photons will be emitted in all directions. In this study, Monte Carlo simulations for a three-dimensional layered sample with finite geometry have been performed to confirm the detectability of Raman signals at all angles and at all sides of the object. We considered a non-transparent cuboid container (high density polyethylene) with explosive material (ammonium nitrate) inside. The simulation results were validated with experimental Raman intensities. Monte Carlo simulation results reveal that the ratio of sub-surface to surface signals improves at geometries other than backscattering. In addition, we demonstrate through simulations the effects of the absorption and scattering coefficients of the layers, and that of the diameter of the excitation beam. The advantage of collecting light from all possible 4 angles, over other collection modes, is that this technique is not geometry specific and molecular identification of layers underneath non-transparent surfaces can be obtained with minimal interference from the surface layer. To what extent all sides of the object will contribute to the total signal will depend on the absorption and scattering coefficients and the physical dimensions. Copyright (c) 2015 John Wiley & Sons, Ltd.

Effect of viscosity contrast on plume formation in density stratified fluids

We perform two and three dimensional numerical simulations of plume formation in density and viscosity stratified fluid systems. We show that the ambient to plume fluid viscosity ratio strongly affects the near wall plume structures (line or sheet plumes) such as plume spacing and shape of plumes. We observe that where mushroom-like plumes are observed for lower viscosity ratios, taller plumes with bulbous heads form for high viscosity ratios. Plume structure and spacing are in good agreement with experimental results. By studying the geometry of the line plumes and the flow in the circulation cells, we discuss the mechanisms of their formation and the dynamics of merging. We show that an increase in the viscosity ratio decreases the total length of line plumes in the planform which indicates a decreased mixing at higher viscosity ratios. (C) 2015 Elsevier Ltd. All rights reserved.

Solid-state optical absorption from optimally tuned time-dependent range-separated hybrid density functional theory

We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from optimally tuned range-separated hybrid density functional theory. The approach, which is fully couched within the formal framework of generalized Kohn-Sham theory, allows for the accurate prediction of exciton binding energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments and prior computations. We further show that with one adjustable parameter, set to produce the known band gap, this method accurately predicts band structures and optical spectra of silicon and lithium fluoride, prototypical covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to studies of materials of increasing size and complexity.

A finite variant of the Toom Model

We present results for a finite variant of the one-dimensional Toom model with closed boundaries. We show that the steady state distribution is not of product form, but is nonetheless simple. In particular, we give explicit formulas for the densities and some nearest neighbour correlation functions. We also give exact results for eigenvalues and multiplicities of the transition matrix using the theory of R-trivial monoids in joint work with A. Schilling, B. Steinberg and N. M. Thiery.

S center dot center dot center dot O chalcogen bonding in sulfa drugs: insights from multipole charge density and X-ray wavefunction of acetazolamide

Experimental charge density analysis combined with the quantum crystallographic technique of X-ray wavefunction refinement (XWR) provides quantitative insights into the intra-and intermolecular interactions formed by acetazolamide, a diuretic drug. Firstly, the analysis of charge density topology at the intermolecular level shows the presence of exceptionally strong interaction motifs such as a DDAA-AADD (D-donor, A-acceptor) type quadruple hydrogen bond motif and a sulfonamide dimer synthon. The nature and strength of intra-molecular S center dot center dot center dot O chalcogen bonding have been characterized using descriptors from the multipole model (MM) and XWR. Although pure geometrical criteria suggest the possibility of two intra-molecular S center dot center dot center dot O chalcogen bonded ring motifs, only one of them satisfies the ``orbital geometry'' so as to exhibit an interaction in terms of an electron density bond path and a bond critical point. The presence of `s-holes' on the sulfur atom leading to the S center dot center dot center dot O chalcogen bond has been visualized on the electrostatic potential surface and Laplacian isosurfaces close to the `reactive surface'. The electron localizability indicator (ELI) and Roby bond orders derived from the `experimental wave function' provide insights into the nature of S center dot center dot center dot O chalcogen bonding.

ABSOLUTE/CONVECTIVE INSTABILITY TRANSITION IN A BACKWARD FACING STEP COMBUSTOR: FUNDAMENTAL MECHANISM AND INFLUENCE OF DENSITY GRADIENT

Hydrodynamic instabilities of the flow field in lean premixed gas turbine combustors can generate velocity perturbations that wrinkle and distort the flame sheet over length scales that are smaller than the flame length. The resultant heat release oscillations can then potentially result in combustion instability. Thus, it is essential to understand the hydrodynamic instability characteristics of the combustor flow field in order to understand its overall influence on combustion instability characteristics. To this end, this paper elucidates the role of fluctuating vorticity production from a linear hydrodynamic stability analysis as the key mechanism promoting absolute/convective instability transitions in shear layers occurring in the flow behind a backward facing step. These results are obtained within the framework of an inviscid, incompressible, local temporal and spatio-temporal stability analysis. Vorticity fluctuations in this limit result from interaction between two competing mechanisms - (1) production from interaction between velocity perturbations and the base flow vorticity gradient and (2) baroclinic torque in the presence of base flow density gradients. This interaction has a significant effect on hydrodynamic instability characteristics when the base flow density and velocity gradients are co-located. Regions in the space of parameters characterizing the base flow velocity profile, i.e. shear layer thickness and ratio of forward to reverse flow velocity, corresponding to convective and absolute instability are identified. The implications of the present results on prior observations of flow instability in other flows such as heated jets and bluff-body stabilized flames is discussed.

Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.