A natural extension of the conventional finite volume method into polygonal unstructured meshes for CFD application.

A new general cell-centered solution procedure based upon the conventional control or finite volume (CV or FV) approach has been developed for numerical heat transfer and fluid flow which encompasses both structured and unstructured meshes for any kind of mixed polygon cell. Unlike conventional FV methods for structured and block structured meshes and both FV and FE methods for unstructured meshes, the irregular control volume (ICV) method does not require the shape of the element or cell to be predefined because it simply exploits the concept of fluxes across cell faces. That is, the ICV method enables meshes employing mixtures of triangular, quadrilateral, and any other higher order polygonal cells to be exploited using a single solution procedure. The ICV approach otherwise preserves all the desirable features of conventional FV procedures for a structured mesh; in the current implementation, collocation of variables at cell centers is used with a Rhie and Chow interpolation (to suppress pressure oscillation in the flow field) in the context of the SIMPLE pressure correction solution procedure. In fact all other FV structured mesh-based methods may be perceived as a subset of the ICV formulation. The new ICV formulation is benchmarked using two standard computational fluid dynamics (CFD) problems i.e., the moving lid cavity and the natural convection driven cavity. Both cases were solved with a variety of structured and unstructured meshes, the latter exploiting mixed polygonal cell meshes. The polygonal mesh experiments show a higher degree of accuracy for equivalent meshes (in nodal density terms) using triangular or quadrilateral cells; these results may be interpreted in a manner similar to the CUPID scheme used in structured meshes for reducing numerical diffusion for flows with changing direction.

Modelling of industrial multi-physics processes—a key role for computational mechanics.

In the analysis of industrial processes, there is an increasing emphasis on systems governed by interacting continuum phenomena. Mathematical models of such multi-physics processes can only be achieved for practical simulations through computational solution procedures—computational mechanics. Examples of such multi-physics systems in the context of metals processing are used to explore some of the key issues. Finite-volume methods on unstructured meshes are proposed as a means to achieve efficient rapid solutions to such systems. Issues associated with the software design, the exploitation of high performance computers, and the concept of the virtual computational-mechanics modelling laboratory are also addressed in this context.

Molecular dynamics simulation of observed c(4x4) and c(4x3) C60 alignments on the Si(100) reconstructed surface.

We have implemented a large-scale classical molecular dynamics simulation at constant temperature to provide a theoretical insight into the results of a recently performed experiment on the monolayer and multi-layer formations of molecular films on the Si(100) reconstructed dimerized surface. Our simulation has successfully reproduced all of the morphologies observed on the monolayer film by this experiment. We have obtained the formation of both c(4 4) and c(4 3) structures of the molecules and have also obtained phase transitions of the former into the latter.

[Book review] Symmetry - A Unifying Concept, István Hargittai and Magdolna Hargittai

Tony Mann provides a review of the book: István Hargittai and Magdolna Hargittai, Symmetry - A Unifying Concept, 1994, Ten Speed Press, 222pp. ISBN: 089815-590-8. £14.99

[Book review] Iterations: The New Image

Tony Mann provides a review of the book: Timothy Druckrey (ed.) Iterations: The New Image. MIT Press, (1994), 199pp. ISBN: 0-2620-4143-X. £35.95

[Book review] Colour: Art and Science

Tony Mann provides a review of the book: Trevor Lamb & Janine Bourriau (Eds.) Colour: Art and Science, 1995. (Darwin College Lectures), Cambridge University Press, 237pp. ISBN: 0521-49645-4 (hbk.) 0521-49963-1 (pbk.)

Does an incremental filling technique reduce polymerization shrinkage stresses

It is widely accepted that volumetric contraction and solidification during the polymerization process of restorative composites in combination with bonding to the hard tissue result in stress transfer and inward deformation of the cavity walls of the restored tooth. Deformation of the walls decreases the size of the cavity during the filling process. This fact has a profound influence on the assumption-raised and discussed in this paper-that an incremental filling technique reduces the stress effect of composite shrinkage on the tooth. Developing stress fields for different incremental filling techniques are simulated in a numerical analysis. The analysis shows that, in a restoration with a well-established bond to the tooth-as is generally desired-incremental filling techniques increase the deformation of the restored tooth. The increase is caused by the incremental deformation of the preparation, which effectively decreases the total amount of composite needed to fill the cavity. This leads to a higher-stressed tooth-composite structure. The study also shows that the assessment of intercuspal distance measurements as well as simplifications based on generalization of the shrinkage stress state cannot be sufficient to characterize the effect of polymerization shrinkage in a tooth-restoration complex. Incremental filling methods may need to be retained for reasons such as densification, adaptation, thoroughness of cure, and bond formation. However, it is very difficult to prove that incrementalization needs to be retained because of the abatement of shrinkage effects.

Computer aided parallelisation tools (CAPTools) — conceptual overview and performance on the parallelisation of structured mesh codes

Computer Aided Parallelisation Tools (CAPTools) is a toolkit designed to automate as much as possible of the process of parallelising scalar FORTRAN 77 codes. The toolkit combines a very powerful dependence analysis together with user supplied knowledge to build an extremely comprehensive and accurate dependence graph. The initial version has been targeted at structured mesh computational mechanics codes (eg. heat transfer, Computational Fluid Dynamics (CFD)) and the associated simple mesh decomposition paradigm is utilised in the automatic code partition, execution control mask generation and communication call insertion. In this, the first of a series of papers [1–3] the authors discuss the parallelisations of a number of case study codes showing how the various component tools may be used to develop a highly efficient parallel implementation in a few hours or days. The details of the parallelisation of the TEAMKE1 CFD code are described together with the results of three other numerical codes. The resulting parallel implementations are then tested on workstation clusters using PVM and an i860-based parallel system showing efficiencies well over 80%.

Exploitation of symbolic information in interprocedural dependence analysis

The requirement for a very accurate dependence analysis to underpin software tools to aid the generation of efficient parallel implementations of scalar code is argued. The current status of dependence analysis is shown to be inadequate for the generation of efficient parallel code, causing too many conservative assumptions to be made. This paper summarises the limitations of conventional dependence analysis techniques, and then describes a series of extensions which enable the production of a much more accurate dependence graph. The extensions include analysis of symbolic variables, the development of a symbolic inequality disproof algorithm and its exploitation in a symbolic Banerjee inequality test; the use of inference engine proofs; the exploitation of exact dependence and dependence pre-domination attributes; interprocedural array analysis; conditional variable definition tracing; integer array tracing and division calculations. Analysis case studies on typical numerical code is shown to reduce the total dependencies estimated from conventional analysis by up to 50%. The techniques described in this paper have been embedded within a suite of tools, CAPTools, which combines analysis with user knowledge to produce efficient parallel implementations of numerical mesh based codes.

Automatic parallel code generation for message passing on distributed memory systems

The availability of a very accurate dependence graph for a scalar code is the basis for the automatic generation of an efficient parallel implementation. The strategy for this task which is encapsulated in a comprehensive data partitioning code generation algorithm is described. This algorithm involves the data partition, calculation of assignment ranges for partitioned arrays, addition of a comprehensive set of execution control masks, altering loop limits, addition and optimisation of communications for all data. In this context, the development and implementation of strategies to merge communications wherever possible has proved an important feature in producing efficient parallel implementations for numerical mesh based codes. The code generation strategies described here are embedded within the Computer Aided Parallelisation tools (CAPTools) software as a key part of a toolkit for automating as much as possible of the parallelisation process for mesh based numerical codes. The algorithms used enables parallelisation of real computational mechanics codes with only minor user interaction and without any prior manual customisation of the serial code to suit the parallelisation tool.