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.     

Stability testing in an integrated scheme

ISO Guide 35 deals with RM stability issues and scrutinizes the evaluation of stability testing results under the assumption that either there is no trend at all (a rather rare situation), or any observed deterministic change is insignificant and thus can be neglected. However, market demands for reliable reference materials are obviously not limited to stable or at least seemingly stable materials. In many analytical applications, analytes and measurands under consideration are known, or at least suspected, to be unstable on time scales that may vary widely from measurand to measurand. The Federal Institute for Materials Research and Testing (BAM) has developed (and successfully uses) an integrated approach in its certification practice. The approach is based on an initial stability study and subsequent post-certification monitoring. Data evaluation is model-based and takes advantage of all information collected in the stability testing scheme(s). It thus allows one to deal with any kind of instability observed, to assess limiting time intervals at any stress condition in the range tested, to estimate a final expiry date for materials with detected instabilities or the maximum admissible re-testing interval for seemingly stable materials, and to assess maximum admissible stress loads during delivery of the material to the customer. The article describes (and exemplifies) typical study layout, the model selection, and the integrated data assessment.     

Molecular flow and wall collision age distributions

The use of storage cells has become a standard technique for internal gas targets in conjunction with high energy storage rings. In case of spin-polarized hydrogen and deuterium gas targets the interaction of the injected atoms with the walls of the storage cell can lead to depolarization and recombination. Thus the number of wall collisions of the atoms in the target gas is important for modeling the processes of spin relaxation and recombination. It is shown in this article that the diffusion process of rarefied gases in long tubes or storage cells can be described with the help of the one-dimensional diffusion equation. Mathematical methods are presented that allow one to calculate collision age distributions (CAD) and their moments analytically. These methods provide a better understanding of the different aspects of diffusion than Monte Carlo calculations. Additionally it is shown that measurements of the atomic density or polarization of a gas sample taken from the center of the tube allow one to determine the possible range of the corresponding density weighted average values along the tube. The calculations are applied to the storage cell geometry of the HERMES internal polarized hydrogen and deuterium gas target. Received 9 July 2001 and Received in final form 18 September 2001     

Dynamic NMR study of the flexibility of 2-amino-3-aroyl-4,6-diaryl-pyrylium salts

Dynamic NMR investigations of a number of 2-amino-3-aroyl-4,6-diaryl-pyrylium salts were carried out. The barrier to rotation of the partial C, N double bond was determined and proved to be in the range of 62 to 63 kJ/mol. Quantum chemical calculations of bond orders and electron densities of the different atoms in the molecules show the distinct double bond character of the exocyclic C, N bond. This is in agreement with the relatively high barrier to rotation. By quantum chemical ab initio 3-21G calculations, the dynamic behaviour of this kind of compounds was simulated; two pairs (image and mirror image) of ground state conformations, in coincidence with the experiment, were obtained. Received: 10 May 1996 / Revised: 1 July 1996 / Accepted: 4 July 1996