Parallel implementation of the recurrence method for computing the power-spectral density of thin avalanche photodiodes

A simulation program has been developed to calculate the power-spectral density of thin avalanche photodiodes, which are used in optical networks. The program extends the time-domain analysis of the dead-space multiplication model to compute the autocorrelation function of the APD impulse response. However, the computation requires a large amount of memory space and is very time consuming. We describe our experiences in parallelizing the code using both MPI and OpenMP. Several array partitioning schemes and scheduling policies are implemented and tested Our results show that the OpenMP code is scalable up to 64 processors on an SGI Origin 2000 machine and has small average errors.

Parallel bandwidth characteristics calculations for thin avalanche photodiodes on a SGI Origin 2000 supercomputer

An important factor for high-speed optical communication is the availability of ultrafast and low-noise photodetectors. Among the semiconductor photodetectors that are commonly used in today’s long-haul and metro-area fiber-optic systems, avalanche photodiodes (APDs) are often preferred over p-i-n photodiodes due to their internal gain, which significantly improves the receiver sensitivity and alleviates the need for optical pre-amplification. Unfortunately, the random nature of the very process of carrier impact ionization, which generates the gain, is inherently noisy and results in fluctuations not only in the gain but also in the time response. Recently, a theory characterizing the autocorrelation function of APDs has been developed by us which incorporates the dead-space effect, an effect that is very significant in thin, high-performance APDs. The research extends the time-domain analysis of the dead-space multiplication model to compute the autocorrelation function of the APD impulse response. However, the computation requires a large amount of memory space and is very time consuming. In this research, we describe our experiences in parallelizing the code in MPI and OpenMP using CAPTools. Several array partitioning schemes and scheduling policies are implemented and tested. Our results show that the code is scalable up to 64 processors on a SGI Origin 2000 machine and has small average errors.

The ParaWise expert assistant - widening accessibility to efficient and scalable tool Generated OpenMP code

Despite the apparent simplicity of the OpenMP directive shared memory programming model and the sophisticated dependence analysis and code generation capabilities of the ParaWise/CAPO tools, experience shows that a level of expertise is required to produce efficient parallel code. In a real world application the investigation of a single loop in a generated parallel code can soon become an in-depth inspection of numerous dependencies in many routines. The additional understanding of dependencies is also needed to effectively interpret the information provided and supply the required feedback. The ParaWise Expert Assistant has been developed to automate this investigation and present questions to the user about, and in the context of, their application code. In this paper, we demonstrate that knowledge of dependence information and OpenMP are no longer essential to produce efficient parallel code with the Expert Assistant. It is hoped that this will enable a far wider audience to use the tools and subsequently, exploit the benefits of large parallel systems.

Tree diameter distribution modelling: introducing the logitlogistic distribution

Johnson's SB distribution is a four-parameter distribution that is transformed into a normal distribution by a logit transformation. By replacing the normal distribution of Johnson's SB with the logistic distribution, we obtain a new distributional model that approximates SB. It is analytically tractable, and we name it the "logitlogistic" (LL) distribution. A generalized four-parameter Weibull model and the Burr XII model are also introduced for comparison purposes. Using the distribution "shape plane" (with axes skew and kurtosis) we compare the "coverage" properties of the LL, the generalized Weibull, and the Burr XII with Johnson's SB, the beta, and the three-parameter Weibull, the main distributions used in forest modelling. The LL is found to have the largest range of shapes. An empirical case study of the distributional models is conducted on 107 sample plots of Chinese fir. The LL performs best among the four-parameter models.

Comparison between weight loss of bends in a pneumatic conveyor and erosion rate obtained in a centrifugal erosion tester for the same materials

In order to find a link between results obtained from a laboratory erosion tester and tests carried out on a pneumatic conveyor, a comparison has been made between weight loss from bends on an industrial-scale pneumatic conveyor and erosion rates obtained in a small centrifugal erosion tester, for the same materials. Identical test conditions have been applied to both experiments so that comparable test results have been obtained. The erosion rate of mild steel commonly used as the wall material of conveyor pipes and pipe bends was determined individually on both test rigs. A relationship between weight loss from the bends and erosion rate determined from the tester has been developed. A discussion based on the results and their applicability to the prediction of wear in pneumatic conveyors concludes the paper.

Effect of particle concentration on erosion rate of mild steel bends in a pneumatic conveyor

Particle concentration is known as a main factor that affects erosion rate of pipe bends in pneumatic conveyors. With consideration of different bend radii, the effect of particle concentration on weight loss of mild steel bends has been investigated in an industrial scale test rig. Experimental results show that there was a significant reduction of the specific erosion rate for high particle concentrations. This reduction was considered to be as a result of the shielding effect during the particle impacts. An empirical model is given. Also a theoretical study of scaling on the shielding effect, and comparisons with some existing models, are presented. It is found that the reduction in specific erosion rate (relative to particle concentration) has a stronger relationship in conveying pipelines than has been found in the erosion tester.

Effect of bend orientation on life and puncture point location due to solid particle erosion of a high concentration flow in pneumatic conveyors

An investigation into predicting failure of pneumatic conveyor pipe bends due to hard solid particle impact erosion has been carried out on an industrial scale test rig. The bend puncture point locations may vary with many factors. However, bend orientation was suspected of being a main factor due to the biased particle distribution pattern of a high concentration flow. In this paper, puncture point locations have been studied with different pipe bend orientations and geometry (a solids loading ratio of 10 being used for the high concentration flow). Test results confirmed that the puncture point location is indeed most significantly influenced by the bend orientation (especially for a high concentration flow) due to the biased particle distribution and biased particle flux distribution.

Assessing the scalability of multiphysics tools for modeling solidification and melting processes on parallel clusters

A comprehensive solution of solidification/melting processes requires the simultaneous representation of free surface fluid flow, heat transfer, phase change, nonlinear solid mechanics and, possibly, electromagnetics together with their interactions, in what is now known as multiphysics simulation. Such simulations are computationally intensive and the implementation of solution strategies for multiphysics calculations must embed their effective parallelization. For some years, together with our collaborators, we have been involved in the development of numerical software tools for multiphysics modeling on parallel cluster systems. This research has involved a combination of algorithmic procedures, parallel strategies and tools, plus the design of a computational modeling software environment and its deployment in a range of real world applications. One output from this research is the three-dimensional parallel multiphysics code, PHYSICA. In this paper we report on an assessment of its parallel scalability on a range of increasingly complex models drawn from actual industrial problems, on three contemporary parallel cluster systems.

Spatial variability of air pollution in the vicinity of a permanent monitoring station in central Paris

The strong spatial and temporal variability of traffic-related air pollution detected at roadside locations in a number of European cities has raised the question of how representative the site and time period of air quality measurements actually can be. To address this question, a 7-month sampling campaign was carried out on a major road axis (Avenue Leclerc) leading to a very busy intersection (Place Basch) in central Paris, covering the surroundings of a permanent air quality monitoring station. This station has recorded the highest CO and NOx concentrations during recent years in the region of Paris. Diffusive BTX samplers as well as a mobile monitoring unit equipped with real-time CO, NOx and O3 analysers and meteorological instruments were used to reveal the small-scale pollution gradients and their temporal trends near the permanent monitoring station. The diffusive measurements provided 7-day averages of benzene, toluene, xylene and other hydrocarbons at different heights above the ground and distances from the kerb covering summer and winter periods. Relevant traffic and meteorological data were also obtained on an hourly basis. Furthermore, three semiempirical dispersion models (STREET-SRI, OSPM and AEOLIUS) were tested for an asymmetric canyon location in Av. Leclerc. The analysis of this comprehensive data set has helped to assess the representativeness of air quality monitoring information.

Pseudo-spectral solutions for fluid flow and heat transfer in electro-metallurgical applications

The pseudo-spectral solution method offers a flexible and fast alternative to the more usual finite element/volume/difference methods, particularly when the long-time transient behaviour of a system is of interest. Since the exact solution is obtained at the grid collocation points superior accuracy can be achieved on modest grid resolution. Furthermore, the grid can be freely adapted with time and in space, to particular flow conditions or geometric variations. This is especially advantageous where strongly coupled, time-dependent, multi-physics solutions are investigated. Examples include metallurgical applications involving the interaction of electromagnetic fields and conducting liquids with a free sutface. The electromagnetic field then determines the instantaneous liquid volume shape and the liquid shape affects in turn the electromagnetic field. In AC applications a thin "skin effect" region results on the free surface that dominates grid requirements. Infinitesimally thin boundary cells can be introduced using Chebyshev polynomial expansions without detriment to the numerical accuracy. This paper presents a general methodology of the pseudo-spectral approach and outlines the solution procedures used. Several instructive example applications are given: the aluminium electrolysis MHD problem, induction melting and stirring and the dynamics of magnetically levitated droplets in AC and DC fields. Comparisons to available analytical solutions and to experimental measurements will be discussed.

Nonlinear MHD stability of aluminium reduction cells

An industrial electrolysis cell used to produce primary aluminium is sensitive to waves at the interface of liquid aluminium and electrolyte. The interface waves are similar to stratified sea layers [1], but the penetrating electric current and the associated magnetic field are intricately involved in the oscillation process, and the observed wave frequencies are shifted from the purely hydrodynamic ones [2]. The interface stability problem is of great practical importance because the electrolytic aluminium production is a major electrical energy consumer, and it is related to environmental pollution rate. The stability analysis was started in [3] and a short summary of the main developments is given in [2]. Important aspects of the multiple mode interaction have been introduced in [4], and a widely used linear friction law first applied in [5]. In [6] a systematic perturbation expansion is developed for the fluid dynamics and electric current problems permitting reduction of the three-dimensional problem to a two dimensional one. The procedure is more generally known as “shallow water approximation” which can be extended for the case of weakly non-linear and dispersive waves. The Boussinesq formulation permits to generalise the problem for non-unidirectionally propagating waves accounting for side walls and for a two fluid layer interface [1]. Attempts to extend the electrolytic cell wave modelling to the weakly nonlinear case have started in [7] where the basic equations are derived, including the nonlinearity and linear dispersion terms. An alternative approach for the nonlinear numerical simulation for an electrolysis cell wave evolution is attempted in [8 and references there], yet, omitting the dispersion terms and without a proper account for the dissipation, the model can predict unstable waves growth only. The present paper contains a generalisation of the previous non linear wave equations [7] by accounting for the turbulent horizontal circulation flows in the two fluid layers. The inclusion of the turbulence model is essential in order to explain the small amplitude self-sustained oscillations of the liquid metal surface observed in real cells, known as “MHD noise”. The fluid dynamic model is coupled to the extended electromagnetic simulation including not only the fluid layers, but the whole bus bar circuit and the ferromagnetic effects [9].

Maximising heat transfer efficiency in the cold crucible induction melting process

Induction heating is an efficient method used to melt electrically conductive materials, particularly if melting takes place in a ceramic crucible. This form of melting is particularly good for alloys, as electromagnetic forces set up by the induction coil lead to vigorous stirring of the melt ensuring homogeneity and uniformity in temperature. However, for certain reactive alloys, or where high purity is required, ceramic crucibles cannot be used, but a water-cooled segmented copper crucible is employed instead. Water cooling prevents meltdown or distortion of the metal wall, but much of the energy goes into the coolant. To reduce this loss, the electromagnetic force generated by the coil is used to push the melt away from the walls and so minimise contact with water-cooled surfaces. Even then, heat is lost through the crucible base where contact is inevitable. In a collaborative programme between Greenwich and Birmingham Universities, computer modelling has been used in conjunction with experiments to improve the superheat attainable in the melt for a,number of alloys, especially for y-TiAl intermetallics to cast aeroengine turbine blades. The model solves the discretised form of the turbulent Navier-Stokes, thermal energy conservation and Maxwell equations using a Spectral Collocation technique. The time-varying melt envelope is followed explicitly during the computation using an adaptive mesh. This paper briefly describes the mathematical model used to represent the interaction between the magnetic field, fluid flow, heat transfer and change of phase in the crucible and identifies the proportions of energy used in the melt, lost in the crucible base and in the crucible walls. The role of turbulence is highlighted as important in controlling heat losses and turbulence damping is introduced as a means of improving superheat. Model validation is against experimental results and shows good agreement with measured temperatures and energy losses in the cooling fluid throughout the melting cycle.

Approximation Algorithms for a Symmetric Quadratic Knapsack Problem and its Scheduling Applications

We consider a knapsack problem to minimize a symmetric quadratic function. We demonstrate that this symmetric quadratic knapsack problem is relevant to two problems of single machine scheduling: the problem of minimizing the weighted sum of the completion times with a single machine non-availability interval under the non-resumable scenario; and the problem of minimizing the total weighted earliness and tardiness with respect to a common small due date. We develop a polynomial-time approximation algorithm that delivers a constant worst-case performance ratio for a special form of the symmetric quadratic knapsack problem. We adapt that algorithm to our scheduling problems and achieve a better performance. For the problems under consideration no fixed-ratio approximation algorithms have been previously known.

An Analysis of Human Behavior during Evacuation

The World Trade Center Evacuation: The evacuation of the WTC complex represents one of the largest full-scale evacuations of people in modern times.