A geometric mass-preserving redistancing scheme for the level set function

In this paper we describe and evaluate a geometric mass-preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element-based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)-based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE-based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well-suited for the highly stretched curvilinear grids used in CFD simulations. Copyright (C) 2010 John Wiley & Sons, Ltd.

Crop rotation scheduling with adjacency constraints

In this article we propose a 0-1 optimization model to determine a crop rotation schedule for each plot in a cropping area. The rotations have the same duration in all the plots and the crops are selected to maximize plot occupation. The crops may have different production times and planting dates. The problem includes planting constraints for adjacent plots and also for sequences of crops in the rotations. Moreover, cultivating crops for green manuring and fallow periods are scheduled into each plot. As the model has, in general, a great number of constraints and variables, we propose a heuristics based on column generation. To evaluate the performance of the model and the method, computational experiments using real-world data were performed. The solutions obtained indicate that the method generates good results.

Sustainable vegetable crop supply problem

We consider an agricultural production problem, in which one must meet a known demand of crops while respecting ecologically-based production constraints. The problem is twofold: in order to meet the demand, one must determine the division of the available heterogeneous arable areas in plots and, for each plot, obtain an appropriate crop rotation schedule. Rotation plans must respect ecologically-based constraints such as the interdiction of certain crop successions, and the regular insertion of fallows and green manures. We propose a linear formulation for this problem, in which each variable is associated with a crop rotation schedule. The model may include a large number of variables and it is, therefore, solved by means of a column-generation approach. We also discuss some extensions to the model, in order to incorporate additional characteristics found in field conditions. A set of computational tests using instances based on real-world data confirms the efficacy of the proposed methodology. (C) 2009 Elsevier B.V. All rights reserved.

Partial spectral projected gradient method with active-set strategy for linearly constrained optimization

A method for linearly constrained optimization which modifies and generalizes recent box-constraint optimization algorithms is introduced. The new algorithm is based on a relaxed form of Spectral Projected Gradient iterations. Intercalated with these projected steps, internal iterations restricted to faces of the polytope are performed, which enhance the efficiency of the algorithm. Convergence proofs are given and numerical experiments are included and commented. Software supporting this paper is available through the Tango Project web page: http://www.ime.usp.br/similar to egbirgin/tango/.

Job rotation in assembly lines employing disabled workers

In this paper we consider the programming of job rotation in the assembly line worker assignment and balancing problem. The motivation for this study comes from the designing of assembly lines in sheltered work centers for the disabled, where workers have different task execution times. In this context, the well-known training aspects associated with job rotation are particularly desired. We propose a metric along with a mixed integer linear model and a heuristic decomposition method to solve this new job rotation problem. Computational results show the efficacy of the proposed heuristics. (C) 2009 Elsevier B.V. All rights reserved.

A Parallel Hardware Architecture for Scale and Rotation Invariant Feature Detection

This paper proposes a parallel hardware architecture for image feature detection based on the Scale Invariant Feature Transform algorithm and applied to the Simultaneous Localization And Mapping problem. The work also proposes specific hardware optimizations considered fundamental to embed such a robotic control system on-a-chip. The proposed architecture is completely stand-alone; it reads the input data directly from a CMOS image sensor and provides the results via a field-programmable gate array coupled to an embedded processor. The results may either be used directly in an on-chip application or accessed through an Ethernet connection. The system is able to detect features up to 30 frames per second (320 x 240 pixels) and has accuracy similar to a PC-based implementation. The achieved system performance is at least one order of magnitude better than a PC-based solution, a result achieved by investigating the impact of several hardware-orientated optimizations oil performance, area and accuracy.

Rotation effect on geodesic and zonal flow modes in tokamak plasmas with isothermal magnetic surfaces

The electrostatic geodesic mode oscillations are investigated in rotating large aspect ratio tokamak plasmas with circular isothermal magnetic surfaces. The analysis is carried out within the magnetohydrodynamic model including heat flux to compensate for the non-adiabatic pressure distribution along the magnetic surfaces in plasmas with poloidal rotation. Instead of two standard geodesic modes, three geodesic continua are found. The two higher branches of the geodesic modes have a small frequency up-shift from ordinary geodesic acoustic and sonic modes due to rotation. The lower geodesic continuum is a newzonal flowmode (geodesic Doppler mode) in plasmas with mainly poloidal rotation. Limits to standard geodesic modes are found. Bifurcation of Alfven continuum by geodesic modes at the rational surfaces is also discussed. Due to that, the frequency of combined geodesic continuum extends from the poloidal rotation frequency to the ion-sound band that can have an important role in suppressing plasma turbulence.

Temporal behaviour of toroidal rotation velocity in the TCABR tokamak

A new method for determining the temporal evolution of plasma rotation is reported in this work. The method is based upon the detection of two different portions of the spectral profile of a plasma impurity line, using a monochromator with two photomultipliers installed at the exit slits. The plasma rotation velocity is determined by the ratio of the two detected signals. The measured toroidal rotation velocities of C III (4647.4 angstrom) and C VI (5290.6 angstrom), at different radial positions in TCABR discharges, show good agreement, within experimental uncertainty, with previous results (Severo et al 2003 Nucl. Fusion 43 1047). In particular, they confirm that the plasma core rotates in the direction opposite to the plasma current, while near the plasma edge (r/a > 0.9) the rotation is in the same direction. This technique was also used to investigate the dependence of toroidal rotation on the poloidal position of gas puffing. The results show that there is no dependence for the plasma core, while for plasma edge (r/a > 0.9) some dependence is observed.

Systematics of nuclear densities, deformations and excitation energies within the context of the generalized rotation-vibration model

We present a large-scale systematics of charge densities, excitation energies and deformation parameters For hundreds of heavy nuclei The systematics is based on a generalized rotation vibration model for the quadrupole and octupole modes and takes into account second-order contributions of the deformations as well as the effects of finite diffuseness values for the nuclear densities. We compare our results with the predictions of classical surface vibrations in the hydrodynamical approximation. (C) 2010 Elsevier B V All rights reserved.

A complex network-based approach for boundary shape analysis

This paper introduces a novel methodology to shape boundary characterization, where a shape is modeled into a small-world complex network. It uses degree and joint degree measurements in a dynamic evolution network to compose a set of shape descriptors. The proposed shape characterization method has all efficient power of shape characterization, it is robust, noise tolerant, scale invariant and rotation invariant. A leaf plant classification experiment is presented on three image databases in order to evaluate the method and compare it with other descriptors in the literature (Fourier descriptors, Curvature, Zernike moments and multiscale fractal dimension). (C) 2008 Elsevier Ltd. All rights reserved.

Tunneling among rotation tori

We consider an integrable Hamiltonian system generated by the resonant normal form in order to study a particular mechanism of tunneling. We isolated near doublets of energy corresponding to rotation tori of the classical dynamics counterpart and the degeneracies breakdown is attributed to rotation-rotation tunneling. (C) 2008 Elsevier B.V. All rights reserved.

Paleomagnetic data and K-Ar ages from Mesozoic units of the Santa Marta massif: A preliminary interpretation for block rotation and translations

We report 6 K-Ar ages and paleomagnetic data from 28 sites collected in Jurassic, Lower Cretaceous and Paleocene rocks of the Santa Marta massif, to test previous hypothesis of rotations and translations of this massif, whose rock assemblage differs from other basement-cored ranges adjacent to the Guyana margin. Three magnetic components were identified in this study. A first component has a direction parallel to the present magnetic field and was uncovered in all units (D 352, I = 25.6, k = 57.35, a95 = 5.3, N = 12). A second component was isolated in Cretaceous limestone and Jurassic volcaniclastic rocks (D = 8.8, I = 8.3, k = 24.71, a95 = 13.7, N = 6), and it was interpreted as of Early Cretaceous age. In Jurassic sites with this component, Early Cretaceous K-Ar ages obtained from this and previous studies are interpreted as reset ages. The third component was uncovered in eight sites of Jurassic volcaniclastic rocks, and its direction indicates negative shallow to moderate inclinations and northeastward declinations. K-Ar ages in these sites are of Early (196.5 +/- 4.9 Ma) to early Late Jurassic age (156.6 +/- 8.9 Ma). Due to local structural complexity and too few Cretaceous outcrops to perform a reliable unconformity test, we only used two sites with (1) K-Ar ages, (2) less structural complexity, and (3) reliable structural data for Jurassic and Cretaceous rocks. The mean direction of the Jurassic component is (D = 20.4, I = -18.2, k = 46.9, a95 = 5.1, n = 18 specimens from two sites). These paleomagnetic data support previous models of northward along-margin translations of Grenvillian-cored massifs. Additionally, clockwise vertical-axis rotation of this massif, with respect to the stable craton, is also documented; the sense of rotation is similar to that proposed for the Perija Range and other ranges of the southern Caribbean margin. More data is needed to confirm the magnitudes of rotations and translations. (C) 2009 Elsevier Ltd. All rights reserved.

Continuous GRASP with a local active-set method for bound-constrained global optimization

Global optimization seeks a minimum or maximum of a multimodal function over a discrete or continuous domain. In this paper, we propose a hybrid heuristic-based on the CGRASP and GENCAN methods-for finding approximate solutions for continuous global optimization problems subject to box constraints. Experimental results illustrate the relative effectiveness of CGRASP-GENCAN on a set of benchmark multimodal test functions.

Algorithms for Maximum Independent Set in Convex Bipartite Graphs

A bipartite graph G = (V, W, E) is convex if there exists an ordering of the vertices of W such that, for each v. V, the neighbors of v are consecutive in W. We describe both a sequential and a BSP/CGM algorithm to find a maximum independent set in a convex bipartite graph. The sequential algorithm improves over the running time of the previously known algorithm and the BSP/CGM algorithm is a parallel version of the sequential one. The complexity of the algorithms does not depend on |W|.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.