To provide a compact high-brightness heavy-ion beam source for Heavy Ion Fusion (HIF), we have been experimenting with merging multi-beamlets in an injector which uses an RF plasma source. An array of converging beamlets was use to produce a beam with the envelope radius, convergence, and ellipticity matched to an electrostatic quadrupole (ESQ) channel. Experimental results were in good quantitative agreement with simulation and have demonstrated the feasibility of this concept. The size of a driver-scale injector system using this approach will be several times smaller than one designed using traditional single large-aperture beams. The success of this experiment has possible significant economical and technical impacts on the architecture of HIF drivers.

Numerical issues with modeling transport of chemicals or solute in realistic large-scale subsurface systems have been a serious concern, even with the continual progress made in both simulation algorithms and computer hardware in the past few decades. The problem remains and becomes even more difficult when dealing with chemical transport in a multiphase flow system using coarse, multidimensional regular or irregular grids, because of the known effects of numerical dispersion associated with moving plume fronts. We have investigated several total-variation-diminishing (TVD) or flux-limiter schemes by implementing and testing them in the T2R3D code, one of the TOUGH2 family of codes. The objectives of this paper are (1) to investigate the possibility of applying these TVD schemes, using multi-dimensional irregular unstructured grids, and (2) to help select more accurate spatial averaging methods for simulating chemical transport given a numerical grid or spatial discretization. We present an application example to show that such TVD schemes are able to effectively reduce numerical dispersion.

None

We have modeled magma flow in a dike rising in a crack whose strike runs from a highland or ridge to an adjacent lowland to determine the effect of topography on the flow, using a 3D hydromechanical code, FLAC3D (http://www.itascacg.com). The aperture, a, is calculated as a variable in a sheet of zones of fixed width d during the simulation as a function of model deformation. The permeability tensor of each zone is adjusted at each time step in response to the pressure in the cell according to the relationship k{sub ij} = {delta}{sub ij} {alpha}{sup 3}/12{mu}d, which is obtained by equating the flow through the layer of permeable zones from Darcy's law with Poiseuille's law under the same gradient. The fluid viscosity is {mu}, and the crack width is a We found a distinct tendency for the flow to be diverted away from the highland end of the strike toward the lowland. For the 4-km long strike length we modeled, eruption was offset between 500 and 1250 m toward the lowland from the center of the strike length. Separation of the geometric effect of the topography from the topographic overburden effect on lateral confining stresses at the crack indicates …