BI-OMP supports DOE's mission in Climate Change Research. The program provides the fundamental understanding of the linkages between carbon and nitrogen cycles in ocean margins. Researchers are providing a mechanistic understanding of these cycles, using the tools of modern molecular biology. The models that will allow policy makers to determine safe levels of greenhouse gases for the Earth System.

The existing superconducting linac ATLAS is in many respects an ideal secondary beam accelerator for an ISOL (Isotope separator on-line) type radioactive beam facility. Such a facility would require the addition of two major accelerator elements: a low charge state injector for the existing heavy ion linac, and a primary beam accelerator providing 220 MV of acceleration for protons and light ions. Development work for both of these elements, including the option of superconducting cavities for the primary beam accelerator is discussed.

Using surface micromachined samples, we demonstrate the accurate measurement of cantilever beam adhesion by using test structures which are adhered over long attachment lengths. We show that this configuration has a deep energy well, such that a fracture equilibrium is easily reached. When compared to the commonly used method of determining the shortest attached beam, the present method is much less sensitive to variations in surface topography or to details of capillary drying.

Adaptive mesh refinement is complicated by both the algorithms and the dynamic nature of the computations. In parallel the complexity of getting good performance is dependent upon the architecture and the application. Most attempts to address the complexity of AMR have lead to the development of library solutions, most have developed object-oriented libraries or frameworks. All attempts to date have made numerous and sometimes conflicting assumptions which make the evaluation of performance of AMR across different applications and architectures difficult or impracticable. The evaluation of different approaches can alternatively be accomplished through simulation of the different AMR processes. In this paper we outline our research work to simulate the processing of adaptive mesh refinement grids using a distributed array class library (P++). This paper presents a combined analytic and empirical approach, since details of the algorithms can be readily predicted (separated into specific phases), while the performance associated with the dynamic behavior must be studied empirically. The result, AMRSim, provides a simple way to develop bounds on the expected performance of AMR calculations subject to constraints given by the algorithms, frameworks, and architecture.

Simulation of nonisothermal, multiphase flow through fractured geothermal reservoirs involves the solution of a system of strongly nonlinear algebraic equations. The Newton-Raphson method used to solve such a nonlinear system of equations requires the evaluation of a Jacobian matrix. In this paper we discuss automatic differentiation (AD) as a method for analytically computing the Jacobian matrix of derivatives. Robustness and efficiency of the AD-generated derivative codes are compared with a conventional derivative computation approach based on first-order finite differences.

An alternative methodology is described for Large-Eddy Simulation of flows involving shocks, turbulence and mixing. In lieu of filtering the governing equations, it is postulated that the large-scale behavior of an ''LES'' fluid, i.e., a fluid with artificial properties, will be similar to that of a real fluid, provided the artificial properties obey certain constraints. The artificial properties consist of modifications to the shear viscosity, bulk viscosity, thermal conductivity and species diffusivity of a fluid. The modified transport coefficients are designed to damp out high wavenumber modes, close to the resolution limit, without corrupting lower modes. Requisite behavior of the artificial properties is discussed and results are shown for a variety of test problems, each designed to exercise different aspects of the models. When combined with a 10th-order compact scheme, the overall method exhibits excellent resolution characteristics for turbulent mixing, while capturing shocks and material interfaces in crisp fashion.

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.

The authors discuss the decay of axion walls bounded by strings and present numerical simulations of the decay process. In these simulations, the decay happens immediately, in a time scale of order the light travel time, and the average energy of the radiated axions is <w{sub a}> {approx_equal} 7m{sub a} for v{sub a}/m{sub a} {approx_equal} 500. <w{sub a}> is found to increase approximately linearly with ln(v{sub a}/m{sub a}). Extrapolation of this behavior yields <w{sub a}> {approx_equal} 60 m{sub a} in axion models of interest.

We consider backgrounds to the detection of the two-body hadronic decay modes of neutral B mesons and baryons. The largest background is due to the correlated production of pairs of high-p/sub T/ hadrons in the target, but this can be adequately rejected provided the experimental apparatus has sufficient resolution in mass and decay vertex. Another possible source of background arises from the production and decay of charmed and strange particles. Since these particles can travel considerable distances before decaying, they can give rise to backgrounds which may not be rejectable by means of vertex cut. We have simulated several backgrounds from charm, and we find them to be small compared to the expected level of signal. 8 refs., 1 fig., 2 tabs

None

We describe an approach for characterizing selective binding between oppositely charged ionic functional groups under biologically relevant conditions. Relative shifts in K-shell x-ray absorption spectra of aqueous cations and carboxylate anions indicate the corresponding binding strengths via perturbations of carbonyl antibonding orbitals. XAS spectra measured for aqueous formate and acetate solutions containing lithium, sodium, and potassium cations reveal monotonically stronger binding of the lighter metals, supporting recent results from simulations and other experiments. The carbon K-edge spectra of the acetate carbonyl feature centered near 290 eV clearly indicate a preferential interaction of sodium versus potassium, which was less apparent with formate. These results are in accord with the Law of Matching Water Affinities, relating relative hydration strengths of ions to their respective tendencies to form contact ion pairs. Density functional theory calculations of K-shell spectra support the experimental findings.

The concept of a disturbed zone surrounding the mined openings of a potential geologic repository for high-level radioactive waste was introduced by the US Nuclear Regulatory Commission (NRC) as a region to be excluded for determining groundwater travel time to the accessible environment, but to be included for determining the impact of underground construction and radioactive decay heat on groundwater movement and radionuclide transport for total system performance analysis. This paper explores both the regulatory and technical necessity for characterizing and modeling a larger region -- the altered zone -- within which the temperature is increased significantly by heat from the high-level waste. Particular attention is given to addressing the effects of heterogeneity on groundwater flux and travel time, showing how these effects might be modeled simply on a macroscopic scale, and outlining its parameters. The effect of uncertainty in the parameter values on the performance of a potential repository can then be easily handled by probabilistic analysis.

Fluids have been sampled from 9 wells and 2 fumaroles fromthe East Flank of the Coso hydrothermal system with a view toidentifying, if possible, the location and characteristics of the heatsource inflows into this portion of the geothermal field. Preliminaryresults show that there has been extensive vapor loss in the system, mostprobably in response to production. Wells 38A-9, 51-16 and 83A-16 showthe highest CO2-CO-CH4-H2 chemical equilibration temperatures, rangingbetween 300-340oC, and apart from 38A-9, the values are generally inaccordance with the measured temperatures in the wells. Calculatedtemperatures for the fractionation of 13C between CO2 and CH4 are inexcess of 400oC in fluids from wells 38A-9, 64-16-RD2 and 51A-16,obviously pointing to equilibrium conditions from deeper portions of thereservoir. Given that the predominant reservoir rock lithologies in theCoso system are relatively silicic (granitic to dioritic), the isotopicsignatures appear to reflect convective circulation and equilibrationwithin rocks close to the plastic-brittle transition. 3He/4He signatures,in conjunction with relative volatile abundances in the Coso fluids,point to a possibly altered mantle source for the heat sourcefluids.

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.

Realistic modeling of coherent synchrotron radiation (CSR) and the space charge force in single-pass systems and rings usually requires at least a two-dimensional (2-D) description of the charge/current density of the bunch. Since that leads to costly computations, one often resorts to a 1-D model of the bunch for first explorations. This paper provides several improvements to previous 1-D theories, eliminating unnecessary approximations and physical restrictions.

Composite Analyses (CA's) are required per DOE Order 435.1 [1], in order to provide a reasonable expectation that DOE low-level waste (LLW) disposal, high-level waste tank closure, and transuranic (TRU) waste disposal in combination with Comprehensive Environmental Response, Compensation, and Liability Act (CERCLA), Resource Conservation and Recovery Act (RCRA), and deactivation and decommissioning (D&D) actions, will not result in the need for future remedial actions in order to ensure radiological protection of the public and environment. This Order requires that an accounting of all sources of DOE man-made radionuclides and DOE enhanced natural radionuclides that are projected to remain on the site after all DOE site operations have ceased. This CA updates the previous CA that was developed in 1997. As part of this CA, an inventory of expected radionuclide residuals was conducted, exposure pathways were screened and a model was developed such that a dose to the MOP at the selected points of exposure might be evaluated.

We report on a general method to fabricate transition metal related defects in diamond. Controlled incorporation of Mo and W in synthetic CVD diamond was achieved by adding volatile metal precursors to the diamond chemical vapor deposition (CVD) growth process. Effects of deposition temperature, grain structure and precursor exposure on the doping level were systematically studied, and doping levels of up to 0.25 at.% have been achieved. The metal atoms are uniformly distributed throughout the diamond grains without any indication of inclusion formation. These results are discussed in context of the kinetically controlled growth process of CVD diamond.

In 2008, a new Multi-Probe Corrosion Monitoring System (MPCMS) was installed in double-shell tank 241-AN-102 on the U.S. Department of Energy's Hanford Site in Washington State. Developmental design work included laboratory testing in simulated tank 241-AN-102 waste to evaluate metal performance for installation on the MPCMS as secondary metal reference electrodes. The MPCMS design includes coupon arrays as well as a wired probe which facilitates measurement of tank potential as well as corrosion rate using electrical resistance (ER) sensors. This paper presents the MPCMS design, field data obtained following installation of the MPCMS in tank 241-AN-102, and a comparison between laboratory potential data obtained using simulated waste and tank potential data obtained following field installation.

Several new and updated BABAR measurements of sin 2{beta} are presented, together with the latest constraints on the Unitarity Triangle angles {alpha} and {gamma}. The higher statistics now available allow more sophisticated analysis techniques, such as time-dependent Dalitz plot fitting. Combined world-average results place tight constraints on the Unitarity Triangle. There is good agreement among the measurements and with the unitarity of the CKM matrix. This represents an impressive verification of the Standard Model description of the quark-flavor sector and of CP violation.

Lawrence Livermore National Laboratory is in the process of constructing the National Ignition Facility, a half million square foot facility which will house a 192 beam laser system capable of generating the 2 million joules of ultraviolet light energy necessary to achieve fusion ignition with inertial targets by 2004. More than 7,000 meter class optics will need to be manufactured by LLNL`s industrial partners to construct the laser system. The components will be manufactured starting in 1998 and will be finished by 2003. In 1994 it became clear through a series of funded cost studies that, in order to fabricate such an unprecedented number of large precision optics in so short a time for the lowest possible cost, new technologies would need to be developed and new factories constructed based on those technologies. At that time, LLNL embarked on an ambitious optics finishing technology development program costing more than $6M over 3 years to develop these technologies, working with three suppliers of large precision optics. While each development program centered upon the specialties and often proprietary technologies already existing in the suppliers facility, many of the technologies required for manufacturing large precision optics at the lowest cost possible are common …

In this paper we describe use of the Aquila active well neutron coincidence counter for nuclear material assays of {sup 235}U in multiple analytical techniques at Savannah River Site (SRS), at the Savannah River National Laboratory (SRNL), and at Argonne West National Laboratory (AWNL). The uses include as a portable passive neutron counter for field measurements searching for evidence of {sup 252}Cf deposits and storage; as a portable active neutron counter using an external activation source for field measurements searching for trace {sup 235}U deposits and holdup; for verification measurements of U-Al reactor fuel elements; for verification measurements of uranium metal; and for verification measurements of process waste of impure uranium in a challenging cement matrix. The wide variety of uses described demonstrate utility of the technique for neutron coincidence verification measurements over the dynamic ranges of 100 g-5000 g for U metal, 200 g-1300 g for U-Al, and 8 g-35 g for process waste. In addition to demonstrating use of the instrument in both the passive and active modes, we also demonstrate its use in both the fast and thermal neutron modes.

DIAS, the Dynamic Information Architecture System, is an object-oriented simulation system that was designed to provide an integrating framework in which new or legacy software applications can operate in a context-driven frame of reference. DIAS provides a flexible and extensible mechanism to allow disparate, and mixed language, software applications to interoperate. DIAS captures the dynamic interplay between different processes or phenomena in the same frame of reference. Finally, DIAS accommodates a broad range of analysis contexts, with widely varying spatial and temporal resolutions and fidelity.

The paper contains over eighty rules for editing graphs, arranged under nine major headings in a logical sequence for editing all the graphs in a manuscript. It is excerpted from a monograph used at the Lawrence Livermore National Laboratory to train beginning technical editors in editing graphs; a corresponding Hypercard stack is also used in this training. 6 refs., 4 figs.

Experiments in a rapid compression machine have examined the influences of variations in pressure, temperature, and equivalence ratio on the autoignition of n-pentane. Equivalence ratios included values from 0.5 to � 2.0, compressed gas initial temperatures were varied between 675K and 980K, and compresed gas initial pressures varied from 8 to 20 bar. Numerical simulations of the same experiments were carried out using a detailed chemical kinetic reaction mechanism. The results are interpreted in terms of a low temperature oxidation mechanism involving addition of molecular oxygen to alkyl and hydroperoxyalkyl radicals. Idealized calculations are reported which identify the major reaction paths at each temperature. Results indicate that in most cases, the reactive gases experience a two-stage autoigni tion. The first stage follows a low temperature alkylperoxy radical isomerization pathway that is effectively quenched when the temperature reaches a level where dissociation reactions of alkylperoxy and hydroperoxyalkylperoxy radicals are more rapid than the reverse addition steps. The second stage is controlled by the onset of dissociation of hydrogen peroxide. Results also show that in some cases, the first stage ignition takes place during the compression stroke in the rapid compression machine, making the interpretation of the experiments somewhat more complex than …