Faster computation of the Karhunen-Loeve expansion using its domain independence property

The goal of this work is to reduce the cost of computing the coefficients in the Karhunen-Loeve (KL) expansion. The KL expansion serves as a useful and efficient tool for discretizing second-order stochastic processes with known covariance function. Its applications in engineering mechanics include discretizing random field models for elastic moduli, fluid properties, and structural response. The main computational cost of finding the coefficients of this expansion arises from numerically solving an integral eigenvalue problem with the covariance function as the integration kernel. Mathematically this is a homogeneous Fredholm equation of second type. One widely used method for solving this integral eigenvalue problem is to use finite element (FE) bases for discretizing the eigenfunctions, followed by a Galerkin projection. This method is computationally expensive. In the current work it is first shown that the shape of the physical domain in a random field does not affect the realizations of the field estimated using KL expansion, although the individual KL terms are affected. Based on this domain independence property, a numerical integration based scheme accompanied by a modification of the domain, is proposed. In addition to presenting mathematical arguments to establish the domain independence, numerical studies are also conducted to demonstrate and test the proposed method. Numerically it is demonstrated that compared to the Galerkin method the computational speed gain in the proposed method is of three to four orders of magnitude for a two dimensional example, and of one to two orders of magnitude for a three dimensional example, while retaining the same level of accuracy. It is also shown that for separable covariance kernels a further cost reduction of three to four orders of magnitude can be achieved. Both normal and lognormal fields are considered in the numerical studies. (c) 2014 Elsevier B.V. All rights reserved.

Molecules support morphology: species status of South Indian populations of the widely distributed Hanuman langur

The taxonomy of the Hanuman langur (Semnopithecus spp.), a widely distributed Asian colobine monkey, has been in a flux for a long time due to much disagreement between various classification schemes. However, results from a recent field-based morphological study were consistent with Hill's (Ceylon J Sci 21:277-305, 1939) species level classification scheme. Here we tested the validity of S. hypoleucos and S. priam, the two South Indian species recognized by Hill. To this end, one mitochondrial and four nuclear markers were sequenced from over 72 non-invasive samples of Hanuman langurs and S. johnii collected from across India. The molecular data were subjected to various tree building methods. The nuclear data was also used in a Bayesian structure analysis and to determine the genealogical sorting index of each hypothesized species. Results from nuclear data suggest that the South Indian population of Hanuman langur consists of two units that correspond to the species recognized by Hill. However in the mitochondrial tree S. johnii and S. priam were polyphyletic probably due to retention of ancestral polymorphism and/or low levels of hybridization. Implications of these results on conservation of Hanuman langurs are also discussed.

Interpreting the upper level structure of the Madden-Julian oscillation

The nonlinear response of a spherical shallow water model to an imposed heat source in the presence of realistic zonal mean zonal winds is investigated numerically. The solutions exhibit elongated, meridionally tilted ridges and troughs indicative of a poleward dispersion of wave activity. As the speed of the jets is increased, the equatorial Kelvin wave is unaffected but the global Rossby wave train coalesces to form a compact, amplified quadrupole structure that bears a striking resemblance to the observed upper level structure of the Madden-Julian oscillation. In the presence of strong subtropical westerly jets, the advection of planetary vorticity by the meridional flow and relative vorticity by the zonally averaged background flow conspire to create the distinctive quadrupole configuration of flanking Rossby waves.

Boundary Tracking using a Hybrid Cyclic Pursuit Scheme

There has been a lot of work in the literature, related to the mapping of boundaries of regions, using multiple agents. Most of these are based on optimization techniques or rely on potential fields to drive the agents towards the boundary and then retain them there while they space out evenly along the perimeter or surface (in two-dimensional and three-dimensional cases, respectively). In this paper an algorithm to track the boundary of a region in space is provided based on the cyclic pursuit scheme. This enables the agents to constantly move along the perimeter in a cluster, thereby tracking a dynamically changing boundary. The trajectories of the agents provide a sketch of the boundary. The use of multiple agents may facilitate minimization of tracking error by providing accurate estimates of points on the boundary, besides providing redundancy. Simulation results are provided to highlight the performance of the proposed scheme.

ADAPTATION OF MW LEVEL NATURAL GAS LEAN BURN ENGINE FOR PRODUCER GAS OPERATION A GRID CONNECTED POWER PLANT

Availability of producer gas engines at MW being limited necessitates to adapt engine from natural gas operation. The present work focus on the development of necessary kit for adapting a 12 cylinder lean burn turbo-charged natural gas engine rated at 900 kWe (Waukesha make VHP5904LTD) to operate on producer and set up an appropriate capacity biomass gasification system for grid linked power generation in Thailand. The overall plant configuration had fuel processing, drying, reactor, cooling and cleaning system, water treatment, engine generator and power evacuation. The overall project is designed for evacuation of 1.5 MWe power to the state grid and had 2 gasification system with the above configuration and 3 engines. Two gasification system each designed for about 1100 kg/hr of woody biomass was connected to the engine using a producer gas carburetor for the necessary Air to fuel ratio control. In the use of PG to fuel IC engines, it has been recognized that the engine response will differ as compared to the response with conventional fueled operation due to the differences in the thermo-physical properties of PG. On fuelling a conventional engine with PG, power de-rating can be expected due to the lower calorific value (LCV), lower adiabatic flame temperature (AFT) and the lower than unity product to reactant more ratio. Further the A/F ratio for producer gas is about 1/10th that of natural gas and requires a different carburetor for engine operation. The research involved in developing a carburetor for varying load conditions. The patented carburetor is based on area ratio control, consisting of a zero pressure regulator and a separate gas and air line along with a mixing zone. The 95 litre engine at 1000 rpm has an electrical efficiency of 33.5 % with a heat input of 2.62 MW. Each engine had two carburetors designed for producer gas flow each capable of handling about 1200 m3/hr in order to provide similar engine heat input at a lower conversion efficiency. Cold flow studies simulating the engine carburetion system results showed that the A/F was maintained in the range of 1.3 +/- 0.1 over the entire flow range. Initially, the gasification system was tested using woody biomass and the gas composition was found to be CO 15 +/- 1.5 % H-2 22 +/- 2% CH4 2.2 +/- 0.5 CO2 11.25 +/- 1.4 % and rest N-2, with the calorific value in the range of 5.0 MJ/kg. After initial trials on the engine to fine tune the control system and adjust various engine operating parameter a peak load of 800 kWe was achieved, while a stable operating conditions was found to be at 750 kWe which is nearly 85 % of the natural gas rating. The specific fuel consumption was found to be 0.9 kg of biomass per kWh.

An Implementation of Partitioned Scheduling Scheme for Hard Real-Time Tasks in Multicore Linux with Fair Share for Linux Tasks

The correctness of a hard real-time system depends its ability to meet all its deadlines. Existing real-time systems use either a pure real-time scheduler or a real-time scheduler embedded as a real-time scheduling class in the scheduler of an operating system (OS). Existing implementations of schedulers in multicore systems that support real-time and non-real-time tasks, permit the execution of non-real-time tasks in all the cores with priorities lower than those of real-time tasks, but interrupts and softirqs associated with these non-real-time tasks can execute in any core with priorities higher than those of real-time tasks. As a result, the execution overhead of real-time tasks is quite large in these systems, which, in turn, affects their runtime. In order that the hard real-time tasks can be executed in such systems with minimal interference from other Linux tasks, we propose, in this paper, an integrated scheduler architecture, called SchedISA, which aims to considerably reduce the execution overhead of real-time tasks in these systems. In order to test the efficacy of the proposed scheduler, we implemented partitioned earliest deadline first (P-EDF) scheduling algorithm in SchedISA on Linux kernel, version 3.8, and conducted experiments on Intel core i7 processor with eight logical cores. We compared the execution overhead of real-time tasks in the above implementation of SchedISA with that in SCHED_DEADLINE's P-EDF implementation, which concurrently executes real-time and non-real-time tasks in Linux OS in all the cores. The experimental results show that the execution overhead of real-time tasks in the above implementation of SchedISA is considerably less than that in SCHED_DEADLINE. We believe that, with further refinement of SchedISA, the execution overhead of real-time tasks in SchedISA can be reduced to a predictable maximum, making it suitable for scheduling hard real-time tasks without affecting the CPU share of Linux tasks.

Seventeen-Level Inverter Formed by Cascading Flying Capacitor and Floating Capacitor H-Bridges

A multilevel inverter for generating 17 voltage levels using a three-level flying capacitor inverter and cascaded H-bridge modules with floating capacitors has been proposed. Various aspects of the proposed inverter like capacitor voltage balancing have been presented in the present paper. Experimental results are presented to study the performance of the proposed converter. The stability of the capacitor balancing algorithm has been verified both during transients and steady-state operation. All the capacitors in this circuit can be balanced instantaneously by using one of the pole voltage combinations. Another advantage of this topology is its ability to generate all the voltages from a single dc-link power supply which enables back-to-back operation of converter. Also, the proposed inverter can be operated at all load power factors and modulation indices. Additional advantage is, if one of the H-bridges fail, the inverter can still be operated at full load with reduced number of levels. This configuration has very low dv/dt and common-mode voltage variation.

Parametric Studies on Saltwater Intrusion into Coastal Aquifers for Anticipate Sea Level Rise

Saltwater intrusion into coastal aquifers is a global issue, exacerbated by increasing demands for freshwater in coastal regions. This study investigates into the parametric analysis on saltwater intrusion in a conceptual, coastal, unconfined aquifer considering wide range of freshwater draft and anticipated sea level rise. The saltwater intrusion under various circumstances is simulated through parametric studies using MODFLOW, MT3DMS and SEAWAT. The MODFLOW is used to simulate the groundwater flow system under changing hydro-dynamics in coastal aquifer. To simulate solute transport MT3DMS and SEAWAT is used. The saltwater intrusion process has direct bearing on hydraulic conductivity and inversely related to porosity. It may also be noted that increase in recharge rate considered in the study does not have much influence on saltwater intrusion. Effect of freshwater draft at locations beyond half of the width of the aquifer considered has marginal effect and hence can be considered as safe zone for freshwater withdrawals. Due to the climate change effect, the anticipated rise in sea level of 0.88 m over a century is considered in the investigation. This causes increase in salinity intrusion by about 25%. The combined effect of sea level rise and freshwater draft (C) 2015 The Authors. Published by Elsevier B.V.

Reduction in DNA topoisomerase I level affects growth, phenotype and nucleoid architecture of Mycobacterium smegmatis

The steady-state negative supercoiling of eubacterial genomes is maintained by the action of DNA topoisomerases. Topoisomerase distribution varies in different species of mycobacteria. While Mycobacterium tuberculosis (Mtb) contains a single type I (Topol) and a single type II (Gyrase) enzyme, Mycobacterium smegmatis (Msm) and other members harbour additional relaxases. Topol is essential for Mtb survival. However, the necessity of Topol or other relaxases in Msm has not been investigated. To recognize the importance of Topol for growth, physiology and gene expression of Msm, we have developed a conditional knock-down strain of Topol in Msm. The Topol-depleted strain exhibited extremely slow growth and drastic changes in phenotypic characteristics. The cessation of growth indicates the essential requirement of the enzyme for the organism in spite of having additional DNA relaxation enzymes in the cell. Notably, the imbalance in Topol level led to the altered expression of topology modulatory proteins, resulting in a diffused nucleoid architecture. Proteomic and transcript analysis of the mutant indicated reduced expression of the genes involved in central metabolic pathways and core DNA transaction processes. RNA polymerase (RNAP) distribution on the transcription units was affected in the Topol-depleted cells, suggesting global alteration in transcription. The study thus highlights the essential requirement of Topol in the maintenance of cellular phenotype, growth characteristics and gene expression in mycobacteria. A decrease in Topol level led to altered RNAP occupancy and impaired transcription elongation, causing severe downstream effects.

Inner Bound on the GDOF of the K-User MIMO Gaussian Symmetric Interference Channel

The K-user multiple input multiple output (MIMO) Gaussian symmetric interference channel where each transmitter has M antennas and each receiver has N antennas is studied from a generalized degrees of freedom (GDOF) perspective. An inner bound on the GDOF is derived using a combination of techniques such as treating interference as noise, zero forcing (ZF) at the receivers, interference alignment (IA), and extending the Han-Kobayashi (HK) scheme to K users, as a function of the number of antennas and the log INR/log SNR level. Several interesting conclusions are drawn from the derived bounds. It is shown that when K > N/M + 1, a combination of the HK and IA schemes performs the best among the schemes considered. When N/M < K <= N/M + 1, the HK-scheme outperforms other schemes and is found to be GDOF optimal in many cases. In addition, when the SNR and INR are at the same level, ZF-receiving and the HK-scheme have the same GDOF performance.

Performance of New Dynamic Benefit-Weighted Scheduling Scheme in DiffServ Networks

In this paper, we design a new dynamic packet scheduling scheme suitable for differentiated service (DiffServ) network. Designed dynamic benefit weighted scheduling (DBWS) uses a dynamic weighted computation scheme loosely based on weighted round robin (WRR) policy. It predicts the weight required by expedited forwarding (EF) service for the current time slot (t) based on two criteria; (i) previous weight allocated to it at time (t-1), and (ii) the average increase in the queue length of EF buffer. This prediction provides smooth bandwidth allocation to all the services by avoiding overbooking of resources for EF service and still providing guaranteed services for it. The performance is analyzed for various scenarios at high, medium and low traffic conditions. The results show that packet loss is minimized, end to end delay is minimized and jitter is reduced and therefore meet quality of service (QoS) requirement of a network.

Macro-Level Assessment of Seismically Induced Landslide Hazard for the State of Sikkim, India Based On GIS Technique

This paper presents a macro-level seismic landslide hazard assessment for the entire state of Sikkim, India, based on the Newmark's methodology. The slope map of Sikkim was derived from ASTER Global Digital Elevation Model (GDEM). Seismic shaking in terms of peak horizontal acceleration (PHA) at bedrock level was estimated from deterministic seismic hazard analysis (DSHA), considering point source model. Peak horizontal acceleration at the surface level for the study area was estimated based on nonlinear site amplification technique, considering B-type NEHRP site class. The PHA at surface was considered to induce driving forces on slopes, thus causing landslides. Knowing the surface level PHA and slope angle, the seismic landslide hazard assessment for each grid point was carried out using Newmark's analysis. The critical static factor of safety required to resist landslide for the PHA (obtained from deterministic analysis) was evaluated and its spatial variation throughout the study area is presented. For any slope in the study area, if the in-situ (available) static factor of safety is greater than the static factor of safety required to resist landslide as predicted in the present study, that slope is considered to be safe.

Coupling geomorphology parameters with lithology for micro-level groundwater resource assessment: a case study from semi-arid hard rock terrain in Tamil Nadu, India

The standard procedure of groundwater resource estimation in India till date is based on the specific yield parameters of each rock type (lithology) derived through pumping test analysis. Using the change in groundwater level, specific yield, and area of influence, groundwater storage change could be estimated. However, terrain conditions in the form of geomorphological variations have an important bearing on the net groundwater recharge. In this study, an attempt was made to use both lithology and geomorphology as input variables to estimate the recharge from different sources in each lithology unit influenced by the geomorphic conditions (lith-geom), season wise separately. The study provided a methodological approach for an evaluation of groundwater in a semi-arid hard rock terrain in Tirunelveli, Tamil Nadu, India. While characterizing the gneissic rock, it was found that the geomorphologic variations in the gneissic rock due to weathering and deposition behaved differently with respect to aquifer recharge. The three different geomorphic units identified in gneissic rock (pediplain shallow weathered (PPS), pediplain moderate weathered (PPM), and buried pediplain moderate (BPM)) showed a significant variation in recharge conditions among themselves. It was found from the study that Peninsular gneiss gives a net recharge value of 0.13 m/year/unit area when considered as a single unit w.r.t. lithology, whereas the same area considered with lith-geom classes gives recharge values between 0.1 and 0.41 m/year presenting a different assessment. It is also found from this study that the stage of development (SOD) for each lith-geom unit in Peninsular gneiss varies from 168 to 230 %, whereas the SOD is 223 % for the lithology as a single unit.

Electronic structure origin of conductivity and oxygen reduction activity changes in low-level Cr-substituted (La, Sr)MnO3

The electronic structure of the (La0.8Sr0.2)(0.98)Mn1-xCrxO3 model series (x = 0, 0.05, or 0.1) was measured using soft X-ray synchrotron radiation at room and elevated temperature. O K-edge near-edge X-ray absorption fine structure (NEXAFS) spectra showed that low-level chromium substitution of (La, Sr)MnO3 resulted in lowered hybridisation between O 2p orbitals and M 3d and M 4sp valance orbitals. Mn L-3-edge resonant photoemission spectroscopy measurements indicated lowered Mn 3d-O 2p hybridisation with chromium substitution. Deconvolution of O K-edge NEXAFS spectra took into account the effects of exchange and crystal field splitting and included a novel approach whereby the pre-peak region was described using the nominally filled t(2g) up arrow state. 10% chromium substitution resulted in a 0.17 eV lowering in the energy of the t(2g) up arrow state, which appears to provide an explanation for the 0.15 eV rise in activation energy for the oxygen reduction reaction, while decreased overlap between hybrid O 2p-Mn 3d states was in qualitative agreement with lowered electronic conductivity. An orbital-level understanding of the thermodynamically predicted solid oxide fuel cell cathode poisoning mechanism involving low-level chromium substitution on the B-site of (La, Sr)MnO3 is presented. (C) 2015 AIP Publishing LLC.

Implicit scheme for meshless compressible Euler solver

In this paper, an implicit scheme is presented for a meshless compressible Euler solver based on the Least Square Kinetic Upwind Method (LSKUM). The Jameson and Yoon's split flux Jacobians formulation is very popular in finite volume methodology, which leads to a scalar diagonal dominant matrix for an efficient implicit procedure (Jameson & Yoon, 1987). However, this approach leads to a block diagonal matrix when applied to the LSKUM meshless method. The above split flux Jacobian formulation, along with a matrix-free approach, has been adopted to obtain a diagonally dominant, robust and cheap implicit time integration scheme. The efficacy of the scheme is demonstrated by computing 2D flow past a NACA 0012 airfoil under subsonic, transonic and supersonic flow conditions. The results obtained are compared with available experiments and other reliable computational fluid dynamics (CFD) results. The present implicit formulation shows good convergence acceleration over the RK4 explicit procedure. Further, the accuracy and robustness of the scheme in 3D is demonstrated by computing the flow past an ONERA M6 wing and a clipped delta wing with aileron deflection. The computed results show good agreement with wind tunnel experiments and other CFD computations.