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.     

Investigation of solvability of the multidimensional two-phase Stefan and the nonstationary filtration Florin problems for second order parabolic equations in weighted Hölder spaces of functions

Multidimensional two-phase Stefan (k=1) and nonstationary filtration Florin (k=0) problems for second order parabolic equations in the case when the free boundary is a graph of a functionx _n =Ψ _k (x′t),x′∈ ^n?1 ,n≥2,t∈(0,T) are studied. A unique solvability theorem in weighted Hölder spaces of functions with time power weight is proved, coercive estimates for solutions are obtained. Bibliography: 30 titles.     

Logarithmic temperature dependence of the conductivity of the two-dimensional metal

We report on the first observation and studies of a weak delocalizing logarithmic temperature dependence of the conductivity, which causes the conductivity of the 2D metal to increase as T decreases down to 16 mK. The prefactor of the logarithmic dependence is found to decrease gradually with density, to vanish at a critical density n _c , 2∼2×10¹² cm⁻², and then to have the opposite sign at n>n _c ,2. The second critical density sets the upper limit on the existence region of the 2D metal, whereas the conductivity at the critical point, G _c ,2∼120e ²/h, sets an upper (low-temperature) limit on its conductivity. Pis’ma Zh. éksp. Teor. Fiz. 68, No. 6, 497–501 (25 September 1998) Published in English in the original Russian journal. Edited by Steve Torstveit.     

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