A key challenge in tracking moving objects is the correspondence problem, that is, the correct propagation of object labels from one time step to another. This is especially true when the objects are non-rigid structures, changing shape, and merging and splitting over time. In this work, we describe a general approach to tracking thousands of non-rigid structures in an image sequence. We show how we can minimize memory requirements and generate accurate results while working with only two frames of the sequence at a time. We demonstrate our results using data from computer simulations of a fluimix problem.

This paper describes the mulitblade coordinate transformation (MBC) modeling process that integrates the dynamics of individual wind turbine blades and expresses them as fixed frames.

We present a detailed angle-resolved photoemission spectroscopy (ARPES) investigation of the RTe3 family, which sets this system as an ideal"textbook" example for the formation of a nesting driven charge density wave (CDW). This family indeed exhibits the full range of phenomena that can be associated to CDWinstabilities, from the opening of large gaps on the best nested parts of Fermi surface (up to 0.4 eV), to the existence of residual metallic pockets. ARPES is the best suited technique to characterize these features, thanks to its unique ability to resolve the electronic structure in k space. An additional advantage of RTe3 is that theband structure can be very accurately described by a simple two dimensional tight-binding (TB) model, which allows one to understand and easily reproduce many characteristics of the CDW. In this paper, we first establish the main features of the electronic structure by comparing our ARPES measurements with the linear muffin-tinorbital band calculations. We use this to define the validity and limits of the TB model. We then present a complete description of the CDW properties and of their strong evolution as a function of R. Using simple models, we are able to reproduce perfectly the evolution of gaps …

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.

The last decade has witnessed a rapid proliferation of superscalar cache-based microprocessors to build high-end capability and capacity computers primarily because of their generality, scalability, and cost effectiveness. However, the constant degradation of superscalar sustained performance, has become a well-known problem in the scientific computing community. This trend has been widely attributed to the use of superscalar-based commodity components who's architectural designs offer a balance between memory performance, network capability, and execution rate that is poorly matched to the requirements of large-scale numerical computations. The recent development of massively parallel vector systems offers the potential to increase the performance gap for many important classes of algorithms. In this study we examine four diverse scientific applications with the potential to run at ultrascale, from the areas of plasma physics, material science, astrophysics, and magnetic fusion. We compare performance between the vector-based Earth Simulator (ES) and Cray X1, with leading superscalar-based platforms: the IBM Power3/4 and the SGI Altix. Results demonstrate that the ES vector systems achieve excellent performance on our application suite - the highest of any architecture tested to date.

Properties of finite nuclei are evaluated with two-nucleon (NN) and three-nucleon (NNN) interactions derived within chiral effective field theory (EFT). The nuclear Hamiltonian is fixed by properties of the A = 2 system, except for two low-energy constants (LECs) that parameterize the short range NNN interaction. We constrain those two LECs by a fit to the A = 3 system binding energy and investigate sensitivity of {sup 4}He, {sup 6}Li, {sup 10,11}B and {sup 12,13}C properties to the variation of the constrained LECs. We identify a preferred choice that gives globally the best description. We demonstrate that the NNN interaction terms significantly improve the binding energies and spectra of mid-p-shell nuclei not just with the preferred choice of the LECs but even within a wide range of the constrained LECs. At the same time, we find that a very high quality description of these nuclei requires further improvements to the chiral Hamiltonian.

Large-scale water injection at The Geysers, California, hasgenerated substantial benefits in terms of sustaining reservoir pressuresand production rates, as well as improving steam composition by reducingthe content of non-condensible gases (NCGs). Two effects have beenrecognized and discussed in the literature as contributing to improvedsteam composition, (1) boiling of injectate provides a source of "clean"steam to production wells, and (2) pressurization effects induced byboiling of injected water reduce upflow of native steam with large NCGconcentrations from depth. In this paper we focus on a possibleadditional effect that could reduce NCGs in produced steam by dissolutionin a condensed aqueous phase.Boiling of injectate causes pressurizationeffects that will fairly rapidly migrate outward, away from the injectionpoint. Pressure increases will cause an increase in the saturation ofcondensed phase due to vapor adsorption on mineral surfaces, andcapillary condensation in small pores. NCGs will dissolve in theadditional condensed phase which, depending upon their solubility, mayreduce NCG concentrations in residual steam.We have analyzed thepartitioning of HCl between vapor and aqueous phases, and have performednumerical simulations of injection into superheated vapor zones. Oursimulations provide evidence that dissolution in the condensed phase canindeed reduce NCG concentrations in produced steam.

We study the effect of interference between the Standard Model Higgs boson resonance and the continuum background in the process {gamma}{gamma} {yields} H {yields} b{bar b} at a photon collider. Taking into account virtual gluon exchange between the final-state quarks, we calculate the leading corrections to the height of the resonance for the case of a light (m{sub H} < 160 GeV) Higgs boson. We find that the interference is destructive and around 0.1-0.2% of the peak height, depending on the mass of the Higgs and the scattering angle. This suppression is smaller by an order of magnitude than the anticipated experimental accuracy at a photon collider. However, the fractional suppression can be significantly larger if the Higgs coupling to b quarks is increased by physics beyond the Standard Model.

A new iterative method is developed to numerically calculate the periodic, matched beam envelope solution of the coupled Kapchinskij-Vladimirskij (KV) equations describing the transverse evolution of a beam in a periodic, linear focusing lattice of arbitrary complexity. Implementation of the method is straightforward. It is highly convergent and can be applied to all usual parameterizations of the matched envelope solutions. The method is applicable to all classes of linear focusing lattices without skew couplings, and also applies to parameters where the matched beam envelope is strongly unstable. Example applications are presented for periodic solenoidal and quadrupole focusing lattices. Convergence properties are summarized over a wide range of system parameters.

The XRD-31 is a fast, windowless X-ray vacuum photodiode developed by EG&G. It is currently the primary fast X-ray detector used to diagnose the X-rays on NIF and OMEGA on the multichannel DANTE spectrometer. The XRD-31 has a dynamic range of less than 1e-12 amps to more than 10 amps. A technique is described to measure the impulse response of the diodes to a 150 fs pulse of 200 nm laser light and a method to calculate the “risetime” for a square pulse and compare it with the computed electron transit time from the photocathode to the anode. Measured response time for 5 XRD-31s assembled in early 2004 was 149.7 ps +-2.75 ps.

In this study, we investigated {gamma}H2AX{sup ser139} and 53BP1{sup ser25}, DNA damage pathway markers, to observe responses to radiation insult. Two Human Mammary Epithelial Cell (HMEC) lines were utilized to research the role of immortalization in DNA damage marker expression, HMEC HMT-3522 (S1) with an infinite lifespan, and a subtype of HMEC 184 (184V) with a finite lifespan. Cells were irradiated with 50 cGy X-rays, fixed with 4% paraformaldehyde after 1 hour repair at 37 C, and processed through immunofluorescence. Cells were visualized with a fluorescent microscope and images were digitally captured using Image-Pro Plus software. The 184V irradiated cells exhibited a more positive punctate response within the nucleus for both DNA damage markers compared to the S1 irradiated cells. We will expand the dose and time course in future studies to augment the preliminary data from this research. It is important to understand whether the process of transformation to immortalization compromises the DNA damage sensor and repair process proteins of HMECs in order to understand what is 'normal' and to evaluate the usefulness of cell lines as experimental models.

EUV optics play a key role in attosecond science since only with higher photon energies is it possible to achieve the wide spectral bandwidth required for ultrashort pulses. Multilayer EUV mirrors have been proposed and are being developed to temporally shape (compress) attosecond pulses. To fully characterize a multilayer optic for pulse applications requires not only knowledge of the reflectivity, as a function of photon energy, but also the reflected phase of the mirror. This work develops the metrologies to determine the reflected phase of an EUV multilayer mirror using the photoelectric effect. The proposed method allows one to determine the optic's impulse response and hence its pulse characteristics.

CH2M HILL Hanford Group, Inc. (the Hanford Tank Farm Operations contractor) and the Department of Energy's Office of River Protection have just completed the first phase of the Hanford Single-Shell Tank RCRA Corrective Action Program. The focus of this first phase was to characterize the nature and extent of past Hanford single-shell tank releases and to characterize the resulting fate and transport of the released contaminants. Most of these plumes are below 20 meters, with some reaching groundwater (at 60 to 120 meters below ground surface [bgs]).

In 2005, the Federal Energy Regulatory Commission (FERC) began discussing regulatory policy for reactive-power procurement and pricing in competitive electricity markets. This paper summarizes findings from a unique, interdisciplinary program of public-interest research that lays a formal foundation for evaluating aspects of FERC staff recommendations and offers early insights that should be useful in guiding policy implementation, specifically by: (1) clarifying the consumers and economic characteristics of reactive power as a basis for creating incentives to appropriately price it, (2) defining specific challenges in creating a competitive market for reactive power as well as new tools needed to help ensure such a market functions efficiently, and (3) demonstrating the importance of accounting for the physical characteristics of the transmission network in planning for reactive power and avoiding the exercise of market power by suppliers.

Previous research has identified links between changes in sea surface temperature (SST) and hurricane intensity. We use climate models to study the possible causes of SST changes in Atlantic and Pacific tropical cyclogenesis regions. The observed SST increases in these regions range from 0.32 to 0.67 C over the 20th century. The 22 climate models examined here suggest that century-timescale SST changes of this magnitude cannot be explained solely by unforced variability of the climate system, even under conservative assumptions regarding the magnitude of this variability. Model simulations that include external forcing by combined anthropogenic and natural factors are generally capable of replicating observed SST changes in both tropical cyclogenesis regions.

Computational grids have the potential for solving large-scale scientific problems using heterogeneous and geographically distributed resources. However, a number of major technical hurdles must be overcome before this goal can be fully realized. One problem critical to the effective utilization of computational grids is efficient job scheduling. Our prior work addressed this challenge by defining a grid scheduling architecture and several job migration strategies. The focus of this study is to explore the impact of data migration under a variety of demanding grid conditions. We evaluate our grid scheduling algorithms by simulating compute servers, various groupings of servers into sites, and inter-server networks, using real workloads obtained from leading supercomputing centers. Several key performance metrics are used to compare the behavior of our algorithms against reference local and centralized scheduling schemes. Results show the tremendous benefits of grid scheduling, even in the presence of input/output data migration - while highlighting the importance of utilizing communication-aware scheduling schemes.

Cylinders and rings fabricated from AerMet{reg_sign} 100 alloy and AISI 1018 steel have been explosively driven to fragmentation in order to determine the fracture strains for these materials under plane strain and uniaxial stress conditions. The phenomena associated with the dynamic expansion and subsequent break up of the cylinders are monitored with high-speed diagnostics. In addition, complementary experiments are performed in which fragments from the explosively driven cylinder are recovered and analyzed to determine the statistical distribution associated with the fragmentation process as well as to determine failure mechanisms. The data are used to determine relevant coefficients for the Johnson-Cook (Hancock-McKenzie) fracture model. Metallurgical analysis of the fragments provides information on damage and failure mechanisms.

In recent work (Nourgaliev, Liou, Theofanous, JCP in press) we demonstrated that numerical simulations of interfacial flows in the presence of strong shear must be cast in dynamically sharp terms (sharp interface treatment or SIM), and that moreover they must meet stringent resolution requirements (i.e., resolving the critical layer). The present work is an outgrowth of that work aiming to overcome consequent limitations on the temporal treatment, which become still more severe in the presence of phase change. The key is to avoid operator splitting between interface motion, fluid convection, viscous/heat diffusion and reactions; instead treating all these non-linear operators fully-coupled within a Newton iteration scheme. To this end, the SIM’s cut-cell meshing is combined with the high-orderaccurate implicit Runge-Kutta and the “recovery” Discontinuous Galerkin methods along with a Jacobian-free, Krylov subspace iteration algorithm and its physics-based preconditioning. In particular, the interfacial geometry (i.e., marker’s positions and volumes of cut cells) is a part of the Newton-Krylov solution vector, so that the interface dynamics and fluid motions are fully-(non-linearly)-coupled. We show that our method is: (a) robust (L-stable) and efficient, allowing to step over stability time steps at will while maintaining high-(up to the 5th)-order temporal accuracy; (b) fully conservative, …

The MINOS experiment uses the intense NuMI beam created at Fermilab and two magnetized tracking calorimeters, one located at Fermilab and one located 735 km away at the Soudan Mine in Minnesota, to make precise measurements of {nu}{sub {mu}} disappearance oscillation parameters. We present recent results from the first two years of NuMI beam operations, including the precise measurement of the atmospheric neutrino oscillation parameters and the search for sterile neutrinos. Future prospects for MINOS will also be discussed, including an improved limit on the {theta}{sub 13} mixing angle by searching for {nu}{sub e} appearance in the {nu}{sub {mu}} beam.

The U.S. Department of Energy's Office of Environmental Management (DOE-EM) was established to achieve the safe and compliant disposition of legacy wastes and facilities from defense nuclear applications. A large majority of these wastes and facilities are 'one-of-a-kind' and unique to DOE. Many of the programs to treat these wastes have been 'first-of-a-kind' and unprecedented in scope and complexity. This has meant that many of the technologies needed to successfully disposition these wastes were not yet developed or required significant re-engineering to be adapted for DOE-EM's needs. The DOE-EM program believes strongly in reducing the technical risk of its projects and has initiated several efforts to reduce those risks: (1) Technology Readiness Assessments to reduce the risks of deployment of new technologies; (2) External Technical Reviews as one of several steps to ensure the timely resolution of engineering and technology issues; and (3) Technical Risk Ratings as a means to monitor and communicate information about technical risks. This paper will present examples of how Technology Readiness Assessments, External Technical Reviews, and Technical Risk Ratings are being used by DOE-EM to reduce technical risks.

Fragmentation is a fundamental process that naturally spans micro to macroscopic scales. Recent advances in algorithms, computer simulations, and hardware enable us to connect the continuum to microstructural regimes in a real simulation through a heterogeneous multiscale mathematical model. We apply this model to the problem of predicting how targets in the NIF chamber dismantle, so that optics and diagnostics can be protected from damage. The mechanics of the initial material fracture depend on the microscopic grain structure. In order to effectively simulate the fragmentation, this process must be modeled at the subgrain level with computationally expensive crystal plasticity models. However, there are not enough computational resources to model the entire NIF target at this microscopic scale. In order to accomplish these calculations, a hierarchical material model (HMM) is being developed. The HMM will allow fine-scale modeling of the initial fragmentation using computationally expensive crystal plasticity, while the elements at the mesoscale can use polycrystal models, and the macroscopic elements use analytical flow stress models. The HMM framework is built upon an adaptive mesh refinement (AMR) capability. We present progress in implementing the HMM in the NIF-ALE-AMR code. Additionally, we present test simulations relevant to NIF targets.

Realistic engineering simulation data often have features that are not optimally resolved due to practical limitations on mesh resolution. To be useful to application engineers, vortex characterization techniques must be sufficiently robust to handle realistic data with complex vortex topologies. In this paper, we present enhancements to the vortex topology identification component of an existing vortex characterization algorithm. The modified techniques are demonstrated by application to three realistic data sets that illustrate the strengths and weaknesses of our approach.

Recent claims of a low energy enhancement in the photon strength function of {sup 96}Mo are investigated. Using the DANCE detector the gamma-ray spectra following resonance neutron capture was measured. The spectrum fitting method was used to indirectly extract a photon strength function from the gamma-ray spectra. No strong low energy enhancement in the photon strength function was found.

I'd like to welcome you to the Agent 2004 conference. As most of you are aware, this conference is the fifth in a series of meetings that began in 1999. A conference followed the next year in 2000. The 2001 conference was skipped because of some conflicts with other conferences, and the conferences have proceeded annually since then. We have the proceedings of the previous conferences available here on CDs. One CD has the proceedings from 1999, 2000, and 2002; the other contains last year's proceedings. The purpose of these conferences is to advance the state of the computational social sciences and to integrate the social sciences with the decision sciences and something that is traditionally known as the management sciences. Those of you in the operations/research area are familiar with the traditional school of modeling simulation that emerged from that scientific area. This conference will bring together a different group of people to talk about the topic of agent-based theories and simulations. This fifth agent conference is one of a group of conferences held annually around the country. Most of you are probably aware of the CASOS Conference held at Carnegie Mellon University, usually in July. UCLA holds the …