On a class of preconditioners for solving the Helmholtz equation

In 1983, a preconditioner was proposed [J. Comput. Phys. 49 (1983) 443] based on the Laplace operator for solving the discrete Helmholtz equation efficiently with CGNR. The preconditioner is especially effective for low wavenumber cases where the linear system is slightly indefinite. Laird [Preconditioned iterative solution of the 2D Helmholtz equation, First Year's Report, St. Hugh's College, Oxford, 2001] proposed a preconditioner where an extra term is added to the Laplace operator. This term is similar to the zeroth order term in the Helmholtz equation but with reversed sign. In this paper, both approaches are further generalized to a new class of preconditioners, the so-called “shifted Laplace” preconditioners of the form Δφ−k²φ with . Numerical experiments for various wavenumbers indicate the effectiveness of the preconditioner. The preconditioner is evaluated in combination with GMRES, Bi-CGSTAB, and CGNR.     

Error estimates for Hermite interpolation on spheres

In this paper, we prove convergence rates for spherical spline Hermite interpolation on the sphere S^d−1 via an error estimate given in a technical report by Luo and Levesley. The functionals in the Hermite interpolation are either point evaluations of pseudodifferential operators or rotational differential operators, the desirable feature of these operators being that they map polynomials to polynomials. Convergence rates for certain derivatives are given in terms of maximum point separation.     

Heavy Majorana neutrinos in e + e − → W + W−

In this paper we investigate the effects of heavy Majorana neutrinos in the reaction e ⁺ e ^? → W ⁺ W ^?. We consider neutrino masses in the 1–10 TeV region. We show that at LEP II and NLC energies it is possible to use this processes to verify indirect evidence of heavy neutral particles with mixing angles of the order sin² α = 0.01. We discuss the unitarity restrictions that can be obtained for vector singlet and fermion-mirror-fermion models.     

Transient binary nucleation from ethanol-hexanol vapour

We study condensation of ethanol-hexanol vapour by numerical solution of kinetic equations. The number of droplets formed in unit volume is computed within self-consistent classical model. It is shown that formation of ethanol-rich droplets prevails at the initial stage of nucleation process, but in the stationary state formation of droplets near the saddle point (on cluster formation energy surface) plays the dominant role. Presented at the 6th Joint Seminar “Development of Materials Science in Research and Education”, Karlštejn, Czech Republic, 17–19 September 1996. This work was supported by Grant No. A1010615 of the Grant Agency of the Academy of Sciences of the Czech Republic.     

Dissipative Nonlinear Evolution Equations and Chaos

In this article we have studied the nonlinear interaction between ellipticity and dissipation in a set of model equations (1.1) and established the relation between this interaction and chaos. In addition to theoretical investigations, extensive numerical simulations with these equations have been made, and different routes to chaos have been found. The numerical studies have revealed the chaotic nature of the solutions.     

Rarefied Flow Computations Using Nonlinear Model Boltzmann Equations

High resolution finite difference schemes for solving the nonlinear model Boltzmann equations are presented for the computations of rarefied gas flows. The discrete ordinate method is first applied to remove the velocity space dependency of the distribution function which renders the model Boltzmann equation in phase space to a set of hyperbolic conservation laws with source terms in physical space. Then a high order essentially nonoscillatory method due to Harten et al. (J. Comput. Phys. 71, 231, 1987) is adapted and extended to solve them. Explicit methods using operator splitting and implicit methods using the lower-upper factorization are described to treat multidimensional problems. The methods are tested for both steady and unsteady rarefied gas flows to illustrate its potential use. The computed results using model Boltzmann equations are found to compare well both with those using the direct simulation Monte Carlo results in the transitional regime flows and those with the continuum Navier-Stokes calculations in near continuum regime flows.