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. 相似文献
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 相似文献
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. 相似文献
We consider the problem of maximizing a linear fractional function on the Pareto efficient frontier of two other linear fractional functions. We present a finite pivoting-type algorithm that solves the maximization problem while computing simultaneously the efficient frontier. Application to multistage efficiency analysis is discussed. An example demonstrating the computational procedure is included. 相似文献
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. 相似文献
RSMR experiments with57Fe radiation were performed on myoglobin. An areasensitive detector was employed for simultaneous angular dependent collection of the scattered quanta up to a maximum angle 2θ of 17‡. Experimental data of polycrystalline and lyophilized myoglobin are compared with computer calculations of the scattering which are based on the atomic coordinates determined by X-ray structure analysis. Special attention has been paid to the influence of coherence effects from collectively moving parts of the protein. A simple model is introduced in order to take into account these segmental motions. Our first results indicate that the sizes of collectively moving segments are comparable with spheres of about 6 å in diameter in dry myoglobin. In myoglobin crystals, where the molecules are surrounded by large hydration shells, the movements appear to be correlated in segments with sizes comparable to helices.