Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

On acoustic and electric waves diffraction in piezoelectric medium by a permeable half-plane crack

The problem of electric and acoustic waves diffraction by a half-plane crack in a transversal isotropic piezoelectric medium is investigated. The crack is assumed to be electric permeable and free of tractions. The so-called “quasi-hyperbolic approximation” [15] is adopted. Applying Laplace transformations and Wiener–Hopf technique a closed form solution is obtained. By the means of Cagniard–de Hoop method a detailed dynamic full electroacoustic wavefield’s investigation is conducted. Mode conversion between electric and acoustic waves, effect of electroacoustic head wave, Bleustein–Gulyaev surface wave and the wavefield structure depending on the type of the incident wave (acoustic or electric) and its angle of incidence are analyzed in details. The dynamic field intensity factors at the crack tip depending on the angle of incidence and on time are derived explicitly. Numerical analysis is presented.     

Essentially optimal explicit Runge–Kutta methods with application to hyperbolic–parabolic equations

Optimal explicit Runge–Kutta methods consider more stages in order to include a particular spectrum in their stability domain and thus reduce time-step restrictions. This idea, so far used mostly for real-line spectra, is generalized to more general spectra in the form of a thin region. In thin regions the eigenvalues may extend away from the real axis into the imaginary plane. We give a direct characterization of optimal stability polynomials containing a maximal thin region and calculate these polynomials for various cases. Semi-discretizations of hyperbolic–parabolic equations are a relevant application which exhibit a thin region spectrum. As a model, linear, scalar advection–diffusion is investigated. The second-order-stabilized explicit Runge–Kutta methods derived from the stability polynomials are applied to advection–diffusion and compressible, viscous fluid dynamics in numerical experiments. Due to the stabilization the time step can be controlled solely from the hyperbolic CFL condition even in the presence of viscous fluxes.     

An algorithm for decomposing a polynomial system into normal ascending sets

We present an algorithm to decompose a polynomial system into a finite set of normal ascending sets such that the set of the zeros of the polynomial system is the union of the sets of the regular zeros of the normal ascending sets.If the polynomial system is zero dimensional,the set of the zeros of the polynomials is the union of the sets of the zeros of the normal ascending sets.     

Estimates from Below for the Spectral Function and for the Remainder in Local Weyl’s Law

We obtain asymptotic lower bounds for the spectral function of the Laplacian and for the remainder in local Weyl’s law on manifolds. In the negatively curved case, thermodynamic formalism is applied to improve the estimates. Key ingredients of the proof include the wave equation parametrix, a pretrace formula and the Dirichlet box principle. Our results develop and extend the unpublished thesis of A. Karnaukh [Ka]. The first author was supported by NSERC, FQRNT, Alfred P. Sloan Foundation Fellowship and Dawson Fellowship. The second author was supported by NSERC and FQRNT. Received: April 2006 Revision: October 2006 Accepted: October 2006     

Machine assignment in a nonlinear multi-product flowline

A group of machines for processing a set of jobs in a manufacturing system is often located in a serial line. An efficient strategy for locating these machines such that the total travel distance or the cost of transporting the jobs is minimized is desired. In this research, the assumption of a linear line with equally spaced machine location is relaxed. This research addressed problems of locating unique machines. It is found that the machine distances possess unique properties in this type of a problem. Utilizing these properties, heuristic strategies are proposed to obtain efficient solution where optimal methods are expected to be computationally prohibitive. A lower bound for the optimum solution is also proposed. Results are encouraging.     

Some remarks on general commutators theorems

The paper proves some general facts on commutators that refer to Fuglede-Putnam classical theorem in the spectral theory of not necessarily selfadjoint operators.     

Constructing Space for Science at the Royal Institution of Great Britain

Those who have worked in the Royal Institution of Great Britain have, since its foundation in 1799, made significant contributions to scientific knowledge, to its practical application, and to its communication to a wide variety of audiences. Such work cannot be carried out in an architectural vacuum, and in this paper we examine how the buildings of the Royal Institution, 20 and 21 Albemarle Street in central London, have shaped the work undertaken within its walls and how, on a number of occasions, the buildings have been reconfigured to take account of the evolving needs of scientific research and communication. This paper is based on the Conservation Plan of the Royal Institution that we wrote during 2003. The Conservation Plan did not examine the land owned by the Royal Institution to the north (i.e., 22 and 23 Albemarle Street; for this area see Richard Garnier, “Grafton Street, Mayfair,” Georgian Group Journal 13 (2003), 210–272), but it did discuss 18 and 19 Albemarle Street. In this paper we concentrate on the core Royal Institution buildings at 20 and 21 Albemarle Street. Other studies of the relationship of architecture,space, and science include Crosbie Smith and Jon Agar, ed., Making Space for Science: Territorial Themes in the Shaping of Knowledge (Basingstoke: Macmillan, 1997); Peter Galison and Emily Thompson, ed., The Architecture of Science (Cambridge, Mass.: MIT Press, 1999); and Sophie Forgan,“The architecture of science and the idea of a university,” Studies in History and Philosophy of Science 20 (1989), 405–434. Frank A.J.L. James is Professor of the History of Science at the Royal Institution; he has written widely on the history of nineteenth-century science in its social and cultural contexts and is editor of the Correspondence of Michael Faraday. He is President of the British Society for the History of Science. Anthony Peers is an Associate of Rodney Melville and Partners where he works in the field of building conservation as an architectural historian. He is a Council member of the Ancient Monument Society.     

Symplectic Yang–Mills theory,Ricci tensor,and connections

A Yang–Mills theory in a purely symplectic framework is developed. The corresponding Euler–Lagrange equations are derived and first integrals are given. We relate the results to the work of Bourgeois and Cahen on preferred symplectic connections.     

Hospital Waiting List Management Models

This paper describes work carried out with hospital consultants in Lancaster District Health Authority to develop and use computerized waiting list management models. The ‘what if?’ models produced have been found to be a very useful tool in the constructive debate about problems faced by consultants and managers, and of the options available to them.