Finite element method was used to analyze the three-point bend experimental data of A533B-1 pressure vessel steel obtained by Sherry, Lidbury, and Beardsmore [1] from -160 to -45 C within the ductile-brittle transition regime. As many researchers have shown, the failure stress ({sigma}{sub f}) of the material could be approximated as a constant. The characteristic length, or the critical distance (r{sub c}) from the crack tip, at which {sigma}{sub f} is reached, is shown to be temperature dependent based on the crack tip stress field calculated by the finite element method. With the J-A{sub 2} two-parameter constraint theory in fracture mechanics, the fracture toughness (J{sub C} or K{sub JC}) can be expressed as a function of the constraint level (A{sub 2}) and the critical distance r{sub c}. This relationship is used to predict the fracture toughness of A533B-1 in the ductile-brittle transition regime with a constant {sigma}{sub f} and a set of temperature-dependent r{sub c}. It can be shown that the prediction agrees well with the test data for wide range of constraint levels from shallow cracks (a/W= 0.075) to deep cracks (a/W= 0.5), where a is the crack length and W is the specimen width.

Sparse parallel factorization is among the most complicated and irregular algorithms to analyze and optimize. Performance depends both on system characteristics such as the floating point rate, the memory hierarchy, and the interconnect performance, as well as input matrix characteristics such as such as the number and location of nonzeros. We present LUsim, a simulation framework for modeling the performance of sparse LU factorization. Our framework uses micro-benchmarks to calibrate the parameters of machine characteristics and additional tools to facilitate real-time performance modeling. We are using LUsim to analyze an existing parallel sparse LU factorization code, and to explore a latency tolerant variant. We developed and validated a model of the factorization in SuperLU_DIST, then we modeled and implemented a new variant of slud, replacing a blocking collective communication phase with a non-blocking asynchronous point-to-point one. Our strategy realized a mean improvement of 11percent over a suite of test matrices.

The ATLAS experiment at the Large Hadron Collider at CERN has a substantial collaborative Education and Outreach project. This article describes its activities and how it promotes physics to students around the world. With the extraordinary possibility to make groundbreaking discoveries, the ATLAS Experiment [1] at the Large Hadron Collider at CERN can play an important role in promoting contemporary physics at school. For many years ATLAS has had a substantial collaborative Education and Outreach (E&O) project in which physicists from various parts of the world take part. When the experiment begins in 2007, students from around the world will be analyzing data using cutting-edge technology. The unprecedented collision energies of the Large Hadron Collider allow ATLAS to decode the 'events' that unfold after the head-on collisions of protons (Fig. 1). The scientific results from these events will reveal much about the basic nature of matter, energy, space, and time. Students and others will be excited as they try to find events that may be signs for dark matter, extra dimensions of space, mini-black holes, string theory, and other fundamental discoveries. Science education and outreach and the promotion of awareness and appreciation of physics research have become important tasks for …

We will provide an overview of the physics which can be addressed by studying fluctuations and correlations in heavy ion collisions. Observables, which have been discussed in the literature will be briefly reviewed and put in context with experiment and information from Lattice QCD. Special attention will be given to the QCD critical point and the first order co-existence region.

As is becoming increasingly clear, the human species evolvedin the East African savannah. Details of the precise evolutionary chainremain unresolved however it appears that the process lasted severalmillion years, culminating with the emergence of modern Homo sapiensroughly 200,000 years ago. Following that final evolutionary developmentmodern Homo sapiens relatively quickly populated the entire world.Clearly modern Homo sapiens is a successful, resourceful and adaptablespecies. In the developed societies, modern humans live an existence farremoved from our evolutionary ancestors. As we have learned over the lastcentury, this "new" lifestyle can often result in unintendedconsequences. Clearly, our modern access to food, shelter, transportationand healthcare has resulted in greatly expanded expected lifespan butthis new lifestyle can also result in the emergence of different kinds ofdiseases and health problems. The environment in modern buildings haslittle resemblance to the environment of the savannah. We strive tocreate environments with little temperature, air movement and lightvariation. Building occupants often express great dissatisfaction withthese modern created environments and a significant fraction even developsomething akin to allergies to specific buildings (sick buildingsyndrome). Are the indoor environments we are creating fundamentallyunhealthy -- when examined from an evolutionary perspective?

The structural approach to joint inversion, entailing common boundaries or gradients, offers a flexible way to invert diverse types of surface-based and/or crosshole geophysical data. The cross-gradients function has been introduced as a means to construct models in which spatial changes in two models are parallel or anti-parallel. Inversion methods that use such structural constraints also provide estimates of non-linear and non-unique field-scale relationships between model parameters. Here, we invert jointly crosshole radar and seismic traveltimes for structurally similar models using an iterative non-linear traveltime tomography algorithm. Application of the inversion scheme to synthetic data demonstrates that it better resolves lithological boundaries than the individual inversions. Tests of the scheme on observed radar and seismic data acquired within a shallow aquifer illustrate that the resultant models have improved correlations with flowmeter data than with models based on individual inversions. The highest correlation with the flowmeter data is obtained when the joint inversion is combined with a stochastic regularization operator, where the vertical integral scale is estimated from the flowmeter data. Point-spread functions shows that the most significant resolution improvements of the joint inversion is in the horizontal direction.

Injection of carbon dioxide (CO{sub 2}) into saline aquifers has been proposed as a means to reduce greenhouse gas emissions (geological carbon sequestration). The injection process can be classified as immiscible displacement of an aqueous phase by a less dense and less viscous gas phase. Under disposal conditions (supercritical CO{sub 2}) the viscosity of carbon dioxide can be less than the viscosity of the aqueous phase by a factor of 15. Because of the lower viscosity, the CO{sub 2} displacement front will have a tendency towards instability so that waves or rounded lobes of saturation may appear and grow into fingers that lead to enhanced dissolution, bypassing, and possibly poor sweep efficiency. This paper presents an analysis, through high-resolution numerical simulations, of the onset of instabilities (viscous fingering) during injection of CO{sub 2} into saline aquifers. We explore the influence of viscosity ratio, relative permeability functions, and capillary pressure on finger growth and spacing. In addition, we address the issues of finger triggering, convergence under grid refinement and boundary condition effects. Simulations were carried out on scalar machines, and on an IBM RS/6000 SP (a distributed-memory parallel computer with 6080 processors) with a parallelized version of TOUGH2.

Simulation studies for transmission of microwaves through electron clouds show good agreement with analytic results. The electron cloud produces a shift in phase of the microwave. Experimental observation of this phenomena would lead to a useful diagnostic tool for accessing the local density of electron clouds in an accelerator. These experiments are being carried out at the CERN SPS and the PEP-II LER at SLAC and is proposed to be done at the Fermilab main injector. In this study, a brief analysis of the phase shift is provided and the results are compared with that obtained from simulations.

Studies of submicron dust using Mie scattering from Nd:YAG lasers and video data of micron to sub-millimeter sized dust on DIII-D tokamak have provided the first data of dust sources and transport during tokamak discharges. During normal operation on DIII-D dust observation rates are low, a few events per discharge or less. The net carbon content of the dust corresponds to a carbon atom density a few orders of magnitude below the core impurity density. Statistical analysis of Mie data collected over months of operation reveal correlation of increased dust rate with increased heating power and impulsive wall loading due to edge localized modes (ELMs) and disruptions. Generation of significant amounts of dust by disruptions is confirmed by the camera data. However, dust production by disruptions alone is insufficient to account for estimated in-vessel dust inventory in DIII-D. After an extended entry vent, thousands of dust particles are observed by cameras in the first 2-3 plasma discharges. Individual particles moving at velocities up to {approx}300 m/s, breakup of larger particles into pieces, and collisions of particles with walls are observed. After {approx}70 discharges, dust levels are reduced to a few events per discharge. In order to calibrate diagnostics and benchmark …

The amounts of hydrocarbon gases trapped in natural hydrate accumulations are enormous, leading to a recent interest in the evaluation of their potential as an energy source. Earlier studies have demonstrated that large volumes of gas can be readily produced at high rates for long times from gas hydrate accumulations by means of depressurization-induced dissociation, using conventional technology and vertical wells. The results of this numerical study indicate that the use of horizontal wells does not confer any practical advantages to gas production from Class 1 deposits. This is because of the large disparity in permeabilities between the hydrate layer (HL) and the underlying free gas zone, leading to a hydrate dissociation that proceeds in a horizontally dominant direction and is uniform along the length of the reservoir. When horizontal wells are placed near the base of the HL in Class 2 deposits, the delay in the evolution of a significant gas production rate outweighs their advantages, which include higher rates and the prevention of flow obstruction problems that often hamper the performance of vertical wells. Conversely, placement of a horizontal well near to top of the HL can lead to dramatic increases in gas production from Class 2 and …

The electronic contribution to friction at semiconductor surfaces was investigated by using a Pt-coated tip with 50nm radius in an atomic force microscope sliding against an n-type GaAs(100) substrate. The GaAs surface was covered by an approximately 1 nm thick oxide layer. Charge accumulation or depletion was induced by the application of forward or reverse bias voltages. We observed a substantial increase in friction force in accumulation (forward bias) with respect to depletion (reverse bias). We propose a model based on the force exerted by the trapped charges that quantitatively explains the experimental observations of excess friction.

Paleooceanographic evidence has been used to postulate that methane from oceanic hydrates may have had a significant role in regulating global climate, implicating global oceanic deposits of methane gas hydrate as the main culprit in instances of rapid climate change that have occurred in the past. However, the behavior of contemporary oceanic methane hydrate deposits subjected to rapid temperature changes, like those predicted under future climate change scenarios, is poorly understood. To determine the fate of the carbon stored in these hydrates, we performed simulations of oceanic gas hydrate accumulations subjected to temperature changes at the seafloor and assessed the potential for methane release into the ocean. Our modeling analysis considered the properties of benthic sediments, the saturation and distribution of the hydrates, the ocean depth, the initial seafloor temperature, and for the first time, estimated the effect of benthic biogeochemical activity. The results show that shallow deposits--such as those found in arctic regions or in the Gulf of Mexico--can undergo rapid dissociation and produce significant methane fluxes of 2 to 13 mol/yr/m{sup 2} over a period of decades, and release up to 1,100 mol of methane per m{sup 2} of seafloor in a century. These fluxes may exceed the …

We report the deflagration behavior of PBX 9501 at pressures up to 300 MPa and temperatures of 150-180 C where the sample has been held at the test temperature for several hours before ignition. The purpose is to determine the effect on the deflagration behavior of material damage caused by prolonged exposure to high temperature. This conditioning is similar to that experienced by an explosive while it being heated to eventual explosion. The results are made more complicated by the presence of a significant thermal gradient along the sample during the temperature ramp and soak. Three major conclusions are: the presence of nitroplasticizer makes PBX 9501 more thermally sensitive than LX-04 with an inert Viton binder; the deflagration behavior of PBX 9501 is more extreme and more inconsistent than that of LX-04; and something in PBX 9501 causes thermal damage to 'heal' as the deflagration proceeds, resulting in a decelerating deflagration front as it travels along the sample.

The 'Criticality HazOp' technique, as developed at Hanford's Plutonium Finishing Plant (PFP), has allowed for efficiencies enabling shortening of the time necessary to complete new or revised criticality safety evaluation reports (CSERs). For example, in the last half of 2007 at PFP, CSER revisions undergoing the 'Criticality HazOp' process were completed at a higher rate than previously achievable. The efficiencies gained through use of the 'Criticality HazOp' process come from the preliminary narrowing of potential scenarios for the Criticality analyst to fully evaluate in preparation of the new or revised CSER, and from the use of a systematized 'Criticality HazOp' group assessment of the relevant conditions to show which few parameter/condition/deviation combinations actually require analytical effort. The 'Criticality HazOp' has not only provided efficiencies of time, but has brought to criticality safety evaluation revisions the benefits of a structured hazard evaluation method and the enhanced insight that may be gained from direct involvement of a team in the process. In addition, involved personnel have gained a higher degree of confidence and understanding of the resulting CSER product.

The electrostatic particle-in-cell codeWARP is currently being expanded in order to study electron cloud effects on the dynamics of the beam in storage rings. Results for the Fermilab main injector (MI) show the existence of a threshold in the electron density beyond which there is rapid emittance growth. The Fermilab MI is being considered for an upgrade as part of the high intensity neutrino source (HINS) effort, which will result in a significant increasing of the bunch intensity relative to its present value, placing it in a regime where electron-cloud effects are expected to become important. Various results from the simulations using WARP are discussed here.

Over the last several years we have predicted and observed plasmas with an index of refraction greater than one in the soft X-ray regime. These plasmas are usually a few times ionized and have ranged from low-Z carbon plasmas to mid-Z tin plasmas. Our main calculational tool has been the average atom code. We have recently observed C{sup 2+} plasmas with an index of refraction greater than one at a wavelength of 46.9 nm (26.44 eV). In this paper we compare the average atom method, AVATOMKG, against two more detailed methods, OPAL and CAK, for calculating the index of refraction for the carbon plasmas and discuss the different approximations used. We present experimental measurements of carbon plasmas that display this anomalous dispersion phenomenon. It is shown that the average atom calculation is a good approximation when the strongest lines dominate the dispersion. However, when weaker lines make a significant contribution, the more detailed calculations such as OPAL and CAK are essential. During the next decade X-ray free electron lasers and other X-ray sources will be available to probe a wider variety of plasmas at higher densities and shorter wavelengths so understanding the index of refraction in plasmas will be even …

The Fermilab main injector (MI) is being considered for an upgrade as part of the high intensity neutrino source (HINS) effort. This upgrade will involve a significant increasing of the bunch intensity relative to its present value. Such an increase will place the MI in a regime in which electron-cloud effects are expected to become important. We have used the electrostatic particle-in-cell code WARP, recently augmented with new modeling capabilities and simulation techniques, to study the dynamics of beam-electron cloud interaction. This work in progress involves a systematic assessment of beam instabilities due to the presence of electron clouds.

The Savannah River National Laboratory (SRNL) has developed a new medium for storage of hydrogen and other gases. This involves fabrication of thin, Porous Walled, Hollow Glass Microspheres (PW-HGMs), with diameters generally in the range of 1 to several hundred microns. What is unique about the glass microballons is that porosity has been induced and controlled within the thin, one micron thick walls, on the scale of 10 to several thousand Angstroms. This porosity results in interesting properties including the ability to use these channels to fill the microballons with special absorbents and other materials, thus providing a contained environment even for reactive species. Gases can now enter the microspheres and be retained on the absorbents, resulting in solid-state and contained storage of even reactive species. Also, the porosity can be altered and controlled in various ways, and even used to filter mixed gas streams within a system. SRNL is involved in about a half dozen different programs involving these PW-HGMs and an overview of some of these activities and results emerging are presented.

Low-defect multilayer coatings are required to fabricate mask blanks for Extreme Ultraviolet Lithography (EUVL). The mask blanks consist of high reflectance E W multilayers on low thermal expansion substrates. A defect density of 0.0025 printable defects/cm{sup 2} for both the mask substrate and the multilayer is required to provide a mask blank yield of 60%. Current low defect multilayer coating technology allows repeated coating-added defect levels of 0.05/cm{sup 2} for defects greater than 90 nm polystyrene latex sphere (PSL) equivalent size for lots of 20 substrates. Extended clean operation of the coating system at levels below 0.08/cm{sup 2} for 3 months of operation has also been achieved. Two substrates with zero added defects in the quality area have been fabricated, providing an existence proof that ultra low defect coatings are possible. Increasing the ion source-to-target distance from 410 to 560 mm to reduce undesired coating of the ion source caused the defect density to increase to 0.2/cm{sup 2}. Deposition and etching diagnostic witness substrates and deposition pinhole cameras showed a much higher level of ion beam spillover (ions missing the sputter target) than expected. Future work will quantify beam spillover, and test designs to reduce spillover, if it is confirmed …

A low-energy {gamma}{gamma} collider has been discussed in the context of a testbed for a {gamma}{gamma} interaction region at the Next Linear Collider(NLC). We consider the production of heavy mesons at such a testbed using Compton-backscattered photons and demonstrate that their production rivals or exceeds those by BELLE, BABAR or LEP where they are produced indirectly via virtual {gamma}{gamma} luminosities.

Edward Teller was one of the great physicists of the twentieth century. His career began just after the key ideas of the quantum revolution of the 1920's were completed, opening vast areas of physics and chemistry to detailed understanding. Thus, his early work in theoretical physics focused on applying the new quantum theory to the understanding of diverse phenomena. These topics included chemical physics, diamagnetism, and nuclear physics. Later, he made key contributions to statistical mechanics, surface physics, solid state, and plasma physics. In many cases, the ideas in these papers are still rich with important ramifications.

We have designed an experimental technique to use on the National Ignition Facility (NIF) laser to achieve very high pressure (P{sub max} > 10 Mbar = 1000 GPa), dense states of matter at moderate temperatures (kT < 0.5 eV = 6000 K), relevant to the core conditions of the giant planets. A discussion of the conditions in the interiors of the giant planets is given, and an experimental design that can approach those conditions is described.

Results are presented on event-by-event fluctuations of transverse momenta p{sub T} in central Pb+Pb interactions at 20A, 30A, 40A, 80A, and 158A GeV. The analysis was performed for charged particles at forward center-of-mass rapidity (1.1 < y*{sub {pi}} < 2.6). Three fluctuation measures were studied: the distribution of average transverse momentum M(p{sub T}) in the event, the {phi}{sub p{sub T}} fluctuation measure, and two-particle transverse momentum correlations. Fluctuations of p{sub T} are small and show no significant energy dependence in the energy range of the CERN Super Proton Synchrotron. Results are compared with QCD-inspired predictions for the critical point, and with the UrQMD model. Transverse momentum fluctuations, similar to multiplicity fluctuations, do not show the increase expected for freeze-out near the critical point of QCD.

We validated one year of Global Forecast System (GFS) predictions of surface meteorological variables (wind speed, air temperature, dewpoint temperature, air pressure) over the entire planet for forecasts extending from zero hours into the future (an analysis) to 36 hours. Approximately 12,000 surface stations world-wide were included in this analysis. Root-Mean-Square- Errors (RMSE) increased as the forecast period increased from zero to 36 hours, but the initial RMSE were almost as large as the 36 hour forecast RMSE for all variables. Typical RMSE were 3 C for air temperature, 2-3mb for sea-level pressure, 3.5 C for dewpoint temperature and 2.5 m/s for wind speed. Approximately 20-40% of the GFS errors can be attributed to a lack of resolution of local features. We attribute the large initial RMSE for the zero hour forecasts to the inability of the GFS to resolve local terrain features that often dominate local weather conditions, e.g., mountain- valley circulations and sea and land breezes. Since the horizontal resolution of the GFS (about 1{sup o} of latitude and longitude) prevents it from simulating these locally-driven circulations, its performance will not improve until model resolution increases by a factor of 10 or more (from about 100 km to …