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.     

On Series of Ordinals and Combinatorics

This paper deals mainly with generalizations of results in finitary combinatorics to infinite ordinals. It is well-known that for finite ordinals ∑_bT<αβ is the number of 2-element subsets of an α-element set. It is shown here that for any well-ordered set of arbitrary infinite order type α, ∑_bT<αβ is the ordinal of the set M of 2-element subsets, where M is ordered in some natural way. The result is then extended to evaluating the ordinal of the set of all n-element subsets for each natural number n ≥ 2. Moreover, series ∑_β<αf(β) are investigated and evaluated, where α is a limit ordinal and the function f belongs to a certain class of functions containing polynomials with natural number coefficients. The tools developed for this result can be extended to cover all infinite α, but the case of finite α appears to be quite problematic.     

Spectrophotometric analysis, thermal analysis and gravimetric determination of some metal ions with oxime and Schiff's base derivatives of N-furoylphenylhydroxylamine

The oxime and Schiff's base of N-furoylphenylhydroxylamine form complexes with Co²⁺, Cu²⁺, Zn²⁺ and Fe³⁺. These complexes were studied by means of the UV — Vis. spectra of nujol mulls, electron spin resonance (ESR) spectra and magnetic susceptibility to determine their stereochemistry. The complexes were characterized via elemental analyses, molar conductivities and thermogravimetric analyses. Organic reagents were used for the gravimetric determination of Co²⁺, Cu²⁺ and Fe³⁺ through the precipitation of their complexes. The compounds were also used for separation of a binary mixture of Cu²⁺ and Cu²⁺. Interferences were studied.     

Linearization of column responses by correlation chromatography. Separation improvement and reduction of undesired adsorption effects

Summary Experiments and simulations prove that correlation chromatography can greatly reduce the disadvantage of a non-linear response of the chromatographic column. A factor that has been accepted as being an important source of error in correlation or multiplex chromatography, has been shown not to be. Separations affected, improve dramatically when correlation chromatography is used, and a substantial amount of correlation noise only arises when there is a large difference in separation between a conventional chromatogram and a correlogram.A model has been developed for simulating these nonlinearities. It is shown that, especially for correlation chromatography, the simulation results match the practical measurements very well.     

Recursive partition and amalgamation with the exponential family: Theory and applications

The theory of tree-growing (RECPAM approach) is developed for outcome variables which are distributed as the canonical exponential family. The general RECPAM approach (consisting of three steps: recursive partition, pruning and amalgamation), is reviewed. This is seen as constructing a partition with maximal information content about a parameter to be predicted, followed by simplification by the elimination of ‘negligible’ information. The measure of information is defined for an exponential family outcome as a deviance difference, and appropriate modifications of pruning and amalgamation rules are discussed. It is further shown how the proposed approach makes it possible to develop tree-growing for situations usually treated by generalized linear models (GLIM). In particular, Poisson and logistic regression can be tree-structured. Moreover, censored survival data can be treated, as in GLIM, by observing a formal equivalence of the likelihood under random censoring and an appropriate Poisson model. Three examples are given of application to Poisson, binary and censored survival data.     

Dominant modes of submerged thin cylindrical shells

The relationship between the dominant mode of the submerged thin cylindrical shell and the flexural wave velocity is investigated. The natural frequency corresponding to the vibration mode is obtained as the solution of characteristic equation of thin cylindrical shell. However, it is difficult to estimate the dominant mode, especially if two or more vibration modes are involved. To estimate the dominant mode of a thin shell in vacuo, the concept of “modified bending stiffness” has been introduced. In this paper, the concept of modified bending stiffness is extended to estimate the dominant mode of a submerged thin cylindrical shell. The dominant mode of a submerged thin cylindrical shell is theoretically discriminated from the other mode based on the smallness of the modified bending stiffness of the submerged shell. The validity of our theory is confirmed by a good agreement between theoretical and experimental results on flexural wave velocity.     

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.     

EXISTENCE OF PERIODIC SOLUTIONS OF SYSTEM OF DIFFERENCE EQUATIONS

This paper studies the nonautonomous nonlinear system of difference equationsΔx(n)=A(n)x(n)+f(n,x(n)),n∈Z,(*) where x(n)∈R~N,A(n)=(a_(ij)(n))N×N is an N×N matrix,with a-(ij)∈C(R,R) for i,j= 1,2,3,...,N,and f=(f_1,f_2,...,f_N)~T∈C(R×R~N,R~N),satisfying A(t+ω)=A(t),f(t+ω,z)=f(t,z) for any t∈R,(t,z)∈R×R~N andωis a positive integer.Sufficient conditions for the existence ofω-periodic solutions to equations (*) are obtained.     

Mass Power Spectrum in a Universe Dominated by the Chaplygin Gas

The mass power spectrum for a Universe dominated by the Chaplygin gas is evaluated numerically from scales of the order of the Hubble horizon to 100 Mpc. The results are compared with a pure baryonic Universe and a cosmological constant model. In all three cases, the spectrum increases with k, the wavenumber of the perturbations. The slope of the spectrum is higher for the baryonic model and smaller for the cosmological constant model, the Chaplygin gas interpolating these two models. The results are analyzed in terms of the sound velocity of the Chaplygin gas and the moment the Universe begins to accelerate.     

Method of lines with boundary elements for 1‐D transient diffusion‐reaction problems

Time‐dependent differential equations can be solved using the concept of method of lines (MOL) together with the boundary element (BE) representation for the spatial linear part of the equation. The BE method alleviates the need for spatial discretization and casts the problem in an integral format. Hence errors associated with the numerical approximation of the spatial derivatives are totally eliminated. An element level local cubic approximation is used for the variable at each time step to facilitate the time marching and the nonlinear terms are represented in a semi‐implicit manner by a local linearization at each time step. The accuracy of the method has been illustrated on a number of test problems of engineering significance. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006