Field-of-values analysis of preconditioned iterative methods for nonsymmetric elliptic problems

Summary. The convergence rate of Krylov subspace methods for the solution of nonsymmetric systems of linear equations, such as GMRES or FOM, is studied. Bounds on the convergence rate are presented which are based on the smallest real part of the field of values of the coefficient matrix and of its inverse. Estimates for these quantities are available during the iteration from the underlying Arnoldi process. It is shown how these bounds can be used to study the convergence properties, in particular, the dependence on the mesh-size and on the size of the skew-symmetric part, for preconditioners for finite element discretizations of nonsymmetric elliptic boundary value problems. This is illustrated for the hierarchical basis and multilevel preconditioners which constitute popular preconditioning strategies for such problems. Received May 3, 1996     

The influence of temperature gradients on the kinetic boundary layer problem for a condensing droplet

We extend an earlier method for solving kinetic boundary layer problems to the case of particles moving in aspatially inhomogeneous background. The method is developed for a gas mixture containing a supersaturated vapor and a light carrier gas from which a small droplet condenses. The release of heat of condensation causes a temperature difference between droplet and gas in the quasistationary state; the kinetic equation describing the vapor is the stationary Klein-Kramers equation for Brownian particles diffusing in a temperature gradient. By means of an expansion in Burnett functions, this equation is transformed into a set of coupled algebrodifferential equations. By numerical integration we construct fundamental solutions of this equation that are subsequently combined linearly to fulfill appropriate mesoscopic boundary conditions for particles leaving the droplet surface. In view of the intrinsic numerical instability of the system of equations, a novel procedure is developed to remove the admixture of fast growing solutions to the solutions of interest. The procedure is tested for a few model problems and then applied to a slightly simplified condensation problem with parameters corresponding to the condensation of mercury in a background of neon. The effects of thermal gradients and thermodiffusion on the growth rate of the droplet are small (of the order of 1%), but well outside of the margin of error of the method.     

Intrinsic thermal stability and quenching recovery of thin-film superconductors with thermal boundary resistance

The stability behaviour of a thin-film superconductor under a localized release of thermal disturbance is investigated. Two-dimensional conjugate film/substrate conduction equation with anisotropic thermal conductivity of the film, and Joule heat are employed to investigate effects of substrate and thermal properties on the intrinsic stability and quenching recovery. To consider the thermal boundary resistance between film and substrate, an interfacial-layer model (ILM) with very low diffusivity and an acoustic mismatch model (AMM) are employed. Results show that the thermal boundary resistance influences strongly the intrinsic stability. Thermal boundary resistance increases intrinsic stability if the thermal conductivity of the substrate or the disturbance energy is large. Higher Biot numbers and thermal conductivity ratios of film to substrate in longitudinal direction influence stability favorably. We demonstrate also that operation of a film/substrate system, such as YBCO/MgO, is either intrinsically stable or irrecoverably unstable.The authors wish to express their sincere appreciation to Dr. R. C. Chen for his invaluable advice and suggestions during the course of this paper. This research was supported by the National Science Council of the R. O. C. through grant NSC 83-0401-E-009-006. The computations were performed on the IBM ES/9000 at the National Center For High-Performance Computing.     

Viscosity solution theory of a class of nonlinear degenerate parabolic equations II. Lipschitz continuity of free boundary

In [1] we construct a unique bounded Hölder continuous viscosity solution for the nonlinear PDEs with the evolutionp-Laplacian equation and its anisotropic version as typical examples. In this part, we investigate the Lipschitz continuity of the free boundary of viscosity solution and its asymptotic spherical symmetricity, however, this result does not include the anisotropic case.This research is supported by the National Natural Sciences Foundation of China.     

Some results from 50 years' research on surface forces

A review is presented about research on surface forces and surface interactions conducted over the past half-century, with some emphasis on the pioneering contributions of the Department of Surface Phenomena at the Institute of Physical Chemistry of the USSR Academy of Sciences.     

具有非局部边界条件的奇摄动反应扩散问题

本文研究了一类具有非局部边界条件的奇摄动反应扩散初始边值问题．在适当的条件下，利用比较定理讨论了问题解的渐近性态．     

Cyclic reduction and FACR methods for piecewise hermite bicubic orthogonal spline collocation

Cyclic reduction and Fourier analysis-cyclic reduction (FACR) methods are presented for the solution of the linear systems which arise when orthogonal spline collocation with piecewise Hermite bicubics is applied to boundary value problems for certain separable partial differential equations on a rectangle. On anN×N uniform partition, the cyclic reduction and Fourier analysis-cyclic reduction methods requireO(N ²log₂ N) andO(N ²log₂log₂ N) arithmetic operations, respectively.     

Riabouchinsky flow with partially penetrable obstacle

An explicitly solvable Riabouchinsky model with a partially penetrable obstacle is introduced. This model applied to the estimation of the efficiency of free flow turbines allows us to take into account the pressure drop past the lamina.     

Trefftz,collocation, and other boundary methods—A comparison

In this article we survey the Trefftz method (TM), the collocation method (CM), and the collocation Trefftz method (CTM). We also review the coupling techniques for the interzonal conditions, which include the indirect Trefftz method, the original Trefftz method, the penalty plus hybrid Trefftz method, and the direct Trefftz method. Other boundary methods are also briefly described. Key issues in these algorithms, including the error analysis, are addressed. New numerical results are reported. Comparisons among TMs and other numerical methods are made. It is concluded that the CTM is the simplest algorithm and provides the most accurate solution with the best numerical stability. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A nodal coupling of finite and boundary elements

This article presents and analyzes a simple method for the exterior Laplace equation through the coupling of finite and boundary element methods. The main novelty is the use of a smooth parametric artificial boundary where boundary elements fit without effort together with a straight approximate triangulation in the bounded area, with the coupling done only in nodes. A numerically integrated version of the algorithm is also analyzed. Finally, an isoparametric variant with higher order is proposed. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 555–570, 2003