Blind equalization and sequence estimation using Markov chain Monte Carlo methods

This paper discusses the problem of restoring a digital input signal that has been degraded by an unknown FIR filter in noise, using the Gibbs sampler. A method for drawing a random sample of a sequence of bits is presented; this is shown to have faster convergence than a scheme by Chen and Li, which draws bits independently. ©1998 IEEE.

Simulation of solitary wave impact on coastal structures using weakly compressible and incompressible SPH methods

In this paper, we engage a Lagrangian, particle-based CFD method, named Smoothed Particle Hydrodynamic (SPH) to study the solitary wave motion and its impact on coastal structures. Two-dimensional weakly compressible and incompressible SPH models were applied to simulate wave impacting on seawall and schematic coastal house. The results confirmed the accuracy of both models for predicting the wave surface profiles. The incompressible SPH model performed better in predicting the pressure field and impact loadings on coastal structures than the weakly compressible SPH model. The results are in qualitatively agreement with experimental results. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).

Energy scattering in weakly non-linear systems.

The Chinese Tam-Tam exhibits non-linear behavior in its vibro-acoustic response. The frequency content of the response during free, unforced vibration smoothly changes, with energy being progressively smeared out over a greater bandwidth with time. This is used as a motivating case for the general study of the phenomenon of energy cascading through weak nonlinearity. Numerical models based upon the Fermi-Pasta-Ulam system of non-linearly coupled oscillators, modified with the addition of damping, have been developed. These were used to study the response of ensembles of systems with randomized natural frequencies. Results from simulations will be presented here. For un-damped systems, individual ensemble members exhibit cyclical energy exchange between linear modes, but the ensemble average displays a steady state. For the ensemble response of damped systems, lightly damped modes can exhibit an effective damping which is higher than predicated by linear theory. The presence of a non-linearity provides a path for energy flow to other modes, increasing the apparent damping spectrum at some frequencies and reducing it at others. The target of this work is a model revealing the governing parameters of a generic system of this type and leading to predictions of the ensemble response.

BarraCUDA – a Fast Sequence Mapping Software using Graphics Processing Units

High-throughput DNA sequencing (HTS) instruments today are capable of generating millions of sequencing reads in a short period of time, and this represents a serious challenge to current bioinformatics pipeline in processing such an enormous amount of data in a fast and economical fashion. Modern graphics cards are powerful processing units that consist of hundreds of scalar processors in parallel in order to handle the rendering of high-definition graphics in real-time. It is this computational capability that we propose to harness in order to accelerate some of the time-consuming steps in analyzing data generated by the HTS instruments. We have developed BarraCUDA, a novel sequence mapping software that utilizes the parallelism of NVIDIA CUDA graphics cards to map sequencing reads to a particular location on a reference genome. While delivering a similar mapping fidelity as other mainstream programs , BarraCUDA is a magnitude faster in mapping throughput compared to its CPU counterparts. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the mapping throughput. BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the mapping of millions of sequencing reads generated by HTS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology. BarraCUDA is currently available at http://seqbarracuda.sf.net

Image sequence segmentation using mixture-models with uniqueness and spatio-temporal consistency constraints

The use of mixture-model techniques for motion estimation and image sequence segmentation was discussed. The issues such as modeling of occlusion and uncovering, determining the relative depth of the objects in a scene, and estimating the number of objects in a scene were also investigated. The segmentation algorithm was found to be computationally demanding, but the computational requirements were reduced as the motion parameters and segmentation of the frame were initialized. The method provided a stable description, in whichthe addition and removal of objects from the description corresponded to the entry and exit of objects from the scene.

BarraCUDA - a fast short read sequence aligner using graphics processing units.

BACKGROUND: With the maturation of next-generation DNA sequencing (NGS) technologies, the throughput of DNA sequencing reads has soared to over 600 gigabases from a single instrument run. General purpose computing on graphics processing units (GPGPU), extracts the computing power from hundreds of parallel stream processors within graphics processing cores and provides a cost-effective and energy efficient alternative to traditional high-performance computing (HPC) clusters. In this article, we describe the implementation of BarraCUDA, a GPGPU sequence alignment software that is based on BWA, to accelerate the alignment of sequencing reads generated by these instruments to a reference DNA sequence. FINDINGS: Using the NVIDIA Compute Unified Device Architecture (CUDA) software development environment, we ported the most computational-intensive alignment component of BWA to GPU to take advantage of the massive parallelism. As a result, BarraCUDA offers a magnitude of performance boost in alignment throughput when compared to a CPU core while delivering the same level of alignment fidelity. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the alignment throughput. CONCLUSIONS: BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the alignment of millions of sequencing reads generated by NGS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology.BarraCUDA is currently available from http://seqbarracuda.sf.net.

Modeling of SFR cores with Serpent-DYN3D codes sequence

DYN3D reactor dynamics nodal diffusion code was originally developed for the analysis of Light Water Reactors. In this paper, we demonstrate the feasibility of using DYN3D for modeling of fast spectrum reactors. A homogenized cross sections data library was generated using continuous energy Monte-Carlo code Serpent which provides significant modeling flexibility compared with traditional deterministic lattice transport codes and tolerable execution time. A representative sodium cooled fast reactor core was modeled with the Serpent-DYN3D code sequence and the results were compared with those produced by ERANOS code and with a 3D full core Monte-Carlo solution. Very good agreement between the codes was observed for the core integral parameters and power distribution suggesting that the DYN3D code with cross section library generated using Serpent can be reliably used for the analysis of fast reactors. © 2012 Elsevier Ltd. All rights reserved.

Numerical study of the impact of process variations on the motional resistance of weakly coupled MEMS resonators

This paper presents a numerical study of the impact of process-induced variations on the achievable motional resistance Rx of one-dimensional, cyclic and cross-coupled architectures of electrostatically transduced MEMS resonators operating in the 250 kHz range. Monte Carlo numerical simulations which accounted for up to 0.75% variation in critical resonator feature sizes were initiated on 1, 2, 3, 4, 5 and 9 coupled MEMS resonators for three distinct coupling architectures. Improvements of 100X in the spread of Rx and 2.7X in mean achievable Rx are reported for the case of 9 resonators when implemented in the cross-coupled topology, as opposed to the traditional one-dimensional chain. © 2013 IEEE.