Increasingly in industry, software design and implementation is object-oriented, developed in C++ or Java, and relies heavily on pre-existing software libraries (e.g. the Microsoft Foundation Classes for C++, the Java API for Java). A similar but more tentative trend is developing in high-performance parallel scientific computing. The transition from serial to parallel application development considerably increases the need for library support: task creation and management, data distribution and dynamic redistribution, and inter-process and inter-processor communication and synchronization must be supported. PADRE is a library to support the interoperability of parallel applications. We feel there is significant need for just such a tool to compliment the many domain-specific application frameworks presently available today, but which are generally not interoperable.

We present a new approach for investigating the kinetics of defect migration during annealing of low-energy, ion-implanted silicon, employing a combination of computer simulations and atomic-resolution tunneling microscopy. Using atomically-clean Si(111)-7x7 as a sink for bulk point defects created by 5 keV Xe and Ar irradiation, we observe distinct, temperature-dependent surface arrival rates for vacancies and interstitials. A combination of simulation tools provides a detailed description of the processes that underly the observed temperature-dependence of defect segregation, and the predictions of the simulations agree closely with the experimental observations.

Modern antennas especially arrays are being placed in layers of materials on complex environments. This technique produces aesthetically pleasing structures if necessary, allows for more freedom in structure planning, and can improve antenna performance. In the past, buried antennas have been studied by numerous authors such as in Reference. Recent work on this subject uses spectral and/or numerical moment method formulations. For high frequency analysis it is important to find efficient and accurate methods for design purposes. A rigorous recursive method for plane waves reflection and transmission coefficients by Richmond has been used in the past for dipoles above multilayer slabs. This solution is modified in this paper to account for forward and backward traveling rays with appropriate spread factors for a dipole in the media. Extensive validation for this approximate method shows good agreement with a Method of Moments code. This code is developed at Lawrence Livermore National Laboratory. The geometry for these comparisons uses a dipole in nontruncated dielectric multilayer slabs.

From the history of the Commission, it began with a large number of members. It was found that it was difficult to gain agreement by correspondence between such a large group. A smaller group was elected to operate by correspondence and make decisions. It operated successfully for a half century in this manner. With funding available, the Commission membership grew larger but they.discussed all matters face to face at Commission meetings. Subcommittees were appointed to pursue specialized topics and members reported and discussed their subcommittee results directly to the Commission at the face to face meetings. With the change in the bylaws, future face to face meetings will no longer be an option for the members of the Commission and its subcommittees, unless all members provide their own funds or those of their host institutions. The funding and membership restrictions are all serious topics, which require a thorough discussion.

The probability density function for wave propagating in a straight perfect electrical conductor (PEC) rough wall tunnel is deduced from the mathematical models of the random electromagnetic fields. The field propagating in caves or tunnels is a complex-valued Gaussian random processing by the Central Limit Theorem. The probability density function for single modal field amplitude in such structure is Ricean. Since both expected value and standard deviation of this field depend only on radial position, the probability density function, which gives what is the power distribution, is a radially dependent function. The radio channel places fundamental limitations on the performance of wireless communication systems in tunnels and caves. The transmission path between the transmitter and receiver can vary from a simple direct line of sight to one that is severely obstructed by rough walls and corners. Unlike wired channels that are stationary and predictable, radio channels can be extremely random and difficult to analyze. In fact, modeling the radio channel has historically been one of the more challenging parts of any radio system design; this is often done using statistical methods. In this contribution, we present the most important statistic property, the field probability density function, of wave propagating in …

We present a mathematical model describing the thermodynamic behavior of shape memory alloy wires, as well as a computational technique to solve the resulting system of partial differential equations. The model consists of conservation equations based on a new Helmholtz free energy potential. The computational technique introduces a viscosity-based continuation method, which allows the model to handle dynamic applications where the temporally local behavior of solutions is desired. Computational experiments document that this combination of modeling and solution techniques appropriately predicts the thermally- and stress-induced martensitic phase transitions, as well as the hysteretic behavior and production of latent heat associated with such materials.

We compare inexact Newton and coordinate descent methods for optimizing the quality of a mesh by repositioning the vertices, where quality is measured by the harmonic mean of the mean-ratio metric. The effects of problem size, element size heterogeneity, and various vertex displacement schemes on the performance of these algorithms are assessed for a series of tetrahedral meshes.

The Kinetic Stabilizer (K-S) concept [1] represents a means for stabilizing axisymmetric mirror and tandem-mirror (T-M) magnetic fusion systems against MHD interchange instability modes. Magnetic fusion research has given us examples of axisymmetric mirror confinement devices in which radial transport rates approach the classical ''Spitzer'' level, i.e. situations in which turbulence if present at all, is at too low a level to adversely affect the radial transport [2,3,4]. If such a low-turbulence condition could be achieved in a T-M system it could lead to a fusion power system that would be simpler, smaller, and easier to develop than one based on closed-field confinement, e.g., the tokamak, where the transport is known to be dominated by turbulence. However, since conventional axisymmetric mirror systems suffer from the MHD interchange instability, the key to exploiting this new opportunity is to find a practical way to stabilize this mode. The K-S represents one avenue to achieving this goal. The starting point for the K-S concept is a theoretical analysis by Ryutov [5]. He showed that a MHD-unstable plasma contained in an axisymmetric mirror cell can be MHD-stabilized by the presence of a low-density plasma on the expanding field lines outside the mirrors. If this …

The response of materials to high-dose x-ray exposures needs to be understood for inertial fusion energy (IFE) and inertial confinement fusion applications, where the requirements for IFE are considerably more stringent. In the IFE context, x-ray damage and/or small levels of ablation are of importance for component survivability, generation of debris, and contamination. Ablation quantities of even 1 angstrom per shot would result in material removal of more than 1 cm per year of operation. If even one part in a million of this material made its way to the final optics, it would coat them with a thickness equivalent to several waves of the laser light. Also, small-scale melting and thermomechanical effects, such as fatigue, can result from x-ray heating. These effects potentially become important when multiple shots are considered, and thus, their study requires use of rep-rated experiments. As a part of the High-Average Power Laser Program, the XAPPER experiment has been initiated at Lawrence Livermore National Laboratory. XAPPER produces high doses of low-energy x-rays at repetition rates of up to 10 Hz. Study of x-ray damage is underway. An overview of facility capabilities, results to date, and future plans are provided.

None

We report results from the third and final year of our project (ROA0101-35) to collect seismic event and waveform data recorded in and around the Arabian Peninsula. This effort involves several elements. We are working with King Abdulaziz City for Science and Technology to collect data from the Saudi National Seismic Network, that consists of 38 digital three-component stations (27 broadband and 11 short-period). We have an ongoing collaboration with the Kuwait Institute for Scientific Research, which runs the eight station Kuwait National Seismic Network. We installed two temporary broadband stations in the United Arab Emirates (funded by NNSA NA-24 Office of Non-Proliferation & International Security). In this paper we present a summary of data collected under these efforts including integration of the raw data into LLNL's Seismic Research Database and preliminary analysis of souce parameters and earth structure.

The DUV-FEL facility has been in operation in High Gain Harmonic Generation (HGHG) mode for one year producing 266 nm output from 177 MeV electrons. In this paper we present preliminary results of the Chirped Pulse Amplification (CPA) of HGHG radiation. In the normal HGHG process, a 1 ps electron beam is seeded by chirped 9 ps long 800 nm Ti:Sapphire laser. The electron beam sees only a narrow fraction of the seed laser bandwidth. However, in the CPA case the seed laser pulse length is reduced to 1 ps, and the electron beam sees the full bandwidth. We introduce an energy chirp on electron beam to match the chirp of the seed pulse, enabling the resonant condition for the whole beam. We present measurements of the spectrum bandwidth for various chirp conditions.

The ab initio No-Core Shell Model (NCSM) has recently been expanded to include nucleon-nucleon (NN) and three-nucleon (3N) interactions at the three-body cluster level. Here it is used to predict binding energies and spectra of p-shell nuclei based on realistic NN and 3N interactions. It is shown that 3N force (3NF) properties can be studied in these nuclear systems. First results show that interactions based on chiral perturbation theory lead to a realistic description of {sup 6}Li.

In earlier work, we have described our experiences with the use of decision tree classifiers to identify radio-emitting galaxies with a bent-double morphology in the FIRST astronomical survey. We now extend this work to include ensembles of decision tree classifiers, including two algorithms developed by us. These algorithms randomize the decision at each node of the tree, and because they consider fewer candidate splitting points, are faster than other methods for creating ensembles. The experiments presented in this paper with our astronomy data show that our algorithms are competitive in accuracy, but faster than other ensemble techniques such as Boosting, Bagging, and Arcx4 with different split criteria.

By producing compact representations of hyperspectral image cubes (hypercubes), image storage requirements and the amount of time it takes to extract essential elements of information can both be dramatically reduced. However, these compact representations must preserve the important spectral features within hypercube pixels and the spatial structure associated with background and objects or phenomena of interest. This paper describes a novel approach for automatically and efficiently generating a particular type of compact hypercube representation, referred to as a supercube. The hypercube is segmented into regions that contain pixels with similar spectral shapes that are spatially connected, and the pixel connectivity constraint can be relaxed. Thresholds of similarity in spectral shape between pairs of pixels are derived directly from the hypercube data. One superpixel is generated for each region as some linear combination of pixels belonging to that region. The superpixels are optimal in the sense that the linear combination coefficients are computed so as to minimize the level of noise. Each hypercube pixel is represented in the supercube by applying a gain and bias to the superpixel assigned to the region containing that pixel. Examples are provided.

Fluid flow experiments are being conducted on core specimens of quartz monzonite retrieved from depths of about 1 km at the Desert Peak East EGS site in Churchill County, Nevada. Our immediate goal is to observe permeability evolution in fractures at pressure and temperature conditions appropriate to the Desert Peak geothermal site. Longer term, we aim to evaluate mechanisms that control the evolution of fracture permeability. In the experiments saline water is flowed through an artificial fracture at a constant rate of 0.02 ml/min over a period of several weeks. The constant flow tests are interrupted at selected times for shorter tests in which flow is either stopped or varied between 0 and 2.0 ml/min. The experiments to date were conducted at a confining pressure of 5.5 MPa, pore pressures of 1.38 MPa or 2.07 MPa and temperatures of 167- 169 C. Measurements include differential pressure and electrical resistance across the specimen. The short-term variable flow rate experiments allow us to calculate the effective hydraulic aperture of the fracture at various times during the experiment. Changes in electrical resistivity provide indirect evidence of ongoing mineral dissolution and precipitation processes that are expected to change fracture permeability over time. The early …

We present equation of state results from impulsively stimulated light scattering (ISLS) experiments conducted in diamond anvil cells on pure supercritical fluids. We have made measurements on fluid H{sub 2}O (water), and CH{sub 3}OH (methanol). Sound speeds measured through ISLS have allowed us to refine existing potential models used in the exponential-6 (EXP-6) detonation product library [Fried, L. E., and Howard, W. M., J. Chem. Phys. 109 (17): 7338-7348 (1998).]. The refined models allow us to more accurately assess the chemical composition at the Chapman-Jouget (C-J) state of common energetic materials. We predict that water is present in appreciable quantities at the C-J state of energetic materials HMX, RDX, and nitro methane.

Laboratory measurements are required to establish relationships between the physical properties of unconsolidated sediments and P- and S-wave propagation through them. Previous work has either focused on measurements of compressional wave properties at depths greater than 500 m for oil industry applications or on measurements of dynamic shear properties at pressures corresponding to depths of less than 50 m for geotechnical applications. Therefore, the effects of lithology, fluid saturation, and compaction on impedance and P- and S-wave velocities of shallow soils are largely unknown. We describe two state-of-the-art laboratory experiments. One setup allows us to measure ultrasonic P-wave velocities at very low pressures in unconsolidated sediments (up to 0.1 MPa). The other experiment allows P- and S-wave velocity measurements at low to medium pressures (up to 20 MPa). We summarize the main velocity and attenuation results on sands and sand - clay mixtures under partially saturated and fully saturated conditions in two ranges of pressures (0 - 0.1 MPa and 0.1 - 20 MPa) representative of the top few meters and the top 1 km, respectively. Under hydrostatic pressures of 0.1 to 20 MPa, our measurements demonstrate a P- and S-wave velocity-dependence in dry sands around a fourth root (0.23 …

This paper focuses on the design consideration, fabrication processes and preliminary testing of the stretchable micro-electrode array. We are developing an implantable, stretchable micro-electrode array using polymer-based microfabrication techniques. The device will serve as the interface between an electronic imaging system and the human eye, directly stimulating retinal neurons via thin film conducting traces and electroplated electrodes. The metal features are embedded within a thin ({approx}50 micron) substrate fabricated using poly (dimethylsiloxane) (PDMS), a biocompatible elastomeric material that has very low water permeability. The conformable nature of PDMS is critical for ensuring uniform contact with the curved surface of the retina. To fabricate the device, we developed unique processes for metalizing PDMS to produce robust traces capable of maintaining conductivity when stretched (5%, SD 1.5), and for selectively passivating the conductive elements. An in situ measurement of residual strain in the PDMS during curing reveals a tensile strain of 10%, explaining the stretchable nature of the thin metalized devices.

The role of grain boundary constraint in strain localization and concomitant constitutive response was examined by performing a series of uniaxial compression tests on a tantalum bicrystal. Tantalum single crystals were diffusion bonded to form a (011) 90 twist boundary that was compressed along the common [011] direction. The plastic deformation resulted in the creation of deformation bands away from the highly constraining grain boundary, resembling those bands known from single crystal plastic deformation. Near the grain boundary, such deformation band formation could not be detected. Instead a distinctive pattern of crystal lattice rotation was observed that filled a rather large volume (several millimeters in size) around the bicrystal grain boundary. The internal deformation band structure as well as the crystal lattice rotation pattern near the bicrystal grain boundary were characterized and found to give greater rates of work hardening in the neighborhood of the grain boundary.

The Performance Confirmation program of the Yucca Mountain Repository Development Project needs to employ remotely operated robots to work inside the emplacement drifts which will have an environment unsuitable for humans (radiation environment of up to 200 rad/hour (mostly gamma rays, some neutrons)) and maximum temperatures of 180 C. The robots will be required to operate inside the drifts for up to 8 hours per mission. Based on available functional requirements, we have developed the following specifications for the power needed by the robots:

The Rhode Island Nuclear Science Center (RINSC), located on the Narragansett Bay Campus of the University of Rhode Island, is a state-owned and US NRC-licensed nuclear facility constructed for educational and industrial applications. The main building of RINSC houses a two-megawatt (2 MW) thermal power critical reactor immersed in demineralized water within a shielded tank. As its original design in 1958 by the Rhode Island Atomic Energy Commission focused on the teaching and research use of the facility, only a minimum of 3.85 kg fissile uranium-235 was maintained in the fuel elements to allow the reactor to reach a critical state. In 1986 when RINSC was temporarily shutdown to start US DOE-directed core conversion project for national security reasons, all the U-Al based Highly-Enriched Uranium (HEU, 93% uranium-235 in the total uranium) fuel elements were replaced by the newly developed U{sub 3}Si{sub 2}-Al based Low Enriched Uranium (LEU, {le}20% uranium-235 in the total uranium) elements. The reactor first went critical after the core conversion was achieved in 1993, and feasibility study on the core upgrade to accommodate Boron Neutron-Captured Therapy (BNCT) was completed in 2000 [3]. The 2-MW critical reactor at RINSC which includes six beam tubes, a thermal column, …

Livermore Computing is an early and aggressive adopter of parallel file systems including, for example, GPFS from IBM and Lustre for our present Linux systems. As such, we have acquired more than our share of battle scars from encountering bugs in 'bleeding edge' file systems that we have pressed into production to serve our customers' massive I/O requirements. A major role of the Scalable I/O Project is to detect errors before our end users do. In order to do this, we have developed highly parallel test codes to stress and probe potentially weak areas of file system behavior. This paper describes those test programs and how we make use of them.

In applications where jets propagate through energetic materials, they have been observed to become sufficiently perturbed to reduce their ability to effectively penetrate subsequent material. Analytical calculations of the jet Bernoulli flow provides an estimate of the onset and extent of such perturbations. Although two-dimensional calculations show the back-flow interaction pressure pulses, the symmetry dictates that the flow remains axial. In three dimensions the same pressure impulses can be asymmetrical if the jet is asymmetrical. The 3D calculations thus show parts of the jet having a significant component of radial velocity. On the average the downstream effects of this radial flow can be estimated and calculated by a 2D code by applying a symmetrical radial component to the jet at the appropriate position as the jet propagates through the energetic material. We have calculated the 3D propagation of a radio graphed TOW2 jet with measured variations in straightness and diameter. The resultant three-dimensional perturbations on the jet result in radial flow, which eventually tears apart the coherent jet flow. This calculated jet is compared with jet radiographs after passage through the energetic material for various material thickness and plate thicknesses. We noted that confinement due to a bounding metal plate …