Novel Relay Selection Rules for Average Interference-Constrained Cognitive AF Relay Networks

Cooperative relaying combined with selection exploits spatial diversity to significantly improve the performance of interference-constrained secondary users in an underlay cognitive radio (CR) network. However, unlike conventional relaying, the state of the links between the relay and the primary receiver affects the choice of the relay. Further, while the optimal amplify-and-forward (AF) relay selection rule for underlay CR is well understood for the peak interference-constraint, this is not so for the less conservative average interference constraint. For the latter, we present three novel AF relay selection (RS) rules, namely, symbol error probability (SEP)-optimal, inverse-of-affine (IOA), and linear rules. We analyze the SEPs of the IOA and linear rules and also develop a novel, accurate approximation technique for analyzing the performance of AF relays. Extensive numerical results show that all the three rules outperform several RS rules proposed in the literature and generalize the conventional AF RS rule.

Three-dimensional numerical analysis of geocell-reinforced soft clay beds by considering the actual geometry of geocell pockets

Due to its complex honeycomb structure, the numerical modeling of the geocell has always been a big challenge. Generally, the equivalent composite approach is used to model the geocells. In the equivalent composite approach, the geocellsoil composite is treated as the soil layer with improved strength and stiffness values. Though this approach is very simple, it is unrealistic to model the geocells as the soil layer. This paper presents a more realistic approach of modeling the geocells in three-dimensional (3D) framework by considering the actual curvature of the geocell pocket. A square footing resting on geocell reinforced soft clay bed was modeled using the ``fast Lagrangian analysis of continua in 3D'' (FLAC(3D)) finite difference package. Three different material models, namely modified Cam-clay, Mohr-Coulomb, and linear elastic were used to simulate the behaviour of foundation soil, infill soil and the geocell, respectively. It was found that the geocells distribute the load laterally to the wider area below the footing as compared to the unreinforced case. More than 50% reduction in the stress was observed in the clay bed in the presence of geocells. In addition to geocells, two other cases, namely, only geogrid and geocell with additional basal geogrid cases were also simulated. The numerical model was systematically validated with the results of the physical model tests. Using the validated numerical model, parametric studies were conducted to evaluate the influence of various geocell properties on the performance of reinforced clay beds.

Numerical study of laminar, standing hydraulic jumps in a planar geometry

We solve the two-dimensional, planar Navier-Stokes equations to simulate a laminar, standing hydraulic jump using a Volume-of-Fluid method. The geometry downstream of the jump has been designed to be similar to experimental conditions by including a pit at the edge of the platform over which liquid film flows. We obtain jumps with and without separation. Increasing the inlet Froude number pushes the jump downstream and makes the slope of the jump weaker, consistent with experimental observations of circular jumps, and decreasing the Reynolds number brings the jump upstream while making it steeper. We study the effect of the length of the domain and that of a downstream obstacle on the structure and location of the jump. The transient flow which leads to a final steady jump is described for the first time to our knowledge. In the moderate Reynolds number regime, we obtain steady undular jumps with a separated bubble underneath the first few undulations. Interestingly, surface tension leads to shortening of wavelength of these undulations. We show that the undulations can be explained using the inviscid theory of Benjamin and Lighthill (Proc. R. Soc. London, Ser. A, 1954). We hope this new finding will motivate experimental verification.

A hybrid method for stochastic response analysis of a vibrating structure

Response analysis of a linear structure with uncertainties in both structural parameters and external excitation is considered here. When such an analysis is carried out using the spectral stochastic finite element method (SSFEM), often the computational cost tends to be prohibitive due to the rapid growth of the number of spectral bases with the number of random variables and the order of expansion. For instance, if the excitation contains a random frequency, or if it is a general random process, then a good approximation of these excitations using polynomial chaos expansion (PCE) involves a large number of terms, which leads to very high cost. To address this issue of high computational cost, a hybrid method is proposed in this work. In this method, first the random eigenvalue problem is solved using the weak formulation of SSFEM, which involves solving a system of deterministic nonlinear algebraic equations to estimate the PCE coefficients of the random eigenvalues and eigenvectors. Then the response is estimated using a Monte Carlo (MC) simulation, where the modal bases are sampled from the PCE of the random eigenvectors estimated in the previous step, followed by a numerical time integration. It is observed through numerical studies that this proposed method successfully reduces the computational burden compared with either a pure SSFEM of a pure MC simulation and more accurate than a perturbation method. The computational gain improves as the problem size in terms of degrees of freedom grows. It also improves as the timespan of interest reduces.

A FRAMEWORK FOR THE ERROR ANALYSIS OF DISCONTINUOUS FINITE ELEMENT METHODS FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS AND APPLICATIONS TO C-0 IP METHODS

In this article, an abstract framework for the error analysis of discontinuous Galerkin methods for control constrained optimal control problems is developed. The analysis establishes the best approximation result from a priori analysis point of view and delivers a reliable and efficient a posteriori error estimator. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. Subsequently, the applications of C-0 interior penalty methods for a boundary control problem as well as a distributed control problem governed by the biharmonic equation subject to simply supported boundary conditions are discussed through the abstract analysis. Numerical experiments illustrate the theoretical findings.

A modified approach to numerical solution of Fredholm integral equations of the second kind

A modified approach to obtain approximate numerical solutions of Fredholin integral equations of the second kind is presented. The error bound is explained by the aid of several illustrative examples. In each example, the approximate solution is compared with the exact solution, wherever possible, and an excellent agreement is observed. In addition, the error bound in each example is compared with the one obtained by the Nystrom method. It is found that the error bound of the present method is smaller than the ones obtained by the Nystrom method. Further, the present method is successfully applied to derive the solution of an integral equation arising in a special Dirichlet problem. (C) 2015 Elsevier Inc. All rights reserved.

QoS Constrained Optimal Sink and Relay Placement in Planned Wireless Sensor Networks

We are given a set of sensors at given locations, a set of potential locations for placing base stations (BSs, or sinks), and another set of potential locations for placing wireless relay nodes. There is a cost for placing a BS and a cost for placing a relay. The problem we consider is to select a set of BS locations, a set of relay locations, and an association of sensor nodes with the selected BS locations, so that the number of hops in the path from each sensor to its BS is bounded by h(max), and among all such feasible networks, the cost of the selected network is the minimum. The hop count bound suffices to ensure a certain probability of the data being delivered to the BS within a given maximum delay under a light traffic model. We observe that the problem is NP-Hard, and is hard to even approximate within a constant factor. For this problem, we propose a polynomial time approximation algorithm (SmartSelect) based on a relay placement algorithm proposed in our earlier work, along with a modification of the greedy algorithm for weighted set cover. We have analyzed the worst case approximation guarantee for this algorithm. We have also proposed a polynomial time heuristic to improve upon the solution provided by SmartSelect. Our numerical results demonstrate that the algorithms provide good quality solutions using very little computation time in various randomly generated network scenarios.

Weakly corrected numerical solutions to stochastically driven nonlinear dynamical systems

A method to weakly correct the solutions of stochastically driven nonlinear dynamical systems, herein numerically approximated through the Eule-Maruyama (EM) time-marching map, is proposed. An essential feature of the method is a change of measures that aims at rendering the EM-approximated solution measurable with respect to the filtration generated by an appropriately defined error process. Using Ito's formula and adopting a Monte Carlo (MC) setup, it is shown that the correction term may be additively applied to the realizations of the numerically integrated trajectories. Numerical evidence, presently gathered via applications of the proposed method to a few nonlinear mechanical oscillators and a semi-discrete form of a 1-D Burger's equation, lends credence to the remarkably improved numerical accuracy of the corrected solutions even with relatively large time step sizes. (C) 2015 Elsevier Inc. All rights reserved.

Numerical Study of Counterflow Diffusion Flame and Water Spray Interaction

This paper reports numerical investigation concerning the interaction of a laminar methane-air counterflow diffusion flame with monodisperse and polydisperse water spray. Commercial code ANSYS FLUENT with reduced chemistry has been used for investigation. Effects of strain rate, Sauter mean diameter (SMD), and droplet size distribution on the temperature along stagnation streamline have been studied. Flame extinction using polydisperse water spray has also been explored. Comparison of monodisperse and polydisperse droplet distribution on flame properties reveals suitability of polydisperse spray in flame temperature reduction beyond a particular SMD. This study also provides a numerical framework to study flame-spray interaction and extinction.

Numerical analysis of segmented-electrode Orbitraps

Simple geometries which are possible alternatives for the Orbitrap are studied in this paper. We have taken up for numerical investigation two segmented-electrode structures, ORB1 and ORB2, to mimic the electric field of the Orbitrap. In the ORB1, the inner spindle-like electrode and the outer barrel-like electrode of the Orbitrap have been replaced by 35 rings and 35 discs of fixed radii, respectively. In this structure two segmented end cap electrodes have been added. In this geometry, different potentials are applied to the different electrodes keeping top-bottom symmetry intact. In the second geometry, ORB2, the inner and outer electrodes of the Orbitrap were replaced by an approximate step structure which follows the profile of the Orbitrap electrodes. In the present study 45 steps have been used. In the ORB2, like the Orbitrap, the inner electrode is held at a negative potential and the outer electrode is at ground potential. For the purpose of comparing the performance of ORB1 and ORB2 with that of the Orbitrap, the following studies have been undertaken: (1) variation of electric potential, (2) computation of ion trajectories, (3) simulation of image currents. These studies have been carried out using both 2D and 3D Boundary Element Method (BEM), the 3D BEM was developed specifically for this study. It has been seen in these investigations that ORB1 and ORB2 have performance similar to that of the Orbitrap, with the performance of the ORB1 being seen to be marginally superior to that of the ORB2. It has been shown that with proper optimization, geometries containing far fewer electrodes can be used as mass analyzers. A novel technique of optimization of the electric field has been proposed with the objective of minimizing the dependence of axial frequency of ion motion on the initial position of an ion. The results on the optimization of 9 and 15 segmented-electrode traps having the same design as ORB1 show that it can provide accurate mass analysis. (C) 2015 Elsevier B.V. All rights reserved.

Scaling analysis and numerical studies on water vapour adsorption in a columnar porous silica gel bed

This article presents a theoretical analysis of heat and mass transfer in a silica gel + water adsorption process using scaling principles. A two-dimensional columnar packed adsorber domain is chosen for the study, with side and bottom walls cooled and vapour inlet from the top. The adsorption process is initiated from the cold walls with a temperature jump of 15 K, whereas the water vapour supply is maintained at a constant inlet pressure of 1 kPa. The first part of the study is dedicated to deriving relevant scales for the adsorption process by an order of magnitude analysis of energy, continuity and momentum equations. In the latter part, the derived scales are compared with the outcome of numerical studies performed for various domain widths and aspect ratio of bed. A good correlation between scaling and simulation results is observed, thereby validating the scaling approach. (C) 2015 Elsevier Ltd. All rights reserved.

Numerical analysis of segmented-electrode Orbitraps

Simple geometries which are possible alternatives for the Orbitrap are studied in this paper. We have taken up for numerical investigation two segmented-electrode structures, ORB1 and ORB2, to mimic the electric field of the Orbitrap. In the ORB1, the inner spindle-like electrode and the outer barrel-like electrode of the Orbitrap have been replaced by 35 rings and 35 discs of fixed radii, respectively. In this structure two segmented end cap electrodes have been added. In this geometry, different potentials are applied to the different electrodes keeping top-bottom symmetry intact. In the second geometry, ORB2, the inner and outer electrodes of the Orbitrap were replaced by an approximate step structure which follows the profile of the Orbitrap electrodes. In the present study 45 steps have been used. In the ORB2, like the Orbitrap, the inner electrode is held at a negative potential and the outer electrode is at ground potential. For the purpose of comparing the performance of ORB1 and ORB2 with that of the Orbitrap, the following studies have been undertaken: (1) variation of electric potential, (2) computation of ion trajectories, (3) simulation of image currents. These studies have been carried out using both 2D and 3D Boundary Element Method (BEM), the 3D BEM was developed specifically for this study. It has been seen in these investigations that ORB1 and ORB2 have performance similar to that of the Orbitrap, with the performance of the ORB1 being seen to be marginally superior to that of the ORB2. It has been shown that with proper optimization, geometries containing far fewer electrodes can be used as mass analyzers. A novel technique of optimization of the electric field has been proposed with the objective of minimizing the dependence of axial frequency of ion motion on the initial position of an ion. The results on the optimization of 9 and 15 segmented-electrode traps having the same design as ORB1 show that it can provide accurate mass analysis. (C) 2015 Elsevier B.V. All rights reserved.

Numerical Investigations for Leakage and Windage Heating in Straight-Through Labyrinth Seals

The ability to quantify leakage flow and windage heating for labyrinth seals with honeycomb lands is critical in understanding gas turbine engine system performance and predicting its component life. Variety of labyrinth seal configurations (number of teeth, stepped or straight, honeycomb cell size) are in use in gas turbines, and for each configuration, there are many geometric factors that can impact a seal's leakage and windage characteristics. This paper describes the development of a numerical methodology aimed at studying the effect of honeycomb lands on leakage and windage heating. Specifically, a three-dimensional computational fluid dynamics (CFD) model is developed utilizing commercial finite volume-based software incorporating the renormalization group (RNG) k-epsilon turbulence model with modified Schmidt number. The modified turbulence model is benchmarked and fine-tuned based on several experiments. Using this model, a broad parametric study is conducted by varying honeycomb cell size, pressure ratio (PR), and radial clearance for a four-tooth straight-through labyrinth seal. The results show good agreement with available experimental data. They further indicate that larger honeycomb cells predict higher seal leakage and windage heating at tighter clearances compared to smaller honeycomb cells and smooth lands. However, at open seal clearances larger honeycomb cells have lower leakage compared to smaller honeycomb cells.

Numerical prediction of nozzle flow separation: Issue of turbulence modeling

Numerical simulation of separated flows in rocket nozzles is challenging because existing turbulence models are unable to predict it correctly. This paper addresses this issue with the Spalart-Allmaras and Shear Stress Transport (SST) eddy-viscosity models, which predict flow separation with moderate success. Their performances have been compared against experimental data for a conical and two contoured subscale nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the nozzle pressure ratio (NPR) and nozzle type. A careful assessment indicated how the model had to be tuned for better, consistent prediction. It is learnt that SST model's failure is caused by limiting of the shear stress inside boundary layer according to Bradshaw's assumption, and by over prediction of jet spreading rate. Accordingly, SST's coefficients were empirically modified to match the experimental wall pressure data. Results confirm that accurate RANS prediction of separation depends on the correct capture of the jet spreading rate, and that it is feasible over a wide range of NPRs by modified values of the diffusion coefficients in the turbulence model. (C) 2015 Elsevier Masson SAS. All rights reserved.

A Stochastic Approximation Algorithm for Quantile Estimation

In this paper, we present two new stochastic approximation algorithms for the problem of quantile estimation. The algorithms uses the characterization of the quantile provided in terms of an optimization problem in 1]. The algorithms take the shape of a stochastic gradient descent which minimizes the optimization problem. Asymptotic convergence of the algorithms to the true quantile is proven using the ODE method. The theoretical results are also supplemented through empirical evidence. The algorithms are shown to provide significant improvement in terms of memory requirement and accuracy.