The optimisation of higher order resampling methods in a multiprocessor environment

The processing of images obtained from satellites often involves highly repetitive calculations on very large amounts of data. This processing is extremely time consuming when these calculations are performed on sequential machines. Parallel computers are well suited to handling computationally expensive operations such as higher order interpolations on large data sets. This paper decribes work undertaken to develop parallel implementations of a set of resampling procedures on an Alliant VFX/4. Each resampling procedure implemented has been optimised in three stages. First, the algorithm has been restructured, where two-dimensional resampling is performed by two one-dimensional resampling operations. Second, each procedure has been reprogrammed in such a way that the autoparallelisation provided by the FX/Fortran compiler has been exploited. Thirdly, data dependent analysis of each procedure has been performed in order to achieve full optimization of each procedure; each procedure has been restructured where appropriate to circumvent vectorisation and concurrency inhibiting data dependencies. The nature and extent of the code optimization achieved for each procedure is presented in this paper. The original code for the most computationally expensive procedure, as targeted at a sequential machine, was found to have an execution time of 4900 seconds on the Alliant VFX/4 when compiled with regular compiler optimization options. Following algorithmic redesigning and reprogramming of the code, as indicated in stage 1 and stage 2, the execution time was reduced to 248 seconds. Restructuring of the code following data dependency analysis indicated in stage 3 in order to avoid data dependencies and allow concurrency and vectorisation, further reduced execution time to 162 seconds. The consequence of this work is that higher-order resampling methods which had not previous been practical are now routinely performed on the Alliant VFX/4 at the University of Dundee.     

Calcium phosphate/calcium oxalate crystal association in urinary stones: implications for heterogeneous nucleation of calcium oxalate

PURPOSE: Most stones contain more than one type of crystals, and some combinations, such as calcium phosphate/calcium oxalate, are more common than others. Epitaxy between the crystals has been suggested to play a role in growth of such stones. The specific aim of this study is to investigate the involvement of calcium phosphate in crystallization of calcium oxalate. MATERIALS AND METHODS: Twenty calcium oxalate stones or stone fragments were examined using various microscopic techniques, including scanning, transmission and back-scattered electron microscopy. Similarly, calcium oxalate stones induced on a plastic foreign body implanted inside urinary bladders of laboratory rats were also investigated. Examination of the interface between calcium phosphate and calcium oxalate crystals was emphasized. RESULTS: Close association between crystals of calcium phosphate and calcium oxalate were found in both the human and rat stones. All crystals examined were associated with an organic matrix on the surface and contained copious amounts of organic material within the crystalline entities. Interface between the crystals also appeared to be occupied by organic matrix. CONCLUSIONS: Results of this and other studies from our laboratory indicate that epitaxy between various crystals, even though theoretically possible, appears unlikely in vivo. The appearance of specific crystalline combinations in stones is probably a result of the urinary environment being conducive for crystallization of those components. Heterogeneous nucleation of calcium oxalate is most probably induced by biological elements, including membranous cellular degradation products.     

Compact WLAN disc antennas

A novel conical beam patch antenna design, suitable for local area network applications, uses a central cylindrical connection from the disc to ground to achieve a large reduction in resonant size. Results on efficiency, bandwidth, and pattern are reported for several prototypes.     

BATCH SEQUENCING

Consider the single machine scheduling problem where there are a number of part types to be processed. A part type is defined as follows: Two parts are of the same part type if the machine does not require a setup in between the processing of these parts. The problem investigated in this paper is to find a sequence of batches of parts (if there are any) where all the requirements for parts are met. A heuristic and an exact algorithm are developed, and computational analysis is performed to measure the performance of the heuristic. The time complexity function of the heuristic is O(n²), and the exact algorithm runs in polynomial time given a fixed upper bound on the number of setups.     

Rheology of concentrated solutions of hyperbranched polyesters

The solution rheology of different generations of hyperbranched polyesters in N‐methyl‐2‐pyrrolidinone (NMP) solvent was examined in this study. The solutions exhibited Newtonian behavior over a wide range of polyester concentrations. Also, the relative viscosities of poly(amidoamine) (PAMAM) dendrimers in ethylenediamine were compared with those of the hyperbranched polyesters in NMP. Both types of dendritic polymers have relative viscosities that are exponential functions of their molar fraction in solution. The slopes of these relative viscosity curves show a linear relationship with respect to the generation number. PAMAM dendrimers have the greater slopes for each generation, reflecting their relatively larger intrinsic viscosity values.     

Validation of a Three-Dimensional Computational Fluid Dynamics Model of a Contact Tank

A three-dimensional (3D) computational fluid dynamics (CFD) model of a contact tank is presented in this paper. The model results are compared against 3D velocities and flow through curve (FTC) data, representing a tracer concentration profile, from a 1:8 scale physical model. The objective is to demonstrate that CFD models can simulate both the FTC and the 3D velocity field quite well. Simultaneous validation of velocities and FTC is important in ascertaining the predictive capabilities of CFD models, as physical model studies indicate that different baffle arrangements can lead to similar FTCs. Therefore, a good prediction of only FTC, as presented in previous 3D CFD model studies, does not necessarily imply a correct simulation of the flow field.     

Determination of gas dispersion in vapor extraction of heavy oil and bitumen

In this work, a mathematical model is developed and simulated to determine gas dispersion along with solubility during the vapor extraction (Vapex) of live oil from a laboratory scale physical model. The physical model is a rectangular block of homogenous porous medium saturated with heavy oil and bitumen. At a given temperature and pressure, the block is initially exposed on its side to a solvent gas, which diffuses into the medium and gets absorbed. The absorption of gas reduces the viscosity of heavy oil and bitumen causing it to drain under gravity. The low-viscosity “live oil” is produced at the bottom of the porous block. The production of live oil with time is accompanied by the shrinkage of oil in the block as well as its increased exposure to gas from top. These phenomena of Vapex are described by the mathematical model, which is used to calculate live oil production with various values of gas solubility and dispersion. Their optimal values are determined for the vapor extraction of Cold Lake bitumen with butane by matching calculated live oil production with its experimental values published earlier.     

Specimen preparation and image processing and analysis techniques for automated quantification of concrete microcracks and voids

Specimen preparation and image processing/analysis techniques were developed for use in automated quantitative microstructural investigation of concrete, focusing on concrete microcracks and voids. Different specimen preparation techniques were developed for use in fluorescent and scanning electron microscopy (SEM) of concrete; then techniques produce a sharp contrast between microcracks/voids and the body of concrete. The image processing/analysis techniques developed specifically for use with concrete address the following usages: automatic threshold; development of intersecting microcracks/voids and connected voids; distinction of microcracks form voids based on geometric attributes; and noise filtration.