The dissertation, titled ''Genetic Algorithms Applied to Nonlinear and Complex Domains'', describes and then applies a new class of powerful search algorithms (GAS) to certain domains. GAS are capable of solving complex and nonlinear problems where many parameters interact to produce a final result such as the optimization of the laser pulse in the interaction of an atom with an intense laser field. GAS can very efficiently locate the global maximum by searching parameter space in problems which are unsuitable for a search using traditional methods. In particular, the dissertation contains new scientific findings in two areas. First, the dissertation examines the interaction of an ultra-intense short laser pulse with atoms. GAS are used to find the optimal frequency for stabilizing atoms in the ionization process. This leads to a new theoretical formulation, to explain what is happening during the ionization process and how the electron is responding to finite (real-life) laser pulse shapes. It is shown that the dynamics of the process can be very sensitive to the ramp of the pulse at high frequencies. The new theory which is formulated, also uses a novel concept (known as the (t,t') method) to numerically solve the time-dependent Schrodinger equation Second, …

Radioactive metal waste makes up a significant portion of the waste currently being sent for disposal. Recovery of this metal as a valuable resource is possible through the use of decontamination technologies. Through the development and use of expert systems a comparison can be made of laser decontamination, a technology currently under development at Ames Laboratory, with currently available decontamination technologies for applicability to the types of metal waste being generated and the effectiveness of these versus simply disposing of the waste. These technologies can be technically and economically evaluated by the use of expert systems techniques to provide a waste management decision making tool that generates, given an identified metal waste, waste management recommendations. The user enters waste characteristic information as input and the system then recommends decontamination technologies, determines residual contamination levels and possible waste management strategies, carries out a cost analysis and then ranks, according to cost, the possibilities for management of the waste. The expert system was developed using information from literature and personnel experienced in the use of decontamination technologies and requires validation by human experts and assignment of confidence factors to the knowledge represented within.

Heavy ion fusion requires injection, transport and acceleration of high current beams. Detailed simulation of such beams requires fully self-consistent space charge fields and three dimensions. WARP3D, developed for this purpose, is a particle-in-cell plasma simulation code optimized to work within the framework of an accelerator`s lattice of accelerating, focusing, and bending elements. The code has been used to study several test problems and for simulations and design of experiments. Two applications are drift compression experiments on the MBE-4 facility at LBL and design of the electrostatic quadrupole injector for the proposed ILSE facility. With aggressive drift compression on MBE-4, anomalous emittance growth was observed. Simulations carried out to examine possible causes showed that essentially all the emittance growth is result of external forces on the beam and not of internal beam space-charge fields. Dominant external forces are the dodecapole component of focusing fields, the image forces on the surrounding pipe and conductors, and the octopole fields that result from the structure of the quadrupole focusing elements. Goal of the design of the electrostatic quadrupole injector is to produce a beam of as low emittance as possible. The simulations show that the dominant effects that increase the emittance are the …

Smooth particle methods are relatively new methods for simulating solid and fluid flows through they have a 20-year history of solving complex hydrodynamic problems in astrophysics, such as colliding planets and stars, for which correct answers are unknown. The results presented in this thesis evaluate the adaptability or fitness of the method for typical hydrocode production problems. For finite hydrodynamic systems, boundary conditions are important. A reflective boundary condition with image particles is a good way to prevent a density anomaly at the boundary and to keep the fluxes continuous there. Boundary values of temperature and velocity can be separately controlled. The gradient algorithm, based on differentiating the smooth particle expression for (u{rho}) and (T{rho}), does not show numerical instabilities for the stress tensor and heat flux vector quantities which require second derivatives in space when Fourier`s heat-flow law and Newton`s viscous force law are used. Smooth particle methods show an interesting parallel linking to them to molecular dynamics. For the inviscid Euler equation, with an isentropic ideal gas equation of state, the smooth particle algorithm generates trajectories isomorphic to those generated by molecular dynamics. The shear moduli were evaluated based on molecular dynamics calculations for the three weighting functions, …

The low temperature and flux dependent growth of ultrathin Ag films on the Si(111)7x7 surface is studied with Reflection High-Energy Electron Diffraction (RHEED). The grazing incidence geometry of RHEED allows for an incident molecular beam normal to the surface, and makes it an ideal surface probe for studying ultrathin film growth in real time. Short-lived oscillations in the diffracted intensity are observed during Ag deposition at 150 K, indicating quasi-layer-by-layer growth mediated by adatom mobility. When the 150 K growth is performed over a wide range of deposition rates F, the peak intensity is observed to scale, i.e. I(Ft) depends only on the total amount deposited, which implies thermally activated diffusion is absent at 150 K. Scaling is not obeyed at higher temperatures (T{ge}473 K) for the growth of the {radical}3{times}{radical}3 R30{degrees} ({radical}3) superstructure. Testing for scaling of the diffracted intensity constitutes a new experimental method which can be applied generally to determine if thermal diffusion is active at a particular temperature. Scaling is consistent with a constant diffusion length R{sub 0}, independent of substrate temperature and deposition rate. The presence of a non-thermal diffusion mechanism (responsible for the constant diffusion length R{sub 0}) is confirmed by monitoring the flux …

This dissertation describes a new algorithm for simulating low frequency, kinetic phenomena in plasmas. DArwin Direct Implicit Particle-in-Cell (DADIPIC), as its name implies, is a combination of the Darwin and direct implicit methods. One of the difficulties in simulating plasmas lies in the enormous disparity between the fundamental scale lengths of a plasma and the scale lengths of the phenomena of interest. The objective is to create models which can ignore the fundamental constraints without eliminating relevant plasma properties. Over the past twenty years several PIC methods have been investigated for overcoming the constraints on explicit electrodynamic PIC. These models eliminate selected high frequency plasma phenomena while retaining kinetic phenomena at low frequency. This dissertation shows that the combination of Darwin and Direct Implicit allows them to operate better than they have been shown to operate in the past. Through the Darwin method the hyperbolic Maxwell`s equations are reformulated into a set of elliptic equations. Propagating light waves do not exist in the formulation so the Courant constraint on the time step is eliminated. The Direct Implicit method is applied only to the electrostatic field with the result that electrostatic plasma oscillations do not have to be resolved for stability. …

Research into laser-driven inertial confinement fusion is now entering a critical juncture with the construction of the National Ignition Facility (NIF) at Lawrence Livermore National Laboratory (LLNL). Many of the remaining unanswered questions concerning NIF involve interactions between lasers and plasmas. With the eventual goal of fusion power in mind, laser-plasma interactions relevant to laser fusion schemes is an important topic in need of further research. This work experimentally addresses some potential shortcuts and pitfalls on the road to laser-driven fusion power. Current plans on NIF have 192 laser beams directed into a small cylindrical cavity which will contain the fusion fuel; to accomplish this the beams must cross in the entrance holes, and this intersection will be in the presence of outward-flowing plasma. To investigate the physics involved, interactions of crossing laser beams in flowing plasmas are investigated with experiments on the Nova laser facility at LLNL. It was found that in a flowing plasma, energy is transferred between two crossing laser beams, and this may have deleterious consequences for energy balance and ignition in NIF. Possible solutions to this problem are presented. A recently-proposed alternative to standard laser-driven fusion, the ''fast ignitor'' concept, is also experimentally addressed in …

This paper describes methods for solid phase extraction and high performance liquid chromatography (HPLC). The following are described: Effects of Resin Sulfonation on the Retention of Polar Organic Compounds in Solid Phase Extraction; Ion-Chromatographic Separation of Alkali Metals In Non-Aqueous Solvents; Cation-Exchange Chromatography in Non-Aqueous Solvents; and Silicalite As a Stationary Phase For HPLC.

Objective is to determine if conventional slow-growth methods will be successful for growing single quasicrystals in this system. Reaction schemes, isopleths, and a ternary liquidus diagram were constructed using DTA, XRD, microscopy, and energy dispersive spectroscopy. Results suggest a peak in the liquidus near the predicted compositions of Al{sub 65}Cu{sub 23}Ru{sub 12} or Al{sub 65}Cu{sub 20}Ru{sub 15} in the face-centered icosahedral phase field. Further study is needed to verify the existence of a congruently melting phase. There appears to be composition for which slow-growth methods will produce a single-quasicrystal. Directional solidification studies were done on 65Al25Cu10Ru.

To overcome the extreme difficulties encountered in the experimental decay studies of proton drip line nuclei, several techniques have been utilized, including a helium-jet transport system, particle identification detectors and mass separation. Improvements to the ion source/extraction region of the He-jet coupled on-line Recoil Atom Mass Analyzer (RAMA) and its target/ion source coupling resulted in significant increases in RAMA efficiencies and its mass resolution, as well as reductions in the overall transit time. At the 88-Inch Cyclotron at LBNL, the decays of {sup 31}Cl, {sup 27}P and {sup 28}P, with half-lives of 150 msec, 260 msec and 270.3 msec, respectively, were examined using a he-jet and low-energy gas {Delta}E-gas {Delta}E-silicon E detector telescopes. Total beta-delayed proton branches of 0.3% and 0.07% in {sup 31}Cl and {sub 27}P, respectively, were estimated. Several proton peaks that had been previously assigned to the decay of {sup 31}Cl were shown to be from the decay of {sup 25}Si. In {sup 27}P, two proton groups at 459 {+-} 14 keV and 610 {+-} 11 keV, with intensities of 7 {+-} 3% and 92 {+-} 4% relative to the main (100%) group were discovered. The Gamow-Teller component of the preceding beta-decay of each observed proton transition …

Indoor air exposures to gaseous contaminants originating in soil can cause large human health risks. To predict and control these exposures, the mechanisms that affect vapor transport in near-surface soils need to be understood. In particular, radon exposure is a concern since average indoor radon concentrations lead to much higher risks than are generally accepted for exposure to other environmental contaminants. This dissertation examines an important component of the indoor radon problem: the impacts of wind on soil-gas and radon transport and entry into buildings. The research includes experimental and modeling studies of wind`s interactions with a building`s superstructure and the resulting soil-gas and radon flows in the surrounding soil. In addition to exploring the effects of steady winds, a novel modeling technique is developed to examine the impacts of fluctuating winds on soil-gas and radon transport.

A one-dimensional reactivity calculation code is developed using first order perturbation theory. The reactivity equation is based on the multi-group transport equation using the discrete ordinates method for angular dependence. In addition to the first order perturbation approximations, the reactivity code uses only the isotropic scattering data, but cross section libraries with higher order scattering data can still be used with this code. The reactivity code obtains all the flux, cross section, and geometry data from the standard interface files created by ONEDANT, a discrete ordinates transport code. Comparisons between calculated and experimental reactivities were done with the central reactivity worth data for Lady Godiva, a bare uranium metal assembly. Good agreement is found for isotopes that do not violate the assumptions in the first order approximation. In general for cases where there are large discrepancies, the discretized cross section data is not accurately representing certain resonance regions that coincide with dominant flux groups in the Godiva assembly. Comparing reactivities calculated with first order perturbation theory and a straight {Delta}k/k calculation shows agreement within 10% indicating the perturbation of the calculated fluxes is small enough for first order perturbation theory to be applicable in the modeled system. Computation time comparisons …

The transport of atoms or molecules over surfaces has been an important area of study for several decades now, with its progress generally limited by the available experimental techniques to characterize the phenomena. A number of methods have been developed over the years to measure surface diffusion yet only very few systems have been characterized to this day mainly due to the physical limitations inherent in these available methods. Even the STM with its astonishing atomically-resolved images of the surface has been limited in terms of its capability to determine mass transport properties. This is because the STM is inherently a ``slow`` instrument, i.e., a finite time is needed for signal averaging in order to produce the image. A need exists for additional surface diffusion measurement techniques, ideally ones which are able to study varied systems and measure a wide range of diffusion rates. The STM (especially because of its highly local nature) presents itself as a promising tool to conduct dynamical studies if its poor time resolution during ``normal operation`` can somehow be overcome. The purpose of this dissertation is to introduce a new technique of using the STM to measure adatom mobility on surfaces -- one with a …

Volume rendering allows us to represent a density cloud with ideal properties (single scattering, no self-shadowing, etc.). Scientific visualization utilizes this technique by mapping an abstract variable or property in a computer simulation to a synthetic density cloud. This thesis extends volume rendering from its limitation of isotropic density clouds to anisotropic and/or noisy density clouds. Design aspects of these techniques are discussed that aid in the comprehension of scientific information. Anisotropic volume rendering is used to represent vector based quantities in scientific visualization. Velocity and vorticity in a fluid flow, electric and magnetic waves in an electromagnetic simulation, and blood flow within the body are examples of vector based information within a computer simulation or gathered from instrumentation. Understand these fields can be crucial to understanding the overall physics or physiology. Three techniques for representing three-dimensional vector fields are presented: Line Bundles, Textured Splats and Hair Splats. These techniques are aimed at providing a high-level (qualitative) overview of the flows, offering the user a substantial amount of information with a single image or animation. Non-homogenous volume rendering is used to represent multiple variables. Computer simulations can typically have over thirty variables, which describe properties whose understanding are useful to …

Trac-II is a two dimensional axisymmetric resistive MHD code. It simulates all three spatial components (r, z, φ) of the magnetic field and fluid velocity vectors, and the plasma is treated as a single fluid with two temperatures (T_e,T_i). In addition, it can optionally include a self-consistent external circuit. Recent modifications to the code include the addition of the 3-T radiation model, a 4-phase (solid-liquid-vapor-plasma) equation of state model (QEOS), a 4-phase electrical/thermal conductivity model, and an implicit solution of poloidal B_z,B_r) magnetic field diffusion. These changes permit a detailed study of fiber and wire array Z-pinches. Specifically, Trac-II is used to study the wire array Z-pinch at the PBFA-Z pulse power generator at Sandia National Laboratory. First, in 1-D we examine the behavior of a single wire in the Z-pinch. Then, using these results as initial radial conditions in 2-D, we investigate the dynamics of wire array configurations in the r-z and r-θ plane. In the r-z plane we examine the growth of the m=0 or ��sausage�� instability in single wires within the array. In the r-θ plane we examine the merging behavior between neighboring wires. Special emphasis is placed on trying to explain how instability growth affects the …

The potential risk of lung cancer has evoked interest in the properties of radon decay products. There are two forms of this progeny: either attached to ambient aerosols, or still in the status of ions/molecules/small clusters. This ``unattached`` activity would give a higher dose per unit of airborne activity than the ``attached`` progeny that are rather poorly deposited. In this thesis, a system for determining unattached radon decay products electrical mobility size distribution by measuring their electrical mobilities was developed, based on the fact that about 88% of {sup 218}Po atoms have unit charge at the end of their recoil after decay from {sup 222}Rn, while the remainder are neutral. Essential part of the setup is the radon-aerosol chamber with the Circular Electrical Mobility Spectrometer (CEMS) inside. CEMS is used for sampling and classifying the charged radioactive clusters produced in the chamber. An alpha- sensitive plastic, CR-39 disk, is placed in CEMS as an inlaid disk electrode and the alpha particle detector. CEMS showed good performance in fine inactive particles` classification. If it also works well for radon decay products, it can offer a convenient size distribution measurement for radioactive ultrafine particles. However, the experiments did not obtain an acceptable …

A group of compounds referred to as metal hydrides, when used as electrode materials, is a less toxic alternative to the cadmium hydroxide electrode found in nickel/cadmium secondary battery systems. For this and other reasons, the nickel/metal hydride battery system is becoming a popular rechargeable battery for electric vehicle and consumer electronics applications. A model of this battery system is presented. Specifically the metal hydride material, LaNi{sub 5}H{sub 6}, is chosen for investigation due to the wealth of information available in the literature on this compound. The model results are compared to experiments found in the literature. Fundamental analyses as well as engineering optimizations are performed from the results of the battery model. In order to examine diffusion limitations in the nickel oxide electrode, a ``pseudo 2-D model`` is developed. This model allows for the theoretical examination of the effects of a diffusion coefficient that is a function of the state of charge of the active material. It is found using present data from the literature that diffusion in the solid phase is usually not an important limitation in the nickel oxide electrode. This finding is contrary to the conclusions reached by other authors. Although diffusion in the nickel oxide …

The highly heterogeneous nature of most geologic media, coupled with the restricted view of the subsurface available through boreholes, makes it difficult to determine the spatial distribution of subsurface hydrologic properties. Without such a description one cannot predict how fluid flow or solute transport will occur through permeable geologic media, and these predictions are critically needed to address many important environmental problems, including toxic chemical spills, leaking underground storage tanks, and long-term radioactive waste isolation. A common concern of these problems is the possible existence of high-permeability pathways connecting the problem to the biosphere. An understanding of flow and transport behavior is also necessary to optimize energy extraction from petroleum or geothermal reservoirs, where identifying low-permeability barriers that compartmentalize reservoirs and hamper efficient resource utilization is a key problem. The present work describes the development and application of a new inverse method for determining the spatial distribution of hydrologic properties (permeability and specific storage) in heterogeneous geologic media, using pressure transients from interference well tests. The method employs fractal concepts to improve efficiency and reliability. It is applicable to any sort of heterogeneous geologic medium in which wells communicate with each other, whether it be porous, fractured, or a combination …

Two-phase cross flow over heat exchanger tubes creates vibrations which contribute greatly to the wear on the tubes. Fluidelastic instability is a major mechanism by which tubes can fail. In this work, the fluidelastic instability of a tube placed in an array subjected to two-phase cross flow has been studied. For the determination of fluidelastic instability, a triangular tube array was used. The tubes were made of acrylic and were 2.2 cm or 2.37 cm in diameter and 20 cm in length. Eighteen tubes and 4 half tubes formed 5 rows with a pitch to diameter ratio of 1.4. All of the tubes except the test tube were rigidly supported at the text section wall. The test tube was flexibly supported with two cantilever beams. By installing cantilever beams horizontally and vertically, drag and lift direction tube vibration were studied. Parameters of tube mass, structural stiffness, natural frequency, and pitch to diameter ratio were varied. The drag coefficients on a rigidly held tube in an array subjected to two-phase cross flow were measured. The tube in an array was located at displaced positions as well as at the normal position in order to study the variation of fluid force as …

Currently one of the most widely used high {Tc} superconductors is the Bi-based compounds Bi{sub 2}Sr{sub 2}CaCu{sub 2}O{sub z} and Bi{sub 2}Sr{sub 2}Ca{sub 2}Cu{sub 3}O{sub z} (known as BSCCO 2212 and 2223 compounds) with {Tc} values of about 85 K and 110 K respectively. Lengths of high performance conductors ranging from 100 to 1000 m long are routinely fabricated and some test magnets have been wound. An additional difficulty here is that although Bi-2212 and Bi-2223 phases exist over a wide range of stoichiometries, neither has been prepared in phase-pure form. So far the most successful method of constructing reliable and robust wires or tapes is the so called powder-in-tube (PIT) technique [1, 2, 3, 4, 5, 6, 7] in which oxide powder of the appropriate stoichiometry and phase content is placed inside a metal tube, deformed into the desired geometry (round wire or flat tape), and annealed to produce the desired superconducting properties. Intermediate anneals are often incorporated between successive deformation steps. Silver is the metal used in this process because it is the most compatible with the reacting phase. In all of the commercial processes for BSCCO, Ag seems to play a special catalytic role promoting the growth …

The chemical-properties of rutherfordium (Rf) and its group 4 homologs were studied by sorption on glass support surfaces coated with cobalt(II)ferrocyanide and by solvent extraction with tributylphosphate (TBP) and triisooctylamine (TIOA). The surface studies showed that the hydrolysis trend in the group 4 elements and the pseudogroup 4 element, lb, decreases in the order Rf>Zr{approx}Hf>Th. This trend was attributed to relativistic effects which predicted that Rf would be more prone to having a coordination number of 6 than 8 in most aqueous solutions due to a destabilization of the 6d{sub 5/2} shell and a stabilization of the 7p{sub l/2} shell. This hydrolysis trend was confirmed in the TBP/HBr solvent extraction studies which showed that the extraction trend decreased in the order Zr>Hf>Rf?Ti for HBr, showing that Rf and Ti did not extract as well because they hydrolyzed more easily than Zr and Hf. The TIOA/HF solvent extraction studies showed that the extraction trend for the group 4 elements decreased in the order Ti>Zr{approx}Hf>Rf, in inverse order from the trend of ionic radii Rf>Zr{approx}Hf>Ti. An attempt was made to produce {sup 263}Rf (a) via the {sup 248}Cm({sup 22}Ne, {alpha}3n) reaction employing thenoyltrifluoroacetone (TTA) solvent extraction chemistry and (b) via the {sup …

Detailed studies of the growth and structure of thin films of metal oxides grown on metal single crystal surfaces using Scanning Tunneling Microscopy (STM) are presented. The oxide overlayer systems studied are iron oxide and titanium oxide on the Pt(III) surface. The complexity of the metal oxides and large lattice mismatches often lead to surface structures with large unit cells. These are particularly suited to a local real space technique such as scanning tunneling microscopy. In particular, the symmetry that is directly observed with the STM elucidates the relationship of the oxide overlayers to the substrate as well as distinguishing, the structures of different oxides.

This dissertation proposal covers research work that started in the spring of 1992. The aim of the research has been to study the structural performance and stability of proposed nuclear waste containers and the enclosed fuel rods to be used in the long term storage of High Level Nuclear Waste (HLNW). This research is in two phases, computational and experimental. The computational phase deals with the linear and nonlinear Finite Element Analysis of the different containers due to various loading conditions during normal handling conditions and due to the effect of long term corrosion while the canister is stored in the drift of a backfilled geological repository. The elastoplastic stability of the nuclear fuel rods were studied under body forces resulting from acceleration vectors at varying angles, resulting from a sudden drop of the canister at an angle onto a hard surface.

A parallel program consists of sets of concurrent and sequential tasks. Often, a reduction (such as array sum) sequentially combines values produced by a parallel computation. Because reductions occur so frequently in otherwise parallel programs, they are good candidates for optimization. Since reductions may introduce dependencies, most languages separate computation and reduction. The Sisal functional language is unique in that reduction is a natural consequence of loop expressions; the parallelism is implicit in the language. Unfortunately, the original language supports only seven reduction operations. To generalize these expressions, the Sisal 90 definition adds user-defined reductions at the language level. Applicable optimizations depend upon the mathematical properties of the reduction. Compilation and execution speed, synchronization overhead, memory use and maximum size influence the final implementation. This paper (1) Defines reduction syntax and compares with traditional concurrent methods; (2) Defines classes of reduction operations; (3) Develops analysis of classes for optimized concurrency; (4) Incorporates reductions into Sisal 1.2 and Sisal 90; (5) Evaluates performance and size of the implementations.