Numerical experiments on unstructured PIC stability.

Energy Technology Data Exchange (ETDEWEB)

Day, David Minot

2011-04-01

Particle-In-Cell (PIC) is a method for plasmas simulation. Particles are pushed with Verlet time integration. Fields are modeled using finite differences on a tensor product mesh (cells). The Unstructured PIC methods studied here use instead finite element discretizations on unstructured (simplicial) meshes. PIC is constrained by stability limits (upper bounds) on mesh and time step sizes. Numerical evidence (2D) and analysis will be presented showing that similar bounds constrain unstructured PIC.

Response of plasma facing components in Tokamaks due to intense energy deposition using Particle-In-Cell (PIC) methods

Science.gov (United States)

Genco, Filippo

Damage to plasma-facing components (PFC) due to various plasma instabilities is still a major concern for the successful development of fusion energy and represents a significant research obstacle in the community. It is of great importance to fully understand the behavior and lifetime expectancy of PFC under both low energy cycles during normal events and highly energetic events as disruptions, Edge-Localized Modes (ELM), Vertical Displacement Events (VDE), and Run-away electron (RE). The consequences of these high energetic dumps with energy fluxes ranging from 10 MJ/m2 up to 200 MJ/m 2 applied in very short periods (0.1 to 5 ms) can be catastrophic both for safety and economic reasons. Those phenomena can cause a) large temperature increase in the target material b) consequent melting, evaporation and erosion losses due to the extremely high heat fluxes c) possible structural damage and permanent degradation of the entire bulk material with probable burnout of the coolant tubes; d) plasma contamination, transport of target material into the chamber far from where it was originally picked. The modeling of off-normal events such as Disruptions and ELMs requires the simultaneous solution of three main problems along time: a) the heat transfer in the plasma facing component b) the interaction of the produced vapor from the surface with the incoming plasma particles c) the transport of the radiation produced in the vapor-plasma cloud. In addition the moving boundaries problem has to be considered and solved at the material surface. Considering the carbon divertor as target, the moving boundaries are two since for the given conditions, carbon doesn't melt: the plasma front and the moving eroded material surface. The current solution methods for this problem use finite differences and moving coordinates system based on the Crank-Nicholson method and Alternating Directions Implicit Method (ADI). Currently Particle-In-Cell (PIC) methods are widely used for solving

Modeling an emittance-dominated elliptical sheet beam with a 212-dimensional particle-in-cell code

International Nuclear Information System (INIS)

Carlsten, Bruce E.

2005-01-01

Modeling a 3-dimensional (3-D) elliptical beam with a 212-D particle-in-cell (PIC) code requires a reduction in the beam parameters. The 212-D PIC code can only model the center slice of the sheet beam, but that can still provide useful information about the beam transport and distribution evolution, even if the beam is emittance dominated. The reduction of beam parameters and resulting interpretation of the simulation is straightforward, but not trivial. In this paper, we describe the beam parameter reduction and emittance issues related to the initial beam distribution. As a numerical example, we use the case of a sheet beam designed for use with a planar traveling-wave amplifier for high power generator for RF ranging from 95 to 300GHz [Carlsten et al., IEEE Trans. Plasma Sci. 33 (2005) 85]. These numerical techniques also apply to modeling high-energy elliptical bunches in RF accelerators

Program Package for 3d PIC Model of Plasma Fiber

Science.gov (United States)

Kulhánek, Petr; Břeň, David

2007-08-01

A fully three dimensional Particle in Cell model of the plasma fiber had been developed. The code is written in FORTRAN 95, implementation CVF (Compaq Visual Fortran) under Microsoft Visual Studio user interface. Five particle solvers and two field solvers are included in the model. The solvers have relativistic and non-relativistic variants. The model can deal both with periodical and non-periodical boundary conditions. The mechanism of the surface turbulences generation in the plasma fiber was successfully simulated with the PIC program package.

A Fokker-Planck-Landau collision equation solver on two-dimensional velocity grid and its application to particle-in-cell simulation

Energy Technology Data Exchange (ETDEWEB)

Yoon, E. S.; Chang, C. S., E-mail: cschang@pppl.gov [Princeton Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08543 (United States); Korea Advanced Institute of Science and Technology, Yuseong-gu, DaeJeon 305-701 (Korea, Republic of)

2014-03-15

An approximate two-dimensional solver of the nonlinear Fokker-Planck-Landau collision operator has been developed using the assumption that the particle probability distribution function is independent of gyroangle in the limit of strong magnetic field. The isotropic one-dimensional scheme developed for nonlinear Fokker-Planck-Landau equation by Buet and Cordier [J. Comput. Phys. 179, 43 (2002)] and for linear Fokker-Planck-Landau equation by Chang and Cooper [J. Comput. Phys. 6, 1 (1970)] have been modified and extended to two-dimensional nonlinear equation. In addition, a method is suggested to apply the new velocity-grid based collision solver to Lagrangian particle-in-cell simulation by adjusting the weights of marker particles and is applied to a five dimensional particle-in-cell code to calculate the neoclassical ion thermal conductivity in a tokamak plasma. Error verifications show practical aspects of the present scheme for both grid-based and particle-based kinetic codes.

The use of electromagnetic particle-in-cell codes in accelerator applications

International Nuclear Information System (INIS)

Eppley, K.

1988-12-01

The techniques developed for the numerical simulation of plasmas have numerous applications relevant to accelerators. The operation of many accelerator components involves transients, interactions between beams and rf fields, and internal plasma oscillations. These effects produce non-linear behavior which can be represented accurately by particle in cell (PIC) simulations. We will give a very brief overview of the algorithms used in PIC Codes. We will examine the range of parameters over which they are useful. We will discuss the factors which determine whether a two or three dimensional simulation is most appropriate. PIC codes have been applied to a wide variety of diverse problems, spanning many of the systems in a linear accelerator. We will present a number of practical examples of the application of these codes to areas such as guns, bunchers, rf sources, beam transport, emittance growth and final focus. 8 refs., 8 figs., 2 tabs

Towards the optimization of a gyrokinetic Particle-In-Cell (PIC) code on large-scale hybrid architectures

International Nuclear Information System (INIS)

Ohana, N; Lanti, E; Tran, T M; Brunner, S; Hariri, F; Villard, L; Jocksch, A; Gheller, C

2016-01-01

With the aim of enabling state-of-the-art gyrokinetic PIC codes to benefit from the performance of recent multithreaded devices, we developed an application from a platform called the “PIC-engine” [1, 2, 3] embedding simplified basic features of the PIC method. The application solves the gyrokinetic equations in a sheared plasma slab using B-spline finite elements up to fourth order to represent the self-consistent electrostatic field. Preliminary studies of the so-called Particle-In-Fourier (PIF) approach, which uses Fourier modes as basis functions in the periodic dimensions of the system instead of the real-space grid, show that this method can be faster than PIC for simulations with a small number of Fourier modes. Similarly to the PIC-engine, multiple levels of parallelism have been implemented using MPI+OpenMP [2] and MPI+OpenACC [1], the latter exploiting the computational power of GPUs without requiring complete code rewriting. It is shown that sorting particles [3] can lead to performance improvement by increasing data locality and vectorizing grid memory access. Weak scalability tests have been successfully run on the GPU-equipped Cray XC30 Piz Daint (at CSCS) up to 4,096 nodes. The reduced time-to-solution will enable more realistic and thus more computationally intensive simulations of turbulent transport in magnetic fusion devices. (paper)

A gridding method for object-oriented PIC codes

International Nuclear Information System (INIS)

Gisler, G.; Peter, W.; Nash, H.; Acquah, J.; Lin, C.; Rine, D.

1993-01-01

A simple, rule-based gridding method for object-oriented PIC codes is described which is not only capable of dealing with complicated structures such as multiply-connected regions, but is also computationally faster than classical gridding techniques. Using, these smart grids, vacant cells (e.g., cells enclosed by conductors) will never have to be stored or calculated, thus avoiding the usual situation of having to zero electromagnetic fields within conductors after valuable cpu time has been spent in calculating the fields within these cells in the first place. This object-oriented gridding technique makes use of encapsulating characteristics of actual physical objects (particles, fields, grids, etc.) in C ++ classes and supporting software reuse of these entities through C ++ class inheritance relations. It has been implemented in the form of a simple two-dimensional plasma particle-in-cell code, and forms the initial effort of an AFOSR research project to develop a flexible software simulation environment for particle-in-cell algorithms based on object-oriented technology

Particle-in-cell vs straight-line airflow Gaussian calculations of concentration and deposition of airborne emissions out to 70 km for two sites of differing meteorological and topographical character

International Nuclear Information System (INIS)

Lange, R.; Dickerson, M.A.; Peterson, K.R.; Sherman, C.A.; Sullivan, T.J.

1976-01-01

Two numerical models for the calculation of air concentration and ground deposition of airborne effluent releases are compared. The Particle-in-Cell (PIC) model and the Straight-Line Airflow Gaussian model were used for the simulation. Two sites were selected for comparison: the Hudson River Valley, New York, and the area around the Savannah River Plant, South Carolina. Input for the models was synthesized from meteorological data gathered in previous studies by various investigators. It was found that the PIC model more closely simulated the three-dimensional effects of the meteorology and topography. Overall, the Gaussian model calculated higher concentrations under stable conditions with better agreement between the two methods during neutral to unstable conditions. In addition, because of its consideration of exposure from the returning plume after flow reversal, the PIC model calculated air concentrations over larger areas than did the Gaussian model

Particle-in-cell simulations of Hall plasma thrusters

Science.gov (United States)

Miranda, Rodrigo; Ferreira, Jose Leonardo; Martins, Alexandre

2016-07-01

Hall plasma thrusters can be modelled using particle-in-cell (PIC) simulations. In these simulations, the plasma is described by a set of equations which represent a coupled system of charged particles and electromagnetic fields. The fields are computed using a spatial grid (i.e., a discretization in space), whereas the particles can move continuously in space. Briefly, the particle and fields dynamics are computed as follows. First, forces due to electric and magnetic fields are employed to calculate the velocities and positions of particles. Next, the velocities and positions of particles are used to compute the charge and current densities at discrete positions in space. Finally, these densities are used to solve the electromagnetic field equations in the grid, which are interpolated at the position of the particles to obtain the acting forces, and restart this cycle. We will present numerical simulations using software for PIC simulations to study turbulence, wave and instabilities that arise in Hall plasma thrusters. We have sucessfully reproduced a numerical simulation of a SPT-100 Hall thruster using a two-dimensional (2D) model. In addition, we are developing a 2D model of a cylindrical Hall thruster. The results of these simulations will contribute to improve the performance of plasma thrusters to be used in Cubesats satellites currenty in development at the Plasma Laboratory at University of Brasília.

Particle-in-cell vs straight line Gaussian calculations for an area of complex topography

International Nuclear Information System (INIS)

Lange, R.; Sherman, C.

1977-01-01

Two numerical models for the calculation of time integrated air concentraton and ground deposition of airborne radioactive effluent releases are compared. The time dependent Particle-in-Cell (PIC) model and the steady state Gaussian plume model were used for the simulation. The area selected for the comparison was the Hudson River Valley, New York. Input for the models was synthesized from meteorological data gathered in previous studies by various investigators. It was found that the PIC model more closely simulated the three-dimensional effects of the meteorology and topography. Overall, the Gaussian model calculated higher concentrations under stable conditions. In addition, because of its consideration of exposure from the returning plume after flow reversal, the PIC model calculated air concentrations over larger areas than did the Gaussian model

LPIC++. A parallel one-dimensional relativistic electromagnetic particle-in-cell code for simulating laser-plasma-interaction

International Nuclear Information System (INIS)

Lichters, R.; Pfund, R.E.W.; Meyer-ter-Vehn, J.

1997-08-01

The code LPIC++ presented here, is based on a one-dimensional, electromagnetic, relativistic PIC code that has originally been developed by one of the authors during a PhD thesis at the Max-Planck-Institut fuer Quantenoptik for kinetic simulations of high harmonic generation from overdense plasma surfaces. The code uses essentially the algorithm of Birdsall and Langdon and Villasenor and Bunemann. It is written in C++ in order to be easily extendable and has been parallelized to be able to grow in power linearly with the size of accessable hardware, e.g. massively parallel machines like Cray T3E. The parallel LPIC++ version uses PVM for communication between processors. PVM is public domain software, can be downloaded from the world wide web. A particular strength of LPIC++ lies in its clear program and data structure, which uses chained lists for the organization of grid cells and enables dynamic adjustment of spatial domain sizes in a very convenient way, and therefore easy balancing of processor loads. Also particles belonging to one cell are linked in a chained list and are immediately accessable from this cell. In addition to this convenient type of data organization in a PIC code, the code shows excellent performance in both its single processor and parallel version. (orig.)

Dynamic load balancing in a concurrent plasma PIC code on the JPL/Caltech Mark III hypercube

International Nuclear Information System (INIS)

Liewer, P.C.; Leaver, E.W.; Decyk, V.K.; Dawson, J.M.

1990-01-01

Dynamic load balancing has been implemented in a concurrent one-dimensional electromagnetic plasma particle-in-cell (PIC) simulation code using a method which adds very little overhead to the parallel code. In PIC codes, the orbits of many interacting plasma electrons and ions are followed as an initial value problem as the particles move in electromagnetic fields calculated self-consistently from the particle motions. The code was implemented using the GCPIC algorithm in which the particles are divided among processors by partitioning the spatial domain of the simulation. The problem is load-balanced by partitioning the spatial domain so that each partition has approximately the same number of particles. During the simulation, the partitions are dynamically recreated as the spatial distribution of the particles changes in order to maintain processor load balance

Parallel pic plasma simulation through particle decomposition techniques

International Nuclear Information System (INIS)

Briguglio, S.; Vlad, G.; Di Martino, B.; Naples, Univ. 'Federico II'

1998-02-01

Particle-in-cell (PIC) codes are among the major candidates to yield a satisfactory description of the detail of kinetic effects, such as the resonant wave-particle interaction, relevant in determining the transport mechanism in magnetically confined plasmas. A significant improvement of the simulation performance of such codes con be expected from parallelization, e.g., by distributing the particle population among several parallel processors. Parallelization of a hybrid magnetohydrodynamic-gyrokinetic code has been accomplished within the High Performance Fortran (HPF) framework, and tested on the IBM SP2 parallel system, using a 'particle decomposition' technique. The adopted technique requires a moderate effort in porting the code in parallel form and results in intrinsic load balancing and modest inter processor communication. The performance tests obtained confirm the hypothesis of high effectiveness of the strategy, if targeted towards moderately parallel architectures. Optimal use of resources is also discussed with reference to a specific physics problem [it

Comparison of several algorithms of the electric force calculation in particle plasma models

International Nuclear Information System (INIS)

Lachnitt, J; Hrach, R

2014-01-01

This work is devoted to plasma modelling using the technique of molecular dynamics. The crucial problem of most such models is the efficient calculation of electric force. This is usually solved by using the particle-in-cell (PIC) algorithm. However, PIC is an approximative algorithm as it underestimates the short-range interactions of charged particles. We propose a hybrid algorithm which adds these interactions to PIC. Then we include this algorithm in a set of algorithms which we test against each other in a two-dimensional collisionless magnetized plasma model. Besides our hybrid algorithm, this set includes two variants of pure PIC and the direct application of Coulomb's law. We compare particle forces, particle trajectories, total energy conservation and the speed of the algorithms. We find out that the hybrid algorithm can be a good replacement of direct Coulomb's law application (quite accurate and much faster). It is however probably unnecessary to use it in practical 2D models.

Plasma non-uniformity in a symmetric radiofrequency capacitively-coupled reactor with dielectric side-wall: a two dimensional particle-in-cell/Monte Carlo collision simulation

Science.gov (United States)

Liu, Yue; Booth, Jean-Paul; Chabert, Pascal

2018-02-01

A Cartesian-coordinate two-dimensional electrostatic particle-in-cell/Monte Carlo collision (PIC/MCC) plasma simulation code is presented, including a new treatment of charge balance at dielectric boundaries. It is used to simulate an Ar plasma in a symmetric radiofrequency capacitively-coupled parallel-plate reactor with a thick (3.5 cm) dielectric side-wall. The reactor size (12 cm electrode width, 2.5 cm electrode spacing) and frequency (15 MHz) are such that electromagnetic effects can be ignored. The dielectric side-wall effectively shields the plasma from the enhanced electric field at the powered-grounded electrode junction, which has previously been shown to produce locally enhanced plasma density (Dalvie et al 1993 Appl. Phys. Lett. 62 3207-9 Overzet and Hopkins 1993 Appl. Phys. Lett. 63 2484-6 Boeuf and Pitchford 1995 Phys. Rev. E 51 1376-90). Nevertheless, enhanced electron heating is observed in a region adjacent to the dielectric boundary, leading to maxima in ionization rate, plasma density and ion flux to the electrodes in this region, and not at the reactor centre as would otherwise be expected. The axially-integrated electron power deposition peaks closer to the dielectric edge than the electron density. The electron heating components are derived from the PIC/MCC simulations and show that this enhanced electron heating results from increased Ohmic heating in the axial direction as the electron density decreases towards the side-wall. We investigated the validity of different analytical formulas to estimate the Ohmic heating by comparing them to the PIC results. The widespread assumption that a time-averaged momentum transfer frequency, v m , can be used to estimate the momentum change can cause large errors, since it neglects both phase and amplitude information. Furthermore, the classical relationship between the total electron current and the electric field must be used with caution, particularly close to the dielectric edge where the (neglected

Enhanced quasi-static particle-in-cell simulation of electron cloud instabilities in circular accelerators

Science.gov (United States)

Feng, Bing

Electron cloud instabilities have been observed in many circular accelerators around the world and raised concerns of future accelerators and possible upgrades. In this thesis, the electron cloud instabilities are studied with the quasi-static particle-in-cell (PIC) code QuickPIC. Modeling in three-dimensions the long timescale propagation of beam in electron clouds in circular accelerators requires faster and more efficient simulation codes. Thousands of processors are easily available for parallel computations. However, it is not straightforward to increase the effective speed of the simulation by running the same problem size on an increasingly number of processors because there is a limit to domain size in the decomposition of the two-dimensional part of the code. A pipelining algorithm applied on the fully parallelized particle-in-cell code QuickPIC is implemented to overcome this limit. The pipelining algorithm uses multiple groups of processors and optimizes the job allocation on the processors in parallel computing. With this novel algorithm, it is possible to use on the order of 102 processors, and to expand the scale and the speed of the simulation with QuickPIC by a similar factor. In addition to the efficiency improvement with the pipelining algorithm, the fidelity of QuickPIC is enhanced by adding two physics models, the beam space charge effect and the dispersion effect. Simulation of two specific circular machines is performed with the enhanced QuickPIC. First, the proposed upgrade to the Fermilab Main Injector is studied with an eye upon guiding the design of the upgrade and code validation. Moderate emittance growth is observed for the upgrade of increasing the bunch population by 5 times. But the simulation also shows that increasing the beam energy from 8GeV to 20GeV or above can effectively limit the emittance growth. Then the enhanced QuickPIC is used to simulate the electron cloud effect on electron beam in the Cornell Energy Recovery Linac

Wavenumber spectrum of whistler turbulence: Particle-in-cell simulation

International Nuclear Information System (INIS)

Saito, S.; Gary, S. Peter; Narita, Y.

2010-01-01

The forward cascade of decaying whistler turbulence is studied in low beta plasma to understand essential properties of the energy spectrum at electron scales, by using a two-dimensional electromagnetic particle-in-cell (PIC) simulation. This simulation demonstrates turbulence in which the energy cascade rate is greater than the dissipation rate at the electron inertial length. The PIC simulation shows that the magnetic energy spectrum of forward-cascaded whistler turbulence at electron inertial scales is anisotropic and develops a very steep power-law spectrum which is consistent with recent solar wind observations. A comparison of the simulated spectrum with that predicted by a phenomenological turbulence scaling model suggests that the energy cascade at the electron inertial scale depends on both magnetic fluctuations and electron velocity fluctuations, as well as on the whistler dispersion relation. Thus, not only kinetic Alfven turbulence but also whistler turbulence may explain recent solar wind observations of very steep magnetic spectra at short scales.

Resolution of the Vlasov-Maxwell system by PIC discontinuous Galerkin method on GPU with OpenCL

Directory of Open Access Journals (Sweden)

Crestetto Anaïs

2013-01-01

Full Text Available We present an implementation of a Vlasov-Maxwell solver for multicore processors. The Vlasov equation describes the evolution of charged particles in an electromagnetic field, solution of the Maxwell equations. The Vlasov equation is solved by a Particle-In-Cell method (PIC, while the Maxwell system is computed by a Discontinuous Galerkin method. We use the OpenCL framework, which allows our code to run on multicore processors or recent Graphic Processing Units (GPU. We present several numerical applications to two-dimensional test cases.

FLIP-MHD: A particle-in-cell mehtod for magnetohydrodynamics

International Nuclear Information System (INIS)

Brackbill, J.U.

1990-01-01

A particle-in-cell (PIC) method, FLIP is extended to magnetohydrodynamic (MHD) flow in two dimensions. Particles are used to reduce computational diffusion of the magnetic field. FLIP is an extension of ''classical'' PIC, where particles have mass, but every other property of the fluid is stored on a grid. In FLIP, particles have every property of the fluid, so that they provide a complete Lagrangian description not only to resolve contact discontinuities but also to reduce computational diffusion of linear and angular momentum. The interactions among the particles are calculated on a grid, for convenience and economy. The present study extends FLIP to MHD, by including information about the magnetic field among the attributes of the particles. 6 refs

Discrete particle noise in particle-in-cell simulations of plasma microturbulence

International Nuclear Information System (INIS)

Nevins, W.M.; Hammett, G.W.; Dimits, A.M.; Dorland, W.; Shumaker, D.E.

2005-01-01

Recent gyrokinetic simulations of electron temperature gradient (ETG) turbulence with the global particle-in-cell (PIC) code GTC [Z. Lin et al., Proceedings of the 20th Fusion Energy Conference, Vilamoura, Portugal, 2004 (IAEA, Vienna, 2005)] yielded different results from earlier flux-tube continuum code simulations [F. Jenko and W. Dorland, Phys. Rev. Lett. 89, 225001 (2002)] despite similar plasma parameters. Differences between the simulation results were attributed to insufficient phase-space resolution and novel physics associated with global simulation models. The results of the global PIC code are reproduced here using the flux-tube PIC code PG3EQ [A. M. Dimits et al., Phys. Rev. Lett. 77, 71 (1996)], thereby eliminating global effects as the cause of the discrepancy. The late-time decay of the ETG turbulence and the steady-state heat transport observed in these PIC simulations are shown to result from discrete particle noise. Discrete particle noise is a numerical artifact, so both these PG3EQ simulations and, by inference, the GTC simulations that they reproduced have little to say about steady-state ETG turbulence and the associated anomalous heat transport. In the course of this work several diagnostics are developed to retrospectively test whether a particular PIC simulation is dominated by discrete particle noise

ELECTRON ACCELERATIONS AT HIGH MACH NUMBER SHOCKS: TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS IN VARIOUS PARAMETER REGIMES

Energy Technology Data Exchange (ETDEWEB)

Matsumoto, Yosuke [Department of Physics, Chiba University, Yayoi-cho 1-33, Inage-ku, Chiba 263-8522 (Japan); Amano, Takanobu; Hoshino, Masahiro, E-mail: ymatumot@astro.s.chiba-u.ac.jp [Department of Earth and Planetary Science, University of Tokyo, Hongo 1-33, Bunkyo-ku, Tokyo 113-0033 (Japan)

2012-08-20

Electron accelerations at high Mach number collisionless shocks are investigated by means of two-dimensional electromagnetic particle-in-cell simulations with various Alfven Mach numbers, ion-to-electron mass ratios, and the upstream electron {beta}{sub e} (the ratio of the thermal pressure to the magnetic pressure). We find electrons are effectively accelerated at a super-high Mach number shock (M{sub A} {approx} 30) with a mass ratio of M/m = 100 and {beta}{sub e} = 0.5. The electron shock surfing acceleration is an effective mechanism for accelerating the particles toward the relativistic regime even in two dimensions with a large mass ratio. Buneman instability excited at the leading edge of the foot in the super-high Mach number shock results in a coherent electrostatic potential structure. While multi-dimensionality allows the electrons to escape from the trapping region, they can interact with the strong electrostatic field several times. Simulation runs in various parameter regimes indicate that the electron shock surfing acceleration is an effective mechanism for producing relativistic particles in extremely high Mach number shocks in supernova remnants, provided that the upstream electron temperature is reasonably low.

Deploying electromagnetic particle-in-cell (EM-PIC) codes on Xeon Phi accelerators boards

Science.gov (United States)

Fonseca, Ricardo

2014-10-01

The complexity of the phenomena involved in several relevant plasma physics scenarios, where highly nonlinear and kinetic processes dominate, makes purely theoretical descriptions impossible. Further understanding of these scenarios requires detailed numerical modeling, but fully relativistic particle-in-cell codes such as OSIRIS are computationally intensive. The quest towards Exaflop computer systems has lead to the development of HPC systems based on add-on accelerator cards, such as GPGPUs and more recently the Xeon Phi accelerators that power the current number 1 system in the world. These cards, also referred to as Intel Many Integrated Core Architecture (MIC) offer peak theoretical performances of >1 TFlop/s for general purpose calculations in a single board, and are receiving significant attention as an attractive alternative to CPUs for plasma modeling. In this work we report on our efforts towards the deployment of an EM-PIC code on a Xeon Phi architecture system. We will focus on the parallelization and vectorization strategies followed, and present a detailed performance evaluation of code performance in comparison with the CPU code.

Charge-conserving FEM-PIC schemes on general grids

International Nuclear Information System (INIS)

Campos Pinto, M.; Jund, S.; Salmon, S.; Sonnendruecker, E.

2014-01-01

Particle-In-Cell (PIC) solvers are a major tool for the understanding of the complex behavior of a plasma or a particle beam in many situations. An important issue for electromagnetic PIC solvers, where the fields are computed using Maxwell's equations, is the problem of discrete charge conservation. In this article, we aim at proposing a general mathematical formulation for charge-conserving finite-element Maxwell solvers coupled with particle schemes. In particular, we identify the finite-element continuity equations that must be satisfied by the discrete current sources for several classes of time-domain Vlasov-Maxwell simulations to preserve the Gauss law at each time step, and propose a generic algorithm for computing such consistent sources. Since our results cover a wide range of schemes (namely curl-conforming finite element methods of arbitrary degree, general meshes in two or three dimensions, several classes of time discretization schemes, particles with arbitrary shape factors and piecewise polynomial trajectories of arbitrary degree), we believe that they provide a useful roadmap in the design of high-order charge-conserving FEM-PIC numerical schemes. (authors)

Electromagnetic particle-in-cell (PIC) method for modeling the formation of metal surface structures induced by femtosecond laser radiation

Energy Technology Data Exchange (ETDEWEB)

Djouder, M. [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria); Lamrous, O., E-mail: omarlamrous@mail.ummto.dz [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria); Mitiche, M.D. [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria); Itina, T.E. [Laboratoire Hubert Curien, UMR CNRS 5516/Université Jean Monnet, 18 rue de Professeur Benoît Lauras, 42000 Saint-Etienne (France); Zemirli, M. [Laboratoire de Physique et Chimie Quantique, Université Mouloud Mammeri de Tizi-ouzou, BP 17 RP, 15000 Tizi-Ouzou (Algeria)

2013-09-01

The particle in cell (PIC) method coupled to the finite-difference time-domain (FDTD) method is used to model the formation of laser induced periodic surface structures (LIPSS) at the early stage of femtosecond laser irradiation of smooth metal surface. The theoretical results were analyzed and compared with experimental data taken from the literature. It was shown that the optical properties of the target are not homogeneous and the ejection of electrons is such that ripples in the electron density were obtained. The Coulomb explosion mechanism was proposed to explain the ripples formation under the considered conditions.

Electromagnetic particle-in-cell (PIC) method for modeling the formation of metal surface structures induced by femtosecond laser radiation

International Nuclear Information System (INIS)

Djouder, M.; Lamrous, O.; Mitiche, M.D.; Itina, T.E.; Zemirli, M.

2013-01-01

The particle in cell (PIC) method coupled to the finite-difference time-domain (FDTD) method is used to model the formation of laser induced periodic surface structures (LIPSS) at the early stage of femtosecond laser irradiation of smooth metal surface. The theoretical results were analyzed and compared with experimental data taken from the literature. It was shown that the optical properties of the target are not homogeneous and the ejection of electrons is such that ripples in the electron density were obtained. The Coulomb explosion mechanism was proposed to explain the ripples formation under the considered conditions.

Apar-T: code, validation, and physical interpretation of particle-in-cell results

Science.gov (United States)

Melzani, Mickaël; Winisdoerffer, Christophe; Walder, Rolf; Folini, Doris; Favre, Jean M.; Krastanov, Stefan; Messmer, Peter

2013-10-01

We present the parallel particle-in-cell (PIC) code Apar-T and, more importantly, address the fundamental question of the relations between the PIC model, the Vlasov-Maxwell theory, and real plasmas. First, we present four validation tests: spectra from simulations of thermal plasmas, linear growth rates of the relativistic tearing instability and of the filamentation instability, and nonlinear filamentation merging phase. For the filamentation instability we show that the effective growth rates measured on the total energy can differ by more than 50% from the linear cold predictions and from the fastest modes of the simulation. We link these discrepancies to the superparticle number per cell and to the level of field fluctuations. Second, we detail a new method for initial loading of Maxwell-Jüttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. The formulation of the relativistic Harris equilibrium is generalized to arbitrary temperature and mass ratios. Both are required for the tearing instability setup. Third, we turn to the key point of this paper and scrutinize the question of what description of (weakly coupled) physical plasmas is obtained by PIC models. These models rely on two building blocks: coarse-graining, i.e., grouping of the order of p ~ 1010 real particles into a single computer superparticle, and field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e., quantities depending on p. They derive from the PIC plasma parameter ΛPIC, which we show to behave as ΛPIC ∝ 1/p. We explore two important implications. One is that PIC collision- and fluctuation-induced thermalization times are expected to scale with the number of superparticles per grid cell, and thus to be a factor p ~ 1010 smaller than in real plasmas, a fact that we confirm with

On the elimination of numerical Cerenkov radiation in PIC simulations

International Nuclear Information System (INIS)

Greenwood, Andrew D.; Cartwright, Keith L.; Luginsland, John W.; Baca, Ernest A.

2004-01-01

Particle-in-cell (PIC) simulations are a useful tool in modeling plasma in physical devices. The Yee finite difference time domain (FDTD) method is commonly used in PIC simulations to model the electromagnetic fields. However, in the Yee FDTD method, poorly resolved waves at frequencies near the cut off frequency of the grid travel slower than the physical speed of light. These slowly traveling, poorly resolved waves are not a problem in many simulations because the physics of interest are at much lower frequencies. However, when high energy particles are present, the particles may travel faster than the numerical speed of their own radiation, leading to non-physical, numerical Cerenkov radiation. Due to non-linear interaction between the particles and the fields, the numerical Cerenkov radiation couples into the frequency band of physical interest and corrupts the PIC simulation. There are two methods of mitigating the effects of the numerical Cerenkov radiation. The computational stencil used to approximate the curl operator can be altered to improve the high frequency physics, or a filtering scheme can be introduced to attenuate the waves that cause the numerical Cerenkov radiation. Altering the computational stencil is more physically accurate but is difficult to implement while maintaining charge conservation in the code. Thus, filtering is more commonly used. Two previously published filters by Godfrey and Friedman are analyzed and compared to ideally desired filter properties

Particle-in-cell simulation of Trichel pulses in pure oxygen

International Nuclear Information System (INIS)

Soria-Hoyo, C; Pontiga, F; Castellanos, A

2007-01-01

The development and propagation of Trichel pulses in oxygen have been numerically simulated using an improved fluid particle-in-cell (PIC) method. The numerical method has been optimized to accurately reproduce sequences of about 100 Trichel pulses (∼1 ms). A classical one-dimensional model of negative corona in sphere-to-plane geometry has been used to formulate the continuity equations for electrons and ions. The effects of ionization, attachment and secondary-electron emission from the cathode have all been considered. The electric field has been obtained from the solution of Poisson's equation in two dimensions. Using this model, the temporal and electrical characteristics of Trichel pulses have been investigated, in particular, the relation between applied voltage, pulse frequency and time-averaged current intensity and charge

Acceleration of PIC simulation with GPU

International Nuclear Information System (INIS)

Suzuki, Junya; Shimazu, Hironori; Fukazawa, Keiichiro; Den, Mitsue

2011-01-01

Particle-in-cell (PIC) is a simulation technique for plasma physics. The large number of particles in high-resolution plasma simulation increases the volume computation required, making it vital to increase computation speed. In this study, we attempt to accelerate computation speed on graphics processing units (GPUs) using KEMPO, a PIC simulation code package. We perform two tests for benchmarking, with small and large grid sizes. In these tests, we run KEMPO1 code using a CPU only, both a CPU and a GPU, and a GPU only. The results showed that performance using only a GPU was twice that of using a CPU alone. While, execution time for using both a CPU and GPU is comparable to the tests with a CPU alone, because of the significant bottleneck in communication between the CPU and GPU. (author)

Verification of frequency scaling laws for capacitive radio-frequency discharges using two-dimensional simulations

International Nuclear Information System (INIS)

Vahedi, V.; Birdsall, C.K.; Lieberman, M.A.; DiPeso, G.; Rognlien, T.D.

1993-01-01

Weakly ionized processing plasmas are studied in two dimensions using a bounded particle-in-cell (PIC) simulation code with a Monte Carlo collision (MCC) package. The MCC package models the collisions between charged and neutral particles, which are needed to obtain a self-sustained plasma and the proper electron and ion energy loss mechanisms. A two-dimensional capacitive radio-frequency (rf) discharge is investigated in detail. Simple frequency scaling laws for predicting the behavior of some plasma parameters are derived and then compared with simulation results, finding good agreements. It is found that as the drive frequency increases, the sheath width decreases, and the bulk plasma becomes more uniform, leading to a reduction of the ion angular spread at the target and an improvement of ion dose uniformity at the driven electrode

Particle-in-cell simulation of Trichel pulses in pure oxygen

Energy Technology Data Exchange (ETDEWEB)

Soria-Hoyo, C [Dpto. Electronica y Electromagnetismo, Universidad de Sevilla, Avda. Reina Mercedes s/n, Sevilla 41012 (Spain); Pontiga, F [Dpto. Fisica Aplicada II, Universidad de Sevilla, Avda. Reina Mercedes s/n, Sevilla 41012 (Spain); Castellanos, A [Dpto. Electronica y Electromagnetismo, Universidad de Sevilla, Avda. Reina Mercedes s/n, Sevilla 41012 (Spain)

2007-08-07

The development and propagation of Trichel pulses in oxygen have been numerically simulated using an improved fluid particle-in-cell (PIC) method. The numerical method has been optimized to accurately reproduce sequences of about 100 Trichel pulses ({approx}1 ms). A classical one-dimensional model of negative corona in sphere-to-plane geometry has been used to formulate the continuity equations for electrons and ions. The effects of ionization, attachment and secondary-electron emission from the cathode have all been considered. The electric field has been obtained from the solution of Poisson's equation in two dimensions. Using this model, the temporal and electrical characteristics of Trichel pulses have been investigated, in particular, the relation between applied voltage, pulse frequency and time-averaged current intensity and charge.

SHARP: A Spatially Higher-order, Relativistic Particle-in-cell Code

Energy Technology Data Exchange (ETDEWEB)

Shalaby, Mohamad; Broderick, Avery E. [Department of Physics and Astronomy, University of Waterloo, 200 University Avenue West, Waterloo, ON, N2L 3G1 (Canada); Chang, Philip [Department of Physics, University of Wisconsin-Milwaukee, 1900 E. Kenwood Boulevard, Milwaukee, WI 53211 (United States); Pfrommer, Christoph [Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam (Germany); Lamberts, Astrid [Theoretical Astrophysics, California Institute of Technology, Pasadena, CA 91125 (United States); Puchwein, Ewald, E-mail: mshalaby@live.ca [Institute of Astronomy and Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge, CB3 0HA (United Kingdom)

2017-05-20

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to the fifth order). We validate our algorithm against several test problems—thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.

Global Particle-in-Cell Simulations of Mercury's Magnetosphere

Science.gov (United States)

Schriver, D.; Travnicek, P. M.; Lapenta, G.; Amaya, J.; Gonzalez, D.; Richard, R. L.; Berchem, J.; Hellinger, P.

2017-12-01

Spacecraft observations of Mercury's magnetosphere have shown that kinetic ion and electron particle effects play a major role in the transport, acceleration, and loss of plasma within the magnetospheric system. Kinetic processes include reconnection, the breakdown of particle adiabaticity and wave-particle interactions. Because of the vast range in spatial scales involved in magnetospheric dynamics, from local electron Debye length scales ( meters) to solar wind/planetary magnetic scale lengths (tens to hundreds of planetary radii), fully self-consistent kinetic simulations of a global planetary magnetosphere remain challenging. Most global simulations of Earth's and other planet's magnetosphere are carried out using MHD, enhanced MHD (e.g., Hall MHD), hybrid, or a combination of MHD and particle in cell (PIC) simulations. Here, 3D kinetic self-consistent hybrid (ion particle, electron fluid) and full PIC (ion and electron particle) simulations of the solar wind interaction with Mercury's magnetosphere are carried out. Using the implicit PIC and hybrid simulations, Mercury's relatively small, but highly kinetic magnetosphere will be examined to determine how the self-consistent inclusion of electrons affects magnetic reconnection, particle transport and acceleration of plasma at Mercury. Also the spatial and energy profiles of precipitating magnetospheric ions and electrons onto Mercury's surface, which can strongly affect the regolith in terms of space weathering and particle outflow, will be examined with the PIC and hybrid codes. MESSENGER spacecraft observations are used both to initiate and validate the global kinetic simulations to achieve a deeper understanding of the role kinetic physics play in magnetospheric dynamics.

Particle-in-Cell Codes for plasma-based particle acceleration

CERN Document Server

Pukhov, Alexander

2016-01-01

Basic principles of particle-in-cell (PIC ) codes with the main application for plasma-based acceleration are discussed. The ab initio full electromagnetic relativistic PIC codes provide the most reliable description of plasmas. Their properties are considered in detail. Representing the most fundamental model, the full PIC codes are computationally expensive. The plasma-based acceler- ation is a multi-scale problem with very disparate scales. The smallest scale is the laser or plasma wavelength (from one to hundred microns) and the largest scale is the acceleration distance (from a few centimeters to meters or even kilometers). The Lorentz-boost technique allows to reduce the scale disparity at the costs of complicating the simulations and causing unphysical numerical instabilities in the code. Another possibility is to use the quasi-static approxi- mation where the disparate scales are separated analytically.

Rise time of proton cut-off energy in 2D and 3D PIC simulations

Science.gov (United States)

Babaei, J.; Gizzi, L. A.; Londrillo, P.; Mirzanejad, S.; Rovelli, T.; Sinigardi, S.; Turchetti, G.

2017-04-01

The Target Normal Sheath Acceleration regime for proton acceleration by laser pulses is experimentally consolidated and fairly well understood. However, uncertainties remain in the analysis of particle-in-cell simulation results. The energy spectrum is exponential with a cut-off, but the maximum energy depends on the simulation time, following different laws in two and three dimensional (2D, 3D) PIC simulations so that the determination of an asymptotic value has some arbitrariness. We propose two empirical laws for the rise time of the cut-off energy in 2D and 3D PIC simulations, suggested by a model in which the proton acceleration is due to a surface charge distribution on the target rear side. The kinetic energy of the protons that we obtain follows two distinct laws, which appear to be nicely satisfied by PIC simulations, for a model target given by a uniform foil plus a contaminant layer that is hydrogen-rich. The laws depend on two parameters: the scaling time, at which the energy starts to rise, and the asymptotic cut-off energy. The values of the cut-off energy, obtained by fitting 2D and 3D simulations for the same target and laser pulse configuration, are comparable. This suggests that parametric scans can be performed with 2D simulations since 3D ones are computationally very expensive, delegating their role only to a correspondence check. In this paper, the simulations are carried out with the PIC code ALaDyn by changing the target thickness L and the incidence angle α, with a fixed a0 = 3. A monotonic dependence, on L for normal incidence and on α for fixed L, is found, as in the experimental results for high temporal contrast pulses.

Implementing particle-in-cell plasma simulation code on the BBN TC2000

International Nuclear Information System (INIS)

Sturtevant, J.E.; Maccabe, A.B.

1990-01-01

The BBN TC2000 is a multiple instruction, multiple data (MIMD) machine that combines a physically distributed memory with a logically shared memory programming environment using the unique Butterfly switch. Particle-In-Cell (PIC) plasma simulations model the interaction of charged particles with electric and magnetic fields. This paper describes the implementation of both a 1-D electrostatic and a 2 1/2-D electromagnetic PIC (particle-in-cell) plasma simulation code on a BBN TC2000. Performance is compared to implementations of the same code on the shared memory Sequent Balance and distributed memory Intel iPSC hypercube

Output power fluctuations due to different weights of macro particles used in particle-in-cell simulations of Cerenkov devices

International Nuclear Information System (INIS)

Bao, Rong; Li, Yongdong; Liu, Chunliang; Wang, Hongguang

2016-01-01

The output power fluctuations caused by weights of macro particles used in particle-in-cell (PIC) simulations of a backward wave oscillator and a travelling wave tube are statistically analyzed. It is found that the velocities of electrons passed a specific slow-wave structure form a specific electron velocity distribution. The electron velocity distribution obtained in PIC simulation with a relative small weight of macro particles is considered as an initial distribution. By analyzing this initial distribution with a statistical method, the estimations of the output power fluctuations caused by different weights of macro particles are obtained. The statistical method is verified by comparing the estimations with the simulation results. The fluctuations become stronger with increasing weight of macro particles, which can also be determined reversely from estimations of the output power fluctuations. With the weights of macro particles optimized by the statistical method, the output power fluctuations in PIC simulations are relatively small and acceptable.

Two-dimensional particle-in-cell plasma source ion implantation of a prolate spheroid target

International Nuclear Information System (INIS)

Cheng-Sen, Liu; Hong-Ying, Han; Xiao-Qing, Peng; Ye, Chang; De-Zhen, Wang

2010-01-01

A two-dimensional particle-in-cell simulation is used to study the time-dependent evolution of the sheath surrounding a prolate spheroid target during a high voltage pulse in plasma source ion implantation. Our study shows that the potential contour lines pack more closely in the plasma sheath near the vertex of the major axis, i.e. where a thinner sheath is formed, and a non-uniform total ion dose distribution is incident along the surface of the prolate spheroid target due to the focusing of ions by the potential structure. Ion focusing takes place not only at the vertex of the major axis, where dense potential contour lines exist, but also at the vertex of the minor axis, where sparse contour lines exist. This results in two peaks of the received ion dose, locating at the vertices of the major and minor axes of the prolate spheroid target, and an ion dose valley, staying always between the vertices, rather than at the vertex of the minor axis

Two-dimensional particle-in-cell plasma source ion implantation of a prolate spheroid target

Science.gov (United States)

Liu, Cheng-Sen; Han, Hong-Ying; Peng, Xiao-Qing; Chang, Ye; Wang, De-Zhen

2010-03-01

A two-dimensional particle-in-cell simulation is used to study the time-dependent evolution of the sheath surrounding a prolate spheroid target during a high voltage pulse in plasma source ion implantation. Our study shows that the potential contour lines pack more closely in the plasma sheath near the vertex of the major axis, i.e. where a thinner sheath is formed, and a non-uniform total ion dose distribution is incident along the surface of the prolate spheroid target due to the focusing of ions by the potential structure. Ion focusing takes place not only at the vertex of the major axis, where dense potential contour lines exist, but also at the vertex of the minor axis, where sparse contour lines exist. This results in two peaks of the received ion dose, locating at the vertices of the major and minor axes of the prolate spheroid target, and an ion dose valley, staying always between the vertices, rather than at the vertex of the minor axis.

Role of cell deformability in the two-dimensional melting of biological tissues

Science.gov (United States)

Li, Yan-Wei; Ciamarra, Massimo Pica

2018-04-01

The size and shape of a large variety of polymeric particles, including biological cells, star polymers, dendrimes, and microgels, depend on the applied stresses as the particles are extremely soft. In high-density suspensions these particles deform as stressed by their neighbors, which implies that the interparticle interaction becomes of many-body type. Investigating a two-dimensional model of cell tissue, where the single particle shear modulus is related to the cell adhesion strength, here we show that the particle deformability affects the melting scenario. On increasing the temperature, stiff particles undergo a first-order solid/liquid transition, while soft ones undergo a continuous solid/hexatic transition followed by a discontinuous hexatic/liquid transition. At zero temperature the melting transition driven by the decrease of the adhesion strength occurs through two continuous transitions as in the Kosterlitz, Thouless, Halperin, Nelson, and Young scenario. Thus, there is a range of adhesion strength values where the hexatic phase is stable at zero temperature, which suggests that the intermediate phase of the epithelial-to-mesenchymal transition could be hexatic type.

Particle-in-Cell Modeling of Magnetized Argon Plasma Flow Through Small Mechanical Apertures

Energy Technology Data Exchange (ETDEWEB)

Adam B. Sefkow and Samuel A. Cohen

2009-04-09

Motivated by observations of supersonic argon-ion flow generated by linear helicon-heated plasma devices, a three-dimensional particle-in-cell (PIC) code is used to study whether stationary electrostatic layers form near mechanical apertures intersecting the flow of magnetized plasma. By self-consistently evaluating the temporal evolution of the plasma in the vicinity of the aperture, the PIC simulations characterize the roles of the imposed aperture and applied magnetic field on ion acceleration. The PIC model includes ionization of a background neutral-argon population by thermal and superthermal electrons, the latter found upstream of the aperture. Near the aperture, a transition from a collisional to a collisionless regime occurs. Perturbations of density and potential, with mm wavelengths and consistent with ion acoustic waves, propagate axially. An ion acceleration region of length ~ 200-300 λD,e forms at the location of the aperture and is found to be an electrostatic double layer, with axially-separated regions of net positive and negative charge. Reducing the aperture diameter or increasing its length increases the double layer strength.

Particle-in-Cell Modeling of Magnetized Argon Plasma Flow Through Small Mechanical Apertures

International Nuclear Information System (INIS)

Sefkow, Adam B.; Cohen, Samuel A.

2009-01-01

Motivated by observations of supersonic argon-ion flow generated by linear helicon-heated plasma devices, a three-dimensional particle-in-cell (PIC) code is used to study whether stationary electrostatic layers form near mechanical apertures intersecting the flow of magnetized plasma. By self-consistently evaluating the temporal evolution of the plasma in the vicinity of the aperture, the PIC simulations characterize the roles of the imposed aperture and applied magnetic field on ion acceleration. The PIC model includes ionization of a background neutral-argon population by thermal and superthermal electrons, the latter found upstream of the aperture. Near the aperture, a transition from a collisional to a collisionless regime occurs. Perturbations of density and potential, with mm wavelengths and consistent with ion acoustic waves, propagate axially. An ion acceleration region of length ∼ 200-300 λ D,e forms at the location of the aperture and is found to be an electrostatic double layer, with axially-separated regions of net positive and negative charge. Reducing the aperture diameter or increasing its length increases the double layer strength

Laser-plasma interactions with a Fourier-Bessel particle-in-cell method

Energy Technology Data Exchange (ETDEWEB)

Andriyash, Igor A., E-mail: igor.andriyash@gmail.com [Synchrotron SOLEIL, L' Orme des Merisiers, Saint Aubin, 91192 Gif-sur-Yvette (France); LOA, ENSTA ParisTech, CNRS, Ecole polytechnique, Université Paris-Saclay, 828 bd des Maréchaux, 91762 Palaiseau cedex (France); Lehe, Remi [Lawrence Berkeley National Laboratory, Berkeley, California 94720 (United States); Lifschitz, Agustin [LOA, ENSTA ParisTech, CNRS, Ecole polytechnique, Université Paris-Saclay, 828 bd des Maréchaux, 91762 Palaiseau cedex (France)

2016-03-15

A new spectral particle-in-cell (PIC) method for plasma modeling is presented and discussed. In the proposed scheme, the Fourier-Bessel transform is used to translate the Maxwell equations to the quasi-cylindrical spectral domain. In this domain, the equations are solved analytically in time, and the spatial derivatives are approximated with high accuracy. In contrast to the finite-difference time domain (FDTD) methods, that are used commonly in PIC, the developed method does not produce numerical dispersion and does not involve grid staggering for the electric and magnetic fields. These features are especially valuable in modeling the wakefield acceleration of particles in plasmas. The proposed algorithm is implemented in the code PLARES-PIC, and the test simulations of laser plasma interactions are compared to the ones done with the quasi-cylindrical FDTD PIC code CALDER-CIRC.

2D arc-PIC code description: methods and documentation

CERN Document Server

Timko, Helga

2011-01-01

Vacuum discharges are one of the main limiting factors for future linear collider designs such as that of the Compact LInear Collider. To optimize machine efficiency, maintaining the highest feasible accelerating gradient below a certain breakdown rate is desirable; understanding breakdowns can therefore help us to achieve this goal. As a part of ongoing theoretical research on vacuum discharges at the Helsinki Institute of Physics, the build-up of plasma can be investigated through the particle-in-cell method. For this purpose, we have developed the 2D Arc-PIC code introduced here. We present an exhaustive description of the 2D Arc-PIC code in two parts. In the first part, we introduce the particle-in-cell method in general and detail the techniques used in the code. In the second part, we provide a documentation and derivation of the key equations occurring in the code. The code is original work of the author, written in 2010, and is therefore under the copyright of the author. The development of the code h...

Electromagnetic direct implicit PIC simulation

International Nuclear Information System (INIS)

Langdon, A.B.

1983-01-01

Interesting modelling of intense electron flow has been done with implicit particle-in-cell simulation codes. In this report, the direct implicit PIC simulation approach is applied to simulations that include full electromagnetic fields. The resulting algorithm offers advantages relative to moment implicit electromagnetic algorithms and may help in our quest for robust and simpler implicit codes

PIC simulations of conical magnetically insulated transmission line with LTD generator: Transition from self-limited to load-limited flow

Science.gov (United States)

Liu, Laqun; Wang, Huihui; Guo, Fan; Zou, Wenkang; Liu, Dagang

2017-04-01

Based on the 3-dimensional Particle-In-Cell (PIC) code CHIPIC3D, with a new circuit boundary algorithm we developed, a conical magnetically insulated transmission line (MITL) with a 1.0-MV linear transformer driver (LTD) is explored numerically. The values of switch jitter time of LTD are critical parameters for the system, which are difficult to be measured experimentally. In this paper, these values are obtained by comparing the PIC results with experimental data of large diode-gap MITL. By decreasing the diode gap, we find that all PIC results agree well with experimental data only if MITL works on self-limited flow no matter how large the diode gap is. However, when the diode gap decreases to a threshold, the self-limited flow would transfer to a load-limited flow. In this situation, PIC results no longer agree with experimental data anymore due to the anode plasma expansion in the diode load. This disagreement is used to estimate the plasma expansion speed.

Multi-Band Light Curves from Two-Dimensional Simulations of Gamma-Ray Burst Afterglows

Science.gov (United States)

MacFadyen, Andrew

2010-01-01

The dynamics of gamma-ray burst outflows is inherently multi-dimensional. 1.) We present high resolution two-dimensional relativistic hydrodynamics simulations of GRBs in the afterglow phase using adaptive mesh refinement (AMR). Using standard synchrotron radiation models, we compute multi-band light curves, from the radio to X-ray, directly from the 2D hydrodynamics simulation data. We will present on-axis light curves for both constant density and wind media. We will also present off-axis light curves relevant for searches for orphan afterglows. We find that jet breaks are smoothed due to both off-axis viewing and wind media effects. 2.) Non-thermal radiation mechanisms in GRB afterglows require substantial magnetic field strengths. In turbulence driven by shear instabilities in relativistic magnetized gas, we demonstrate that magnetic field is naturally amplified to half a percent of the total energy (epsilon B = 0.005). We will show high resolution three dimensional relativistic MHD simulations of this process as well as particle in cell (PIC) simulations of mildly relativistic collisionless shocks.

A two-dimensional analytical model of laminar flame in lycopodium dust particles

Energy Technology Data Exchange (ETDEWEB)

Rahbari, Alireza [Shahid Rajaee Teacher Training University, Tehran (Iran, Islamic Republic of); Shakibi, Ashkan [Iran University of Science and Technology, Tehran (Iran, Islamic Republic of); Bidabadi, Mehdi [Combustion Research Laboratory, Narmak, Tehran (Iran, Islamic Republic of)

2015-09-15

A two-dimensional analytical model is presented to determine the flame speed and temperature distribution of micro-sized lycopodium dust particles. This model is based on the assumptions that the particle burning rate in the flame front is controlled by the process of oxygen diffusion and the flame structure consists of preheat, reaction and post flame zones. In the first step, the energy conservation equations for fuel-lean condition are expressed in two dimensions, and then these differential equations are solved using the required boundary condition and matching the temperature and heat flux at the interfacial boundaries. Consequently, the obtained flame temperature and flame speed distributions in terms of different particle diameters and equivalence ratio for lean mixture are compared with the corresponding experimental data for lycopodium dust particles. Consequently, it is shown that this two-dimensional model demonstrates better agreement with the experimental results compared to the previous models.

A two-dimensional analytical model of laminar flame in lycopodium dust particles

International Nuclear Information System (INIS)

Rahbari, Alireza; Shakibi, Ashkan; Bidabadi, Mehdi

2015-01-01

A two-dimensional analytical model is presented to determine the flame speed and temperature distribution of micro-sized lycopodium dust particles. This model is based on the assumptions that the particle burning rate in the flame front is controlled by the process of oxygen diffusion and the flame structure consists of preheat, reaction and post flame zones. In the first step, the energy conservation equations for fuel-lean condition are expressed in two dimensions, and then these differential equations are solved using the required boundary condition and matching the temperature and heat flux at the interfacial boundaries. Consequently, the obtained flame temperature and flame speed distributions in terms of different particle diameters and equivalence ratio for lean mixture are compared with the corresponding experimental data for lycopodium dust particles. Consequently, it is shown that this two-dimensional model demonstrates better agreement with the experimental results compared to the previous models.

Current-Sheet Formation and Reconnection at a Magnetic X Line in Particle-in-Cell Simulations

Science.gov (United States)

Black, C.; Antiochos, S. K.; Hesse, M.; Karpen, J. T.; Kuznetsova, M. M.; Zenitani, S.

2011-01-01

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the heliophysics community, particularly in the context of magnetic reconnection. Reconnection governs the large-scale energy release and topological rearrangement of magnetic fields in a wide variety of laboratory, heliophysical, and astrophysical systems. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high-resolution particle-in-cell (PIC) simulations. The initial minimum-energy, potential magnetic field is perturbed by excess thermal pressure introduced into the particle distribution function far from the X line. Subsequently, the relaxation of this added stress leads self-consistently to the development of a current sheet that reconnects for imposed stress of sufficient strength. We compare the time-dependent evolution and final state of our PIC simulations with macroscopic magnetohydrodynamic simulations assuming both uniform and localized electrical resistivities (C. R. DeVore et al., this meeting), as well as with force-free magnetic-field equilibria in which the amount of reconnection across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for understanding magnetic-reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the solar corona and terrestrial magnetotail.

Review of particle-in-cell modeling for the extraction region of large negative hydrogen ion sources for fusion

Science.gov (United States)

Wünderlich, D.; Mochalskyy, S.; Montellano, I. M.; Revel, A.

2018-05-01

Particle-in-cell (PIC) codes are used since the early 1960s for calculating self-consistently the motion of charged particles in plasmas, taking into account external electric and magnetic fields as well as the fields created by the particles itself. Due to the used very small time steps (in the order of the inverse plasma frequency) and mesh size, the computational requirements can be very high and they drastically increase with increasing plasma density and size of the calculation domain. Thus, usually small computational domains and/or reduced dimensionality are used. In the last years, the available central processing unit (CPU) power strongly increased. Together with a massive parallelization of the codes, it is now possible to describe in 3D the extraction of charged particles from a plasma, using calculation domains with an edge length of several centimeters, consisting of one extraction aperture, the plasma in direct vicinity of the aperture, and a part of the extraction system. Large negative hydrogen or deuterium ion sources are essential parts of the neutral beam injection (NBI) system in future fusion devices like the international fusion experiment ITER and the demonstration reactor (DEMO). For ITER NBI RF driven sources with a source area of 0.9 × 1.9 m2 and 1280 extraction apertures will be used. The extraction of negative ions is accompanied by the co-extraction of electrons which are deflected onto an electron dump. Typically, the maximum negative extracted ion current is limited by the amount and the temporal instability of the co-extracted electrons, especially for operation in deuterium. Different PIC codes are available for the extraction region of large driven negative ion sources for fusion. Additionally, some effort is ongoing in developing codes that describe in a simplified manner (coarser mesh or reduced dimensionality) the plasma of the whole ion source. The presentation first gives a brief overview of the current status of the ion

Wavelet-based Poisson Solver for use in Particle-In-Cell Simulations

International Nuclear Information System (INIS)

Terzic, B.; Mihalcea, D.; Bohn, C.L.; Pogorelov, I.V.

2005-01-01

We report on a successful implementation of a wavelet based Poisson solver for use in 3D particle-in-cell (PIC) simulations. One new aspect of our algorithm is its ability to treat the general(inhomogeneous) Dirichlet boundary conditions (BCs). The solver harnesses advantages afforded by the wavelet formulation, such as sparsity of operators and data sets, existence of effective preconditioners, and the ability simultaneously to remove numerical noise and further compress relevant data sets. Having tested our method as a stand-alone solver on two model problems, we merged it into IMPACT-T to obtain a fully functional serial PIC code. We present and discuss preliminary results of application of the new code to the modeling of the Fermilab/NICADD and AES/JLab photoinjectors

Thrust calculation of electric solar wind sail by particle-in-cell simulation

Energy Technology Data Exchange (ETDEWEB)

Hoshi, Kento [Kyoto Univ. (Japan). Dept. of Electrical Engineering; Kojima, Hirotsugu; Yamakawa, Hiroshi [Kyoto Univ. (Japan). Research Inst. for Sustainable Humanosphere; Muranaka, Takanobu [Chukyo Univ., Nagoya (Japan). Dept. of Electrical Engineering

2016-07-01

In this study, thrust characteristics of an electric solar wind sail were numerically evaluated using full threedimensional particle-in-cell (PIC) simulation. The thrust obtained from the PIC simulation was lower than the thrust estimations obtained in previous studies. The PIC simulation indicated that ambient electrons strongly shield the electrostatic potential of the tether of the sail, and the strong shield effect causes a greater thrust reduction than has been obtained in previous studies. Additionally, previous expressions of the thrust estimation were modified by using the shielded potential structure derived from the present simulation results. The modified thrust estimation agreed very well with the thrust obtained from the PIC simulation.

Thrust calculation of electric solar wind sail by particle-in-cell simulation

International Nuclear Information System (INIS)

Hoshi, Kento; Kojima, Hirotsugu; Yamakawa, Hiroshi; Muranaka, Takanobu

2016-01-01

In this study, thrust characteristics of an electric solar wind sail were numerically evaluated using full threedimensional particle-in-cell (PIC) simulation. The thrust obtained from the PIC simulation was lower than the thrust estimations obtained in previous studies. The PIC simulation indicated that ambient electrons strongly shield the electrostatic potential of the tether of the sail, and the strong shield effect causes a greater thrust reduction than has been obtained in previous studies. Additionally, previous expressions of the thrust estimation were modified by using the shielded potential structure derived from the present simulation results. The modified thrust estimation agreed very well with the thrust obtained from the PIC simulation.

Particle-in-cell Simulations with Kinetic Electrons

International Nuclear Information System (INIS)

Lewandowski, J.L.V.

2004-01-01

A new scheme, based on an exact separation between adiabatic and nonadiabatic electron responses, for particle-in-cell (PIC) simulations of drift-type modes is presented. The (linear and nonlinear) elliptic equations for the scalar fields are solved using a multi-grid solver. The new scheme yields linear growth rates in excellent agreement with theory and it is shown to conserve energy well into the nonlinear regime. It is also demonstrated that simulations with few electrons are reliable and accurate, suggesting that large-scale, PIC simulations with electron dynamics in toroidal geometry (e.g., tokamaks and stellarators plasmas) are within reach of present-day massively parallel supercomputers

Two-dimensional particle-in-cell simulation of the expansion of a plasma into a rarefied medium

Energy Technology Data Exchange (ETDEWEB)

Sarri, G; Quinn, K; Kourakis, I; Borghesi, M [Centre for Plasma Physics, The Queens University of Belfast, Belfast BT7 1NN (United Kingdom); Murphy, G C; Drury, L O C [Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2 (Ireland); Dieckmann, M E; Ynnerman, A [Department of Science and Technology (ITN), Linkoeping University, 60174 Norrkoping (Sweden); Bret, A, E-mail: gsarri01@qub.ac.uk [ETSI Industriales, Universidad de Castilla-La Mancha, 13071 Ciudad Real (Spain)

2011-07-15

The expansion of a dense plasma through a more rarefied ionized medium has been studied by means of two-dimensional particle-in-cell simulations. The initial conditions involve a density jump by a factor of 100, located in the middle of an otherwise equally dense electron-proton plasma with uniform proton and electron temperatures of 10 eV and 1 keV, respectively. Simulations show the creation of a purely electrostatic collisionless shock together with an ion-acoustic soliton tied to its downstream region. The shock front is seen to evolve in filamentary structures consistently with the onset of the ion-ion instability. Meanwhile, an un-magnetized drift instability is triggered in the core part of the dense plasma. Such results explain recent experimental laser-plasma experiments, carried out in similar conditions, and are of intrinsic relevance to non-relativistic shock scenarios in the solar and astrophysical systems.

Two-dimensional particle-in-cell simulation of the expansion of a plasma into a rarefied medium

International Nuclear Information System (INIS)

Sarri, G; Quinn, K; Kourakis, I; Borghesi, M; Murphy, G C; Drury, L O C; Dieckmann, M E; Ynnerman, A; Bret, A

2011-01-01

The expansion of a dense plasma through a more rarefied ionized medium has been studied by means of two-dimensional particle-in-cell simulations. The initial conditions involve a density jump by a factor of 100, located in the middle of an otherwise equally dense electron-proton plasma with uniform proton and electron temperatures of 10 eV and 1 keV, respectively. Simulations show the creation of a purely electrostatic collisionless shock together with an ion-acoustic soliton tied to its downstream region. The shock front is seen to evolve in filamentary structures consistently with the onset of the ion-ion instability. Meanwhile, an un-magnetized drift instability is triggered in the core part of the dense plasma. Such results explain recent experimental laser-plasma experiments, carried out in similar conditions, and are of intrinsic relevance to non-relativistic shock scenarios in the solar and astrophysical systems.

Enhanced stopping of macro-particles in particle-in-cell simulations

International Nuclear Information System (INIS)

May, J.; Tonge, J.; Ellis, I.; Mori, W. B.; Fiuza, F.; Fonseca, R. A.; Silva, L. O.; Ren, C.

2014-01-01

We derive an equation for energy transfer from relativistic charged particles to a cold background plasma appropriate for finite-size particles that are used in particle-in-cell simulation codes. Expressions for one-, two-, and three-dimensional particles are presented, with special attention given to the two-dimensional case. This energy transfer is due to the electric field of the wake set up in the background plasma by the relativistic particle. The enhanced stopping is dependent on the q 2 /m, where q is the charge and m is the mass of the relativistic particle, and therefore simulation macro-particles with large charge but identical q/m will stop more rapidly. The stopping power also depends on the effective particle shape of the macro-particle. These conclusions are verified in particle-in-cell simulations. We present 2D simulations of test particles, relaxation of high-energy tails, and integrated fast ignition simulations showing that the enhanced drag on macro-particles may adversely affect the results of these simulations in a wide range of high-energy density plasma scenarios. We also describe a particle splitting algorithm which can potentially overcome this problem and show its effect in controlling the stopping of macro-particles

Load management strategy for Particle-In-Cell simulations in high energy particle acceleration

Energy Technology Data Exchange (ETDEWEB)

Beck, A., E-mail: beck@llr.in2p3.fr [Laboratoire Leprince-Ringuet, École polytechnique, CNRS-IN2P3, Palaiseau 91128 (France); Frederiksen, J.T. [Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 København Ø (Denmark); Dérouillat, J. [CEA, Maison de La Simulation, 91400 Saclay (France)

2016-09-01

In the wake of the intense effort made for the experimental CILEX project, numerical simulation campaigns have been carried out in order to finalize the design of the facility and to identify optimal laser and plasma parameters. These simulations bring, of course, important insight into the fundamental physics at play. As a by-product, they also characterize the quality of our theoretical and numerical models. In this paper, we compare the results given by different codes and point out algorithmic limitations both in terms of physical accuracy and computational performances. These limitations are illustrated in the context of electron laser wakefield acceleration (LWFA). The main limitation we identify in state-of-the-art Particle-In-Cell (PIC) codes is computational load imbalance. We propose an innovative algorithm to deal with this specific issue as well as milestones towards a modern, accurate high-performance PIC code for high energy particle acceleration.

Electrostatic and electromagnetic instabilities associated with electrostatic shocks: Two-dimensional particle-in-cell simulation

International Nuclear Information System (INIS)

Kato, Tsunehiko N.; Takabe, Hideaki

2010-01-01

A two-dimensional electromagnetic particle-in-cell simulation with the realistic ion-to-electron mass ratio of 1836 is carried out to investigate the electrostatic collisionless shocks in relatively high-speed (∼3000 km s -1 ) plasma flows and also the influence of both electrostatic and electromagnetic instabilities, which can develop around the shocks, on the shock dynamics. It is shown that the electrostatic ion-ion instability can develop in front of the shocks, where the plasma is under counterstreaming condition, with highly oblique wave vectors as was shown previously. The electrostatic potential generated by the electrostatic ion-ion instability propagating obliquely to the shock surface becomes comparable with the shock potential and finally the shock structure is destroyed. It is also shown that in front of the shock the beam-Weibel instability gradually grows as well, consequently suggesting that the magnetic field generated by the beam-Weibel instability becomes important in long-term evolution of the shock and the Weibel-mediated shock forms long after the electrostatic shock vanished. It is also observed that the secondary electrostatic shock forms in the reflected ions in front of the primary electrostatic shock.

On the confinement of a Dirac particle to a two-dimensional ring

International Nuclear Information System (INIS)

Bakke, K.; Furtado, C.

2012-01-01

In this contribution, we propose a new model for studying the confinement of a spin-half particle to a two-dimensional quantum ring for systems described by the Dirac equation by introducing a new coupling into the Dirac equation. We show that the introduction of this new coupling into the Dirac equation yields a generalization of the two-dimensional quantum ring model proposed by Tan and Inkson [W.-C. Tan, J.C. Inkson, Semicond. Sci. Technol. 11 (1996) 1635] for relativistic spin-half quantum particles. -- Highlights: ► Two-dimensional ring model for condensed matter systems described by the Dirac equation. ► Exact solutions of the Dirac equation. ► Persistent currents for Dirac-like systems confined to a two-dimensional quantum ring.

Low frequency, electrodynamic simulation of kinetic plasmas with the DArwin Direct Implicit Particle-In-Cell (DADIPIC) method

International Nuclear Information System (INIS)

Gibbons, M.R.

1995-06-01

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. With the elimination of these constraints spatial and temporal discretization can be much larger than that possible with explicit, electrodynamic PIC. The code functions in a two dimensional Cartesian region and has been implemented with all components of the particle velocities, the E-field, and the B-field. Internal structures, conductors or dielectrics, may be placed in the simulation region, can be set at desired potentials, and driven with specified currents

The three-dimensional particle-in-cell simulation analysis of cavity of high power subterahertz pulsed gyrotron

International Nuclear Information System (INIS)

Ito, Koyu; Jiang, Weihua

2013-01-01

High power sub-terahertz pulsed gyrotrons for Collective Thomson Scattering (CTS) diagnostics of fusion plasmas are being developed. The typical target parameters are: output power of 100-200 kW, operation frequency of 300 GHz, and pulsed length > 10 us. In order to support experimental development, numerical simulations were carried out by using Particle-In-Cell (PIC) code MAGIC. The oscillation mode of the electromagnetic radiation was selected as TE_1_5_,_2, for which the beam parameters and cavity dimensions were determined accordingly. The simulation results have showed maximum power of 144 kW at oscillation frequency of 292.80 GHz, with oscillation efficiency of 22.15%. (author)

Coupling motion of colloidal particles in quasi-two-dimensional confinement

International Nuclear Information System (INIS)

Ma, Jun; Jing, Guangyin

2014-01-01

The Brownian motion of colloidal particles in quasi-two-dimensional (q2D) confinement displays a distinct kinetic character from that in bulk. Here we experimentally report dynamic coupling motion of Brownian particles in a relatively long process (∼100 h), which displays a quasi-equilibrium state in the q2D system. In the quasi-equilibrium state, the q2D confinement results in the coupling of particle motions, which slowly damps the motion and interaction of particles until the final equilibrium state is reached. The process of approaching the equilibrium is a random relaxation of a many-body interaction system of Brownian particles. As the relaxation proceeds for ∼100 h, the system reaches the equilibrium state in which the energy gained by the particles from the stochastic collision in the whole system is counteracted by the dissipative energy resulting from the collision. The relaxation time of this stochastic q2D system is 17.7 h. The theory is developed to explain coupling motions of Brownian particles in q2D confinement. (paper)

SELF-CONSISTENT LANGEVIN SIMULATION OF COULOMB COLLISIONS IN CHARGED-PARTICLE BEAMS

International Nuclear Information System (INIS)

QIANG, J.; RYNE, R.; HABIB, S.

2000-01-01

In many plasma physics and charged-particle beam dynamics problems, Coulomb collisions are modeled by a Fokker-Planck equation. In order to incorporate these collisions, we present a three-dimensional parallel Langevin simulation method using a Particle-In-Cell (PIC) approach implemented on high-performance parallel computers. We perform, for the first time, a fully self-consistent simulation, in which the FR-iction and diffusion coefficients are computed FR-om first principles. We employ a two-dimensional domain decomposition approach within a message passing programming paradigm along with dynamic load balancing. Object oriented programming is used to encapsulate details of the communication syntax as well as to enhance reusability and extensibility. Performance tests on the SGI Origin 2000 and the Cray T3E-900 have demonstrated good scalability. Work is in progress to apply our technique to intrabeam scattering in accelerators

Exactly energy conserving semi-implicit particle in cell formulation

International Nuclear Information System (INIS)

Lapenta, Giovanni

2017-01-01

We report a new particle in cell (PIC) method based on the semi-implicit approach. The novelty of the new method is that unlike any of its semi-implicit predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. Recent research has presented fully implicit methods where energy conservation is obtained as part of a non-linear iteration procedure. The new method (referred to as Energy Conserving Semi-Implicit Method, ECSIM), instead, does not require any non-linear iteration and its computational cycle is similar to that of explicit PIC. The properties of the new method are: i) it conserves energy exactly to round-off for any time step or grid spacing; ii) it is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency and allowing the user to select any desired time step; iii) it eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length; iv) the particle mover has a computational complexity identical to that of the explicit PIC, only the field solver has an increased computational cost. The new ECSIM is tested in a number of benchmarks where accuracy and computational performance are tested. - Highlights: • We present a new fully energy conserving semi-implicit particle in cell (PIC) method based on the implicit moment method (IMM). The new method is called Energy Conserving Implicit Moment Method (ECIMM). • The novelty of the new method is that unlike any of its predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. • The new method is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency. • The new method eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length. • These

Exactly energy conserving semi-implicit particle in cell formulation

Energy Technology Data Exchange (ETDEWEB)

Lapenta, Giovanni, E-mail: giovanni.lapenta@kuleuven.be

2017-04-01

We report a new particle in cell (PIC) method based on the semi-implicit approach. The novelty of the new method is that unlike any of its semi-implicit predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. Recent research has presented fully implicit methods where energy conservation is obtained as part of a non-linear iteration procedure. The new method (referred to as Energy Conserving Semi-Implicit Method, ECSIM), instead, does not require any non-linear iteration and its computational cycle is similar to that of explicit PIC. The properties of the new method are: i) it conserves energy exactly to round-off for any time step or grid spacing; ii) it is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency and allowing the user to select any desired time step; iii) it eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length; iv) the particle mover has a computational complexity identical to that of the explicit PIC, only the field solver has an increased computational cost. The new ECSIM is tested in a number of benchmarks where accuracy and computational performance are tested. - Highlights: • We present a new fully energy conserving semi-implicit particle in cell (PIC) method based on the implicit moment method (IMM). The new method is called Energy Conserving Implicit Moment Method (ECIMM). • The novelty of the new method is that unlike any of its predecessors at the same time it retains the explicit computational cycle and conserves energy exactly. • The new method is unconditionally stable in time, freeing the user from the need to resolve the electron plasma frequency. • The new method eliminates the constraint of the finite grid instability, allowing the user to select any desired resolution without being forced to resolve the Debye length. • These

Electrostatic plasma simulation by Particle-In-Cell method using ANACONDA package

International Nuclear Information System (INIS)

Blandón, J S; Grisales, J P; Riascos, H

2017-01-01

Electrostatic plasma is the most representative and basic case in plasma physics field. One of its main characteristics is its ideal behavior, since it is assumed be in thermal equilibrium state. Through this assumption, it is possible to study various complex phenomena such as plasma oscillations, waves, instabilities or damping. Likewise, computational simulation of this specific plasma is the first step to analyze physics mechanisms on plasmas, which are not at equilibrium state, and hence plasma is not ideal. Particle-In-Cell (PIC) method is widely used because of its precision for this kind of cases. This work, presents PIC method implementation to simulate electrostatic plasma by Python, using ANACONDA packages. The code has been corroborated comparing previous theoretical results for three specific phenomena in cold plasmas: oscillations, Two-Stream instability (TSI) and Landau Damping(LD). Finally, parameters and results are discussed. (paper)

NAMMA TWO-DIMENSIONAL STEREO PROBE AND CLOUD PARTICLE IMAGER V1

Data.gov (United States)

National Aeronautics and Space Administration — This Cloud Microphysics dataset consists of data from two probes used to measure the size, shape, and concentration of cloud particles; the two-dimensional stereo...

Development of in-situ visualization tool for PIC simulation

International Nuclear Information System (INIS)

Ohno, Nobuaki; Ohtani, Hiroaki

2014-01-01

As the capability of a supercomputer is improved, the sizes of simulation and its output data also become larger and larger. Visualization is usually carried out on a researcher's PC with interactive visualization software after performing the computer simulation. However, the data size is becoming too large to do it currently. A promising answer is in-situ visualization. For this case a simulation code is coupled with the visualization code and visualization is performed with the simulation on the same supercomputer. We developed an in-situ visualization tool for particle-in-cell (PIC) simulation and it is provided as a Fortran's module. We coupled it with a PIC simulation code and tested the coupled code on Plasma Simulator supercomputer, and ensured that it works. (author)

Numerical heating in Particle-In-Cell simulations with Monte Carlo binary collisions

Science.gov (United States)

Alves, E. Paulo; Mori, Warren; Fiuza, Frederico

2017-10-01

The binary Monte Carlo collision (BMCC) algorithm is a robust and popular method to include Coulomb collision effects in Particle-in-Cell (PIC) simulations of plasmas. While a number of works have focused on extending the validity of the model to different physical regimes of temperature and density, little attention has been given to the fundamental coupling between PIC and BMCC algorithms. Here, we show that the coupling between PIC and BMCC algorithms can give rise to (nonphysical) numerical heating of the system, that can be far greater than that observed when these algorithms operate independently. This deleterious numerical heating effect can significantly impact the evolution of the simulated system particularly for long simulation times. In this work, we describe the source of this numerical heating, and derive scaling laws for the numerical heating rates based on the numerical parameters of PIC-BMCC simulations. We compare our theoretical scalings with PIC-BMCC numerical experiments, and discuss strategies to minimize this parasitic effect. This work is supported by DOE FES under FWP 100237 and 100182.

Multi-grid Particle-in-cell Simulations of Plasma Microturbulence

International Nuclear Information System (INIS)

Lewandowski, J.L.V.

2003-01-01

A new scheme to accurately retain kinetic electron effects in particle-in-cell (PIC) simulations for the case of electrostatic drift waves is presented. The splitting scheme, which is based on exact separation between adiabatic and on adiabatic electron responses, is shown to yield more accurate linear growth rates than the standard df scheme. The linear and nonlinear elliptic problems that arise in the splitting scheme are solved using a multi-grid solver. The multi-grid particle-in-cell approach offers an attractive path, both from the physics and numerical points of view, to simulate kinetic electron dynamics in global toroidal plasmas

Nonequilibrium Gyrokinetic Fluctuation Theory and Sampling Noise in Gyrokinetic Particle-in-cell Simulations

International Nuclear Information System (INIS)

Krommes, John A.

2007-01-01

The present state of the theory of fluctuations in gyrokinetic (GK) plasmas and especially its application to sampling noise in GK particle-in-cell (PIC) simulations is reviewed. Topics addressed include the Δf method, the fluctuation-dissipation theorem for both classical and GK many-body plasmas, the Klimontovich formalism, sampling noise in PIC simulations, statistical closure for partial differential equations, the theoretical foundations of spectral balance in the presence of arbitrary noise sources, and the derivation of Kadomtsev-type equations from the general formalism

Nonequilibrium Gyrokinetic Fluctuation Theory and Sampling Noise in Gyrokinetic Particle-in-cell Simulations

Energy Technology Data Exchange (ETDEWEB)

John A. Krommes

2007-10-09

The present state of the theory of fluctuations in gyrokinetic GK plasmas and especially its application to sampling noise in GK particle-in-cell PIC simulations is reviewed. Topics addressed include the Δf method, the fluctuation-dissipation theorem for both classical and GK many-body plasmas, the Klimontovich formalism, sampling noise in PIC simulations, statistical closure for partial differential equations, the theoretical foundations of spectral balance in the presence of arbitrary noise sources, and the derivation of Kadomtsev-type equations from the general formalism.

Two-dimensional PIC simulations of ion beam instabilities in Supernova-driven plasma flows

Energy Technology Data Exchange (ETDEWEB)

Dieckmann, M E; Shukla, P K [Institut fuer Theoretische Physik IV, Ruhr-Universitaet Bochum, D-44780 Bochum (Germany); Meli, A; Mastichiadis, A [Department of Physics, National University of Athens, Panepistimiopolis, Zografos 15783 (Greece); Drury, L O C [Dublin Institute for Advanced Studies, Dublin 2 (Ireland)], E-mail: markd@tp4.rub.de

2008-06-15

Supernova remnant blast shells can reach the flow speed v{sub s} = 0.1c and shocks form at its front. Instabilities driven by shock-reflected ion beams heat the plasma in the foreshock, which may inject particles into diffusive acceleration. The ion beams can have the speed v{sub b} {approx} v{sub s}. For v{sub b} << v{sub s} the Buneman or upper-hybrid instabilities dominate, while for v{sub b} >> v{sub s} the filamentation and mixed modes grow faster. Here the relevant waves for v{sub b} {approx} v{sub s} are examined and how they interact nonlinearly with the particles. The collision of two plasma clouds at the speed v{sub s} is modelled with particle-in-cell simulations, which convect with them magnetic fields oriented perpendicular to their flow velocity vector. One simulation models equally dense clouds and the other one uses a density ratio of 2. Both simulations show upper-hybrid waves that are planar over large spatial intervals and that accelerate electrons to {approx}10 keV. The symmetric collision yields only short oscillatory wave pulses, while the asymmetric collision also produces large-scale electric fields, probably through a magnetic pressure gradient. The large-scale fields destroy the electron phase space holes and they accelerate the ions, which facilitates the formation of a precursor shock.

Fusion PIC code performance analysis on the Cori KNL system

Energy Technology Data Exchange (ETDEWEB)

Koskela, Tuomas S. [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Deslippe, Jack [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Friesen, Brian [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Raman, Karthic [INTEL Corp. (United States)

2017-05-25

We study the attainable performance of Particle-In-Cell codes on the Cori KNL system by analyzing a miniature particle push application based on the fusion PIC code XGC1. We start from the most basic building blocks of a PIC code and build up the complexity to identify the kernels that cost the most in performance and focus optimization efforts there. Particle push kernels operate at high AI and are not likely to be memory bandwidth or even cache bandwidth bound on KNL. Therefore, we see only minor benefits from the high bandwidth memory available on KNL, and achieving good vectorization is shown to be the most beneficial optimization path with theoretical yield of up to 8x speedup on KNL. In practice we are able to obtain up to a 4x gain from vectorization due to limitations set by the data layout and memory latency.

A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions

Science.gov (United States)

Chacón, L.; Chen, G.

2016-07-01

We extend a recently proposed fully implicit PIC algorithm for the Vlasov-Darwin model in multiple dimensions (Chen and Chacón (2015) [1]) to curvilinear geometry. As in the Cartesian case, the approach is based on a potential formulation (ϕ, A), and overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. Conservation theorems for local charge and global energy are derived in curvilinear representation, and then enforced discretely by a careful choice of the discretization of field and particle equations. Additionally, the algorithm conserves canonical-momentum in any ignorable direction, and preserves the Coulomb gauge ∇ ṡ A = 0 exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. We demonstrate the accuracy and efficiency properties of the algorithm with numerical experiments in mapped meshes in 1D-3V and 2D-3V.

Two-dimensional position sensitive silicon photodiode as a charged particle detector

International Nuclear Information System (INIS)

Kovacevic, K.; Zadro, M.

1999-01-01

A two-dimensional position sensitive silicon photodiode has been tested for measurement of position and energy of charged particles. Position nonlinearity and resolution, as well as energy resolution and ballistic deficit were measured for 5.486 MeV α-particles. The results obtained for different pulse shaping time constants are presented

Quantum mechanical treatment of a constrained particle on two dimensional sphere

Energy Technology Data Exchange (ETDEWEB)

Jahangiri, L., E-mail: laleh.jahangiry@yahoo.com; Panahi, H., E-mail: t-panahi@guilan.ac.ir

2016-12-15

In this work, we study the motion of a particle on two dimensional sphere. By writing the Schrodinger equation, we obtain the wave function and energy spectra for three dimensional harmonic oscillator potential plus trigonometric Rosen–Morse non-central potential. By letting three special cases for intertwining operator, we investigate the energy spectra and wave functions for Smorodinsky–Winternitz potential model.

The Particle-in-Cell and Kinetic Simulation Software Center

Science.gov (United States)

Mori, W. B.; Decyk, V. K.; Tableman, A.; Fonseca, R. A.; Tsung, F. S.; Hu, Q.; Winjum, B. J.; An, W.; Dalichaouch, T. N.; Davidson, A.; Hildebrand, L.; Joglekar, A.; May, J.; Miller, K.; Touati, M.; Xu, X. L.

2017-10-01

The UCLA Particle-in-Cell and Kinetic Simulation Software Center (PICKSC) aims to support an international community of PIC and plasma kinetic software developers, users, and educators; to increase the use of this software for accelerating the rate of scientific discovery; and to be a repository of knowledge and history for PIC. We discuss progress towards making available and documenting illustrative open-source software programs and distinct production programs; developing and comparing different PIC algorithms; coordinating the development of resources for the educational use of kinetic software; and the outcomes of our first sponsored OSIRIS users workshop. We also welcome input and discussion from anyone interested in using or developing kinetic software, in obtaining access to our codes, in collaborating, in sharing their own software, or in commenting on how PICKSC can better serve the DPP community. Supported by NSF under Grant ACI-1339893 and by the UCLA Institute for Digital Research and Education.

Activation of AMP-Activated Protein Kinase α and Extracelluar Signal-Regulated Kinase Mediates CB-PIC-Induced Apoptosis in Hypoxic SW620 Colorectal Cancer Cells

Directory of Open Access Journals (Sweden)

Sung-Yun Cho

2013-01-01

Full Text Available Here, antitumor mechanism of cinnamaldehyde derivative CB-PIC was elucidated in human SW620 colon cancer cells. CB-PIC significantly exerted cytotoxicity, increased sub-G1 accumulation, and cleaved PARP with apoptotic features, while it enhanced the phosphorylation of AMPK alpha and ACC as well as activated the ERK in hypoxic SW620 cells. Furthermore, CB-PIC suppressed the expression of HIF1 alpha, Akt, and mTOR and activated the AMPK phosphorylation in hypoxic SW620 cells. Conversely, silencing of AMPKα blocked PARP cleavage and ERK activation induced by CB-PIC, while ERK inhibitor PD 98059 attenuated the phosphorylation of AMPKα in hypoxic SW620 cells, implying cross-talk between ERK and AMPKα. Furthermore, cotreatment of CB-PIC and metformin enhanced the inhibition of HIF1α and Akt/mTOR and the activation of AMPKα and pACC in hypoxic SW620 cells. In addition, CB-PIC suppressed the growth of SW620 cells inoculated in BALB/c athymic nude mice, and immunohistochemistry revealed that CB-PIC treatment attenuated the expression of Ki-67, CD34, and CAIX and increased the expression of pAMPKα in CB-PIC-treated group. Interestingly, CP-PIC showed better antitumor activity in SW620 colon cancer cells under hypoxia than under normoxia, since it may be applied to chemoresistance. Overall, our findings suggest that activation of AMPKα and ERK mediates CB-PIC-induced apoptosis in hypoxic SW620 colon cancer cells.

Concurrent particle-in-cell plasma simulation on a multi-transputer parallel computer

International Nuclear Information System (INIS)

Khare, A.N.; Jethra, A.; Patel, Kartik

1992-01-01

This report describes the parallelization of a Particle-in-Cell (PIC) plasma simulation code on a multi-transputer parallel computer. The algorithm used in the parallelization of the PIC method is described. The decomposition schemes related to the distribution of the particles among the processors are discussed. The implementation of the algorithm on a transputer network connected as a torus is presented. The solutions of the problems related to global communication of data are presented in the form of a set of generalized communication functions. The performance of the program as a function of data size and the number of transputers show that the implementation is scalable and represents an effective way of achieving high performance at acceptable cost. (author). 11 refs., 4 figs., 2 tabs., appendices

QUICKSILVER - A general tool for electromagnetic PIC simulation

International Nuclear Information System (INIS)

Seidel, David B.; Coats, Rebecca S.; Johnson, William A.; Kiefer, Mark L.; Mix, L. Paul; Pasik, Michael F.; Pointon, Timothy D.; Quintenz, Jeffrey P.; Riley, Douglas J.; Turner, C. David

1997-01-01

The dramatic increase in computational capability that has occurred over the last ten years has allowed fully electromagnetic simulations of large, complex, three-dimensional systems to move progressively from impractical, to expensive, and recently, to routine and widespread. This is particularly true for systems that require the motion of free charge to be self-consistently treated. The QUICKSILVER electromagnetic Particle-In-Cell (EM-PIC) code has been developed at Sandia National Laboratories to provide a general tool to simulate a wide variety of such systems. This tool has found widespread use for many diverse applications, including high-current electron and ion diodes, magnetically insulated power transmission systems, high-power microwave oscillators, high-frequency digital and analog integrated circuit packages, microwave integrated circuit components, antenna systems, radar cross-section applications, and electromagnetic interaction with biological material. This paper will give a brief overview of QUICKSILVER and provide some thoughts on its future development

A particle-in-cell method for modeling small angle Coulomb collisions in plasmas

International Nuclear Information System (INIS)

Parker, S.E.

1989-01-01

We propose a computational method to self-consistently model small angle collisional effects. This method may be added to standard Particle-In-Cell (PIC) plasma simulations to include collisions, or as an alternative to solving the Fokker-Planck (FP) equation using finite difference methods. The distribution function is represented by a large number of particles. The particle velocities change due to the drag force, and the diffusion in velocity is represented by a random process. This is similar to previous Monte-Carlo methods except we calculate the drag force and diffusion tensor self- consistently. The particles are weighted to a grid in velocity space and associated ''Poisson equations'' are solved for the Rosenbluth potentials. The motivation is to avoid the very time consuming method of Coulomb scattering pair by pair. First the approximation for small angle Coulomb collisions is discussed. Next, the FP-PIC collision method is outlined. Then we show a test of the particle advance modeling an electron beam scattering off a fixed ion background. 4 refs

Performance of particle in cell methods on highly concurrent computational architectures

International Nuclear Information System (INIS)

Adams, M.F.; Ethier, S.; Wichmann, N.

2009-01-01

Particle in cell (PIC) methods are effective in computing Vlasov-Poisson system of equations used in simulations of magnetic fusion plasmas. PIC methods use grid based computations, for solving Poisson's equation or more generally Maxwell's equations, as well as Monte-Carlo type methods to sample the Vlasov equation. The presence of two types of discretizations, deterministic field solves and Monte-Carlo methods for the Vlasov equation, pose challenges in understanding and optimizing performance on today large scale computers which require high levels of concurrency. These challenges arises from the need to optimize two very different types of processes and the interactions between them. Modern cache based high-end computers have very deep memory hierarchies and high degrees of concurrency which must be utilized effectively to achieve good performance. The effective use of these machines requires maximizing concurrency by eliminating serial or redundant work and minimizing global communication. A related issue is minimizing the memory traffic between levels of the memory hierarchy because performance is often limited by the bandwidths and latencies of the memory system. This paper discusses some of the performance issues, particularly in regard to parallelism, of PIC methods. The gyrokinetic toroidal code (GTC) is used for these studies and a new radial grid decomposition is presented and evaluated. Scaling of the code is demonstrated on ITER sized plasmas with up to 16K Cray XT3/4 cores.

Performance of particle in cell methods on highly concurrent computational architectures

International Nuclear Information System (INIS)

Adams, M F; Ethier, S; Wichmann, N

2007-01-01

Particle in cell (PIC) methods are effective in computing Vlasov-Poisson system of equations used in simulations of magnetic fusion plasmas. PIC methods use grid based computations, for solving Poisson's equation or more generally Maxwell's equations, as well as Monte-Carlo type methods to sample the Vlasov equation. The presence of two types of discretizations, deterministic field solves and Monte-Carlo methods for the Vlasov equation, pose challenges in understanding and optimizing performance on today large scale computers which require high levels of concurrency. These challenges arises from the need to optimize two very different types of processes and the interactions between them. Modern cache based high-end computers have very deep memory hierarchies and high degrees of concurrency which must be utilized effectively to achieve good performance. The effective use of these machines requires maximizing concurrency by eliminating serial or redundant work and minimizing global communication. A related issue is minimizing the memory traffic between levels of the memory hierarchy because performance is often limited by the bandwidths and latencies of the memory system. This paper discusses some of the performance issues, particularly in regard to parallelism, of PIC methods. The gyrokinetic toroidal code (GTC) is used for these studies and a new radial grid decomposition is presented and evaluated. Scaling of the code is demonstrated on ITER sized plasmas with up to 16K Cray XT3/4 cores

Energy transport in a shear flow of particles in a two-dimensional dusty plasma.

Science.gov (United States)

Feng, Yan; Goree, J; Liu, Bin

2012-11-01

A shear flow of particles in a laser-driven two-dimensional (2D) dusty plasma is observed in a study of viscous heating and thermal conduction. Video imaging and particle tracking yields particle velocity data, which we convert into continuum data, presented as three spatial profiles: mean particle velocity (i.e., flow velocity), mean-square particle velocity, and mean-square fluctuations of particle velocity. These profiles and their derivatives allow a spatially resolved determination of each term in the energy and momentum continuity equations, which we use for two purposes. First, by balancing these terms so that their sum (i.e., residual) is minimized while varying viscosity η and thermal conductivity κ as free parameters, we simultaneously obtain values for η and κ in the same experiment. Second, by comparing the viscous heating and thermal conduction terms, we obtain a spatially resolved characterization of the viscous heating.

Electromagnetic ''particle-in-cell'' plasma simulation

International Nuclear Information System (INIS)

Langdon, A.B.

1985-01-01

''PIC'' simulation tracks particles through electromagnetic fields calculated self-consistently from the charge and current densities of the particles themselves, external sources, and boundaries. Already used extensively in plasma physics, such simulations have become useful in the design of accelerators and their r.f. sources. 5 refs

Electromagnetic particle-in-cell simulations of the solar wind interaction with lunar magnetic anomalies.

Science.gov (United States)

Deca, J; Divin, A; Lapenta, G; Lembège, B; Markidis, S; Horányi, M

2014-04-18

We present the first three-dimensional fully kinetic and electromagnetic simulations of the solar wind interaction with lunar crustal magnetic anomalies (LMAs). Using the implicit particle-in-cell code iPic3D, we confirm that LMAs may indeed be strong enough to stand off the solar wind from directly impacting the lunar surface forming a mini-magnetosphere, as suggested by spacecraft observations and theory. In contrast to earlier magnetohydrodynamics and hybrid simulations, the fully kinetic nature of iPic3D allows us to investigate the space charge effects and in particular the electron dynamics dominating the near-surface lunar plasma environment. We describe for the first time the interaction of a dipole model centered just below the lunar surface under plasma conditions such that only the electron population is magnetized. The fully kinetic treatment identifies electromagnetic modes that alter the magnetic field at scales determined by the electron physics. Driven by strong pressure anisotropies, the mini-magnetosphere is unstable over time, leading to only temporal shielding of the surface underneath. Future human exploration as well as lunar science in general therefore hinges on a better understanding of LMAs.

Resonant and non-resonant whistlers-particle interaction in the radiation belts

NARCIS (Netherlands)

E. Camporeale (Enrico)

2015-01-01

htmlabstractWe study the wave-particle interactions between lower band chorus whistlers and an anisotropic tenuous population of relativistic electrons. We present the first direct comparison of first-principle Particle-in-Cell (PIC) simulations with a quasi-linear diffusion code. In the PIC

Loading relativistic Maxwell distributions in particle simulations

Science.gov (United States)

Zenitani, S.

2015-12-01

In order to study energetic plasma phenomena by using particle-in-cell (PIC) and Monte-Carlo simulations, we need to deal with relativistic velocity distributions in these simulations. However, numerical algorithms to deal with relativistic distributions are not well known. In this contribution, we overview basic algorithms to load relativistic Maxwell distributions in PIC and Monte-Carlo simulations. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are newly proposed in a physically transparent manner. Their acceptance efficiencies are 􏰅50% for generic cases and 100% for symmetric distributions. They can be combined with arbitrary base algorithms.

Continuation of full-scale three-dimensional numerical experiments on high-intensity particle and laser beam-matter interactions

Energy Technology Data Exchange (ETDEWEB)

Mori, Warren, B.

2012-12-01

We present results from the grant entitled, Continuation of full-scale three-dimensional numerical experiments on high-intensity particle and laser beam-matter interactions. The research significantly advanced the understanding of basic high-energy density science (HEDS) on ultra intense laser and particle beam plasma interactions. This advancement in understanding was then used to to aid in the quest to make 1 GeV to 500 GeV plasma based accelerator stages. The work blended basic research with three-dimensions fully nonlinear and fully kinetic simulations including full-scale modeling of ongoing or planned experiments. The primary tool was three-dimensional particle-in-cell simulations. The simulations provided a test bed for theoretical ideas and models as well as a method to guide experiments. The research also included careful benchmarking of codes against experiment. High-fidelity full-scale modeling provided a means to extrapolate parameters into regimes that were not accessible to current or near term experiments, thereby allowing concepts to be tested with confidence before tens to hundreds of millions of dollars were spent building facilities. The research allowed the development of a hierarchy of PIC codes and diagnostics that is one of the most advanced in the world.

Particle-in-cell plasma simulation codes on the connection machine

International Nuclear Information System (INIS)

Walker, D.W.

1991-01-01

Methods for implementing three-dimensional, electromagnetic, relativistic PIC plasma simulation codes on the Connection Machine (CM-2) are discussed. The gather and scatter phases of the PIC algorithm involve indirect indexing of data, which results in a large amount of communication on the CM-2. Different data decompositions are described that seek to reduce the amount of communication while maintaining good load balance. These methods require the particles to be spatially sorted at the start of each time step, which introduced another form of overhead. The different methods are implemented in CM Fortran on the CM-2 and compared. It was found that the general router is slow in performing the communication in the gather and scatter steps, which precludes an efficient CM Fortran implementation. An alternative method that uses PARIS calls and the NEWS communication network to pipeline data along the axes of the VP set is suggested as a more efficient algorithm

Particle-in-cell simulations of fast magnetic field penetration into plasmas due to the Hall electric field

International Nuclear Information System (INIS)

Swanekamp, S.B.; Grossmann, J.M.; Fruchtman, A.; Oliver, B.V.; Ottinger, P.F.

1996-01-01

Particle-in-cell (PIC) simulations are used to study the penetration of magnetic field into plasmas in the electron-magnetohydrodynamic (EMHD) regime. These simulations represent the first definitive verification of EMHD with a PIC code. When ions are immobile, the PIC results reproduce many aspects of fluid treatments of the problem. However, the PIC results show a speed of penetration that is between 10% and 50% slower than predicted by one-dimensional fluid treatments. In addition, the PIC simulations show the formation of vortices in the electron flow behind the EMHD shock front. The size of these vortices is on the order of the collisionless electron skin depth and is closely coupled to the effects of electron inertia. An energy analysis shows that one-half the energy entering the plasma is stored as magnetic field energy while the other half is shared between internal plasma energy (thermal motion and electron vortices) and electron kinetic energy loss from the volume to the boundaries. The amount of internal plasma energy saturates after an initial transient phase so that late in time the rate that magnetic energy increases in the plasma is the same as the rate at which kinetic energy flows out through the boundaries. When ions are mobile it is observed that axial magnetic field penetration is followed by localized thinning in the ion density. The density thinning is produced by the large electrostatic fields that exist inside the electron vortices which act to reduce the space-charge imbalance necessary to support the vortices. This mechanism may play a role during the opening process of a plasma opening switch. copyright 1996 American Institute of Physics

The Particle inside a Ring: A Two-Dimensional Quantum Problem Visualized by Scanning Tunneling Microscopy

Science.gov (United States)

Ellison, Mark D.

2008-01-01

The one-dimensional particle-in-a-box model used to introduce quantum mechanics to students suffers from a tenuous connection to a real physical system. This article presents a two-dimensional model, the particle confined within a ring, that directly corresponds to observations of surface electrons in a metal trapped inside a circular barrier.…

A three-dimensional electrostatic particle-in-cell methodology on unstructured Delaunay-Voronoi grids

International Nuclear Information System (INIS)

Gatsonis, Nikolaos A.; Spirkin, Anton

2009-01-01

The mathematical formulation and computational implementation of a three-dimensional particle-in-cell methodology on unstructured Delaunay-Voronoi tetrahedral grids is presented. The method allows simulation of plasmas in complex domains and incorporates the duality of the Delaunay-Voronoi in all aspects of the particle-in-cell cycle. Charge assignment and field interpolation weighting schemes of zero- and first-order are formulated based on the theory of long-range constraints. Electric potential and fields are derived from a finite-volume formulation of Gauss' law using the Voronoi-Delaunay dual. Boundary conditions and the algorithms for injection, particle loading, particle motion, and particle tracking are implemented for unstructured Delaunay grids. Error and sensitivity analysis examines the effects of particles/cell, grid scaling, and timestep on the numerical heating, the slowing-down time, and the deflection times. The problem of current collection by cylindrical Langmuir probes in collisionless plasmas is used for validation. Numerical results compare favorably with previous numerical and analytical solutions for a wide range of probe radius to Debye length ratios, probe potentials, and electron to ion temperature ratios. The versatility of the methodology is demonstrated with the simulation of a complex plasma microsensor, a directional micro-retarding potential analyzer that includes a low transparency micro-grid.

Parallel Finite Element Particle-In-Cell Code for Simulations of Space-charge Dominated Beam-Cavity Interactions

International Nuclear Information System (INIS)

Candel, A.; Kabel, A.; Ko, K.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Prudencio, E.; Schussman, G.; Uplenchwar, R.

2007-01-01

Over the past years, SLAC's Advanced Computations Department (ACD) has developed the parallel finite element (FE) particle-in-cell code Pic3P (Pic2P) for simulations of beam-cavity interactions dominated by space-charge effects. As opposed to standard space-charge dominated beam transport codes, which are based on the electrostatic approximation, Pic3P (Pic2P) includes space-charge, retardation and boundary effects as it self-consistently solves the complete set of Maxwell-Lorentz equations using higher-order FE methods on conformal meshes. Use of efficient, large-scale parallel processing allows for the modeling of photoinjectors with unprecedented accuracy, aiding the design and operation of the next-generation of accelerator facilities. Applications to the Linac Coherent Light Source (LCLS) RF gun are presented

Plume expansion of a laser-induced plasma studied with the particle-in-cell method

DEFF Research Database (Denmark)

Ellegaard, O.; Nedelea, T.; Schou, Jørgen

2002-01-01

energy as well as electron energy. We have estimated the time constant for energy transfer between the electrons and the ions. The scaling of these processes is given by a single parameter determined by the Debye length obtained from the electron density in the plasma outside the surface. (C) 2002......The initial stage of laser-induced plasma plume expansion from a solid in vacuum and the effect of the Coulomb field have been studied. We have performed a one-dimensional numerical calculation by mapping the charge on a computational grid according to the particle-in-cell (PIC) method of Birdsall...... et al. It is assumed that the particle ablation from a surface with a fixed temperature takes place as a pulse, i.e. within a finite period of time. A number of characteristic quantities for the plasma plume are compared with similar data for expansion of neutrals as well as fluid models: Density...

New apparatus for direct counting of β particles from two-dimensional gels and an application to changes in protein synthesis due to cell density

International Nuclear Information System (INIS)

Anderson, H.L.; Puck, T.T.; Shera, E.B.

1987-07-01

A new method is described for scanning two-dimensional gels by the direct counting of β particles instead of autoradiography. The methodology is described; results are compared with autoradiographic results; and data are presented demonstrating changed patterns of protein synthesis accompanying changes in cell density. The method is rapid and permits identification of differences in protein abundance of approximately 10% for a substantial fraction of the more prominent proteins. A modulation effect of more than 5 standard deviations, accompanying contact inhibition of cell growth, is shown to occur for an appreciable number of these proteins. The method promises to be applicable to a variety of biochemical and genetic experiments designed to delineate changes in protein synthesis accompanying changes in genome, molecular environment, history, and state of differentiation of the cell populations studied. 13 refs., 8 figs., 4 tabs

Energy dispersion of charged particles decelerated in a two-dimensional electrostatic field of the type x1/n

International Nuclear Information System (INIS)

Zashkvara, V.V.; Bok, A.A.

1992-01-01

Two components of the spatial dispersion of particles with respect to kinetic energy can be distinguished of the motion of charged particle beams in electrostatic mirros with a two-dimensional field φ(x,y) ans xz symmetry plane. The first is the longitudinal dispersion, which is along the z axis perpendicular to the field; the second is the transverse dispersion, along the x axis parallel to the field vector in the plane of symmetry. The longitudinal dispersion is a basic characteristic of electrostatic mirrors used as energy analyzers. It has been shown that for first-order angular focusing, the longitudinal dispersion, divided by the focal length, is independent of the structure of the two-dimensional field and is a function only of the angle at which the charged particle beam enters the mirror. The transverse dispersion stems from the energy dependence of the penetration depth of the beam as it is decelerated, and it plays an important role when the energy of a charged particle beam is analyzed by the filtering principle, making use of the property of an electrostatic mirror to transmit or reflect charged particles with kinetic energy in a specified interval. This type of dispersion in electrostatic mirrors with two-dimensional fields has not been analyzed systematically. In the present note the authors consider a particular type of two-dimensional electrostatic field which is characterized by a large transverse dispersion, many times larger than in existing electrostatic reflecting filters employing planar and cylindrical fields

On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

International Nuclear Information System (INIS)

Meyers, M.D.; Huang, C.-K.; Zeng, Y.; Yi, S.A.; Albright, B.J.

2015-01-01

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTD scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models

On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

Science.gov (United States)

Meyers, M. D.; Huang, C.-K.; Zeng, Y.; Yi, S. A.; Albright, B. J.

2015-09-01

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTD scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models.

PIC Simulations of Hypersonic Plasma Instabilities

Science.gov (United States)

Niehoff, D.; Ashour-Abdalla, M.; Niemann, C.; Decyk, V.; Schriver, D.; Clark, E.

2013-12-01

The plasma sheaths formed around hypersonic aircraft (Mach number, M > 10) are relatively unexplored and of interest today to both further the development of new technologies and solve long-standing engineering problems. Both laboratory experiments and analytical/numerical modeling are required to advance the understanding of these systems; it is advantageous to perform these tasks in tandem. There has already been some work done to study these plasmas by experiments that create a rapidly expanding plasma through ablation of a target with a laser. In combination with a preformed magnetic field, this configuration leads to a magnetic "bubble" formed behind the front as particles travel at about Mach 30 away from the target. Furthermore, the experiment was able to show the generation of fast electrons which could be due to instabilities on electron scales. To explore this, future experiments will have more accurate diagnostics capable of observing time- and length-scales below typical ion scales, but simulations are a useful tool to explore these plasma conditions theoretically. Particle in Cell (PIC) simulations are necessary when phenomena are expected to be observed at these scales, and also have the advantage of being fully kinetic with no fluid approximations. However, if the scales of the problem are not significantly below the ion scales, then the initialization of the PIC simulation must be very carefully engineered to avoid unnecessary computation and to select the minimum window where structures of interest can be studied. One method of doing this is to seed the simulation with either experiment or ion-scale simulation results. Previous experiments suggest that a useful configuration for studying hypersonic plasma configurations is a ring of particles rapidly expanding transverse to an external magnetic field, which has been simulated on the ion scale with an ion-hybrid code. This suggests that the PIC simulation should have an equivalent configuration

One-dimensional hybrid-direct kinetic simulation of the discharge plasma in a Hall thruster

International Nuclear Information System (INIS)

Hara, Kentaro; Boyd, Iain D.; Kolobov, Vladimir I.

2012-01-01

In order to model the non-equilibrium plasma within the discharge region of a Hall thruster, the velocity distribution functions (VDFs) must be obtained accurately. A direct kinetic (DK) simulation method that directly solves the plasma Boltzmann equation can achieve better resolution of VDFs in comparison to particle simulations, such as the particle-in-cell (PIC) method that inherently include statistical noise. In this paper, a one-dimensional hybrid-DK simulation, which uses a DK simulation for heavy species and a fluid model for electrons, is developed and compared to a hybrid-PIC simulation. Time-averaged results obtained from the hybrid-DK simulation are in good agreement with hybrid-PIC results and experimental data. It is shown from a comparison of using a kinetic simulation and solving the continuity equation that modeling of the neutral atoms plays an important role for simulations of the Hall thruster discharge plasma. In addition, low and high frequency plasma oscillations are observed. Although the kinetic nature of electrons is not resolved due to the use of a fluid model, the hybrid-DK model provides spatially and temporally well-resolved plasma properties and an improved resolution of VDFs for heavy species with less statistical noise in comparison to the hybrid-PIC method.

An exploratory study of three-dimensional MP-PIC-based simulation of bubbling fluidized beds with and without baffles

DEFF Research Database (Denmark)

Yang, Shuai; Wu, Hao; Lin, Weigang

2018-01-01

In this study, the flow characteristics of Geldart A particles in a bubbling fluidized bed with and without perforated plates were simulated by the multiphase particle-in-cell (MP-PIC)-based Eulerian-Lagrangian method. A modified structure-based drag model was developed based on our previous work....... Other drag models including the Parker and Wen-Yu-Ergun drag models were also employed to investigate the effects of drag models on the simulation results. Although the modified structure-based drag model better predicts the gas-solid flow dynamics of a baffle-free bubbling fluidized bed in comparison...... with the experimental data, none of these drag models predict the gas-solid flow in a baffled bubbling fluidized bed sufficiently well because of the treatment of baffles in the Barracuda software. To improve the simulation accuracy, future versions of Barracuda should address the challenges of incorporating the bed...

Swelling of two-dimensional polymer rings by trapped particles.

Science.gov (United States)

Haleva, E; Diamant, H

2006-09-01

The mean area of a two-dimensional Gaussian ring of N monomers is known to diverge when the ring is subject to a critical pressure differential, p c ~ N -1. In a recent publication (Eur. Phys. J. E 19, 461 (2006)) we have shown that for an inextensible freely jointed ring this divergence turns into a second-order transition from a crumpled state, where the mean area scales as [A]~N-1, to a smooth state with [A]~N(2). In the current work we extend these two models to the case where the swelling of the ring is caused by trapped ideal-gas particles. The Gaussian model is solved exactly, and the freely jointed one is treated using a Flory argument, mean-field theory, and Monte Carlo simulations. For a fixed number Q of trapped particles the criticality disappears in both models through an unusual mechanism, arising from the absence of an area constraint. In the Gaussian case the ring swells to such a mean area, [A]~ NQ, that the pressure exerted by the particles is at p c for any Q. In the freely jointed model the mean area is such that the particle pressure is always higher than p c, and [A] consequently follows a single scaling law, [A]~N(2) f (Q/N), for any Q. By contrast, when the particles are in contact with a reservoir of fixed chemical potential, the criticality is retained. Thus, the two ensembles are manifestly inequivalent in these systems.

Quasi-two-dimensional complex plasma containing spherical particles and their binary agglomerates.

Science.gov (United States)

Chaudhuri, M; Semenov, I; Nosenko, V; Thomas, H M

2016-05-01

A unique type of quasi-two-dimensional complex plasma system was observed which consisted of monodisperse microspheres and their binary agglomerations (dimers). The particles and their dimers levitated in a plasma sheath at slightly different heights and formed two distinct sublayers. The system did not crystallize and may be characterized as a disordered solid. The dimers were identified based on their characteristic appearance in defocused images, i.e., rotating interference fringe patterns. The in-plane and interplane particle separations exhibit nonmonotonic dependence on the discharge pressure.

Potential formation in a one-dimensional bounded plasma system containing a two-electron temperature plasma: Kinetic model and PIC simulation

International Nuclear Information System (INIS)

Gyergyek, T.; Jurcic-Zlobec, B.; Cercek, M.

2008-01-01

Potential formation in a bounded plasma system that contains electrons with a two-temperature velocity distribution and is terminated by a floating, electron emitting electrode (collector) is studied by a one-dimensional kinetic model. A method on how to determine the boundary conditions at the collector for the numerical solution of the Poisson equation is presented. The difference between the regular and the irregular numerical solutions of the Poisson equation is explained. The regular numerical solution of the Poisson equation fulfills the boundary conditions at the source and can be computed for any distance from the collector. The irregular solution does not fulfill the source boundary conditions and the computation breaks down at some distance from the collector. An excellent agreement of the values of the potential at the inflection point found from the numerical solution of the Poisson equation with the values predicted by the analytical model is obtained. Potential, electric field, and particle density profiles found by the numerical solution of the Poisson equation are compared to the profiles obtained with the particle in cell computer simulation. A very good quantitative agreement of the potential and electric field profiles is obtained. For certain values of the parameters the analytical model predicts three possible values of the potential at the inflection point. In such cases always only one of the corresponding numerical solutions of the Poisson equation is regular, while the other two are irregular. The regular numerical solution of the Poisson equation always corresponds to the solution of the model that predicts the largest ion flux to the collector

Nonlinear PIC simulation in a Penning trap

International Nuclear Information System (INIS)

Lapenta, G.; Delzanno, G.L.; Finn, J. M.

2002-01-01

We study the nonlinear dynamics of a Penning trap plasma, including the effect of the finite length and end curvature of the plasma column. A new cylindrical PIC code, called KANDINSKY, has been implemented by using a new interpolation scheme. The principal idea is to calculate the volume of each cell from a particle volume, in the same manner as it is done for the cell charge. With this new method, the density is conserved along streamlines and artificial sources of compressibility are avoided. The code has been validated with a reference Eulerian fluid code. We compare the dynamics of three different models: a model with compression effects, the standard Euler model and a geophysical fluid dynamics model. The results of our investigation prove that Penning traps can really be used to simulate geophysical fluids

Recent advances in the modeling of plasmas with the Particle-In-Cell methods

Science.gov (United States)

Vay, Jean-Luc; Lehe, Remi; Vincenti, Henri; Godfrey, Brendan; Lee, Patrick; Haber, Irv

2015-11-01

The Particle-In-Cell (PIC) approach is the method of choice for self-consistent simulations of plasmas from first principles. The fundamentals of the PIC method were established decades ago but improvements or variations are continuously being proposed. We report on several recent advances in PIC related algorithms, including: (a) detailed analysis of the numerical Cherenkov instability and its remediation, (b) analytic pseudo-spectral electromagnetic solvers in Cartesian and cylindrical (with azimuthal modes decomposition) geometries, (c) arbitrary-order finite-difference and generalized pseudo-spectral Maxwell solvers, (d) novel analysis of Maxwell's solvers' stencil variation and truncation, in application to domain decomposition strategies and implementation of Perfectly Matched Layers in high-order and pseudo-spectral solvers. Work supported by US-DOE Contracts DE-AC02-05CH11231 and the US-DOE SciDAC program ComPASS. Used resources of NERSC, supported by US-DOE Contract DE-AC02-05CH11231.

Terahertz radiation driven by two-color laser pulses at near-relativistic intensities: Competition between photoionization and wakefield effects

Science.gov (United States)

González de Alaiza Martínez, P.; Davoine, X.; Debayle, A.; Gremillet, L.; Bergé, L.

2016-01-01

We numerically investigate terahertz (THz) pulse generation by linearly-polarized, two-color femtosecond laser pulses in highly-ionized argon. Major processes consist of tunneling photoionization and ponderomotive forces associated with transverse and longitudinal field excitations. By means of two-dimensional particle-in-cell (PIC) simulations, we reveal the importance of photocurrent mechanisms besides transverse and longitudinal plasma waves for laser intensities >1015 W/cm2. We demonstrate the following. (i) With two-color pulses, photoionization prevails in the generation of GV/m THz fields up to 1017 W/cm2 laser intensities and suddenly loses efficiency near the relativistic threshold, as the outermost electron shell of ionized Ar atoms has been fully depleted. (ii) PIC results can be explained by a one-dimensional Maxwell-fluid model and its semi-analytical solutions, offering the first unified description of the main THz sources created in plasmas. (iii) The THz power emitted outside the plasma channel mostly originates from the transverse currents. PMID:27255689

Comparison of multi-fluid moment models with particle-in-cell simulations of collisionless magnetic reconnection

International Nuclear Information System (INIS)

Wang, Liang; Germaschewski, K.; Hakim, Ammar H.; Bhattacharjee, A.

2015-01-01

We introduce an extensible multi-fluid moment model in the context of collisionless magnetic reconnection. This model evolves full Maxwell equations and simultaneously moments of the Vlasov-Maxwell equation for each species in the plasma. Effects like electron inertia and pressure gradient are self-consistently embedded in the resulting multi-fluid moment equations, without the need to explicitly solving a generalized Ohm's law. Two limits of the multi-fluid moment model are discussed, namely, the five-moment limit that evolves a scalar pressures for each species and the ten-moment limit that evolves the full anisotropic, non-gyrotropic pressure tensor for each species. We first demonstrate analytically and numerically that the five-moment model reduces to the widely used Hall magnetohydrodynamics (Hall MHD) model under the assumptions of vanishing electron inertia, infinite speed of light, and quasi-neutrality. Then, we compare ten-moment and fully kinetic particle-in-cell (PIC) simulations of a large scale Harris sheet reconnection problem, where the ten-moment equations are closed with a local linear collisionless approximation for the heat flux. The ten-moment simulation gives reasonable agreement with the PIC results regarding the structures and magnitudes of the electron flows, the polarities and magnitudes of elements of the electron pressure tensor, and the decomposition of the generalized Ohm's law. Possible ways to improve the simple local closure towards a nonlocal fully three-dimensional closure are also discussed

Many-particle correlations in quasi-two-dimensional electron-hole systems

International Nuclear Information System (INIS)

Nikolaev, Valentin

2002-01-01

This thesis reports a theoretical investigation of many-particle correlation effects in semiconductor heterostructures containing quantum wells. Particular attention is paid towards quasi-particle pair correlations. Using the Green's function technique and the ladder approximation as a basis, the generalized mass action law, which describes the redistribution of particles between correlated and uncorrelated states in quasi-two-dimensional systems for different temperatures and total densities, is derived. The expression is valid beyond the low-density limit, which allows us to investigate the transition of the system from a dilute exciton gas to a dense electron-hole plasma. A generalized Levinson theorem, which takes k-space filling into account, is formulated. Screening in quasi-two-dimensional systems is analyzed rigorously. Firstly, the qualitatively new mechanism of static local screening by indirect excitons is studied using the simple Thomas-Fermi approximation. Then, a detailed many-body description suitable for a proper account of dynamic screening by a quasi-2D electron-hole plasma, and consistent with the previously derived mass action law, is provided. The generalized Lindhard approximation and excitonic plasmon-pole approximations are also derived. The theory is applied to single and double quantum wells. A self-consistent procedure is developed for numerical investigation of the ionization degree of an electron-hole plasma at different values of temperature/exciton Rydberg ratios. This procedure accounts for screening, k-space filling (exciton bleaching), and the formation of excitons. An abrupt jump in the value of the ionization degree that happens with an increase of the carrier density or temperature (Mott transition) is found in a certain density-temperature region. It has been found that the critical density of the Mott transition for indirect excitons may be much smaller than that for direct excitons. A suggestion has been made that some of the

SPECT3D - A multi-dimensional collisional-radiative code for generating diagnostic signatures based on hydrodynamics and PIC simulation output

Science.gov (United States)

MacFarlane, J. J.; Golovkin, I. E.; Wang, P.; Woodruff, P. R.; Pereyra, N. A.

2007-05-01

SPECT3D is a multi-dimensional collisional-radiative code used to post-process the output from radiation-hydrodynamics (RH) and particle-in-cell (PIC) codes to generate diagnostic signatures (e.g. images, spectra) that can be compared directly with experimental measurements. This ability to post-process simulation code output plays a pivotal role in assessing the reliability of RH and PIC simulation codes and their physics models. SPECT3D has the capability to operate on plasmas in 1D, 2D, and 3D geometries. It computes a variety of diagnostic signatures that can be compared with experimental measurements, including: time-resolved and time-integrated spectra, space-resolved spectra and streaked spectra; filtered and monochromatic images; and X-ray diode signals. Simulated images and spectra can include the effects of backlighters, as well as the effects of instrumental broadening and time-gating. SPECT3D also includes a drilldown capability that shows where frequency-dependent radiation is emitted and absorbed as it propagates through the plasma towards the detector, thereby providing insights on where the radiation seen by a detector originates within the plasma. SPECT3D has the capability to model a variety of complex atomic and radiative processes that affect the radiation seen by imaging and spectral detectors in high energy density physics (HEDP) experiments. LTE (local thermodynamic equilibrium) or non-LTE atomic level populations can be computed for plasmas. Photoabsorption rates can be computed using either escape probability models or, for selected 1D and 2D geometries, multi-angle radiative transfer models. The effects of non-thermal (i.e. non-Maxwellian) electron distributions can also be included. To study the influence of energetic particles on spectra and images recorded in intense short-pulse laser experiments, the effects of both relativistic electrons and energetic proton beams can be simulated. SPECT3D is a user-friendly software package that runs

The limiting current in a one-dimensional situation: Transition from a space charge limited to magnetically limited flow

International Nuclear Information System (INIS)

Kumar, Raghwendra; Biswas, Debabrata

2008-01-01

For a nonrelativistic electron beam propagating in a cylindrical drift tube, it is shown that the limiting current density does not saturate to the electrostatic one-dimensional (1D) estimate with increasing beam radius. Fully electromagnetic particle-in-cell (PIC) simulation studies show that beyond a critical aspect ratio, the limiting current density is lower than the 1D electrostatic prediction. The lowering in the limiting current density is found to be due to the transition from the space charge limited to magnetically limited flow. An adaptation of Alfven's single particle trajectory method is used to estimate the magnetically limited current as well as the critical radius beyond which the flow is magnetically limited in a drift tube. The predictions are found to be in close agreement with PIC simulations

Measurement of correlations between low-frequency vibrational modes and particle rearrangements in quasi-two-dimensional colloidal glasses

NARCIS (Netherlands)

Chen, K.; Manning, M.L.; Yunker, P.J.; Ellenbroek, W.G.; Zhang, Zexin; Liu, Andrea J.; Yodh, A.G.

2011-01-01

We investigate correlations between low-frequency vibrational modes and rearrangements in two-dimensional colloidal glasses composed of thermosensitive microgel particles, which readily permit variation of the sample packing fraction. At each packing fraction, the particle displacement covariance

Tracer particles in two-dimensional elastic networks diffuse logarithmically slow

International Nuclear Information System (INIS)

Lizana, Ludvig; Ambjörnsson, Tobias; Lomholt, Michael A

2017-01-01

Several experiments on tagged molecules or particles in living systems suggest that they move anomalously slow—their mean squared displacement (MSD) increase slower than linearly with time. Leading models aimed at understanding these experiments predict that the MSD grows as a power law with a growth exponent that is smaller than unity. However, in some experiments the growth is so slow (fitted exponent  ∼0.1–0.2) that they hint towards other mechanisms at play. In this paper, we theoretically demonstrate how in-plane collective modes excited by thermal fluctuations in a two dimensional membrane lead to logarithmic time dependence for the the tracer particle’s MSD. (paper)

Magnetic manipulation of particles and cells in ferrofluid flow through straight microchannels using two magnets

Science.gov (United States)

Zeng, Jian

Microfluidic devices have been increasingly used in the past two decades for particle and cell manipulations in many chemical and biomedical applications. A variety of force fields have been demonstrated to control particle and cell transport in these devices including electric, magnetic, acoustic, and optical forces etc. Among these particle handling techniques, the magnetic approach provides clear advantages over others such as low cost, noninvasive, and free of fluid heating issues. However, the current knowledge of magnetic control of particle transport is still very limited, especially lacking is the handling of diamagnetic particle. This thesis is focused on the magnetic manipulation of diamagnetic particles and cells in ferrofluid flow through the use of a pair of permanent magnets. By varying the configuration of the two magnets, diverse operations of particles and cells is implemented in a straight microchannel that can potentially be integrated into lab-on-a-chip devices for various applications. First, an approach for embedding two, symmetrically positioned, repulsive permanent magnets about a straight rectangular microchannel in a PDMS-based microfluidic device is developed for particle focusing. Focusing particles and cells into a tight stream is often required in order for continuous detection, counting, and sorting. The closest distance between the magnets is limited only by the size of the magnets involved in the fabrication process. The device is used to implement and investigate the three-dimensional magnetic focusing of polystyrene particles in ferrofluid microflow with both top-view and side-view visualizations. The effects of flow speed and particle size on the particle focusing effectiveness are studied. This device is also applied to magnetically focus yeast cells in ferrofluid, which proves to be biocompatible as verified by cell viability test. In addition, an analytical model is developed and found to be able to predict the experimentally

Loading relativistic Maxwell distributions in particle simulations

Energy Technology Data Exchange (ETDEWEB)

Zenitani, Seiji, E-mail: seiji.zenitani@nao.ac.jp [National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588 (Japan)

2015-04-15

Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ≈50% for generic cases and 100% for symmetric distributions. They can be combined with arbitrary base algorithms.

Loading relativistic Maxwell distributions in particle simulations

International Nuclear Information System (INIS)

Zenitani, Seiji

2015-01-01

Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ≈50% for generic cases and 100% for symmetric distributions. They can be combined with arbitrary base algorithms

Novel methods in the Particle-In-Cell accelerator Code-Framework Warp

Energy Technology Data Exchange (ETDEWEB)

Vay, J-L [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Grote, D. P. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Cohen, R. H. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Friedman, A. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

2012-12-26

The Particle-In-Cell (PIC) Code-Framework Warp is being developed by the Heavy Ion Fusion Science Virtual National Laboratory (HIFS-VNL) to guide the development of accelerators that can deliver beams suitable for high-energy density experiments and implosion of inertial fusion capsules. It is also applied in various areas outside the Heavy Ion Fusion program to the study and design of existing and next-generation high-energy accelerators, including the study of electron cloud effects and laser wakefield acceleration for example. This study presents an overview of Warp's capabilities, summarizing recent original numerical methods that were developed by the HIFS-VNL (including PIC with adaptive mesh refinement, a large-timestep 'drift-Lorentz' mover for arbitrarily magnetized species, a relativistic Lorentz invariant leapfrog particle pusher, simulations in Lorentz-boosted frames, an electromagnetic solver with tunable numerical dispersion and efficient stride-based digital filtering), with special emphasis on the description of the mesh refinement capability. In addition, selected examples of the applications of the methods to the abovementioned fields are given.

Particle-in-cell plasma simulations of the modified two-stream instability

Directory of Open Access Journals (Sweden)

K. Schlegel

1994-08-01

Full Text Available We model the modified two-stream plasma instability occurring in the ionospheric E-region using a 2.5-dimensional particle-in-cell code. Compared to previous similar work we concentrate on simulated quantities that can easily be measured in the real ionosphere by coherent radars or rockets, such as the Doppler velocity, the backscattered power, backscattered spectra, aspect angle behaviour and electron temperature enhancement. Despite using a relatively small simulation model, we obtain remarkably good agreement between actual observed and simulated plasma parameters. The advantage of such a small system is that we were able to perform (other than in previous related work many simulation runs with different sets of input parameters, thus studying the unstable plasma under various conditions.

Fabricating large two-dimensional single colloidal crystals by doping with active particles

NARCIS (Netherlands)

van der Meer, B; Filion, L; Dijkstra, M

2016-01-01

Using simulations we explore the behaviour of two-dimensional colloidal (poly)crystals doped with active particles. We show that these active dopants can provide an elegant new route to removing grain boundaries in polycrystals. Specifically, we show that active dopants both generate and are

A 3D gyrokinetic particle-in-cell simulation of fusion plasma microturbulence on parallel computers

Science.gov (United States)

Williams, T. J.

1992-12-01

One of the grand challenge problems now supported by HPCC is the Numerical Tokamak Project. A goal of this project is the study of low-frequency micro-instabilities in tokamak plasmas, which are believed to cause energy loss via turbulent thermal transport across the magnetic field lines. An important tool in this study is gyrokinetic particle-in-cell (PIC) simulation. Gyrokinetic, as opposed to fully-kinetic, methods are particularly well suited to the task because they are optimized to study the frequency and wavelength domain of the microinstabilities. Furthermore, many researchers now employ low-noise delta(f) methods to greatly reduce statistical noise by modelling only the perturbation of the gyrokinetic distribution function from a fixed background, not the entire distribution function. In spite of the increased efficiency of these improved algorithms over conventional PIC algorithms, gyrokinetic PIC simulations of tokamak micro-turbulence are still highly demanding of computer power--even fully-vectorized codes on vector supercomputers. For this reason, we have worked for several years to redevelop these codes on massively parallel computers. We have developed 3D gyrokinetic PIC simulation codes for SIMD and MIMD parallel processors, using control-parallel, data-parallel, and domain-decomposition message-passing (DDMP) programming paradigms. This poster summarizes our earlier work on codes for the Connection Machine and BBN TC2000 and our development of a generic DDMP code for distributed-memory parallel machines. We discuss the memory-access issues which are of key importance in writing parallel PIC codes, with special emphasis on issues peculiar to gyrokinetic PIC. We outline the domain decompositions in our new DDMP code and discuss the interplay of different domain decompositions suited for the particle-pushing and field-solution components of the PIC algorithm.

On the size distribution of one-, two- and three-dimensional Voronoi cells

International Nuclear Information System (INIS)

Marthinsen, K.

1994-03-01

The present report gives a presentation of the different cell size distribution obtained by computer simulations of random Voronoi cell structures in one-, two- and three-dimensional space. The random Voronoi cells are constructed from cell centroids randomly distributed along a string, in the plane and in three-dimensional space, respectively. The size distributions are based on 2-3 · 10 4 cells. For the spacial polyhedra both the distribution of volumes, areas and radii are presented, and the two latter quantities are compared to the distributions of areas and radii from a planar section through the three-dimensional structure as well as to the corresponding distributions obtained from a pure two-dimensional cell structure. 11 refs., 11 figs

Development of 2D particle-in-cell code to simulate high current, low ...

Indian Academy of Sciences (India)

Abstract. A code for 2D space-charge dominated beam dynamics study in beam trans- port lines is developed. The code is used for particle-in-cell (PIC) simulation of z-uniform beam in a channel containing solenoids and drift space. It can also simulate a transport line where quadrupoles are used for focusing the beam.

Optimal conclusive teleportation of a d-dimensional two-particle unknown quantum state

Institute of Scientific and Technical Information of China (English)

Yang Yu-Guang; Wen Qiao-Yan; Zhu Fu-Chen

2006-01-01

A conclusive teleportation protocol of a d-dimensional two-particle unknown quantum state using three ddimensional particles in an arbitrary pure state is proposed. A sender teleports the unknown state conclusively to a receiver by using the positive operator valued measure(POVM) and introducing an ancillary qudit to perform the generalized Bell basis measurement. We calculate the optimal teleportation fidelity. We also discuss and analyse the reason why the information on the teleported state is lost in the course of the protocol.

Three-dimensional relativistic pair plasma reconnection with radiative feedback in the Crab Nebula

Energy Technology Data Exchange (ETDEWEB)

Cerutti, B. [Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 (United States); Werner, G. R.; Uzdensky, D. A. [Center for Integrated Plasma Studies, Physics Department, University of Colorado, UCB 390, Boulder, CO 80309-0390 (United States); Begelman, M. C., E-mail: bcerutti@astro.princeton.edu, E-mail: greg.werner@colorado.edu, E-mail: uzdensky@colorado.edu, E-mail: mitch@jila.colorado.edu [JILA, University of Colorado and National Institute of Standards and Technology, UCB 440, Boulder, CO 80309-0440 (United States)

2014-02-20

The discovery of rapid synchrotron gamma-ray flares above 100 MeV from the Crab Nebula has attracted new interest in alternative particle acceleration mechanisms in pulsar wind nebulae. Diffuse shock-acceleration fails to explain the flares because particle acceleration and emission occur during a single or even sub-Larmor timescale. In this regime, the synchrotron energy losses induce a drag force on the particle motion that balances the electric acceleration and prevents the emission of synchrotron radiation above 160 MeV. Previous analytical studies and two-dimensional (2D) particle-in-cell (PIC) simulations indicate that relativistic reconnection is a viable mechanism to circumvent the above difficulties. The reconnection electric field localized at X-points linearly accelerates particles with little radiative energy losses. In this paper, we check whether this mechanism survives in three dimension (3D), using a set of large PIC simulations with radiation reaction force and with a guide field. In agreement with earlier works, we find that the relativistic drift kink instability deforms and then disrupts the layer, resulting in significant plasma heating but few non-thermal particles. A moderate guide field stabilizes the layer and enables particle acceleration. We report that 3D magnetic reconnection can accelerate particles above the standard radiation reaction limit, although the effect is less pronounced than in 2D with no guide field. We confirm that the highest-energy particles form compact bunches within magnetic flux ropes, and a beam tightly confined within the reconnection layer, which could result in the observed Crab flares when, by chance, the beam crosses our line of sight.

FPGA Implementation of one-dimensional and two-dimensional cellular automata

International Nuclear Information System (INIS)

D'Antone, I.

1999-01-01

This report describes the hardware implementation of one-dimensional and two-dimensional cellular automata (CAs). After a general introduction to the cellular automata, we consider a one-dimensional CA used to implement pseudo-random techniques in built-in self test for VLSI. Due to the increase in digital ASIC complexity, testing is becoming one of the major costs in the VLSI production. The high electronics complexity, used in particle physics experiments, demands higher reliability than in the past time. General criterions are given to evaluate the feasibility of the circuit used for testing and some quantitative parameters are underlined to optimize the architecture of the cellular automaton. Furthermore, we propose a two-dimensional CA that performs a peak finding algorithm in a matrix of cells mapping a sub-region of a calorimeter. As in a two-dimensional filtering process, the peaks of the energy clusters are found in one evolution step. This CA belongs to Wolfram class II cellular automata. Some quantitative parameters are given to optimize the architecture of the cellular automaton implemented in a commercial field programmable gate array (FPGA)

PIConGPU - How to build one of the fastest GPU particle-in-cell codes in the world

Energy Technology Data Exchange (ETDEWEB)

Burau, Heiko; Debus, Alexander; Helm, Anton; Huebl, Axel; Kluge, Thomas; Widera, Rene; Bussmann, Michael; Schramm, Ulrich; Cowan, Thomas [HZDR, Dresden (Germany); Juckeland, Guido; Nagel, Wolfgang [TU Dresden (Germany); ZIH, Dresden (Germany); Schmitt, Felix [NVIDIA (United States)

2013-07-01

We present the algorithmic building blocks of PIConGPU, one of the fastest implementations of the particle-in-cell algortihm on GPU clusters. PIConGPU is a highly-scalable, 3D3V electromagnetic PIC code that is used in laser plasma and astrophysical plasma simulations.

Memory-efficient optimization of Gyrokinetic particle-to-grid interpolation for multicore processors

Energy Technology Data Exchange (ETDEWEB)

Madduri, Kamesh [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Williams, Samuel [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Ethier, Stephane [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Oliker, Leonid [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Shalf, John [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Strohmaier, Erich [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Yelicky, Katherine [Univ. of California, Berkeley, CA (United States)

2009-01-01

We present multicore parallelization strategies for the particle-to-grid interpolation step in the Gyrokinetic Toroidal Code (GTC), a 3D particle-in-cell (PIC) application to study turbulent transport in magnetic-confinement fusion devices. Particle-grid interpolation is a known performance bottleneck in several PIC applications. In GTC, this step involves particles depositing charges to a 3D toroidal mesh, and multiple particles may contribute to the charge at a grid point. We design new parallel algorithms for the GTC charge deposition kernel, and analyze their performance on three leading multicore platforms. We implement thirteen different variants for this kernel and identify the best-performing ones given typical PIC parameters such as the grid size, number of particles per cell, and the GTC-specific particle Larmor radius variation. We find that our best strategies can be 2x faster than the reference optimized MPI implementation, and our analysis provides insight into desirable architectural features for high-performance PIC simulation codes.

Wavelet-Based Poisson Solver for Use in Particle-in-Cell Simulations

CERN Document Server

Terzic, Balsa; Mihalcea, Daniel; Pogorelov, Ilya V

2005-01-01

We report on a successful implementation of a wavelet-based Poisson solver for use in 3D particle-in-cell simulations. One new aspect of our algorithm is its ability to treat the general (inhomogeneous) Dirichlet boundary conditions. The solver harnesses advantages afforded by the wavelet formulation, such as sparsity of operators and data sets, existence of effective preconditioners, and the ability simultaneously to remove numerical noise and further compress relevant data sets. Having tested our method as a stand-alone solver on two model problems, we merged it into IMPACT-T to obtain a fully functional serial PIC code. We present and discuss preliminary results of application of the new code to the modelling of the Fermilab/NICADD and AES/JLab photoinjectors.

PIC simulation of a thermal anisotropy-driven Weibel instability in a circular rarefaction wave

International Nuclear Information System (INIS)

Dieckmann, M E; Sarri, G; Kourakis, I; Borghesi, M; Murphy, G C; O'C Drury, L; Bret, A; Romagnani, L; Ynnerman, A

2012-01-01

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here whether the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here that this electric field modification triggers a second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally. (paper)

PIC simulation of a thermal anisotropy-driven Weibel instability in a circular rarefaction wave

Science.gov (United States)

Dieckmann, M. E.; Sarri, G.; Murphy, G. C.; Bret, A.; Romagnani, L.; Kourakis, I.; Borghesi, M.; Ynnerman, A.; O'C Drury, L.

2012-02-01

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here whether the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here that this electric field modification triggers a second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally.

3-D electromagnetic plasma particle simulations on the Intel Delta parallel computer

International Nuclear Information System (INIS)

Wang, J.; Liewer, P.C.

1994-01-01

A three-dimensional electromagnetic PIC code has been developed on the 512 node Intel Touchstone Delta MIMD parallel computer. This code is based on the General Concurrent PIC algorithm which uses a domain decomposition to divide the computation among the processors. The 3D simulation domain can be partitioned into 1-, 2-, or 3-dimensional sub-domains. Particles must be exchanged between processors as they move among the subdomains. The Intel Delta allows one to use this code for very-large-scale simulations (i.e. over 10 8 particles and 10 6 grid cells). The parallel efficiency of this code is measured, and the overall code performance on the Delta is compared with that on Cray supercomputers. It is shown that their code runs with a high parallel efficiency of ≥ 95% for large size problems. The particle push time achieved is 115 nsecs/particle/time step for 162 million particles on 512 nodes. Comparing with the performance on a single processor Cray C90, this represents a factor of 58 speedup. The code uses a finite-difference leap frog method for field solve which is significantly more efficient than fast fourier transforms on parallel computers. The performance of this code on the 128 node Cray T3D will also be discussed

An electrostatic particle-in-cell model for a lower hybrid grill

International Nuclear Information System (INIS)

Rantamaeki, K.

1998-01-01

In recent lower hybrid (LH) current drive experiments, generation of hot spots and impurities in the grill region have been observed on Tore Supra and Tokamak de Varennes (TdeV). A possible explanation is the parasitic absorption of the LH power in front of the grill. In parasitic absorption, the short-wavelength part of the lower hybrid spectrum can resonantly interact with the cold edge electrons. In this work, the absorption of the LH waves and the generation of fast electrons near the waveguide mouth is investigated with a new tool in this context: particle-in-cell (PIC) simulations. The advantage of this new method is that the electric field is calculated self-consistently. The PIC simulations also provide the key parameters for the hot spot problem: the absorbed power, the radial deposition profiles and the absorption length. A grill model has been added to the 2d3v PIC code XPDP2. Two sets of simulations were made. The first simulations used a phenomenological grill model. Strong absorption in the edge plasma was obtained. About 5% of the coupled power was absorbed within 1.7 mm in the case with fairly large amount of power in the modes with large parallel refractive index. Consequently, a rapid generation of fast electrons took place in the same region. In order to model experiments with realistic wave spectra, the PIC code was coupled to the slow wave antenna coupling code SWAN. The absorption within 1.7 mm in front of the grill was found to be between 2 and 5%. In the short time of a few wave periods, part of the initially thermal electrons (T e = 100 eV) were accelerated to velocities corresponding to a few keV. (orig.)

An electrostatic particle-in-cell model for a lower hybrid grill

Energy Technology Data Exchange (ETDEWEB)

Rantamaeki, K

1998-07-01

In recent lower hybrid (LH) current drive experiments, generation of hot spots and impurities in the grill region have been observed on Tore Supra and Tokamak de Varennes (TdeV). A possible explanation is the parasitic absorption of the LH power in front of the grill. In parasitic absorption, the short-wavelength part of the lower hybrid spectrum can resonantly interact with the cold edge electrons. In this work, the absorption of the LH waves and the generation of fast electrons near the waveguide mouth is investigated with a new tool in this context: particle-in-cell (PIC) simulations. The advantage of this new method is that the electric field is calculated self-consistently. The PIC simulations also provide the key parameters for the hot spot problem: the absorbed power, the radial deposition profiles and the absorption length. A grill model has been added to the 2d3v PIC code XPDP2. Two sets of simulations were made. The first simulations used a phenomenological grill model. Strong absorption in the edge plasma was obtained. About 5% of the coupled power was absorbed within 1.7 mm in the case with fairly large amount of power in the modes with large parallel refractive index. Consequently, a rapid generation of fast electrons took place in the same region. In order to model experiments with realistic wave spectra, the PIC code was coupled to the slow wave antenna coupling code SWAN. The absorption within 1.7 mm in front of the grill was found to be between 2 and 5%. In the short time of a few wave periods, part of the initially thermal electrons (T{sub e} = 100 eV) were accelerated to velocities corresponding to a few keV. (orig.)

Parallel treatment of simulation particles in particle-in-cell codes on SUPRENUM

International Nuclear Information System (INIS)

Seldner, D.

1990-02-01

This report contains the program documentation and description of the program package 2D-PLAS, which has been developed at the Nuclear Research Center Karlsruhe in the Institute for Data Processing in Technology (IDT) under the auspices of the BMFT. 2D-PLAS is a parallel program version of the treatment of the simulation particles of the two-dimensional stationary particle-in-cell code BFCPIC which has been developed at the Nuclear Research Center Karlsruhe. This parallel version has been designed for the parallel computer SUPRENUM. (orig.) [de

Flame propagation in two-dimensional solids: Particle-resolved studies with complex plasmas

Science.gov (United States)

Yurchenko, S. O.; Yakovlev, E. V.; Couëdel, L.; Kryuchkov, N. P.; Lipaev, A. M.; Naumkin, V. N.; Kislov, A. Yu.; Ovcharov, P. V.; Zaytsev, K. I.; Vorob'ev, E. V.; Morfill, G. E.; Ivlev, A. V.

2017-10-01

Using two-dimensional (2D) complex plasmas as an experimental model system, particle-resolved studies of flame propagation in classical 2D solids are carried out. Combining experiments, theory, and molecular dynamics simulations, we demonstrate that the mode-coupling instability operating in 2D complex plasmas reveals all essential features of combustion, such as an activated heat release, two-zone structure of the self-similar temperature profile ("flame front"), as well as thermal expansion of the medium and temperature saturation behind the front. The presented results are of relevance for various fields ranging from combustion and thermochemistry, to chemical physics and synthesis of materials.

Analysis of the beam halo in negative ion sources by using 3D3V PIC code

Energy Technology Data Exchange (ETDEWEB)

Miyamoto, K., E-mail: kmiyamot@naruto-u.ac.jp [Naruto University of Education, 748 Nakashima, Takashima, Naruto-cho, Naruto-shi, Tokushima 772-8502 (Japan); Nishioka, S.; Goto, I.; Hatayama, A. [Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522 (Japan); Hanada, M.; Kojima, A.; Hiratsuka, J. [Japan Atomic Energy Agency, 801-1 Mukouyama, Naka 319-0913 (Japan)

2016-02-15

The physical mechanism of the formation of the negative ion beam halo and the heat loads of the multi-stage acceleration grids are investigated with the 3D PIC (particle in cell) simulation. The following physical mechanism of the beam halo formation is verified: The beam core and the halo consist of the negative ions extracted from the center and the periphery of the meniscus, respectively. This difference of negative ion extraction location results in a geometrical aberration. Furthermore, it is shown that the heat loads on the first acceleration grid and the second acceleration grid are quantitatively improved compared with those for the 2D PIC simulation result.

Discrete particle simulation of bubble and slug formation in a two-dimensional gas-fluidised bed: A hard-sphere approach.

NARCIS (Netherlands)

Hoomans, B.P.B.; Kuipers, J.A.M.; Briels, Willem J.; van Swaaij, Willibrordus Petrus Maria

1996-01-01

A discrete particle model of a gas-fluidised bed has been developed and in this the two-dimensional motion of the individual, spherical particles was directly calculated from the forces acting on them, accounting for the interaction between the particles and the interstitial gas phase. Our collision

A Generalized Weight-Based Particle-In-Cell Simulation Scheme

International Nuclear Information System (INIS)

Lee, W.W.; Jenkins, T.G.; Ethier, S.

2010-01-01

A generalized weight-based particle simulation scheme suitable for simulating magnetized plasmas, where the zeroth-order inhomogeneity is important, is presented. The scheme is an extension of the perturbative simulation schemes developed earlier for particle-in-cell (PIC) simulations. The new scheme is designed to simulate both the perturbed distribution ((delta)f) and the full distribution (full-F) within the same code. The development is based on the concept of multiscale expansion, which separates the scale lengths of the background inhomogeneity from those associated with the perturbed distributions. The potential advantage for such an arrangement is to minimize the particle noise by using (delta)f in the linear stage stage of the simulation, while retaining the flexibility of a full-F capability in the fully nonlinear stage of the development when signals associated with plasma turbulence are at a much higher level than those from the intrinsic particle noise.

2D PIC simulations for an EN discharge with magnetized electrons and unmagnetized ions

Science.gov (United States)

Lieberman, Michael A.; Kawamura, Emi; Lichtenberg, Allan J.

2009-10-01

We conducted 2D particle-in-cell (PIC) simulations for an electronegative (EN) discharge with magnetized electrons and unmagnetized ions, and compared the results to a previously developed 1D (radial) analytical model of an EN plasma with strongly magnetized electrons and weakly magnetized ions [1]. In both cases, there is a static uniform applied magnetic field in the axial direction. The 1D radial model mimics the wall losses of the particles in the axial direction by introducing a bulk loss frequency term νL. A special (desired) solution was found in which only positive and negative ions but no electrons escaped radially. The 2D PIC results show good agreement with the 1D model over a range of parameters and indicate that the analytical form of νL employed in [1] is reasonably accurate. However, for the PIC simulations, there is always a finite flux of electrons to the radial wall which is about 10 to 30% of the negative ion flux.[4pt] [1] G. Leray, P. Chabert, A.J. Lichtenberg and M.A. Lieberman, J. Phys. D, accepted for publication 2009.

The one-particle scenario for the metal-insulator transition in two-dimensional systems at T = 0

CERN Document Server

Tarasov, Y V

2003-01-01

The conductance of bounded disordered electron systems is calculated by reducing the original dynamic problem of arbitrary dimensionality to a set of strictly one-dimensional problems for one-particle mode propagators. The metallic ground state of a two-dimensional conductor, which is considered as a limiting case of three-dimensional quantum waveguide, is shown to result from its multi-modeness. As the waveguide thickness is reduced, e.g., by applying a 'pressing' potential, the electron system undergoes a set of continuous phase transitions related to discrete variations of the number of extended modes. The closing of the last current carrying mode is regarded as a phase transition of the electron system from metallic to dielectric state. The obtained results agree qualitatively with the observed 'anomalies' of resistivity of different two-dimensional electron and hole systems.

A Particle-In-Cell approach to particle flux shaping with a surface mask

Directory of Open Access Journals (Sweden)

G. Kawamura

2017-08-01

Full Text Available The Particle-In-Cell simulation code PICS has been developed to study plasma in front of a surface with two types of masks, step-type and roof-type. Parameter scans with regard to magnetic field angle, electron density, and mask height were carried out to understand their influence on ion particle flux distribution on a surface. A roof-type mask with a small mask height yields short decay length in the flux distribution which is consistent with that estimated experimentally. A roof-type mask with a large height yields very long decay length and the flux value does not depend on a mask height or an electron density, but rather on a mask length and a biasing voltage of the surface. Mask height also changes the flux distribution apart from the mask because of the shading effect of the mask. Electron density changes the distribution near the mask edge according to the Debye length. Dependence of distribution on parameters are complicated especially for a roof-type mask, and simulation study with various parameters are useful to understand the physical reasons of dependence and also is useful as a tool for experiment studies.

Particle Control in Phase Space by Global K-Means Clustering

DEFF Research Database (Denmark)

Frederiksen, Jacob Trier; Lapenta, G.; Pessah, M. E.

2015-01-01

We devise and explore an iterative optimization procedure for controlling particle populations in particle-in-cell (PIC) codes via merging and splitting of computational macro-particles. Our approach, is to compute an optimal representation of the global particle phase space structure while decre...

Two-particle correlations in the one-dimensional Hubbard model: a ground-state analytical solution

CERN Document Server

Vallejo, E; Espinosa, J E

2003-01-01

A solution to the extended Hubbard Hamiltonian for the case of two-particles in an infinite one-dimensional lattice is presented, using a real-space mapping method and the Green function technique. This Hamiltonian considers the on-site (U) and the nearest-neighbor (V) interactions. The method is based on mapping the correlated many-body problem onto an equivalent site-impurity tight-binding one in a higher dimensional space. In this new space we obtained the analytical solution for the ground state binding energy. Results are in agreement with the numerical solution obtained previously [1], and with those obtained in the reciprocal space [2]. (Author)

Horizons in 2+1-dimensional collapse of particles

Indian Academy of Sciences (India)

A simple, geometrical construction is given for three-dimensional spacetimes with negative cosmological constant that contain two particles colliding head-on. Depending on parameters like particle masses and distance, the combined geometry will be that of a particle, or of a black hole. In the black hole case the horizon is ...

Recent progress in 3D EM/EM-PIC simulation with ARGUS and parallel ARGUS

International Nuclear Information System (INIS)

Mankofsky, A.; Petillo, J.; Krueger, W.; Mondelli, A.; McNamara, B.; Philp, R.

1994-01-01

ARGUS is an integrated, 3-D, volumetric simulation model for systems involving electric and magnetic fields and charged particles, including materials embedded in the simulation region. The code offers the capability to carry out time domain and frequency domain electromagnetic simulations of complex physical systems. ARGUS offers a boolean solid model structure input capability that can include essentially arbitrary structures on the computational domain, and a modular architecture that allows multiple physics packages to access the same data structure and to share common code utilities. Physics modules are in place to compute electrostatic and electromagnetic fields, the normal modes of RF structures, and self-consistent particle-in-cell (PIC) simulation in either a time dependent mode or a steady state mode. The PIC modules include multiple particle species, the Lorentz equations of motion, and algorithms for the creation of particles by emission from material surfaces, injection onto the grid, and ionization. In this paper, we present an updated overview of ARGUS, with particular emphasis given in recent algorithmic and computational advances. These include a completely rewritten frequency domain solver which efficiently treats lossy materials and periodic structures, a parallel version of ARGUS with support for both shared memory parallel vector (i.e. CRAY) machines and distributed memory massively parallel MIMD systems, and numerous new applications of the code

The CNCSN: one, two- and three-dimensional coupled neutral and charged particle discrete ordinates code package

International Nuclear Information System (INIS)

Voloschenko, A.M.; Gukov, S.V.; Kryuchkov, V.P.; Dubinin, A.A.; Sumaneev, O.V.

2005-01-01

The CNCSN package is composed of the following codes: -) KATRIN-2.0: a three-dimensional neutral and charged particle transport code; -) KASKAD-S-2.5: a two-dimensional neutral and charged particle transport code; -) ROZ-6.6: a one-dimensional neutral and charged particle transport code; -) ARVES-2.5: a preprocessor for the working macroscopic cross-section format FMAC-M for transport calculations; -) MIXERM: a utility code for preparing mixtures on the base of multigroup cross-section libraries in ANISN format; -) CEPXS-BFP: a version of the Sandia National Lab. multigroup coupled electron-photon cross-section generating code CEPXS, adapted for solving the charged particles transport in the Boltzmann-Fokker-Planck formulation with the use of discrete ordinate method; -) SADCO-2.4: Institute for High-Energy Physics modular system for generating coupled nuclear data libraries to provide high-energy particles transport calculations by multigroup method; -) KATRIF: the post-processor for the KATRIN code; -) KASF: the post-processor for the KASKAD-S code; and ROZ6F: the post-processor for the ROZ-6 code. The coding language is Fortran-90

Three-dimensional electromagnetic solitary waves in an underdense plasma in PIC simulations

International Nuclear Information System (INIS)

Bulanov, Sergei; Esirkepov, Timur; Nishihara, Katsunobu; Pegoraro, Francesco

2002-01-01

A three-dimensional sub-cycle relativistic electromagnetic soliton has been observed for the first time in a 3D Particle-in-Cell simulation of the propagation of an intense short laser pulse in an underdense plasma. The structure of the 3D soliton is identified. It resembles an oscillating electric dipole and has a strong charge separation and toroidal magnetic field component. We call this structure a TM-soliton (transverse magnetic). The 3D TM-soliton resembles a 2D P-soliton in the plane of electric field polarization, and a 2D S-soliton in the perpendicular plane. The core of the soliton is positively charged on average in time, and this results in its Coulomb explosion and in ion heating. Then the soliton evolves into a post-soliton, which is a slowly expanding quasi-neutral cavity in the plasma

Analysis of instability growth and collisionless relaxation in thermionic converters using 1-D PIC simulations

International Nuclear Information System (INIS)

Kreh, B.B.

1994-12-01

This work investigates the role that the beam-plasma instability may play in a thermionic converter. The traditional assumption of collisionally dominated relaxation is questioned, and the beam-plasma instability is proposed as a possible dominant relaxation mechanism. Theory is developed to describe the beam-plasma instability in the cold-plasma approximation, and the theory is tested with two common Particle-in-Cell (PIC) simulation codes. The theory is first confirmed using an unbounded plasma PIC simulation employing periodic boundary conditions, ES1. The theoretically predicted growth rates are on the order of the plasma frequencies, and ES1 simulations verify these predictions within the order of 1%. For typical conditions encountered in thermionic converters, the resulting growth period is on the order of 7 x 10 -11 seconds. The bounded plasma simulation PDP1 was used to evaluate the influence of finite geometry and the electrode boundaries. For this bounded plasma, a two-stream interaction was supported and resulting in nearly complete thermalization in approximately 5 x 10 -10 seconds. Since the electron-electron collision rate of 10 9 Hz and the electron atom collision rate of 10 7 Hz are significantly slower than the rate of development of these instabilities, the instabilities appear to be an important relaxation mechanism

Large Scale Earth's Bow Shock with Northern IMF as Simulated by PIC Code in Parallel with MHD Model

Science.gov (United States)

Baraka, Suleiman

2016-06-01

In this paper, we propose a 3D kinetic model (particle-in-cell, PIC) for the description of the large scale Earth's bow shock. The proposed version is stable and does not require huge or extensive computer resources. Because PIC simulations work with scaled plasma and field parameters, we also propose to validate our code by comparing its results with the available MHD simulations under same scaled solar wind (SW) and (IMF) conditions. We report new results from the two models. In both codes the Earth's bow shock position is found to be ≈14.8 R E along the Sun-Earth line, and ≈29 R E on the dusk side. Those findings are consistent with past in situ observations. Both simulations reproduce the theoretical jump conditions at the shock. However, the PIC code density and temperature distributions are inflated and slightly shifted sunward when compared to the MHD results. Kinetic electron motions and reflected ions upstream may cause this sunward shift. Species distributions in the foreshock region are depicted within the transition of the shock (measured ≈2 c/ ω pi for Θ Bn = 90° and M MS = 4.7) and in the downstream. The size of the foot jump in the magnetic field at the shock is measured to be (1.7 c/ ω pi ). In the foreshocked region, the thermal velocity is found equal to 213 km s-1 at 15 R E and is equal to 63 km s -1 at 12 R E (magnetosheath region). Despite the large cell size of the current version of the PIC code, it is powerful to retain macrostructure of planets magnetospheres in very short time, thus it can be used for pedagogical test purposes. It is also likely complementary with MHD to deepen our understanding of the large scale magnetosphere.

Particle simulation of a two-dimensional electrostatic plasma

International Nuclear Information System (INIS)

Patel, K.

1989-01-01

Computer simulation is a growing field of research and plasma physics is one of the important areas where it is being applied today. This report describes the particle method of simulating a two-dimensional electrostatic plasma. The methods used to discretise the plasma equations and integrate the equations of motion are outlined. The algorithm used in building a simulation program is described. The program is applied to simulating the Two-stream Instability occurring within an infinite plasma. The results of the simulation are presented. The growth rate of the instability as simulated is in excellent agreement with the growth rate as calculated using linear theory. Diagnostic techniques used in interpreting the data generated by the simulation program are discussed. A comparison of the computing environment of the ND and PC from a user's viewpoint is presented. It is observed that the PC is an acceptable computing tool for certain (non-trivial) physics problems, and that more extensive use of its computing power should be made. (author). 5 figs

Particle Demagnetization in Collisionless Magnetic Reconnection

Science.gov (United States)

Hesse, Michael

2006-01-01

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Analysis of effects of laser profiles on fast electron generation by two-dimensional Particle-In-Cell simulations

International Nuclear Information System (INIS)

Hata, M.

2010-01-01

Complete text of publication follows. A cone-guided target is used in the Fast Ignition Realization Experiment project phase-I (FIREX-I) and optimization of its design is performed. However a laser profile is not optimized much, because the laser profile that is the best for core heating is not known well. To find that, it is useful to investigate characteristics of generated fast electrons in each condition of different laser profiles. In this research, effects of laser profiles on fast electron generation are investigated on somewhat simple conditions by two-dimensional Particle-In-Cell simulations. In these simulations, a target is made up of Au pre-plasma and Au plasma. The Au pre-plasma has the exponential profile in the x direction with the scale length L = 4.0 μm and the density from 0.10 n cr to 20 n cr . The Au plasma has the flat profile in the x direction with 10 μm width and 20 n cr . Plasma profiles are uniform in the y direction. The ionization degree and the mass number of plasmas are 40 and 197, where the ionization degree is determined by PINOCO simulations. PINOCO is a two-dimensional radiation hydrodynamics simulation code, which simulates formation of the high-density plasma during the compression phase in the fast ignition. A laser is assumed to propagate as plane wave from the negative x direction to the positive x direction. Laser profiles are supposed to be uniform in the y direction. Three different laser profiles, namely flat one with t flat = 100 fs, Gaussian one with t rise/fall = 47.0 fs and flat + Gaussian one with t rise/fall = 23.5 fs and t flat = 50 fs are used. The energy and the peak intensity are constant with E = 10 7 J/cm 2 and I L = 10 20 W/cm 2 in all cases of different laser profiles. We compare results in each condition of three different laser profiles and investigate effects of laser profiles on fast electron generation. Time-integrated energy spectra are similar in all cases of three different laser profiles. In the

Collision probability in two-dimensional lattice by ray-trace method and its applications to cell calculations

International Nuclear Information System (INIS)

Tsuchihashi, Keichiro

1985-03-01

A series of formulations to evaluate collision probability for multi-region cells expressed by either of three one-dimensional coordinate systems (plane, sphere and cylinder) or by the general two-dimensional cylindrical coordinate system is presented. They are expressed in a suitable form to have a common numerical process named ''Ray-Trace'' method. Applications of the collision probability method to two optional treatments for the resonance absorption are presented. One is a modified table-look-up method based on the intermediate resonance approximation, and the other is a rigorous method to calculate the resonance absorption in a multi-region cell in which nearly continuous energy spectra of the resonance neutron range can be solved and interaction effect between different resonance nuclides can be evaluated. Two works on resonance absorption in a doubly heterogeneous system with grain structure are presented. First, the effect of a random distribution of particles embedded in graphite diluent on the resonance integral is studied. Next, the ''Accretion'' method proposed by Leslie and Jonsson to define the collision probability in a doubly heterogeneous system is applied to evaluate the resonance absorption in coated particles dispersed in fuel pellet of the HTGR. Several optional models are proposed to define the collision rates in the medium with the microscopic heterogeneity. By making use of the collision probability method developed by the present study, the JAERI thermal reactor standard nuclear design code system SRAC has been developed. Results of several benchmark tests for the SRAC are presented. The analyses of critical experiments of the SHE, DCA, and FNR show good agreement of critical masses with their experimental values. (J.P.N.)

PIC 16 F84

International Nuclear Information System (INIS)

Jung, Gi Cheol; Min, Han Sik

2001-11-01

The contents of this book are introduction of microprocessor, basic for microcomputer practice, introduction of one chip micro computer, basic command of PIC, instructions simulator and in circuit emulator ; what a simulator of PIC is, and MPLAB direction, making PIC rom writer and instructions, of La's PIC Micro Programmer, PIC programming ; learning Command with examples, and controlling hardware with C-language, practical task for PIC application, a line tracer automobile and making ultrasonic radar ; circuit, source program and monitor program.

PIC 16 F84

Energy Technology Data Exchange (ETDEWEB)

Jung, Gi Cheol; Min, Han Sik

2001-11-15

The contents of this book are introduction of microprocessor, basic for microcomputer practice, introduction of one chip micro computer, basic command of PIC, instructions simulator and in circuit emulator ; what a simulator of PIC is, and MPLAB direction, making PIC rom writer and instructions, of La's PIC Micro Programmer, PIC programming ; learning Command with examples, and controlling hardware with C-language, practical task for PIC application, a line tracer automobile and making ultrasonic radar ; circuit, source program and monitor program.

PIC microcomputer guide for beginner

International Nuclear Information System (INIS)

Shin, Chulho

2001-03-01

This book comprised of four parts. The first part deals with computer one chip, voltage current, resistance, electronic components, logical element, TTL and CMOS, memory and I/O and MDS. The second part is about PIC16C84 which describes its memory structure, registers and PIC16C84 command. The third part deals with LED control program, jet car LED, quiz buzzer program, LED spectrum, digital dice, two digital dices and time bomb. The last part introduces PIC16C71 and temperature controller.

Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations

Energy Technology Data Exchange (ETDEWEB)

Qin, Hong; Liu, Jian; Xiao, Jianyuan; Zhang, Ruili; He, Yang; Wang, Yulei; Sun, Yajuan; Burby, Joshua W.; Ellison, Leland; Zhou, Yao

2015-12-14

Particle-in-cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretising its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root searching or global matrix inversion, enabling applications of the proposed method to very large-scale plasma simulations with many, e.g. 10(9), degrees of freedom. The long-term accuracy and fidelity of the algorithm enables us to numerically confirm Mouhot and Villani's theory and conjecture on nonlinear Landau damping over several orders of magnitude using the PIC method, and to calculate the nonlinear evolution of the reflectivity during the mode conversion process from extraordinary waves to Bernstein waves.

A parallel code named NEPTUNE for 3D fully electromagnetic and pic simulations

International Nuclear Information System (INIS)

Dong Ye; Yang Wenyuan; Chen Jun; Zhao Qiang; Xia Fang; Ma Yan; Xiao Li; Sun Huifang; Chen Hong; Zhou Haijing; Mao Zeyao; Dong Zhiwei

2010-01-01

A parallel code named NEPTUNE for 3D fully electromagnetic and particle-in-cell (PIC) simulations is introduced, which could run on the Linux system with hundreds to thousand CPUs. NEPTUNE is suitable to simulate entire 3D HPM devices; many HPM devices are simulated and designed by using it. In NEPTUNE code, the electromagnetic fields are updated by using the finite-difference in time domain (FDTD) method of solving Maxwell equations and the particles are moved by using Buneman-Boris advance method of solving relativistic Newton-Lorentz equation. Electromagnetic fields and particles are coupled by using liner weighing interpolation PIC method, and the electric filed components are corrected by using Boris method of solve Poisson equation in order to ensure charge-conservation. NEPTUNE code could construct many complicated geometric structures, such as arbitrary axial-symmetric structures, plane transforming structures, slow-wave-structures, coupling holes, foils, and so on. The boundary conditions used in NEPTUNE code are introduced in brief, including perfectly electric conductor boundary, external wave boundary, and particle boundary. Finally, some typical HPM devices are simulated and test by using NEPTUNE code, including MILO, RBWO, VCO, and RKA. The simulation results are with correct and credible physical images, and the parallel efficiencies are also given. (authors)

Programming 16-Bit PIC Microcontrollers in C Learning to Fly the PIC 24

CERN Document Server

Di Jasio, Lucio

2011-01-01

New in the second edition: * MPLAB X support and MPLAB C for the PIC24F v3 and later libraries * I2C™ interface * 100% assembly free solutions * Improved video, PAL/NTSC * Improved audio, RIFF files decoding * PIC24F GA1, GA2, GB1 and GB2 support   Most readers will associate Microchip's name with the ubiquitous 8-bit PIC microcontrollers but it is the new 16-bit PIC24F family that is truly stealing the scene. Orders of magnitude increases of performance, memory size and the rich peripheral set make programming these devices in C a must. This new guide by Microchip insid

A One-Dimensional Particle-in-Cell Model of Plasma Build-Up in Vacuum Arcs

CERN Document Server

Timko, H; Kovermann, J; Taborelli, M; Nordlund, K; Descoeudres, A; Schneider, R; Calatroni, S; Matyash, K; Wuensch, W; Hansen, A; Grudiev, A

2011-01-01

Understanding the mechanism of plasma build-up in vacuum arcs is essential in many fields of physics. A one-dimensional particle-in-cell computer simulation model is presented, which models the plasma developing from a field emitter tip under electrical breakdown conditions, taking into account the relevant physical phenomena. As a starting point, only an external electric field and an initial enhancement factor of the tip are assumed. General requirements for plasma formation have been identified and formulated in terms of the initial local field and a critical neutral density. The dependence of plasma build-up on tip melting current, the evaporation rate of neutrals and external circuit time constant has been investigated for copper and simulations imply that arcing involves melting currents around 0.5-1 A/mu m(2),evaporation of neutrals to electron field emission ratios in the regime 0.01 - 0.05, plasma build-up timescales in the order of similar to 1 - 10 ns and two different regimes depending on initial ...

Two-dimensional beam profiles and one-dimensional projections

Science.gov (United States)

Findlay, D. J. S.; Jones, B.; Adams, D. J.

2018-05-01

One-dimensional projections of improved two-dimensional representations of transverse profiles of particle beams are proposed for fitting to data from harp-type monitors measuring beam profiles on particle accelerators. Composite distributions, with tails smoothly matched on to a central (inverted) parabola, are shown to give noticeably better fits than single gaussian and single parabolic distributions to data from harp-type beam profile monitors all along the proton beam transport lines to the two target stations on the ISIS Spallation Neutron Source. Some implications for inferring beam current densities on the beam axis are noted.

Role of positive ions on the surface production of negative ions in a fusion plasma reactor type negative ion source--Insights from a three dimensional particle-in-cell Monte Carlo collisions model

Science.gov (United States)

Fubiani, G.; Boeuf, J. P.

2013-11-01

Results from a 3D self-consistent Particle-In-Cell Monte Carlo Collisions (PIC MCC) model of a high power fusion-type negative ion source are presented for the first time. The model is used to calculate the plasma characteristics of the ITER prototype BATMAN ion source developed in Garching. Special emphasis is put on the production of negative ions on the plasma grid surface. The question of the relative roles of the impact of neutral hydrogen atoms and positive ions on the cesiated grid surface has attracted much attention recently and the 3D PIC MCC model is used to address this question. The results show that the production of negative ions by positive ion impact on the plasma grid is small with respect to the production by atomic hydrogen or deuterium bombardment (less than 10%).

Role of positive ions on the surface production of negative ions in a fusion plasma reactor type negative ion source—Insights from a three dimensional particle-in-cell Monte Carlo collisions model

International Nuclear Information System (INIS)

Fubiani, G.; Boeuf, J. P.

2013-01-01

Results from a 3D self-consistent Particle-In-Cell Monte Carlo Collisions (PIC MCC) model of a high power fusion-type negative ion source are presented for the first time. The model is used to calculate the plasma characteristics of the ITER prototype BATMAN ion source developed in Garching. Special emphasis is put on the production of negative ions on the plasma grid surface. The question of the relative roles of the impact of neutral hydrogen atoms and positive ions on the cesiated grid surface has attracted much attention recently and the 3D PIC MCC model is used to address this question. The results show that the production of negative ions by positive ion impact on the plasma grid is small with respect to the production by atomic hydrogen or deuterium bombardment (less than 10%)

A parallel 3D particle-in-cell code with dynamic load balancing

International Nuclear Information System (INIS)

Wolfheimer, Felix; Gjonaj, Erion; Weiland, Thomas

2006-01-01

A parallel 3D electrostatic Particle-In-Cell (PIC) code including an algorithm for modelling Space Charge Limited (SCL) emission [E. Gjonaj, T. Weiland, 3D-modeling of space-charge-limited electron emission. A charge conserving algorithm, Proceedings of the 11th Biennial IEEE Conference on Electromagnetic Field Computation, 2004] is presented. A domain decomposition technique based on orthogonal recursive bisection is used to parallelize the computation on a distributed memory environment of clustered workstations. For problems with a highly nonuniform and time dependent distribution of particles, e.g., bunch dynamics, a dynamic load balancing between the processes is needed to preserve the parallel performance. The algorithm for the detection of a load imbalance and the redistribution of the tasks among the processes is based on a weight function criterion, where the weight of a cell measures the computational load associated with it. The algorithm is studied with two examples. In the first example, multiple electron bunches as occurring in the S-DALINAC [A. Richter, Operational experience at the S-DALINAC, Proceedings of the Fifth European Particle Accelerator Conference, 1996] accelerator are simulated in the absence of space charge fields. In the second example, the SCL emission and electron trajectories in an electron gun are simulated

A parallel 3D particle-in-cell code with dynamic load balancing

Energy Technology Data Exchange (ETDEWEB)

Wolfheimer, Felix [Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Schlossgartenstr.8, 64283 Darmstadt (Germany)]. E-mail: wolfheimer@temf.de; Gjonaj, Erion [Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Schlossgartenstr.8, 64283 Darmstadt (Germany); Weiland, Thomas [Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Schlossgartenstr.8, 64283 Darmstadt (Germany)

2006-03-01

A parallel 3D electrostatic Particle-In-Cell (PIC) code including an algorithm for modelling Space Charge Limited (SCL) emission [E. Gjonaj, T. Weiland, 3D-modeling of space-charge-limited electron emission. A charge conserving algorithm, Proceedings of the 11th Biennial IEEE Conference on Electromagnetic Field Computation, 2004] is presented. A domain decomposition technique based on orthogonal recursive bisection is used to parallelize the computation on a distributed memory environment of clustered workstations. For problems with a highly nonuniform and time dependent distribution of particles, e.g., bunch dynamics, a dynamic load balancing between the processes is needed to preserve the parallel performance. The algorithm for the detection of a load imbalance and the redistribution of the tasks among the processes is based on a weight function criterion, where the weight of a cell measures the computational load associated with it. The algorithm is studied with two examples. In the first example, multiple electron bunches as occurring in the S-DALINAC [A. Richter, Operational experience at the S-DALINAC, Proceedings of the Fifth European Particle Accelerator Conference, 1996] accelerator are simulated in the absence of space charge fields. In the second example, the SCL emission and electron trajectories in an electron gun are simulated.

Particle in cell simulation of peaking switch for breakdown evaluation

Energy Technology Data Exchange (ETDEWEB)

Umbarkar, Sachin B.; Bindu, S.; Mangalvedekar, H.A.; Saxena, A.; Singh, N.M., E-mail: sachin.b.umbarkar@gmail.com [Department of Electric Engineering, Veermata Jijabai Technological Institute, Mumbai (India); Sharma, Archana; Saroj, P.C.; Mittal, K.C. [Accelerator Pulse Power Division, Bhabha Atomic Research Centre, Mumbai (India)

2014-07-01

Marx generator connected to peaking capacitor and peaking switch can generate Ultra-Wideband (UWB) radiation. A new peaking switch is designed for converting the existing nanosecond Marx generator to a UWB source. The paper explains the particle in cell (PIC) simulation for this peaking switch, using MAGIC 3D software. This peaking switch electrode is made up of copper tungsten material and is fixed inside the hermitically sealed derlin material. The switch can withstand a gas pressure up to 13.5 kg/cm{sup 2}. The lower electrode of the switch is connected to the last stage of the Marx generator. Initially Marx generator (without peaking stage) in air; gives the output pulse with peak amplitude of 113.75 kV and pulse rise time of 25 ns. Thus, we design a new peaking switch to improve the rise time of output pulse and to pressurize this peaking switch separately (i.e. Marx and peaking switch is at different pressure). The PIC simulation gives the particle charge density, current density, E counter plot, emitted electron current, and particle energy along the axis of gap between electrodes. The charge injection and electric field dependence on ionic dissociation phenomenon are briefly analyzed using this simulation. The model is simulated with different gases (N{sub 2}, H{sub 2}, and Air) under different pressure (2 kg/cm{sup 2}, 5 kg/cm{sup 2}, 10 kg/cm{sup 2}). (author)

Statistics, distillation, and ordering emergence in a two-dimensional stochastic model of particles in counterflowing streams

Science.gov (United States)

Stock, Eduardo Velasco; da Silva, Roberto; Fernandes, H. A.

2017-07-01

In this paper, we propose a stochastic model which describes two species of particles moving in counterflow. The model generalizes the theoretical framework that describes the transport in random systems by taking into account two different scenarios: particles can work as mobile obstacles, whereas particles of one species move in the opposite direction to the particles of the other species, or particles of a given species work as fixed obstacles remaining in their places during the time evolution. We conduct a detailed study about the statistics concerning the crossing time of particles, as well as the effects of the lateral transitions on the time required to the system reaches a state of complete geographic separation of species. The spatial effects of jamming are also studied by looking into the deformation of the concentration of particles in the two-dimensional corridor. Finally, we observe in our study the formation of patterns of lanes which reach the steady state regardless of the initial conditions used for the evolution. A similar result is also observed in real experiments involving charged colloids motion and simulations of pedestrian dynamics based on Langevin equations, when periodic boundary conditions are considered (particles counterflow in a ring symmetry). The results obtained through Monte Carlo simulations and numerical integrations are in good agreement with each other. However, differently from previous studies, the dynamics considered in this work is not Newton-based, and therefore, even artificial situations of self-propelled objects should be studied in this first-principles modeling.

Electron acceleration in the Solar corona - 3D PiC code simulations of guide field reconnection

Science.gov (United States)

Alejandro Munoz Sepulveda, Patricio

2017-04-01

The efficient electron acceleration in the solar corona detected by means of hard X-ray emission is still not well understood. Magnetic reconnection through current sheets is one of the proposed production mechanisms of non-thermal electrons in solar flares. Previous works in this direction were based mostly on test particle calculations or 2D fully-kinetic PiC simulations. We have now studied the consequences of self-generated current-aligned instabilities on the electron acceleration mechanisms by 3D magnetic reconnection. For this sake, we carried out 3D Particle-in-Cell (PiC) code numerical simulations of force free reconnecting current sheets, appropriate for the description of the solar coronal plasmas. We find an efficient electron energization, evidenced by the formation of a non-thermal power-law tail with a hard spectral index smaller than -2 in the electron energy distribution function. We discuss and compare the influence of the parallel electric field versus the curvature and gradient drifts in the guiding-center approximation on the overall acceleration, and their dependence on different plasma parameters.

Cell-centered particle weighting algorithm for PIC simulations in a non-uniform 2D axisymmetric mesh

Science.gov (United States)

Araki, Samuel J.; Wirz, Richard E.

2014-09-01

Standard area weighting methods for particle-in-cell simulations result in systematic errors on particle densities for a non-uniform mesh in cylindrical coordinates. These errors can be significantly reduced by using weighted cell volumes for density calculations. A detailed description on the corrected volume calculations and cell-centered weighting algorithm in a non-uniform mesh is provided. The simple formulas for the corrected volume can be used for any type of quadrilateral and/or triangular mesh in cylindrical coordinates. Density errors arising from the cell-centered weighting algorithm are computed for radial density profiles of uniform, linearly decreasing, and Bessel function in an adaptive Cartesian mesh and an unstructured mesh. For all the density profiles, it is shown that the weighting algorithm provides a significant improvement for density calculations. However, relatively large density errors may persist at outermost cells for monotonically decreasing density profiles. A further analysis has been performed to investigate the effect of the density errors in potential calculations, and it is shown that the error at the outermost cell does not propagate into the potential solution for the density profiles investigated.

Two-dimensional simulation of sintering process

International Nuclear Information System (INIS)

Vasconcelos, Vanderley de; Pinto, Lucio Carlos Martins; Vasconcelos, Wander L.

1996-01-01

The results of two-dimensional simulations are directly applied to systems in which one of the dimensions is much smaller than the others, and to sections of three dimensional models. Moreover, these simulations are the first step of the analysis of more complex three-dimensional systems. In this work, two basic features of the sintering process are studied: the types of particle size distributions related to the powder production processes and the evolution of geometric parameters of the resultant microstructures during the solid-state sintering. Random packing of equal spheres is considered in the sintering simulation. The packing algorithm does not take into account the interactive forces between the particles. The used sintering algorithm causes the densification of the particle set. (author)

Cold-Fluid Equilibrium of a Large-Aspect-Ratio Ellipse-Shaped Charged-Particle Beam in a Non-Axisymmetric Periodic Permanent Magnet Focusing Field

CERN Document Server

Zhou, Jing; Chen Chi Ping

2005-01-01

A new class of equilibrium is discovered for a large-aspect-ratio ellipse-shaped charged-particle beam in a non-axisymmetric periodic permanent magnet focusing field. A paraxial cold-fluid model is employed to derive the equilibrium flow properties and generalized envelope equations with negligibly small emittance. A periodic beam equilibrium solution is obtained numerically from the generalized envelope equations. It is shown that the beam edges are well confined in both transverse directions, and that the equilibrium beam exhibits a small-angle periodic wobble as it propagates. A two-dimensional particle-in-cell (PIC) code, PFB2D, is used to verify the theoretical predictions in the paraxial limit, and to establish validity under non-paraxial situations and the influence of the conductor walls of the beam tunnel.

Two dimensional kinetic analysis of electrostatic harmonic plasma waves

Energy Technology Data Exchange (ETDEWEB)

Fonseca-Pongutá, E. C.; Ziebell, L. F.; Gaelzer, R. [Instituto de Física, UFRGS, 91501-970 Porto Alegre, RS (Brazil); Yoon, P. H. [IPST, University of Maryland, College Park, Maryland 20742 (United States); SSR, Kyung Hee University, Yongin, Gyeonggi 446-701 (Korea, Republic of)

2016-06-15

Electrostatic harmonic Langmuir waves are virtual modes excited in weakly turbulent plasmas, first observed in early laboratory beam-plasma experiments as well as in rocket-borne active experiments in space. However, their unequivocal presence was confirmed through computer simulated experiments and subsequently theoretically explained. The peculiarity of harmonic Langmuir waves is that while their existence requires nonlinear response, their excitation mechanism and subsequent early time evolution are governed by essentially linear process. One of the unresolved theoretical issues regards the role of nonlinear wave-particle interaction process over longer evolution time period. Another outstanding issue is that existing theories for these modes are limited to one-dimensional space. The present paper carries out two dimensional theoretical analysis of fundamental and (first) harmonic Langmuir waves for the first time. The result shows that harmonic Langmuir wave is essentially governed by (quasi)linear process and that nonlinear wave-particle interaction plays no significant role in the time evolution of the wave spectrum. The numerical solutions of the two-dimensional wave spectra for fundamental and harmonic Langmuir waves are also found to be consistent with those obtained by direct particle-in-cell simulation method reported in the literature.

Noise-induced drift in two-dimensional anisotropic systems

Science.gov (United States)

Farago, Oded

2017-10-01

We study the isothermal Brownian dynamics of a particle in a system with spatially varying diffusivity. Due to the heterogeneity of the system, the particle's mean displacement does not vanish even if it does not experience any physical force. This phenomenon has been termed "noise-induced drift," and has been extensively studied for one-dimensional systems. Here, we examine the noise-induced drift in a two-dimensional anisotropic system, characterized by a symmetric diffusion tensor with unequal diagonal elements. A general expression for the mean displacement vector is derived and presented as a sum of two vectors, depicting two distinct drifting effects. The first vector describes the tendency of the particle to drift toward the high diffusivity side in each orthogonal principal diffusion direction. This is a generalization of the well-known expression for the noise-induced drift in one-dimensional systems. The second vector represents a novel drifting effect, not found in one-dimensional systems, originating from the spatial rotation in the directions of the principal axes. The validity of the derived expressions is verified by using Langevin dynamics simulations. As a specific example, we consider the relative diffusion of two transmembrane proteins, and demonstrate that the average distance between them increases at a surprisingly fast rate of several tens of micrometers per second.

Effect of Substrate Friction in a Two-Dimensional Granular Couette Shearing Cell

Science.gov (United States)

Templeman, Chris; Garg, Shila

2001-03-01

An investigation of the effect of substrate friction on the kinematics of rigid granular material in a two-dimensional granular Couette shearing cell was conducted. Cylindrical disks resting on a substrate were packed between a stationary outer ring and a rotating inner wheel. Previous work reports the velocity and particle rotation rates as a function of packing fraction and shearing rates [1]. The authors report the existence of a stick-slip condition of the disks in contact with the shearing wheel. The focus of our study is to investigate the impact of the substrate friction on the stick-slip condition as well as the kinematics of the system in general. [1] C.T. Veje, Daniel W. Howell, and R.P Behringer, Phys. Rev. E 59, 739 (1999). This research was partially supported by the Copeland Fund, administered by The College of Wooster. C.T. received support from NASA GRC LERCIP internship program.

Two-Dimensional IIR Filter Design Using Simulated Annealing Based Particle Swarm Optimization

Directory of Open Access Journals (Sweden)

Supriya Dhabal

2014-01-01

Full Text Available We present a novel hybrid algorithm based on particle swarm optimization (PSO and simulated annealing (SA for the design of two-dimensional recursive digital filters. The proposed method, known as SA-PSO, integrates the global search ability of PSO with the local search ability of SA and offsets the weakness of each other. The acceptance criterion of Metropolis is included in the basic algorithm of PSO to increase the swarm’s diversity by accepting sometimes weaker solutions also. The experimental results reveal that the performance of the optimal filter designed by the proposed SA-PSO method is improved. Further, the convergence behavior as well as optimization accuracy of proposed method has been improved significantly and computational time is also reduced. In addition, the proposed SA-PSO method also produces the best optimal solution with lower mean and variance which indicates that the algorithm can be used more efficiently in realizing two-dimensional digital filters.

A parallel implementation of particle tracking with space charge effects on an INTEL iPSC/860

International Nuclear Information System (INIS)

Chang, L.; Bourianoff, G.; Cole, B.; Machida, S.

1993-05-01

Particle-tracking simulation is one of the scientific applications that is well-suited to parallel computations. At the Superconducting Super Collider, it has been theoretically and empirically demonstrated that particle tracking on a designed lattice can achieve very high parallel efficiency on a MIMD Intel iPSC/860 machine. The key to such success is the realization that the particles can be tracked independently without considering their interaction. The perfectly parallel nature of particle tracking is broken if the interaction effects between particles are included. The space charge introduces an electromagnetic force that will affect the motion of tracked particles in 3-D space. For accurate modeling of the beam dynamics with space charge effects, one needs to solve three-dimensional Maxwell field equations, usually by a particle-in-cell (PIC) algorithm. This will require each particle to communicate with its neighbor grids to compute the momentum changes at each time step. It is expected that the 3-D PIC method will degrade parallel efficiency of particle-tracking implementation on any parallel computer. In this paper, we describe an efficient scheme for implementing particle tracking with space charge effects on an INTEL iPSC/860 machine. Experimental results show that a parallel efficiency of 75% can be obtained

The effects of particle recycling on the divertor plasma: A particle-in-cell with Monte Carlo collision simulation

Science.gov (United States)

Chang, Mingyu; Sang, Chaofeng; Sun, Zhenyue; Hu, Wanpeng; Wang, Dezhen

2018-05-01

A Particle-In-Cell (PIC) with Monte Carlo Collision (MCC) model is applied to study the effects of particle recycling on divertor plasma in the present work. The simulation domain is the scrape-off layer of the tokamak in one-dimension along the magnetic field line. At the divertor plate, the reflected deuterium atoms (D) and thermally released deuterium molecules (D2) are considered. The collisions between the plasma particles (e and D+) and recycled neutral particles (D and D2) are described by the MCC method. It is found that the recycled neutral particles have a great impact on divertor plasma. The effects of different collisions on the plasma are simulated and discussed. Moreover, the impacts of target materials on the plasma are simulated by comparing the divertor with Carbon (C) and Tungsten (W) targets. The simulation results show that the energy and momentum losses of the C target are larger than those of the W target in the divertor region even without considering the impurity particles, whereas the W target has a more remarkable influence on the core plasma.

Resonant and non-resonant whistlers-particle interaction in the radiation belts

NARCIS (Netherlands)

E. Camporeale (Enrico)

2014-01-01

htmlabstractWe study the wave-particle interactions between lower band chorus whistlers and an anisotropic tenuous population of relativistic electrons. We present the first direct comparison of first-principle Particle-in-Cell (PIC) simulations with a quasi-linear diffusion code, in this context.

Real-Time Two-Dimensional Magnetic Particle Imaging for Electromagnetic Navigation in Targeted Drug Delivery

Science.gov (United States)

Le, Tuan-Anh; Zhang, Xingming; Hoshiar, Ali Kafash; Yoon, Jungwon

2017-01-01

Magnetic nanoparticles (MNPs) are effective drug carriers. By using electromagnetic actuated systems, MNPs can be controlled noninvasively in a vascular network for targeted drug delivery (TDD). Although drugs can reach their target location through capturing schemes of MNPs by permanent magnets, drugs delivered to non-target regions can affect healthy tissues and cause undesirable side effects. Real-time monitoring of MNPs can improve the targeting efficiency of TDD systems. In this paper, a two-dimensional (2D) real-time monitoring scheme has been developed for an MNP guidance system. Resovist particles 45 to 65 nm in diameter (5 nm core) can be monitored in real-time (update rate = 2 Hz) in 2D. The proposed 2D monitoring system allows dynamic tracking of MNPs during TDD and renders magnetic particle imaging-based navigation more feasible. PMID:28880220

AP-Cloud: Adaptive Particle-in-Cloud method for optimal solutions to Vlasov–Poisson equation

International Nuclear Information System (INIS)

Wang, Xingyu; Samulyak, Roman; Jiao, Xiangmin; Yu, Kwangmin

2016-01-01

We propose a new adaptive Particle-in-Cloud (AP-Cloud) method for obtaining optimal numerical solutions to the Vlasov–Poisson equation. Unlike the traditional particle-in-cell (PIC) method, which is commonly used for solving this problem, the AP-Cloud adaptively selects computational nodes or particles to deliver higher accuracy and efficiency when the particle distribution is highly non-uniform. Unlike other adaptive techniques for PIC, our method balances the errors in PDE discretization and Monte Carlo integration, and discretizes the differential operators using a generalized finite difference (GFD) method based on a weighted least square formulation. As a result, AP-Cloud is independent of the geometric shapes of computational domains and is free of artificial parameters. Efficient and robust implementation is achieved through an octree data structure with 2:1 balance. We analyze the accuracy and convergence order of AP-Cloud theoretically, and verify the method using an electrostatic problem of a particle beam with halo. Simulation results show that the AP-Cloud method is substantially more accurate and faster than the traditional PIC, and it is free of artificial forces that are typical for some adaptive PIC techniques.

Effects of Temperature and Residence Time on the Emissions of PIC and Fine Particles during Fixed Bed Combustion of Conifer Stemwood Pellets

Energy Technology Data Exchange (ETDEWEB)

Boman, Christoffer; Lindmark, Fredrik; Oehman, Marcus; Nordin, Anders [Umeaa Univ. (Sweden). Energy Technology and Thermal Process Chemistry; Pettersson, Esbjoern [Energy Technology Centre, Piteaa (Sweden); Westerholm, Roger [Stockholm Univ., Arrhenius Laboratory (Sweden). Dept. of Analytical Chemistry

2006-07-15

The use of wood fuel Pellets has proved to be well suited for the small-scale market enabling controlled and efficient combustion with low emission of products of incomplete combustion (PIC). Still a potential for further emission reduction exists and a thorough understanding of the influence of combustion conditions on the emission characteristics of air pollutants like PAH and particulate matter (PM) is important. The objective was to determine the effects of temperature and residence time on the emission performance and characteristics with focus on hydrocarbons and PM during combustion of conifer stemwood Pellets in a laboratory fixed bed reactor (<5 kW). Temperature and residence time after the bed section were varied according to statistical experimental designs (650-970 deg C and 0.5-3.5 s) with the emission responses; CO, organic gaseous carbon, NO, 20 VOC compounds, 43 PAH compounds, PM{sub tot}, fine particle mass/count median diameter (MMD and CMD) and number concentration. Temperature was negatively correlated with the emissions of all studied PIC with limited effects of residence time. The PM{sub tot} emissions of 15-20 mg/MJ was in all cases dominated by fine (<1 {mu}m) particles of K, Na, S, Cl, C, O and Zn. Increased residence time resulted in increased fine particle sizes (i.e. MMD and CMD) and decreased number concentrations. The importance of high temperature (>850 deg C) in the bed zone with intensive, air rich and well mixed isothermal conditions for 0.5-1.0 s in the post combustion zone was illustrated for wood Pellets combustion with almost a total depletion of all studied PIC. The results emphasize the need for further verification studies and technology development work.

Three-dimensional particle image velocimetry measurement technique

International Nuclear Information System (INIS)

Hassan, Y.A.; Seeley, C.H.; Henderson, J.A.; Schmidl, W.D.

2004-01-01

The experimental flow visualization tool, Particle Image Velocimetry (PIV), is being used to determine the velocity field in two-dimensional fluid flows. In the past few years, the technique has been improved to allow the capture of flow fields in three dimensions. This paper describes changes which were made to two existing two-dimensional tracking algorithms to enable them to track three-dimensional PIV data. Results of the tests performed on these three-dimensional routines with synthetic data are presented. Experimental data was also used to test the tracking algorithms. The test setup which was used to acquire the three-dimensional experimental data is described, along with the results from both of the tracking routines which were used to analyze the experimental data. (author)

Measurement of correlations between low-frequency vibrational modes and particle rearrangements in quasi-two-dimensional colloidal glasses.

Science.gov (United States)

Chen, Ke; Manning, M L; Yunker, Peter J; Ellenbroek, Wouter G; Zhang, Zexin; Liu, Andrea J; Yodh, A G

2011-09-02

We investigate correlations between low-frequency vibrational modes and rearrangements in two-dimensional colloidal glasses composed of thermosensitive microgel particles, which readily permit variation of the sample packing fraction. At each packing fraction, the particle displacement covariance matrix is measured and used to extract the vibrational spectrum of the "shadow" colloidal glass (i.e., the particle network with the same geometry and interactions as the sample colloid but absent damping). Rearrangements are induced by successive, small reductions in the packing fraction. The experimental results suggest that low-frequency quasilocalized phonon modes in colloidal glasses, i.e., modes that present low energy barriers for system rearrangements, are spatially correlated with rearrangements in this thermal system.

Model independent search for new particles in two-dimensional mass space using events with missing energy, two jets and two leptons with the CMS detector

CERN Document Server

AUTHOR|(CDS)2080070; Hebbeker, Thomas

2017-07-07

The discovery of a new particle consistent with the standard model Higgs boson at the Large Hadron Collider in 2012 completed the standard model of particle physics (SM). Despite its remarkable success many questions remain unexplained. Numerous theoretical models, predicting the existence of new heavy particles, provide answers to these unresolved questions and are tested at high energy experiments such as the Compact Muon Solenoid (CMS) detector at the Large Hadron Collider (LHC). In this thesis a model independent search method for new particles in two-dimensional mass space in events with missing transverse energy is presented using 19.7 $\\mbox{fb}^{-1}$ of proton-proton collision data recorded by the CMS detector at a centre of mass energy $\\sqrt{s}$ = 8 TeV at the LHC. The analysis searches for signatures of pair-produced new heavy particles $\\mbox{T}^\\prime$ which decay further into unknown heavy particles $\\mbox{W}^\\prime$ and SM quarks $q$ ($\\mbox{T}^\\prime\\overline{\\mbox{T}^\\prime} \\rightarrow {...

Anomalous electron transport in Hall-effect thrusters: Comparison between quasi-linear kinetic theory and particle-in-cell simulations

Science.gov (United States)

Lafleur, T.; Martorelli, R.; Chabert, P.; Bourdon, A.

2018-06-01

Kinetic drift instabilities have been implicated as a possible mechanism leading to anomalous electron cross-field transport in E × B discharges, such as Hall-effect thrusters. Such instabilities, which are driven by the large disparity in electron and ion drift velocities, present a significant challenge to modelling efforts without resorting to time-consuming particle-in-cell (PIC) simulations. Here, we test aspects of quasi-linear kinetic theory with 2D PIC simulations with the aim of developing a self-consistent treatment of these instabilities. The specific quantities of interest are the instability growth rate (which determines the spatial and temporal evolution of the instability amplitude), and the instability-enhanced electron-ion friction force (which leads to "anomalous" electron transport). By using the self-consistently obtained electron distribution functions from the PIC simulations (which are in general non-Maxwellian), we find that the predictions of the quasi-linear kinetic theory are in good agreement with the simulation results. By contrast, the use of Maxwellian distributions leads to a growth rate and electron-ion friction force that is around 2-4 times higher, and consequently significantly overestimates the electron transport. A possible method for self-consistently modelling the distribution functions without requiring PIC simulations is discussed.

One and two dimensional simulations on beat wave acceleration

International Nuclear Information System (INIS)

Mori, W.; Joshi, C.; Dawson, J.M.; Forslund, D.W.; Kindel, J.M.

1984-01-01

Recently there has been considerable interest in the use of fast-large-amplitude plasma waves as the basis for a high energy particle accelerator. In these schemes, lasers are used to create the plasma wave. To date the few simulation studies on this subject have been limited to one-dimensional, short rise time simulations. Here the authors present results from simulations in which more realistic parameters are used. In addition, they present the first two dimensional simulations on this subject. One dimensional simulations on a 2 1/2-D relativistic electromagnetic particle code, in which only a few cells were used in one direction, on colinear optical mixing are presented. In these simulations the laser rise time, laser intensity, plasma density, plasma temperature and system size were varied. The simulations indicate that the theory of Rosenbluth and Liu is applicable over a wide range of parameters. In addition, simulations with a DC magnetic field are presented in order to study the ''Surfatron'' concept

Correlated particle dynamics in concentrated quasi-two-dimensional suspensions

International Nuclear Information System (INIS)

Diamant, H; Cui, B; Lin, B; Rice, S A

2005-01-01

We investigate theoretically and experimentally how the hydrodynamically correlated lateral motion of particles in a suspension confined between two surfaces is affected by the suspension concentration. Despite the long range of the correlations (decaying as 1/r 2 with the inter-particle distance r), the concentration effect is present only at short inter-particle distances for which the static pair correlation is nonuniform. This is in sharp contrast with the effect of hydrodynamic screening in unconfined suspensions, where increasing the concentration changes the prefactor of the large-distance correlation

Plasma and BIAS Modeling: Self-Consistent Electrostatic Particle-in-Cell with Low-Density Argon Plasma for TiC

Directory of Open Access Journals (Sweden)

Jürgen Geiser

2011-01-01

processes. In this paper we present a new model taken into account a self-consistent electrostatic-particle in cell model with low density Argon plasma. The collision model are based of Monte Carlo simulations is discussed for DC sputtering in lower pressure regimes. In order to simulate transport phenomena within sputtering processes realistically, a spatial and temporal knowledge of the plasma density and electrostatic field configuration is needed. Due to relatively low plasma densities, continuum fluid equations are not applicable. We propose instead a Particle-in-cell (PIC method, which allows the study of plasma behavior by computing the trajectories of finite-size particles under the action of an external and self-consistent electric field defined in a grid of points.

Boltzmann electron PIC simulation of the E-sail effect

Directory of Open Access Journals (Sweden)

P. Janhunen

2015-12-01

Full Text Available The solar wind electric sail (E-sail is a planned in-space propulsion device that uses the natural solar wind momentum flux for spacecraft propulsion with the help of long, charged, centrifugally stretched tethers. The problem of accurately predicting the E-sail thrust is still somewhat open, however, due to a possible electron population trapped by the tether. Here we develop a new type of particle-in-cell (PIC simulation for predicting E-sail thrust. In the new simulation, electrons are modelled as a fluid, hence resembling hybrid simulation, but in contrast to normal hybrid simulation, the Poisson equation is used as in normal PIC to calculate the self-consistent electrostatic field. For electron-repulsive parts of the potential, the Boltzmann relation is used. For electron-attractive parts of the potential we employ a power law which contains a parameter that can be used to control the number of trapped electrons. We perform a set of runs varying the parameter and select the one with the smallest number of trapped electrons which still behaves in a physically meaningful way in the sense of producing not more than one solar wind ion deflection shock upstream of the tether. By this prescription we obtain thrust per tether length values that are in line with earlier estimates, although somewhat smaller. We conclude that the Boltzmann PIC simulation is a new tool for simulating the E-sail thrust. This tool enables us to calculate solutions rapidly and allows to easily study different scenarios for trapped electrons.

Three-dimensional labeling program for elucidation of the geometric properties of biological particles in three-dimensional space.

Science.gov (United States)

Nomura, A; Yamazaki, Y; Tsuji, T; Kawasaki, Y; Tanaka, S

1996-09-15

For all biological particles such as cells or cellular organelles, there are three-dimensional coordinates representing the centroid or center of gravity. These coordinates and other numerical parameters such as volume, fluorescence intensity, surface area, and shape are referred to in this paper as geometric properties, which may provide critical information for the clarification of in situ mechanisms of molecular and cellular functions in living organisms. We have established a method for the elucidation of these properties, designated the three-dimensional labeling program (3DLP). Algorithms of 3DLP are so simple that this method can be carried out through the use of software combinations in image analysis on a personal computer. To evaluate 3DLP, it was applied to a 32-cell-stage sea urchin embryo, double stained with FITC for cellular protein of blastomeres and propidium iodide for nuclear DNA. A stack of optical serial section images was obtained by confocal laser scanning microscopy. The method was found effective for determining geometric properties and should prove applicable to the study of many different kinds of biological particles in three-dimensional space.

Single-particle density matrix and superfluidity in the two-dimensional Bose Coulomb fluid

International Nuclear Information System (INIS)

Minguzzi, A.; Tosi, M.P.; Davoudi, B.

2002-01-01

A study by Magro and Ceperley [Phys. Rev. Lett. 73, 826 (1994)] has shown that the ground state of the two-dimensional fluid of charged bosons with logarithmic interactions is not Bose condensed, but exhibits algebraic off-diagonal order in the single-particle density matrix ρ(r). We use a hydrodynamic Hamiltonian expressed in terms of density and phase operators, in combination with an f-sum rule on the superfluid fraction, to reproduce these results and to extend the evaluation of the density matrix to finite temperature T. This approach allows us to treat the liquid as a superfluid in the absence of a condensate. The algebraic decay of the one-body density matrix is due to correlations between phase fluctuations, and we find that the exponent in the power law is determined by the superfluid density n s (T). We also find that the plasmon gap in the single-particle energy spectrum at long wavelengths decreases with increasing T and closes at the critical temperature for the onset of superfluidity

Progress in the study of mesh refinement for particle-in-cell plasma simulations and its application to heavy ion fusion

International Nuclear Information System (INIS)

Vay, J.-L.; Friedman, A.; Grote, D.P.

2002-01-01

The numerical simulation of the driving beams in a heavy ion fusion power plant is a challenging task, and, despite rapid progress in computer power, one must consider the use of the most advanced numerical techniques. One of the difficulties of these simulations resides in the disparity of scales in time and in space which must be resolved. When these disparities are in distinctive zones of the simulation region, a method which has proven to be effective in other areas (e.g. fluid dynamics simulations) is the Adaptive-Mesh-Refinement (AMR) technique. We follow in this article the progress accomplished in the last few months in the merging of the AMR technique with Particle-In-Cell (PIC) method. This includes a detailed modeling of the Lampel-Tiefenback solution for the one-dimensional diode using novel techniques to suppress undesirable numerical oscillations and an AMR patch to follow the head of the particle distribution. We also report new results concerning the modeling of ion sources using the axisymmetric WARPRZ-AMR prototype showing the utility of an AMR patch resolving the emitter vicinity and the beam edge

Topological aspect of disclinations in two-dimensional crystals

International Nuclear Information System (INIS)

Wei-Kai, Qi; Tao, Zhu; Yong, Chen; Ji-Rong, Ren

2009-01-01

By using topological current theory, this paper studies the inner topological structure of disclinations during the melting of two-dimensional systems. From two-dimensional elasticity theory, it finds that there are topological currents for topological defects in homogeneous equation. The evolution of disclinations is studied, and the branch conditions for generating, annihilating, crossing, splitting and merging of disclinations are given. (the physics of elementary particles and fields)

The discrete cones methods for two-dimensional neutral particle transport problems with voids

International Nuclear Information System (INIS)

Watanabe, Y.; Maynard, C.W.

1983-01-01

One of the most widely applied deterministic methods for time-independent, two-dimensional neutron transport calculations is the discrete ordinates method (DSN). The DSN solution, however, fails to be accurate in a void due to the ray effect. In order to circumvent this drawback, the authors have been developing a novel approximation: the discrete cones method (DCN), where a group of particles in a cone are simultaneously traced instead of particles in discrete directions for the DSN method. Programs, which apply to the DSN method in a non-vacuum region and the DCN method in a void, have been written for transport calculations in X-Y coordinates. The solutions for test problems demonstrate mitigation of the ray effect in voids without loosing the computational efficiency of the DSN method

THE IBEX RIBBON AND THE PICKUP ION RING STABILITY IN THE OUTER HELIOSHEATH. II. MONTE-CARLO AND PARTICLE-IN-CELL MODEL RESULTS

Energy Technology Data Exchange (ETDEWEB)

Niemiec, J. [Institute of Nuclear Physics PAN, Radzikowskiego 152 31-342 Krakow (Poland); Florinski, V.; Heerikhuisen, J. [Department of Space Science and Center for Space Plasma and Aeronomic Research, University of Alabama in Huntsville, Huntsville, AL 35899 (United States); Nishikawa, K.-I., E-mail: jacek.niemiec@ifj.edu.pl, E-mail: vaf0001@uah.edu, E-mail: jh0004@uah.edu, E-mail: ken-ichi.nishikawa-1@nasa.gov [Department of Physics, University of Alabama in Huntsville, Huntsville, AL 35899 (United States)

2016-08-01

The nearly circular ribbon of energetic neutral atom (ENA) emission discovered by NASA’s Interplanetary Boundary EXplorer satellite ( IBEX ), is most commonly attributed to the effect of charge exchange of secondary pickup ions (PUIs) gyrating about the magnetic field in the outer heliosheath (OHS) and the interstellar space beyond. The first paper in the series (Paper I) presented a theoretical analysis of the pickup process in the OHS and hybrid-kinetic simulations, revealing that the kinetic properties of freshly injected proton rings depend sensitively on the details of their velocity distribution. It was demonstrated that only rings that are not too narrow (parallel thermal spread above a few km s{sup −1}) and not too wide (parallel temperature smaller than the core plasma temperature) could remain stable for a period of time long enough to generate ribbon ENAs. This paper investigates the role of electron dynamics and the extra spatial degree of freedom in the ring ion scattering process with the help of two-dimensional full particle-in-cell (PIC) kinetic simulations. A good agreement is observed between ring evolution under unstable conditions in hybrid and PIC models, and the dominant modes are found to propagate parallel to the magnetic field. We also present more realistic ribbon PUI distributions generated using Monte Carlo simulations of atomic hydrogen in the global heliosphere and examine the effect of both the cold ring-like and the hot “halo” PUIs produced from heliosheath ENAs on the ring stability. It is shown that the second PUI population enhances the fluctuation growth rate, leading to faster isotropization of the solar-wind-derived ring ions.

In flight calibrations of Ibis/PICsIT

International Nuclear Information System (INIS)

Malaguti, G.; Di Cocco, G.; Foschini, L.; Stephen, J.B.; Bazzano, A.; Ubertini, P.; Bird, A.J.; Laurent, P.; Segreto, A.

2003-01-01

PICsIT (Pixellated Imaging Caesium Iodide Telescope) is the high energy detector of the IBIS telescope on-board the INTEGRAL satellite. It consists of 4096 independent detection units, ∼ 0.7 cm 2 in cross-section, operating in the energy range between 175 keV and 10 MeV. The intrinsically low signal to noise ratio in the gamma-ray astronomy domain implies very long observations, lasting 10 5 - 10 6 s. Moreover, the image formation principle on which PICsIT works is that of coded imaging in which the entire detection plane contributes to each decoded sky pixel. For these two main reasons, the monitoring, and possible correction, of the spatial and temporal non-uniformity of pixel performances, especially in terms of gain and energy resolution, is of paramount importance. The IBIS on-board 22 Na calibration source allows the calibration of each pixel at an accuracy of <0.5% by integrating the data from a few revolutions at constant temperature. The two calibration lines, at 511 and 1275 keV, allow also the measurement and monitoring of the PICsIT energy resolution which proves to be very stable at ∼ 19% and ∼ 9% (FWHM) respectively, and consistent with the values expected analytical predictions checked against pre-launch tests. (authors)

Evaluation of the osteogenic differentiation of gingiva-derived stem cells grown on culture plates or in stem cell spheroids: Comparison of two- and three-dimensional cultures.

Science.gov (United States)

Lee, Sung-Il; Ko, Youngkyung; Park, Jun-Beom

2017-09-01

Three-dimensional cell culture systems provide a convenient in vitro model for the study of complex cell-cell and cell-matrix interactions in the absence of exogenous substrates. The current study aimed to evaluate the osteogenic differentiation potential of gingiva-derived stem cells cultured in two-dimensional or three-dimensional systems. To the best of our knowledge, the present study is the first to compare the growth of gingiva-derived stem cells in monolayer culture to a three-dimensional culture system with microwells. For three-dimensional culture, gingiva-derived stem cells were isolated and seeded into polydimethylsiloxane-based concave micromolds. Alkaline phosphatase activity and alizarin red S staining assays were then performed to evaluate osteogenesis and the degree of mineralization, respectively. Stem cell spheroids had a significantly increased level of alkaline phosphatase activity and mineralization compared with cells from the two-dimensional culture. In addition, an increase in mineralized deposits was observed with an increase in the loading cell number. The results of present study indicate that gingiva-derived stem cell spheroids exhibit an increased osteogenic potential compared with stem cells from two-dimensional culture. This highlights the potential of three-dimensional culture systems using gingiva-derived stem cells for regenerative medicine applications requiring stem cells with osteogenic potential.

Introducing a distributed unstructured mesh into gyrokinetic particle-in-cell code, XGC

Science.gov (United States)

Yoon, Eisung; Shephard, Mark; Seol, E. Seegyoung; Kalyanaraman, Kaushik

2017-10-01

XGC has shown good scalability for large leadership supercomputers. The current production version uses a copy of the entire unstructured finite element mesh on every MPI rank. Although an obvious scalability issue if the mesh sizes are to be dramatically increased, the current approach is also not optimal with respect to data locality of particles and mesh information. To address these issues we have initiated the development of a distributed mesh PIC method. This approach directly addresses the base scalability issue with respect to mesh size and, through the use of a mesh entity centric view of the particle mesh relationship, provides opportunities to address data locality needs of many core and GPU supported heterogeneous systems. The parallel mesh PIC capabilities are being built on the Parallel Unstructured Mesh Infrastructure (PUMI). The presentation will first overview the form of mesh distribution used and indicate the structures and functions used to support the mesh, the particles and their interaction. Attention will then focus on the node-level optimizations being carried out to ensure performant operation of all PIC operations on the distributed mesh. Partnership for Edge Physics Simulation (EPSI) Grant No. DE-SC0008449 and Center for Extended Magnetohydrodynamic Modeling (CEMM) Grant No. DE-SC0006618.

Two-dimensional Monte Carlo simulations of structures of a suspension comprised of magnetic and nonmagnetic particles in uniform magnetic fields

International Nuclear Information System (INIS)

Peng Xiaoling; Min Yong; Ma Tianyu; Luo Wei; Yan Mi

2009-01-01

The structures of suspensions comprised of magnetic and nonmagnetic particles in magnetic fields are studied using two-dimensional Monte Carlo simulations. The magnetic interaction among magnetic particles, magnetic field strength, and concentrations of both magnetic and nonmagnetic particles are considered as key influencing factors in the present work. The results show that chain-like clusters of magnetic particles are formed along the field direction. The size of the clusters increases with increasing magnetic interaction between magnetic particles, while it keeps nearly unchanged as the field strength increases. As the concentration of magnetic particles increases, both the number and size of the clusters increase. Moreover, nonmagnetic particles are found to hinder the migration of magnetic ones. As the concentration of nonmagnetic particles increases, the hindrance on migration of magnetic particles is enhanced

Three-dimensional single-particle tracking in live cells: news from the third dimension

International Nuclear Information System (INIS)

Dupont, A; Wehnekamp, F; Katayama, Y; Lamb, D C; Gorelashvili, M; Schüller, V; Arcizet, D; Heinrich, D

2013-01-01

Single-particle tracking (SPT) is of growing importance in the biophysical community. It is used to investigate processes such as drug and gene delivery, viral uptake, intracellular trafficking or membrane-bound protein mobility. Traditionally, SPT is performed in two dimensions (2D) because of its technical simplicity. However, life occurs in three dimensions (3D) and many methods have been recently developed to track particles in 3D. Now, is the third dimension worth the effort? Here we investigate the differences between the 2D and 3D analyses of intracellular transport with the 3D development of a time-resolved mean square displacement (MSD) analysis introduced previously. The 3D trajectories, and the 2D projections, of fluorescent nanoparticles were obtained with an orbital tracking microscope in two different cell types: in Dictyostelium discoideum ameba and in adherent, more flattened HuH-7 human cells. As expected from the different 3D organization of both cells’ cytoskeletons, a third of the active transport was lost upon projection in the ameba whereas the identification of the active phases was barely affected in the HuH-7 cells. In both cell types, we found intracellular diffusion to be anisotropic and the diffusion coefficient values derived from the 2D analysis were therefore biased. (paper)

Fourier analysis of cell-wise Block-Jacobi splitting in two-dimensional geometry

International Nuclear Information System (INIS)

Rosa, M.; Warsa, J. S.; Kelley, T. M.

2009-01-01

A Fourier analysis is conducted in two-dimensional (2D) geometry for the discrete ordinates (S N ) approximation of the neutron transport problem solved with Richardson iteration (Source Iteration) using the cell-wise Block-Jacobi (BJ) algorithm. The results of the Fourier analysis show that convergence of cell-wise BJ can degrade, leading to a spectral radius equal to 1, in problems containing optically thin cells. For problems containing cells that are optically thick, instead, the spectral radius tends to 0. Hence, in the optically thick-cell regime, cell-wise BJ is rapidly convergent even for problems that are scattering dominated, with a scattering ratio c close to 1. (authors)

Comparisons of angularly and spectrally resolved Bremsstrahlung measurements to two-dimensional multi-stage simulations of short-pulse laser-plasma interactions

Energy Technology Data Exchange (ETDEWEB)

Chen, C. D.; Kemp, A. J.; Pérez, F.; Link, A.; Key, M. H.; McLean, H.; Ping, Y.; Patel, P. K. [Lawrence Livermore National Laboratory (United States); Beg, F. N.; Chawla, S.; Sorokovikova, A.; Westover, B. [University of California, San Diego (United States); Morace, A. [University of Milan (Italy); Stephens, R. B. [General Atomics (United States); Streeter, M. [Imperial College London (United Kingdom)

2013-05-15

A 2-D multi-stage simulation model incorporating realistic laser conditions and a fully resolved electron distribution handoff has been developed and compared to angularly and spectrally resolved Bremsstrahlung measurements from high-Z planar targets. For near-normal incidence and 0.5-1 × 10{sup 20} W/cm{sup 2} intensity, particle-in-cell (PIC) simulations predict the existence of a high energy electron component consistently directed away from the laser axis, in contrast with previous expectations for oblique irradiation. Measurements of the angular distribution are consistent with a high energy component when directed along the PIC predicted direction, as opposed to between the target normal and laser axis as previously measured.

Approximate solutions for the two-dimensional integral transport equation. Solution of complex two-dimensional transport problems

International Nuclear Information System (INIS)

Sanchez, Richard.

1980-11-01

This work is divided into two parts: the first part deals with the solution of complex two-dimensional transport problems, the second one (note CEA-N-2166) treats the critically mixed methods of resolution. A set of approximate solutions for the isotropic two-dimensional neutron transport problem has been developed using the interface current formalism. The method has been applied to regular lattices of rectangular cells containing a fuel pin, cladding, and water, or homogenized structural material. The cells are divided into zones that are homogeneous. A zone-wise flux expansion is used to formulate a direct collision probability problem within a cell. The coupling of the cells is effected by making extra assumptions on the currents entering and leaving the interfaces. Two codes have been written: CALLIOPE uses a cylindrical cell model and one or three terms for the flux expansion, and NAUSICAA uses a two-dimensional flux representation and does a truly two-dimensional calculation inside each cell. In both codes, one or three terms can be used to make a space-independent expansion of the angular fluxes entering and leaving each side of the cell. The accuracies and computing times achieved with the different approximations are illustrated by numerical studies on two benchmark problems and by calculations performed in the APOLLO multigroup code [fr

One-dimensional integral equations for a system of three identical particles in the boundary condition models and the possibility of changing the off-shell behaviour of the two-particle t-matrix

International Nuclear Information System (INIS)

Efimov, V.N.; Schulz, H.

1976-01-01

It is shown that in the framework of the boundary condition models (BCM) for the two-particle interaction the Schroedinger equation for the system of three identical bosons can be reduced to the one-dimensional integral equation in an exact way. The method used for obtaining such an equation is based on a special consideration of the two-particle off-shell wave functions. The binding energy of the simple three-particle system is calculated. It is indicated that by means of the equation obtained it is possible to change the off-shell behaviour of the two-particle t-matrix and therefore to simulate three particle effects. (Auth.)

Monte Carlo simulation of ion-neutral charge exchange collisions and grid erosion in an ion thruster

Science.gov (United States)

Peng, Xiaohang; Ruyten, Wilhelmus M.; Keefer, Dennis

1991-01-01

A combined particle-in-cell (PIC)/Monte Carlo simulation model has been developed in which the PIC method is used to simulate the charge exchange collisions. It is noted that a number of features were reproduced correctly by this code, but that its assumption of two-dimensional axisymmetry for a single set of grid apertures precluded the reproduction of the most characteristic feature of actual test data; namely, the concentrated grid erosion at the geometric center of the hexagonal aperture array. The first results of a three-dimensional code, which takes into account the hexagonal symmetry of the grid, are presented. It is shown that, with this code, the experimentally observed erosion patterns are reproduced correctly, demonstrating explicitly the concentration of sputtering between apertures.

Realistic PIC modelling of laser-plasma interaction: a direct implicit method with adjustable damping and high order weight functions

International Nuclear Information System (INIS)

Drouin, M.

2009-11-01

This research thesis proposes a new formulation of the relativistic implicit direct method, based on the weak formulation of the wave equation which is solved by means of a Newton algorithm. The first part of this thesis deals with the properties of the explicit particle-in-cell (PIC) methods: properties and limitations of an explicit PIC code, linear analysis of a numerical plasma, numerical heating phenomenon, interest of a higher order interpolation function, and presentation of two applications in high density relativistic laser-plasma interaction. The second and main part of this report deals with adapting the direct implicit method to laser-plasma interaction: presentation of the state of the art, formulating of the direct implicit method, resolution of the wave equation. The third part concerns various numerical and physical validations of the ELIXIRS code: case of laser wave propagation in vacuum, demonstration of the adjustable damping which is a characteristic of the proposed algorithm, influence of space-time discretization on energy conservation, expansion of a thermal plasma in vacuum, two cases of plasma-beam unsteadiness in relativistic regime, and then a case of the overcritical laser-plasma interaction

Quasi-integrability and two-dimensional QCD

International Nuclear Information System (INIS)

Abdalla, E.; Mohayaee, R.

1996-10-01

The notion of integrability in two-dimensional QCD is discussed. We show that in spite of an infinite number of conserved charges, particle production is not entirely suppressed. This phenomenon, which we call quasi-integrability, is explained in terms of quantum corrections to the combined algebra of higher-conserved and spectrum-generating currents. We predict the qualitative form of particle production probabilities and verify that they are in agreement with numerical data. We also discuss four-dimensional self-dual Yang-Mills theory in the light of our results. (author). 25 refs, 4 figs, 1 tab

Numerical study on the two-dimensional flows of plasma and ionizing gas using trial particles

International Nuclear Information System (INIS)

Brushlinskij, K.V.; Kozlov, A.N.; Morozov, A.I.

1985-01-01

Two-dimensional flows of plasma and ionized αs in a channel between two coaxial electrodes are considered in the MHD-model with account of Hall effect. Stationary solutions of the problem on the flow are obtained either analytically in approximation of a ''smooth'' channel - for ideal conducting plasma, or numerically using the methos of establishment - in the ge-neral case of finite conductivity. A method of further numerical analysis of some peculiarities of flow is suggested in the paper. It is based on studying dynamics of single ''test'' particles in fields of the main MHD plasma flow. Trajectory of the test ion is calculated with account for interaction forces with earlier determined electromagentic field and friction responsible for Coulomb collisions with particles of the background flow. The calculations display trajectories of test particles with different masses, initial positions and initial rates. They are shown to be dose to current lines of background medium in plasma of finite conductivity, that testified to the virtue of effectiveness of the MHD-model. In case of ideal conductivity trajectories of test and background particles can noticeably differ from one another. Stabilization effects of motion of particles accidentally knocked out from the flow and separation of pariticles of different mass by electromao.netic forces are considered

Post-arc current simulation based on measurement in vacuum circuit breaker with a one-dimensional particle-in-cell model

Science.gov (United States)

Jia, Shenli; Mo, Yongpeng; Shi, Zongqian; Li, Junliang; Wang, Lijun

2017-10-01

The post-arc dielectric recovery process has a decisive effect on the current interruption performance in a vacuum circuit breaker. The dissipation of residual plasma at the moment of current zero under the transient recovery voltage, which is the first stage of the post-arc dielectric recovery process and forms the post-arc current, has attracted many concerns. A one-dimensional particle-in-cell model is developed to simulate the measured post-arc current in the vacuum circuit breaker in this paper. At first, the parameters of the residual plasma are estimated roughly by the waveform of the post-arc current which is taken from measurements. After that, different components of the post-arc current, which are formed by the movement of charged particles in the residual plasma, are discussed. Then, the residual plasma density is adjusted according to the proportion of electrons and ions absorbed by the post-arc anode derived from the particle-in-cell simulation. After this adjustment, the post-arc current waveform obtained from the simulation is closer to that obtained from measurements.

Particle-in-cell simulation of two-dimensional electron velocity shear driven instability in relativistic domain

Energy Technology Data Exchange (ETDEWEB)

Shukla, Chandrasekhar, E-mail: chandrasekhar.shukla@gmail.com; Das, Amita, E-mail: amita@ipr.res.in [Institute for Plasma Research, Bhat, Gandhinagar 382428 (India); Patel, Kartik [Bhabha Atomic Research Centre, Trombay, Mumbai 400 085 (India)

2016-08-15

We carry out particle-in-cell simulations to study the instabilities associated with a 2-D sheared electron flow configuration against a neutralizing background of ions. Both weak and strong relativistic flow velocities are considered. In the weakly relativistic case, we observe the development of electromagnetic Kelvin-Helmholtz instability with similar characteristics as that predicted by the electron Magnetohydrodynamic (EMHD) model. On the contrary, in a strong relativistic case, the compressibility effects of electron fluid dominate and introduce upper hybrid electrostatic oscillations transverse to the flow which are very distinct from EMHD fluid behavior. In the nonlinear regime, both weak and strong relativistic cases lead to turbulence with broad power law spectrum.

Fluctuations and symmetries in two-dimensional active gels.

Science.gov (United States)

Sarkar, N; Basu, A

2011-04-01

Motivated by the unique physical properties of biological active matter, e.g., cytoskeletal dynamics in eukaryotic cells, we set up effective two-dimensional (2d) coarse-grained hydrodynamic equations for the dynamics of thin active gels with polar or nematic symmetries. We use the well-known three-dimensional (3d) descriptions (K. Kruse et al., Eur. Phys. J. E 16, 5 (2005); A. Basu et al., Eur. Phys. J. E 27, 149 (2008)) for thin active-gel samples confined between parallel plates with appropriate boundary conditions to derive the effective 2d constitutive relations between appropriate thermodynamic fluxes and generalised forces for small deviations from equilibrium. We consider three distinct cases, characterised by spatial symmetries and boundary conditions, and show how such considerations dictate the structure of the constitutive relations. We use these to study the linear instabilities, calculate the correlation functions and the diffusion constant of a small tagged particle, and elucidate their dependences on the activity or nonequilibrium drive.

A fully kinetic, self-consistent particle simulation model of the collisionless plasma--sheath region

International Nuclear Information System (INIS)

Procassini, R.J.; Birdsall, C.K.; Morse, E.C.

1990-01-01

A fully kinetic particle-in-cell (PIC) model is used to self-consistently determine the steady-state potential profile in a collisionless plasma that contacts a floating, absorbing boundary. To balance the flow of particles to the wall, a distributed source region is used to inject particles into the one-dimensional system. The effect of the particle source distribution function on the source region and collector sheath potential drops, and particle velocity distributions is investigated. The ion source functions proposed by Emmert et al. [Phys. Fluids 23, 803 (1980)] and Bissell and Johnson [Phys. Fluids 30, 779 (1987)] (and various combinations of these) are used for the injection of both ions and electrons. The values of the potential drops obtained from the PIC simulations are compared to those from the theories of Emmert et al., Bissell and Johnson, and Scheuer and Emmert [Phys. Fluids 31, 3645 (1988)], all of which assume that the electron density is related to the plasma potential via the Boltzmann relation. The values of the source region and total potential drop are found to depend on the choice of the electron source function, as well as the ion source function. The question of an infinite electric field at the plasma--sheath interface, which arises in the analyses of Bissell and Johnson and Scheuer and Emmert, is also addressed

Diffusion and sorption in particles and two-dimensional dispersion in a porous media

International Nuclear Information System (INIS)

Rasmuson, A.

1980-01-01

A solution of the two-dimensional differential equation of dispersion from a disk source, coupled with a differential equation of diffusion and sorption in particles, is developed. The solution is obtained by the successive use of the Laplace and the Hankel transforms and is given in the form of an infinite double-integral. If the lateral dispersion is negligible, the solution is shown to simplify to a solution presented earlier. Dimensionless quantities are introduced. A steady-state condition is obtained after long time. This is investigated in some detail. An expression is derived for the highest concentration which may be expected at a point in space. An important relation is obtained when longitudinal dispersion is neglected. The solution for any value of the lateral dispersion coefficient and radial distance from the source is then obtained by simple multiplication of a solution for no lateral dispersion with the steady-state value. A method for integrating the infinite double integral is given. Some typical examples are shown. (Auth.)

Two-dimensional electron flow in pulsed power transmission lines and plasma opening switches

International Nuclear Information System (INIS)

Church, B.W.; Longcope, D.W.; Ng, C.K.; Sudan, R.N.

1991-01-01

The operation of magnetically insulated transmission lines (MITL) and the interruption of current in a plasma opening switch (POS) are determined by the physics of the electrons emitted by the cathode surface. A mathematical model describes the self-consistent two-dimensional flow of an electron fluid. A finite element code, FERUS, has been developed to solve the two equations which describe Poisson's and Ampere's law in two dimensions. The solutions from this code are obtained for parameters where the electron orbits are considerably modified by the self-magnetic field of the current. Next, the self-insulated electron flow in a MITL with a step change in cross-section is studied using a conventional two-dimensional fully electromagnetic particle-in-cell code, MASK. The equations governing two-dimensional quasi-static electron flow are solved numerically by a third technique which is suitable for predicting current interruption in a POS. The object of the study is to determine the critical load impedance, Z CL , required for current interruption for a given applied voltage, cathode voltage and plasma length. (author). 9 refs, 5 figs

Metal Detector By Using PIC Microcontroller Interfacing With PC

Directory of Open Access Journals (Sweden)

Yin Min Theint

2015-06-01

Full Text Available Abstract This system proposes metal detector by using PIC microcontroller interfacing with PC. The system uses PIC microcontroller as the main controller whether the detected metal is ferrous metal or non-ferrous metal. Among various types of metal sensors and various types of metal detecting technologies concentric type induction coil sensor and VLF very low frequency metal detecting technology are used in this system. This system consists of two configurations Hardware configuration and Software configuration. The hardware components include induction coil sensors which senses the frequency changes of metal a PIC microcontroller personal computer PC buzzer light emitting diode LED and webcam. The software configuration includes a program controller interface. PIC MikroCprogramming language is used to implement the control system. This control system is based on the PIC 16F887 microcontroller.This system is mainly used in mining and high security places such as airport plaza shopping mall and governmental buildings.

One-particle versus two-particle crossover in weakly coupled Hubbard chains and ladders: perturbative renormalization group approach

International Nuclear Information System (INIS)

Kishine, Jun-Ichiro; Yonemitsu, Kenji

1998-01-01

Physical nature of dimensional crossovers in weakly coupled Hubbard chains and ladders has been discussed within the framework of the perturbative renormalization-group (PRG) approach. The difference between these two cases originates from different universality classes which the corresponding isolated systems belong to. In the present work, we discuss the nature of the dimensional crossovers in the weakly coupled chains and ladders, with emphasis on the difference between the two cases within the framework of the PRG approach. The difference of the universality class of the isolated chain and ladder profoundly affects the relevance or irrelevance of the inter-chain/ladder one-particle hopping. The strong coupling phase of the isolated ladder makes the one-particle process irrelevant so that the d-wave superconducting transition can be induced via the two-particle crossover in the weakly coupled ladders. The weak coupling phase of the isolated chain makes the one-particle process relevant so that the two-particle crossover can hardly be realized in the coupled chains. (Copyright (1998) World Scientific Publishing Co. Pte. Ltd)

Dynamic Load Balancing for PIC code using Eulerian/Lagrangian partitioning

OpenAIRE

Sauget, Marc; Latu, Guillaume

2017-01-01

This document presents an analysis of different load balance strategies for a Plasma physics code that models high energy particle beams with PIC method. A comparison of different load balancing algorithms is given: static or dynamic ones. Lagrangian and Eulerian partitioning techniques have been investigated.

Performance analysis of three-dimensional-triple-level cell and two-dimensional-multi-level cell NAND flash hybrid solid-state drives

Science.gov (United States)

Sakaki, Yukiya; Yamada, Tomoaki; Matsui, Chihiro; Yamaga, Yusuke; Takeuchi, Ken

2018-04-01

In order to improve performance of solid-state drives (SSDs), hybrid SSDs have been proposed. Hybrid SSDs consist of more than two types of NAND flash memories or NAND flash memories and storage-class memories (SCMs). However, the cost of hybrid SSDs adopting SCMs is more expensive than that of NAND flash only SSDs because of the high bit cost of SCMs. This paper proposes unique hybrid SSDs with two-dimensional (2D) horizontal multi-level cell (MLC)/three-dimensional (3D) vertical triple-level cell (TLC) NAND flash memories to achieve higher cost-performance. The 2D-MLC/3D-TLC hybrid SSD achieves up to 31% higher performance than the conventional 2D-MLC/2D-TLC hybrid SSD. The factors of different performance between the proposed hybrid SSD and the conventional hybrid SSD are analyzed by changing its block size, read/write/erase latencies, and write unit of 3D-TLC NAND flash memory, by means of a transaction-level modeling simulator.

Few helium atoms in quasi two-dimensional space

International Nuclear Information System (INIS)

Kilic, Srecko; Vranjes, Leandra

2003-01-01

Two, three and four 3 He and 4 He atoms in quasi two-dimensional space above graphite and cesium surfaces and in 'harmonic' potential perpendicular to the surface have been studied. Using some previously examined variational wave functions and the Diffusion Monte Carlo procedure, it has been shown that all molecules: dimers, trimers and tetramers, are bound more strongly than in pure two- and three-dimensional space. The enhancement of binding with respect to unrestricted space is more pronounced on cesium than on graphite. Furthermore, for 3 He 3 ( 3 He 4 ) on all studied surfaces, there is an indication that the configuration of a dimer and a 'free' particle (two dimers) may be equivalently established

SD card projects using the PIC microcontroller

CERN Document Server

Ibrahim, Dogan

2010-01-01

PIC Microcontrollers are a favorite in industry and with hobbyists. These microcontrollers are versatile, simple, and low cost making them perfect for many different applications. The 8-bit PIC is widely used in consumer electronic goods, office automation, and personal projects. Author, Dogan Ibrahim, author of several PIC books has now written a book using the PIC18 family of microcontrollers to create projects with SD cards. This book is ideal for those practicing engineers, advanced students, and PIC enthusiasts that want to incorporate SD Cards into their devices. SD cards are che

The Darwin direct implicit particle-in-cell (DADIPIC) method for simulation of low frequency plasma phenomena

International Nuclear Information System (INIS)

Gibbons, M.R.; Hewett, D.W.

1995-01-01

We describe a new algorithm for simulating low frequency, kinetic phenomena in plasma. Darwin direct implicit particle-in-cell (DADIPIC), as its name implies, is a combination of the Darwin and direct implicit methods. 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. With the elimination of these constraints spatial and temporal discretization can be much larger than that possible with explicit, electrodynamic PIC. We discuss the algorithms for pushing the particles and solving the fields in 2D cartesian geometry. We also detail boundary conditions for conductors and dielectrics. Finally, we present two test cases, electron cyclotron waves and collisionless heating in inductively coupled plasmas. For these test cases DADIPIC shows agreement with analytic kinetic theory and good energy conservation characteristics. 33 refs., 7 figs., 2 tabs

Efficient particle-in-cell simulation of auroral plasma phenomena using a CUDA enabled graphics processing unit

Science.gov (United States)

Sewell, Stephen

This thesis introduces a software framework that effectively utilizes low-cost commercially available Graphic Processing Units (GPUs) to simulate complex scientific plasma phenomena that are modeled using the Particle-In-Cell (PIC) paradigm. The software framework that was developed conforms to the Compute Unified Device Architecture (CUDA), a standard for general purpose graphic processing that was introduced by NVIDIA Corporation. This framework has been verified for correctness and applied to advance the state of understanding of the electromagnetic aspects of the development of the Aurora Borealis and Aurora Australis. For each phase of the PIC methodology, this research has identified one or more methods to exploit the problem's natural parallelism and effectively map it for execution on the graphic processing unit and its host processor. The sources of overhead that can reduce the effectiveness of parallelization for each of these methods have also been identified. One of the novel aspects of this research was the utilization of particle sorting during the grid interpolation phase. The final representation resulted in simulations that executed about 38 times faster than simulations that were run on a single-core general-purpose processing system. The scalability of this framework to larger problem sizes and future generation systems has also been investigated.

Three-dimensional particle-in-cell simulation on gain saturation effect of microchannel plate

Energy Technology Data Exchange (ETDEWEB)

Wang, Qiangqiang; Yuan, Zheng; Cao, Zhurong, E-mail: cao33jin@aliyun.com; Deng, Bo; Chen, Tao; Deng, Keli [Research Center of Laser Fusion, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900 (China)

2016-07-15

We present here the results of the simulation work, using the three-dimensional particle-in-cell method, on the performance of the lead glass microchannel plate under saturated state. We calculated the electron cascade process with different DC bias voltages under both self-consistent condition and non-self-consistent condition. The comparative results have demonstrated that the strong self-consistent field can suppress the cascade process and make the microchannel plate saturated. The simulation results were also compared to the experimental data and good agreement was obtained. The simulation results also show that the electron multiplication process in the channel is accompanied by the buildup process of positive charges in the channel wall. Though the interactions among the secondary electron cloud in the channel, the positive charges in the channel wall, and the external acceleration field can make the electron-surface collision more frequent, the collision energy will be inevitably reduced, thus the electron gain will also be reduced.

Plasmoid statistics in relativistic magnetic reconnection

Science.gov (United States)

Petropoulou, M.; Christie, I. M.; Sironi, L.; Giannios, D.

2018-04-01

Plasmoids, overdense blobs of plasma containing magnetic fields and high-energy particles, are a self-consistent outcome of the reconnection process in the relativistic regime. Recent two-dimensional particle-in-cell (PIC) simulations have shown that plasmoids can undergo a variety of processes (e.g. mergers, bulk acceleration, growth, and advection) within the reconnection layer. We developed a Monte Carlo code, benchmarked with the recent PIC simulations, to examine the effects of these processes on the steady-state size and momentum distributions of the plasmoid chain. The differential plasmoid size distribution is shown to be a power law, ranging from a few plasma skin depths to ˜0.1 of the reconnection layer's length. The power-law slope is shown to be linearly dependent upon the ratio of the plasmoid acceleration and growth rates, which slightly decreases with increasing plasma magnetization. We perform a detailed comparison of our results with those of recent PIC simulations and briefly discuss the astrophysical implications of our findings through the representative case of flaring events from blazar jets.

Electron and ion heating by whistler turbulence: Three-dimensional particle-in-cell simulations

International Nuclear Information System (INIS)

Hughes, R. Scott; Gary, S. Peter; Wang, Joseph

2014-01-01

Three-dimensional particle-in-cell simulations of decaying whistler turbulence are carried out on a collisionless, homogeneous, magnetized, electron-ion plasma model. In addition, the simulations use an initial ensemble of relatively long wavelength whistler modes with a broad range of initial propagation directions with an initial electron beta β e = 0.05. The computations follow the temporal evolution of the fluctuations as they cascade into broadband turbulent spectra at shorter wavelengths. Three simulations correspond to successively larger simulation boxes and successively longer wavelengths of the initial fluctuations. The computations confirm previous results showing electron heating is preferentially parallel to the background magnetic field B o , and ion heating is preferentially perpendicular to B o . The new results here are that larger simulation boxes and longer initial whistler wavelengths yield weaker overall dissipation, consistent with linear dispersion theory predictions of decreased damping, stronger ion heating, consistent with a stronger ion Landau resonance, and weaker electron heating

Pair Interaction of Dislocations in Two-Dimensional Crystals

Science.gov (United States)

Eisenmann, C.; Gasser, U.; Keim, P.; Maret, G.; von Grünberg, H. H.

2005-10-01

The pair interaction between crystal dislocations is systematically explored by analyzing particle trajectories of two-dimensional colloidal crystals measured by video microscopy. The resulting pair energies are compared to Monte Carlo data and to predictions derived from the standard Hamiltonian of the elastic theory of dislocations. Good agreement is found with respect to the distance and temperature dependence of the interaction potential, but not regarding the angle dependence where discrete lattice effects become important. Our results on the whole confirm that the dislocation Hamiltonian allows a quantitative understanding of the formation and interaction energies of dislocations in two-dimensional crystals.

Storage of Maize in Purdue Improved Crop Storage (PICS) Bags.

Science.gov (United States)

Williams, Scott B; Murdock, Larry L; Baributsa, Dieudonne

2017-01-01

Interest in using hermetic technologies as a pest management solution for stored grain has risen in recent years. One hermetic approach, Purdue Improved Crop Storage (PICS) bags, has proven successful in controlling the postharvest pests of cowpea. This success encouraged farmers to use of PICS bags for storing other crops including maize. To assess whether maize can be safely stored in PICS bags without loss of quality, we carried out laboratory studies of maize grain infested with Sitophilus zeamais (Motshulsky) and stored in PICS triple bags or in woven polypropylene bags. Over an eight month observation period, temperatures in the bags correlated with ambient temperature for all treatments. Relative humidity inside PICS bags remained constant over this period despite the large changes that occurred in the surrounding environment. Relative humidity in the woven bags followed ambient humidity closely. PICS bags containing S. zeamais-infested grain saw a significant decline in oxygen compared to the other treatments. Grain moisture content declined in woven bags, but remained high in PICS bags. Seed germination was not significantly affected over the first six months in all treatments, but declined after eight months of storage when infested grain was held in woven bags. Relative damage was low across treatments and not significantly different between treatments. Overall, maize showed no signs of deterioration in PICS bags versus the woven bags and PICS bags were superior to woven bags in terms of specific metrics of grain quality.

Few quantum particles on one dimensional lattices

Energy Technology Data Exchange (ETDEWEB)

Valiente Cifuentes, Manuel

2010-06-18

There is currently a great interest in the physics of degenerate quantum gases and low-energy few-body scattering due to the recent experimental advances in manipulation of ultracold atoms by light. In particular, almost perfect periodic potentials, called optical lattices, can be generated. The lattice spacing is fixed by the wavelength of the laser field employed and the angle betwen the pair of laser beams; the lattice depth, defining the magnitude of the different band gaps, is tunable within a large interval of values. This flexibility permits the exploration of different regimes, ranging from the ''free-electron'' picture, modified by the effective mass for shallow optical lattices, to the tight-binding regime of a very deep periodic potential. In the latter case, effective single-band theories, widely used in condensed matter physics, can be implemented with unprecedent accuracy. The tunability of the lattice depth is nowadays complemented by the use of magnetic Feshbach resonances which, at very low temperatures, can vary the relevant atom-atom scattering properties at will. Moreover, optical lattices loaded with gases of effectively reduced dimensionality are experimentally accessible. This is especially important for one spatial dimension, since most of the exactly solvable models in many-body quantum mechanics deal with particles on a line; therefore, experiments with one-dimensional gases serve as a testing ground for many old and new theories which were regarded as purely academic not so long ago. The physics of few quantum particles on a one-dimensional lattice is the topic of this thesis. Most of the results are obtained in the tight-binding approximation, which is amenable to exact numerical or analytical treatment. For the two-body problem, theoretical methods for calculating the stationary scattering and bound states are developed. These are used to obtain, in closed form, the two-particle solutions of both the Hubbard and

Few quantum particles on one dimensional lattices

International Nuclear Information System (INIS)

Valiente Cifuentes, Manuel

2010-01-01

There is currently a great interest in the physics of degenerate quantum gases and low-energy few-body scattering due to the recent experimental advances in manipulation of ultracold atoms by light. In particular, almost perfect periodic potentials, called optical lattices, can be generated. The lattice spacing is fixed by the wavelength of the laser field employed and the angle betwen the pair of laser beams; the lattice depth, defining the magnitude of the different band gaps, is tunable within a large interval of values. This flexibility permits the exploration of different regimes, ranging from the ''free-electron'' picture, modified by the effective mass for shallow optical lattices, to the tight-binding regime of a very deep periodic potential. In the latter case, effective single-band theories, widely used in condensed matter physics, can be implemented with unprecedent accuracy. The tunability of the lattice depth is nowadays complemented by the use of magnetic Feshbach resonances which, at very low temperatures, can vary the relevant atom-atom scattering properties at will. Moreover, optical lattices loaded with gases of effectively reduced dimensionality are experimentally accessible. This is especially important for one spatial dimension, since most of the exactly solvable models in many-body quantum mechanics deal with particles on a line; therefore, experiments with one-dimensional gases serve as a testing ground for many old and new theories which were regarded as purely academic not so long ago. The physics of few quantum particles on a one-dimensional lattice is the topic of this thesis. Most of the results are obtained in the tight-binding approximation, which is amenable to exact numerical or analytical treatment. For the two-body problem, theoretical methods for calculating the stationary scattering and bound states are developed. These are used to obtain, in closed form, the two-particle solutions of both the Hubbard and extended Hubbard models

SMILEI: A collaborative, open-source, multi-purpose PIC code for the next generation of super-computers

Science.gov (United States)

Grech, Mickael; Derouillat, J.; Beck, A.; Chiaramello, M.; Grassi, A.; Niel, F.; Perez, F.; Vinci, T.; Fle, M.; Aunai, N.; Dargent, J.; Plotnikov, I.; Bouchard, G.; Savoini, P.; Riconda, C.

2016-10-01

Over the last decades, Particle-In-Cell (PIC) codes have been central tools for plasma simulations. Today, new trends in High-Performance Computing (HPC) are emerging, dramatically changing HPC-relevant software design and putting some - if not most - legacy codes far beyond the level of performance expected on the new and future massively-parallel super computers. SMILEI is a new open-source PIC code co-developed by both plasma physicists and HPC specialists, and applied to a wide range of physics-related studies: from laser-plasma interaction to astrophysical plasmas. It benefits from an innovative parallelization strategy that relies on a super-domain-decomposition allowing for enhanced cache-use and efficient dynamic load balancing. Beyond these HPC-related developments, SMILEI also benefits from additional physics modules allowing to deal with binary collisions, field and collisional ionization and radiation back-reaction. This poster presents the SMILEI project, its HPC capabilities and illustrates some of the physics problems tackled with SMILEI.

One- and Two-dimensional Solitary Wave States in the Nonlinear Kramers Equation with Movement Direction as a Variable

Science.gov (United States)

Sakaguchi, Hidetsugu; Ishibashi, Kazuya

2018-06-01

We study self-propelled particles by direct numerical simulation of the nonlinear Kramers equation for self-propelled particles. In our previous paper, we studied self-propelled particles with velocity variables in one dimension. In this paper, we consider another model in which each particle exhibits directional motion. The movement direction is expressed with a variable ϕ. We show that one-dimensional solitary wave states appear in direct numerical simulations of the nonlinear Kramers equation in one- and two-dimensional systems, which is a generalization of our previous result. Furthermore, we find two-dimensionally localized states in the case that each self-propelled particle exhibits rotational motion. The center of mass of the two-dimensionally localized state exhibits circular motion, which implies collective rotating motion. Finally, we consider a simple one-dimensional model equation to qualitatively understand the formation of the solitary wave state.

Two-dimensional ion effects in relativistic diodes

International Nuclear Information System (INIS)

Poukey, J.W.

1975-01-01

In relativistic diodes, ions are emitted from the anode plasma. The effects and properties of these ions are studied via a two-dimensional particle simulation code. The space charge of these ions enhances the electron emission, and this additional current (including that of the ions, themselves) aids in obtaining superpinched electron beams for use in pellet fusion studies. (U.S.)

Particle-in-cell simulation of helical structure onset in plasma fiber with dust grains

International Nuclear Information System (INIS)

Kulhanek, Petr; Bren, David; Kaizr, Vaclav; Pasek, Jan

2002-01-01

Fully three dimensional PIC program package for the helical pinch numerical simulation was developed in our department. Both electromagnetic and gravitational interactions are incorporated into the model. Collisions are treated via Monte Carlo methods. The program package enabled to prove the conditions of onset of spiral and helical structures in the pinch

Two dimensional neutral transport analysis in tokamak plasma

International Nuclear Information System (INIS)

Shimizu, Katsuhiro; Azumi, Masafumi

1987-02-01

Neutral particle influences the particle and energy balance, and play an important role on sputtering impurity and the charge exchange loss of neutral beam injection. In order to study neutral particle behaviour including the effects of asymmetric source and divertor configuration, the two dimensional neutral transport code has been developed using the Monte-Carlo techniques. This code includes the calculation of the H α radiation intensity based on the collisional-radiation model. The particle confinement time of the joule heated plasma in JT-60 tokamak is evaluated by comparing the calculated H α radiation intensity with the experimental data. The effect of the equilibrium on the neutral density profile in high-β plasma is also investigated. (author)

Fractional calculus phenomenology in two-dimensional plasma models

Science.gov (United States)

Gustafson, Kyle; Del Castillo Negrete, Diego; Dorland, Bill

2006-10-01

Transport processes in confined plasmas for fusion experiments, such as ITER, are not well-understood at the basic level of fully nonlinear, three-dimensional kinetic physics. Turbulent transport is invoked to describe the observed levels in tokamaks, which are orders of magnitude greater than the theoretical predictions. Recent results show the ability of a non-diffusive transport model to describe numerical observations of turbulent transport. For example, resistive MHD modeling of tracer particle transport in pressure-gradient driven turbulence for a three-dimensional plasma reveals that the superdiffusive (2̂˜t^α where α> 1) radial transport in this system is described quantitatively by a fractional diffusion equation Fractional calculus is a generalization involving integro-differential operators, which naturally describe non-local behaviors. Our previous work showed the quantitative agreement of special fractional diffusion equation solutions with numerical tracer particle flows in time-dependent linearized dynamics of the Hasegawa-Mima equation (for poloidal transport in a two-dimensional cold-ion plasma). In pursuit of a fractional diffusion model for transport in a gyrokinetic plasma, we now present numerical results from tracer particle transport in the nonlinear Hasegawa-Mima equation and a planar gyrokinetic model. Finite Larmor radius effects will be discussed. D. del Castillo Negrete, et al, Phys. Rev. Lett. 94, 065003 (2005).

Relativistic effects in the energy loss of a fast charged particle moving parallel to a two-dimensional electron gas

Science.gov (United States)

Mišković, Zoran L.; Akbari, Kamran; Segui, Silvina; Gervasoni, Juana L.; Arista, Néstor R.

2018-05-01

We present a fully relativistic formulation for the energy loss rate of a charged particle moving parallel to a sheet containing two-dimensional electron gas, allowing that its in-plane polarization may be described by different longitudinal and transverse conductivities. We apply our formulation to the case of a doped graphene layer in the terahertz range of frequencies, where excitation of the Dirac plasmon polariton (DPP) in graphene plays a major role. By using the Drude model with zero damping we evaluate the energy loss rate due to excitation of the DPP, and show that the retardation effects are important when the incident particle speed and its distance from graphene both increase. Interestingly, the retarded energy loss rate obtained in this manner may be both larger and smaller than its non-retarded counterpart for different combinations of the particle speed and distance.

Young gamma-ray pulsar: from modeling the gamma-ray emission to the particle-in-cell simulations of the global magnetosphere

Science.gov (United States)

Brambilla, Gabriele; Kalapotharakos, Constantions; Timokhin, Andrey; Kust Harding, Alice; Kazanas, Demosthenes

2016-04-01

Accelerated charged particles flowing in the magnetosphere produce pulsar gamma-ray emission. Pair creation processes produce an electron-positron plasma that populates the magnetosphere, in which the plasma is very close to force-free. However, it is unknown how and where the plasma departs from the ideal force-free condition, which consequently inhibits the understanding of the emission generation. We found that a dissipative magnetosphere outside the light cylinder effectively reproduces many aspects of the young gamma-ray pulsar emission as seen by the Fermi Gamma-ray Space Telescope, and through particle-in-cell simulations (PIC), we started explaining this configuration self-consistently. These findings show that, together, a magnetic field structure close to force-free and the assumption of gamma-ray curvature radiation as the emission mechanism are strongly compatible with the observations. Two main issues from the previously used models that our work addresses are the inability to explain luminosity, spectra, and light curve features at the same time and the inconsistency of the electrodynamics. Moreover, using the PIC simulations, we explore the effects of different pair multiplicities on the magnetosphere configurations and the locations of the accelerating regions. Our work aims for a self-consistent modeling of the magnetosphere, connecting the microphysics of the pair-plasma to the global magnetosphere macroscopic quantities. This direction will lead to a greater understanding of pulsar emission at all wavelengths, as well as to concrete insights into the physics of the magnetosphere.

Review of multi-dimensional large-scale kinetic simulation and physics validation of ion acceleration in relativistic laser-matter interaction

International Nuclear Information System (INIS)

Wu, Hui-Chun; Hegelich, B.M.; Fernandez, J.C.; Shah, R.C.; Palaniyappan, S.; Jung, D.; Yin, L.; Albright, B.J.; Bowers, K.; Kwan, T.J.

2012-01-01

Two new experimental technologies enabled realization of Break-out afterburner (BOA) - High quality Trident laser and free-standing C nm-targets. VPIC is an powerful tool for fundamental research of relativistic laser-matter interaction. Predictions from VPIC are validated - Novel BOA and Solitary ion acceleration mechanisms. VPIC is a fully explicit Particle In Cell (PIC) code: models plasma as billions of macro-particles moving on a computational mesh. VPIC particle advance (which typically dominates computation) has been optimized extensively for many different supercomputers. Laser-driven ions lead to realization promising applications - Ion-based fast ignition; active interrogation, hadron therapy.

PICs in the injector complex - what are we talking about?

International Nuclear Information System (INIS)

Hanke, K

2014-01-01

This presentation will identify PIC activities for the LHC injector chain, and point out borderline cases to pure consolidation and upgrade. The most important PIC items will be listed for each LIU project (PSB, PS, SPS) and categorized by a) the risk if not performed and b) the implications of doing them. This will in particular address the consequences on performance, schedule, reliability, commissioning time, operational complexity etc. The additional cost of PICs with regard to pure consolidation will be estimated and possible time lines for the implementation of the PICs will be discussed. In this context, it will be evaluated if the PICs can be implemented over several machine stops

Particle-in-a-box model of one-dimensional excitons in conjugated polymers

Science.gov (United States)

Pedersen, Thomas G.; Johansen, Per M.; Pedersen, Henrik C.

2000-04-01

A simple two-particle model of excitons in conjugated polymers is proposed as an alternative to usual highly computationally demanding quantum chemical methods. In the two-particle model, the exciton is described as an electron-hole pair interacting via Coulomb forces and confined to the polymer backbone by rigid walls. Furthermore, by integrating out the transverse part, the two-particle equation is reduced to one-dimensional form. It is demonstrated how essentially exact solutions are obtained in the cases of short and long conjugation length, respectively. From a linear combination of these cases an approximate solution for the general case is obtained. As an application of the model the influence of a static electric field on the electron-hole overlap integral and exciton energy is considered.

Focusing of sub-micrometer particles and bacteria enabled by two-dimensional acoustophoresis

DEFF Research Database (Denmark)

Antfolk, M.; Muller, Peter Barkholt; Augustsson, P.

2014-01-01

Handling of sub-micrometer bioparticles such as bacteria are becoming increasingly important in the biomedical field and in environmental and food analysis. As a result, there is an increased need for less labor-intensive and time-consuming handling methods. Here, an acoustophoresis-based microfl......Handling of sub-micrometer bioparticles such as bacteria are becoming increasingly important in the biomedical field and in environmental and food analysis. As a result, there is an increased need for less labor-intensive and time-consuming handling methods. Here, an acoustophoresis......-based microfluidic chip that uses ultrasound to focus sub-micrometer particles and bacteria, is presented. The ability to focus sub-micrometer bioparticles in a standing one-dimensional acoustic wave is generally limited by the acoustic-streaming-induced drag force, which becomes increasingly significant the smaller...... particles as small as 0.5 μm in diameter in microchannels of square or rectangular cross sections, is demonstrated. Numerical analysis was used to determine generic transverse particle trajectories in the channels, which revealed spiral-shaped trajectories of the sub-micrometer particles towards the center...

Studies of Particle Wake Potentials in Plasmas

Science.gov (United States)

Ellis, Ian; Graziani, Frank; Glosli, James; Strozzi, David; Surh, Michael; Richards, David; Decyk, Viktor; Mori, Warren

2011-10-01

Fast Ignition studies require a detailed understanding of electron scattering, stopping, and energy deposition in plasmas with variable values for the number of particles within a Debye sphere. Presently there is disagreement in the literature concerning the proper description of these processes. Developing and validating proper descriptions requires studying the processes using first-principle electrostatic simulations and possibly including magnetic fields. We are using the particle-particle particle-mesh (PPPM) code ddcMD and the particle-in-cell (PIC) code BEPS to perform these simulations. As a starting point in our study, we examine the wake of a particle passing through a plasma in 3D electrostatic simulations performed with ddcMD and with BEPS using various cell sizes. In this poster, we compare the wakes we observe in these simulations with each other and predictions from Vlasov theory. Prepared by LLNL under Contract DE-AC52-07NA27344 and by UCLA under Grant DE-FG52-09NA29552.

Integrated Work Management: PIC, Course 31884

Energy Technology Data Exchange (ETDEWEB)

Simpson, Lewis Edward [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2017-09-08

The person-in-charge (PIC) plays a key role in the integrated work management (IWM) process at Los Alamos National Laboratory (LANL, or the Laboratory) because the PIC is assigned responsibility and authority by the responsible line manager (RLM) for the overall validation, coordination, release, execution, and closeout of a work activity in accordance with IWM. This course, Integrated Work Management: PIC (Course 31884), describes the PIC’s IWM roles and responsibilities. This course also discusses IWM requirements that the PIC must meet. For a general overview of the IWM process, see self-study Course 31881, Integrated Work Management: Overview. For instruction on the preparer’s role, see self-study Course 31883, Integrated Work Management: Preparer.

A deformable particle-in-cell method for advective transport in geodynamic modeling

Science.gov (United States)

Samuel, Henri

2018-06-01

This paper presents an improvement of the particle-in-cell method commonly used in geodynamic modeling for solving pure advection of sharply varying fields. Standard particle-in-cell approaches use particle kernels to transfer the information carried by the Lagrangian particles to/from the Eulerian grid. These kernels are generally one-dimensional and non-evolutive, which leads to the development of under- and over-sampling of the spatial domain by the particles. This reduces the accuracy of the solution, and may require the use of a prohibitive amount of particles in order to maintain the solution accuracy to an acceptable level. The new proposed approach relies on the use of deformable kernels that account for the strain history in the vicinity of particles. It results in a significant improvement of the spatial sampling by the particles, leading to a much higher accuracy of the numerical solution, for a reasonable computational extra cost. Various 2D tests were conducted to compare the performances of the deformable particle-in-cell method with the particle-in-cell approach. These consistently show that at comparable accuracy, the deformable particle-in-cell method was found to be four to six times more efficient than standard particle-in-cell approaches. The method could be adapted to 3D space and generalized to cases including motionless transport.

Application of a method for comparing one-dimensional and two-dimensional models of a ground-water flow system

International Nuclear Information System (INIS)

Naymik, T.G.

1978-01-01

To evaluate the inability of a one-dimensional ground-water model to interact continuously with surrounding hydraulic head gradients, simulations using one-dimensional and two-dimensional ground-water flow models were compared. This approach used two types of models: flow-conserving one-and-two dimensional models, and one-dimensional and two-dimensional models designed to yield two-dimensional solutions. The hydraulic conductivities of controlling features were varied and model comparison was based on the travel times of marker particles. The solutions within each of the two model types compare reasonably well, but a three-dimensional solution is required to quantify the comparison

Laser-driven magnetic reconnection in the multi-plasmoid regime

Science.gov (United States)

Totorica, Samuel; Abel, Tom; Fiuza, Frederico

2017-10-01

Magnetic reconnection is a promising candidate mechanism for accelerating the nonthermal particles associated with explosive astrophysical phenomena. Laboratory experiments are starting to probe multi-plasmoid regimes of relevance for particle acceleration. We have performed two- and three-dimensional particle-in-cell (PIC) simulations to explore particle acceleration for parameters relevant to laser-driven reconnection experiments. We have extended our previous work to explore particle acceleration in larger system sizes. Our results show the transition to plasmoid-dominated acceleration associated with the merging and contraction of plasmoids that further extend the maximum energy of the power-law tail of the particle distribution. Furthermore, we have modeled Coulomb collisions and will discuss the influence of collisionality on the plasmoid formation, dynamics, and particle acceleration.

Development of a relativistic Particle In Cell code PARTDYN for linear accelerator beam transport

Energy Technology Data Exchange (ETDEWEB)

Phadte, D., E-mail: deepraj@rrcat.gov.in [LPD, Raja Ramanna Centre for Advanced Technology, Indore 452013 (India); Patidar, C.B.; Pal, M.K. [MAASD, Raja Ramanna Centre for Advanced Technology, Indore (India)

2017-04-11

A relativistic Particle In Cell (PIC) code PARTDYN is developed for the beam dynamics simulation of z-continuous and bunched beams. The code is implemented in MATLAB using its MEX functionality which allows both ease of development as well higher performance similar to a compiled language like C. The beam dynamics calculations carried out by the code are compared with analytical results and with other well developed codes like PARMELA and BEAMPATH. The effect of finite number of simulation particles on the emittance growth of intense beams has been studied. Corrections to the RF cavity field expressions were incorporated in the code so that the fields could be calculated correctly. The deviations of the beam dynamics results between PARTDYN and BEAMPATH for a cavity driven in zero-mode have been discussed. The beam dynamics studies of the Low Energy Beam Transport (LEBT) using PARTDYN have been presented.

A dynamic neutral fluid model for the PIC scheme

Science.gov (United States)

Wu, Alan; Lieberman, Michael; Verboncoeur, John

2010-11-01

Fluid diffusion is an important aspect of plasma simulation. A new dynamic model is implemented using the continuity and boundary equations in OOPD1, an object oriented one-dimensional particle-in-cell code developed at UC Berkeley. The model is described and compared with analytical methods given in [1]. A boundary absorption parameter can be adjusted from ideal absorption to ideal reflection. Simulations exhibit good agreement with analytic time dependent solutions for the two ideal cases, as well as steady state solutions for mixed cases. For the next step, fluid sources and sinks due to particle-particle or particle-fluid collisions within the simulation volume and to surface reactions resulting in emission or absorption of fluid species will be implemented. The resulting dynamic interaction between particle and fluid species will be an improvement to the static fluid in the existing code. As the final step in the development, diffusion for multiple fluid species will be implemented. [4pt] [1] M.A. Lieberman and A.J. Lichtenberg, Principles of Plasma Discharges and Materials Processing, 2nd Ed, Wiley, 2005.

Particle-in-cell simulations of high energy electron production by intense laser pulses in underdense plasmas

International Nuclear Information System (INIS)

Susumu, Kato; Eisuke, Miura; Kazuyoshi, Koyama; Mitsumori, Tanimoto; Masahiro, Adachi

2004-01-01

The propagation of intense laser pulses and the generation of high energy electrons from underdense plasmas are investigated using two dimensional particle-in-cell simulations. When the ratio of the laser power to the critical power of relativistic self-focusing gets the optimal value, the laser pulse propagates in a steady way and electrons have maximum energies. (author)

Particle-in-cell simulations of high energy electron production by intense laser pulses in underdense plasmas

Energy Technology Data Exchange (ETDEWEB)

Susumu, Kato; Eisuke, Miura; Kazuyoshi, Koyama [National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki (Japan); Mitsumori, Tanimoto [Meisei Univ., Dept. of Electrical Engineering, Hino, Tokyo (Japan); Masahiro, Adachi [Hiroshima Univ., Graduate school of Advanced Science of Matter, Higashi-Hiroshima, Hiroshima (Japan)

2004-07-01

The propagation of intense laser pulses and the generation of high energy electrons from underdense plasmas are investigated using two dimensional particle-in-cell simulations. When the ratio of the laser power to the critical power of relativistic self-focusing gets the optimal value, the laser pulse propagates in a steady way and electrons have maximum energies. (author)

User's guide for TWOHEX: a code package for two-dimensional, neutral-particle transport in equilateral triangular meshes

International Nuclear Information System (INIS)

Walters, W.F.; Brinkley, F.W.; Marr, D.R.

1984-10-01

TWOHEX solves the two-dimensional multigroup transport equation on an equilateral triangular mesh in the x,y plane. Both regular and adjoint, inhomogeneous (fixed source) and homogeneous problems are solved. Three problem domains are treated by TWOHEX. The whole core domain is a 60 0 parallelogram with vacuum boundary conditions on each face. The third core domain is a 120 0 parallelogram with two vacuum and two rotational boundary conditions. The sixth core domain is a 60 0 parallelogram with two vacuum and two rotational boundary conditions. General anisotropic scattering is allowed and an anisotropic inhomogeneous source may be input as cell averages

Kawasaki dynamics with two types of particles : critical droplets

NARCIS (Netherlands)

Hollander, den W.Th.F.; Nardi, F.R.; Troiani, A.

2012-01-01

This is the third in a series of three papers in which we study a two-dimensional lattice gas consisting of two types of particles subject to Kawasaki dynamics at low temperature in a large finite box with an open boundary. Each pair of particles occupying neighboring sites has a negative binding

Kawasaki dynamics with two types of particles : critical droplets

NARCIS (Netherlands)

Hollander, den W.Th.F.; Nardi, F.R.; Troiani, A.

2012-01-01

This is the third in a series of three papers in which we study a two-dimensional lattice gas consisting of two types of particles subject to Kawasaki dynamics at low temperature in a large ¿nite box with an open boundary. Each pair of particles occupying neighboring sites has a negative binding

Two-dimensional particle simulation of negative ion extraction from a volume source

International Nuclear Information System (INIS)

Naitou, H.; Fukumasa, O.; Sakachou, K.; Mutou, K.

1995-01-01

Two-dimensional electrostatic particle simulation was done to study the extraction of negative ions from a volume plasma source. The simulation model is a rectangular system which consists of an extraction grid, a plasma grid, and a grounded wall. Full dynamics of electrons, ions, and negative ions are followed. Negative ions are extracted from the plasma region to the extraction grid through a slit in the plasma grid. For the lower value of extraction grid potential, the simulation results agree with the Child-Langumuir law, where the extracted negative ion current is proportional to the three-halves power of the potential of the extraction grid. For the higher value of extraction grid potential, the space charge effect of negative ions, which enter into the beamline at the top of the concavity of the positive ion boundary, reduces the negative ion current from the prediction of the Child-Langumuir law. ((orig.))

Impenetrable Mass-Imbalanced Particles in One-Dimensional Harmonic Traps

DEFF Research Database (Denmark)

Salami Dehkharghani, Amin; Volosniev, A. G.; Zinner, N. T.

2016-01-01

Strongly interacting particles in one dimension subject to external confinement have become a topic of considerable interest due to recent experimental advances and the development of new theoretical methods to attack such systems. In the case of equal mass fermions or bosons with two or more...... internal degrees of freedom, one can map the problem onto the well-known Heisenberg spin models. However, many interesting physical systems contain mixtures of particles with different masses. Therefore, a generalization of the recent strong-coupling techniques would be highly desirable....... This is particularly important since such problems are generally considered non-integrable and thus the hugely successful Bethe ansatz approach cannot be applied. Here we discuss some initial steps towards this goal by investigating small ensembles of one-dimensional harmonically trapped particles where pairwise...

X-ray studies on a two-dimensional diffusion limited system for cell growth

International Nuclear Information System (INIS)

Hlatky, L.R.; Alpen, E.L.

1985-01-01

X-ray studies were performed on cells grown in a new type of in vitro multicellular system, the ''sandwich'' system. This system is a two-dimensional array of cells, sandwiched between transparent slides which are impermeable to oxygen. The cell system is subject to self-created diffusion gradients of nutrients, metabolic products and, most importantly, oxygen. Sandwiches are analogous to living cross sections of multicellular spheroids or of poorly vascularized tumors. They contain a necrotic center, which the authors show to be due to diffusion limitations, an intermediate region which has a large fraction of quiescent cells and a cycling outer rim. One advantage sandwiches have over three-dimensional tumor models (sheproids) is that can control the amount of cell to cell contact and thereby separate effects due to oxygen or other gradients from effects due to contact. The authors present x-ray survival curves for sandwiches of various cell densities and compare them to x-ray survival curves done for spheroids and monolayers of the same cell line

Effects of increased vertebral number on carcass weight in PIC pigs.

Science.gov (United States)

Huang, Jieping; Zhang, Mingming; Ye, Runqing; Ma, Yun; Lei, Chuzhao

2017-12-01

Variation of the vertebral number is associated with carcass traits in pigs. However, results from different populations do not match well with others, especially for carcass weight. Therefore, effects of increased vertebral number on carcass weight were investigated by analyzing the relationship between two loci multi-vertebra causal loci (NR6A1 g.748 C > T and VRTN g.20311_20312ins291) and carcass weight in PIC pigs. Results from the association study between vertebral number and carcass weight showed that increased thoracic number had negative effects on carcass weight, but the results were not statistically significant. Further, VRTN Ins/Ins genotype increased more than one thoracic than that of Wt/Wt genotype on average in this PIC population. Meanwhile, there was a significant negative effect of VRTN Ins on carcass weight (P carcass weight in PIC pigs. © 2017 Japanese Society of Animal Science.

Metal Detector By Using PIC Microcontroller Interfacing With PC

OpenAIRE

Yin Min Theint; Myo Maung Maung; Hla Myo Tun

2015-01-01

Abstract This system proposes metal detector by using PIC microcontroller interfacing with PC. The system uses PIC microcontroller as the main controller whether the detected metal is ferrous metal or non-ferrous metal. Among various types of metal sensors and various types of metal detecting technologies concentric type induction coil sensor and VLF very low frequency metal detecting technology are used in this system. This system consists of two configurations Hardware configuration and Sof...

The MICHELLE 2D/3D ES PIC Code Advances and Applications

CERN Document Server

Petillo, John; De Ford, John F; Dionne, Norman J; Eppley, Kenneth; Held, Ben; Levush, Baruch; Nelson, Eric M; Panagos, Dimitrios; Zhai, Xiaoling

2005-01-01

MICHELLE is a new 2D/3D steady-state and time-domain particle-in-cell (PIC) code* that employs electrostatic and now magnetostatic finite-element field solvers. The code has been used to design and analyze a wide variety of devices that includes multistage depressed collectors, gridded guns, multibeam guns, annular-beam guns, sheet-beam guns, beam-transport sections, and ion thrusters. Latest additions to the MICHELLE/Voyager tool are as follows: 1) a prototype 3D self magnetic field solver using the curl-curl finite-element formulation for the magnetic vector potential, employing edge basis functions and accumulating current with MICHELLE's new unstructured grid particle tracker, 2) the electrostatic field solver now accommodates dielectric media, 3) periodic boundary conditions are now functional on all grids, not just structured grids, 4) the addition of a global optimization module to the user interface where both electrical parameters (such as electrode voltages)can be optimized, and 5) adaptive mesh ref...

Dispersion characteristics of two-dimensional unmagnetized dielectric plasma photonic crystal

International Nuclear Information System (INIS)

Li-Mei, Qi; Zi-Qiang, Yang; Feng, Lan; Xi, Gao; Da-Zhi, Li

2010-01-01

This paper studies dispersion characteristics of the transverse magnetic (TM) mode for two-dimensional unmagnetized dielectric plasma photonic crystal by a modified plane wave method. First, the cutoff behaviour is made clear by using the Maxwell–Garnett effective medium theory, and the influences of dielectric filling factor and dielectric constant on effective plasma frequency are analysed. Moreover, the occurence of large gaps in dielectric plasma photonic crystal is demonstrated by comparing the skin depth with the lattice constant, and the influence of plasma frequency on the first three gaps is also studied. Finally, by using the particle-in-cell simulation method, a transmission curve in the Γ – X direction is obtained in dielectric plasma photonic crystal, which is in accordance with the dispersion curves calculated by the modified plane wave method, and the large gap between the transmission points of 27 GHz and 47 GHz is explained by comparing the electric field patterns in particle-in-cell simulation

Differentiation Potential of Human Chorion-Derived Mesenchymal Stem Cells into Motor Neuron-Like Cells in Two- and Three-Dimensional Culture Systems.

Science.gov (United States)

Faghihi, Faezeh; Mirzaei, Esmaeil; Ai, Jafar; Lotfi, Abolfazl; Sayahpour, Forough Azam; Barough, Somayeh Ebrahimi; Joghataei, Mohammad Taghi

2016-04-01

Many people worldwide suffer from motor neuron-related disorders such as amyotrophic lateral sclerosis and spinal cord injuries. Recently, several attempts have been made to recruit stem cells to modulate disease progression in ALS and also regenerate spinal cord injuries. Chorion-derived mesenchymal stem cells (C-MSCs), used to be discarded as postpartum medically waste product, currently represent a class of cells with self renewal property and immunomodulatory capacity. These cells are able to differentiate into mesodermal and nonmesodermal lineages such as neural cells. On the other hand, gelatin, as a simply denatured collagen, is a suitable substrate for cell adhesion and differentiation. It has been shown that electrospinning of scaffolds into fibrous structure better resembles the physiological microenvironment in comparison with two-dimensional (2D) culture system. Since there is no report on potential of human chorion-derived MSCs to differentiate into motor neuron cells in two- and three-dimensional (3D) culture systems, we set out to determine the effect of retinoic acid (RA) and sonic hedgehog (Shh) on differentiation of human C-MSCs into motor neuron-like cells cultured on tissue culture plates (2D) and electrospun nanofibrous gelatin scaffold (3D).

Direct numerical simulations of agglomeration of circular colloidal particles in two-dimensional shear flow

International Nuclear Information System (INIS)

Choi, Young Joon; Djilali, Ned

2016-01-01

Colloidal agglomeration of nanoparticles in shear flow is investigated by solving the fluid-particle and particle-particle interactions in a 2D system. We use an extended finite element method in which the dynamics of the particles is solved in a fully coupled manner with the flow, allowing an accurate description of the fluid-particle interfaces without the need of boundary-fitted meshes or of empirical correlations to account for the hydrodynamic interactions between the particles. Adaptive local mesh refinement using a grid deformation method is incorporated with the fluid-structure interaction algorithm, and the particle-particle interaction at the microscopic level is modeled using the Lennard-Jones potential. Motivated by the process used in fabricating fuel cell catalysts from a colloidal ink, the model is applied to investigate agglomeration of colloidal particles under external shear flow in a sliding bi-periodic Lees-Edwards frame with varying shear rates and particle fraction ratios. Both external shear and particle fraction are found to have a crucial impact on the structure formation of colloidal particles in a suspension. Segregation intensity and graph theory are used to analyze the underlying agglomeration patterns and structures, and three agglomeration regimes are identified

Voro++: a three-dimensional Voronoi cell library in C++

Energy Technology Data Exchange (ETDEWEB)

Rycroft, Chris

2009-01-15

Voro++ is a free software library for the computation of three dimensional Voronoi cells. It is primarily designed for applications in physics and materials science, where the Voronoi tessellation can be a useful tool in the analysis of densely-packed particle systems, such as granular materials or glasses. The software comprises of several C++ classes that can be modified and incorporated into other programs. A command-line utility is also provided that can use most features of the code. Voro++ makes use of a direct cell-by-cell construction, which is particularly suited to handling special boundary conditions and walls. It employs algorithms which are tolerant for numerical precision errors, and it has been successfully employed on very large particle systems.

Ring accelerators

International Nuclear Information System (INIS)

Gisler, G.; Faehl, R.

1983-01-01

We present two-dimensional simulations in (r-z) and r-theta) cylinderical geometries of imploding-liner-driven accelerators of rings of charged particles. We address issues of azimuthal and longitudinal stability of the rings. We discuss self-trapping designs in which beam injection and extraction is aided by means of external cusp fields. Our simulations are done with the 2-1/2-D particle-in-cell plasma simulation code CLINER, which combines collisionless, electromagnetic PIC capabilities with a quasi-MHD finite element package

Waterlike anomalies in a two-dimensional core-softened potential

Science.gov (United States)

Bordin, José Rafael; Barbosa, Marcia C.

2018-02-01

We investigate the structural, thermodynamic, and dynamic behavior of a two-dimensional (2D) core-corona system using Langevin dynamics simulations. The particles are modeled by employing a core-softened potential which exhibits waterlike anomalies in three dimensions. In previous studies in a quasi-2D system a new region in the pressure versus temperature phase diagram of structural anomalies was observed. Here we show that for the two-dimensional case two regions in the pressure versus temperature phase diagram with structural, density, and diffusion anomalies are observed. Our findings indicate that, while the anomalous region at lower densities is due the competition between the two length scales in the potential at higher densities, the anomalous region is related to the reentrance of the melting line.

Monrelativistic particle in a magnetic field in two-dimensional Lobachevsky space, the cylindrical coordinates and the Poincare half-plane

International Nuclear Information System (INIS)

Ovsiyu, E.M.

2012-01-01

Exact solutions of the Schrodinger equation in the two-dimensional Riemannian space of negative curvature, the hyperbolic Lobachevsky plane, in the presence of an external magnetic field, which is an analog of a uniform magnetic field in the Minkowski space, are constructed. The description uses the cylindrical and quasi-Cartesian coordinates. The quasi-Cartesian coordinates determine the Poincare half-plane. In the both coordinate systems, the Schrodinger equation is solved exactly, the wave functions are constructed. A generalized formula for energy levels is found, which describes the quantized motion of a particle in a magnetic field in the Lobachevsky plane. (authors)

Expression of recombinant myostatin propeptide pPIC9K-Msp plasmid in Pichia pastoris.

Science.gov (United States)

Du, W; Xia, J; Zhang, Y; Liu, M J; Li, H B; Yan, X M; Zhang, J S; Li, N; Zhou, Z Y; Xie, W Z

2015-12-28

Myostatin propeptide can inhibit the biological activity of myostatin protein and promote muscle growth. To express myostatin propeptide in vitro with a higher biological activity, we performed codon optimization on the sheep myostatin propeptide gene sequence, and mutated aspartic acid-76 to alanine based on the codon usage bias of Pichia pastoris and the enhanced biological activity of myostatin propeptide mutant. Modified myostatin propeptide gene was cloned into the pPIC9K plasmid to form the recombinant plasmid pPIC9K-Msp. Recombinant plasmid pPIC9K-Msp was transformed into Pichia pastoris GS115 by electrotransformation. Transformed cells were screened, and methanol was used to induce expression. SDS-PAGE and western blotting were used to verify the successful expression of myostatin propeptide with biological activity in Pichia pastoris, providing the basis for characterization of this protein.

Particle acceleration in near critical density plasma

International Nuclear Information System (INIS)

Gu, Y.J.; Kong, Q.; Kawata, S.; Izumiyama, T.; Nagashima, T.

2013-01-01

Charged particle acceleration schemes driven by ultra intense laser and near critical density plasma interactions are presented. They include electron acceleration in a plasma channel, ion acceleration by the Coulomb explosion and high energy electron beam driven ion acceleration. It is found that under the near critical density plasma both ions and electrons are accelerated with a high acceleration gradient. The electron beam containing a large charge quantity is accelerated well with 23 GeV/cm. The collimated ion bunch reaches 1 GeV. The investigations and discussions are based on 2.5D PIC (particle-in-cell) simulations. (author)

Two-dimensional nano-lattice in Fe-Co-Ni-Al-Cu alloys

International Nuclear Information System (INIS)

Kalanov, M.U.; Ibragimova, E.M.; Khamraeva, R.N.; Rustamova, V.M.; Ummatov, H.D.

2007-01-01

Full text: The high coercive strength of the dispersionally solidified alloys on the base of Fe-Co-Ni-Al-Cu system appears as a result of the special thermomagnetic annealing, when particles of the strong magnetic phase are distinguished in non-magnetic matrix along an external magnetic field direction. The neutron studying allows one to reveal the correlation between magnetization and inclusion axes, and also existence of magnetic microcell and perfectness of the lattice. This work presents results of neutron diffraction study with a double-crystal spectrometer (0.145 nm). Plate like samples of size 18 12 4 mm 3 were cut from a single crystal of alloy UNDK35 T5 along (100) plane. Magnetic field of 6 kOe was applied perpendicular to the neutron beam. Zero-field spectrum had only random variation of the background. Under the applied magnetic field two maxima appeared at the angles of 12 and 24 minute. In the case of the magnetic field directed in parallel to the scattering vector, the two maxima disappeared as expected. It is evidence that nuclear scattering is less than magnetic one and the observed maxima correspond to (10) and (20) reflections from a two dimensional ferro-magnetic microcell. The cell parameter of the magnetic microcell was found 40.6 nm. The coherent scattering region size was 120-160 nm. The ferro-magnetic rod diameter estimated from the peak widths was 16 nm. The diffraction pattern for the demagnetized sample strongly differs from the initial magnetized sample, where a diffuse reflection was observed near Bragg reflection and related with residual magnetization. So, the magnetic inclusions created in the Fe-Co-Ni-Al-Cu system at the thermomagnetic annealing by means of disintegration of the solid solution are strong ferro-magnetic and one-domain. These particles form the two-dimensional magnetic microcell and interact each to other within 3-4 periods of the cell. (authors)

MEASURED DIAMETERS OF TWO F STARS IN THE β PIC MOVING GROUP

International Nuclear Information System (INIS)

Simon, M.; Schaefer, G. H.

2011-01-01

We report angular diameters of HIP 560 and HIP 21547, two F spectral-type pre-main-sequence members of the β Pic Moving Group. We used the east-west 314 m long baseline of the CHARA Array. The measured limb-darkened angular diameters of HIP 560 and HIP 21547 are 0.492 ± 0.032 and 0.518 ± 0.009 mas, respectively. The corresponding stellar radii are 2.1 and 1.6 R ☉ for HIP 560 and HIP 21547, respectively. These values indicate that the stars are truly young. Analyses using the evolutionary tracks calculated by Siess, Dufour, and Forestini and the tracks of the Yonsei-Yale group yield consistent results. Analyzing the measurements on an angular diameter versus color diagram we find that the ages of the two stars are indistinguishable; their average value is 13 ± 2 Myr. The masses of HIP 560 and HIP 21547 are 1.65 ± 0.02 and 1.75 ± 0.05 M ☉ , respectively. However, analysis of the stellar parameters on a Hertzsprung-Russell diagram yields ages at least 5 Myr older. Both stars are rapid rotators. The discrepancy between the two types of analyses has a natural explanation in gravitational darkening. Stellar oblateness, however, does not affect our measurements of angular diameters.

Photon management in two-dimensional disordered media.

Science.gov (United States)

Vynck, Kevin; Burresi, Matteo; Riboli, Francesco; Wiersma, Diederik S

2012-12-01

Elaborating reliable and versatile strategies for efficient light coupling between free space and thin films is of crucial importance for new technologies in energy efficiency. Nanostructured materials have opened unprecedented opportunities for light management, notably in thin-film solar cells. Efficient coherent light trapping has been accomplished through the careful design of plasmonic nanoparticles and gratings, resonant dielectric particles and photonic crystals. Alternative approaches have used randomly textured surfaces as strong light diffusers to benefit from their broadband and wide-angle properties. Here, we propose a new strategy for photon management in thin films that combines both advantages of an efficient trapping due to coherent optical effects and broadband/wide-angle properties due to disorder. Our approach consists of the excitation of electromagnetic modes formed by multiple light scattering and wave interference in two-dimensional random media. We show, by numerical calculations, that the spectral and angular responses of thin films containing disordered photonic patterns are intimately related to the in-plane light transport process and can be tuned through structural correlations. Our findings, which are applicable to all waves, are particularly suited for improving the absorption efficiency of thin-film solar cells and can provide a new approach for high-extraction-efficiency light-emitting diodes.

Collisional particle-in-cell modeling for energy transport accompanied by atomic processes in dense plasmas

Energy Technology Data Exchange (ETDEWEB)

Mishra, R.; Beg, F. N. [Center for Energy Research, University of California, San Diego, California 92093 (United States); Leblanc, P.; Sentoku, Y. [Department of Physics, University of Nevada, Reno, Nevada 89557 (United States); Wei, M. S. [General Atomics, San Diego, California 92121 (United States)

2013-07-15

Fully relativistic collisional Particle-in-Cell (PIC) code, PICLS, has been developed to study extreme energy density conditions produced in intense laser-solid interaction. Recent extensions to PICLS, such as the implementation of dynamic ionization, binary collisions in a partially ionized plasma, and radiative losses, enhance the efficacy of simulating intense laser plasma interaction and subsequent energy transport in resistive media. Different ionization models are introduced and benchmarked against each other to check the suitability of the model. The atomic physics models are critical to determine the energy deposition and transport in dense plasmas, especially when they consist of high Z (atomic number) materials. Finally we demonstrate the electron transport simulations to show the importance of target material on fast electron dynamics.

Four-dimensional optical manipulation of colloidal particles

DEFF Research Database (Denmark)

Rodrigo, P.J.; Daria, V.R.; Glückstad, J.

2005-01-01

We transform a TEM00 laser mode into multiple counterpropagating optical traps to achieve four-dimensional simultaneous manipulation of multiple particles. Efficient synthesis and dynamic control of the counterpropagating-beam traps is carried out via the generalized phase contrast method......, and a spatial polarization-encoding scheme. Our experiments genuinely demonstrate real-time, interactive particle-position control for forming arbitrary volumetric constellations and complex three-dimensional trajectories of multiple particles. This opens up doors for cross-disciplinary cutting-edge research...

[PICS: pharmaceutical inspection cooperation scheme].

Science.gov (United States)

Morénas, J

2009-01-01

The pharmaceutical inspection cooperation scheme (PICS) is a structure containing 34 participating authorities located worldwide (October 2008). It has been created in 1995 on the basis of the pharmaceutical inspection convention (PIC) settled by the European free trade association (EFTA) in1970. This scheme has different goals as to be an international recognised body in the field of good manufacturing practices (GMP), for training inspectors (by the way of an annual seminar and experts circles related notably to active pharmaceutical ingredients [API], quality risk management, computerized systems, useful for the writing of inspection's aide-memoires). PICS is also leading to high standards for GMP inspectorates (through regular crossed audits) and being a room for exchanges on technical matters between inspectors but also between inspectors and pharmaceutical industry.

Pseudo three-dimensional modeling of particle-fuel packing using distinct element method

International Nuclear Information System (INIS)

Yuki, Daisuke; Takata, Takashi; Yamaguchi, Akira

2007-01-01

Vibration-based packing of sphere-pac fuel is a key technology in a nuclear fuel manufacturing. In the production process of sphere-pac fuel, a Mixed Oxide (MOX) fuel is formed to spherical form and is packed in a cladding tube by adding a vibration force. In the present study, we have developed a numerical simulation method to investigate the behavior of the particles in a vibrated tube using the Distinct Element Method (DEM). In general, the DEM requires a significant computational cost. Therefore we propose a new approach in which a small particle can move through the space between three larger particles even in the two-dimensional simulation. We take into account an equivalent three-dimensional effect in the equations of motion. Thus it is named pseudo three-dimensional modeling. (author)

BOA, Beam Optics Analyzer A Particle-In-Cell Code

International Nuclear Information System (INIS)

Bui, Thuc

2007-01-01

The program was tasked with implementing time dependent analysis of charges particles into an existing finite element code with adaptive meshing, called Beam Optics Analyzer (BOA). BOA was initially funded by a DOE Phase II program to use the finite element method with adaptive meshing to track particles in unstructured meshes. It uses modern programming techniques, state-of-the-art data structures, so that new methods, features and capabilities are easily added and maintained. This Phase II program was funded to implement plasma simulations in BOA and extend its capabilities to model thermal electrons, secondary emissions, self magnetic field and implement a more comprehensive post-processing and feature-rich GUI. The program was successful in implementing thermal electrons, secondary emissions, and self magnetic field calculations. The BOA GUI was also upgraded significantly, and CCR is receiving interest from the microwave tube and semiconductor equipment industry for the code. Implementation of PIC analysis was partially successful. Computational resource requirements for modeling more than 2000 particles begin to exceed the capability of most readily available computers. Modern plasma analysis typically requires modeling of approximately 2 million particles or more. The problem is that tracking many particles in an unstructured mesh that is adapting becomes inefficient. In particular memory requirements become excessive. This probably makes particle tracking in unstructured meshes currently unfeasible with commonly available computer resources. Consequently, Calabazas Creek Research, Inc. is exploring hybrid codes where the electromagnetic fields are solved on the unstructured, adaptive mesh while particles are tracked on a fixed mesh. Efficient interpolation routines should be able to transfer information between nodes of the two meshes. If successfully developed, this could provide high accuracy and reasonable computational efficiency.

Assembly of positioner of automated two-dimensional scan coupled to X-ray fluorescence spectrometry

International Nuclear Information System (INIS)

Silva, Leonardo Santiago Melgaço

2011-01-01

This work describes the design and assembling of a prototype automated positioner two-dimensional scanning coupled to X-ray fluorescence spectrometry. The work aims to achieve a portable and easy to use, device of broad utility in the analysis of samples by X-ray fluorescence area of expertise and research. The two-dimensional scanning of the positioner is by means of two stepper motors controlled by a microcontroller PIC 16F877A, encoder and optical sensors. The user interacts with the XY table through an interface program for the Windows operating system, which communicates with the microcontroller through the serial port. The system of Fluorescence Spectroscopy incorporated into the positioner consists of a system commercially available system from the company AMPTEK, where the primary source of excitation of the sample was a source of 241 Am of 59.5 KeV emissions. Resolution and accuracy of tests were performed in the XY scanning process and reproducibility of the same kit with the fluorescence spectrometry X-ray. Qualitative tests by X-ray fluorescence spectrometry in samples were performed to demonstrate the applicability and versatility of the project. It follows that the prototype illustrates a possible adequately to portable device for X-ray spectrometry of two-dimensional. (author)

Particle-in-cell analysis of beam-wave interaction in gyrotron cavity with tapered magnetic field

Energy Technology Data Exchange (ETDEWEB)

Kumar, A., E-mail: anil.gyrotron@gmail.com [Gyrotron Lab., Microwave Tube Area, Central Electronics Engineering Research Inst. (CEERI, CSIR), Pilani, Rajasthan (India); Banasthali Univ., Dept. of Physics, Banasthali, Rajasthan (India); Khatun, H.; Kumar, N.; Singh, U.; Sinha, A.K. [Gyrotron Lab., Microwave Tube Area, Central Electronics Engineering Research Inst. (CEERI, CSIR), Pilani, Rajasthan (India); Vyas, V. [Banasthali Univ., Dept. of Physics, Banasthali, Rajasthan (India)

2010-11-15

A commercially available electromagnetic simulator -- MAGIC, a particle-in-cell (PIC) code -- has been used to carry out a comparative study of the beam-wave interaction under uniform and tapered magnetic field profiles of a 42 GHz, 200kW gyrotron. The magnetic field profile across the resonant cavity varies by ±6.5% with a peak value of 1.615 T. The MAGIC simulation shows the desire performance of the gyrotron under both magnetic field conditions with an operating mode TE{sub 03} and a pitch factor of 1.26. The analysis of the simulated results show that stability in the power growth was reached more quickly and achieved higher output power in the case of a tapered magnetic field. (author)

Two-dimensional boundary-value problem for ion-ion diffusion

International Nuclear Information System (INIS)

Tuszewski, M.; Lichtenberg, A.J.

1977-01-01

Like-particle diffusion is usually negligible compared with unlike-particle diffusion because it is two orders higher in spatial derivatives. When the ratio of the ion gyroradius to the plasma transverse dimension is of the order of the fourth root of the mass ratio, previous one-dimensional analysis indicated that like-particle diffusion is significant. A two-dimensional boundary-value problem for ion-ion diffusion is investigated. Numerical solutions are found with models for which the nonlinear partial differential equation reduces to an ordinary fourth-order differential equation. These solutions indicate that the ion-ion losses are higher by a factor of six for a slab geometry, and by a factor of four for circular geometry, than estimated from dimensional analysis. The solutions are applied to a multiple mirror experiment stabilized with a quadrupole magnetic field which generates highly elliptical flux surfaces. It is found that the ion-ion losses dominate the electron-ion losses and that these classical radial losses contribute to a significant decrease of plasma lifetime, in qualitiative agreement with the experimental results

Final Report for grant ER54958, 'Gyrokinetic Particle Simulation of Turbulent Transport in Burning Plasmas'

International Nuclear Information System (INIS)

Decyk, Viktor K.

2011-01-01

The UCLA contribution to this collaborative proposal is in two general areas. One area is to enhance the performance of GTC to perform optimally on the DOE leadership class computers, the other part is to contribute to the overall object-based design for GTC as it evolves. Most of the effort during this grant period has been in the former area. High performance computing is undergoing a revolution of greatly increasing parallelism, which is expected to lead to an exaflop computer by the end of the decade. However, the path there is uncertain. A number of hardware architectures have been proposed as well as new parallel languages. It is a disruptive time. In addition to this grant, this work is also funded by a gift from Northrop Grumman and UCLA's Institute for Digital Research and Education. In spite of the variety of proposed hardware, we feel that there is a hardware abstraction that describes the most likely features of the next generation high performance computers. This abstraction consists of a hierarchy of computational layers. At the lowest layer, we have a SIMD (vector) processor whose computational elements work together in lockstep and have fast synchronization and shared memory. At the next higher layer, we have a collection of such processors which communicate through a slower shared memory. At the highest layer, we have a cluster of such collections, which communicate via message-passing. The bottom two layers of this abstraction has been implemented in a language called OpenCL. Although different hardware implements different features of this abstraction, such as the SIMD or vector length, such an abstraction allows one to design parameterizable algorithms that can adapt to different hardware. We believe the hardware which most closely resembles a future exaflop computer is a collection of Graphical Processing Units (GPUs), and we have focused our attention on developing Particle-in-Cell (PIC) algorithms for this hardware initially. PIC codes are

Finding two-dimensional peaks

International Nuclear Information System (INIS)

Silagadze, Z.K.

2007-01-01

Two-dimensional generalization of the original peak finding algorithm suggested earlier is given. The ideology of the algorithm emerged from the well-known quantum mechanical tunneling property which enables small bodies to penetrate through narrow potential barriers. We merge this 'quantum' ideology with the philosophy of Particle Swarm Optimization to get the global optimization algorithm which can be called Quantum Swarm Optimization. The functionality of the newborn algorithm is tested on some benchmark optimization problems

Confinement limit of a Dirac particle in two and three dimensions

International Nuclear Information System (INIS)

Toyama, F.M.; Nogami, Y.

2010-01-01

Consider a particle that is in a stationary state described by the Dirac equation with a finite-range potential. In two and three dimensions the particle can be confined to an arbitrarily small spatial region. This is in contrast to the one-dimensional case in which the confinement region cannot be much narrower than the Compton wavelength.

Pattern formation in two-dimensional square-shoulder systems

International Nuclear Information System (INIS)

Fornleitner, Julia; Kahl, Gerhard

2010-01-01

Using a highly efficient and reliable optimization tool that is based on ideas of genetic algorithms, we have systematically studied the pattern formation of the two-dimensional square-shoulder system. An overwhelming wealth of complex ordered equilibrium structures emerge from this investigation as we vary the shoulder width. With increasing pressure three structural archetypes could be identified: cluster lattices, where clusters of particles occupy the sites of distorted hexagonal lattices, lane formation, and compact particle arrangements with high coordination numbers. The internal complexity of these structures increases with increasing shoulder width.

Pattern formation in two-dimensional square-shoulder systems

Energy Technology Data Exchange (ETDEWEB)

Fornleitner, Julia [Institut fuer Festkoerperforschung, Forschungsszentrum Juelich, D-52425 Juelich (Germany); Kahl, Gerhard, E-mail: fornleitner@cmt.tuwien.ac.a [Institut fuer Theoretische Physik and Centre for Computational Materials Science (CMS), Technische Universitaet Wien, Wiedner Hauptstrasse 8-10, A-1040 Wien (Austria)

2010-03-17

Using a highly efficient and reliable optimization tool that is based on ideas of genetic algorithms, we have systematically studied the pattern formation of the two-dimensional square-shoulder system. An overwhelming wealth of complex ordered equilibrium structures emerge from this investigation as we vary the shoulder width. With increasing pressure three structural archetypes could be identified: cluster lattices, where clusters of particles occupy the sites of distorted hexagonal lattices, lane formation, and compact particle arrangements with high coordination numbers. The internal complexity of these structures increases with increasing shoulder width.

Two-dimensional transport of tokamak plasmas

International Nuclear Information System (INIS)

Hirshman, S.P.; Jardin, S.C.

1979-01-01

A reduced set of two-fluid transport equations is obtained from the conservation equations describing the time evolution of the differential particle number, entropy, and magnetic fluxes in an axisymmetric toroidal plasma with nested magnetic surfaces. Expanding in the small ratio of perpendicular to parallel mobilities and thermal conductivities yields as solubility constraints one-dimensional equations for the surface-averaged thermodynamic variables and magnetic fluxes. Since Ohm's law E +u x B =R', where R' accounts for any nonideal effects, only determines the particle flow relative to the diffusing magnetic surfaces, it is necessary to solve a single two-dimensional generalized differential equation, (partial/partialt) delpsi. (delp - J x B) =0, to find the absolute velocity of a magnetic surface enclosing a fixed toroidal flux. This equation is linear but nonstandard in that it involves flux surface averages of the unknown velocity. Specification of R' and the cross-field ion and electron heat fluxes provides a closed system of equations. A time-dependent coordinate transformation is used to describe the diffusion of plasma quantities through magnetic surfaces of changing shape

Four-dimensional image display for associated particle imaging

International Nuclear Information System (INIS)

Headley, G.; Beyerle, A.; Durkee, R.; Hurley, P.; Tunnell, L.

1994-01-01

Associated particle imaging (API) is a three-dimensional neutron gamma imaging technique which provides both spatial and spectral information about an unknown. A local area network consisting of a UNIX fileserver and multiple DOS workstations has been chosen to perform the data acquisition and display functions. The data are acquired with a CAMAC system, stored in list mode, and sorted on the fileserver for display on the DOS workstations. Three of the display PCs, interacting with the fileserver, provide coordinated views as the operator ''slices'' the image. The operator has a choice of: a one-dimensional shadowgram from any side, two-dimensional shadowgrams from any side; a three-dimensional view (either perspective projection or stereoscopic). A common color scheme is used to carry energy information into the spatial images. ((orig.))

Two-dimensional model of a freely expanding plasma

International Nuclear Information System (INIS)

Khalid, Q.

1975-01-01

The free expansion of an initially confined plasma is studied by the computer experiment technique. The research is an extension to two dimensions of earlier work on the free expansion of a collisionless plasma in one dimension. In the two-dimensional rod model, developed in this research, the plasma particles, electrons and ions are modeled as infinitely long line charges or rods. The line charges move freely in two dimensions normal to their parallel axes, subject only to a self-consistent electric field. Two approximations, the grid approximation and the periodic boundary condition are made in order to reduce the computation time. In the grid approximation, the space occupied by the plasma at a given time is divided into boxes. The particles are subject to an average electric field calculated for that box assuming that the total charge within each box is located at the center of the box. However, the motion of each particle is exactly followed. The periodic boundary condition allows us to consider only one-fourth of the total number of particles of the plasma, representing the remaining three-fourths of the particles as symmetrically placed images of those whose positions are calculated. This approximation follows from the expected azimuthal symmetry of the plasma. The dynamics of the expansion are analyzed in terms of average ion and electron positions, average velocities, oscillation frequencies and relative distribution of energy between thermal, flow and electric field energies. Comparison is made with previous calculations of one-dimensional models which employed plane, spherical or cylindrical sheets as charged particles. In order to analyze the effect of the grid approximation, the model is solved for two different grid sizes and for each grid size the plasma dynamics is determined. For the initial phase of expansion, the agreement for the two grid sizes is found to be good

GAP--a PIC-type fluid code

International Nuclear Information System (INIS)

Marder, B.M.

1975-01-01

GAP, a PIC-type fluid code for computing compressible flows, is described and demonstrated. While retaining some features of PIC, it is felt that the GAP approach is conceptually and operationally simpler. 9 figures

Boundary effects in a quasi-two-dimensional driven granular fluid.

Science.gov (United States)

Smith, N D; Smith, M I

2017-12-01

The effect of a confining boundary on the spatial variations in granular temperature of a driven quasi-two-dimensional layer of particles is investigated experimentally. The radial drop in the relative granular temperature ΔT/T exhibits a maximum at intermediate particle numbers which coincides with a crossover from kinetic to collisional transport of energy. It is also found that at low particle numbers, the distributions of radial velocities are increasingly asymmetric as one approaches the boundary. The radial and tangential granular temperatures split, and in the tails of the radial velocity distribution there is a higher population of fast moving particles traveling away rather than towards the boundary.

Large gap plasma display cell with auxiliary electrodes: macro-cell experiments and two-dimensional modelling

International Nuclear Information System (INIS)

Ouyang, J T; Callegari, Th; Caillier, B; Boeuf, J-P

2003-01-01

In this paper we use a two-dimensional fluid model and a 'macroscopic' PDP cell to investigate the possibility of using large gap configurations with auxiliary electrodes to improve the efficiency of PDP discharge cells. The large gap allows operation in a transient positive column regime where energy is more efficiently deposited into xenon excitation, while the auxiliary electrodes are used to keep reasonable values of the operating voltage. Two types of auxiliary electrode configurations (floating and powered) are considered. The discharge characteristics and the discharge efficiency in exciting xenon are studied with simulations and by measuring the intensity of infrared emission from xenon and visible emission from neon in a macroscopic PDP cell. The results show that an efficient positive column regime can be achieved at reasonably low operating voltages when the auxiliary electrode configuration is carefully designed

PIC Activation through Functional Interplay between Mediator and TFIIH.

Science.gov (United States)

Malik, Sohail; Molina, Henrik; Xue, Zhu

2017-01-06

The multiprotein Mediator coactivator complex functions in large part by controlling the formation and function of the promoter-bound preinitiation complex (PIC), which consists of RNA polymerase II and general transcription factors. However, precisely how Mediator impacts the PIC, especially post-recruitment, has remained unclear. Here, we have studied Mediator effects on basal transcription in an in vitro transcription system reconstituted from purified components. Our results reveal a close functional interplay between Mediator and TFIIH in the early stages of PIC development. We find that under conditions when TFIIH is not normally required for transcription, Mediator actually represses transcription. TFIIH, whose recruitment to the PIC is known to be facilitated by the Mediator, then acts to relieve Mediator-induced repression to generate an active form of the PIC. Gel mobility shift analyses of PICs and characterization of TFIIH preparations carrying mutant XPB translocase subunit further indicate that this relief of repression is achieved through expending energy via ATP hydrolysis, suggesting that it is coupled to TFIIH's established promoter melting activity. Our interpretation of these results is that Mediator functions as an assembly factor that facilitates PIC maturation through its various stages. Whereas the overall effect of the Mediator is to stimulate basal transcription, its initial engagement with the PIC generates a transcriptionally inert PIC intermediate, which necessitates energy expenditure to complete the process. Copyright © 2016 Elsevier Ltd. All rights reserved.

Logarithmic Superdiffusion in Two Dimensional Driven Lattice Gases

Science.gov (United States)

Krug, J.; Neiss, R. A.; Schadschneider, A.; Schmidt, J.

2018-03-01

The spreading of density fluctuations in two-dimensional driven diffusive systems is marginally anomalous. Mode coupling theory predicts that the diffusivity in the direction of the drive diverges with time as (ln t)^{2/3} with a prefactor depending on the macroscopic current-density relation and the diffusion tensor of the fluctuating hydrodynamic field equation. Here we present the first numerical verification of this behavior for a particular version of the two-dimensional asymmetric exclusion process. Particles jump strictly asymmetrically along one of the lattice directions and symmetrically along the other, and an anisotropy parameter p governs the ratio between the two rates. Using a novel massively parallel coupling algorithm that strongly reduces the fluctuations in the numerical estimate of the two-point correlation function, we are able to accurately determine the exponent of the logarithmic correction. In addition, the variation of the prefactor with p provides a stringent test of mode coupling theory.

Multichannel scattering amplitudes of microparticles in a quantum well with two-dimensional -potential

International Nuclear Information System (INIS)

Sedrakian, D.M.; Badalyan, D.H.; Sedrakian, L.R.

2015-01-01

Quasi-one-dimensional quantum particle scattering on two-dimensional δ-potential is considered. Analytical expressions for the amplitudes of the multi-channel transmission and reflection are given. The problem for the case when the number of channels is finite and equal N, and the particle falls on the potential moving through the channel l is solved. The case of a three channel scattering is studied in details. It is shown that under conditions k 2 → 0 and k 3 → 0 'overpopulation' of particles on the second and third channels occurs. The points of δ-potential location which provide a full 'overpopulation' of particles is also found

Designing Embedded Systems with PIC Microcontrollers Principles and Applications

CERN Document Server

Wilmshurst, Tim

2009-01-01

PIC microcontrollers are used worldwide in commercial and industrial devices. The 8-bit PIC which this book focuses on is a versatile work horse that completes many designs. An engineer working with applications that include a microcontroller will no doubt come across the PIC sooner rather than later. It is a must to have a working knowledge of this 8-bit technology. This book takes the novice from introduction of embedded systems through to advanced development techniques for utilizing and optimizing the PIC family of microcontrollers in your device. To truly understand the PIC, assembly and

Hybrid-PIC Computer Simulation of the Plasma and Erosion Processes in Hall Thrusters

Science.gov (United States)

Hofer, Richard R.; Katz, Ira; Mikellides, Ioannis G.; Gamero-Castano, Manuel

2010-01-01

HPHall software simulates and tracks the time-dependent evolution of the plasma and erosion processes in the discharge chamber and near-field plume of Hall thrusters. HPHall is an axisymmetric solver that employs a hybrid fluid/particle-in-cell (Hybrid-PIC) numerical approach. HPHall, originally developed by MIT in 1998, was upgraded to HPHall-2 by the Polytechnic University of Madrid in 2006. The Jet Propulsion Laboratory has continued the development of HPHall-2 through upgrades to the physical models employed in the code, and the addition of entirely new ones. Primary among these are the inclusion of a three-region electron mobility model that more accurately depicts the cross-field electron transport, and the development of an erosion sub-model that allows for the tracking of the erosion of the discharge chamber wall. The code is being developed to provide NASA science missions with a predictive tool of Hall thruster performance and lifetime that can be used to validate Hall thrusters for missions.

Stable high efficiency two-dimensional perovskite solar cells via cesium doping

KAUST Repository

Zhang, Xu

2017-08-15

Two-dimensional (2D) organic-inorganic perovskites have recently emerged as one of the most important thin-film solar cell materials owing to their excellent environmental stability. The remaining major pitfall is their relatively poor photovoltaic performance in contrast to 3D perovskites. In this work we demonstrate cesium cation (Cs) doped 2D (BA)(MA)PbI perovskite solar cells giving a power conversion efficiency (PCE) as high as 13.7%, the highest among the reported 2D devices, with excellent humidity resistance. The enhanced efficiency from 12.3% (without Cs) to 13.7% (with 5% Cs) is attributed to perfectly controlled crystal orientation, an increased grain size of the 2D planes, superior surface quality, reduced trap-state density, enhanced charge-carrier mobility and charge-transfer kinetics. Surprisingly, it is found that the Cs doping yields superior stability for the 2D perovskite solar cells when subjected to a high humidity environment without encapsulation. The device doped using 5% Cs degrades only ca. 10% after 1400 hours of exposure in 30% relative humidity (RH), and exhibits significantly improved stability under heating and high moisture environments. Our results provide an important step toward air-stable and fully printable low dimensional perovskites as a next-generation renewable energy source.

Size-sensitive particle trajectories in three-dimensional micro-bubble acoustic streaming flows

Science.gov (United States)

Volk, Andreas; Rossi, Massimiliano; Hilgenfeldt, Sascha; Rallabandi, Bhargav; Kähler, Christian; Marin, Alvaro

2015-11-01

Oscillating microbubbles generate steady streaming flows with interesting features and promising applications for microparticle manipulation. The flow around oscillating semi-cylindrical bubbles has been typically assumed to be independent of the axial coordinate. However, it has been recently revealed that particle motion is strongly three-dimensional: Small tracer particles follow vortical trajectories with pronounced axial displacements near the bubble, weaving a toroidal stream-surface. A well-known consequence of bubble streaming flows is size-dependent particle migration, which can be exploited for sorting and trapping of microparticles in microfluidic devices. In this talk, we will show how the three-dimensional toroidal topology found for small tracer particles is modified as the particle size increases up to 1/3 of the bubble radius. Our results show size-sensitive particle positioning along the axis of the semi-cylindrical bubble. In order to analyze the three-dimensional sorting and trapping capabilities of the system, experiments with an imposed flow and polydisperse particle solutions are also shown.

Comparison of Yeast Cell Protein Solubilization Procedures for Two-dimensional Electrophoresis

DEFF Research Database (Denmark)

Harder, A; Wildgruber, R; Nawrocki, A

1999-01-01

Three different procedures for the solubilization of yeast (S. cerevisiae) cell proteins were compared on the basis of the obtained two-dimensional (2-D) polypeptide patterns. Major emphasis was laid on minimizing handling steps, protein modification or degradation, and quantitative loss of high...... with sodium dodecyl sulfate (SDS) buffer, consisting of 1% SDS and 100 mM tris(hydroxymethyl)aminomethane (Tris)-HCl, pH 7.0, followed by dilution with "standard" lysis buffer, and (iii) boiling the sample with SDS during cell lysis, followed by dilution with thiourea/urea lysis buffer (2 M thiourea/ 7 M urea...

Aspects of jamming in two-dimensional athermal frictionless systems.

Science.gov (United States)

Reichhardt, C; Reichhardt, C J Olson

2014-05-07

In this work we provide an overview of jamming transitions in two dimensional systems focusing on the limit of frictionless particle interactions in the absence of thermal fluctuations. We first discuss jamming in systems with short range repulsive interactions, where the onset of jamming occurs at a critical packing density and where certain quantities show a divergence indicative of critical behavior. We describe how aspects of the dynamics change as the jamming density is approached and how these dynamics can be explored using externally driven probes. Different particle shapes can produce jamming densities much lower than those observed for disk-shaped particles, and we show how jamming exhibits fragility for some shapes while for other shapes this is absent. Next we describe the effects of long range interactions and jamming behavior in systems such as charged colloids, vortices in type-II superconductors, and dislocations. We consider the effect of adding obstacles to frictionless jamming systems and discuss connections between this type of jamming and systems that exhibit depinning transitions. Finally, we discuss open questions such as whether the jamming transition in all these different systems can be described by the same or a small subset of universal behaviors, as well as future directions for studies of jamming transitions in two dimensional systems, such as jamming in self-driven or active matter systems.

PICS bags safely store unshelled and shelled groundnuts in Niger.

Science.gov (United States)

Baributsa, D; Baoua, I B; Bakoye, O N; Amadou, L; Murdock, L L

2017-05-01

We conducted an experiment in Niger to evaluate the performance of hermetic triple layer (Purdue Improved Crop Storage- PICS) bags for the preservation of shelled and unshelled groundnut Arachis hypogaea L. Naturally-infested groundnut was stored in PICS bags and woven bags for 6.7 months. After storage, the average oxygen level in the PICS bags fell from 21% to 18% (v/v) and 21%-15% (v/v) for unshelled and shelled groundnut, respectively. Identified pests present in the stored groundnuts were Tribolium castaneum (Herbst), Corcyra cephalonica (Stainton) and Cryptolestes ferrugineus (Stephens). After 6.7 months of storage, in the woven bag, there was a large increase in the pest population accompanied by a weight loss of 8.2% for unshelled groundnuts and 28.7% for shelled groundnut. In PICS bags for both shelled and unshelled groundnuts, by contrast, the density of insect pests did not increase, there was no weight loss, and the germination rate was the same compared to that recorded at the beginning of the experiment. Storing shelled groundnuts in PICS bags is the most cost-effective way as it increases the quantity of grain stored.

Three-dimensional investigation of the two-phase flow structure in a bubbly pipe flow

International Nuclear Information System (INIS)

Schmidl, W.; Hassan, Y.A.; Ortiz-Villafuerte, J.

1996-01-01

Particle image velocimetry (PIV) is a nonintrusive measurement technique that can be used to study the structure of various fluid flows. PIV is used to measure the time-varying, full-field velocity data of a particle-seeded flow field within either a two-dimensional plane or three-dimensional volume. PIV is a very efficient measurement technique since it can obtain both qualitative and quantitative spatial information about the flow field being studied. The quantitative spatial velocity information can be further processed into information of flow parameters such as vorticity and turbulence over extended areas. The objective of this study was to apply recent advances and improvements in the PIV flow measurement technique to the full-field, nonintrusive analysis of a three-dimensional, two-phase fluid flow system in such a manner that both components of the two-phase system could be experimentally quantified

Blended particle filters for large-dimensional chaotic dynamical systems

Science.gov (United States)

Majda, Andrew J.; Qi, Di; Sapsis, Themistoklis P.

2014-01-01

A major challenge in contemporary data science is the development of statistically accurate particle filters to capture non-Gaussian features in large-dimensional chaotic dynamical systems. Blended particle filters that capture non-Gaussian features in an adaptively evolving low-dimensional subspace through particles interacting with evolving Gaussian statistics on the remaining portion of phase space are introduced here. These blended particle filters are constructed in this paper through a mathematical formalism involving conditional Gaussian mixtures combined with statistically nonlinear forecast models compatible with this structure developed recently with high skill for uncertainty quantification. Stringent test cases for filtering involving the 40-dimensional Lorenz 96 model with a 5-dimensional adaptive subspace for nonlinear blended filtering in various turbulent regimes with at least nine positive Lyapunov exponents are used here. These cases demonstrate the high skill of the blended particle filter algorithms in capturing both highly non-Gaussian dynamical features as well as crucial nonlinear statistics for accurate filtering in extreme filtering regimes with sparse infrequent high-quality observations. The formalism developed here is also useful for multiscale filtering of turbulent systems and a simple application is sketched below. PMID:24825886

Fast ignition: Dependence of the ignition energy on source and target parameters for particle-in-cell-modelled energy and angular distributions of the fast electrons

Energy Technology Data Exchange (ETDEWEB)

Bellei, C.; Divol, L.; Kemp, A. J.; Key, M. H.; Larson, D. J.; Strozzi, D. J.; Marinak, M. M.; Tabak, M.; Patel, P. K. [Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550 (United States)

2013-05-15

The energy and angular distributions of the fast electrons predicted by particle-in-cell (PIC) simulations differ from those historically assumed in ignition designs of the fast ignition scheme. Using a particular 3D PIC calculation, we show how the ignition energy varies as a function of source-fuel distance, source size, and density of the pre-compressed fuel. The large divergence of the electron beam implies that the ignition energy scales with density more weakly than the ρ{sup −2} scaling for an idealized beam [S. Atzeni, Phys. Plasmas 6, 3316 (1999)], for any realistic source that is at some distance from the dense deuterium-tritium fuel. Due to the strong dependence of ignition energy with source-fuel distance, the use of magnetic or electric fields seems essential for the purpose of decreasing the ignition energy.

Relativistic two-and three-particle scattering equations using instant and light-front dynamics

International Nuclear Information System (INIS)

Adhikari, S.K.; Tomio, L.; Frederico, T.

1992-01-01

Starting from the Bethe-Salpeter equation for two particles in the ladder approximation and integrating over the time component of momentum we derive three dimensional scattering integral equations satisfying constraints of unitarity and relativity, both employing the light-front and instant-form variables. The equations we arrive at are those first derived by Weinberg and by Blankenbecler and Sugar, and are shown to be related by a transformation of variables. Hence we show how to perform and relate identical dynamical calculation using these two equations. We extends this procedure to the case of three particles interacting via two-particle separable potentials. Using light-front and instant form variables we suggest a couple of three dimensional three-particle scattering equations satisfying constraints of two and three-particle unitarity and relativity. The three-particle light-front equation is shown to be approximately related by a transformation of variables to one of the instant-form three-particle equations. (author)

Particle-two particle interaction in configuration space

International Nuclear Information System (INIS)

Kuzmichev, V.E.

1982-07-01

The problem if three indentical particles with zero-range two-particle interaction is considered. An explicit expression for the effective potential between one particle and the remaining two-particle system is obtained in the coordinate representation. It is shown that for arbitrary energies, at small and, for zero energy, at large distances rho between the one particle and centre of mass of the other two particles the diagonal matrix element of the effective potential is attractive and proportional to 1/rho 2 . This property of the effective potenial explains both the Thomas singularity and the Efimov effect. In the case of zero total energy of the system the general form of the solution of the three-particle integral equation is found in configuration space. (orig.)

Engineering two-photon high-dimensional states through quantum interference

Science.gov (United States)

Zhang, Yingwen; Roux, Filippus S.; Konrad, Thomas; Agnew, Megan; Leach, Jonathan; Forbes, Andrew

2016-01-01

Many protocols in quantum science, for example, linear optical quantum computing, require access to large-scale entangled quantum states. Such systems can be realized through many-particle qubits, but this approach often suffers from scalability problems. An alternative strategy is to consider a lesser number of particles that exist in high-dimensional states. The spatial modes of light are one such candidate that provides access to high-dimensional quantum states, and thus they increase the storage and processing potential of quantum information systems. We demonstrate the controlled engineering of two-photon high-dimensional states entangled in their orbital angular momentum through Hong-Ou-Mandel interference. We prepare a large range of high-dimensional entangled states and implement precise quantum state filtering. We characterize the full quantum state before and after the filter, and are thus able to determine that only the antisymmetric component of the initial state remains. This work paves the way for high-dimensional processing and communication of multiphoton quantum states, for example, in teleportation beyond qubits. PMID:26933685

Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

International Nuclear Information System (INIS)

Li, Fei; Yu, Peicheng; Xu, Xinlu; Fiuza, Frederico; Decyk, Viktor K.

2017-01-01

In this study we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift (1^ direction). We show that this eliminates the main NCI modes with moderate |k_1|, while keeps additional main NCI modes well outside the range of physical interest with higher |k_1|. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along 1^ which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss’ Law is satisfied. Lastly, we present simulation examples of one relativistically drifting plasma, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.

Repulsion of polarized particles from two-dimensional materials

Science.gov (United States)

Rodríguez-Fortuño, Francisco J.; Picardi, Michela F.; Zayats, Anatoly V.

2018-05-01

Repulsion of nanoparticles, molecules, and atoms from surfaces can have important applications in nanomechanical devices, microfluidics, optical manipulation, and atom optics. Here, through the solution of a classical scattering problem, we show that a dipole source oscillating at a frequency ω can experience a robust and strong repulsive force when its near-field interacts with a two-dimensional material. As an example, the case of graphene is considered, showing that a broad bandwidth of repulsion can be obtained at frequencies for which propagation of plasmon modes is allowed 0 chemical potential tunable electrically or by chemical doping.

Two-dimensional nanoparticle self-assembly using plasma-induced Ostwald ripening

International Nuclear Information System (INIS)

Tang, J; Photopoulos, P; Tsoukalas, D; Tserepi, A

2011-01-01

In this work, a novel Ag nanoparticle self-assembly process based on plasma-induced two-dimensional Ostwald ripening is demonstrated. Ag nanoparticles are deposited on p-doped Si substrates using a DC magnetron sputtering process. With the assistance of O 2 /Ar plasma treatment, different sizes and patterns of Ag nanoparticles are formed, due to the Ostwald ripening. The evolution of plasma-induced nanoparticle ripening is studied and a clear increase in particle size and a decrease in particle density are observed with increasing plasma treatment. From the experiments, it is concluded that the initial nanoparticle density and the plasma gas mixture (Ar/O 2 ratio) are important factors that affect the ripening process. The proposed plasma-directed Ag nanoparticle self-assembly provides a rapid method of tailoring the nanoparticle distribution on substrates, with potential applications in the fields of solar cells, biosensors, and catalysis.

Low-noise electromagnetic δf particle-in-cell simulation of electron Bernstein waves

International Nuclear Information System (INIS)

Xiang Nong; Cary, John R.; Barnes, Daniel C.; Carlsson, John

2006-01-01

The conversion of the extraordinary (X) mode to an electron Bernstein wave (EBW) is one way to get rf energy into an overdense plasma. Analysis of this is complex, as the EBW is a fully kinetic wave, and so its linear propagation is described by an intractable integro-differential equation. Nonlinear effects cannot be calculated within this rubric at all. Full particle-in-cell (PIC) simulations cannot be used for these analyses, as the noise levels for reasonable simulation parameters are much greater than the typical rf amplitudes. It is shown that the delta-f computations are effective for this analysis. In particular, the accuracy of those computations has been verified by comparison with full PIC, cold plasma theory, and small gyroradius theory. This computational method is then used to analyze mode conversion in different frequency regimes. In particular, reasonable agreement with the theoretical predictions of Ram and Schultz [Phys. Plasmas 7, 4084 (2000)] in the linear regime is found, where 100% X-B mode conversion has been obtained when the driving frequency is less than twice the electron gyrofrequency. The results show that cold-plasma theory well predicts the mode conversion efficiency, as is consistent with the phase-space picture of mode conversion. From this it can be shown that nearly 100% X-B mode conversion cannot be obtained when the frequency is higher than the electron second harmonic cyclotron frequency

Self-organized defect strings in two-dimensional crystals.

Science.gov (United States)

Lechner, Wolfgang; Polster, David; Maret, Georg; Keim, Peter; Dellago, Christoph

2013-12-01

Using experiments with single-particle resolution and computer simulations we study the collective behavior of multiple vacancies injected into two-dimensional crystals. We find that the defects assemble into linear strings, terminated by dislocations with antiparallel Burgers vectors. We show that these defect strings propagate through the crystal in a succession of rapid one-dimensional gliding and rare rotations. While the rotation rate decreases exponentially with the number of defects in the string, the diffusion constant is constant for large strings. By monitoring the separation of the dislocations at the end points, we measure their effective interactions with high precision beyond their spontaneous formation and annihilation, and we explain the double-well form of the dislocation interaction in terms of continuum elasticity theory.

Three Dimensional Simulation of Ion Thruster Plume-Spacecraft Interaction Based on a Graphic Processor Unit

International Nuclear Information System (INIS)

Ren Junxue; Xie Kan; Qiu Qian; Tang Haibin; Li Juan; Tian Huabing

2013-01-01

Based on the three-dimensional particle-in-cell (PIC) method and Compute Unified Device Architecture (CUDA), a parallel particle simulation code combined with a graphic processor unit (GPU) has been developed for the simulation of charge-exchange (CEX) xenon ions in the plume of an ion thruster. Using the proposed technique, the potential and CEX plasma distribution are calculated for the ion thruster plume surrounding the DS1 spacecraft at different thrust levels. The simulation results are in good agreement with measured CEX ion parameters reported in literature, and the GPU's results are equal to a CPU's. Compared with a single CPU Intel Core 2 E6300, 16-processor GPU NVIDIA GeForce 9400 GT indicates a speedup factor of 3.6 when the total macro particle number is 1.1×10 6 . The simulation results also reveal how the back flow CEX plasma affects the spacecraft floating potential, which indicates that the plume of the ion thruster is indeed able to alleviate the extreme negative floating potentials of spacecraft in geosynchronous orbit

New relativistic particle-in-cell simulation studies of prompt and early afterglows from GRBs

International Nuclear Information System (INIS)

Ken-Ichi Nishikawa

2008-01-01

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electro-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the collisionless relativistic shock particle acceleration is due to plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electrons' transverse deflection behind the jet head. The '' jitter '' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. (author)

Particle simulation algorithms with short-range forces in MHD and fluid flow

International Nuclear Information System (INIS)

Cable, S.; Tajima, T.; Umegaki, K.

1992-07-01

Attempts are made to develop numerical algorithms for handling fluid flows involving liquids and liquid-gas mixtures. In these types of systems, the short-range intermolecular interactions are important enough to significantly alter behavior predicted on the basis of standard fluid mechanics and magnetohydrodynamics alone. We have constructed a particle-in-cell (PIC) code for the purpose of studying the effects of these interactions. Of the algorithms considered, the one which has been successfully implemented is based on a MHD particle code developed by Brunel et al. In the version presented here, short range forces are included in particle motion by, first, calculating the forces between individual particles and then, to prevent aliasing, interpolating these forces to the computational grid points, then interpolating the forces back to the particles. The code has been used to model a simple two-fluid Rayleigh-Taylor instability. Limitations to the accuracy of the code exist at short wavelengths, where the effects of the short-range forces would be expected to be most pronounced

(Nearly) portable PIC code for parallel computers

International Nuclear Information System (INIS)

Decyk, V.K.

1993-01-01

As part of the Numerical Tokamak Project, the author has developed a (nearly) portable, one dimensional version of the GCPIC algorithm for particle-in-cell codes on parallel computers. This algorithm uses a spatial domain decomposition for the fields, and passes particles from one domain to another as the particles move spatially. With only minor changes, the code has been run in parallel on the Intel Delta, the Cray C-90, the IBM ES/9000 and a cluster of workstations. After a line by line translation into cmfortran, the code was also run on the CM-200. Impressive speeds have been achieved, both on the Intel Delta and the Cray C-90, around 30 nanoseconds per particle per time step. In addition, the author was able to isolate the data management modules, so that the physics modules were not changed much from their sequential version, and the data management modules can be used as open-quotes black boxes.close quotes

Plasma density enhancement in atmospheric-pressure dielectric-barrier discharges by high-voltage nanosecond pulse in the pulse-on period: a PIC simulation

International Nuclear Information System (INIS)

Sang Chaofeng; Sun Jizhong; Wang Dezhen

2010-01-01

A particle-in-cell (PIC) plus Monte Carlo collision simulation is employed to investigate how a sustainable atmospheric pressure single dielectric-barrier discharge responds to a high-voltage nanosecond pulse (HVNP) further applied to the metal electrode. The results show that the HVNP can significantly increase the plasma density in the pulse-on period. The ion-induced secondary electrons can give rise to avalanche ionization in the positive sheath, which widens the discharge region and enhances the plasma density drastically. However, the plasma density stops increasing as the applied pulse lasts over certain time; therefore, lengthening the pulse duration alone cannot improve the discharge efficiency further. Physical reasons for these phenomena are then discussed.

Plasma density enhancement in atmospheric-pressure dielectric-barrier discharges by high-voltage nanosecond pulse in the pulse-on period: a PIC simulation

Science.gov (United States)

Sang, Chaofeng; Sun, Jizhong; Wang, Dezhen

2010-02-01

A particle-in-cell (PIC) plus Monte Carlo collision simulation is employed to investigate how a sustainable atmospheric pressure single dielectric-barrier discharge responds to a high-voltage nanosecond pulse (HVNP) further applied to the metal electrode. The results show that the HVNP can significantly increase the plasma density in the pulse-on period. The ion-induced secondary electrons can give rise to avalanche ionization in the positive sheath, which widens the discharge region and enhances the plasma density drastically. However, the plasma density stops increasing as the applied pulse lasts over certain time; therefore, lengthening the pulse duration alone cannot improve the discharge efficiency further. Physical reasons for these phenomena are then discussed.

Three-dimensional two-phase mass transport model for direct methanol fuel cells

International Nuclear Information System (INIS)

Yang, W.W.; Zhao, T.S.; Xu, C.

2007-01-01

A three-dimensional (3D) steady-state model for liquid feed direct methanol fuel cells (DMFC) is presented in this paper. This 3D mass transport model is formed by integrating five sub-models, including a modified drift-flux model for the anode flow field, a two-phase mass transport model for the porous anode, a single-phase model for the polymer electrolyte membrane, a two-phase mass transport model for the porous cathode, and a homogeneous mist-flow model for the cathode flow field. The two-phase mass transport models take account the effect of non-equilibrium evaporation/ condensation at the gas-liquid interface. A 3D computer code is then developed based on the integrated model. After being validated against the experimental data reported in the literature, the code was used to investigate numerically transport behaviors at the DMFC anode and their effects on cell performance

Solving the two-dimensional Schrödinger equation using basis ...

Indian Academy of Sciences (India)

Ihab H Naeim

2017-10-19

Oct 19, 2017 ... We shall study the case of a two-dimensional quantum system .... Solving (6) for ck,l is tantamount to pro- ... case, the final computational problem becomes quite ..... matrix approach fails in the case of two particles con-.

Two-dimensional analytical model of a proton exchange membrane fuel cell

International Nuclear Information System (INIS)

Liu, Jia Xing; Guo, Hang; Ye, Fang; Ma, Chong Fang

2017-01-01

In this study, a two-dimensional full cell analytical model of a proton exchange membrane fuel cell is developed. The analytical model describes electrochemical reactions on the anode and cathode catalyst layer, reactants diffusion in the gas diffusion layer, and gases flow in the gas channel, etc. The analytical solution is derived according to the basic physical equations. The performance predicted by the model is in good agreement with the experimental data. The results show that the polarization mainly occurs in the cathode side of the proton exchange membrane fuel cell. The anodic overpotential cannot be neglected. The hydrogen and oxygen concentrations decrease along the channel flow direction. The hydrogen and oxygen concentrations in the catalyst layer decrease with the current density. As predicted by the model, concentration polarization mainly occurs in the cathode side. - Highlights: • A 2D full cell analytical model of a proton exchange membrane fuel cell is developed. • The analytical solution is deduced according to the basic equations. • The anode overpotential is not so small that it cannot be neglected. • Species concentration distributions in the fuel cell is obtained and analyzed.

Three-Dimensional Simulation of Ultrasound-Induced Microalgal Cell Disruption.

Science.gov (United States)

Wang, M; Yuan, W; Hale, Andy

2016-03-01

The three-dimensional distribution (x, y, and z) of ultrasound-induced microalgal cell disruption in a sonochemical reactor was predicted by solving the Helmholtz equation using a three-dimensional acoustic module in the COMSOL Multiphysics software. The simulated local ultrasound pressure at any given location (x, y, and z) was found to correlate with cell disruption of a freshwater alga, Scenedesmus dimorphus, represented by the change of algal cell particle/debris concentration, chlorophyll-a fluorescence density (CAFD), and Nile red stained lipid fluorescence density (LFD), which was also validated by the model reaction of potassium iodide oxidation (the Weissler reaction). Furthermore, the effect of ultrasound power intensity and processing duration on algal cell disruption was examined to address the limitation of the model.

Molecular Dynamic Studies of Particle Wake Potentials in Plasmas

Science.gov (United States)

Ellis, Ian; Graziani, Frank; Glosli, James; Strozzi, David; Surh, Michael; Richards, David; Decyk, Viktor; Mori, Warren

2010-11-01

Fast Ignition studies require a detailed understanding of electron scattering, stopping, and energy deposition in plasmas with variable values for the number of particles within a Debye sphere. Presently there is disagreement in the literature concerning the proper description of these processes. Developing and validating proper descriptions requires studying the processes using first-principle electrostatic simulations and possibly including magnetic fields. We are using the particle-particle particle-mesh (P^3M) code ddcMD to perform these simulations. As a starting point in our study, we examined the wake of a particle passing through a plasma. In this poster, we compare the wake observed in 3D ddcMD simulations with that predicted by Vlasov theory and those observed in the electrostatic PIC code BEPS where the cell size was reduced to .03λD.

Characterizing a New Candidate Benchmark Brown Dwarf Companion in the β Pic Moving Group

Science.gov (United States)

Phillips, Caprice; Bowler, Brendan; Liu, Michael C.; Mace, Gregory N.; Sokal, Kimberly R.

2018-01-01

Benchmark brown dwarfs are objects that have at least two measured fundamental quantities such as luminosity and age, and therefore can be used to test substellar atmospheric and evolutionary models. Nearby, young, loose associations such as the β Pic moving group represent some of the best regions in which to identify intermediate-age benchmark brown dwarfs due to their well-constrained ages and metallicities. We present a spectroscopic study of a new companion at the hydrogen-burning limit orbiting a low-mass star at a separation of 9″ (650 AU) in the 23 Myr old β Pic moving group. The medium-resolution near-infrared spectrum of this companion from IRTF/SpeX shows clear signs of low surface gravity and yields an index-based spectral type of M6±1 with a VL-G gravity on the Allers & Liu classification system. Currently, there are four known brown dwarf and giant planet companions in the β Pic moving group: HR 7329 B, PZ Tel B, β Pic b, and 51 Eri b. Depending on its exact age and accretion history, this new object may represent the third brown dwarf companion and fifth substellar companion in this association.

Performance of PICS bags under extreme conditions in the sahel zone of Niger.

Science.gov (United States)

Baoua, Ibrahim B; Bakoye, Ousmane; Amadou, Laouali; Murdock, Larry L; Baributsa, Dieudonne

2018-03-01

Experiments in Niger assessed whether extreme environmental conditions including sunlight exposure affect the performance of triple-layer PICS bags in protecting cowpea grain against bruchids. Sets of PICS bags and woven polypropylene bags as controls containing 50 kg of naturally infested cowpea grain were held in the laboratory or outside with sun exposure for four and one-half months. PICS bags held either inside or outside exhibited no significant increase in insect damage and no loss in weight after 4.5 months of storage compared to the initial values. By contrast, woven bags stored inside or outside side by side with PICS bags showed several-fold increases in insects present in or on the grain and significant losses in grain weight. Grain stored inside in PICS bags showed no reduction in germination versus the initial value but there was a small but significant drop in germination of grain in PICS bags held outside (7.6%). Germination rates dropped substantially more in grain stored in woven bags inside (16.1%) and still more in woven bags stored outside (60%). PICS bags held inside and outside retained their ability to maintain internal reduced levels of oxygen and elevated levels of carbon dioxide. Exposure to extreme environmental conditions degraded the external polypropylene outer layer of the PICS triple-layer bag. Even so, the internal layers of polyethylene were more slowly degraded. The effects of exposure to sunlight, temperature and humidity variation within the sealed bags are described.

Four-dimensional Hall mechanics as a particle on CP3

International Nuclear Information System (INIS)

Bellucci, Stefano; Casteill, Pierre-Yves; Nersessian, Armen

2003-01-01

In order to establish an explicit connection between four-dimensional Hall effect on S 4 and six-dimensional Hall effect on CP 3 , we perform the Hamiltonian reduction of a particle moving on CP 3 in a constant magnetic field to the four-dimensional Hall mechanics (i.e., a-bar particle on S 4 in a SU(2) instanton field). This reduction corresponds to fixing the isospin of the latter system

Two-dimensional impurity transport calculations for a high recycling divertor

International Nuclear Information System (INIS)

Brooks, J.N.

1986-04-01

Two dimensional analysis of impurity transport in a high recycling divertor shows asymmetric particle fluxes to the divertor plate, low helium pumping efficiency, and high scrapeoff zone shielding for sputtered impurities

Segregation in quasi-two-dimensional granular systems

International Nuclear Information System (INIS)

Rivas, Nicolas; Cordero, Patricio; Soto, Rodrigo; Risso, Dino

2011-01-01

Segregation for two granular species is studied numerically in a vertically vibrated quasi-two-dimensional (quasi-2D) box. The height of the box is smaller than two particle diameters so that particles are limited to a submonolayer. Two cases are considered: grains that differ in their density but have equal size, and grains that have equal density but different diameters, while keeping the quasi-2D condition. It is observed that in both cases, for vibration frequencies beyond a certain threshold-which depends on the density or diameter ratios-segregation takes place in the lateral directions. In the quasi-2D geometry, gravity does not play a direct role in the in-plane dynamics and gravity does not point to the segregation directions; hence, several known segregation mechanisms that rely on gravity are discarded. The segregation we observe is dominated by a lack of equipartition between the two species; the light particles exert a larger pressure than the heavier ones, inducing the latter to form clusters. This energy difference in the horizontal direction is due to the existence of a fixed point characterized by vertical motion and hence vanishing horizontal energy. Heavier and bigger grains are more rapidly attracted to the fixed point and the perturbations are less efficient in taking them off the fixed point when compared to the lighter grains. As a consequence, heavier and bigger grains have less horizontal agitation than lighter ones. Although limited by finite size effects, the simulations suggest that the two cases we consider differ in the transition character: one is continuous and the other is discontinuous. In the cases where grains differ in mass on varying the control parameter, partial segregation is first observed, presenting many clusters of heavier particles. Eventually, a global cluster is formed with impurities; namely lighter particles are present inside. The transition looks continuous when characterized by several segregation order

Equilibrium spherically curved two-dimensional Lennard-Jones systems

NARCIS (Netherlands)

Voogd, J.M.; Sloot, P.M.A.; van Dantzig, R.

2005-01-01

To learn about basic aspects of nano-scale spherical molecular shells during their formation, spherically curved two-dimensional N-particle Lennard-Jones systems are simulated, studying curvature evolution paths at zero-temperature. For many N-values (N < 800) equilibrium configu- rations are traced

Resistive-strips micromegas detectors with two-dimensional readout

Science.gov (United States)

Byszewski, M.; Wotschack, J.

2012-02-01

Micromegas detectors show very good performance for charged particle tracking in high rate environments as for example at the LHC. It is shown that two coordinates can be extracted from a single gas gap in these detectors. Several micromegas chambers with spark protection by resistive strips and two-dimensional readout have been tested in the context of the R&D work for the ATLAS Muon System upgrade.

Magnetoresistance oscillations of two-dimensional electron systems in lateral superlattices with structured unit cells

Science.gov (United States)

Gerhardts, Rolf R.

2015-11-01

Model calculations for commensurability oscillations of the low-field magnetoresistance of two-dimensional electron systems (2DES) in lateral superlattices, consisting of unit cells with an internal structure, are compared with recent experiments. The relevant harmonics of the effective modulation potential depend not only on the geometrical structure of the modulated unit cell, but also strongly on the nature of the modulation. While higher harmonics of an electrostatically generated surface modulation are exponentially damped at the position of the 2DES about 90 nm below the surface, no such damping appears for strain-induced modulation generated, e.g., by the deposition of stripes of calixarene resist on the surface before cooling down the sample.

Optical image encryption based on phase retrieval combined with three-dimensional particle-like distribution

International Nuclear Information System (INIS)

Chen, Wen; Chen, Xudong; Sheppard, Colin J R

2012-01-01

We propose a new phase retrieval algorithm for optical image encryption in three-dimensional (3D) space. The two-dimensional (2D) plaintext is considered as a series of particles distributed in 3D space, and an iterative phase retrieval algorithm is developed to encrypt the series of particles into phase-only masks. The feasibility and effectiveness of the proposed method are demonstrated by a numerical experiment, and the advantages and security of the proposed optical cryptosystems are also analyzed and discussed. (paper)

Merging for Particle-Mesh Complex Particle Kinetic Modeling of the Multiple Plasma Beams

Science.gov (United States)

Lipatov, Alexander S.

2011-01-01

We suggest a merging procedure for the Particle-Mesh Complex Particle Kinetic (PMCPK) method in case of inter-penetrating flow (multiple plasma beams). We examine the standard particle-in-cell (PIC) and the PMCPK methods in the case of particle acceleration by shock surfing for a wide range of the control numerical parameters. The plasma dynamics is described by a hybrid (particle-ion-fluid-electron) model. Note that one may need a mesh if modeling with the computation of an electromagnetic field. Our calculations use specified, time-independent electromagnetic fields for the shock, rather than self-consistently generated fields. While a particle-mesh method is a well-verified approach, the CPK method seems to be a good approach for multiscale modeling that includes multiple regions with various particle/fluid plasma behavior. However, the CPK method is still in need of a verification for studying the basic plasma phenomena: particle heating and acceleration by collisionless shocks, magnetic field reconnection, beam dynamics, etc.

Generation of two-dimensional binary mixtures in complex plasmas

Science.gov (United States)

Wieben, Frank; Block, Dietmar

2016-10-01

Complex plasmas are an excellent model system for strong coupling phenomena. Under certain conditions the dust particles immersed into the plasma form crystals which can be analyzed in terms of structure and dynamics. Previous experiments focussed mostly on monodisperse particle systems whereas dusty plasmas in nature and technology are polydisperse. Thus, a first and important step towards experiments in polydisperse systems are binary mixtures. Recent experiments on binary mixtures under microgravity conditions observed a phase separation of particle species with different radii even for small size disparities. This contradicts several numerical studies of 2D binary mixtures. Therefore, dedicated experiments are required to gain more insight into the physics of polydisperse systems. In this contribution first ground based experiments on two-dimensional binary mixtures are presented. Particular attention is paid to the requirements for the generation of such systems which involve the consideration of the temporal evolution of the particle properties. Furthermore, the structure of these two-component crystals is analyzed and compared to simulations. This work was supported by the Deutsche Forschungsgemeinschaft DFG in the framework of the SFB TR24 Greifswald Kiel, Project A3b.