The initial multi-mode interfacial velocity and density perturbations present at the onset of a small Atwood number, incompressible, miscible, Rayleigh-Taylor instability-driven mixing layer have been quantified using a combination of experimental techniques. The streamwise interfacial and spanwise interfacial perturbations were measured using high-resolution thermocouples and planar laser-induced fluorescence (PLIF), respectively. The initial multi-mode streamwise velocity perturbations at the two-fluid density interface were measured using particle-image velocimetry (PIV). It was found that the measured initial conditions describe an initially anisotropic state, in which the perturbations in the streamwise and spanwise directions are independent of one another. The evolution of various fluctuating velocity and density statistics, together with velocity and density variance spectra, were measured using PIV and high-resolution thermocouple data. The evolution of the velocity and density statistics is used to investigate the early-time evolution and the onset of strongly-nonlinear, transitional dynamics within the mixing layer. The early-time evolution of the density and vertical velocity variance spectra indicate that velocity fluctuations are the dominant mechanism driving the instability development. The implications of the present experimental measurements on the initialization of Reynolds-averaged turbulent transport and mixing models and of direct and large-eddy simulations of Rayleigh-Taylor instability-induced turbulence are discussed.

The force fields used in molecular computational biology are not mathematically defined in such a way that their mathematical representation would facilitate the straightforward application of volume visualization techniques. To visualize energy, it is necessary to define a spatial mapping for these fields. Equipped with such a mapping, we can generate volume renderings of the internal energy states in a molecule. We describe our force field, the spatial mapping that we used for energy, and the visualizations that we produced from this mapping. We provide images and animations that offer insight into the computational behavior of the energy optimization algorithms that we employ.

Most physics analysis jobs involve multiple selection steps on the input data. These selection steps are called ''cuts'' or ''queries''. A common strategy to implement these queries is to read all input data from files and then process the queries in memory. In many applications the number of variables used to define these queries is a relative small portion of the overall data set therefore reading all variables into memory takes unnecessarily long time. In this paper we describe an integration effort that can significantly reduce this unnecessary reading by using an efficient compressed bitmap index technology. The primary advantage of this index is that it can process arbitrary combinations of queries very efficiently, while most other indexing technologies suffer from the ''curse of dimensionality'' as the number of queries increases. By integrating this index technology with the ROOT analysis framework, the end-users can benefit from the added efficiency without having to modify their analysis programs. Our performance results show that for multi-dimensional queries, bitmap indices outperform the traditional analysis method up to a factor of 10.

This paper reports on experiments aimed at developing a new high-intensity H{sup -} ion source with long lifetime whose concept had recently been introduced. Starting from the motivation for this effort, several steps of the earlier development work are recapitulated, and the performance of the latest design variant is discussed in detail. The basic concept consists in coupling an ECR ion source to a standard SNS multi-cusp H{sup -} ion source that is driven by pulsed dc, rather than rf, power. As a key result, an electron beam of 1.5 A current has been extracted from the ECR discharge operating at 1.9 kW c. w. power, and a maximum discharge current of 17.5 A was achieved in the H{sup -} ion source. Production of H{sup -} ions, however could not yet been demonstrated in the one, preliminary, experiment conducted so far. The paper concludes by outlining further envisaged development steps for the plasma generator and an expansion towards a novel extraction system.

Steven Impact Tests were performed at low velocity on the explosives TNT, Comp B, and C-4 in attempts to obtain a threshold for reaction. A 76 mm helium driven gas gun was used to accelerate the Steven Test projectiles up to approximately 200 m/s in attempts to react (ignite) the explosive samples. Blast overpressure gauges, acoustic microphones, standard video and high-speed photography were used to characterize the level of any high explosive reaction violence. No bulk reactions were observed in the TNT, Composition B, or C-4 explosive samples impacted up to velocities in the range of 190-200 m/s. This work will outline the experimental details and discuss the lack of reaction when compared to the reaction thresholds of other common explosives.

Measurements made in a laser heated diamond-anvil cell are reported that extend the melting curve of Xe to 80 GPa and 3350 K. The steep lowering of the melting slope (dT/dP) that occurs near 17 GPa and 2750 K results from the hybridization of the p-like valence and d-like conduction states with the formation of clusters in the liquid having Icosahedral Short-Range Order (ISRO).

For the FAIR-project at GSI a model dipole was built at BNL with the nominal field of 4 T and a nominal ramp rate of 1 T/S. The magnet design was similar to the RHIC dipole, with some changes for loss reduction and better cooling. The magnet was already successfully tested in a vertical cryostat, with good training behavior. Cryogenic losses were measured and first results of field harmonics were published. However, for a better understanding of the cooling process, quench currents at several ramp rates were investigated. Detailed measurements of the field harmonics at 2 T/S between 0 and 4 T were performed.

The possibility of Cherenkov-like gluon bremsstrahlung in dense matter is studied. We point out that the occurrence of Cherenkov radiation in dense matter is sensitive to the presence of partonic bound states. This is illustrated by a calculation of the dispersion relation of a massless particle in a simple model in which it couples to two different massive resonance states. We further argue that detailed spectroscopy of jet correlations can directly probe the index of refraction of this matter, which in turn will provide information about the mass scale of these partonic bound states.

A second order accurate embedded boundary method for the two-dimensional wave equation with discontinuous wave propagation speed is described. The wave equation is discretized on a Cartesian grid with constant grid size and the interface (across which the wave speed is discontinuous) is allowed to intersect the mesh in an arbitrary fashion. By using ghost points on either side of the interface, previous embedded boundary techniques for the Neumann and Dirichlet problems are generalized to satisfy the jump conditions across the interface to second order accuracy. The resulting discretization of the jump conditions has the desirable property that each ghost point can be updated independently of all other ghost points, resulting in a fully explicit time-integration method. Numerical examples are given where the method is used to study electro-magnetic scattering of a plane wave by a dielectric cylinder. The numerical solutions are evaluated against the analytical solution due to Mie, and point-wise second order accuracy is confirmed.

Crop yield forecasts provide useful information to a range of users. Yields for several crops in California are currently forecast based on field surveys and farmer interviews, while for many crops official forecasts do not exist. As broad-scale crop yields are largely dependent on weather, measurements from existing meteorological stations have the potential to provide a reliable, timely, and cost-effective means to anticipate crop yields. We developed weather-based models of state-wide yields for 12 major California crops (wine grapes, lettuce, almonds, strawberries, table grapes, hay, oranges, cotton, tomatoes, walnuts, avocados, and pistachios), and tested their accuracy using cross-validation over the 1980-2003 period. Many crops were forecast with high accuracy, as judged by the percent of yield variation explained by the forecast, the number of yields with correctly predicted direction of yield change, or the number of yields with correctly predicted extreme yields. The most successfully modeled crop was almonds, with 81% of yield variance captured by the forecast. Predictions for most crops relied on weather measurements well before harvest time, allowing for lead times that were longer than existing procedures in many cases.

Recent BABAR results on electroweak penguin and leptonic decays are reviewed. In particular, the measurements of B {yields} K{sup (*)}l{sup +}l{sup -} and the preliminary results on B {yields} X{sub s}l{sup +}l{sup -} are presented. Also summarized are the preliminary limits on B{sup +} {yields} l{sup +}{nu} (l = e,{mu}) and B{sup +} {yields} K{sup +}{nu}{bar {nu}}.

Collecting debris samples following a nuclear event requires that operations be conducted from a considerable stand-off distance. An ultra-wide range gamma detector system has been constructed to accomplish both long range radiation search and close range hot sample collection functions. Constructed and tested on a REMOTEC Andros platform, the system has demonstrated reliable operation over six orders of magnitude of gamma dose from 100's of uR/hr to over 100 R/hr. Functional elements include a remotely controlled variable collimator assembly, a NaI(Tl)/photomultiplier tube detector, a proprietary digital radiation instrument, a coaxially mounted video camera, a digital compass, and both local and remote control computers with a user interface designed for long range operations. Long range sensitivity and target location, as well as close range sample selection performance are presented.

The National Ignition Facility (NIF) is a 192 beam Nd-glass laser facility presently under construction at LLNL. When completed, NIF will produce 1.8 MJ, 500 TW of ultraviolet light making it the world's largest and most powerful laser system. NIF will be the world's preeminent facility for performing experiments for Inertial Confinement Fusion (ICF) and High Energy Density Science (HEDS). The Project, begun in 1995, is over 80% complete. The building and the beam path are essentially complete. Nearly all of the functionality of the laser subsystems has been demonstrated. NIF has demonstrated on a single beam basis that it meets its performance goals and shown the laser's precision and flexibility for pulse shaping, pointing, and timing. Beam conditioning techniques, important for target performance, were also demonstrated. The focal spot can be tailored to user specifications using phase plates. Temporal smoothing using smoothing by spectral dispersion (SSD) as well as polarization smoothing was demonstrated. The remaining work is mostly to complete the optics and install them in the beam path and complete the utilities. Presently, eight beams have been activated through the amplifiers and spatial filters to the switchyard wall. Over 150 kJ of 1 {omega} light has been produced …

A new direct constrained optimization algorithm forminimizing the Kohn-Sham (KS) total energy functional is presented inthis paper. The key ingredients of this algorithm involve projecting thetotal energy functional into a sequences of subspaces of small dimensionsand seeking the minimizer of total energy functional within eachsubspace. The minimizer of a subspace energy functional not only providesa search direction along which the KS total energy functional decreasesbut also gives an optimal "step-length" to move along this searchdirection. A numerical example is provided to demonstrate that this newdirect constrained optimization algorithm can be more efficient than theself-consistent field (SCF) iteration.

We present a new technique to extend the embedded-atom method (EAM) for the simulations of non-bulk systems down to the atomic cluster level. To overcome the limitation of the traditional bulk-fit EAM interatomic potentials, bond characteristics from first-principles calculations are systematically included by introducing a local structure dependent prefactor with three additional parameters to the conventional EAM many-body term. The additional parameters improve the local potential landscape virtually for the entire range of atomic configuration space in a quantitative sense. The proposed scheme is applied to two different EAM function sets and validated for both bulk and non-bulk environments in elemental platinum. The obtained material properties, including the binding energies of Pt particles and the Pt adatom diffusion barrier on the Pt (111) surface, show a significant improvement over the conventional EAM formalism.

A diagnostic was developed for the determination of temporal history of an X-ray spot. A pair of thin (0.5 mm) slits image the x-ray spot to a fast scintillator which is coupled to a fast detector, thus sampling a slice of the X-Ray spot. Two other scintillator/detectors are used to determine the position of the spot and total forward dose. The slit signal is normalized to the dose and the resulting signal is analyzed to get the spot size. The position information is used to compensate for small changes due to spot motion and misalignment. The time resolution of the diagnostic is about 1 ns and measures spots from 0.5 mm to over 3 mm. The theory and equations used to calculate spot size and position are presented, as well as data. The calculations assume a symmetric, Gaussian spot. The spot data is generated by the ETA II accelerator, a 2kA, 5.5 MeV, 60 ns electron beam focused on a Tantalum target. The spot generated is typically about 1 mm FWHM. Comparisons are made to an X-ray pinhole camera which images the X-Ray spot (in 2D) at four time slices.

Previous X-ray Diffraction (XRD) and Nuclear Magnetic Resonance (NMR) studies on Ti-doped NaAlH{sub 4} revealed the reaction products of two heavily doped (33.3 at.%) samples that were solvent-mixed and mechanically-milled. This investigation revealed that nano-crystalline or amorphous Al{sub 2}O{sub 3} forms from the possible coordination of aluminum with oxygen atom of the furan ring system from added tetrahydrofuran (THF) in the solvent-mixed sample, and that TiAl{sub 3} forms in mechanically-milled samples. The present paper provides a more sophisticated NMR investigation of the these materials. On heavily doped (33.3 at.%) solvent-mixed samples, {sup 27}Al Magic Angle Spinning (MAS) NMR {sup 27}Al multiple quantum MAS (MQMAS) indicates the presence of an oxide layer of Al{sub 2}O{sub 3} on the surfaces of potentially bulk nanocrystalline Ti, nanocrystalline TiAl{sub 3}, and/or metallic aluminum. The {sup 1}H MAS NMR data also indicate the possible coordination of aluminum with the oxygen atom in the THF. On heavily doped samples that were mechanically milled, {sup 27}Al MAS NMR and static NMR confirms the presence of TiAl{sub 3}. In addition, the {sup 1}H MAS NMR and {sup 1}H spin-lattice relaxation (T{sub 1}) measurements are consistent with the presence of TiH{sub 2}. These results are in agreement with recent …

We have developed a novel insulator concept that involves the use of alternating layers of conductors and insulators with periods less than 1 mm. We have demonstrated that these structures perform 2 to 5 times better than conventional insulators in long pulse, short pulse, and alternating polarity applications. We present new testing results showing exceptional behavior at DC, with gradients in excess of 110kV/cm in vacuum.

We compare several different parallel implementation approaches for the clustering operations performed during adaptive gridding operations in patch-based structured adaptive mesh refinement (SAMR) applications. Specifically, we target the clustering algorithm of Berger and Rigoutsos (BR91), which is commonly used in many SAMR applications. The baseline for comparison is a simplistic parallel extension of the original algorithm that works well for up to O(10{sup 2}) processors. Our goal is a clustering algorithm for machines of up to O(10{sup 5}) processors, such as the 64K-processor IBM BlueGene/Light system. We first present an algorithm that avoids the unneeded communications of the simplistic approach to improve the clustering speed by up to an order of magnitude. We then present a new task-parallel implementation to further reduce communication wait time, adding another order of magnitude of improvement. The new algorithms also exhibit more favorable scaling behavior for our test problems. Performance is evaluated on a number of large scale parallel computer systems, including a 16K-processor BlueGene/Light system.

We study dynamical heterogeneity and growing dynamical length scales in two kinetically constrained models, namely, the one- and two-vacancy assisted triangular lattice gases. One of the models is a strong glassformer and the other is a fragile glassformer. Both exhibit heterogeneous dynamics with broadly distributed timescales as seen in the distribution of persistence times. We show that the Stokes-Einstein relation is violated, to a greater degree in the fragile glassformer, and show how this violation is related to dynamic heterogeneity. Finally, we extract dynamical length scales from structure factors of mobile particles and show, quantitatively, the growth of this lengthscale as density increases. Our results indicate that the one-vacancy assisted lattice gas is in the same universality class as the one-spin facilitated Fredrickson-Andersen model, while the two-vacancy assisted model is in a different universality class. Our findings support the idea that the dynamics of supercooled liquids can be understood from the existence of a zero temperature dynamic critical point.

A moderate (M{approx}5) earthquake struck the northeastern United Arab Emirates (UAE) and northern Oman on March 11, 2002. The event was felt over a wide area of the northern Emirates and was accompanied by smaller (felt) events before and after the March 11 main shock. The event was large enough to be detected and located by global networks at teleseismic distances. We estimated focal mechanism and depth from broadband complete regional waveform modeling. We report a normal mechanism with a slight right-lateral strike-slip component consistent with the large-scale tectonics. The normal component suggests relaxation of obducted crust of the Semail Ophilite (specifically, the Khor Fakkan Block) while the right-lateral strike-slip component of the mechanism is consistent with shear across the Oman Line. Felt earthquakes are rare in the region, however no regional seismic network exists in the UAE to determine local seismicity. This event offers a unique opportunity to study the active tectonics of the region as well as inform future studies of seismic hazard in the UAE and northern Oman.

Residential single family building practice currently ignores the losses of energy and water caused by the poor design of hot water systems. These losses include; the waste of water while waiting for hot water to get to the point of use; the wasted heat as water cools down in the distribution system after a draw; and the energy needed to reheat water that was already heated once before. Average losses of water are estimated to be 6.35 gallons (24.0 L) per day. (This is water that is rundown the drain without being used while waiting for hot water.) The amount of wasted hot water has been calculated to be 10.9 gallons (41.3L) per day. (This is water that was heated, but either is not used or issued after it has cooled off.) A check on the reasonableness of this estimate is made by showing that total residential hot water use averages about 52.6 gallons (199 L) per day. This indicates about 20 percent of average daily hot water is wasted.

None

SciDAC has had a major impact on computational beam dynamics and the design of particle accelerators. Particle accelerators--which account for half of the facilities in the DOE Office of Science Facilities for the Future of Science 20 Year Outlook--are crucial for US scientific, industrial, and economic competitiveness. Thanks to SciDAC, accelerator design calculations that were once thought impossible are now carried routinely, and new challenging and important calculations are within reach. SciDAC accelerator modeling codes are being used to get the most science out of existing facilities, to produce optimal designs for future facilities, and to explore advanced accelerator concepts that may hold the key to qualitatively new ways of accelerating charged particle beams. In this poster we present highlights from the SciDAC Accelerator Science and Technology (AST) project Beam Dynamics focus area in regard to algorithm development, software development, and applications.