This paper describes capabilities, evolution, performance, and applications of the Global Arrays (GA) toolkit. GA was created to provide application programmers with an inteface that allows them to distribute data while maintaining the type of global index space and programming syntax similar to that available when programming on a single processor. The goal of GA is to free the programmer from the low level management of communication and allow them to deal with their problems at the level at which they were originally formulated. At the same time, compatibility of GA with MPI enables the programmer to take advatage of the existing MPI software/libraries when available and appropriate. The variety of applications that have been implemented using Global Arrays attests to the attractiveness of using higher level abstractions to write parallel code.

This paper describes a computerized clavicle identification system primarily designed to resolve the identities of unaccounted-for U.S. soldiers who fought in the Korean War. Elliptical Fourier analysis is used to quantify the clavicle outline shape from skeletons and postero-anterior antemortem chest radiographs to rank individuals in terms of metric distance. Similar to leading fingerprint identification systems, shortlists of the top matching candidates are extracted for subsequent human visual assessment. Two independent tests of the computerized system using 17 field-recovered skeletons and 409 chest radiographs demonstrate that true-positive matches are captured within the top 5% of the sample 75% of the time. These results are outstanding given the eroded state of some field-recovered skeletons and the faintness of the 1950's photofluorographs. These methods enhance the capability to resolve several hundred cold cases for which little circumstantial information exists and current DNA and dental record technologies cannot be applied.

Mesh-based PDE simulation codes are becoming increasingly sophisticated and rely on advanced meshing and discretization tools. Unfortunately, it is still difficult to interchange or interoperate tools developed by different communities to experiment with various technologies or to develop new capabilities. To address these difficulties, we have developed component interfaces designed to support the information flow of mesh-based PDE simulations. We describe this information flow and discuss typical roles and services provided by the geometry, mesh, and field components of the simulation. Based on this delineation for the roles of each component, we give a high-level description of the abstract data model and set of interfaces developed by the Department of Energy's Interoperable Tools for Advanced Petascale Simulation (ITAPS) center. These common interfaces are critical to our interoperability goal, and we give examples of several services based upon these interfaces including mesh adaptation and mesh improvement.

Advances in computational, and imaging techniques have enabled the rapid development of three-dimensional (3-D) models of biological systems in unprecedented detail. Using these advances, 3-D models of the lungs and nasal passages of the rat and human are being developed to ultimately improve predictions of airborne pollutant dosimetry. Techniques for imaging the respiratory tract by magnetic resonance imaging (MRI) were developed to improve the speed and accuracy of geometric data collection for mesh reconstruction. The MRI resolution is comparable to that obtained by manual measurements but at much greater speed and accuracy. Newly developed software (NWGrid) was utilized to translate imaging data from MR into 3-D mesh structures. Together, these approaches significantly reduced the time to develop a 3-D model. This more robust airway structure will ultimately facilitate modeling gas or vapor exchange between the respiratory tract and vasculature as well as enable linkages of dosimetry with cell response models. The 3-D, finite volume, viscoelastic mesh structures form the geometric basis for computational fluid dynamics modeling of inhalation, exhalation and the delivery of individual particles (or concentrations of gas or vapors) to discrete regions of the respiratory tract. The ability of these 3-D models to resolve dosimetry at such a high level of detail will require new techniques to measure regional airflows and particulate deposition for model validation.