RESUMO
Rf-driven ion sources for accelerators and many industrial applications benefit from detailed numerical modeling and simulation of plasma characteristics. For instance, modeling of the Spallation Neutron Source (SNS) internal antenna H(-) source has indicated that a large plasma velocity is induced near bends in the antenna where structural failures are often observed. This could lead to improved designs and ion source performance based on simulation and modeling. However, there are significant separations of time and spatial scales inherent to Rf-driven plasma ion sources, which makes it difficult to model ion sources with explicit, kinetic Particle-In-Cell (PIC) simulation codes. In particular, if both electron and ion motions are to be explicitly modeled, then the simulation time step must be very small, and total simulation times must be large enough to capture the evolution of the plasma ions, as well as extending over many Rf periods. Additional physics processes such as plasma chemistry and surface effects such as secondary electron emission increase the computational requirements in such a way that even fully parallel explicit PIC models cannot be used. One alternative method is to develop fluid-based codes coupled with electromagnetics in order to model ion sources. Time-domain fluid models can simulate plasma evolution, plasma chemistry, and surface physics models with reasonable computational resources by not explicitly resolving electron motions, which thereby leads to an increase in the time step. This is achieved by solving fluid motions coupled with electromagnetics using reduced-physics models, such as single-temperature magnetohydrodynamics (MHD), extended, gas dynamic, and Hall MHD, and two-fluid MHD models. We show recent results on modeling the internal antenna H(-) ion source for the SNS at Oak Ridge National Laboratory using the fluid plasma modeling code USim. We compare demonstrate plasma temperature equilibration in two-temperature MHD models for the SNS source and present simulation results demonstrating plasma evolution over many Rf periods for different plasma temperatures. We perform the calculations in parallel, on unstructured meshes, using finite-volume solvers in order to obtain results in reasonable time.
RESUMO
We examine the appearance of power-law behavior in rooted tree graphs in the context of river networks. It has long been observed that the tails of statistical distributions of upstream areas in river networks, measured above every link, obey a power-law relationship over a range of scales. We examine this behavior by considering a subset of all links, defined as those links which drain complete Strahler basins, where the Strahler order defines a discrete measure of scale, for self-similar networks with both deterministic and random topologies. We find an excellent power-law structure in the tail probabilities for complete Strahler basin areas, over many ranges of scale. We show analytically that the tail probabilities converge to a power law under the assumptions of (1) simple scaling of the distributions of complete Strahler basin areas and (2) application of Horton's law of stream numbers. The convergence to a power law does not occur for all underlying distributions, but for a large class of statistical distributions which have specific limiting properties. For example, underlying distributions which are exponential and gamma distributed, while not power-law scaling, produce power laws in the tail probabilities when rescaled and sampled according to Horton's law of stream numbers. The power-law exponent is given by the expression phi=ln(R(b))/ln(R(A)), where R(b) is the bifurcation ratio and R(A) is the Horton area ratio. It is commonly observed that R(b) approximately equal R(A) in many river basins, implying that the tail probability exponent for complete Strahler basins is close to 1.0.
