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.     

A dual graph-norm refinement indicator for finite volume approximations of the Euler equations

Summary. We generalise and apply a refinement indicator of the type originally designed by Mackenzie, Süli and Warnecke in [15] and [16] for linear Friedrichs systems to the Euler equations of inviscid, compressible fluid flow. The Euler equations are symmetrized by means of entropy variables and locally linearized about a constant state to obtain a symmetric hyperbolic system to which an a posteriori error analysis of the type introduced in [15] can be applied. We discuss the details of the implementation of the refinement indicator into the DLR--Code which is based on a finite volume method of box type on an unstructured grid and present numerical results. Received May 15, 1995 / Revised version received April 17, 1996     

Study of the charge-exchange reactions\pi ^ - p \to (\pi ^0 ,\eta ,\eta \prime )n at 63 GeV

Muon electron pairs were detected in an Al multiplate spark chamber, exposed to a neutrino beam from the CERN PS. The leptons were not accompanied by other particles, except occasionally by protons. The background came mainly from muon associated π⁰ production, with one decay gamma lost. It was determined empirically, together with the small contribution from υ _e reactions. For electron energies above 2 GeV the background is 5.7±1.5 events, whereas 18 (μe)-candidates have been observed. Hence the effect is established, with a rate of about 10^?4 as compared to the muonic reactions above 3 GeV. Charm creation as the origin of this (μe)-production process is excluded; heavy neutral lepton production does not fit the kinematics observed. Instead the events are compatible with the two-body decay of an object with variable invariant mass of order 1 GeV, possibly resulting from axion interactions.     

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.     

Sharp estimates in terms of moduli of smoothness for compound cubature rules

Summary In 1980 Dahmen-DeVore-Scherer introduced a modulus of continuity which turns out to reflect invariance properties of compound cubature rules effectively. Accordingly, sharp error bounds are derived, the existence of relevant counterexamples being a consequence of a quantitative resonance principle, established previously.     

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.     

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 [M. Büttiker,Z. Phys. B 68:161 (1987)]. A specific example of this behavior [N. G. van Kampen,IBM J. Res. Dev. 32:107 (1988); R. Landauer,J. Stat. Phys. 53:233 (1988).] is the motion of overdamped particles in a ring subject to a nonuniform temperature field. When the drift and the temperature, which are periodic with a period equal to the ring circumference, are not in phase a noise-induced circulating current results.This paper will appear in a forthcoming issue of theJournal of Statistical Physics.