Three parallel programming paradigms: Comparisons on an archetypal PDE computation

Three paradigms for distributed-memory parallel computation that free the application programmer from the details of message passing are compared for an archetypal structured scientific computation -- a nonlinear, structured-grid partial differential equation boundary value problem -- using the same algorithm on the same hardware. All of the paradigms -- parallel languages represented by the Portland Group's HPF, (semi-)automated serial-to-parallel source-to-source translation represented by CAP-Tools from the University of Greenwich, and parallel libraries represented by Argonne's PETSc -- are found to be easy to use for this problem class, and all are reasonably effective in exploiting concurrency after a short learning curve. The level of involvement required by the application programmer under any paradigm includes specification of the data partitioning, corresponding to a geometrically simple decomposition of the domain of the PDE. Programming in SPMD style for the PETSc library requires writing only the routines that discretize the PDE and its Jacobian, managing subdomain-to-processor mappings (affine global-to-local index mappings), and interfacing to library solver routines. Programming for HPF requires a complete sequential implementation of the same algorithm as a starting point, introduction of concurrency through subdomain blocking (a task similar to the index mapping), and modest experimentation with rewriting loops to elucidate to the compiler the latent concurrency. Programming with CAPTools involves feeding the same sequential implementation to the CAPTools interactive parallelization system, and guiding the source-to-source code transformation by responding to various queries about quantities knowable only at runtime. Results representative of "the state of the practice" for a scaled sequence of structured grid problems are given on three of the most important contemporary high-performance platforms: the IBM SP, the SGI Origin 2000, and the CRAYY T3E.

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.

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 and volume methods, particularly when the long-time transient behaviour of a system is of interest. The exact solution is obtained at grid collocation points leading to superior accuracy on modest grids. Furthermore, the grid can be freely adapted in time and space to particular flow conditions or geometric variations, especially useful 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 surface. The electromagnetic field determines the instantaneous liquid volume shape, which then affects the electromagnetic field. A general methodology of the pseudo-spectral approach is presented, with several instructive example applications: the aluminium electrolysis MHD problem, induction melting in a cold crucible and the dynamics of AC/DC magnetically levitated droplets. Finally, comparisons with available analytical solutions and to experimental measurements are discussed.

Model-based pilot and data power adaptation in psam with periodic delayed feedback

We consider the optimum design of pilot-symbol-assisted modulation (PSAM) schemes with feedback. The received signal is periodically fed back to the transmitter through a noiseless delayed link and the time-varying channel is modeled as a Gauss-Markov process. We optimize a lower bound on the channel capacity which incorporates the PSAM parameters and Kalman-based channel estimation and prediction. The parameters available for the capacity optimization are the data power adaptation strategy, pilot spacing and pilot power ratio, subject to an average power constraint. Compared to the optimized open-loop PSAM (i.e., the case where no feedback is provided from the receiver), our results show that even in the presence of feedback delay, the optimized power adaptation provides higher information rates at low signal-to-noise ratios (SNR) in medium-rate fading channels. However, in fast fading channels, even the presence of modest feedback delay dissipates the advantages of power adaptation.

Reinforcement corrosion initiation and activation times in concrete structures exposed to severe marine environments

The corrosion of steel reinforcement bars in reinforced concrete structures exposed to severe marine environments usually is attributed to the aggressive nature of chloride ions. In some cases in practice corrosion has been observed to commence already within a few years of exposure even with considerable concrete cover to the reinforcement and apparently high quality concretes. However, there are a number of other cases in practice for which corrosion initiation took much longer, even in cases with quite modest concrete cover and modest concrete quality. Many of these structures show satisfactory long-term structural performance, despite having high levels of localized chloride concentrations at the reinforcement. This disparity was noted already more than 50 years ago, but appears still not fully explained. This paper presents a systematic overview of cases reported in the engineering and corrosion literature and considers possible reasons for these differences. Consistent with observations by others, the data show that concretes made from blast furnace cements have better corrosion durability properties. The data also strongly suggest that concretes made with limestone or non-reactive dolomite aggregates or sufficiently high levels of other forms of calcium carbonates have favourable reinforcement corrosion properties. Both corrosion initiation and the onset of significant damage are delayed. Some possible reasons for this are explored briefly.

Inhibition of prandial and waterspray-induced rat grooming by 8-OH-DPAT

In addition to its hyperphagic effect in rats, 8-OH-DPAT also reduces grooming, but it is uncertain whether the inhibition of grooming is a specific effect or a consequence of response competition from eating. The present experiments explored the effects of 8-OH-DPAT on periprandial grooming and grooming elicited by spraying rats with water. Momentary time sampling over 30 or 60 min, with behaviour scored in one of 6 or 7 (depending on food availability) mutually exclusive categories (feeding, active, scratching, face grooming, body grooming, genital grooming and resting) at 15s intervals, was used for data collection. Non-deprived rats were tested in the presence and absence of food and baseline grooming levels were manipulated by spraying the dorsal surface of the back with water. Data were submitted to ANOVA. The first experiment confirmed that 8-OH-DPAT increased food intake and that this was associated with a parallel increase in feeding observations; active observations were also increased, but resting and total grooming observations were reduced: scratching was reduced even at 0.003mg/kg, face- and body-grooming were reduced at doses > 0.03mg/kg and genital-grooming was least sensitive, only being reduced at 0.1mg/kg. The second experiment revealed that spraying with water had no effect on food intake, feeding or resting observations, but increased total grooming (largely due to increased body-grooming) and reduced activity observations. In rats sprayed with water, 8-OH-DPAT increased food intake (0.1mg/kg) and observations of feeding (0.003 & 0.1mg/kg), but total grooming was dose-dependently inhibited, with genital-grooming most sensitive(> 0.003mg/kg), followed by face-grooming (> 0.01mg/kg) and body-grooming (>0.03mg/kg), whilst low levels of scratching were unaffected. The final experiment tested water-sprayed rats in the absence of food: 8-OH-DPAT increased resting and reduced total grooming, mostly as a consequence of reductions in face- and body-grooming, but there were also modest reductions in scratching. These results confirm that 8-OH-DPAT has a suppressant effect on all aspects of grooming, except where there are probable floor effects, and that this is independent of response competition from increased eating.