Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

Transport as a consequence of state-dependent diffusion

Overdamped particles subject to a drift in a force field with sinusoidal space dependence and also a sinusoidally modulated space-dependent diffusion, with the same period as the drift, experience a net driving force. The resulting current depends on the amplitude of the modulation of the diffusion and is a periodic function of the phase difference between the sinusoidal drift and the sinusoidal modulation of the diffusion. For small modulation amplitudes a particle subject to state-dependent noise behaves the same way as a particle subject to thermal noise but with a drift which, in addition to the sinusoidal term, contains a net force term.Dedicated to Professor Harry Thomas on the occasion of his 60th birthday     

Analytic calculation of scaling dimensions: Tricritical hard squares and critical hard hexagons

The finite-size corrections, central chargesc, and scaling dimensionsx of tricritical hard squares and critical hard hexagons are calculated analytically. This is achieved by solving the special functional equation or inversion identity satisfied by the commuting row transfer matrices of these lattice models at criticality. The results are expressed in terms of Rogers dilogarithms. For tricritical hard squares we obtainc=7/10,x=3/40, 1/5, 7/8, 6/5 and for hard hexagons we obtainc=4/5,x=2/15, 4/5, 17/15, 4/3, 9/5, in accord with the predictions of conformal and modular invariance.     

Portable and Automatic Mössbauer Analysis

A portable Mössbauer spectrometer, developed for extraterrestrial applications, opens up new industrial applications of MBS. But for industrial applications, an available tool for fast data analysis is also required, and it should be easy to handle. The analysis of Mössbauer spectra and their parameters is a barrier for the popularity of this wide-applicable spectroscopic technique in industry. Based on experience, the analysis of a Mössbauer spectrum is time-consuming and requires the dedication of a specialist. However, the analysis of Mössbauer spectra, from the fitting to the identification of the sample phases, can be faster using by genetic algorithms, fuzzy logic and artificial neural networks. Industrial applications are very specific ones and the data analysis can be performed using these algorithms. In combination with an automatic analysis, the Mössbauer spectrometer can be used as a probe instrument which covers the main industrial needs for an on-line monitoring of its products, processes and case studies. Some of these real industrial applications will be discussed.     

Changes in the kurtosis parameter of super-Gaussian beams passing through a spherically aberrated lens

Based on the Collins formula and irradiance moments definition, the propagation of the kurtosis parameter of super-Gaussian beams through a spherically aberrated lens is studied. Detailed numerical results are given. It is shown that, as compared with aberration-free super-Gaussian beams, the profile of spherically aberrated super-Gaussian beams can be leptokurtic, mesokurtic or platykurtic, depending on the aberration and propagation distance. The results for the spherically aberrated Gaussian beams can be dealt with as a special case treated in this paper. Finally, the advantage of our method is pointed out.     

From a Nonlinear, Nonconvex Variational Problem to a Linear, Convex Formulation

   Abstract. We propose a general approach to deal with nonlinear, nonconvex variational problems based on a reformulation of the problem resulting in an optimization problem with linear cost functional and convex constraints. As a first step we explicitly explore these ideas to some one-dimensional variational problems and obtain specific conclusions of an analytical and numerical nature.