Piecewise-linearized methods for oscillators with limit cycles

A piecewise linearization method based on the linearization of nonlinear ordinary differential equations in small intervals, that provides piecewise analytical solutions in each interval and smooth solutions everywhere, is developed for the study of the limit cycles of smooth and non-smooth, conservative and non-conservative, nonlinear oscillators. It is shown that this method provides nonlinear maps for the displacement and velocity which depend on the previous values through the nonlinearity and its partial derivatives with respect to time, displacement and velocity, and yields non-standard finite difference formulae. It is also shown by means of five examples that the piecewise linearization method presented here is more robust and yields more accurate (in terms of displacement, energy and frequency) solutions than the harmonic balance procedure, the method of slowly varying amplitude and phase, and other non-standard finite difference equations.     

Affine Arithmetic: Concepts and Applications

Affine arithmetic is a model for self-validated numerical computation that keeps track of first-order correlations between computed and input quantities. We explain the main concepts in affine arithmetic and how it handles the dependency problem in standard interval arithmetic. We also describe some of its applications.     

A modular integer GCD algorithm

This paper describes the first algorithm to compute the greatest common divisor (GCD) of two n-bit integers using a modular representation for intermediate values U, V and also for the result. It is based on a reduction step, similar to one used in the accelerated algorithm [T. Jebelean, A generalization of the binary GCD algorithm, in: ISSAC '93: International Symposium on Symbolic and Algebraic Computation, Kiev, Ukraine, 1993, pp. 111–116; K. Weber, The accelerated integer GCD algorithm, ACM Trans. Math. Softw. 21 (1995) 111–122] when U and V are close to the same size, that replaces U by (U−bV)/p, where p is one of the prime moduli and b is the unique integer in the interval (−p/2,p/2) such that . When the algorithm is executed on a bit common CRCW PRAM with O(nlognlogloglogn) processors, it takes O(n) time in the worst case. A heuristic model of the average case yields O(n/logn) time on the same number of processors.     

Amorphous chalcogenide AgSbS2 films prepared by pulsed laser deposition

The pulsed laser deposition technique has been applied to prepare amorphous ternary AgSbS films. The films were prepared from AgSbS₂ bulk glass using a KrF excimer laser. The composition of prepared films according to the results of the energy dispersive X-ray analysis and the Rutherford backscattering was close to bulk one. Optical transmission and spectral dependence of the refractive index proved the good optical quality of the films. The Ar⁺ ion laser dot exposures of the films show a potential applicability of the films as a new type of optical recording material. PACS 78.66.Jg; 81.15.Fg; 81.40.Wx     

SASEGASA: A New Generic Parallel Evolutionary Algorithm for Achieving Highest Quality Results

This paper presents a new generic Evolutionary Algorithm (EA) for retarding the unwanted effects of premature convergence. This is accomplished by a combination of interacting generic methods. These generalizations of a Genetic Algorithm (GA) are inspired by population genetics and take advantage of the interactions between genetic drift and migration. In this regard a new selection scheme is introduced, which is designed to directedly control genetic drift within the population by advantageous self-adaptive selection pressure steering. Additionally this new selection model enables a quite intuitive heuristics to detect premature convergence. Based upon this newly postulated basic principle the new selection mechanism is combined with the already proposed Segregative Genetic Algorithm (SEGA), an advanced Genetic Algorithm (GA) that introduces parallelism mainly to improve global solution quality. As a whole, a new generic evolutionary algorithm (SASEGASA) is introduced. The performance of the algorithm is evaluated on a set of characteristic benchmark problems. Computational results show that the new method is capable of producing highest quality solutions without any problem-specific additions.     

Mössbauer Study of Ceramic Finds from the Galería de las Ofrendas, Chavín de Huántar

Ceramic finds from the Galería de las Ofrendas at Chavín de Huántar and surface finds from the settlement of Chavín were characterised by combining the results of archaeological typology with archaeometric studies using neutron activation analysis, Mössbauer spectroscopy, X-ray diffraction and thin-section microscopy. Sherds from the pyramid Tello are included in the study as representative of local material. The analyses show that the vessels were made from different raw materials and that different firing procedures were used in their production. Sherds of certain styles largely exhibit similar types of Mössbauer patterns and in many instances also have similar element compositions. This supports the archaeological notion that the vessels were brought to Chavín from the provinces, perhaps on the occasion of a festivity.     

Kinetics of oxygen exchange in strontium titanate

The kinetics of oxygen exchange are of primary importance for the application of titanates as fast resistive oxygen sensors. The sensor’s conductivity is correlated with the oxygen partial pressure pO₂ of the surrounding atmosphere: Due to oxygen surface transfer and subsequent diffusion of oxygen vacancies V _O ^·· , a pO₂ change gives rise to a conductivity change of the sample. While bulk diffusion usually occurs very fast, the surface transfer reaction becomes the rate determining step for thin samples and for low temperatures. We have shown that in the case of acceptor doped SrTiO₃ the kinetics of the surface transfer reaction can be strongly influenced through stoichiometric changes brought about by thin coatings of alkaline earth metal oxides (e.g. SrO). In contrast to the commonly used jump method (conductivity response to a sudden pO₂ change in the time domain), a model is presented which is based on the frequency-domain analysis of amplitude and phase shift of the response signal obtained from a pO₂ modulation in a fast kinetic measurement set-up. This method allows not only for measuring response times in the sub-millisecond range but also for distinguishing between behaviour either controlled by volume diffusion or by surface transfer reaction. Paper presented at the 9th EuroConference on Ionics, Ixia, Rhodes, Greece, Sept. 15 – 21, 2002.