Solution of a class of surface water wave scattering problems involving non-reflecting vertical barriers

A large class of scattering problems of surface water waves by vertical barriers lead to mixed boundary value problems for Laplace equation. Specific attentions are paid, in the present article, to highlight an analytical method to handle this class of problems of surface water wave scattering, when the barriers in question are non-reflecting in nature. A new set of boundary conditions is proposed for such non-reflecting barriers and tile resulting boundary value problems are handled in the linearized theory of water waves. Three basic poblems of scattering by vertical barriers are solved. The present new theory of non-reflecting vertical barriers predict new transmission coefficients and tile solutions of tile mathematical problems turn out to be extremely simple and straight forward as compared to the solution for other types of barriers handled previously.

On the solution of bivariate population balance equations for aggregation: Pivotwise expansion of solution domain

The solution of a bivariate population balance equation (PBE) for aggregation of particles necessitates a large 2-d domain to be covered. A correspondingly large number of discretized equations for particle populations on pivots (representative sizes for bins) are solved, although at the end only a relatively small number of pivots are found to participate in the evolution process. In the present work, we initiate solution of the governing PBE on a small set of pivots that can represent the initial size distribution. New pivots are added to expand the computational domain in directions in which the evolving size distribution advances. A self-sufficient set of rules is developed to automate the addition of pivots, taken from an underlying X-grid formed by intersection of the lines of constant composition and constant particle mass. In order to test the robustness of the rule-set, simulations carried out with pivotwise expansion of X-grid are compared with those obtained using sufficiently large fixed X-grids for a number of composition independent and composition dependent aggregation kernels and initial conditions. The two techniques lead to identical predictions, with the former requiring only a fraction of the computational effort. The rule-set automatically reduces aggregation of particles of same composition to a 1-d problem. A midway change in the direction of expansion of domain, effected by the addition of particles of different mean composition, is captured correctly by the rule-set. The evolving shape of a computational domain carries with it the signature of the aggregation process, which can be insightful in complex and time dependent aggregation conditions. (c) 2012 Elsevier Ltd. All rights reserved.

Quasi-based hierarchical clustering for land cover mapping using satellite images

This paper presents an improved hierarchical clustering algorithm for land cover mapping problem using quasi-random distribution. Initially, Niche Particle Swarm Optimization (NPSO) with pseudo/quasi-random distribution is used for splitting the data into number of cluster centers by satisfying Bayesian Information Criteria (BIC). Themain objective is to search and locate the best possible number of cluster and its centers. NPSO which highly depends on the initial distribution of particles in search space is not been exploited to its full potential. In this study, we have compared more uniformly distributed quasi-random with pseudo-random distribution with NPSO for splitting data set. Here to generate quasi-random distribution, Faure method has been used. Performance of previously proposed methods namely K-means, Mean Shift Clustering (MSC) and NPSO with pseudo-random is compared with the proposed approach - NPSO with quasi distribution(Faure). These algorithms are used on synthetic data set and multi-spectral satellite image (Landsat 7 thematic mapper). From the result obtained we conclude that use of quasi-random sequence with NPSO for hierarchical clustering algorithm results in a more accurate data classification.

Automatic compilation of MATLAB programs for synergistic execution on heterogeneous processors

MATLAB is an array language, initially popular for rapid prototyping, but is now being increasingly used to develop production code for numerical and scientific applications. Typical MATLAB programs have abundant data parallelism. These programs also have control flow dominated scalar regions that have an impact on the program's execution time. Today's computer systems have tremendous computing power in the form of traditional CPU cores and throughput oriented accelerators such as graphics processing units(GPUs). Thus, an approach that maps the control flow dominated regions to the CPU and the data parallel regions to the GPU can significantly improve program performance. In this paper, we present the design and implementation of MEGHA, a compiler that automatically compiles MATLAB programs to enable synergistic execution on heterogeneous processors. Our solution is fully automated and does not require programmer input for identifying data parallel regions. We propose a set of compiler optimizations tailored for MATLAB. Our compiler identifies data parallel regions of the program and composes them into kernels. The problem of combining statements into kernels is formulated as a constrained graph clustering problem. Heuristics are presented to map identified kernels to either the CPU or GPU so that kernel execution on the CPU and the GPU happens synergistically and the amount of data transfer needed is minimized. In order to ensure required data movement for dependencies across basic blocks, we propose a data flow analysis and edge splitting strategy. Thus our compiler automatically handles composition of kernels, mapping of kernels to CPU and GPU, scheduling and insertion of required data transfer. The proposed compiler was implemented and experimental evaluation using a set of MATLAB benchmarks shows that our approach achieves a geometric mean speedup of 19.8X for data parallel benchmarks over native execution of MATLAB.

Bipolar Poisson Solution for Independent Double-Gate MOSFET

We propose a new set of input voltage equations (IVEs) for independent double-gate MOSFET by solving the governing bipolar Poisson equation (PE) rigorously. The proposed IVEs, which involve the Legendre's incomplete elliptic integral of the first kind and Jacobian elliptic functions and are valid from accumulation to inversion regimes, are shown to have good agreement with the numerical solution of the same PE for all bias conditions.

An extended finite element model coupled with level set method for stable pitting corrosion

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution), which requires only three inputs, namely the solid metal concentration, saturation concentration of the dissolved metal ions and diffusion coefficient. A combined eXtended Finite Element Model (XFEM) and level set method is developed in this paper. The extended finite element model handles the jump discontinuity in the metal concentrations at the interface, by using discontinuous-derivative enrichment formulation for concentration discontinuity at the interface. This eliminates the requirement of using front conforming mesh and re-meshing after each time step as in conventional finite element method. A numerical technique known as level set method tracks the position of the moving interface and updates it over time. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed method is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions.

An extended finite-element model coupled with level set method for analysis of growth of corrosion pits in metallic structures

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution). The interface experiences a jump discontinuity in metal concentration. The extended finite-element model (XFEM) handles this jump discontinuity by using discontinuous-derivative enrichment formulation, eliminating the requirement of using front conforming mesh and re-meshing after each time step as in the conventional finite-element method. However, prior interface location is required so as to solve the governing equations for concentration field for which a numerical technique, the level set method, is used for tracking the interface explicitly and updating it over time. The level set method is chosen as it is independent of shape and location of the interface. Thus, a combined XFEM and level set method is developed in this paper. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed model is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions. An empirical model for pitting potential is also derived based on the finite-element results. Studies show that pitting profile depends on factors such as ion concentration, solution pH and temperature to a large extent. Studying the individual and combined effects of these factors on pitting potential is worth knowing, as pitting potential directly influences corrosion rate.

QEP Based Solution for Power Allocation in Amplify and Forward Networks

The problem of cooperative beamforming for maximizing the achievable data rate of an energy constrained two-hop amplify-and-forward (AF) network is considered. Assuming perfect channel state information (CSI) of all the nodes, we evaluate the optimal scaling factor for the relay nodes. Along with individual power constraint on each of the relay nodes, we consider a weighted sum power constraint. The proposed iterative algorithm initially solves a set of relaxed problems with weighted sum power constraint and then updates the solution to accommodate individual constraints. These relaxed problems in turn are solved using a sequence of Quadratic Eigenvalue Problems (QEP). The key contribution of this letter is the generalization of cooperative beamforming to incorporate both the individual and weighted sum constraint. Furthermore, we have proposed a novel algorithm based on Quadratic Eigenvalue Problem (QEP) and discussed its convergence.

A new solution for the roommate problem: The Q-stable matchings

The aim of this paper is to propose a new solution for the roommate problem with strict preferences. We introduce the solution of maximum irreversibility and consider almost stable matchings (Abraham et al. [2])and maximum stable matchings (Ta [30] [32]). We find that almost stable matchings are incompatible with the other two solutions. Hence, to solve the roommate problem we propose matchings that lie at the intersection of the maximum irreversible matchings and maximum stable matchings, which are called Q-stable matchings. These matchings are core consistent and we offer an effi cient algorithm for computing one of them. The outcome of the algorithm belongs to an absorbing set.

Intramolecular conflict : conformation and self-assembly of architecturally complex macromolecules in solution

The solution behavior of linear polymer chains is well understood, having been the subject of intense study throughout the previous century. As plastics have become ubiquitous in everyday life, polymer science has grown into a major field of study. The conformation of a polymer in solution depends on the molecular architecture and its interactions with the surroundings. Developments in synthetic techniques have led to the creation of precision-tailored polymeric materials with varied topologies and functionalities. In order to design materials with the desired properties, it is imperative to understand the relationships between polymer architecture and their conformation and behavior. To meet that need, this thesis investigates the conformation and self-assembly of three architecturally complex macromolecular systems with rich and varied behaviors driven by the resolution of intramolecular conflicts. First we describe the development of a robust and facile synthetic approach to reproducible bottlebrush polymers (Chapter 2). The method was used to produce homologous series of bottlebrush polymers with polynorbornene backbones, which revealed the effect of side-chain and backbone length on the overall conformation in both good and theta solvent conditions (Chapter 3). The side-chain conformation was obtained from a series of SANS experiments and determined to be indistinguishable from the behavior of free linear polymer chains. Using deuterium-labeled bottlebrushes, we were able for the first time to directly observe the backbone conformation of a bottlebrush polymer which showed self-avoiding walk behavior. Secondly, a series of SANS experiments was conducted on a homologous series of Side Group Liquid Crystalline Polymers (SGLCPs) in a perdeuterated small molecule liquid crystal (5CB). Monodomain, aligned, dilute samples of SGLCP-b-PS block copolymers were seen to self-assemble into complex micellar structures with mutually orthogonally oriented anisotropies at different length scales (Chapter 4). Finally, we present the results from the first scattering experiments on a set of fuel-soluble, associating telechelic polymers. We observed the formation of supramolecular aggregates in dilute (≤0.5wt%) solutions of telechelic polymers and determined that the choice of solvent has a significant effect on the strength of association and the size of the supramolecules (Chapter 5). A method was developed for the direct estimation of supramolecular aggregation number from SANS data. The insight into structure-property relationships obtained from this work will enable the more targeted development of these molecular architectures for their respective applications.

Active reservoir management: a model solution

Steady-state procedures, of their very nature, cannot deal with dynamic situations. Statistical models require extensive calibration, and predictions often have to be made for environmental conditions which are often outside the original calibration conditions. In addition, the calibration requirement makes them difficult to transfer to other lakes. To date, no computer programs have been developed which will successfully predict changes in species of algae. The obvious solution to these limitations is to apply our limnological knowledge to the problem and develop functional models, so reducing the requirement for such rigorous calibration. Reynolds has proposed a model, based on fundamental principles of algal response to environmental events, which has successfully recreated the maximum observed biomass, the timing of events and a fair simulation of the species succession in several lakes. A forerunner of this model was developed jointly with Welsh Water under contract to Messrs. Wallace Evans and Partners, for use in the Cardiff Bay Barrage study. In this paper the authors test a much developed form of this original model against a more complex data-set and, using a simple example, show how it can be applied as an aid in the choice of management strategy for the reduction of problems caused by eutrophication. Some further developments of the model are indicated.

Technical and mapping support for marine spatial planning: Final report

Washington depends on a healthy coastal and marine ecosystem to maintain a thriving economy and vibrant communities. These ecosystems support critical habitats for wildlife and a growing number of often competing ocean activities, such as fishing, transportation, aquaculture, recreation, and energy production. Planners, policy makers and resource managers are being challenged to sustainably balance ocean uses, and environmental conservation in a finite space and with limited information. This balancing act can be supported by spatial planning. Marine spatial planning (MSP) is a planning process that enables integrated, forward looking, and consistent decision making on the human uses of the oceans and coasts. It can improve marine resource management by planning for human uses in locations that reduce conflict, increase certainty, and support a balance among social, economic, and ecological benefits we receive from ocean resources. In March 2010, the Washington state legislature enacted a marine spatial planning law (RCW §43.372) to address resource use conflicts in Washington waters. In 2011, a report to the legislature and a workshop on human use data provided guidance for the marine spatial planning process. The report outlines a set of recommendations for the State to effectively undertake marine spatial planning and this work plan will support some of these recommendations, such as: federal integration, regional coordination, developing mechanisms to integrate scientific and technical expertise, developing data standards, and accessing and sharing spatial data. In 2012 the Governor amended the existing law to focus funding on mapping and ecosystem assessments for Washington’s Pacific coast and the legislature provided $2.1 million in funds to begin marine spatial planning off Washington’s coast. The funds are appropriated through the Washington Department of Natural Resources Marine Resources Stewardship Account with coordination among the State Ocean Caucus, the four Coastal Treaty Tribes, four coastal Marine Resource Committees and the newly formed stakeholder body, the Washington Coastal Marine Advisory Council.

Integration of fisheries acoustics surveys and bathymetric mapping to characterize midwater-seafloor habitats of US Virgin Islands and Puerto Rico (2008-2010)

NOAA’s Coral Reef Conservation program (CRCP) develops coral reef management priorities by bringing together various partners to better understand threats to coral reef ecosystems with the goal of conserving, protecting and restoring these resources. Place-based and ecosystem-based management approaches employed by CRCP require that spatially explicit information about benthic habitats and fish utilization are available to characterize coral reef ecosystems and set conservation priorities. To accomplish this, seafloor habitat mapping of coral reefs around the U.S. Virgin Islands (USVI) and Puerto Rico has been ongoing since 2004. In 2008, fishery acoustics surveys were added to NOAA survey missions in the USVI and Puerto Rico to assess fish distribution and abundance in relation to benthic habitats in high priority conservation areas. NOAA’s National Centers for Coastal Ocean Science (NCCOS) have developed fisheries acoustics survey capabilities onboard the NOAA ship Nancy Foster to complement the CRCP seafloor habitat mapping effort spearheaded by the Center for Coastal Monitoring and Assessment Biogeography Branch (CCMA-BB). The integration of these activities has evolved on the Nancy Foster over the three years summarized in this report. A strategy for improved operations and products has emerged over that time. Not only has the concurrent operation of multibeam and fisheries acoustics surveys been beneficial in terms of optimizing ship time and resources, this joint effort has advanced an integrated approach to characterizing bottom and mid-water habitats and the fishes associated with them. CCMA conducts multibeam surveys to systematically map and characterize coral reef ecosystems, resulting in products such as high resolution bathymetric maps, backscatter information, and benthic habitat classification maps. These products focus on benthic features and live bottom habitats associated with them. NCCOS Centers (the Center for Coastal Fisheries and Habitat Research and the Center for Coastal Environmental Health and Biomolecular Research) characterize coral reef ecosystems by using fisheries acoustics methods to capture biological information through the entire water column. Spatially-explicit information on marine resources derived from fisheries acoustics surveys, such as maps of fish density, supports marine spatial planning strategies and decision making by providing a biological metric for evaluating coral reef ecosystems and assessing impacts from pollution, fishing pressure, and climate change. Data from fisheries acoustics surveys address management needs by providing a measure of biomass in management areas, detecting spatial and temporal responses in distribution relative to natural and anthropogenic impacts, and identifying hotspots that support high fish abundance or fish aggregations. Fisheries acoustics surveys conducted alongside multibeam mapping efforts inherently couple water column data with information on benthic habitats and provide information on the heterogeneity of both benthic habitats and biota in the water column. Building on this information serves to inform resource managers regarding how fishes are organized around habitat structure and the scale at which these relationships are important. Where resource managers require place-based assessments regarding the location of critical habitats along with high abundances of fish, concurrent multibeam and fisheries acoustics surveys serve as an important tool for characterizing and prioritizing coral reef ecosystems. This report summarizes the evolution of fisheries acoustics surveys onboard the NOAA ship Nancy Foster from 2008 to 2010, in conjunction with multibeam data collection, aimed at characterizing benthic and mid-water habitats in high priority conservation areas around the USVI and Puerto Rico. It also serves as a resource for the continued development of consistent data products derived from acoustic surveys. By focusing on the activities of 2010, this report highlights the progress made to date and illustrates the potential application of fisheries data derived from acoustic surveys to the management of coral reef ecosystems.

Comparative molecular cytogenetic studies in the order Carnivora: mapping chromosomal rearrangements onto the phylogenetic tree

We have made a set of chromosome-specific painting probes for the American mink by degenerate oligonucleotide primed-PCR (DOP-PCR) amplification of flow-sorted chromosomes. The painting probes were used to delimit homologous chromosomal segments among human, red fox, dog, cat and eight species of the family Mustelidae, including the European mink, steppe and forest polecats, least weasel, mountain weasel, Japanese sable, striped polecat, and badger. Based on the results of chromosome painting and G-banding, comparative maps between these species have been established. The integrated map demonstrates a high level of karyotype conservation among mustelid species. Comparative analysis of the conserved chromosomal segments among mustelids and outgroup species revealed 18 putative ancestral autosomal segments that probably represent the ancestral chromosomes, or chromosome arms, in the karyotype of the most recent ancestor of the family Mustelidae. The proposed 2n = 38 ancestral Mustelidae karyotype appears to have been retained in some modern mustelids, e.g., Martes, Lutra, ktonyx, and Vormela. The derivation of the mustelid karyotypes from the putative ancestral state resulted from centric fusions, fissions, the addition of heterochromatic arms, and occasional pericentric inversions. Our results confirm many of the evolutionary conclusions suggested by other data and strengthen the topology of the carnivore phylogenetic tree through the inclusion of genome-wide chromosome rearrangements. Copyright (C) 2002 S. KargerAG, Basel.

Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.