Modelica is an object-oriented acausal modeling language that is well positioned to become a de-facto standard for expressing models of complex physical systems. To simulate a model expressed in Modelica, it needs to be translated into executable code. For generating run-time efficient code, such a translation needs to employ algebraic formula manipulations. As the SPARK solver has been shown to be competitive for generating such code but currently cannot be used with the Modelica language, we report in this paper how SPARK's symbolic and numerical algorithms can be implemented in OpenModelica, an open-source implementation of a Modelica modeling and simulation environment. We also report benchmark results that show that for our air flow network simulation benchmark, the SPARK solver is competitive with Dymola, which is believed to provide the best solver for Modelica.

None

The balancing domain decomposition methods by constraints are extended to solving both nonsymmetric, positive definite and symmetric, indefinite linear systems. In both cases, certain nonstandard primal constraints are included in the coarse problems of BDDC algorithms to accelerate the convergence. Under the assumption that the subdomain size is small enough, a convergence rate estimate for the GMRES iteration is established that the rate is independent of the number of subdomains and depends only slightly on the subdomain problem size. Numerical experiments for several two-dimensional examples illustrate the fast convergence of the proposed algorithms.

The Superconducting Magnet group at Lawrence Berkeley National Laboratory has been developing Nb{sub 3}Sn high-field accelerator magnet technology for the last fifteen years. In order to support the magnet R&D effort, we are developing a diagnostic system that can help identify the causes of performance limiting quenches by recording small flux-changes within the magnet prior to quench-onset. These analysis techniques were applied to the test results from recent Nb{sub 3}Sn magnets. This paper will examine various types of events and their distinguishing characteristics. The present measurement techniques are discussed along with the design of a new data acquisition system that will substantially improve the quality of the recorded signals.

An integrated laboratory scale, 15 kW high-temperature electrolysis facility has been developed at the Idaho National Laboratory under the U.S. Department of Energy Nuclear Hydrogen Initiative. Initial operation of this facility resulted in over 400 hours of operation with an average hydrogen production rate of approximately 0.9 Nm3/hr. The integrated laboratory scale facility is designed to address larger-scale issues such as thermal management (feed-stock heating, high-temperature gas handling), multiple-stack hot-zone design, multiple-stack electrical configurations, and other “integral” issues. This paper documents the initial operation of the ILS, with experimental details about heat-up, initial stack performance, as well as long-term operation and stack degradation.

Here we describe a simple micromolding method to construct three-dimensional arrays of organotypic epithelial tissue structures that approximate in vivo histology. An elastomeric stamp containing an array of posts of defined geometry and spacing is used to mold microscale cavities into the surface of type I collagen gels. Epithelial cells are seeded into the cavities and covered with a second layer of collagen. The cells reorganize into hollow tissues corresponding to the geometry of the cavities. Patterned tissue arrays can be produced in 3-4 h and will undergo morphogenesis over the following one to three days. The protocol can easily be adapted to study a variety of tissues and aspects of normal and neoplastic development.

None

The authors show how models of neutrino masses and mixings can be differentiated on the basis of their predictions for {theta}{sub 13} and lepton flavor violation in radiative charged lepton decays and {mu} - e conversion. They illustrate the lepton flavor violation results for five predictive SO(10) SUSY GUT models and point out the relative importance of their heavy right-handed neutrino mass spectra and {theta}{sub 13} predictions.

We have developed complementary soft x-ray scattering and reflectometry techniques that allow for the morphological analysis of thin polymer films without resorting to chemical modification or isotopic 2 labeling. With these techniques, we achieve significant, x-ray energy-dependent contrast between carbon atoms in different chemical environments using soft x-ray resonance at the carbon edge. Because carbon-containing samples absorb strongly in this region, the scattering length density depends on both the real and imaginary parts of the atomic scattering factors. Using a model polymer film of poly(styrene-b-methyl methacrylate), we show that the soft x-ray reflectivity data is much more sensitive to these atomic scattering factors than the soft x-ray scattering data. Nevertheless, fits to both types of data yield useful morphological details on the polymer?slamellar structure that are consistent with each other and with literature values.

For nearly 40 years, aerial radiological search and survey missions have been performed by the United States Department of Energy's (USDOE) Remote Sensing Laboratory (RSL). Originally created in 1967 as Aerial Measurement Operations (AMO), the AMS mission has expanded to include acquiring baseline measurements, performing periodic area monitoring, and responding to radiological emergencies. In an accident scenario, AMS fixed-wing and/or rotary-wing systems can be deployed to map radiological deposition. A fixed-wing system is on standby twenty-fours per day, seven days per week and can be deployed within four hours of notification. It can quickly evaluate high levels of radiation which may constitute immediate health risks. To accomplish its mission the fixed-wing aircraft utilizes the Spectral Aerial Radiological Computer System (SPARCS) which records gross count and spectral information. Data from SPARCS is telemetered to ground stations and secure websites where it can be viewed and evaluated in near-real time. The rotary-wing system deploys following the critical phase of an accident and supports the DOE's Consequence Management Response Team (CMRT) in determining long term consequences of the accident. The rotary wing aircraft utilizes the Radiation and Environmental Data Acquisition and Recording System (REDAR). A 25-liter sodium iodide (NaI) spectral system and precise …

We will start with a brief overview of neutrino oscillation physics with emphasis on the remaining unanswered questions. Next, after mentioning near future reactor and accelerator experiments searching for a non zero {theta}{sub 13}, we will introduce the plans for the next generation of long-baseline accelerator neutrino oscillation experiments. We will focus on experiments utilizing powerful (0.7-2.1 MW) Fermilab neutrino beams, either existing or in the design phase.

This paper describes the status of the stadium-sized National Ignition Facility (NIF), the world's largest laser system and first operational multi-megajoule laser. The 192-beam NIF, located at Lawrence Livermore National Laboratory (LLNL), is 96% complete and scheduled for completion in March 2009. The NIF laser will produce nanosecond laser pulses with energies up to approximately 4 MJ in the infrared (laser wavelength = 1.053-{micro}m) and 2MJ in the ultraviolet (laser wavelength = 0.35-{micro}m). With these energies NIF will access conditions of pressure and temperature not previously available on earth, allowing it to conduct experiments in support of the nation's national security, energy, and fundamental science goals. First ignition experiments at NIF are scheduled for FY2010. This paper will provide an overview of the NIF laser and the ignition, energy, and fundamental science activities at NIF.

Vacuum studies of metal single crystal surfaces using electron and molecular beam scattering revealed that the surface atoms relocate when the surface is clean (reconstruction) and when it is covered by adsorbates (adsorbate induced restructuring). It was also discovered that atomic steps and other low coordination surface sites are active for breaking chemical bonds (H-H, O=O, C-H, C=O and C-C) with high reaction probability. Investigations at high reactant pressures using sum frequency generation (SFG)--vibrational spectroscopy and high pressure scanning tunneling microscopy (HPSTM) revealed bond breaking at low reaction probability sites on the adsorbate-covered metal surface, and the need for adsorbate mobility for continued turnover. Since most catalysts (heterogeneous, enzyme and homogeneous) are nanoparticles, colloid synthesis methods were developed to produce monodispersed metal nanoparticles in the 1-10 nm range and controlled shapes to use them as new model catalyst systems in two-dimensional thin film form or deposited in mesoporous three-dimensional oxides. Studies of reaction selectivity in multipath reactions (hydrogenation of benzene, cyclohexene and crotonaldehyde) showed that reaction selectivity depends on both nanoparticle size and shape. The oxide-metal nanoparticle interface was found to be an important catalytic site because of the hot electron flow induced by exothermic reactions like carbon monoxide oxidation.

Different methods for the calculation of cross sections with many QCD particles are compared. To this end, CSW vertex rules, Berends-Giele recursion and Feynman-diagram based techniques are implemented as well as various methods for the treatment of colours and phase space integration. We find that typically there is only a small window of jet multiplicities, where the CSW technique has efficiencies comparable or better than both of the other two methods.

This paper provides the status of the offshore wind energy project proposals in the United States and describes strategic issues faced by the U.S. wind industry.

Interest in high-quality integral benchmark data is increasing as efforts to quantify and reduce calculational uncertainties accelerate to meet the demands of next generation reactor and advanced fuel cycle concepts. The International Reactor Physics Experiment Evaluation Project (IRPhEP) and the International Criticality Safety Benchmark Evaluation Project (ICSBEP) continue to expand their efforts and broaden their scope to identify, evaluate, and provide integral benchmark data for method and data validation. Benchmark model specifications provided by these two projects are used heavily by the international reactor physics, nuclear data, and criticality safety communities. Thus far, 14 countries have contributed to the IRPhEP, and 20 have contributed to the ICSBEP. The status of the IRPhEP and ICSBEP is discussed in this paper, and the future of the two projects is outlined and discussed. Selected benchmarks that have been added to the IRPhEP and ICSBEP handbooks since PHYSOR’06 are highlighted, and the future of the two projects is discussed.

This paper reports the interim status in identifying and reviewing photovoltaic (PV) codes and standards (C&S) and related electrical activities for grid-connected, high-penetration PV systems with a focus on U.S. electric utility distribution grid interconnection.

Cloud droplet size distributions hence the key microphysical quantities (e.g., radar reflectivity, droplet concentration, liquid water content, relative dispersion, and mean-volume radius) are determined by different physical mechanisms, including pre-cloud aerosols as CCNs, cloud updraft, and various turbulent entrainment-mixing processes. Therefore, different relationships among these microphysical properties are expected in response to these various mechanisms. The effect of turbulent entrainment-mixing processes is particularly vexing, with different entrainment-mixing processes likely leading to different microphysical relationships. Cloud radar has been widely used to infer the cloud liquid water content (L) from the measurement of radar reflectivity (Z) using a Z-L relationship. Existing Z-L expressions have been often obtained empirically, and differ substantially (Khain et al. 2008). The discrepancy among Z-L relations, which has been hindering the application of cloud radar in measuring cloud properties, likely stems from the different relationships between the relevant microphysical properties caused by different physical processes. This study first analyzes the Z-L relationship theoretically, and identify the key microphysical properties that affect this relationship, and then address the effects of various processes on the Z-L relationship by discerning the characteristics of the relationships between the relative dispersion, droplet concentration, liquid water content, and mean-volume radius calculated from in-situ …

Computational science research components were vital parts of the SciDAC-1 accelerator project and are continuing to play a critical role in newly-funded SciDAC-2 accelerator project, the Community Petascale Project for Accelerator Science and Simulation (ComPASS). Recent advances and achievements in the area of computational science research in support of petascale electromagnetic modeling for accelerator design analysis are presented, which include shape determination of superconducting RF cavities, mesh-based multilevel preconditioner in solving highly-indefinite linear systems, moving window using h- or p- refinement for time-domain short-range wakefield calculations, and improved scalable application I/O.

The authors report a measurement of resonance parameters of the orbitally excited (L = 1) narrow B{sup 0} mesons in decays to B{sup (*)+}{pi}{sup -} using 1.7 fb{sup -1} of data collected by the CDF II detector at the Fermilab Tevatron. The mass and width of the B*{sub 2}{sup 0} state are measured to be m(B*{sub 2}{sup 0}) = 5740.2{sub -1.8}{sup +1.7}(stat.){sub -0.8}{sup +0.9}(syst.) MeV/c{sup 2} and {Lambda}(B*{sub 2}{sup 0}) = 22.7{sub -3.2}{sup +3.8}(stat.){sub -10.2}{sup +3.2}(syst.) MeV/c{sub 2}. The mass difference between the B*{sub 2}{sup 0} and B{sub 1}{sup 0} states is measured to be 14.9{sub -2.5}{sup +2.2}(stat.){sub -1.4}{sup +1.2}(syst.) MeV/c{sup 2}, resulting in a B{sub 1}{sup 0} mass of 5725.3{sub -2.2}{sup +1.6}(stat.){sub -1.5}{sup +1.4}(syst.) MeV/c{sup 2}. This is currently the most precise measurement of the masses of these states and the first measurement of the B*{sub 2}{sup 0} width.

Bitmap indexes are known to be efficient for ad-hoc range queries that are common in data warehousing and scientific applications. However, they suffer from the curse of cardinality, that is, their efficiency deteriorates as attribute cardinalities increase. A number of strategies have been proposed, but none of them addresses the problem adequately. In this paper, we propose a novel binned bitmap index that greatly reduces the cost to answer queries, and therefore breaks the curse of cardinality. The key idea is to augment the binned index with an Order-preserving Bin-based Clustering (OrBiC) structure. This data structure significantly reduces the I/O operations needed to resolve records that cannot be resolved with the bitmaps. To further improve the proposed index structure, we also present a strategy to create single-valued bins for frequent values. This strategy reduces index sizes and improves query processing speed. Overall, the binned indexes with OrBiC great improves the query processing speed, and are 3 - 25 times faster than the best available indexes for high-cardinality data.

The adsorption and dissociation of water on a Ru(0001) surface containing a small amount ({le} 3 %) of carbon impurities was studied by scanning tunneling microscopy (STM). Various surface species are formed depending on the temperature. These include molecular H{sub 2}O, H{sub 2}O-C complexes, H, O, OH and CH. Clusters of either pure H{sub 2}O or mixed H{sub 2}O-OH species are also formed. Each of these species produces a characteristic contrast in the STM images and can be identified by experiment and by ab initio total energy calculations coupled with STM image simulations. Manipulation of individual species via excitation of vibrational modes with the tunneling electrons has been used as supporting evidence.

An RF power balance method is used to measure the synchrotron radiation losses and the wake field losses. We present the history of the losses in the Low Energy Ring (LER) and the High Energy Ring (HER) during the last several runs of PEP-II.

This paper assesses the accuracy of the simplified frame cavity conduction/convection and radiation models presented in ISO 15099 and used in software for rating and labeling window products. Temperatures and U-factors for typical horizontal window frames with internal cavities are compared; results from Computational Fluid Dynamics (CFD) simulations with detailed radiation modeling are used as a reference. Four different frames were studied. Two were made of polyvinyl chloride (PVC) and two of aluminum. For each frame, six different simulations were performed, two with a CFD code and four with a building-component thermal-simulation tool using the Finite Element Method (FEM). This FEM tool addresses convection using correlations from ISO 15099; it addressed radiation with either correlations from ISO 15099 or with a detailed, view-factor-based radiation model. Calculations were performed using the CFD code with and without fluid flow in the window frame cavities; the calculations without fluid flow were performed to verify that the CFD code and the building-component thermal-simulation tool produced consistent results. With the FEM-code, the practice of subdividing small frame cavities was examined, in some cases not subdividing, in some cases subdividing cavities with interconnections smaller than five millimeters (mm) (ISO 15099) and in some cases subdividing cavities …