This thesis is in the form of a background research. It contains pictures, schemes and colored graphs. The aim is to chart basic findings in the area of progress in computer physics and the posibility of their subsequent usage or the simulation of real problems in plasma physics. The component of my thesis is a comparison of particular accesses in the course of plasma simulation.

No description available.

We use Particle-in-Cell Method¡]PIC¡^and Monte-Carlo-Collision Method¡]MCC¡^to model a plasma colliding system .By this way¡Awe can realize the behaviors and effects in the motion of plasma collision . Our mainly discuss is that the system will produce a thin layer¡]Sheath¡^between plasma and a workpiece ¡Aand the effects of colliding phenomenon in this thin layer¡]Sheath¡^ .

Computational simulations of air glow discharge plasma in the presence of supersonic flow are presented. The glow discharge model is based on a self-consistent, multi-species, continuum description of the plasma with finite-rate chemistry effects. The glow discharge model is coupled to a compressible Navier-Stokes solver to study the effect of the plasma on the flow and the counter-effect of the flow on the plasma. A finite-rate air chemistry model is presented and validated against experiments from the literature at a pressure of 600 mTorr. Computational results are compared with experimentally measured V-I characteristics, axial positive ion densities and electron temperature, and reasonably good qualitative and quantitative agreement is observed. The validated air plasma model is then used to study the effect of the surface plasma discharge on M=3 supersonic flow at freestream pressure 18 Torr and the corresponding effects of the flow on the discharge structure in two dimensions. The species concentrations and the gas temperature are examined in the absence and presence of bulk supersonic flow. The peak gas temperature from the computations is found to be 1180 K with the surface plasma alone in the absence of flow, and 830 K in the presence of supersonic flow. Results indicate that O- ions can have comparable densities to electrons in the pressure range 1-20 Torr, and that O2- ion densities are at least two orders of magnitude smaller over the pressure range considered. Different ion species are found to be dominant in the absence and presence of supersonic flow, highlighting the importance of including finite-rate chemistry effects in discharge models for understanding plasma actuator physical phenomena. Electrode polarity effects are investigated, and the cathode upstream actuation is found to be stronger than the actuation strength with the cathode downstream, which is consistent with experimental findings of several groups. A parallel computing implementation of the plasma and flow simulation tools has been developed and is used to study the three-dimensional plasma actuator configuration with circular pin electrodes. / text

Plasma-based accelerators offer the potential to achieve accelerating gradients orders of magnitude higher than are typical in conventional accelerators. A Plasma Accelerator Research Station has been proposed using the CLARA accelerator at Daresbury Laboratory. In this thesis, theory and the results of particle-in-cell simulations are presented investigating experiments that could be conducted using CLARA as well as the preceding VELA and CLARA Front End. Plasma wakefield acceleration was found to be viable with both CLARA and CLARA Front End, with accelerating gradients of GV/m and 100 MV/m scale respectively. Drive-witness and tailored bunch structures based on the CLARA bunch were also investigated. Plasma focus- ing of the VELA and CLARA Front End bunches was studied in simulations, showing that substantial focusing gradient could be achieved using a passive plasma lens. A plasma beam dump scheme using varying plasma density is also presented. This scheme allows the performance of a passive plasma beam dump to be maintained as the bunch is decelerated and has some advantages over a previously proposed method.

Recently, experimental observations have shown that Polar Mesospheric Summer Echoes PMSE may be modulated by radio wave heating the irregularity source region with a ground-based ionospheric heating facilities. It is clear from these past investigations that the temporal behavior of PMSE during ionospheric heating shows promise as a diagnostic for the associated dust layer. To investigate the temporal behavior of plasma irregularities thought to produce PMSE, this work describes a new model that incorporates both finite diffusion time effects as well as dust charging. The hybrid model utilizes fluid ions described by continuity and momentum equations, electrons whose behavior is determined from quasi-neutrality, and charged dust described by the standard Particle-In-Cell PIC method. The model has been used to investigate the temporal behavior of charged dust associated electron irregularities during electron temperature enhancement associated with radio wave heating. The model predicts that the temporal behavior of the irregularities depends on the ratio of the electron-ion ambipolar diffusion time to the dust particle charging time Td/Tc. The results indicate that typically for Td/Tc << 1, an enhancement in electron irregularity amplitude occurs for a period after turn-off of the radio wave heating. The work also predicts that for Td/Tc >> 1, an enhancement in electron irregularity amplitude occurs for a time period after the turn-on of the radio wave heating. Due to the dependence of Td on irregularity scale-size, these results have important implications for observations of PMSE modification at different radar frequencies. Both continuous and discrete charging model were embedded into this computational model, the results were compared and analyzed. It is evident that significant diagnostic information may be available about the dust layer from the temporal behavior of the electron irregularities during the heating process which modifies the background electron temperature. Particularly interesting and important periods of the temporal behavior are during the turn-on and turn-off of the radio wave heating. Although a number of past theoretical and experimental investigations have considered both these on and off period, this dissertation considers further possibilities for diagnostic information available as well as the underlying physical processes. Approximate analytical models are developed and compared to a more accurate full computational model as a reference. Then from the temporal behavior of the electron irregularities during the turn-on and turn-off of the radio wave heating, the analytical models are used to obtain possible diagnostic information for various charged dust and background plasma quantities. Finally, two experiment campaigns have been performed at HAARP, Gakona, Alaska. Preliminary observation results look promising for the existence of PMSE turn-on overshoot. However, more careful experiments need to be done before firm conclusions can be drawn. The new designed Echotek digital receiver is ready for use now. It will be much superior to the experimental setup used for measurements in the previous campaign.Therefore, future experimental campaigns are planning next year to support the theoretical research. / Ph. D.

This thesis describes development of the DRACO plasma simulation code. DRACO is an electro-static (ES) code which uses the particle-in-cell (PIC) formulation to track plasma particles through a computational domain, and operates within the Air Force COLISEUM framework. The particles are tracked on a non-standard mesh, which combines the benefits of a Cartesian mesh with the surface-resolving power of an unstructured mesh. DRACO contains its own mesher, called VOLCAR, which is also described in this work. DRACO was applied to a fully kinetic simulation of an ion-beam neutralization. The thruster configuration and running parameters were based on the NASA's 40cm NEXT ion thruster. The neutralization process was divided into three steps. Electron dynamics was studied by assuming an initial beam neutralization, which was accomplished by injecting both electrons and ions from the optics. Performing the simulation on a full-sized domain with cell size much greater than the Debye length resulted in a formation of a virtual anode. Decrease of the cell size to match the Debye length was not feasible, since it would require a million-fold increase in the number of simulation nodes. Instead, a scaling scheme was devised. Simulations were performed on thruster scaled down by a factor of 100, but its operating parameters were also adjusted such that the produced plasma environment did not change. Loss of electrons at the boundary of the finite simulation domain induced a numerical instability. The instability resulted in a strong axial electric field which sucked out electrons from the beam. It was removed by introducing an energy based particle boundary condition. Combination of surface scaling and energy boundary resulted in physically sound simulation results. Comparison were made between the Maxwellian and polytropic temperatures, as well as between simulation electron density and one predicted by the Boltzmann relationship. The cathode was modeled individually from the beam by introducing a positively charged collector plate at a distance corresponding to the beam edge. The local Debye length at the cathode tip was too small to be resolved by the mesh, even if mesh-refinement was incorporated. Since the simulation was not concerned with the near-tip region, two modifications were performed. First, the a limiting value of charge density at the tip was imposed. Second, the cathode potential was allowed to float. These two modifications were necessary to prevent development of a strong potential gradient at the cathode tip. The modified cathode model was combined with ion injection from the optics to model the actual beam neutralization. Three configurations were tested: a single thruster, a 2x2 cluster with individual cathodes and a similar cluster with a single large neutralizer. Neither of the cases achieved neutralization comparable to one in the base-line pre-neutralized case. The reason for the discrepancy is not known, but it does not seem to be due a loss of electrons at the walls. The difference could be due to limited extent of the modeled physics. An additional work is required to answer this question. / Master of Science

A two-dimensional, time-dependent model and code have been developed to model electrothermal (ET) plasma discharges. ET plasma discharges are capillary discharges that draw tens of kA of electric current. The current heats the plasma, and the plasma radiates energy to the capillary walls. The capillary walls ablate by melting and vaporizing and by sublimation. The newly developed model and code is called the Three-fluid, 2D Electrothermal Plasma Flow Simulator (THOR). THOR simulates the electron, ion, and neutral species as separate fluids coupled through interaction terms. The two-dimensional modeling capabilities made available in this new code represent a tool for the exploration and analysis of the physics involved in ET plasma discharges that has never before been available. Previous simulation models of ET plasma discharges have relied primarily on a 1D description of the plasma. These models have often had to include a tunable correction factor to account for the vapor shield layer - a layer of cold ablated vapor separating the plasma core from the ablating surface and limiting the radiation heat flux to the capillary wall. Some studies have incorporated a 2D description of the plasma boundary layer and shown that the effects of a vapor shield layer can be modeled using this 2D description. However, these 2D modeling abilities have not been extended to the simulation of pulsed ET plasma discharges. The development of a fully-2D and time-dependent simulation model of an entire ET plasma source has enabled the investigation of the 2D development of the vapor shield layer and direct comparison with experiments. In addition, this model has provided novel insight into the inherently 2D nature of the internal flow characteristics involved within the plasma channel in an ET plasma discharge. The model is also able to capture the effects of inter-species interactions. This work focuses on the development of the THOR model. The model has been implemented using C++ and takes advantage of modern supercomputing resources. The THOR model couples the 2D hydrodynamics and the interactions of the plasma species through joule heating, ionization, recombination, and elastic collisions. The analysis of simulation results focuses on emergent internal flow characteristics, direct simulation of the vapor shield layer, and the investigation of source geometry effects on simulated plasma parameters. The effect of elastic collisions between electrons and heavy species are shown to affect internal flow characteristics and cause the development of back-flow inside the ET plasma source. The development of the vapor shield layer has been captured using the diffusion approximation for radiation heat transfer within the ET plasma source with simulated results matching experimental measurements. The relationship between source radius and peak current density inside ET plasma discharges has also been explored, and the transition away from the ablation-controlled operation of ET plasma discharges has been observed. / Ph. D.

A complete model of helium-like line and continuum emission in astrophysical plasmas has been incorporated into the plasma simulation code CLOUDY. All elements between He and Zn are treated, any number of levels can be considered, and a full treatment of radiative and collisional processes is included. This includes photoionization from all levels, line transfer including continuum pumping and destruction by background opacities, scattering, and collisional processes. The model is calculated self-consistently with the ionization and thermal structure of the surrounding nebula. The result is a complete line and continuum spectrum of the nebula. The model helium atom is described and compared to a second standalone helium atom in the low-density case. The effects of the mixing of singlet and triplet terms, the truncation of the physical system, and the convergence of the predicted line intensities as a function of the number of quantum levels explicitly included are considered. New Case-B emissivities are calculated for the helium atom at a range of electron temperatures and densities common in planetary nebulae. Observations of the Orion Nebula are analyzed and compared with predictions of the model helium atom. Observations of low-metallicity extragalactic objects by other authors are analyzed. The methods and details of the model helium-like ions are described. The standard X-ray diagnostics of these ions are revisited and augmented with semi-analytical and numerical calculations of ultraviolet line diagnostics. Finally, a new interface between CLOUDY and the X-ray spectral analysis tool XSPEC is discussed.

No description available.