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.     

THE QUALITATIVE CONTENT OF RENEWABLE RESOURCE MODELS

This paper analyzes the extent to which standard dynamic renewable resource models possess refutable implications. Both the steady state comparative static and local comparative dynamic properties of the standard model are studied. A unified framework is developed which enables one to analyze the qualitative properties of any standard renewable resource model. This is achieved by explicitly linking the local stability, steady state comparative static, and local comparative dynamic properties of the model.