This paper describes simulation and experimental results for a 1400A, {+-} 900V peak rated, switch mode power supply for SNS Injection Kicker Magnets. For each magnet (13 m{Omega}, 160{micro}H), the power supply must supply controlled pulses at 60 Hz repetition rate. The pulse current must rise from zero to maximum in less than 1 millisec in a controlled manner, flat top for up to 2 millisec, and should fall in a controlled manner to less than 4A within 500{micro}s. The low current performance during fall time is the biggest challenge in this power supply. The simulation results show that to meet the controlled fall of the current and the current ripple requirements, voltage loop bandwidth of at least 10 kHz and switching frequency of at least 100 kHz are required. To achieve high power high frequency switching with IGBT switches, a series connected topology with three phase shifted (O{sup o}, 60{sup o} & 120{sup o}) converters each with 40 kHz switching frequency (IGBT at 20kHz), has been achieved. In this paper, the circuit topology, relevant system specifications and experimental results that meet the requirements of the power supply are described in detail. A unique six pulse SCR rectifier circuit with …

Version 3 of the Hadley Centre Atmospheric Model (HadAM3) has been used to demonstrate one means of comparing a general circulation model with observations for a specific climate perturbation, namely the strong 1997/98 El Nino. This event was characterized by the collapse of the tropical Pacific's Walker circulation, caused by the lack of a zonal sea surface temperature gradient during the El Nino. Relative to normal years, cloud altitudes were lower in the western portion of the Pacific and higher in the eastern portion. HadAM3 likewise produced the observed collapse of the Walker circulation, and it did a reasonable job of reproducing the west/east cloud structure changes. This illustrates that the 1997/98 El Nino serves as a useful means of testing cloud-climate interactions in climate models.

The systematic requirements to divert an object on an earth-impacting course are developed relating the minimum velocity perturbation (both magnitude and direction) to the time available before impact. This, coupled with the accuracy to which orbits can be determined, restricts the time available for any mitigation technology to operate. Because nuclear energy densities are nearly a million times higher than those possible with chemical bonds, it is the most mass efficient means for storing delivering energy with today's technology. The question is how to most effectively apply that energy. This paper will examine the simple case of shattering the body, as well as a more controlled approach in which one or more small velocity increments divert a body. The optimal approach depends on the detailed circumstances, but in either case, already developed technology permits a successful diversion with a few years to decades of notice. The success of nuclear options on relatively short timescales permits consideration of other technologies that while not so well developed might be sufficiently improved to divert small (100 meter) bodies.

Many turbulent premixed flames of practical interest are statistically stationary. They occur in combustors that have anchoring mechanisms to prevent blow-off and flashback. The stabilization devices often introduce a level of geometric complexity that is prohibitive for detailed computational studies of turbulent flame dynamics. As a result, typical detailed simulations are performed in simplified model configurations such as decaying isotropic turbulence or inflowing turbulence. In these configurations, the turbulence seen by the flame either decays or, in the latter case, increases as the flame accelerates toward the turbulent inflow. This limits the duration of the eddy evolutions experienced by the flame at a given level of turbulent intensity, so that statistically valid observations cannot be made. In this paper, we apply a feedback control to computationally stabilize an otherwise unstable turbulent premixed flame in two dimensions. For the simulations, we specify turbulent in flow conditions and dynamically adjust the integrated fueling rate to control the mean location of the flame in the domain. We outline the numerical procedure, and illustrate the behavior of the control algorithm. We use the simulations to study the propagation and the local chemical variability of turbulent flame chemistry.

We present a proposal for lowering the overhead of interface contract checking for science and engineering applications. Run-time enforcement of assertions is a well-known technique for improving the quality of software; however, the performance penalty is often too high for their retention during deployment, especially for long-running applications that depend upon iterative operations. With an efficient adaptive approach the benefits of run-time checking can continue to accrue with minimal overhead. Examples from scientific software interfaces being developed in the high performance computing research community will be used to measure the efficiency and effectiveness of this approach.

The usual approach to deal with noise present in many real-world optimization problems is to take an arbitrary number of samples of the objective function and use the sample average as an estimate of the true objective value. The number of samples is typically chosen arbitrarily and remains constant for the entire optimization process. This paper studies an adaptive sampling technique that varies the number of samples based on the uncertainty of deciding between two individuals. Experiments demonstrate the effect of adaptive sampling on the final solution quality reached by a genetic algorithm and the computational cost required to find the solution. The results suggest that the adaptive technique can effectively eliminate the need to set the sample size a priori, but in many cases it requires high computational costs.

Residual radioactive waste was removed from Tank 18F in the F-Area Tank Farm at Savannah River Site (SRS), using the advanced design mixer pump (ADMP). Known as a slurry pump, the ADMP is a 55 foot long pump with an upper motor mounted to a steel super structure, which spans the top of the waste tank. The motor is connected by a long vertical drive shaft to a centrifugal pump, which is submerged in waste near the tank bottom. The pump mixes, or slurries, the waste within the tank so that it may be transferred out of the tank. Tank 18F is a 1.3 million gallon, 85 foot diameter underground waste storage tank, which has no internal components such as cooling coils or structural supports. The tank contained a residual 47,000 gallons of nuclear waste, consisting of a gelatinous radioactive waste known as sludge and particulate zeolite. The prediction of the ADMP success was based on nearly twenty five years of research and the application of that research to slurry pump technology. Many personnel at SRS and Pacific Northwest National Laboratories (PNNL) have significantly contributed to these efforts. This report summarizes that research which is pertinent to the ADMP performance …

An analysis of the standard system of differential equations describing multi-speed flows of multi-phase media is performed. It is proved that the Cauchy problem, as posed in most best-estimate thermal-hydraulic codes, results in unstable solutions and potentially unreliable description of many physical phenomena. A system of equations, free from instability effects, is developed allowing more rigorous numerical modeling.

We present the first measurement of directed flow (v{sub 1}) at the Relativistic Heavy Ion Collider (RHIC). v{sub 1} is found to be consistent with zero at pseudorapidities {eta} from -1.2 to 1.2, then rises to the level of a couple of percent over the range 2.4 < |{eta}| < 4. The latter observation is similar to that from NA49 if the SPS rapidities are shifted by the difference in beam rapidity between RHIC and SPS. We studied the evolution of elliptic flow from p + p collisions through d + Au collision, and onto Au + Au collisions. Measurements of higher harmonics are presented and discussed.

The addition of the two Forward TPCs to the STAR detector allows one to measure anisotropic flow at forward pseudorapidities. This made possible the first measurement of directed flow at collision energies of {radical}s{sub NN} = 200 GeV. PHOBOS' results on elliptic flow at forward rapidities were confirmed, and the sign of v{sub 2} was determined to be positive for the first time at RHIC energies. The higher harmonic, v{sub 4}, is consistent with the recently suggested v{sub 2}2 scaling behavior.

Optimizing compilers have a long history of applying loop transformations to C and Fortran scientific applications. However, such optimizations are rare in compilers for object-oriented languages such as C++ or Java, where loops operating on user-defined types are left unoptimized due to their unknown semantics. Our goal is to reduce the performance penalty of using high-level object-oriented abstractions. We propose an approach that allows the explicit communication between programmers and compilers. We have extended the traditional Fortran loop optimizations with an open interface. Through this interface, we have developed techniques to automatically recognize and optimize user-defined array abstractions. In addition, we have developed an adapted constant-propagation algorithm to automatically propagate properties of abstractions. We have implemented these techniques in a C++ source-to-source translator and have applied them to optimize several kernels written using an array-class library. Our experimental results show that using our approach, applications using high-level abstractions can achieve comparable, and in cases superior, performance to that achieved by efficient low-level hand-written codes.

Mercury emissions from coal fired plants will be limited by regulations enforced by the Environmental Protection Agency. However, there is still debate over whether the limits should be on a plant specific basis or a nationwide basis. The nationwide basis allows a Cap and Trade program similar to that for other air pollutants. Therefore, a major issue is the magnitude and extent of local deposition. Computer modeling suggests that increased local deposition will occur on a local (2 to 10 Km) to regional scale (20 to 50 Km) with the increase being a small percentage of background deposition on the regional scale. The amount of deposition depends upon many factors including emission rate, chemical form of mercury emitted (with reactive gaseous mercury depositing more readily than elemental mercury), other emission characteristics (stack height, exhaust temperature, etc), and meteorological conditions. Modeling suggests that wet deposition will lead to the highest deposition rates and that these will occur locally. Dry deposition is also predicted to deposit approximately the same amount of mass as wet deposition, but over a much greater area. Therefore, dry deposition rates will contribute a fraction of total deposition on the regional scale. The models have a number of …

The rationale for regulating air emissions of mercury from U.S. coal-fired power plants largely depends on mathematical dispersion modeling, including the atmospheric chemistry processes that affect the partitioning of Hg emissions into elemental (Hg{sub 0}) and the reactive (RGM) forms that may deposit more rapidly near sources. This paper considers and evaluates the empirical support for this paradigm. We consider the extant experimental data at three spatial scales: local (< 30 km), regional (< {approx}300 km), and national (multi-state data). An additional issue involves the finding of excess Hg levels in urban areas.

Sponsored by the US Nuclear Regulatory Commission (NRC), Brookhaven National Laboratory (BNL) is carrying out a research program to develop a technical basis to support the safety evaluation of deeply embedded and/or buried (DEB) structures as proposed for advanced reactor designs. In this program, the methods and computer programs established for the assessment of soil-structure interaction (SSI) effects for the current generation of light water reactors are evaluated to determine their applicability and adequacy in capturing the seismic behavior of DEB structures. This paper presents an assessment of the simplified vs. detailed methodologies for seismic analyses of DEB structures. In this assessment, a lump-mass beam model is used for the simplified approach and a finite element representation is employed for the detailed method. A typical containment structure embedded in a soil profile representative of a typical nuclear power plant site was utilized, considering various embedment depths from shallow to full burial. BNL used the CARES program for the simplified model and the SASSI2000 program for the detailed analyses. The calculated response spectra at the key locations of the DEB structure are used for the performance assessment of the applied methods for different depths of burial. Included in the paper are: …

New or modified models have been proposed for the much-studied Si(111)-(2x1) surface structure, including: a reverse-tilted p-bonded chain model (by Zitzlsperger et al); a three-bond scission model (by Haneman et al); and a p-bonded chain model with enhanced vibrations (present work). These models are compared here to the generally accepted modified p-bonded chain model (by Himpsel et al, 1984), by analyzing low-energy electron diffraction (LEED) I-V curves measured earlier. Using the efficient automated tensor LEED technique, the models can be refined to a much greater degree than with earlier methods of LEED analysis. This study distinctly favors the earlier modified p-bonded chain model, but with strongly enhanced vibrations. To compare models that have different numbers of adjustable free parameters a Hamilton ratio test is used: it can distinguish between improvement due to a better model and improvement due only to more parameters.

We calculate the atmospheric structure of an accretion disk around a Kerr black hole and obtain its X-ray spectrum, which exhibits prominent atomic transitions under certain circumstances. The gravitational and Doppler (red)shifts of the C V, C VI, O VII, O VIII, and Fe I-XXVI emission lines are observable in active galaxies. We quantify the line emissivities as a function of radius, to identify the effects of atmospheric structure, and to determine the usefulness of these lines for probing the disk energetics. The line emissivities do not always scale linearly with the incident radiative energy, as in the case of Fe XXV and Fe XXVI. Our model incorporates photoionization and thermal balance for the plasma, the hydrostatic approximation perpendicular to the plane of the disk, and general relativistic tidal forces. We include radiative recombination rates, fluorescence yields, Compton scattering, and photoelectric opacities for the most abundant elements.

We have developed interatomic potentials for Pt-Ni and Pt-Re alloys within the modified embedded atom method (MEAM). Furthermore, we applied these potentials to study the equilibrium structures of Pt75Ni25 and Pt75Re25 nanoparticles at T=600 K using the Monte Carlo method. In this work, the nanoparticles are assumed to have disordered fcc cubo-octahedral shapes (terminated by {l_brace}111{r_brace} and {l_brace}100{r_brace} facets) and contain from 586 to 4033 atoms (corresponding to a diameter from 2.5 to 5 nm). It was found that, due to surface segregation, (1) the Pt75Ni25 nanoparticles form a surface-sandwich structure: the Pt atoms are enriched in the outermost and third atomic shells, while the Ni atoms are enriched in the second atomic shell; (2) the equilibrium Pt75Re25 nanoparticles adopt a core-shell structure: a Pt-enriched shell surrounding a Pt-deficient core.

President Eisenhower's hopes for nuclear technology still resonate, but the challenges to fulfilling them are much different today. On December 8, 1953, President Eisenhower, returning from his meeting with the leaders of Britain and France at the Bermuda Summit, flew directly to New York to address the United Nations General Assembly. His presentation, known afterwards as the ''Atoms for Peace'' speech, was bold, broad, and visionary. Eisenhower highlighted dangers associated with the further spread of nuclear weapons and the end of the thermonuclear monopoly, but the president also pointed to opportunities. Earlier that year, Stalin had died and the Korean War armistice was signed. Talks on reunification of Austria were about to begin. The speech sought East-West engagement and outlined a framework for reducing nuclear threats to security while enhancing the civilian benefits of nuclear technology. One specific proposal offered to place surplus military fissile material under the control of an ''international atomic energy agency'' to be used for peaceful purposes, especially economic development. Eisenhower clearly recognized the complex interrelationships between different nuclear technologies and the risks and the benefits that accrue from each. The widespread use of civilian nuclear technology and absence of any use of a nuclear weapon …

The Material Control and Accountability (MC and A) Group has developed an electronic scan system to acknowledge issuance and returns of Tamper Indicating Device (TID) seals. Important MC and A features of the system are: 1. The system requires the issuer identification and the identification of applicators to be entered, thereby ensuring that at least two qualified issuers and applicators possess the seals. Also, Operations is prompted to ensure the Two Person Rule is met during issuance and application of the seals. 2. All input is date and time stamped to ease resolution of anomalies. 3. The system requires all information to be input before allowing the user to logoff, thereby eliminating the problem of incomplete information in the records. 4. The input is immediately available to FBL MC and A personnel outside the facility who then know to expect completed procedures and forms regarding these transactions. The paper will describe the application of these features in routine operations as well as the development effort and final configuration of the system.

QR and LU factorizations for dense matrices are important linear algebra computations that are widely used in scientific applications. To efficiently perform these computations on modern computers, the factorization algorithms need to be blocked when operating on large matrices to effectively exploit the deep cache hierarchy prevalent in today's computer memory systems. Because both QR (based on Householder transformations) and LU factorization algorithms contain complex loop structures, few compilers can fully automate the blocking of these algorithms. Though linear algebra libraries such as LAPACK provides manually blocked implementations of these algorithms, by automatically generating blocked versions of the computations, more benefit can be gained such as automatic adaptation of different blocking strategies. This paper demonstrates how to apply an aggressive loop transformation technique, dependence hoisting, to produce efficient blockings for both QR and LU with partial pivoting. We present different blocking strategies that can be generated by our optimizer and compare the performance of auto-blocked versions with manually tuned versions in LAPACK, both using reference BLAS, ATLAS BLAS and native BLAS specially tuned for the underlying machine architectures.

This paper discusses an innovative application of new Markov localization techniques that combat the problem of odometry drift, allowing a novel control architecture developed at the Idaho National Engineering and Environmental Laboratory (INEEL) to be utilized within a sensor characterization facility developed at the Remote Sensing Laboratory (RSL) in Nevada. The new robotic capability provided by the INEEL will allow RSL to test and evaluate a wide variety of sensors including radiation detection systems, machine vision systems, and sensors that can detect and track heat sources (e.g. human bodies, machines, chemical plumes). By accurately moving a target at varying speeds along designated paths, the robotic solution allows the detection abilities of a wide variety of sensors to be recorded and analyzed.

We report the first observations of the fourth harmonic (v{sub 4}) in the azimuthal distribution of particles at RHIC. The measurement was done taking advantage of the large elliptic flow generated at RHIC. The integrated v{sub 4} is about a factor of 10 smaller than v{sub 2}. For the sixth (v{sub 6}) and eighth (v{sub 8}) harmonics upper limits on the magnitudes are reported.

We present the B physics results from the CDF and D0 experiments at the Tevatron Run II at Fermilab and their future prospect. This includes various B mass and lifetime measurements, B mixing, the confirmation of the discovery of the X particle, rare decays, CP violation, and spectroscopy.

The Spallation Neutron Source accelerator systems will provide a 1 GeV, 1.44 MW proton beam to a liquid mercury target for neutron production. Beam tuning dumps are provided at the end of the linac (the Linac Dump) and in the Ring-to-Target transport line (the Extraction Dump) [1]. Thin windows are required to separate the accelerator vacuum from the poor vacuum upstream of the beam dump. There are several challenging engineering issues that have been addressed in the window design. Namely, handling of the high local power density deposited by the stripped electrons from the H-beam accelerated in the linac, and the need for low-exposure removal and replacement of an activated window. The thermal design of the linac dump window is presented, as is the design of a vacuum clamp and mechanism that allows remote removal and replacement of the window.