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.     

Accelerating Effect in Kerr-Newman-Kasuya Field

We have obtained expressions of the accelerating effect in Kerr-Newman Kasuya field. These expressionsinclude four parameters: mass m, angular momentum a, electric charge q, and magnetic charge φ. Furthermore we studyits special case (vi = 0). We get the following conclusion. In the gravitation field of souse mass with electric charge qand magnetic charge b, the acceleration of test particle has not only radial component but also transverse component.When θ = 0, the acceleration is minimum, and when θ = π/2, the acceleration is maximum. Furthermore, we discussthe effects of electric charge q and magnetic charge φ respectively.     

Structures and stability of Al<Subscript><Emphasis Type="Bold">7</Emphasis></Subscript>C and Al<Subscript><Emphasis Type="Bold">7</Emphasis></Subscript>N clusters

We report results of the atomic and electronic structures of Al₇C cluster using ab initio molecular dynamics with ultrasoft pseudopotentials and generalized gradient approximation. The lowest energy structure is found to be the one in which carbon atom occupies an interstitial position in Al₇ cluster. The electronic structure shows that the recent observation [Chem. Phys. Lett. 316, 31 (2000)] of magic behavior of Al₇C^- cluster is due to a large highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) gap which makes Al₇C^- chemically inert. These results have further led us to the finding of a new neutral magic cluster Al₇N which has the same number of valence electrons as in Al₇C^- and a large HOMO-LUMO gap of 1.99 eV. Further, calculations have been carried out on (Al₇N)₂ to study interaction between magic clusters. Received 28 July 2001     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

有限时滞微分系统依照双测度的实际稳定性

本文对具有有限时滞的泛函微分方程建立了关于依照两种测度的实际稳定性的Razumikhin型判定定理，其中未采用通常的辅助函数，且可运用多个含有状态变量x的部分变元的Lyapunov函数，得出部分变远实际稳定性的判定定理，从而改进了已有的结果。     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

Stability of steady nonisothermal flow in a vertical layer with permeable boundaries in the microconvection model

The problem of the stability of steady convective viscous incompressible fluid flow in a vertical layer with boundaries at different temperatures is considered in the presence of transverse blowing through the layer. The complete spectral problem is solved for a silicon melt. The neutral curve is constructed and the critical Grashof number is found. The numerical calculations show that the presence of transverse blowing significantly affects the flow stability. As compared with the Oberbeck-Boussinesq model, in the microconvection model the instability develops at lower wavenumbers.