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.

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.)

Low-noise Collision Operators for Particle-in-cell Simulations

International Nuclear Information System (INIS)

Lewandowski, J.L.V.

2005-01-01

A new method to implement low-noise collision operators in particle-in-cell simulations is presented. The method is based on the fact that relevant collision operators can be included naturally in the Lagrangian formulation that exemplifies the particle-in-cell simulation method. Numerical simulations show that the momentum and energy conservation properties of the simulated plasma associated with the low-noise collision operator are improved as compared with standard collision algorithms based on random numbers

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

Investigation of Rising-Sun Magnetrons Operated at Relativistic Voltages Using Three Dimensional Particle-in-Cell Simulation

International Nuclear Information System (INIS)

Lemke, R.W.; Genoni, T.C.; Spencer, T.A.

1999-01-01

This work is an attempt to elucidate effects that may limit efficiency in magnetrons operated at relativistic voltages (V ∼ 500 kV). Three-dimensional particle-in-cell simulation is used to investigate the behavior of 14 and 22 cavity, cylindrical, rising-sun magnetrons. Power is extracted radially through a single iris located at the end of every other cavity. Numerical results show that in general output power and efficiency increase approximately linearly with increasing iris width (decreasing vacuum Q) until the total Q becomes too low for stable oscillation in the n-mode to be maintained. Beyond this point mode competition and/or switching occur and efficiency decreases. Results reveal that the minimum value of Q (maximum efficiency) that can be achieved prior to the onset of mode competition is significantly affected by the magnitude of the 0-space-harmonic of the π-mode, a unique characteristic of rising-suns, and by the magnitude of the electron current density (space-charge effects). By minimizing these effects, up to 3.7 GW output power has been produced at an efficiency of 40%

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)

Particle-in-cell simulations of anomalous transport in a Penning discharge

Science.gov (United States)

Carlsson, Johan; Kaganovich, Igor; Powis, Andrew; Raitses, Yevgeny; Romadanov, Ivan; Smolyakov, Andrei

2018-06-01

Electrostatic particle-in-cell simulations of a Penning discharge are performed in order to investigate azimuthally asymmetric, spoke-like structures previously observed in experiments. Two-dimensional simulations show that for Penning-discharge conditions, a persistent nonlinear spoke-like structure forms readily and rotates in the direction of E × B and electron diamagnetic drifts. The azimuthal velocity is within about a factor of 2 of the ion acoustic speed. The spoke frequency follows the experimentally observed scaling with ion mass, which indicates the importance of ion inertia in spoke formation. The spoke provides enhanced (anomalous) radial electron transport, and the effective cross-field conductivity is several times larger than the classical (collisional) value. The level of anomalous current obtained in the simulations is in good agreement with the experimental data. The rotating spoke channels most of the radial current, observable by an edge probe as short pulses.

Quasi-One-Dimensional Particle-in-Cell Simulation of Magnetic Nozzles

Science.gov (United States)

Ebersohn, Frans H.; Sheehan, J. P.; Gallimore, Alec D.; Shebalin, John V.

2015-01-01

A method for the quasi-one-dimensional simulation of magnetic nozzles is presented and simulations of a magnetic nozzle are performed. The effects of the density variation due to plasma expansion and the magnetic field forces on ion acceleration are investigated. Magnetic field forces acting on the electrons are found to be responsible for the formation of potential structures which accelerate ions. The effects of the plasma density variation alone are found to only weakly affect ion acceleration. Strongly diverging magnetic fields drive more rapid potential drops.

A general concurrent algorithm for plasma particle-in-cell simulation codes

International Nuclear Information System (INIS)

Liewer, P.C.; Decyk, V.K.

1989-01-01

We have developed a new algorithm for implementing plasma particle-in-cell (PIC) simulation codes on concurrent processors with distributed memory. This algorithm, named the general concurrent PIC algorithm (GCPIC), has been used to implement an electrostatic PIC code on the 33-node JPL Mark III Hypercube parallel computer. To decompose at PIC code using the GCPIC algorithm, the physical domain of the particle simulation is divided into sub-domains, equal in number to the number of processors, such that all sub-domains have roughly equal numbers of particles. For problems with non-uniform particle densities, these sub-domains will be of unequal physical size. Each processor is assigned a sub-domain and is responsible for updating the particles in its sub-domain. This algorithm has led to a a very efficient parallel implementation of a well-benchmarked 1-dimensional PIC code. The dominant portion of the code, updating the particle positions and velocities, is nearly 100% efficient when the number of particles is increased linearly with the number of hypercube processors used so that the number of particles per processor is constant. For example, the increase in time spent updating particles in going from a problem with 11,264 particles run on 1 processor to 360,448 particles on 32 processors was only 3% (parallel efficiency of 97%). Although implemented on a hypercube concurrent computer, this algorithm should also be efficient for PIC codes on other parallel architectures and for large PIC codes on sequential computers where part of the data must reside on external disks. copyright 1989 Academic Press, Inc

ASPEN: A fully kinetic, reduced-description particle-in-cell model for simulating parametric instabilities

International Nuclear Information System (INIS)

Vu, H.X.; Bezzerides, B.; DuBois, D.F.

1999-01-01

A fully kinetic, reduced-description particle-in-cell (RPIC) model is presented in which deviations from quasineutrality, electron and ion kinetic effects, and nonlinear interactions between low-frequency and high-frequency parametric instabilities are modeled correctly. The model is based on a reduced description where the electromagnetic field is represented by three separate temporal envelopes in order to model parametric instabilities with low-frequency and high-frequency daughter waves. Because temporal envelope approximations are invoked, the simulation can be performed on the electron time scale instead of the time scale of the light waves. The electrons and ions are represented by discrete finite-size particles, permitting electron and ion kinetic effects to be modeled properly. The Poisson equation is utilized to ensure that space-charge effects are included. The RPIC model is fully three dimensional and has been implemented in two dimensions on the Accelerated Strategic Computing Initiative (ASCI) parallel computer at Los Alamos National Laboratory, and the resulting simulation code has been named ASPEN. The authors believe this code is the first particle-in-cell code capable of simulating the interaction between low-frequency and high-frequency parametric instabilities in multiple dimensions. Test simulations of stimulated Raman scattering, stimulated Brillouin scattering, and Langmuir decay instability are presented

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

Object-Oriented Parallel Particle-in-Cell Code for Beam Dynamics Simulation in Linear Accelerators

International Nuclear Information System (INIS)

Qiang, J.; Ryne, R.D.; Habib, S.; Decky, V.

1999-01-01

In this paper, we present an object-oriented three-dimensional parallel particle-in-cell code for beam dynamics simulation in linear accelerators. A two-dimensional parallel domain decomposition approach is employed within a message passing programming paradigm along with a dynamic load balancing. Implementing object-oriented software design provides the code with better maintainability, reusability, and extensibility compared with conventional structure based code. This also helps to encapsulate the details of communications syntax. Performance tests on SGI/Cray T3E-900 and SGI Origin 2000 machines show good scalability of the object-oriented code. Some important features of this code also include employing symplectic integration with linear maps of external focusing elements and using z as the independent variable, typical in accelerators. A successful application was done to simulate beam transport through three superconducting sections in the APT linac design

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

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.

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.

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

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.

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

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

Particle-in-Cell Laser-Plasma Simulation on Xeon Phi Coprocessors

OpenAIRE

Surmin, I. A.; Bastrakov, S. I.; Efimenko, E. S.; Gonoskov, A. A.; Korzhimanov, A. V.; Meyerov, I. B.

2015-01-01

This paper concerns development of a high-performance implementation of the Particle-in-Cell method for plasma simulation on Intel Xeon Phi coprocessors. We discuss suitability of the method for Xeon Phi architecture and present our experience of porting and optimization of the existing parallel Particle-in-Cell code PICADOR. Direct porting with no code modification gives performance on Xeon Phi close to 8-core CPU on a benchmark problem with 50 particles per cell. We demonstrate step-by-step...

Three-dimensional simulation of the motion of a single particle under a simulated turbulent velocity field

Science.gov (United States)

Moreno-Casas, P. A.; Bombardelli, F. A.

2015-12-01

A 3D Lagrangian particle tracking model is coupled to a 3D channel velocity field to simulate the saltation motion of a single sediment particle moving in saltation mode. The turbulent field is a high-resolution three dimensional velocity field that reproduces a by-pass transition to turbulence on a flat plate due to free-stream turbulence passing above de plate. In order to reduce computational costs, a decoupled approached is used, i.e., the turbulent flow is simulated independently from the tracking model, and then used to feed the 3D Lagrangian particle model. The simulations are carried using the point-particle approach. The particle tracking model contains three sub-models, namely, particle free-flight, a post-collision velocity and bed representation sub-models. The free-flight sub-model considers the action of the following forces: submerged weight, non-linear drag, lift, virtual mass, Magnus and Basset forces. The model also includes the effect of particle angular velocity. The post-collision velocities are obtained by applying conservation of angular and linear momentum. The complete model was validated with experimental results from literature within the sand range. Results for particle velocity time series and distribution of particle turbulent intensities are presented.

IBSIMU: A three-dimensional simulation software for charged particle optics

Energy Technology Data Exchange (ETDEWEB)

Kalvas, T.; Tarvainen, O.; Ropponen, T.; Steczkiewicz, O.; Aerje, J. [Department of Physics, University of Jyvaeskylae, Jyvaeskylae 40500 (Finland); Clark, H. [Cyclotron Institute, Texas A and M University, College Station, Texas 77843 (United States)

2010-02-15

A general-purpose three-dimensional (3D) simulation code IBSIMU for charged particle optics with space charge is under development at JYFL. The code was originally developed for designing a slit-beam plasma extraction and nanosecond scale chopping for pulsed neutron generator, but has been developed further and has been used for many applications. The code features a nonlinear FDM Poisson's equation solver based on fast stabilized biconjugate gradient method with ILU0 preconditioner for solving electrostatic fields. A generally accepted nonlinear plasma model is used for plasma extraction. Magnetic fields can be imported to the simulations from other programs. The particle trajectories are solved using adaptive Runge-Kutta method. Steady-state and time-dependent problems can be modeled in cylindrical coordinates, two-dimensional (slit) geometry, or full 3D. The code is used via C++ programming language for versatility but it features an interactive easy-to-use postprocessing tool for diagnosing fields and particle trajectories. The open source distribution and public documentation make the code well suited for scientific use. IBSIMU has been used for modeling the 14 GHz ECR ion source extraction and for designing a four-electrode extraction for a 2.45 GHz microwave ion source at Jyvaeskylae. A grid extraction has also been designed for producing large uniform beam for creating conditions similar to solar wind. The code has also been used to design a H{sup -} extraction with electron dumping for the Cyclotron Institute of Texas A and M University.

Three dimensional particle simulation of drift wave fluctuations in a sheared magnetic field

International Nuclear Information System (INIS)

Sydora, R.D.; Leboeuf, J.N.; Thayer, D.R.; Diamond, P.H.; Tajima, T.

1985-08-01

Three dimensional particle simulations of collisionless drift waves in sheared magnetic fields were performed in order to determine the nonlinear behavior of inverse electron resonance dynamics in the presence of thermal fluctuations. It is found that stochastic electron diffusion in the electron resonance overlap region can destabilize the drift wave eigenmodes. Numerical evaluations based on a nonlinear electron resonance broadening theory give predictions in accord with the frequency and growth rates found in the simulation of short wavelength modes (k/sub y/rho/sub s/ greater than or equal to1)

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)

Two dimensional numerical simulation of gas discharges: comparison between particle-in-cell and FCT techniques

Energy Technology Data Exchange (ETDEWEB)

Soria-Hoyo, C; Castellanos, A [Departamento de Electronica y Electromagnetismo, Facultad de Fisica, Universidad de Sevilla, Avda. Reina Mercedes s/n, 41012 Sevilla (Spain); Pontiga, F [Departamento de Fisica Aplicada II, EUAT, Universidad de Sevilla, Avda. Reina Mercedes s/n, 41012 Sevilla (Spain)], E-mail: cshoyo@us.es

2008-10-21

Two different numerical techniques have been applied to the numerical integration of equations modelling gas discharges: a finite-difference flux corrected transport (FD-FCT) technique and a particle-in-cell (PIC) technique. The PIC technique here implemented has been specifically designed for the simulation of 2D electrical discharges using cylindrical coordinates. The development and propagation of a streamer between two parallel electrodes has been used as a convenient test to compare the performance of both techniques. In particular, the phase velocity of the cathode directed streamer has been used to check the internal consistency of the numerical simulations. The results obtained from the two techniques are in reasonable agreement with each other, and both techniques have proved their ability to follow the high gradients of charge density and electric field present in this type of problems. Moreover, the streamer velocities predicted by the simulation are in accordance with the typical experimental values.

Two dimensional numerical simulation of gas discharges: comparison between particle-in-cell and FCT techniques

International Nuclear Information System (INIS)

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

2008-01-01

Two different numerical techniques have been applied to the numerical integration of equations modelling gas discharges: a finite-difference flux corrected transport (FD-FCT) technique and a particle-in-cell (PIC) technique. The PIC technique here implemented has been specifically designed for the simulation of 2D electrical discharges using cylindrical coordinates. The development and propagation of a streamer between two parallel electrodes has been used as a convenient test to compare the performance of both techniques. In particular, the phase velocity of the cathode directed streamer has been used to check the internal consistency of the numerical simulations. The results obtained from the two techniques are in reasonable agreement with each other, and both techniques have proved their ability to follow the high gradients of charge density and electric field present in this type of problems. Moreover, the streamer velocities predicted by the simulation are in accordance with the typical experimental values.

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

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 ...

Particle-in-cell simulation of x-ray wakefield acceleration and betatron radiation in nanotubes

Directory of Open Access Journals (Sweden)

Xiaomei Zhang

2016-10-01

Full Text Available Though wakefield acceleration in crystal channels has been previously proposed, x-ray wakefield acceleration has only recently become a realistic possibility since the invention of the single-cycled optical laser compression technique. We investigate the acceleration due to a wakefield induced by a coherent, ultrashort x-ray pulse guided by a nanoscale channel inside a solid material. By two-dimensional particle-in-cell computer simulations, we show that an acceleration gradient of TeV/cm is attainable. This is about 3 orders of magnitude stronger than that of the conventional plasma-based wakefield accelerations, which implies the possibility of an extremely compact scheme to attain ultrahigh energies. In addition to particle acceleration, this scheme can also induce the emission of high energy photons at ∼O(10–100  MeV. Our simulations confirm such high energy photon emissions, which is in contrast with that induced by the optical laser driven wakefield scheme. In addition to this, the significantly improved emittance of the energetic electrons has been discussed.

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.

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)

On the stimulated Raman sidescattering in inhomogeneous plasmas: revisit of linear theory and three-dimensional particle-in-cell simulations

Science.gov (United States)

Xiao, C. Z.; Zhuo, H. B.; Yin, Y.; Liu, Z. J.; Zheng, C. Y.; Zhao, Y.; He, X. T.

2018-02-01

Stimulated Raman sidescattering (SRSS) in inhomogeneous plasma is comprehensively revisited on both theoretical and numerical aspects due to the increasing concern of its detriments to inertial confinement fusion. Firstly, two linear mechanisms of finite beam width and collisional effects that could suppress SRSS are investigated theoretically. Thresholds for the eigenmode and wave packet in a finite-width beam are derived as a supplement to the theory proposed by Mostrom and Kaufman (1979 Phys. Rev. Lett. 42 644). Collisional absorption of SRSS is efficient at high-density plasma and high-Z material, otherwise, it allows emission of sidescattering. Secondly, we have performed the first three-dimensional particle-in-cell simulations in the context of SRSS to investigate its linear and nonlinear effects. Simulation results are qualitatively agreed with the linear theory. SRSS with the maximum growth gain is excited at various densities, grows to an amplitude that is comparable with the pump laser, and evolutes to lower densities with a large angle of emergence. Competitions between SRSS and other parametric instabilities such as stimulated Raman backscattering, two-plasmon decay, and stimulated Brillouin scattering are discussed. These interaction processes are determined by gains, occurrence sites, scattering geometries of each instability, and will affect subsequent evolutions. Nonlinear effects of self-focusing and azimuthal magnetic field generation are observed to be accompanied with SRSS. In addition, it is found that SRSS is insensitive to ion motion, collision (low-Z material), and electron temperature.

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

A high-order particle-in-cell method for low density plasma flow and the simulation of gyrotron resonator devices

International Nuclear Information System (INIS)

Stock, Andreas

2013-01-01

Within this thesis a parallelized, transient, three-dimensional, high-order discontinuous Galerkin Particle-in-Cell solver is developed and used to simulate the resonant cavity of a gyrotron. The high-order discontinuous Galerkin approach - a Finite-Element type method - provides a fast and efficient algorithm to numerically solve Maxwell's equations used within this thesis. Besides its outstanding dissipation and dispersion properties, the discontinuous Galerkin approach easily allows for using unstructured grids, as required to simulate complex-shaped engineering devices. The discontinuous Galerkin approach approximates a wavelength with significantly less degrees of freedom compared to other methods, e.g. Finite Difference methods. Furthermore, the parallelization capabilities of the discontinuous Galerkin framework are excellent due to the very local dependencies between the elements. These properties are essential for the efficient numerical treatment of the Vlasov-Maxwell system with the Particle-in-Cell method. This system describes the self-consistent interaction of charged particles and the electromagnetic field. As central application within this thesis gyrotron resonators are simulated with the discontinuous Galerkin Particle-in-Cell method on high-performance-computers. The gyrotron is a high-power millimeter wave source, used for the electron cyclotron resonance heating of magnetically confined fusion plasma, e.g. in the Wendelstein 7-X experimental fusion-reactor. Compared to state-of-the-art simulation tools used for the design of gyrotron resonators the Particle-in-Cell method does not use any significant physically simplifications w.r.t. the modelling of the particle-field-interaction, the geometry and the wave-spectrum. Hence, it is the method of choice for validation of current simulation tools being restricted by these simplifications. So far, the Particle-in-Cell method was restricted to be used for demonstration calculations only, because

A high-order particle-in-cell method for low density plasma flow and the simulation of gyrotron resonator devices

Energy Technology Data Exchange (ETDEWEB)

Stock, Andreas

2013-04-26

Within this thesis a parallelized, transient, three-dimensional, high-order discontinuous Galerkin Particle-in-Cell solver is developed and used to simulate the resonant cavity of a gyrotron. The high-order discontinuous Galerkin approach - a Finite-Element type method - provides a fast and efficient algorithm to numerically solve Maxwell's equations used within this thesis. Besides its outstanding dissipation and dispersion properties, the discontinuous Galerkin approach easily allows for using unstructured grids, as required to simulate complex-shaped engineering devices. The discontinuous Galerkin approach approximates a wavelength with significantly less degrees of freedom compared to other methods, e.g. Finite Difference methods. Furthermore, the parallelization capabilities of the discontinuous Galerkin framework are excellent due to the very local dependencies between the elements. These properties are essential for the efficient numerical treatment of the Vlasov-Maxwell system with the Particle-in-Cell method. This system describes the self-consistent interaction of charged particles and the electromagnetic field. As central application within this thesis gyrotron resonators are simulated with the discontinuous Galerkin Particle-in-Cell method on high-performance-computers. The gyrotron is a high-power millimeter wave source, used for the electron cyclotron resonance heating of magnetically confined fusion plasma, e.g. in the Wendelstein 7-X experimental fusion-reactor. Compared to state-of-the-art simulation tools used for the design of gyrotron resonators the Particle-in-Cell method does not use any significant physically simplifications w.r.t. the modelling of the particle-field-interaction, the geometry and the wave-spectrum. Hence, it is the method of choice for validation of current simulation tools being restricted by these simplifications. So far, the Particle-in-Cell method was restricted to be used for demonstration calculations only, because

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

Optimized Loading for Particle-in-cell Gyrokinetic Simulations

International Nuclear Information System (INIS)

Lewandowski, J.L.V.

2004-01-01

The problem of particle loading in particle-in-cell gyrokinetic simulations is addressed using a quadratic optimization algorithm. Optimized loading in configuration space dramatically reduces the short wavelength modes in the electrostatic potential that are partly responsible for the non-conservation of total energy; further, the long wavelength modes are resolved with good accuracy. As a result, the conservation of energy for the optimized loading is much better that the conservation of energy for the random loading. The method is valid for any geometry and can be coupled to optimization algorithms in velocity space

Advanced particle-in-cell simulation techniques for modeling the Lockheed Martin Compact Fusion Reactor

Science.gov (United States)

Welch, Dale; Font, Gabriel; Mitchell, Robert; Rose, David

2017-10-01

We report on particle-in-cell developments of the study of the Compact Fusion Reactor. Millisecond, two and three-dimensional simulations (cubic meter volume) of confinement and neutral beam heating of the magnetic confinement device requires accurate representation of the complex orbits, near perfect energy conservation, and significant computational power. In order to determine initial plasma fill and neutral beam heating, these simulations include ionization, elastic and charge exchange hydrogen reactions. To this end, we are pursuing fast electromagnetic kinetic modeling algorithms including a two implicit techniques and a hybrid quasi-neutral algorithm with kinetic ions. The kinetic modeling includes use of the Poisson-corrected direct implicit, magnetic implicit, as well as second-order cloud-in-cell techniques. The hybrid algorithm, ignoring electron inertial effects, is two orders of magnitude faster than kinetic but not as accurate with respect to confinement. The advantages and disadvantages of these techniques will be presented. Funded by Lockheed Martin.

Sparse grid techniques for particle-in-cell schemes

Science.gov (United States)

Ricketson, L. F.; Cerfon, A. J.

2017-02-01

We propose the use of sparse grids to accelerate particle-in-cell (PIC) schemes. By using the so-called ‘combination technique’ from the sparse grids literature, we are able to dramatically increase the size of the spatial cells in multi-dimensional PIC schemes while paying only a slight penalty in grid-based error. The resulting increase in cell size allows us to reduce the statistical noise in the simulation without increasing total particle number. We present initial proof-of-principle results from test cases in two and three dimensions that demonstrate the new scheme’s efficiency, both in terms of computation time and memory usage.

Two-way coupling of magnetohydrodynamic simulations with embedded particle-in-cell simulations

Science.gov (United States)

Makwana, K. D.; Keppens, R.; Lapenta, G.

2017-12-01

We describe a method for coupling an embedded domain in a magnetohydrodynamic (MHD) simulation with a particle-in-cell (PIC) method. In this two-way coupling we follow the work of Daldorff et al. (2014) [19] in which the PIC domain receives its initial and boundary conditions from MHD variables (MHD to PIC coupling) while the MHD simulation is updated based on the PIC variables (PIC to MHD coupling). This method can be useful for simulating large plasma systems, where kinetic effects captured by particle-in-cell simulations are localized but affect global dynamics. We describe the numerical implementation of this coupling, its time-stepping algorithm, and its parallelization strategy, emphasizing the novel aspects of it. We test the stability and energy/momentum conservation of this method by simulating a steady-state plasma. We test the dynamics of this coupling by propagating plasma waves through the embedded PIC domain. Coupling with MHD shows satisfactory results for the fast magnetosonic wave, but significant distortion for the circularly polarized Alfvén wave. Coupling with Hall-MHD shows excellent coupling for the whistler wave. We also apply this methodology to simulate a Geospace Environmental Modeling (GEM) challenge type of reconnection with the diffusion region simulated by PIC coupled to larger scales with MHD and Hall-MHD. In both these cases we see the expected signatures of kinetic reconnection in the PIC domain, implying that this method can be used for reconnection studies.

Species Entropies in the Kinetic Range of Collisionless Plasma Turbulence: Particle-in-cell Simulations

Science.gov (United States)

Gary, S. Peter; Zhao, Yinjian; Hughes, R. Scott; Wang, Joseph; Parashar, Tulasi N.

2018-06-01

Three-dimensional particle-in-cell simulations of the forward cascade of decaying turbulence in the relatively short-wavelength kinetic range have been carried out as initial-value problems on collisionless, homogeneous, magnetized electron-ion plasma models. The simulations have addressed both whistler turbulence at β i = β e = 0.25 and kinetic Alfvén turbulence at β i = β e = 0.50, computing the species energy dissipation rates as well as the increase of the Boltzmann entropies for both ions and electrons as functions of the initial dimensionless fluctuating magnetic field energy density ε o in the range 0 ≤ ε o ≤ 0.50. This study shows that electron and ion entropies display similar rates of increase and that all four entropy rates increase approximately as ε o , consistent with the assumption that the quasilinear premise is valid for the initial conditions assumed for these simulations. The simulations further predict that the time rates of ion entropy increase should be substantially greater for kinetic Alfvén turbulence than for whistler turbulence.

Particle-in-cell simulations of plasma accelerators and electron-neutral collisions

Directory of Open Access Journals (Sweden)

David L. Bruhwiler

2001-10-01

Full Text Available We present 2D simulations of both beam-driven and laser-driven plasma wakefield accelerators, using the object-oriented particle-in-cell code XOOPIC, which is time explicit, fully electromagnetic, and capable of running on massively parallel supercomputers. Simulations of laser-driven wakefields with low \\(∼10^{16} W/cm^{2}\\ and high \\(∼10^{18} W/cm^{2}\\ peak intensity laser pulses are conducted in slab geometry, showing agreement with theory and fluid simulations. Simulations of the E-157 beam wakefield experiment at the Stanford Linear Accelerator Center, in which a 30 GeV electron beam passes through 1 m of preionized lithium plasma, are conducted in cylindrical geometry, obtaining good agreement with previous work. We briefly describe some of the more significant modifications to XOOPIC required by this work, and summarize the issues relevant to modeling relativistic electron-neutral collisions in a particle-in-cell code.

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

Ef: Software for Nonrelativistic Beam Simulation by Particle-in-Cell Algorithm

Science.gov (United States)

Boytsov, A. Yu.; Bulychev, A. A.

2018-04-01

Understanding of particle dynamics is crucial in construction of electron guns, ion sources and other types of nonrelativistic beam devices. Apart from external guiding and focusing systems, a prominent role in evolution of such low-energy beams is played by particle-particle interaction. Numerical simulations taking into account these effects are typically accomplished by a well-known particle-in-cell method. In practice, for convenient work a simulation program should not only implement this method, but also support parallelization, provide integration with CAD systems and allow access to details of the simulation algorithm. To address the formulated requirements, development of a new open source code - Ef - has been started. It's current features and main functionality are presented. Comparison with several analytical models demonstrates good agreement between the numerical results and the theory. Further development plans are discussed.

Ef: Software for Nonrelativistic Beam Simulation by Particle-in-Cell Algorithm

Directory of Open Access Journals (Sweden)

Boytsov A. Yu.

2018-01-01

Full Text Available Understanding of particle dynamics is crucial in construction of electron guns, ion sources and other types of nonrelativistic beam devices. Apart from external guiding and focusing systems, a prominent role in evolution of such low-energy beams is played by particle-particle interaction. Numerical simulations taking into account these effects are typically accomplished by a well-known particle-in-cell method. In practice, for convenient work a simulation program should not only implement this method, but also support parallelization, provide integration with CAD systems and allow access to details of the simulation algorithm. To address the formulated requirements, development of a new open source code - Ef - has been started. It's current features and main functionality are presented. Comparison with several analytical models demonstrates good agreement between the numerical results and the theory. Further development plans are discussed.

Particle-in-cell simulations on spontaneous thermal magnetic field fluctuations

Energy Technology Data Exchange (ETDEWEB)

Simões, F. J. R. Jr.; Pavan, J. [Instituto de Física e Matemática, UFPel, Pelotas, RS (Brazil); Gaelzer, R.; Ziebell, L. F. [Instituto de Física, UFRGS, Porto Alegre, RS (Brazil); Yoon, P. H. [Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742 (United States)

2013-10-15

In this paper an electromagnetic particle code is used to investigate the spontaneous thermal emission. Specifically we perform particle-in-cell simulations employing a non-relativistic isotropic Maxwellian particle distribution to show that thermal fluctuations are related to the origin of spontaneous magnetic field fluctuation. These thermal fluctuations can become seed for further amplification mechanisms and thus be considered at the origin of the cosmological magnetic field, at microgauss levels. Our numerical results are in accordance with theoretical results presented in the literature.

Particle-in-Cell laser-plasma simulation on Xeon Phi coprocessors

Science.gov (United States)

Surmin, I. A.; Bastrakov, S. I.; Efimenko, E. S.; Gonoskov, A. A.; Korzhimanov, A. V.; Meyerov, I. B.

2016-05-01

This paper concerns the development of a high-performance implementation of the Particle-in-Cell method for plasma simulation on Intel Xeon Phi coprocessors. We discuss the suitability of the method for Xeon Phi architecture and present our experience in the porting and optimization of the existing parallel Particle-in-Cell code PICADOR. Direct porting without code modification gives performance on Xeon Phi close to that of an 8-core CPU on a benchmark problem with 50 particles per cell. We demonstrate step-by-step optimization techniques, such as improving data locality, enhancing parallelization efficiency and vectorization leading to an overall 4.2 × speedup on CPU and 7.5 × on Xeon Phi compared to the baseline version. The optimized version achieves 16.9 ns per particle update on an Intel Xeon E5-2660 CPU and 9.3 ns per particle update on an Intel Xeon Phi 5110P. For a real problem of laser ion acceleration in targets with surface grating, where a large number of macroparticles per cell is required, the speedup of Xeon Phi compared to CPU is 1.6 ×.

Advective isotope transport by mixing cell and particle tracking algorithms

International Nuclear Information System (INIS)

Tezcan, L.; Meric, T.

1999-01-01

The 'mixing cell' algorithm of the environmental isotope data evaluation is integrated with the three dimensional finite difference ground water flow model (MODFLOW) to simulate the advective isotope transport and the approach is compared with the 'particle tracking' algorithm of the MOC3D, that simulates three-dimensional solute transport with the method of characteristics technique

A Particle-in-Cell Simulation for the Traveling Wave Direct Energy Converter (TWDEC) for Fusion Propulsion

Science.gov (United States)

Chap, Andrew; Tarditi, Alfonso G.; Scott, John H.

2013-01-01

A Particle-in-cell simulation model has been developed to study the physics of the Traveling Wave Direct Energy Converter (TWDEC) applied to the conversion of charged fusion products into electricity. In this model the availability of a beam of collimated fusion products is assumed; the simulation is focused on the conversion of the beam kinetic energy into alternating current (AC) electric power. The model is electrostatic, as the electro-dynamics of the relatively slow ions can be treated in the quasistatic approximation. A two-dimensional, axisymmetric (radial-axial coordinates) geometry is considered. Ion beam particles are injected on one end and travel along the axis through ring-shaped electrodes with externally applied time-varying voltages, thus modulating the beam by forming a sinusoidal pattern in the beam density. Further downstream, the modulated beam passes through another set of ring electrodes, now electrically oating. The modulated beam induces a time alternating potential di erence between adjacent electrodes. Power can be drawn from the electrodes by connecting a resistive load. As energy is dissipated in the load, a corresponding drop in beam energy is measured. The simulation encapsulates the TWDEC process by reproducing the time-dependent transfer of energy and the particle deceleration due to the electric eld phase time variations.

Dissipative particle dynamics simulations of deformation and aggregation of healthy and diseased red blood cells in a tube flow

International Nuclear Information System (INIS)

Ye, Ting; Phan-Thien, Nhan; Khoo, Boo Cheong; Lim, Chwee Teck

2014-01-01

In this paper, we report simulation results assessing the deformation and aggregation of mixed healthy and malaria-infected red blood cells (RBCs) in a tube flow. A three dimensional particle model based on Dissipative Particle Dynamics (DPD) is developed to predict the tube flow containing interacting cells. The cells are also modelled by DPD, with a Morse potential to characterize the cell-cell interaction. As validation tests, a single RBC in a tube flow and two RBCs in a static flow are simulated to examine the cell deformation and intercellular interaction, respectively. The study of two cells, one healthy and the other malaria-infected RBCs in a tube flow demonstrates that the malaria-infected RBC (in the leading position along flow direction) has different effects on the healthy RBC (in the trailing position) at the different stage of parasite development or at the different capillary number. With parasitic development, the malaria-infected RBC gradually loses its deformability, and in turn the corresponding trailing healthy RBC also deforms less due to the intercellular interaction. With increasing capillary number, both the healthy and malaria-infected RBCs are likely to undergo an axisymmetric motion. The minimum intercellular distance becomes small enough so that rouleaux is easily formed, i.e., the healthy and malaria-infected RBCs are difficultly disaggregated

Dissipative particle dynamics simulations of deformation and aggregation of healthy and diseased red blood cells in a tube flow

Energy Technology Data Exchange (ETDEWEB)

Ye, Ting; Phan-Thien, Nhan, E-mail: Nhan@nus.edu.sg; Khoo, Boo Cheong; Lim, Chwee Teck [Department of Mechanical Engineering, National University of Singapore, Singapore 119260 (Singapore)

2014-11-15

In this paper, we report simulation results assessing the deformation and aggregation of mixed healthy and malaria-infected red blood cells (RBCs) in a tube flow. A three dimensional particle model based on Dissipative Particle Dynamics (DPD) is developed to predict the tube flow containing interacting cells. The cells are also modelled by DPD, with a Morse potential to characterize the cell-cell interaction. As validation tests, a single RBC in a tube flow and two RBCs in a static flow are simulated to examine the cell deformation and intercellular interaction, respectively. The study of two cells, one healthy and the other malaria-infected RBCs in a tube flow demonstrates that the malaria-infected RBC (in the leading position along flow direction) has different effects on the healthy RBC (in the trailing position) at the different stage of parasite development or at the different capillary number. With parasitic development, the malaria-infected RBC gradually loses its deformability, and in turn the corresponding trailing healthy RBC also deforms less due to the intercellular interaction. With increasing capillary number, both the healthy and malaria-infected RBCs are likely to undergo an axisymmetric motion. The minimum intercellular distance becomes small enough so that rouleaux is easily formed, i.e., the healthy and malaria-infected RBCs are difficultly disaggregated.

Numerical Simulation of a Tumor Growth Dynamics Model Using Particle Swarm Optimization.

Science.gov (United States)

Wang, Zhijun; Wang, Qing

Tumor cell growth models involve high-dimensional parameter spaces that require computationally tractable methods to solve. To address a proposed tumor growth dynamics mathematical model, an instance of the particle swarm optimization method was implemented to speed up the search process in the multi-dimensional parameter space to find optimal parameter values that fit experimental data from mice cancel cells. The fitness function, which measures the difference between calculated results and experimental data, was minimized in the numerical simulation process. The results and search efficiency of the particle swarm optimization method were compared to those from other evolutional methods such as genetic algorithms.

Electron Debye scale Kelvin-Helmholtz instability: Electrostatic particle-in-cell simulations

International Nuclear Information System (INIS)

Lee, Sang-Yun; Lee, Ensang; Kim, Khan-Hyuk; Lee, Dong-Hun; Seon, Jongho; Jin, Ho

2015-01-01

In this paper, we investigated the electron Debye scale Kelvin-Helmholtz (KH) instability using two-dimensional electrostatic particle-in-cell simulations. We introduced a velocity shear layer with a thickness comparable to the electron Debye length and examined the generation of the KH instability. The KH instability occurs in a similar manner as observed in the KH instabilities in fluid or ion scales producing surface waves and rolled-up vortices. The strength and growth rate of the electron Debye scale KH instability is affected by the structure of the velocity shear layer. The strength depends on the magnitude of the velocity and the growth rate on the velocity gradient of the shear layer. However, the development of the electron Debye scale KH instability is mainly determined by the electric field generated by charge separation. Significant mixing of electrons occurs across the shear layer, and a fraction of electrons can penetrate deeply into the opposite side fairly far from the vortices across the shear layer

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.

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.

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.

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

Virtual reality in urban water management: communicating urban flooding with particle-based CFD simulations.

Science.gov (United States)

Winkler, Daniel; Zischg, Jonatan; Rauch, Wolfgang

2018-01-01

For communicating urban flood risk to authorities and the public, a realistic three-dimensional visual display is frequently more suitable than detailed flood maps. Virtual reality could also serve to plan short-term flooding interventions. We introduce here an alternative approach for simulating three-dimensional flooding dynamics in large- and small-scale urban scenes by reaching out to computer graphics. This approach, denoted 'particle in cell', is a particle-based CFD method that is used to predict physically plausible results instead of accurate flow dynamics. We exemplify the approach for the real flooding event in July 2016 in Innsbruck.

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.

Turbulence in Three Dimensional Simulations of Magnetopause Reconnection

Science.gov (United States)

Drake, J. F.; Price, L.; Swisdak, M.; Burch, J. L.; Cassak, P.; Dahlin, J. T.; Ergun, R.

2017-12-01

We present two- and three-dimensional particle-in-cell simulations of the 16 October 2015 MMS magnetopause reconnection event. While the two-dimensional simulation is laminar, turbulence develops at both the x-line and along the magnetic separatrices in the three-dimensional simulation. This turbulence is electromagnetic in nature, is characterized by a wavevector k given by kρ e ˜(m_e/m_i)0.25 with ρ e the electron Larmor radius, and appears to have the ion pressure gradient as its source of free energy. Taken together, these results suggest the instability is a variant of the lower-hybrid drift instability. The turbulence produces electric field fluctuations in the out-of-plane direction (the direction of the reconnection electric field) with an amplitude of around ± 10 mV/m, which is much greater than the reconnection electric field of around 0.1 mV/m. Such large values of the out-of-plane electric field have been identified in the MMS data. The turbulence in the simulation controls the scale lengths of the density profile and current layers in asymmetric reconnection, driving them closer to √ {ρ eρ_i } than the ρ e or de scalings seen in 2D reconnection simulations, where de is the electron inertial length. The turbulence is strong enough to make the magnetic field around the reconnection island chaotic and produces both anomalous resistivity and anomalous viscosity. Each contribute significantly to breaking the frozen-in condition in the electron diffusion region. The crescent-shaped features in velocity space seen both in MMS observations and in two-dimensional simulations survive, even in the turbulent environment of the three-dimensional system. We compare and contrast these results to a three-dimensional simulation of the 8 December 2015 MMS magnetopause reconnection event in which the reconnecting and out-of-plane guide fields are comparable. LHDI is still present in this event, although its appearance is modified by the presence of the guide

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

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.))

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

A comprehensive study of MPI parallelism in three-dimensional discrete element method (DEM) simulation of complex-shaped granular particles

Science.gov (United States)

Yan, Beichuan; Regueiro, Richard A.

2018-02-01

A three-dimensional (3D) DEM code for simulating complex-shaped granular particles is parallelized using message-passing interface (MPI). The concepts of link-block, ghost/border layer, and migration layer are put forward for design of the parallel algorithm, and theoretical scalability function of 3-D DEM scalability and memory usage is derived. Many performance-critical implementation details are managed optimally to achieve high performance and scalability, such as: minimizing communication overhead, maintaining dynamic load balance, handling particle migrations across block borders, transmitting C++ dynamic objects of particles between MPI processes efficiently, eliminating redundant contact information between adjacent MPI processes. The code executes on multiple US Department of Defense (DoD) supercomputers and tests up to 2048 compute nodes for simulating 10 million three-axis ellipsoidal particles. Performance analyses of the code including speedup, efficiency, scalability, and granularity across five orders of magnitude of simulation scale (number of particles) are provided, and they demonstrate high speedup and excellent scalability. It is also discovered that communication time is a decreasing function of the number of compute nodes in strong scaling measurements. The code's capability of simulating a large number of complex-shaped particles on modern supercomputers will be of value in both laboratory studies on micromechanical properties of granular materials and many realistic engineering applications involving granular materials.

Electromagnetic particle-in-cell simulations of Applied-B proton diodes

International Nuclear Information System (INIS)

Slutz, S.A.; Seidel, D.B.; Coats, R.S.

1986-01-01

Fully electromagnetic particle-in-cell simulations of Applied-B ion diodes have been performed using the magic code. These calculations indicate that Applied-B diodes can be nearly 100% efficient. Furthermore, the simulations exhibit an impedance relaxation phenomenon due to the buildup of electron space charge near the anode which causes a time-dependent enhancement of the ion emission above the Child--Langmuir value. This phenomenon may at least partially explain the rapidly decreasing impedance that has been observed in Applied-B ion diode experiments. The results of our numerical simulations will be compared to experimental data on Applied-B ion diodes and to analytic theories of their operation

Vectorization of a particle code used in the simulation of rarefied hypersonic flow

Science.gov (United States)

Baganoff, D.

1990-01-01

A limitation of the direct simulation Monte Carlo (DSMC) method is that it does not allow efficient use of vector architectures that predominate in current supercomputers. Consequently, the problems that can be handled are limited to those of one- and two-dimensional flows. This work focuses on a reformulation of the DSMC method with the objective of designing a procedure that is optimized to the vector architectures found on machines such as the Cray-2. In addition, it focuses on finding a better balance between algorithmic complexity and the total number of particles employed in a simulation so that the overall performance of a particle simulation scheme can be greatly improved. Simulations of the flow about a 3D blunt body are performed with 10 to the 7th particles and 4 x 10 to the 5th mesh cells. Good statistics are obtained with time averaging over 800 time steps using 4.5 h of Cray-2 single-processor CPU time.

Computer-assisted Particle-in-Cell code development

International Nuclear Information System (INIS)

Kawata, S.; Boonmee, C.; Teramoto, T.; Drska, L.; Limpouch, J.; Liska, R.; Sinor, M.

1997-12-01

This report presents a new approach for an electromagnetic Particle-in-Cell (PIC) code development by a computer: in general PIC codes have a common structure, and consist of a particle pusher, a field solver, charge and current density collections, and a field interpolation. Because of the common feature, the main part of the PIC code can be mechanically developed on a computer. In this report we use the packages FIDE and GENTRAN of the REDUCE computer algebra system for discretizations of field equations and a particle equation, and for an automatic generation of Fortran codes. The approach proposed is successfully applied to the development of 1.5-dimensional PIC code. By using the generated PIC code the Weibel instability in a plasma is simulated. The obtained growth rate agrees well with the theoretical value. (author)

Framework to model neutral particle flux in convex high aspect ratio structures using one-dimensional radiosity

Science.gov (United States)

Manstetten, Paul; Filipovic, Lado; Hössinger, Andreas; Weinbub, Josef; Selberherr, Siegfried

2017-02-01

We present a computationally efficient framework to compute the neutral flux in high aspect ratio structures during three-dimensional plasma etching simulations. The framework is based on a one-dimensional radiosity approach and is applicable to simulations of convex rotationally symmetric holes and convex symmetric trenches with a constant cross-section. The framework is intended to replace the full three-dimensional simulation step required to calculate the neutral flux during plasma etching simulations. Especially for high aspect ratio structures, the computational effort, required to perform the full three-dimensional simulation of the neutral flux at the desired spatial resolution, conflicts with practical simulation time constraints. Our results are in agreement with those obtained by three-dimensional Monte Carlo based ray tracing simulations for various aspect ratios and convex geometries. With this framework we present a comprehensive analysis of the influence of the geometrical properties of high aspect ratio structures as well as of the particle sticking probability on the neutral particle flux.

Interactive methods for exploring particle simulation data

Energy Technology Data Exchange (ETDEWEB)

Co, Christopher S.; Friedman, Alex; Grote, David P.; Vay, Jean-Luc; Bethel, E. Wes; Joy, Kenneth I.

2004-05-01

In this work, we visualize high-dimensional particle simulation data using a suite of scatter plot-based visualizations coupled with interactive selection tools. We use traditional 2D and 3D projection scatter plots as well as a novel oriented disk rendering style to convey various information about the data. Interactive selection tools allow physicists to manually classify ''interesting'' sets of particles that are highlighted across multiple, linked views of the data. The power of our application is the ability to correspond new visual representations of the simulation data with traditional, well understood visualizations. This approach supports the interactive exploration of the high-dimensional space while promoting discovery of new particle behavior.

Fully implicit Particle-in-cell algorithms for multiscale plasma simulation

Energy Technology Data Exchange (ETDEWEB)

Chacon, Luis [Los Alamos National Laboratory

2015-07-16

The outline of the paper is as follows: Particle-in-cell (PIC) methods for fully ionized collisionless plasmas, explicit vs. implicit PIC, 1D ES implicit PIC (charge and energy conservation, moment-based acceleration), and generalization to Multi-D EM PIC: Vlasov-Darwin model (review and motivation for Darwin model, conservation properties (energy, charge, and canonical momenta), and numerical benchmarks). The author demonstrates a fully implicit, fully nonlinear, multidimensional PIC formulation that features exact local charge conservation (via a novel particle mover strategy), exact global energy conservation (no particle self-heating or self-cooling), adaptive particle orbit integrator to control errors in momentum conservation, and canonical momenta (EM-PIC only, reduced dimensionality). The approach is free of numerical instabilities: ω_peΔt >> 1, and Δx >> λ_D. It requires many fewer dofs (vs. explicit PIC) for comparable accuracy in challenging problems. Significant CPU gains (vs explicit PIC) have been demonstrated. The method has much potential for efficiency gains vs. explicit in long-time-scale applications. Moment-based acceleration is effective in minimizing N_FE, leading to an optimal algorithm.

Blood–plasma separation in Y-shaped bifurcating microfluidic channels: a dissipative particle dynamics simulation study

International Nuclear Information System (INIS)

Li, Xuejin; Karniadakis, George Em; Popel, Aleksander S

2012-01-01

The motion of a suspension of red blood cells (RBCs) flowing in a Y-shaped bifurcating microfluidic channel is investigated using a validated low-dimensional RBC model based on dissipative particle dynamics. Specifically, the RBC is represented as a closed torus-like ring of ten colloidal particles, which leads to efficient simulations of blood flow in microcirculation over a wide range of hematocrits. Adaptive no-slip wall boundary conditions were implemented to model hydrodynamic flow within a specific wall structure of diverging three-dimensional microfluidic channels, paying attention to controlling density fluctuations. Plasma skimming and the all-or-nothing phenomenon of RBCs in a bifurcating microfluidic channel have been investigated in our simulations for healthy and diseased blood, including the size of a cell-free layer on the daughter branches. The feed hematocrit level in the parent channel has considerable influence on blood–plasma separation. Compared to the blood–plasma separation efficiencies of healthy RBCs, malaria-infected stiff RBCs (iRBCs) have a tendency to travel into the low flow-rate daughter branch because of their different initial distribution in the parent channel. Our simulation results are consistent with previously published experimental results and theoretical predictions. (paper)

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)

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.

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.

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 simulation and auto-stereoscopic 3D display of the battlefield environment based on the particle system algorithm

Science.gov (United States)

Ning, Jiwei; Sang, Xinzhu; Xing, Shujun; Cui, Huilong; Yan, Binbin; Yu, Chongxiu; Dou, Wenhua; Xiao, Liquan

2016-10-01

The army's combat training is very important now, and the simulation of the real battlefield environment is of great significance. Two-dimensional information has been unable to meet the demand at present. With the development of virtual reality technology, three-dimensional (3D) simulation of the battlefield environment is possible. In the simulation of 3D battlefield environment, in addition to the terrain, combat personnel and the combat tool ,the simulation of explosions, fire, smoke and other effects is also very important, since these effects can enhance senses of realism and immersion of the 3D scene. However, these special effects are irregular objects, which make it difficult to simulate with the general geometry. Therefore, the simulation of irregular objects is always a hot and difficult research topic in computer graphics. Here, the particle system algorithm is used for simulating irregular objects. We design the simulation of the explosion, fire, smoke based on the particle system and applied it to the battlefield 3D scene. Besides, the battlefield 3D scene simulation with the glasses-free 3D display is carried out with an algorithm based on GPU 4K super-multiview 3D video real-time transformation method. At the same time, with the human-computer interaction function, we ultimately realized glasses-free 3D display of the simulated more realistic and immersed 3D battlefield environment.

Three-dimensional simulation of super-radiant Smith-Purcell radiation

International Nuclear Information System (INIS)

Li, D.; Imasaki, K.; Yang, Z.; Park, Gun-Sik

2006-01-01

A simulation of coherent and super-radiant Smith-Purcell radiation is performed in the gigahertz regime using a three-dimensional particle-in-cell code. The simulation model supposes a rectangular grating to be driven by a single electron bunch and a train of periodic bunches, respectively. The true Smith-Purcell radiation is distinguished from the evanescent wave, which has an angle independent frequency lower than the minimum allowed Smith-Purcell frequency. We also find that the super-radiant radiations excited by periodic bunches are emitted at higher harmonics of the bunching frequency and at the corresponding Smith-Purcell angles

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.

Plasma physics via particle simulation

International Nuclear Information System (INIS)

Birdsall, C.K.

1981-01-01

Plasmas are studied by following the motion of many particles in applied and self fields, analytically, experimentally and computationally. Plasmas for magnetic fusion energy devices are very hot, nearly collisionless and magnetized, with scale lengths of many ion gyroradii and Debye lengths. The analytic studies of such plasmas are very difficult as the plasma is nonuniform, anisotropic and nonlinear. The experimental studies have become very expensive in time and money, as the size, density and temperature approach fusion reactor values. Computational studies using many particles and/or fluids have complemented both theories and experiments for many years and have progressed to fully three dimensional electromagnetic models, albeit with hours of running times on the fastest largest computers. Particle simulation methods are presented in some detail, showing particle advance from acceleration to velocity to position, followed by calculation of the fields from charge and current densities and then further particle advance, and so on. Limitations due to the time stepping and use of a spatial grid are given, to avoid inaccuracies and instabilities. Examples are given for an electrostatic program in one dimension of an orbit averaging program, and for a three dimensional electromagnetic program. Applications of particle simulations of plasmas in magnetic and inertial fusion devices continue to grow, as well as to plasmas and beams in peripheral devices, such as sources, accelerators, and converters. (orig.)

Extended particle-based simulation for magnetic-aligned compaction of hard magnetic particles

Energy Technology Data Exchange (ETDEWEB)

Soda, Rikio; Takagi, Kenta; Ozaki, Kimihiro, E-mail: r-soda@aist.go.jp

2015-12-15

In order to understand the magnetic-aligned compaction process, we develop a three-dimensional (3D) discrete element method for simulating the motion of hard magnetic particles subjected to strong compression and magnetic fields. The proposed simulation model also considers the exact magnetic force involved via the calculation of the magnetic moment. First, to validate the simulation model, single-action compaction in the absence of a magnetic field was calculated. The calculated compaction curves are in good quantitative agreement with experimental ones. Based on this simulation model, the alignment behavior of Nd–Fe–B particles during compression under the application of a static magnetic field. The developed simulation model enables the visualization of particle behavior including the misorientation of the magnetization easy axis, which provided the quantitative relationships between applied pressure and particle misorientation. - Highlights: • A practical 3D DEM simulation technique for magnetic-aligned compaction was developed. • An extended simulation model was introduced for hard magnetic particles. • Magnetic-aligned compaction was simulated using the developed simulation model.

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 of magnetically driven reconnection using laser-powered capacitor coils

Science.gov (United States)

Huang, Kai; Lu, Quanming; Gao, Lan; Ji, Hantao; Wang, Xueyi; Fan, Feibin

2018-05-01

In this paper, we propose an experimental scheme to fulfill magnetically driven reconnections. Here, two laser beams are focused on a capacitor-coil target and then strong currents are wired in two parallel circular coils. Magnetic reconnection occurs between the two magnetic bubbles created by the currents in the two parallel circular coils. A two-dimensional particle-in-cell simulation model in the cylindrical coordinate is used to investigate such a process, and the simulations are performed in the (r ,z ) plane. The results show that with the increase of the currents in the two coils, the associated magnetic bubbles expand and a current sheet is formed between the two bubbles. Magnetic reconnection occurs when the current sheet is sufficiently thin. A quadrupole structure of the magnetic field in the θ direction ( Bθ ) is generated in the diffusion region and a strong electron current along the r direction ( Je r ) is also formed due to the existence of the high-speed electron flow away from the X line in the center of the outflow region. Because the X line is a circle along the θ direction, the convergence of the plasma flow around r =0 will lead to the asymmetry of Je r and Bθ between the two outflow regions of magnetic reconnection.

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.

Lagrangian Particle Tracking Simulation for Warm-Rain Processes in Quasi-One-Dimensional Domain

Science.gov (United States)

Kunishima, Y.; Onishi, R.

2017-12-01

Conventional cloud simulations are based on the Euler method and compute each microphysics process in a stochastic way assuming infinite numbers of particles within each numerical grid. They therefore cannot provide the Lagrangian statistics of individual particles in cloud microphysics (i.e., aerosol particles, cloud particles, and rain drops) nor discuss the statistical fluctuations due to finite number of particles. We here simulate the entire precipitation process of warm-rain, with tracking individual particles. We use the Lagrangian Cloud Simulator (LCS), which is based on the Euler-Lagrangian framework. In that framework, flow motion and scalar transportation are computed with the Euler method, and particle motion with the Lagrangian one. The LCS tracks particle motions and collision events individually with considering the hydrodynamic interaction between approaching particles with a superposition method, that is, it can directly represent the collisional growth of cloud particles. It is essential for trustworthy collision detection to take account of the hydrodynamic interaction. In this study, we newly developed a stochastic model based on the Twomey cloud condensation nuclei (CCN) activation for the Lagrangian tracking simulation and integrated it into the LCS. Coupling with the Euler computation for water vapour and temperature fields, the initiation and condensational growth of water droplets were computed in the Lagrangian way. We applied the integrated LCS for a kinematic simulation of warm-rain processes in a vertically-elongated domain of, at largest, 0.03×0.03×3000 (m3) with horizontal periodicity. Aerosol particles with a realistic number density, 5×107 (m3), were evenly distributed over the domain at the initial state. Prescribed updraft at the early stage initiated development of a precipitating cloud. We have confirmed that the obtained bulk statistics fairly agree with those from a conventional spectral-bin scheme for a vertical column

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)

Quasi-three-dimensional particle imaging with digital holography.

Science.gov (United States)

Kemppinen, Osku; Heinson, Yuli; Berg, Matthew

2017-05-01

In this work, approximate three-dimensional structures of microparticles are generated with digital holography using an automated focus method. This is done by stacking a collection of silhouette-like images of a particle reconstructed from a single in-line hologram. The method enables estimation of the particle size in the longitudinal and transverse dimensions. Using the discrete dipole approximation, the method is tested computationally by simulating holograms for a variety of particles and attempting to reconstruct the known three-dimensional structure. It is found that poor longitudinal resolution strongly perturbs the reconstructed structure, yet the method does provide an approximate sense for the structure's longitudinal dimension. The method is then applied to laboratory measurements of holograms of single microparticles and their scattering patterns.

Lattice Boltzmann method used to simulate particle motion in a conduit

Directory of Open Access Journals (Sweden)

Dolanský Jindřich

2017-06-01

Full Text Available A three-dimensional numerical simulation of particle motion in a pipe with a rough bed is presented. The simulation based on the Lattice Boltzmann Method (LBM employs the hybrid diffuse bounce-back approach to model moving boundaries. The bed of the pipe is formed by stationary spherical particles of the same size as the moving particles. Particle movements are induced by gravitational and hydrodynamic forces. To evaluate the hydrodynamic forces, the Momentum Exchange Algorithm is used. The LBM unified computational frame makes it possible to simulate both the particle motion and the fluid flow and to study mutual interactions of the carrier liquid flow and particles and the particle–bed and particle–particle collisions. The trajectories of simulated and experimental particles are compared. The Particle Tracking method is used to track particle motion. The correctness of the applied approach is assessed.

A particle-based model to simulate the micromechanics of single-plant parenchyma cells and aggregates

International Nuclear Information System (INIS)

Van Liedekerke, P; Tijskens, E; Smeedts, B; Ramon, H; Ghysels, P; Samaey, G; Roose, D

2010-01-01

This paper is concerned with addressing how plant tissue mechanics is related to the micromechanics of cells. To this end, we propose a mesh-free particle method to simulate the mechanics of both individual plant cells (parenchyma) and cell aggregates in response to external stresses. The model considers two important features in the plant cell: (1) the cell protoplasm, the interior liquid phase inducing hydrodynamic phenomena, and (2) the cell wall material, a viscoelastic solid material that contains the protoplasm. In this particle framework, the cell fluid is modeled by smoothed particle hydrodynamics (SPH), a mesh-free method typically used to address problems with gas and fluid dynamics. In the solid phase (cell wall) on the other hand, the particles are connected by pairwise interactions holding them together and preventing the fluid to penetrate the cell wall. The cell wall hydraulic conductivity (permeability) is built in as well through the SPH formulation. Although this model is also meant to be able to deal with dynamic and even violent situations (leading to cell wall rupture or cell–cell debonding), we have concentrated on quasi-static conditions. The results of single-cell compression simulations show that the conclusions found by analytical models and experiments can be reproduced at least qualitatively. Relaxation tests revealed that plant cells have short relaxation times (1 µs–10 µs) compared to mammalian cells. Simulations performed on cell aggregates indicated an influence of the cellular organization to the tissue response, as was also observed in experiments done on tissues with a similar structure

Three-dimensional simulation of charge collection and multiple-bit upset in Si devices

International Nuclear Information System (INIS)

Dodd, P.E.; Sexton, F.W.; Winokur, P.S.

1994-01-01

In this paper, three-dimensional numerical simulation is used to explore the basic charge-collection mechanisms in silicon n + /p diodes. For diodes on lightly-doped substrates ( 15 cm -3 ) struck by a 100-MeV Fe ion, the funneling effect is very strong and essentially all collection is by funnel-assisted drift. This drift collection may occur as late as several nanoseconds after the strike, later than is usually associated with drift collection. For moderately-doped substrates (∼1 x 10 16 cm -3 ) and epitaxial structures grown on heavily-doped substrates, the funnel effect is weaker and drift and diffusion are of more equal importance. For 5-MeV He (α-particle) strikes with low-density charge tracks, the charge-collection transient exhibits both drift and diffusion regimes regardless of the substrate doping. Simulations of diodes with passive external loads indicate that while the current response is altered considerably by the load, total collected charge is not greatly affected for the simple resistive loads studied. Three-dimensional mixed-mode simulation is performed to investigate charge-collection behavior and upset mechanisms in complete CMOS SRAM cells. Simulations of double SRAM cell structures indicate that only collection by diffusion from ''between-node'' strikes is capable of producing multiple-bit upsets in the simulated technology. Limitations of the simulations, specifically carrier-carrier scattering models and large concentration gradients, are also discussed

Realistic simulations of a cyclotron spiral inflector within a particle-in-cell framework

Science.gov (United States)

Winklehner, Daniel; Adelmann, Andreas; Gsell, Achim; Kaman, Tulin; Campo, Daniela

2017-12-01

We present an upgrade to the particle-in-cell ion beam simulation code opal that enables us to run highly realistic simulations of the spiral inflector system of a compact cyclotron. This upgrade includes a new geometry class and field solver that can handle the complicated boundary conditions posed by the electrode system in the central region of the cyclotron both in terms of particle termination, and calculation of self-fields. Results are benchmarked against the analytical solution of a coasting beam. As a practical example, the spiral inflector and the first revolution in a 1 MeV /amu test cyclotron, located at Best Cyclotron Systems, Inc., are modeled and compared to the simulation results. We find that opal can now handle arbitrary boundary geometries with relative ease. Simulated injection efficiencies and beam shape compare well with measured efficiencies and a preliminary measurement of the beam distribution after injection.

Particle-in-cell simulations of electron transport from plasma turbulence: recent progress in gyrokinetic particle simulations of turbulent plasmas

International Nuclear Information System (INIS)

Lin, Z; Rewoldt, G; Ethier, S; Hahm, T S; Lee, W W; Lewandowski, J L V; Nishimura, Y; Wang, W X

2005-01-01

Recent progress in gyrokinetic particle-in-cell simulations of turbulent plasmas using the gyrokinetic toroidal code (GTC) is surveyed. In particular, recent results for electron temperature gradient (ETG) modes and their resulting transport are presented. Also, turbulence spreading, and the effects of the parallel nonlinearity, are described. The GTC code has also been generalized for non-circular plasma cross-section, and initial results are presented. In addition, two distinct methods of generalizing the GTC code to be electromagnetic are described, along with preliminary results. Finally, a related code, GTC-Neo, for calculating neoclassical fluxes, electric fields, and velocities, are described

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.

High energy gain in three-dimensional simulations of light sail acceleration

Energy Technology Data Exchange (ETDEWEB)

Sgattoni, A., E-mail: andrea.sgattoni@polimi.it [Dipartimento di Energia, Politecnico di Milano, Milano (Italy); CNR, Istituto Nazionale di Ottica, u.o.s. “Adriano Gozzini,” Pisa (Italy); Sinigardi, S. [CNR, Istituto Nazionale di Ottica, u.o.s. “Adriano Gozzini,” Pisa (Italy); Dipartimento di Fisica e Astronomia, Università di Bologna, Bologna (Italy); INFN sezione di Bologna, Bologna (Italy); Macchi, A. [CNR, Istituto Nazionale di Ottica, u.o.s. “Adriano Gozzini,” Pisa (Italy); Dipartimento di Fisica “Enrico Fermi,” Università di Pisa, Pisa (Italy)

2014-08-25

The dynamics of radiation pressure acceleration in the relativistic light sail regime are analysed by means of large scale, three-dimensional (3D) particle-in-cell simulations. Differently to other mechanisms, the 3D dynamics leads to faster and higher energy gain than in 1D or 2D geometry. This effect is caused by the local decrease of the target density due to transverse expansion leading to a “lighter sail.” However, the rarefaction of the target leads to an earlier transition to transparency limiting the energy gain. A transverse instability leads to a structured and inhomogeneous ion distribution.

High energy gain in three-dimensional simulations of light sail acceleration

International Nuclear Information System (INIS)

Sgattoni, A.; Sinigardi, S.; Macchi, A.

2014-01-01

The dynamics of radiation pressure acceleration in the relativistic light sail regime are analysed by means of large scale, three-dimensional (3D) particle-in-cell simulations. Differently to other mechanisms, the 3D dynamics leads to faster and higher energy gain than in 1D or 2D geometry. This effect is caused by the local decrease of the target density due to transverse expansion leading to a “lighter sail.” However, the rarefaction of the target leads to an earlier transition to transparency limiting the energy gain. A transverse instability leads to a structured and inhomogeneous ion distribution.

Self-organization and oscillation of negatively charged dust particles in a 2-dimensional dusty plasma

Energy Technology Data Exchange (ETDEWEB)

Song, Y.L. [College of Science, China Agricultural University, Beijing 100083 (China); Huang, F., E-mail: huangfeng@cau.edu.cn [College of Science, China Agricultural University, Beijing 100083 (China); Chen, Z.Y., E-mail: chenzy@mail.buct.edu.cn [Department of Physics, Beijing University of Chemical Technology, Beijing 100029 (China); State Key Laboratory of Laser Propulsion & Application, Beijing 101416 (China); Liu, Y.H. [School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025 (China); Yu, M.Y. [Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027 (China); Institute for Theoretical Physics I, Ruhr University, D-44801 Bochum (Germany)

2016-02-22

Negatively charged dust particles immersed in 2-dimensional dusty plasma system are investigated by molecular dynamics simulations. The effects of the confinement potential and attraction interaction potential on dust particle self-organization are studied in detail and two typical dust particle distributions are obtained when the system reaches equilibrium. The average radial velocity (ARV), average radial force (ARF) and radial mean square displacement are employed to analyze the dust particles' dynamics. Both ARVs and ARFs exhibit oscillation behaviors when the simulation system reaches equilibrium state. The relationships between the oscillation and confinement potential and attraction potential are studied in this paper. The simulation results are qualitatively similar to experimental results. - Highlights: • Self-organization and oscillation of a 2-dimensional dusty plasma is investigated. • Effect of the confinement potential on dust self-organization and oscillation is given. • Effect of the attraction potential on dust self-organization and oscillation is studied.

Self-organization and oscillation of negatively charged dust particles in a 2-dimensional dusty plasma

International Nuclear Information System (INIS)

Song, Y.L.; Huang, F.; Chen, Z.Y.; Liu, Y.H.; Yu, M.Y.

2016-01-01

Negatively charged dust particles immersed in 2-dimensional dusty plasma system are investigated by molecular dynamics simulations. The effects of the confinement potential and attraction interaction potential on dust particle self-organization are studied in detail and two typical dust particle distributions are obtained when the system reaches equilibrium. The average radial velocity (ARV), average radial force (ARF) and radial mean square displacement are employed to analyze the dust particles' dynamics. Both ARVs and ARFs exhibit oscillation behaviors when the simulation system reaches equilibrium state. The relationships between the oscillation and confinement potential and attraction potential are studied in this paper. The simulation results are qualitatively similar to experimental results. - Highlights: • Self-organization and oscillation of a 2-dimensional dusty plasma is investigated. • Effect of the confinement potential on dust self-organization and oscillation is given. • Effect of the attraction potential on dust self-organization and oscillation is studied.

Modelling and simulation of particle-particle interaction in a magnetophoretic bio-separation chip

Science.gov (United States)

Alam, Manjurul; Golozar, Matin; Darabi, Jeff

2018-04-01

A Lagrangian particle trajectory model is developed to predict the interaction between cell-bead particle complexes and to track their trajectories in a magnetophoretic bio-separation chip. Magnetic flux gradients are simulated in the OpenFOAM CFD software and imported into MATLAB to obtain the trapping lengths and trajectories of the particles. A connector vector is introduced to calculate the interaction force between cell-bead complexes as they flow through a microfluidic device. The interaction force calculations are performed for cases where the connector vector is parallel, perpendicular, and at an angle of 45° with the applied magnetic field. The trajectories of the particles are simulated by solving a system of eight ordinary differential equations using a fourth order Runge-Kutta method. The model is then used to study the effects of geometric positions and angles of the connector vector between the particles as well as the cell size, number of beads per cell, and flow rate on the interaction force and trajectories of the particles. The results show that the interaction forces may be attractive or repulsive, depending on the orientation of the connector vector distance between the particle complexes and the applied magnetic field. When the interaction force is attractive, the particles are observed to merge and trap sooner than a single particle, whereas a repulsive interaction force has little or no effect on the trapping length.

Particle-in-cell simulations of plasma opening switch with external magnetic field

International Nuclear Information System (INIS)

Chen Yulan; Zeng Zhengzhong; Sun Fengju

2003-01-01

Fully electromagnetic particle-in-cell simulations are performed to study the effects of an external magnetic field on coaxial plasma opening switch (POS). The simulation results show that POS opening performance can be significantly improved only when external longitudinal magnetic field coils are placed at the cathode side, and an additional azimuthal magnetic field is effective whether the central electrode is of positive or negative polarity. Voltage multiplication coefficient K rises with the additional magnetic field increasing till the electron current is completely magnetically insulated during the opening of POS

Particle-in-cell numerical simulations of a cylindrical Hall thruster with permanent magnets

Science.gov (United States)

Miranda, Rodrigo A.; Martins, Alexandre A.; Ferreira, José L.

2017-10-01

The cylindrical Hall thruster (CHT) is a propulsion device that offers high propellant utilization and performance at smaller dimensions and lower power levels than traditional Hall thrusters. In this paper we present first results of a numerical model of a CHT. This model solves particle and field dynamics self-consistently using a particle-in-cell approach. We describe a number of techniques applied to reduce the execution time of the numerical simulations. The specific impulse and thrust computed from our simulations are in agreement with laboratory experiments. This simplified model will allow for a detailed analysis of different thruster operational parameters and obtain an optimal configuration to be implemented at the Plasma Physics Laboratory at the University of Brasília.

Axisymmetric particle-in-cell simulations of diamagnetic-cavity formation in vacuum

International Nuclear Information System (INIS)

Gisler, G.

1989-01-01

Axisymmetric simulations of the expansion of a hot plasma suddenly introduced into a vacuum containing a weak magnetic field were performed using an electromagnetic particle-in-cell code. Both uniform and gradient fields have been used, with the simulation axis along the principle field direction. The formation of a diamagnetic cavity requires an initial plasma β > 1; as the expansion proceeds, β diminishes, and the field eventually recovers. The maximum spatial extent of the cavity and its duration can be obtained from simple dynamical considerations. Field-aligned ion acceleration behind the electron front is observed in all field geometries and strengths. In the case of expansion into a divergent field, the plasma is found to move down the field gradient by ambipolar diffusion. These simulations are relevant to active release experiments in the Earth's magnetosphere, to pellet ablation experiments, and to the naturally occurring diamagnetic bubbles observed at the Earth's foreshock

A parallel electrostatic Particle-in-Cell method on unstructured tetrahedral grids for large-scale bounded collisionless plasma simulations

Science.gov (United States)

Averkin, Sergey N.; Gatsonis, Nikolaos A.

2018-06-01

An unstructured electrostatic Particle-In-Cell (EUPIC) method is developed on arbitrary tetrahedral grids for simulation of plasmas bounded by arbitrary geometries. The electric potential in EUPIC is obtained on cell vertices from a finite volume Multi-Point Flux Approximation of Gauss' law using the indirect dual cell with Dirichlet, Neumann and external circuit boundary conditions. The resulting matrix equation for the nodal potential is solved with a restarted generalized minimal residual method (GMRES) and an ILU(0) preconditioner algorithm, parallelized using a combination of node coloring and level scheduling approaches. The electric field on vertices is obtained using the gradient theorem applied to the indirect dual cell. The algorithms for injection, particle loading, particle motion, and particle tracking are parallelized for unstructured tetrahedral grids. The algorithms for the potential solver, electric field evaluation, loading, scatter-gather algorithms are verified using analytic solutions for test cases subject to Laplace and Poisson equations. Grid sensitivity analysis examines the L2 and L∞ norms of the relative error in potential, field, and charge density as a function of edge-averaged and volume-averaged cell size. Analysis shows second order of convergence for the potential and first order of convergence for the electric field and charge density. Temporal sensitivity analysis is performed and the momentum and energy conservation properties of the particle integrators in EUPIC are examined. The effects of cell size and timestep on heating, slowing-down and the deflection times are quantified. The heating, slowing-down and the deflection times are found to be almost linearly dependent on number of particles per cell. EUPIC simulations of current collection by cylindrical Langmuir probes in collisionless plasmas show good comparison with previous experimentally validated numerical results. These simulations were also used in a parallelization

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

Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell

Science.gov (United States)

Osnes, A. N.; Vartdal, M.; Pettersson Reif, B. A.

2018-05-01

The formation of jets from a shock-accelerated cylindrical shell of particles, confined in a Hele-Shaw cell, is studied by means of numerical simulation. A number of simulations have been performed, systematically varying the coupling between the gas and solid phases in an effort to identify the primary mechanism(s) responsible for jet formation. We find that coupling through drag is sufficient for the formation of jets. Including the effect of particle volume fraction and particle collisions did not alter the general behaviour, but had some influence on the length, spacing and number of jets. Furthermore, we find that the jet selection process starts early in the dispersal process, during the initial expansion of the particle layer.

The dynamics of low-β plasma clouds as simulated by a three-dimensional, electromagnetic particle code

International Nuclear Information System (INIS)

Neubert, T.; Miller, R.H.; Buneman, O.; Nishikawa, K.I.

1992-01-01

The dynamics of low-β plasma clouds moving perpendicular to an ambient magnetic field in vacuum and in a background plasma is simulated by means of a three-dimensional, electromagnetic, and relativistic particle simulation code. The simulations show the formation of the space charge sheaths at the sides of the cloud with the associated polarization electric field which facilitate the cross-field propagation, as well as the sheaths at the front and rear end of the cloud caused by the larger ion Larmor radius, which allows ions to move ahead and lag behind the electrons as they gyrate. Results on the cloud dynamics and electromagnetic radiation include the following: (1) In a background plasma, electron and ion sheaths expand along the magnetic field at the same rate, whereas in vacuum the electron sheath expands much faster than the ion sheath. (2) Sheath electrons are accelerated up to relativistic energies. This result indicates that artificial plasma clouds released in the ionosphere or magnetosphere may generate optical emissions (aurora) as energetic sheath electrons scatter in the upper atmosphere. (3) The expansion of the electron sheaths is analogous to the ejection of high-intensity electron beams from spacecraft. (4) Second-order and higher-order sheaths are formed which extend out into the ambient plasma. (5) Formation of the sheaths and the polarization field reduces the forward momentum of the cloud. (6) The coherent component of the particle gyromotion is damped in time as the particles establish a forward directed drift velocity. (7) The coherent particle gyrations generate electromagnetic radiation

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.

Three-dimensional gyrokinetic particle-in-cell simulation of plasmas on a massively parallel computer: Final report on LDRD Core Competency Project, FY 1991--FY 1993

International Nuclear Information System (INIS)

Byers, J.A.; Williams, T.J.; Cohen, B.I.; Dimits, A.M.

1994-01-01

One of the programs of the Magnetic fusion Energy (MFE) Theory and computations Program is studying the anomalous transport of thermal energy across the field lines in the core of a tokamak. We use the method of gyrokinetic particle-in-cell simulation in this study. For this LDRD project we employed massively parallel processing, new algorithms, and new algorithms, and new formal techniques to improve this research. Specifically, we sought to take steps toward: researching experimentally-relevant parameters in our simulations, learning parallel computing to have as a resource for our group, and achieving a 100 x speedup over our starting-point Cray2 simulation code's performance

Direct Lagrangian tracking simulations of particles in vertically-developing atmospheric clouds

Science.gov (United States)

Onishi, Ryo; Kunishima, Yuichi

2017-11-01

We have been developing the Lagrangian Cloud Simulator (LCS), which follows the so-called Euler-Lagrangian framework, where flow motion and scalar transportations (i.e., temperature and humidity) are computed with the Euler method and particle motion with the Lagrangian method. The LCS simulation considers the hydrodynamic interaction between approaching particles for robust collision detection. This leads to reliable simulations of collision growth of cloud droplets. Recently the activation process, in which aerosol particles become tiny liquid droplets, has been implemented in the LCS. The present LCS can therefore consider the whole warm-rain precipitation processes -activation, condensation, collision and drop precipitation. In this talk, after briefly introducing the LCS, we will show kinematic simulations using the LCS for quasi-one dimensional domain, i.e., vertically elongated 3D domain. They are compared with one-dimensional kinematic simulations using a spectral-bin cloud microphysics scheme, which is based on the Euler method. The comparisons show fairly good agreement with small discrepancies, the source of which will be presented. The Lagrangian statistics, obtained for the first time for the vertical domain, will be the center of discussion. This research was supported by MEXT as ``Exploratory Challenge on Post-K computer'' (Frontiers of Basic Science: Challenging the Limits).

Dimensional analysis and prediction of dielectrophoretic crossover frequency of spherical particles

Directory of Open Access Journals (Sweden)

Che-Kai Yeh

2017-06-01

Full Text Available The manipulation of biological cells and micrometer-scale particles using dielectrophoresis (DEP is an indispensable technique for lab-on-a-chip systems for many biological and colloidal science applications. However, existing models, including the dipole model and numerical simulations based on Maxwell stress tensor (MST, cannot achieve high accuracy and high computation efficiency at the same time. The dipole model is widely used and provides adequate predictions on the crossover frequency of submicron particles, but cannot predict the crossover frequency for larger particles accurately; on the other hand, the MST method offers high accuracy for a wide variety of particle sizes and shapes, but is time-consuming and may lack predictive understanding of the interplay between key parameters. Here we present a mathematical model, using dimensional analysis and the Buckingham pi theorem, that permits high accuracy and efficiency in predicting the crossover frequency of spherical particles. The curve fitting and calculation are performed using commercial packages OriginLab and MATLAB, respectively. In addition, through this model we also can predict the conditions in which no crossover frequency exists. Also, we propose a pair of dimensionless parameters, forming a functional relation, that provide physical insights into the dependency of the crossover frequency on five key parameters. The model is verified under several scenarios using comprehensive MST simulations by COMSOL Multiphysics software (COMSOL, Inc. and some published experimental data.

Quantifying the mechanical micro-environment during three-dimensional cell expansion on microbeads by means of individual cell-based modelling.

Science.gov (United States)

Smeets, Bart; Odenthal, Tim; Tijskens, Engelbert; Ramon, Herman; Van Oosterwyck, Hans

2013-10-01

Controlled in vitro three-dimensional cell expansion requires culture conditions that optimise the biophysical micro-environment of the cells during proliferation. In this study, we propose an individual cell-based modelling platform for simulating the mechanics of cell expansion on microcarriers. The lattice-free, particle-based method considers cells as individual interacting particles that deform and move over time. The model quantifies how the mechanical micro-environment of individual cells changes during the time of confluency. A sensitivity analysis is performed, which shows that changes in the cell-specific properties of cell-cell adhesion and cell stiffness cause the strongest change in the mechanical micro-environment of the cells. Furthermore, the influence of the mechanical properties of cells and microbead is characterised. The mechanical micro-environment is strongly influenced by the adhesive properties and the size of the microbead. Simulations show that even in the absence of strong biological heterogeneity, a large heterogeneity in mechanical stresses can be expected purely due to geometric properties of the culture system.

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.

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

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

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.

Combining neural networks and signed particles to simulate quantum systems more efficiently

Science.gov (United States)

Sellier, Jean Michel

2018-04-01

Recently a new formulation of quantum mechanics has been suggested which describes systems by means of ensembles of classical particles provided with a sign. This novel approach mainly consists of two steps: the computation of the Wigner kernel, a multi-dimensional function describing the effects of the potential over the system, and the field-less evolution of the particles which eventually create new signed particles in the process. Although this method has proved to be extremely advantageous in terms of computational resources - as a matter of fact it is able to simulate in a time-dependent fashion many-body systems on relatively small machines - the Wigner kernel can represent the bottleneck of simulations of certain systems. Moreover, storing the kernel can be another issue as the amount of memory needed is cursed by the dimensionality of the system. In this work, we introduce a new technique which drastically reduces the computation time and memory requirement to simulate time-dependent quantum systems which is based on the use of an appropriately tailored neural network combined with the signed particle formalism. In particular, the suggested neural network is able to compute efficiently and reliably the Wigner kernel without any training as its entire set of weights and biases is specified by analytical formulas. As a consequence, the amount of memory for quantum simulations radically drops since the kernel does not need to be stored anymore as it is now computed by the neural network itself, only on the cells of the (discretized) phase-space which are occupied by particles. As its is clearly shown in the final part of this paper, not only this novel approach drastically reduces the computational time, it also remains accurate. The author believes this work opens the way towards effective design of quantum devices, with incredible practical implications.

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.

Vectorization of a particle simulation method for hypersonic rarefied flow

Science.gov (United States)

Mcdonald, Jeffrey D.; Baganoff, Donald

1988-01-01

An efficient particle simulation technique for hypersonic rarefied flows is presented at an algorithmic and implementation level. The implementation is for a vector computer architecture, specifically the Cray-2. The method models an ideal diatomic Maxwell molecule with three translational and two rotational degrees of freedom. Algorithms are designed specifically for compatibility with fine grain parallelism by reducing the number of data dependencies in the computation. By insisting on this compatibility, the method is capable of performing simulation on a much larger scale than previously possible. A two-dimensional simulation of supersonic flow over a wedge is carried out for the near-continuum limit where the gas is in equilibrium and the ideal solution can be used as a check on the accuracy of the gas model employed in the method. Also, a three-dimensional, Mach 8, rarefied flow about a finite-span flat plate at a 45 degree angle of attack was simulated. It utilized over 10 to the 7th particles carried through 400 discrete time steps in less than one hour of Cray-2 CPU time. This problem was chosen to exhibit the capability of the method in handling a large number of particles and a true three-dimensional geometry.

Vectorization of a particle simulation method for hypersonic rarefied flow

International Nuclear Information System (INIS)

Mcdonald, J.D.; Baganoff, D.

1988-01-01

An efficient particle simulation technique for hypersonic rarefied flows is presented at an algorithmic and implementation level. The implementation is for a vector computer architecture, specifically the Cray-2. The method models an ideal diatomic Maxwell molecule with three translational and two rotational degrees of freedom. Algorithms are designed specifically for compatibility with fine grain parallelism by reducing the number of data dependencies in the computation. By insisting on this compatibility, the method is capable of performing simulation on a much larger scale than previously possible. A two-dimensional simulation of supersonic flow over a wedge is carried out for the near-continuum limit where the gas is in equilibrium and the ideal solution can be used as a check on the accuracy of the gas model employed in the method. Also, a three-dimensional, Mach 8, rarefied flow about a finite-span flat plate at a 45 degree angle of attack was simulated. It utilized over 10 to the 7th particles carried through 400 discrete time steps in less than one hour of Cray-2 CPU time. This problem was chosen to exhibit the capability of the method in handling a large number of particles and a true three-dimensional geometry. 14 references

Macroscale implicit electromagnetic particle simulation of magnetized plasmas

International Nuclear Information System (INIS)

Tanaka, Motohiko.

1988-01-01

An electromagnetic and multi-dimensional macroscale particle simulation code (MACROS) is presented which enables us to make a large time and spatial scale kinetic simulation of magnetized plasmas. Particle ions, finite mass electrons with the guiding-center approximation and a complete set of Maxwell equations are employed. Implicit field-particle coupled equations are derived in which a time-decentered (slightly backward) finite differential scheme is used to achieve stability for large time and spatial scales. It is shown analytically that the present simulation scheme suppresses high frequency electromagnetic waves and that it accurately reproduces low frequency waves in the plasma. These properties are verified by numerical examination of eigenmodes in a 2-D thermal equilibrium plasma and by that of the kinetic Alfven wave. (author)

One-dimensional plasma simulation studies

International Nuclear Information System (INIS)

Friberg, Ari; Virtamo, Jorma

1976-01-01

Some basic plasma phenomena are studied by a one-dimensional electrostatic simulation code. A brief description of the code and its application to a test problem is given. The experiments carried out include Landau damping of an excited wave, particle retardation by smoothed field and beam-plasma instability. In each case, the set-up of the experiment is described and the results are compared with theoretical predictions. In the theoretical discussions, the oscillatory behaviour found in the Landau damping experiment is explained, an explicit formula for the particle retardation rate is derived and a rudimentary picture of the beam-plasma instability in terms of quasilinear theory is given. (author)

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

PAM: Particle automata model in simulation of Fusarium graminearum pathogen expansion.

Science.gov (United States)

Wcisło, Rafał; Miller, S Shea; Dzwinel, Witold

2016-01-21

The multi-scale nature and inherent complexity of biological systems are a great challenge for computer modeling and classical modeling paradigms. We present a novel particle automata modeling metaphor in the context of developing a 3D model of Fusarium graminearum infection in wheat. The system consisting of the host plant and Fusarium pathogen cells can be represented by an ensemble of discrete particles defined by a set of attributes. The cells-particles can interact with each other mimicking mechanical resistance of the cell walls and cell coalescence. The particles can move, while some of their attributes can be changed according to prescribed rules. The rules can represent cellular scales of a complex system, while the integrated particle automata model (PAM) simulates its overall multi-scale behavior. We show that due to the ability of mimicking mechanical interactions of Fusarium tip cells with the host tissue, the model is able to simulate realistic penetration properties of the colonization process reproducing both vertical and lateral Fusarium invasion scenarios. The comparison of simulation results with micrographs from laboratory experiments shows encouraging qualitative agreement between the two. Copyright © 2015 Elsevier Ltd. All rights reserved.

Finite element approach to global gyrokinetic particle-in-cell simulations using magnetic coordinate

International Nuclear Information System (INIS)

Fivaz, M.; Brunner, S.; Ridder, G. de; Sauter, O.; Tran, T.M.; Vaclavik, J.; Villard, L.; Appert, K.

1997-08-01

We present a fully-global linear gyrokinetic simulation code (GYGLES) aimed at describing the instable spectrum of the ion-temperature-gradient modes in toroidal geometry. We formulate the Particle-In-Cell method with finite elements defined in magnetic coordinates, which provides excellent numerical convergence properties. The poloidal mode structure corresponding to k // =0 is extracted without approximation from the equations, which reduces drastically the numerical resolution needed. The code can simulate routinely modes with both very long and very short toroidal wavelengths, can treat realistic (MHD) equilibria of any size and runs on a massively parallel computer. (author) 10 figs., 28 refs

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.

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

Effects of dimensionality and laser polarization on kinetic simulations of laser-ion acceleration in the transparency regime

Science.gov (United States)

Stark, David; Yin, Lin; Albright, Brian; Guo, Fan

2017-10-01

The often cost-prohibitive nature of three-dimensional (3D) kinetic simulations of laser-plasma interactions has resulted in heavy use of two-dimensional (2D) simulations to extract physics. However, depending on whether the polarization is modeled as 2D-S or 2D-P (laser polarization in and out of the simulation plane, respectively), different results arise. In laser-ion acceleration in the transparency regime, VPIC particle-in-cell simulations show that 2D-S and 2D-P capture different physics that appears in 3D simulations. The electron momentum distribution is virtually two-dimensional in 2D-P, unlike the more isotropic distributions in 2D-S and 3D, leading to greater heating in the simulation plane. As a result, target expansion time scales and density thresholds for the onset of relativistic transparency differ dramatically between 2D-S and 2D-P. The artificial electron heating in 2D-P exaggerates the effectiveness of target-normal sheath acceleration (TNSA) into its dominant acceleration mechanism, whereas 2D-S and 3D both have populations accelerated preferentially during transparency to higher energies than those of TNSA. Funded by the LANL Directed Research and Development Program.

Lattice-Boltzmann Method with Dynamic Grid Refinement for Simulating Particle Deposition on a Single Fibre

Directory of Open Access Journals (Sweden)

Helmut Schomburg

2013-03-01

Full Text Available In this work a numerical approach to predict the deposition behaviour of nano-scale particles on the surface of a single fibre by resolving the resulting dendrite-like particle structures in detail is presented. The gas flow simulation is carried out by a two-dimensional Lattice-Boltzmann method, which is coupled with a Lagrangian approach for the particle motion. To decrease calculation time and system requirements the Lattice-Boltzmann model is extended to allow for local grid refinement. Because of the a priori unknown location of deposition, the simulation procedure starts on a coarse mesh which is then locally refined in a fully adaptive way in regions of accumulated particles. After each deposition the fluid flow is recalculated in order to resolve the coupling of the flow with the growing particle structures correctly. For the purpose of avoiding unphysical blocking of flow by growing particle dendrites the Lattice-Boltzmann method is extended to permeable cells in these regions using the Brinkmann equation. This extended deposition model is compared to simpler approaches, where the deposit has no retroaction on the flow or is treated as a solid structure. It is clear that the permeable model is most realistic and allows considering the particle deposition on a fibre as two-dimensional problem. Comprehensive simulations were conducted for analysing the importance of different parameters, i.e. free-stream velocity and particle diameter on the deposit structure. The results of this sensitivity analysis agree qualitatively well with former published numerical and experimental results. Finally the structure of the particle deposit was quantitatively characterised by using a modified fractal dimension.

Simulation of high-energy particle production through sausage and kink instabilities in pinched plasma discharges

International Nuclear Information System (INIS)

Haruki, Takayuki; Yousefi, Hamid Reza; Masugata, Katsumi; Sakai, Jun-Ichi; Mizuguchi, Yusuke; Makino, Nao; Ito, Hiroaki

2006-01-01

In an experimental plasma, high-energy particles were observed by using a plasma focus device, to obtain energies of a few hundred keV for electrons, up to MeV for ions. In order to study the mechanism of high-energy particle production in pinched plasma discharges, a numerical simulation was introduced. By use of a three-dimensional relativistic and fully electromagnetic particle-in-cell code, the dynamics of a Z-pinch plasma, thought to be unstable against sausage and kink instabilities, are investigated. In this work, the development of sausage and kink instabilities and subsequent high-energy particle production are shown. In the model used here, cylindrically distributed electrons and ions are driven by an external electric field. The driven particles spontaneously produce a current, which begins to pinch by the Lorentz force. Initially the pinched current is unstable against a sausage instability, and then becomes unstable against a kink instability. As a result high-energy particles are observed

Investigating the settling dynamics of cohesive silt particles with particle-resolving simulations

Science.gov (United States)

Sun, Rui; Xiao, Heng; Sun, Honglei

2018-01-01

The settling of cohesive sediment is ubiquitous in aquatic environments, and the study of the settling process is important for both engineering and environmental reasons. In the settling process, the silt particles show behaviors that are different from non-cohesive particles due to the influence of inter-particle cohesive force. For instance, the flocs formed in the settling process of cohesive silt can loosen the packing, and thus the structural densities of cohesive silt beds are much smaller than that of non-cohesive sand beds. While there is a consensus that cohesive behaviors depend on the characteristics of sediment particles (e.g., Bond number, particle size distribution), little is known about the exact influence of these characteristics on the cohesive behaviors. In addition, since the cohesive behaviors of the silt are caused by the inter-particle cohesive forces, the motions of and the contacts among silt particles should be resolved to study these cohesive behaviors in the settling process. However, studies of the cohesive behaviors of silt particles in the settling process based on particle-resolving approach are still lacking. In the present work, three-dimensional settling process is investigated numerically by using CFD-DEM (Computational Fluid Dynamics-Discrete Element Method). The inter-particle collision force, the van der Waals force, and the fluid-particle interaction forces are considered. The numerical model is used to simulate the hindered settling process of silt based on the experimental setup in the literature. The results obtained in the simulations, including the structural densities of the beds, the characteristic lines, and the particle terminal velocity, are in good agreement with the experimental observations in the literature. To the authors' knowledge, this is the first time that the influences of non-dimensional Bond number and particle polydispersity on the structural densities of silt beds have been investigated separately

An arbitrary curvilinear-coordinate method for particle-in-cell modeling

International Nuclear Information System (INIS)

Fichtl, C A; Finn, J M; Cartwright, K L

2012-01-01

A new approach to kinetic simulation of plasmas in complex geometries, based on the particle-in-cell (PIC) simulation method, is explored. In the two-dimensional (2D) electrostatic version of our method, called the arbitrary curvilinear-coordinate PIC method, all essential PIC operations are carried out in 2D on a uniform grid on the unit square logical domain, and mapped to a nonuniform boundary-fitted grid on the physical domain. As the resulting logical grid equations of motion are not separable, we have developed an extension of the semi-implicit modified leapfrog integration technique to preserve the symplectic nature of the logical grid particle mover. A generalized, curvilinear-coordinate formulation of Poisson's equations to solve for the electrostatic fields on the uniform logical grid is also developed. By our formulation, we compute the plasma charge density on the logical grid based on the particles' positions on the logical domain. That is, the plasma particles are weighted to the uniform logical grid and the self-consistent mean electrostatic fields obtained from the solution of the logical grid Poisson equation are interpolated to the particle positions on the logical grid. This process eliminates the complexity associated with the weighting and interpolation processes on the nonuniform physical grid and allows us to run the PIC method on arbitrary boundary-fitted meshes. (paper)

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.

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)

HIDENEK: an implicit particle simulation of kinetic-MHD phenomena in three-dimensional plasmas

International Nuclear Information System (INIS)

Tanaka, Motohiko.

1993-05-01

An advanced 'kinetic-MHD' simulation method and its applications to plasma physics are given in this lecture. This method is quite suitable for studying strong nonlinear, kinetic processes associated with large space-scale, low-frequency electromagnetic phenomena of plasmas. A full set of the Maxwell equations, and the Newton-Lorentz equations of motion for particle ions and guiding-center electrons are adopted. In order to retain only the low-frequency waves and instabilities, implicit particle-field equations are derived. The present implicit-particle method is proved to reproduce the MHD eigenmodes such as Alfven, magnetosonic and kinetic Alfven waves in a thermally near-equilibrium plasma. In the second part of the lecture, several physics applications are shown. These include not only the growth of the instabilities of beam ions against the background plasmas and helical kink of the current, but they also demonstrate nonlinear results such as pitch-angle scattering of the ions. Recent progress in the simulation of the Kelvin-Helmholtz instability is also presented with a special emphasis on the mixing of plasma particles. (author)

Particle-in-cell Simulations of Raman Laser Amplification in Ionizing Plasmas

International Nuclear Information System (INIS)

Clark, Daniel S.; Fisch, Nathaniel J.

2003-01-01

By using the amplifying laser pulse in a plasma-based backward Raman laser amplifier to generate the plasma by photo-ionization of a gas simultaneous with the amplification process, possible instabilities of the pumping laser pulse can be avoided. Particle-in-cell simulations are used to study this amplification mechanism, and earlier results using more elementary models of the Raman interaction are verified [D.S. Clark and N.J. Fisch, Phys. Plasmas, 9 (6): 2772-2780, 2002]. The effects (unique to amplification in ionizing plasmas and not included in previous simulations) of blue-shifting of the pump and seed laser pulses and the generation of a wake are observed not significantly to impact the amplification process. As expected theoretically, the peak output intensity is found to be limited to I ∼ 10 17 W/cm 2 by forward Raman scattering of the amplifying seed. The integrity of the ionization front of the seed pulse against the development of a possible transverse modulation instability is also demonstrated

A Novel Microgravity Simulator Applicable for Three-Dimensional Cell Culturing

Science.gov (United States)

Wuest, Simon L.; Richard, Stéphane; Walther, Isabelle; Furrer, Reinhard; Anderegg, Roland; Sekler, Jörg; Egli, Marcel

2014-10-01

Random Positioning Machines (RPM) were introduced decades ago to simulate microgravity. Since then numerous experiments have been carried out to study its influence on biological samples. The machine is valued by the scientific community involved in space relevant topics as an excellent experimental tool to conduct pre-studies, for example, before sending samples into space. We have developed a novel version of the traditional RPM to broaden its operative range. This novel version has now become interesting to researchers who are working in the field of tissue engineering, particularly those interested in alternative methods for three-dimensional (3D) cell culturing. The main modifications concern the cell culture condition and the algorithm that controls the movement of the frames for the nullification of gravity. An incubator was integrated into the inner frame of the RPM allowing precise control over the cell culture environment. Furthermore, several feed-throughs now allow a permanent supply of gas like CO 2. All these modifications substantially improve conditions to culture cells; furthermore, the rewritten software responsible for controlling the movement of the frames enhances the quality of the generated microgravity. Cell culture experiments were carried out with human lymphocytes on the novel RPM model to compare the obtained response to the results gathered on an older well-established RPM as well as to data from space flights. The overall outcome of the tests validates this novel RPM for cell cultivation under simulated microgravity conditions.

Multi-Algorithm Particle Simulations with Spatiocyte.

Science.gov (United States)

Arjunan, Satya N V; Takahashi, Koichi

2017-01-01

As quantitative biologists get more measurements of spatially regulated systems such as cell division and polarization, simulation of reaction and diffusion of proteins using the data is becoming increasingly relevant to uncover the mechanisms underlying the systems. Spatiocyte is a lattice-based stochastic particle simulator for biochemical reaction and diffusion processes. Simulations can be performed at single molecule and compartment spatial scales simultaneously. Molecules can diffuse and react in 1D (filament), 2D (membrane), and 3D (cytosol) compartments. The implications of crowded regions in the cell can be investigated because each diffusing molecule has spatial dimensions. Spatiocyte adopts multi-algorithm and multi-timescale frameworks to simulate models that simultaneously employ deterministic, stochastic, and particle reaction-diffusion algorithms. Comparison of light microscopy images to simulation snapshots is supported by Spatiocyte microscopy visualization and molecule tagging features. Spatiocyte is open-source software and is freely available at http://spatiocyte.org .

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

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

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 Simulations of Asymmetric Dual Frequency Capacitive Discharge Physics

Science.gov (United States)

Wu, Alan; Lichtenberg, A. J.; Lieberman, M. A.; Verboncoeur, J. P.

2003-10-01

Dual frequency capacitive discharges are finding increasing use for etching in the microelectronics industry. In the ideal case, the high frequency power (typically 27.1-160 MHz) controls the plasma density and the low frequency power (typically 2-13.56 MHz) controls the ion energy. The electron power deposition and the dynamics of dual frequency rf sheaths are not well understood. We report on particle-in-cell computer simulations of an asymmetric dual frequency argon discharge. The simulations are performed in 1D (radial) geometry using the bounded electrostatic code XPDP1. Operating parameters are 27.1/2 MHz high/low frequencies, 10/13 cm inner/outer radii, 3-200 mTorr pressures, and 10^9-10^11 cm-3 densities. We determine the power deposition and sheath dynamics for the high frequency power alone, and with various added low frequency powers. We compare the simulation results to simple global models of dual frequency discharges. Support provided by Lam Research, NSF Grant ECS-0139956, California industries, and UC-SMART Contract SM99-10051.

photon-plasma: A modern high-order particle-in-cell code

International Nuclear Information System (INIS)

Haugbølle, Troels; Frederiksen, Jacob Trier; Nordlund, Åke

2013-01-01

We present the photon-plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach, the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks

A variational Bayesian multiple particle filtering scheme for large-dimensional systems

KAUST Repository

Ait-El-Fquih, Boujemaa

2016-06-14

This paper considers the Bayesian filtering problem in high-dimensional nonlinear state-space systems. In such systems, classical particle filters (PFs) are impractical due to the prohibitive number of required particles to obtain reasonable performances. One approach that has been introduced to overcome this problem is the concept of multiple PFs (MPFs), where the state-space is split into low-dimensional subspaces and then a separate PF is applied to each subspace. Remarkable performances of MPF-like filters motivated our investigation here into a new strategy that combines the variational Bayesian approach to split the state-space with random sampling techniques, to derive a new computationally efficient MPF. The propagation of each particle in the prediction step of the resulting filter requires generating only a single particle in contrast with standard MPFs, for which a set of (children) particles is required. We present simulation results to evaluate the behavior of the proposed filter and compare its performances against standard PF and a MPF.

A variational Bayesian multiple particle filtering scheme for large-dimensional systems

KAUST Repository

Ait-El-Fquih, Boujemaa; Hoteit, Ibrahim

2016-01-01

This paper considers the Bayesian filtering problem in high-dimensional nonlinear state-space systems. In such systems, classical particle filters (PFs) are impractical due to the prohibitive number of required particles to obtain reasonable performances. One approach that has been introduced to overcome this problem is the concept of multiple PFs (MPFs), where the state-space is split into low-dimensional subspaces and then a separate PF is applied to each subspace. Remarkable performances of MPF-like filters motivated our investigation here into a new strategy that combines the variational Bayesian approach to split the state-space with random sampling techniques, to derive a new computationally efficient MPF. The propagation of each particle in the prediction step of the resulting filter requires generating only a single particle in contrast with standard MPFs, for which a set of (children) particles is required. We present simulation results to evaluate the behavior of the proposed filter and compare its performances against standard PF and a MPF.

Revealing the sub-structures of the magnetic reconnection separatrix via particle-in-cell simulation

International Nuclear Information System (INIS)

Zhou, M.; Deng, X. H.; Pang, Y.; Xu, X. J.; Yao, M.; Huang, S. Y.; Yuan, Z. G.; Li, H. M.; Wang, D. D.; Wang, Y. H.

2012-01-01

Magnetic separatrix is an important boundary layer separating the inflow and outflow regions in magnetic reconnection. In this article, we investigate the sub-structures of the separatrix region by using two-and-half dimensional electromagnetic particle-in-cell simulation. The separatrix region can be divided into two sub-regions in terms of the ion and electron frozen-in conditions. Far from the neutral sheet, ions and electrons are magnetized in magnetic fields. Approaching the neutral sheet, ion frozen-in condition is broken in a narrow region (∼c/ω pi ) at the edge of a density cavity, while electrons are frozen-in to magnetic fields. In this region, electric field E z is around zero, and the convective term –(v i × B) is balanced by the Hall term in the generalized Ohm’s law because ions carry the perpendicular current. Inside the density cavity, both ion and electron frozen-in conditions are broken. The region consists of two sub-ion or electron-scale layers, which contain intense electric fields. Formation of the two sub-layers is due to the complex electron flow pattern around the separatrix region. In the layer, E z is balanced by a combination of Hall term and the divergence of electron pressure tensor, with the Hall term being dominant. Our preliminary simulation result shows that the separatrix region in guide field reconnection also contains two sub-regions: the inner region and the outer region. However, the inner region contains only one current layer in contrast with the case without guide field.

Particle simulation of intense electron cyclotron heating and beat-wave current drive

International Nuclear Information System (INIS)

Cohen, B.I.

1987-01-01

High-power free-electron lasers make new methods possible for heating plasmas and driving current in toroidal plasmas with electromagnetic waves. We have undertaken particle simulation studies with one and two dimensional, relativistic particle simulation codes of intense pulsed electron cyclotron heating and beat-wave current drive. The particle simulation methods here are conventional: the algorithms are time-centered, second-order-accurate, explicit, leap-frog difference schemes. The use of conventional methods restricts the range of space and time scales to be relatively compact in the problems addressed. Nevertheless, experimentally relevant simulations have been performed. 10 refs., 2 figs

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

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.

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

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

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.

PARTICLE-IN-CELL SIMULATION OF A STRONG DOUBLE LAYER IN A NONRELATIVISTIC PLASMA FLOW: ELECTRON ACCELERATION TO ULTRARELATIVISTIC SPEEDS

International Nuclear Information System (INIS)

Dieckmann, Mark E.; Bret, Antoine

2009-01-01

Two charge- and current-neutral plasma beams are modeled with a one-dimensional particle-in-cell simulation. The beams are uniform and unbounded. The relative speed between both beams is 0.4c. One beam is composed of electrons and protons, and the other of protons and negatively charged oxygen (dust). All species have the temperature 9.1 keV. A Buneman instability develops between the electrons of the first beam and the protons of the second beam. The wave traps the electrons, which form plasmons. The plasmons couple energy into the ion acoustic waves, which trap the protons of the second beam. A structure similar to a proton phase-space hole develops, which grows through its interaction with the oxygen and the heated electrons into a rarefaction pulse. This pulse drives a double layer, which accelerates a beam of electrons to about 50 MeV, which is comparable to the proton kinetic energy. The proton distribution eventually evolves into an electrostatic shock. Beams of charged particles moving at such speeds may occur in the foreshock of supernova remnant (SNR) shocks. This double layer is thus potentially relevant for the electron acceleration (injection) into the diffusive shock acceleration by SNR shocks.

Stochastic plasma heating by electrostatic waves: a comparison between a particle-in-cell simulation and a laboratory experiment

International Nuclear Information System (INIS)

Fivaz, M.; Fasoli, A.; Appert, K.; Trans, T.M.; Tran, M.Q.; Skiff, F.

1993-08-01

Dynamical chaos is produced by the interaction between plasma particles and two electrostatic waves. Experiments performed in a linear magnetized plasma and a 1D particle-in-cell simulation agree qualitatively: above a threshold wave amplitude, ion stochastic diffusion and heating occur on a fast time scale. Self-consistency appears to limit the extent of the heating process. (author) 5 figs., 18 refs

Particle-in-cell simulations of the lasertron

International Nuclear Information System (INIS)

Jones, M.E.; Peter, W.K.

1985-01-01

The lasertron is a device (either rf or dc) for producing intense, very short, pulsed electron beams (tens of picoseconds). In the dc lasertron, a laser is pulsed repetitively onto a photoemissive cathode. In general, the current is not space-charge limited and follows the laser intensity. The electron pulse is then accelerated out of the device by a constant voltage. By using the laser the need for a subharmonic buncher is eliminated. In the rf lasertron, the diode becomes an rf cavity. This improves the breakdown characteristics of the device, allowing higher voltages to be applied and hence higher currents to be obtained. The calculations are aimed at producing a 10 nC electron beam with an emittance of less than 40π mm-mrad for use in free-electron laser experiments at Los Alamos. Other applications of the lasertron include efficient microwave or rf generation. A class of electrode shapes has been obtained which in the absence of space charge produces no emittance growth. These shapes have been studied with the particle-in-cell simulation model ISIS, and the electrodes which produce minimum emittance including the effect of space charge have been determined. Unique emittance problems associated with the time dependence of the beam pulse are studied and conditions for reducing these effects are discussed. 5 refs., 5 figs

Particle-In-Cell Simulations of the Solar Wind Interaction with Lunar Crustal Magnetic Anomalies: Magnetic Cusp Regions

Science.gov (United States)

Poppe, A. R.; Halekas, J. S.; Delory, G. T.; Farrell, W. M.

2012-01-01

As the solar wind is incident upon the lunar surface, it will occasionally encounter lunar crustal remanent magnetic fields. These magnetic fields are small-scale, highly non-dipolar, have strengths up to hundreds of nanotesla, and typically interact with the solar wind in a kinetic fashion. Simulations, theoretical analyses, and spacecraft observations have shown that crustal fields can reflect solar wind protons via a combination of magnetic and electrostatic reflection; however, analyses of surface properties have suggested that protons may still access the lunar surface in the cusp regions of crustal magnetic fields. In this first report from a planned series of studies, we use a 1 1/2-dimensional, electrostatic particle-in-cell code to model the self-consistent interaction between the solar wind, the cusp regions of lunar crustal remanent magnetic fields, and the lunar surface. We describe the self-consistent electrostatic environment within crustal cusp regions and discuss the implications of this work for the role that crustal fields may play regulating space weathering of the lunar surface via proton bombardment.

Optimization of Particle Search Algorithm for CFD-DEM Simulations

Directory of Open Access Journals (Sweden)

G. Baryshev

2013-09-01

Full Text Available Discrete element method has numerous applications in particle physics. However, simulating particles as discrete entities can become costly for large systems. In time-driven DEM simulation most computation time is taken by contact search stage. We propose an efficient collision detection method which is based on sorting particles by their coordinates. Using multiple sorting criteria allows minimizing number of potential neighbours and defines fitness of this approach for simulation of massive systems in 3D. This method is compared to a common approach that consists of placing particles onto a grid of cells. Advantage of the new approach is independence of simulation parameters upon particle radius and domain size.

Explicit simulation of ice particle habits in a Numerical Weather Prediction Model

Science.gov (United States)

Hashino, Tempei

2007-05-01

This study developed a scheme for explicit simulation of ice particle habits in Numerical Weather Prediction (NWP) Models. The scheme is called Spectral Ice Habit Prediction System (SHIPS), and the goal is to retain growth history of ice particles in the Eulerian dynamics framework. It diagnoses characteristics of ice particles based on a series of particle property variables (PPVs) that reflect history of microphysieal processes and the transport between mass bins and air parcels in space. Therefore, categorization of ice particles typically used in bulk microphysical parameterization and traditional bin models is not necessary, so that errors that stem from the categorization can be avoided. SHIPS predicts polycrystals as well as hexagonal monocrystals based on empirically derived habit frequency and growth rate, and simulates the habit-dependent aggregation and riming processes by use of the stochastic collection equation with predicted PPVs. Idealized two dimensional simulations were performed with SHIPS in a NWP model. The predicted spatial distribution of ice particle habits and types, and evolution of particle size distributions showed good quantitative agreement with observation This comprehensive model of ice particle properties, distributions, and evolution in clouds can be used to better understand problems facing wide range of research disciplines, including microphysics processes, radiative transfer in a cloudy atmosphere, data assimilation, and weather modification.

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

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

Self-diffusion in monodisperse three-dimensional magnetic fluids by molecular dynamics simulations

Energy Technology Data Exchange (ETDEWEB)

Dobroserdova, A.B. [Ural Federal University, Lenin Av. 51, Ekaterinburg (Russian Federation); Kantorovich, S.S., E-mail: alla.dobroserdova@urfu.ru [Ural Federal University, Lenin Av. 51, Ekaterinburg (Russian Federation); University of Vienna, Sensengasse 8, Vienna (Austria)

2017-06-01

In the present work we study the self-diffusion behaviour in the three-dimensional monodisperse magnetic fluids using the Molecular Dynamics Simulation and Density Functional Theory. The peculiarity of computer simulation is to study two different systems: dipolar and soft sphere ones. In the theoretical method, it is important to choose the approximation for the main structures, which are chains. We compare the theoretical results and the computer simulation data for the self-diffusion coefficient as a function of the particle volume fraction and magnetic dipole-dipole interaction parameter and find the qualitative and quantitative agreement to be good. - Highlights: • The paper deals with the study of the self-diffusion in monodisperse three-dimensional magnetic fluids. • The theoretical approach contains the free energy density functional minimization. • Computer simulations are performed by the molecular dynamics method. • We have a good qualitative and quantitative agreement between the theoretical results and computer simulation data.

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

Three-dimensional simulation of grain mixing in three different rotating drum designs for solid-state fermentation

NARCIS (Netherlands)

Schutyser, M.A.I.; Weber, F.J.; Briels, W.J.; Boom, R.M.; Rinzema, A.

2002-01-01

A previously published two-dimensional discrete particle simulation model for radial mixing behavior of various slowly rotating drums for solid-state fermentation (SSF) has been extended to a three-dimensional model that also predicts axial mixing. Radial and axial mixing characteristics were

Kinetic structures of quasi-perpendicular shocks in global particle-in-cell simulations

International Nuclear Information System (INIS)

Peng, Ivy Bo; Markidis, Stefano; Laure, Erwin; Johlander, Andreas; Vaivads, Andris; Khotyaintsev, Yuri; Henri, Pierre; Lapenta, Giovanni

2015-01-01

We carried out global Particle-in-Cell simulations of the interaction between the solar wind and a magnetosphere to study the kinetic collisionless physics in super-critical quasi-perpendicular shocks. After an initial simulation transient, a collisionless bow shock forms as a result of the interaction of the solar wind and a planet magnetic dipole. The shock ramp has a thickness of approximately one ion skin depth and is followed by a trailing wave train in the shock downstream. At the downstream edge of the bow shock, whistler waves propagate along the magnetic field lines and the presence of electron cyclotron waves has been identified. A small part of the solar wind ion population is specularly reflected by the shock while a larger part is deflected and heated by the shock. Solar wind ions and electrons are heated in the perpendicular directions. Ions are accelerated in the perpendicular direction in the trailing wave train region. This work is an initial effort to study the electron and ion kinetic effects developed near the bow shock in a realistic magnetic field configuration

Propagation of localized structures in relativistic magnetized electron-positron plasmas using particle-in-cell simulations

Energy Technology Data Exchange (ETDEWEB)

López, Rodrigo A. [Departamento de Física, Facultad de Ciencias Físicas y Matemáticas, Universidad de Concepción, Concepción 4070386 (Chile); Muñoz, Víctor [Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago (Chile); Viñas, Adolfo F. [Geospace Physics Laboratory, Heliophysics Science Division, NASA Goddard Space Flight Center, Greenbelt, Maryland 20771 (United States); Valdivia, Juan A. [Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago (Chile); Centro para el Desarrollo de la Nanociencia y la Nanotecnología (CEDENNA), Santiago 9170124 (Chile)

2015-09-15

We use a particle-in-cell simulation to study the propagation of localized structures in a magnetized electron-positron plasma with relativistic finite temperature. We use as initial condition for the simulation an envelope soliton solution of the nonlinear Schrödinger equation, derived from the relativistic two fluid equations in the strongly magnetized limit. This envelope soliton turns out not to be a stable solution for the simulation and splits in two localized structures propagating in opposite directions. However, these two localized structures exhibit a soliton-like behavior, as they keep their profile after they collide with each other due to the periodic boundary conditions. We also observe the formation of localized structures in the evolution of a spatially uniform circularly polarized Alfvén wave. In both cases, the localized structures propagate with an amplitude independent velocity.

Propagation of localized structures in relativistic magnetized electron-positron plasmas using particle-in-cell simulations

International Nuclear Information System (INIS)

López, Rodrigo A.; Muñoz, Víctor; Viñas, Adolfo F.; Valdivia, Juan A.

2015-01-01

We use a particle-in-cell simulation to study the propagation of localized structures in a magnetized electron-positron plasma with relativistic finite temperature. We use as initial condition for the simulation an envelope soliton solution of the nonlinear Schrödinger equation, derived from the relativistic two fluid equations in the strongly magnetized limit. This envelope soliton turns out not to be a stable solution for the simulation and splits in two localized structures propagating in opposite directions. However, these two localized structures exhibit a soliton-like behavior, as they keep their profile after they collide with each other due to the periodic boundary conditions. We also observe the formation of localized structures in the evolution of a spatially uniform circularly polarized Alfvén wave. In both cases, the localized structures propagate with an amplitude independent velocity

Mesh Refinement for Particle-In-Cell Plasma Simulations: Applications to - and benefits for - Heavy-Ion-Fusion

International Nuclear Information System (INIS)

Vay, J.-L.; Colella, P.; McCorquodale, P.; Van Straalen, B.; 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 simulation of the power plant as a whole, or even of the driver, is not yet possible. Despite the rapid progress in computer power, past and anticipated, one must consider the use of the most advanced numerical techniques, if we are to reach our goal expeditiously. 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 mesh refinement technique. They discuss the challenges posed by the implementation of this technique into plasma simulations (due to the presence of particles and electromagnetic waves). They will present the prospects for and projected benefits of its application to heavy ion fusion. In particular to the simulation of the ion source and the final beam propagation in the chamber. A collaboration project is under way at LBNL between the Applied Numerical Algorithms Group (ANAG) and the HIF group to couple the Adaptive Mesh Refinement (AMR) library (CHOMBO) developed by the ANAG group to the Particle-In-Cell accelerator code WARP developed by the HIF-VNL. They describe their progress and present their initial findings

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.

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 ...

Implicit particle simulation of electromagnetic plasma phenomena

International Nuclear Information System (INIS)

Kamimura, T.; Montalvo, E.; Barnes, D.C.; Leboeuf, J.N.; Tajima, T.

1986-11-01

A direct method for the implicit particle simulation of electromagnetic phenomena in magnetized, multi-dimensional plasmas is developed. The method is second-order accurate for ωΔt < 1, with ω a characteristic frequency and time step Δt. Direct time integration of the implicit equations with simplified space differencing allows the consistent inclusion of finite particle size. Decentered time differencing of the Lorentz force permits the efficient simulation of strongly magnetized plasmas. A Fourier-space iterative technique for solving the implicit field corrector equation, based on the separation of plasma responses perpendicular and parallel to the magnetic field and longitudinal and transverse to the wavevector, is described. Wave propagation properties in a uniform plasma are in excellent agreement with theoretical expectations. Applications to collisionless tearing and coalescence instabilities further demonstrate the usefulness of the algorithm. (author)

Scientific Visualization and Simulation for Multi-dimensional Marine Environment Data

Science.gov (United States)

Su, T.; Liu, H.; Wang, W.; Song, Z.; Jia, Z.

2017-12-01

As higher attention on the ocean and rapid development of marine detection, there are increasingly demands for realistic simulation and interactive visualization of marine environment in real time. Based on advanced technology such as GPU rendering, CUDA parallel computing and rapid grid oriented strategy, a series of efficient and high-quality visualization methods, which can deal with large-scale and multi-dimensional marine data in different environmental circumstances, has been proposed in this paper. Firstly, a high-quality seawater simulation is realized by FFT algorithm, bump mapping and texture animation technology. Secondly, large-scale multi-dimensional marine hydrological environmental data is virtualized by 3d interactive technologies and volume rendering techniques. Thirdly, seabed terrain data is simulated with improved Delaunay algorithm, surface reconstruction algorithm, dynamic LOD algorithm and GPU programming techniques. Fourthly, seamless modelling in real time for both ocean and land based on digital globe is achieved by the WebGL technique to meet the requirement of web-based application. The experiments suggest that these methods can not only have a satisfying marine environment simulation effect, but also meet the rendering requirements of global multi-dimension marine data. Additionally, a simulation system for underwater oil spill is established by OSG 3D-rendering engine. It is integrated with the marine visualization method mentioned above, which shows movement processes, physical parameters, current velocity and direction for different types of deep water oil spill particle (oil spill particles, hydrates particles, gas particles, etc.) dynamically and simultaneously in multi-dimension. With such application, valuable reference and decision-making information can be provided for understanding the progress of oil spill in deep water, which is helpful for ocean disaster forecasting, warning and emergency response.

Simulations of super-structure domain walls in two dimensional assemblies of magnetic nanoparticles

DEFF Research Database (Denmark)

Jordanovic, Jelena; Beleggia, Marco; Schiøtz, Jakob

2015-01-01

We simulate the formation of domain walls in two-dimensional assemblies of magnetic nanoparticles. Particle parameters are chosen to match recent electron holography and Lorentz microscopy studies of almost monodisperse cobalt nanoparticles assembled into regular, elongated lattices. As the parti......We simulate the formation of domain walls in two-dimensional assemblies of magnetic nanoparticles. Particle parameters are chosen to match recent electron holography and Lorentz microscopy studies of almost monodisperse cobalt nanoparticles assembled into regular, elongated lattices...... taking the role of the atomic spins. The coupling is, however, different. The superspins interact only by dipolar interactions as exchange coupling between individual nanoparticles may be neglected due to interparticle spacing. We observe that it is energetically favorable to introduce domain walls...... oriented along the long dimension of nanoparticle assemblies rather than along the short dimension. This is unlike what is typically observed in continuous magnetic materials, where the exchange interaction introduces an energetic cost proportional to the area of the domain walls. Structural disorder...

Chromium poisoning in (La,Sr)MnO3 cathode: Three-dimensional simulation of a solid oxide fuel cell

OpenAIRE

Miyoshi, Kota; Iwai, Hiroshi; Kishimoto, Masashi; Saito, Motohiro; Yoshida, Hideo

2016-01-01

A three-dimensional numerical model of a single solid oxide fuel cell (SOFC) considering chromium poisoning on the cathode side has been developed to investigate the evolution of the SOFC performance over long-term operation. The degradation model applied in the simulation describes the loss of the cathode electrochemical activity as a decrease in the active triple-phase boundary (TPB) length. The calculations are conducted for two types of cell: lanthanum strontium manganite (LSM)/yttria-sta...

Effect of Finite Particle Size on Convergence of Point Particle Models in Euler-Lagrange Multiphase Dispersed Flow

Science.gov (United States)

Nili, Samaun; Park, Chanyoung; Haftka, Raphael T.; Kim, Nam H.; Balachandar, S.

2017-11-01

Point particle methods are extensively used in simulating Euler-Lagrange multiphase dispersed flow. When particles are much smaller than the Eulerian grid the point particle model is on firm theoretical ground. However, this standard approach of evaluating the gas-particle coupling at the particle center fails to converge as the Eulerian grid is reduced below particle size. We present an approach to model the interaction between particles and fluid for finite size particles that permits convergence. We use the generalized Faxen form to compute the force on a particle and compare the results against traditional point particle method. We apportion the different force components on the particle to fluid cells based on the fraction of particle volume or surface in the cell. The application is to a one-dimensional model of shock propagation through a particle-laden field at moderate volume fraction, where the convergence is achieved for a well-formulated force model and back coupling for finite size particles. Comparison with 3D direct fully resolved numerical simulations will be used to check if the approach also improves accuracy compared to the point particle model. Work supported by the U.S. Department of Energy, National Nuclear Security Administration, Advanced Simulation and Computing Program, as a Cooperative Agreement under the Predictive Science Academic Alliance Program, under Contract No. DE-NA0002378.

Simulation study of effects of initial particle size distribution on dissolution

International Nuclear Information System (INIS)

Wang, G.; Xu, D.S.; Ma, N.; Zhou, N.; Payton, E.J.; Yang, R.; Mills, M.J.; Wang, Y.

2009-01-01

Dissolution kinetics of γ' particles in binary Ni-Al alloys with different initial particle size distributions (PSD) is studied using a three-dimensional (3D) quantitative phase field model. By linking model inputs directly to thermodynamic and atomic mobility databases, microstructural evolution during dissolution is simulated in real time and length scales. The model is first validated against analytical solution for dissolution of a single γ' particle in 1D and numerical solution in 3D before it is applied to investigate the effects of initial PSD on dissolution kinetics. Four different types of PSD, uniform, normal, log-normal and bimodal, are considered. The simulation results show that the volume fraction of γ' particles decreases exponentially with time, while the temporal evolution of average particle size depends strongly on the initial PSD

A fully-implicit Particle-In-Cell Monte Carlo Collision code for the simulation of inductively coupled plasmas

Science.gov (United States)

Mattei, S.; Nishida, K.; Onai, M.; Lettry, J.; Tran, M. Q.; Hatayama, A.

2017-12-01

We present a fully-implicit electromagnetic Particle-In-Cell Monte Carlo collision code, called NINJA, written for the simulation of inductively coupled plasmas. NINJA employs a kinetic enslaved Jacobian-Free Newton Krylov method to solve self-consistently the interaction between the electromagnetic field generated by the radio-frequency coil and the plasma response. The simulated plasma includes a kinetic description of charged and neutral species as well as the collision processes between them. The algorithm allows simulations with cell sizes much larger than the Debye length and time steps in excess of the Courant-Friedrichs-Lewy condition whilst preserving the conservation of the total energy. The code is applied to the simulation of the plasma discharge of the Linac4 H- ion source at CERN. Simulation results of plasma density, temperature and EEDF are discussed and compared with optical emission spectroscopy measurements. A systematic study of the energy conservation as a function of the numerical parameters is presented.

Hybrid three-dimensional variation and particle filtering for nonlinear systems

International Nuclear Information System (INIS)

Leng Hong-Ze; Song Jun-Qiang

2013-01-01

This work addresses the problem of estimating the states of nonlinear dynamic systems with sparse observations. We present a hybrid three-dimensional variation (3DVar) and particle piltering (PF) method, which combines the advantages of 3DVar and particle-based filters. By minimizing the cost function, this approach will produce a better proposal distribution of the state. Afterwards the stochastic resampling step in standard PF can be avoided through a deterministic scheme. The simulation results show that the performance of the new method is superior to the traditional ensemble Kalman filtering (EnKF) and the standard PF, especially in highly nonlinear systems

Multi-dimensional PIC-simulations of parametric instabilities for shock-ignition conditions

Directory of Open Access Journals (Sweden)

Riconda C.

2013-11-01

Full Text Available Laser-plasma interaction is investigated for conditions relevant for the shock-ignition (SI scheme of inertial confinement fusion using two-dimensional particle-in-cell (PIC simulations of an intense laser beam propagating in a hot, large-scale, non-uniform plasma. The temporal evolution and interdependence of Raman- (SRS, and Brillouin- (SBS, side/backscattering as well as Two-Plasmon-Decay (TPD are studied. TPD is developing in concomitance with SRS creating a broad spectrum of plasma waves near the quarter-critical density. They are rapidly saturated due to plasma cavitation within a few picoseconds. The hot electron spectrum created by SRS and TPD is relatively soft, limited to energies below one hundred keV.

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.

Particle Simulation of Fractional Diffusion Equations

KAUST Repository

Allouch, Samer

2017-07-12

This work explores different particle-based approaches to the simulation of one-dimensional fractional subdiffusion equations in unbounded domains. We rely on smooth particle approximations, and consider four methods for estimating the fractional diffusion term. The first method is based on direct differentiation of the particle representation, it follows the Riesz definition of the fractional derivative and results in a non-conservative scheme. The other three methods follow the particle strength exchange (PSE) methodology and are by construction conservative, in the sense that the total particle strength is time invariant. The first PSE algorithm is based on using direct differentiation to estimate the fractional diffusion flux, and exploiting the resulting estimates in an integral representation of the divergence operator. Meanwhile, the second one relies on the regularized Riesz representation of the fractional diffusion term to derive a suitable interaction formula acting directly on the particle representation of the diffusing field. A third PSE construction is considered that exploits the Green\\'s function of the fractional diffusion equation. The performance of all four approaches is assessed for the case of a one-dimensional diffusion equation with constant diffusivity. This enables us to take advantage of known analytical solutions, and consequently conduct a detailed analysis of the performance of the methods. This includes a quantitative study of the various sources of error, namely filtering, quadrature, domain truncation, and time integration, as well as a space and time self-convergence analysis. These analyses are conducted for different values of the order of the fractional derivatives, and computational experiences are used to gain insight that can be used for generalization of the present constructions.

Particle Simulation of Fractional Diffusion Equations

KAUST Repository

Allouch, Samer; Lucchesi, Marco; Maî tre, O. P. Le; Mustapha, K. A.; Knio, Omar

2017-01-01

This work explores different particle-based approaches to the simulation of one-dimensional fractional subdiffusion equations in unbounded domains. We rely on smooth particle approximations, and consider four methods for estimating the fractional diffusion term. The first method is based on direct differentiation of the particle representation, it follows the Riesz definition of the fractional derivative and results in a non-conservative scheme. The other three methods follow the particle strength exchange (PSE) methodology and are by construction conservative, in the sense that the total particle strength is time invariant. The first PSE algorithm is based on using direct differentiation to estimate the fractional diffusion flux, and exploiting the resulting estimates in an integral representation of the divergence operator. Meanwhile, the second one relies on the regularized Riesz representation of the fractional diffusion term to derive a suitable interaction formula acting directly on the particle representation of the diffusing field. A third PSE construction is considered that exploits the Green's function of the fractional diffusion equation. The performance of all four approaches is assessed for the case of a one-dimensional diffusion equation with constant diffusivity. This enables us to take advantage of known analytical solutions, and consequently conduct a detailed analysis of the performance of the methods. This includes a quantitative study of the various sources of error, namely filtering, quadrature, domain truncation, and time integration, as well as a space and time self-convergence analysis. These analyses are conducted for different values of the order of the fractional derivatives, and computational experiences are used to gain insight that can be used for generalization of the present constructions.

Viscosity of dilute suspensions of rodlike particles: A numerical simulation method

Science.gov (United States)

Yamamoto, Satoru; Matsuoka, Takaaki

1994-02-01

The recently developed simulation method, named as the particle simulation method (PSM), is extended to predict the viscosity of dilute suspensions of rodlike particles. In this method a rodlike particle is modeled by bonded spheres. Each bond has three types of springs for stretching, bending, and twisting deformation. The rod model can therefore deform by changing the bond distance, bond angle, and torsion angle between paired spheres. The rod model can represent a variety of rigidity by modifying the bond parameters related to Young's modulus and the shear modulus of the real particle. The time evolution of each constituent sphere of the rod model is followed by molecular-dynamics-type approach. The intrinsic viscosity of a suspension of rodlike particles is derived from calculating an increased energy dissipation for each sphere of the rod model in a viscous fluid. With and without deformation of the particle, the motion of the rodlike particle was numerically simulated in a three-dimensional simple shear flow at a low particle Reynolds number and without Brownian motion of particles. The intrinsic viscosity of the suspension of rodlike particles was investigated on orientation angle, rotation orbit, deformation, and aspect ratio of the particle. For the rigid rodlike particle, the simulated rotation orbit compared extremely well with theoretical one which was obtained for a rigid ellipsoidal particle by use of Jeffery's equation. The simulated dependence of the intrinsic viscosity on various factors was also identical with that of theories for suspensions of rigid rodlike particles. For the flexible rodlike particle, the rotation orbit could be obtained by the particle simulation method and it was also cleared that the intrinsic viscosity decreased as occurring of recoverable deformation of the rodlike particle induced by flow.

Fractional variational problems and particle in cell gyrokinetic simulations with fuzzy logic approach for tokamaks

Directory of Open Access Journals (Sweden)

Rastović Danilo

2009-01-01

Full Text Available In earlier Rastovic's papers [1] and [2], the effort was given to analyze the stochastic control of tokamaks. In this paper, the deterministic control of tokamak turbulence is investigated via fractional variational calculus, particle in cell simulations, and fuzzy logic methods. Fractional integrals can be considered as approximations of integrals on fractals. The turbulent media could be of the fractal structure and the corresponding equations should be changed to include the fractal features of the media.

Transient Simulation of Accumulating Particle Deposition in Pipe Flow

Science.gov (United States)

Hewett, James; Sellier, Mathieu

2015-11-01

Colloidal particles that deposit in pipe systems can lead to fouling which is an expensive problem in both the geothermal and oil & gas industries. We investigate the gradual accumulation of deposited colloids in pipe flow using numerical simulations. An Euler-Lagrangian approach is employed for modelling the fluid and particle phases. Particle transport to the pipe wall is modelled with Brownian motion and turbulent diffusion. A two-way coupling exists between the fouled material and the pipe flow; the local mass flux of depositing particles is affected by the surrounding fluid in the near-wall region. This coupling is modelled by changing the cells from fluid to solid as the deposited particles exceed each local cell volume. A similar method has been used to model fouling in engine exhaust systems (Paz et al., Heat Transfer Eng., 34(8-9):674-682, 2013). We compare our deposition velocities and deposition profiles with an experiment on silica scaling in turbulent pipe flow (Kokhanenko et al., 19th AFMC, 2014).

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

Optimization of the particle pusher in a diode simulation code

International Nuclear Information System (INIS)

Theimer, M.M.; Quintenz, J.P.

1979-09-01

The particle pusher in Sandia's particle-in-cell diode simulation code has been rewritten to reduce the required run time of a typical simulation. The resulting new version of the code has been found to run up to three times as fast as the original with comparable accuracy. The cost of this optimization was an increase in storage requirements of about 15%. The new version has also been written to run efficiently on a CRAY-1 computing system. Steps taken to affect this reduced run time are described. Various test cases are detailed

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.

PRIAM: A self consistent finite element code for particle simulation in electromagnetic fields

International Nuclear Information System (INIS)

Le Meur, G.; Touze, F.

1990-06-01

A 2 1/2 dimensional, relativistic particle simulation code is described. A short review of the used mixed finite element method is given. The treatment of the driving terms (charge and current densities), initial, boundary conditions are exposed. Graphical results are shown

Three dimensional simulated modelling of diffusion capacitance of ...

African Journals Online (AJOL)

A three dimensional (3-D) simulated modelling was developed to analyse the excess minority carrier density in the base of a polycrystalline bifacial silicon solar cell. The concept of junction recombination velocity was ado-pted to quantify carrier flow through the junction, and to examine the solar cell diffusion capacitance for ...

Two-dimensional Nonlinear Simulations of Temperature-anisotropy Instabilities with a Proton-alpha Drift

Science.gov (United States)

Markovskii, S. A.; Chandran, Benjamin D. G.; Vasquez, Bernard J.

2018-04-01

We present two-dimensional hybrid simulations of proton-cyclotron and mirror instabilities in a proton-alpha plasma with particle-in-cell ions and a neutralizing electron fluid. The instabilities are driven by the protons with temperature perpendicular to the background magnetic field larger than the parallel temperature. The alpha particles with initially isotropic temperature have a nonzero drift speed with respect to the protons. The minor ions are known to influence the relative effect of the proton-cyclotron and mirror instabilities. In this paper, we show that the mirror mode can dominate the power spectrum at the nonlinear stage even if its linear growth rate is significantly lower than that of the proton-cyclotron mode. The proton-cyclotron instability combined with the alpha-proton drift is a possible cause of the nonzero magnetic helicity observed in the solar wind for fluctuations propagating nearly parallel to the magnetic field. Our simulations generally confirm this concept but reveal a complex helicity spectrum that is not anticipated from the linear theory of the instability.

Simulation of bluff-body flows using iterative penalization in a multiresolution particle-mesh vortex method

DEFF Research Database (Denmark)

Spietz, Henrik Juul; Hejlesen, Mads Mølholm; Walther, Jens Honore

in the oncoming flow. This may lead to structural instability e.g. when the shedding frequency aligns with the natural frequency of the structure. Fluid structure interaction must especially be considered when designing long span bridges. A three dimensional vortex-in-cell method is applied for the direct......The ability to predict aerodynamic forces, due to the interaction of a fluid flow with a solid body, is central in many fields of engineering and is necessary to identify error-prone structural designs. In bluff-body flows the aerodynamic forces oscillate due to vortex shedding and variations...... numerical simulation of the flow past a bodies of arbitrary shape. Vortex methods use a simple formulation where only the trajectories of discrete vortex particles are simulated. The Lagrangian formulation eliminates the CFL type condition that Eulerian methods have to satisfy. This allows vortex methods...

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

Thermal Field Analysis and Simulation of an Infrared Belt Furnace Used for Solar Cells

Directory of Open Access Journals (Sweden)

Bai Lu

2014-01-01

Full Text Available During solar cell firing, volatile organic compounds (VOC and a small number of metal particles were removed using the gas flow. When the gas flow was disturbed by the thermal field of infrared belt furnace and structure, the metal particles in the discharging gas flow randomly adhered to the surface of solar cell, possibly causing contamination. Meanwhile, the gas flow also affected the thermal uniformity of the solar cell. In this paper, the heating mechanism of the solar cell caused by radiation, convection, and conduction during firing was analyzed. Afterward, four 2-dimensional (2D models of the furnace were proposed. The transient thermal fields with different gas inlets, outlets, and internal structures were simulated. The thermal fields and the temperature of the solar cell could remain stable and uniform when the gas outlets were installed at the ends and in the middle of the furnace, with the gas inlets being distributed evenly. To verify the results, we produced four types of furnaces according to the four simulated results. The experimental results indicated that the thermal distribution of the furnace and the characteristics of the solar cells were consistent with the simulation. These experiments improved the efficiency of the solar cells while optimizing the solar cell manufacturing equipment.

Graphics processing unit accelerated three-dimensional model for the simulation of pulsed low-temperature plasmas

Energy Technology Data Exchange (ETDEWEB)

Fierro, Andrew, E-mail: andrew.fierro@ttu.edu; Dickens, James; Neuber, Andreas [Center for Pulsed Power and Power Electronics, Department of Electrical and Computer Engineering, Texas Tech University, Lubbock, Texas 79409 (United States)

2014-12-15

A 3-dimensional particle-in-cell/Monte Carlo collision simulation that is fully implemented on a graphics processing unit (GPU) is described and used to determine low-temperature plasma characteristics at high reduced electric field, E/n, in nitrogen gas. Details of implementation on the GPU using the NVIDIA Compute Unified Device Architecture framework are discussed with respect to efficient code execution. The software is capable of tracking around 10 × 10{sup 6} particles with dynamic weighting and a total mesh size larger than 10{sup 8} cells. Verification of the simulation is performed by comparing the electron energy distribution function and plasma transport parameters to known Boltzmann Equation (BE) solvers. Under the assumption of a uniform electric field and neglecting the build-up of positive ion space charge, the simulation agrees well with the BE solvers. The model is utilized to calculate plasma characteristics of a pulsed, parallel plate discharge. A photoionization model provides the simulation with additional electrons after the initial seeded electron density has drifted towards the anode. Comparison of the performance benefits between the GPU-implementation versus a CPU-implementation is considered, and a speed-up factor of 13 for a 3D relaxation Poisson solver is obtained. Furthermore, a factor 60 speed-up is realized for parallelization of the electron processes.

Study on Characteristics of Constricted DC Plasma Using Particle-In-Cell Simulator

International Nuclear Information System (INIS)

Jo, Jong Gap; Park, Yeong Shin; Hwang, Yong Seok

2010-01-01

In dc glow discharge, when anode size is smaller than cathode, very small and bright plasma ball occurs in front of anode. This plasma is called constricted dc plasma and characterized by a high plasma density in positive glow, so called plasma ball, compared to the conventional dc plasma. For the reason, this plasma is utilized to ion or electron beam sources since the beam currents are enhanced by the dense anode glow. However, correlations between characteristics of the plasma (plasma density, electron temperature and space potential) and discharge conditions (anode size, discharge voltage, discharge current, pressure) have been a little investigated definitely clear in previous study because of the trouble of a diagnosis. The plasma ball which is the most essential part of the constricted plasma is too small to diagnose precisely without disturbing plasma. Therefore, we tried to analyze the constricted plasma through computer simulation with Particle-In-Cell (PIC) code. In this study, simulation result of constricted dc plasma as well as conventional dc glow discharge will be addressed and compared with each others

Three-dimensional simulation of flow and combustion for pulverised coal injection

Energy Technology Data Exchange (ETDEWEB)

Guo, B.Y.; Zulli, P.; Rogers, H.; Mathieson, J.G.; Yu, A.B. [BlueScope Steel Research, Port Kembla, NSW (Australia)

2005-07-01

A three-dimensional numerical model of pulverised coal injection has been developed for simulating coal flow and combustion in the tuyere and raceway of a blast furnace. The model has been used to simulate previously reported combustion tests, which feature an inclined co-axial lance with an annular cooling gas. The predicted coal burnout agrees well with that measured for three coals with volatile contents and particle size ranging between 20.2-36.4% and particle sizes 1-200 {mu}m. Many important phenomena including flow asymmetry, recirculating flow and particle dispersion in the combustion chamber have been predicted. The current model can reproduce the experimental observations including the effects on burnout of coal flowrate and the introduction of methane for lance cooling.

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...

Self-assembly of polyhedral metal–organic framework particles into three-dimensional ordered superstructures

NARCIS (Netherlands)

Avci, Civan; Imaz, Inhar; Carné-Sánchez, Arnau; Pariente, Jose Angel; Tasios, Nikos; Pérez-Carvajal, Javier; Alonso, Maria Isabel; Blanco, Alvaro; Dijkstra, M.; López, Cefe; Maspoch, Daniel

Self-assembly of particles into long-range, three-dimensional, ordered superstructures is crucial for the design of a variety of materials, including plasmonic sensing materials, energy or gas storage systems, catalysts and photonic crystals. Here, we have combined experimental and simulation data

Development of a particle method of characteristics (PMOC) for one-dimensional shock waves

Science.gov (United States)

Hwang, Y.-H.

2018-03-01

In the present study, a particle method of characteristics is put forward to simulate the evolution of one-dimensional shock waves in barotropic gaseous, closed-conduit, open-channel, and two-phase flows. All these flow phenomena can be described with the same set of governing equations. The proposed scheme is established based on the characteristic equations and formulated by assigning the computational particles to move along the characteristic curves. Both the right- and left-running characteristics are traced and represented by their associated computational particles. It inherits the computational merits from the conventional method of characteristics (MOC) and moving particle method, but without their individual deficiencies. In addition, special particles with dual states deduced to the enforcement of the Rankine-Hugoniot relation are deliberately imposed to emulate the shock structure. Numerical tests are carried out by solving some benchmark problems, and the computational results are compared with available analytical solutions. From the derivation procedure and obtained computational results, it is concluded that the proposed PMOC will be a useful tool to replicate one-dimensional shock waves.

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure Induced by HZE Particles

Science.gov (United States)

Ponomarev, Artem; Plante, Ianik; George, Kerry; Wu, Honglu

2014-01-01

The formation of double-strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte-Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high-energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation-induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex-type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05-0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta-rays becomes more important.

A Coulomb collision algorithm for weighted particle simulations

Science.gov (United States)

Miller, Ronald H.; Combi, Michael R.

1994-01-01

A binary Coulomb collision algorithm is developed for weighted particle simulations employing Monte Carlo techniques. Charged particles within a given spatial grid cell are pair-wise scattered, explicitly conserving momentum and implicitly conserving energy. A similar algorithm developed by Takizuka and Abe (1977) conserves momentum and energy provided the particles are unweighted (each particle representing equal fractions of the total particle density). If applied as is to simulations incorporating weighted particles, the plasma temperatures equilibrate to an incorrect temperature, as compared to theory. Using the appropriate pairing statistics, a Coulomb collision algorithm is developed for weighted particles. The algorithm conserves energy and momentum and produces the appropriate relaxation time scales as compared to theoretical predictions. Such an algorithm is necessary for future work studying self-consistent multi-species kinetic transport.

Particle-in-cell Simulations of Raman Laser Amplification in Preformed Plasmas

International Nuclear Information System (INIS)

Clark, Daniel S.; Fisch, Nathaniel J.

2003-01-01

Two critical issues in the amplification of laser pulses by backward Raman scattering in plasma slabs are the saturation mechanism of the amplification effect (which determines the maximum attainable output intensity of a Raman amplifier) and the optimal plasma density for amplification. Previous investigations [V.M. Malkin, et al., Phys. Rev. Lett., 82 (22):4448-4451, 1999] identified forward Raman scattering and modulational instabilities of the amplifying seed as the likely saturation mechanisms and lead to an estimated unfocused output intensities of 10 17 W/cm 2 . The optimal density for amplification is determined by the competing constraints of minimizing the plasma density so as to minimize the growth rate of the instabilities leading to saturation but also maintaining the plasma sufficiently dense that the driven Langmuir wave responsible for backscattering does not break prematurely. Here, particle-in-cell code are simulations presented which verify that saturation of backward Raman amplification does occur at intensities of ∼10 17 W/cm 2 by forward Raman scattering and modulational instabilities. The optimal density for amplification in a plasma with the representative temperature of T(sub)e = 200 eV is also shown in these simulations to be intermediate between the cold plasma wave-breaking density and the density limit found by assuming a water bag electron distribution function

Simulation of the relationship between porosity and tortuosity in porous media with cubic particles

International Nuclear Information System (INIS)

Tang Xiao-Wu; Sun Zu-Feng; Cheng Guan-Chu

2012-01-01

Tortuosity is an important parameter used in areas such as vascular medicine, neurobiology, and the field of soil permeability and diffusion to express the mass transport in porous media. It is a function of the porosity and the shape and distribution of particles. In this paper, the tortuosity of cubic particles is calculated. With the assumption that the porous medium is homogeneous, the problem is converted to the micro-level over a unit cell, and geometry models of flow paths are proposed. In three-dimensional (3D) cells, the flow paths are too complicated to define. Hence, the 3D models are converted to two-dimensional (2D) models to simplify the calculation process. It is noticed that the path in the 2D model is shorter than that in the 3D model. As a result, triangular particles and the interaction are also taken into consideration to account for the longer distance respectively. We have proposed quadrate particle and interaction (QI) and quadrate and triangular particle (QT) models with cubic particles. Both models have shown good agreement with the experimental data. It is also found that they can predict the toruosities of some kinds of porous media, like freshwater sediment and Negev chalk

Two-dimensional numerical simulation of boron diffusion for pyramidally textured silicon

International Nuclear Information System (INIS)

Ma, Fa-Jun; Duttagupta, Shubham; Shetty, Kishan Devappa; Meng, Lei; Hoex, Bram; Peters, Ian Marius; Samudra, Ganesh S.

2014-01-01

Multidimensional numerical simulation of boron diffusion is of great relevance for the improvement of industrial n-type crystalline silicon wafer solar cells. However, surface passivation of boron diffused area is typically studied in one dimension on planar lifetime samples. This approach neglects the effects of the solar cell pyramidal texture on the boron doping process and resulting doping profile. In this work, we present a theoretical study using a two-dimensional surface morphology for pyramidally textured samples. The boron diffusivity and segregation coefficient between oxide and silicon in simulation are determined by reproducing measured one-dimensional boron depth profiles prepared using different boron diffusion recipes on planar samples. The established parameters are subsequently used to simulate the boron diffusion process on textured samples. The simulated junction depth is found to agree quantitatively well with electron beam induced current measurements. Finally, chemical passivation on planar and textured samples is compared in device simulation. Particularly, a two-dimensional approach is adopted for textured samples to evaluate chemical passivation. The intrinsic emitter saturation current density, which is only related to Auger and radiative recombination, is also simulated for both planar and textured samples. The differences between planar and textured samples are discussed

Cell structures caused by settling particles in turbulent Rayleigh-Bénard convection

Science.gov (United States)

Lee, Changhoon; Park, Sangro

2016-11-01

Turbulent thermal convection is an important phenomenon frequently found in nature and industrial processes, often with laden particles. In the last several decades, the vast majority of studies have addressed single phase convective flow with focus on the scaling relation of flow parameters associated with heat transfer. Particle-laden Rayleigh-Bénard convection, however, has not been sufficiently studied. In this study, modulation of cell structures by settling particles in turbulent Rayleigh-Bénard convection in a doubly periodic square channel is investigated using direct numerical simulation with a point particle approach. Flow parameters are fixed at Rayleigh number=106, Prandtl number=0.7, the aspect ratio=6, and Froude number=0.19. We report from the simulations that settling heavy particles modulate irregular large-scale thermal plume structures into organized polygonal cell structures. Different shapes of flow structures are obtained for different particle diameters and mass loadings. We found that polygonal cell structures arise due to asymmetric feedback force exerted by particles onto hot and cold plumes. Increasing the number of particles augments the asymmetry and the polygonal cell structures become smaller, eventually going to the hexagonal structures.

Three-dimensional Virtual Simulation of Oil Spill of Yangtze River in Chongqing Area Based on Emergency Decision

Science.gov (United States)

Chen, Shuzhe; Huang, Liwen

the river of Yangtze River in Chongqing area is continuous curved. Hydrology and channel situation is complex, and the transportation is busy. With the increasing of shipments of hazardous chemicals year by year, oil spill accident risk is rising. So establishment of three-dimensional virtual simulation of oil spill and its application in decision-making has become an urgent task. This paper detailed the process of three-dimensional virtual simulation of oil spill and established a system of three-dimensional virtual Simulation of oil spill of Yangtze River in Chongqing area by establishing an oil spill model of the Chongqing area based on oil particles model, and the system has been used in emergency decision to provide assistance for the oil spill response.

Implementation of a 3D plasma particle-in-cell code on a MIMD parallel computer

International Nuclear Information System (INIS)

Liewer, P.C.; Lyster, P.; Wang, J.

1993-01-01

A three-dimensional plasma particle-in-cell (PIC) code has been implemented on the Intel Delta MIMD parallel supercomputer using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors: A processor is assigned a subdomain and all the particles in it. Particles must be exchanged between processors as they move. Results are presented comparing the efficiency for 1-, 2- and 3-dimensional partitions of the three dimensional domain. This algorithm has been found to be very efficient even when a large fraction (e.g. 30%) of the particles must be exchanged at every time step. On the 512-node Intel Delta, up to 125 million particles have been pushed with an electrostatic push time of under 500 nsec/particle/time step

Time-dependent transport of energetic particles in magnetic turbulence: computer simulations versus analytical theory

Science.gov (United States)

Arendt, V.; Shalchi, A.

2018-06-01

We explore numerically the transport of energetic particles in a turbulent magnetic field configuration. A test-particle code is employed to compute running diffusion coefficients as well as particle distribution functions in the different directions of space. Our numerical findings are compared with models commonly used in diffusion theory such as Gaussian distribution functions and solutions of the cosmic ray Fokker-Planck equation. Furthermore, we compare the running diffusion coefficients across the mean magnetic field with solutions obtained from the time-dependent version of the unified non-linear transport theory. In most cases we find that particle distribution functions are indeed of Gaussian form as long as a two-component turbulence model is employed. For turbulence setups with reduced dimensionality, however, the Gaussian distribution can no longer be obtained. It is also shown that the unified non-linear transport theory agrees with simulated perpendicular diffusion coefficients as long as the pure two-dimensional model is excluded.

Dynamic Simulation of Random Packing of Polydispersive Fine Particles

Science.gov (United States)

Ferraz, Carlos Handrey Araujo; Marques, Samuel Apolinário

2018-02-01

In this paper, we perform molecular dynamic (MD) simulations to study the two-dimensional packing process of both monosized and random size particles with radii ranging from 1.0 to 7.0 μm. The initial positions as well as the radii of five thousand fine particles were defined inside a rectangular box by using a random number generator. Both the translational and rotational movements of each particle were considered in the simulations. In order to deal with interacting fine particles, we take into account both the contact forces and the long-range dispersive forces. We account for normal and static/sliding tangential friction forces between particles and between particle and wall by means of a linear model approach, while the long-range dispersive forces are computed by using a Lennard-Jones-like potential. The packing processes were studied assuming different long-range interaction strengths. We carry out statistical calculations of the different quantities studied such as packing density, mean coordination number, kinetic energy, and radial distribution function as the system evolves over time. We find that the long-range dispersive forces can strongly influence the packing process dynamics as they might form large particle clusters, depending on the intensity of the long-range interaction strength.

Simulation of radiation effects on three-dimensional computer optical memories

Science.gov (United States)

Moscovitch, M.; Emfietzoglou, D.

1997-01-01

A model was developed to simulate the effects of heavy charged-particle (HCP) radiation on the information stored in three-dimensional computer optical memories. The model is based on (i) the HCP track radial dose distribution, (ii) the spatial and temporal distribution of temperature in the track, (iii) the matrix-specific radiation-induced changes that will affect the response, and (iv) the kinetics of transition of photochromic molecules from the colored to the colorless isomeric form (bit flip). It is shown that information stored in a volume of several nanometers radius around the particle's track axis may be lost. The magnitude of the effect is dependent on the particle's track structure.

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

Self-consistent Analysis of Three-dimensional Uniformly Charged Ellipsoid with Zero Emittance

International Nuclear Information System (INIS)

Batygin, Yuri K.

2001-01-01

A self-consistent treatment of a three-dimensional ellipsoid with negligible emittance in time-dependent external field is performed. Envelope equations describing the evolution of an ellipsoid boundary are discussed. For a complete model it is required that the initial particle momenta be a linear function of the coordinates. Numerical example and verification of the problem by a 3-dimensional particle-in-cell simulations are given

Unit-cell design for two-dimensional phase-field simulation of microstructure evolution in single-crystal Ni-based superalloys during solidification

Directory of Open Access Journals (Sweden)

Dongjia Cao

2017-12-01

Full Text Available Phase-field simulation serves as an effective tool for quantitative characterization of microstructure evolution in single-crystal Ni-based superalloys during solidification nowadays. The classic unit cell is either limited to γ dendrites along crystal orientation or too ideal to cover complex morphologies for γ dendrites. An attempt to design the unit cell for two-dimensional (2-D phase-field simulations of microstructure evolution in single-crystal Ni-based superalloys during solidification was thus performed by using the MICRESS (MICRostructure Evolution Simulation Software in the framework of the multi-phase-field (MPF model, and demonstrated in a commercial TMS-113 superalloy. The coupling to CALPHAD (CALculation of PHAse Diagram thermodynamic database was realized via the TQ interface and the experimental diffusion coefficients were utilized in the simulation. Firstly, the classic unit cell with a single γ dendrite along crystal orientation was employed for the phase-field simulation in order to reproduce the microstructure features. Then, such simple unit cell was extended into the cases with two other different crystal orientations, i.e., and . Thirdly, for crystal orientations, the effect of γ dendritic orientations and unit cell sizes on microstructure and microsegregation was comprehensively studied, from which a new unit cell with multiple γ dendrites was proposed. The phase-field simulation with the newly proposed unit cell was further performed in the TMS-113 superalloy, and the microstructure features including the competitive growth of γ dendrites, microsegregation of different solutes and distribution of γ′ grains, can be nicely reproduced.

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure iIduced by HZE Particles

Science.gov (United States)

Ponomarev, Artem; Plante, Ianik; George, Kerry; Wu, Honglu

2014-01-01

The formation of double-strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte-Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high-energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation-induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex-type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05-0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta-rays becomes more important.

Particle-in-cell Simulation of Dipolarization Front Associated Whistlers

Science.gov (United States)

Lin, D.; Scales, W.; Ganguli, G.; Crabtree, C. E.

2017-12-01

Dipolarization fronts (DFs) are dipolarized magnetic field embedded in the Earthward propagating bursty bulk flows (BBFs), which separates the hot, tenuous high-speed flow from the cold, dense, and slowly convecting surrounding plasma [Runov et al. 2011]. Broadband fluctuations have been observed at DFs including the electromagnetic whistler waves and electrostatic lower hybrid waves in the Very Low Frequency (VLF) range [e.g., Zhou et al. 2009, Deng et al. 2010]. There waves are suggested to be able heat electrons and play a critical role in the plasma sheet dynamics [Chaston et al., 2012, Angelopoulos et al., 2013]. However, their generation mechanism and role in the energy conversion are still under debate. The gradient scale of magnetic field, plasma density at DFs in the near-Earth magnetotail is comparable to or lower than the ion gyro radius [Runov et al., 2011, Fu et al., 2012, Breuillard et al., 2016]. Such strongly inhomogeneous configuration could be unstable to the electron-ion hybrid (EIH) instability, which arises from strongly sheared transverse flow and is in the VLF range [Ganguli et al. 1988, Ganguli et al. 2014]. The equilibrium of the EIH theory implies an anisotropy of electron temperature, which are likely to drive the whistler waves observed in DFs [Deng et al., 2010, Gary et al., 2011]. In order to better understand how the whistler waves are generated in DFs and whether the EIH theory is applicable, a fully electromagnetic particle-in-cell (EMPIC) model is used to simulate the EIH instability with similar equilibrium configurations in DF observations. The EMPIC model deals with three dimensions in the velocity space and two dimensions in the configuration space, which is quite ready to include the third configuration dimension. Simulation results will be shown in this presentation.

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.

Gyrokinetic particle-in-cell simulations of plasma microturbulence on advanced computing platforms

International Nuclear Information System (INIS)

Ethier, S; Tang, W M; Lin, Z

2005-01-01

Since its introduction in the early 1980s, the gyrokinetic particle-in-cell (PIC) method has been very successfully applied to the exploration of many important kinetic stability issues in magnetically confined plasmas. Its self-consistent treatment of charged particles and the associated electromagnetic fluctuations makes this method appropriate for studying enhanced transport driven by plasma turbulence. Advances in algorithms and computer hardware have led to the development of a parallel, global, gyrokinetic code in full toroidal geometry, the gyrokinetic toroidal code (GTC), developed at the Princeton Plasma Physics Laboratory. It has proven to be an invaluable tool to study key effects of low-frequency microturbulence in fusion plasmas. As a high-performance computing applications code, its flexible mixed-model parallel algorithm has allowed GTC to scale to over a thousand processors, which is routinely used for simulations. Improvements are continuously being made. As the US ramps up its support for the International Tokamak Experimental Reactor (ITER), the need for understanding the impact of turbulent transport in burning plasma fusion devices is of utmost importance. Accordingly, the GTC code is at the forefront of the set of numerical tools being used to assess and predict the performance of ITER on critical issues such as the efficiency of energy confinement in reactors

Computer simulation of phenomena in plasma via particles

International Nuclear Information System (INIS)

Alves, M.V.; Bittencourt, J.A.

1988-06-01

The method of plasma computer simulation via particles has become an efficient tool to investigate the time and spatial evolution of various physical phenomena in plasmas. This method is based on the study of the individual plasma particle motions interacting with one another and with the externally applied fields. Although fairly simple, it allows a non-linear analysis of complex plasma physical phenomena and to obtain diagnostics even for regions of the system where experimental measurements would be difficult to make. In this report, a general view of the electrostatic one-dimensional computer code ES1, originally developed by A. Bruce Langdon, is presented. The main mathematical artifice in this code is the use of a spatial grid in which various plasma particles are represented by ''superparticles'', using a given shape function. The principal characteristics of the model, the approximations made and the mathematical methods used to solve the equations involved, are described. The specification of the input parameters which characterize the system, the initial conditions and the graphic diagnostics which can be utilized, are also described. Results are presented illustrating graphically the behavior of the plasma oscillations, the two-stream instability and the beam-plasma instability. (author) [pt

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.

Load-balancing techniques for a parallel electromagnetic particle-in-cell code

Energy Technology Data Exchange (ETDEWEB)

PLIMPTON,STEVEN J.; SEIDEL,DAVID B.; PASIK,MICHAEL F.; COATS,REBECCA S.

2000-01-01

QUICKSILVER is a 3-d electromagnetic particle-in-cell simulation code developed and used at Sandia to model relativistic charged particle transport. It models the time-response of electromagnetic fields and low-density-plasmas in a self-consistent manner: the fields push the plasma particles and the plasma current modifies the fields. Through an LDRD project a new parallel version of QUICKSILVER was created to enable large-scale plasma simulations to be run on massively-parallel distributed-memory supercomputers with thousands of processors, such as the Intel Tflops and DEC CPlant machines at Sandia. The new parallel code implements nearly all the features of the original serial QUICKSILVER and can be run on any platform which supports the message-passing interface (MPI) standard as well as on single-processor workstations. This report describes basic strategies useful for parallelizing and load-balancing particle-in-cell codes, outlines the parallel algorithms used in this implementation, and provides a summary of the modifications made to QUICKSILVER. It also highlights a series of benchmark simulations which have been run with the new code that illustrate its performance and parallel efficiency. These calculations have up to a billion grid cells and particles and were run on thousands of processors. This report also serves as a user manual for people wishing to run parallel QUICKSILVER.

Load-balancing techniques for a parallel electromagnetic particle-in-cell code

International Nuclear Information System (INIS)

Plimpton, Steven J.; Seidel, David B.; Pasik, Michael F.; Coats, Rebecca S.

2000-01-01

QUICKSILVER is a 3-d electromagnetic particle-in-cell simulation code developed and used at Sandia to model relativistic charged particle transport. It models the time-response of electromagnetic fields and low-density-plasmas in a self-consistent manner: the fields push the plasma particles and the plasma current modifies the fields. Through an LDRD project a new parallel version of QUICKSILVER was created to enable large-scale plasma simulations to be run on massively-parallel distributed-memory supercomputers with thousands of processors, such as the Intel Tflops and DEC CPlant machines at Sandia. The new parallel code implements nearly all the features of the original serial QUICKSILVER and can be run on any platform which supports the message-passing interface (MPI) standard as well as on single-processor workstations. This report describes basic strategies useful for parallelizing and load-balancing particle-in-cell codes, outlines the parallel algorithms used in this implementation, and provides a summary of the modifications made to QUICKSILVER. It also highlights a series of benchmark simulations which have been run with the new code that illustrate its performance and parallel efficiency. These calculations have up to a billion grid cells and particles and were run on thousands of processors. This report also serves as a user manual for people wishing to run parallel QUICKSILVER

Radiation in Particle Simulations

International Nuclear Information System (INIS)

More, R.; Graziani, F.; Glosli, J.; Surh, M.

2010-01-01

Hot dense radiative (HDR) plasmas common to Inertial Confinement Fusion (ICF) and stellar interiors have high temperature (a few hundred eV to tens of keV), high density (tens to hundreds of g/cc) and high pressure (hundreds of megabars to thousands of gigabars). Typically, such plasmas undergo collisional, radiative, atomic and possibly thermonuclear processes. In order to describe HDR plasmas, computational physicists in ICF and astrophysics use atomic-scale microphysical models implemented in various simulation codes. Experimental validation of the models used to describe HDR plasmas are difficult to perform. Direct Numerical Simulation (DNS) of the many-body interactions of plasmas is a promising approach to model validation but, previous work either relies on the collisionless approximation or ignores radiation. We present four methods that attempt a new numerical simulation technique to address a currently unsolved problem: the extension of molecular dynamics to collisional plasmas including emission and absorption of radiation. The first method applies the Lienard-Weichert solution of Maxwell's equations for a classical particle whose motion is assumed to be known. The second method expands the electromagnetic field in normal modes (planewaves in a box with periodic boundary-conditions) and solves the equation for wave amplitudes coupled to the particle motion. The third method is a hybrid molecular dynamics/Monte Carlo (MD/MC) method which calculates radiation emitted or absorbed by electron-ion pairs during close collisions. The fourth method is a generalization of the third method to include small clusters of particles emitting radiation during close encounters: one electron simultaneously hitting two ions, two electrons simultaneously hitting one ion, etc. This approach is inspired by the virial expansion method of equilibrium statistical mechanics. Using a combination of these methods we believe it is possible to do atomic-scale particle simulations of

Discrete Element Simulation of Elastoplastic Shock Wave Propagation in Spherical Particles

Directory of Open Access Journals (Sweden)

M. Shoaib

2011-01-01

Full Text Available Elastoplastic shock wave propagation in a one-dimensional assembly of spherical metal particles is presented by extending well-established quasistatic compaction models. The compaction process is modeled by a discrete element method while using elastic and plastic loading, elastic unloading, and adhesion at contacts with typical dynamic loading parameters. Of particular interest is to study the development of the elastoplastic shock wave, its propagation, and reflection during entire loading process. Simulation results yield information on contact behavior, velocity, and deformation of particles during dynamic loading. Effects of shock wave propagation on loading parameters are also discussed. The elastoplastic shock propagation in granular material has many practical applications including the high-velocity compaction of particulate material.

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

Simulation of dust-acoustic waves

International Nuclear Information System (INIS)

Winske, D.; Murillo, M.S.; Rosenberg, M.

1998-01-01

The authors use molecular dynamics (MD) and particle-in-cell (PIC) simulation methods to investigate the dispersion relation of dust-acoustic waves in a one-dimensional, strongly coupled (Coulomb coupling parameter Λ = ratio of the Coulomb energy to the thermal energy = 120) dusty plasma. They study both cases where the dust is represented by a small number of simulation particles that form into a regular array structure (crystal limit) as well as where the dust is represented by a much larger number of particles (fluid limit)

E × B electron drift instability in Hall thrusters: Particle-in-cell simulations vs. theory

Science.gov (United States)

Boeuf, J. P.; Garrigues, L.

2018-06-01

The E × B Electron Drift Instability (E × B EDI), also called Electron Cyclotron Drift Instability, has been observed in recent particle simulations of Hall thrusters and is a possible candidate to explain anomalous electron transport across the magnetic field in these devices. This instability is characterized by the development of an azimuthal wave with wavelength in the mm range and velocity on the order of the ion acoustic velocity, which enhances electron transport across the magnetic field. In this paper, we study the development and convection of the E × B EDI in the acceleration and near plume regions of a Hall thruster using a simplified 2D axial-azimuthal Particle-In-Cell simulation. The simulation is collisionless and the ionization profile is not-self-consistent but rather is given as an input parameter of the model. The aim is to study the development and properties of the instability for different values of the ionization rate (i.e., of the total ion production rate or current) and to compare the results with the theory. An important result is that the wavelength of the simulated azimuthal wave scales as the electron Debye length and that its frequency is on the order of the ion plasma frequency. This is consistent with the theory predicting destruction of electron cyclotron resonance of the E × B EDI in the non-linear regime resulting in the transition to an ion acoustic instability. The simulations also show that for plasma densities smaller than under nominal conditions of Hall thrusters the field fluctuations induced by the E × B EDI are no longer sufficient to significantly enhance electron transport across the magnetic field, and transit time instabilities develop in the axial direction. The conditions and results of the simulations are described in detail in this paper and they can serve as benchmarks for comparisons between different simulation codes. Such benchmarks would be very useful to study the role of numerical noise (numerical

Three-dimensional growth simulation: A study of substrate oriented films

International Nuclear Information System (INIS)

Besnard, A; Martin, N; Carpentier, L

2010-01-01

Monte Carlo simulations are developed to simulate the growth of three-dimensional columnar microstructure in thin films. We are studying in particular oriented microstructure like those produced with the Glancing Angle Deposition technique (GLAD). Some geometrical characteristics of the particles flux, the organization of defect sites on the substrate surface and the atomic surface diffusion are mainly investigated in order to predict the growth processes and the resulting features of the films. This study reports on simulations of thin film growth exhibiting an oblique and zigzag columnar microstructure. Column angle evolution and density are investigated versus incidence angle α or period number n and compared with experimental measurements.

Fully non-linear multi-species Fokker-Planck-Landau collisions for gyrokinetic particle-in-cell simulations of fusion plasma

Science.gov (United States)

Hager, Robert; Yoon, E. S.; Ku, S.; D'Azevedo, E. F.; Worley, P. H.; Chang, C. S.

2015-11-01

We describe the implementation, and application of a time-dependent, fully nonlinear multi-species Fokker-Planck-Landau collision operator based on the single-species work of Yoon and Chang [Phys. Plasmas 21, 032503 (2014)] in the full-function gyrokinetic particle-in-cell codes XGC1 [Ku et al., Nucl. Fusion 49, 115021 (2009)] and XGCa. XGC simulations include the pedestal and scrape-off layer, where significant deviations of the particle distribution function from a Maxwellian can occur. Thus, in order to describe collisional effects on neoclassical and turbulence physics accurately, the use of a non-linear collision operator is a necessity. Our collision operator is based on a finite volume method using the velocity-space distribution functions sampled from the marker particles. Since the same fine configuration space mesh is used for collisions and the Poisson solver, the workload due to collisions can be comparable to or larger than the workload due to particle motion. We demonstrate that computing time spent on collisions can be kept affordable by applying advanced parallelization strategies while conserving mass, momentum, and energy to reasonable accuracy. We also show results of production scale XGCa simulations in the H-mode pedestal and compare to conventional theory. Work supported by US DOE OFES and OASCR.

Two-dimensional simulation of gas concentration impedance for a planar solid oxide fuel cell

International Nuclear Information System (INIS)

Fadaei, M.; Mohammadi, R.; Ghassemi, M.

2014-01-01

Highlights: • The 2D simulation shows another feature in concentration impedance. • The channel gas transport causes a capacitive behavior. • Anode polarization variation has a significant influence on velocity distribution. • The influence of 2D simulation is important for channel height bigger than 2 mm. - Abstract: This paper presents a two-dimensional model for a planar solid oxide fuel cell (SOFC) anode in order to simulate the steady-state performance characteristics as well as the electrochemical impedance spectra. The developed model couples the mass transport with the electrochemical kinetics. The transient conservation equations (momentum and species equations) are solved numerically and the linear kinetic is used for the anode electrochemistry. In order to solve the system of the nonlinear equations, an in-house code based on the finite volume method is developed and utilized. A parametric study is also carried out and the results are discussed. Results show a capacitive semicircle in the Nyquist plot which is identical to the gas concentration impedance. The simulation results are in good agreement with published data

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 simulation for the effect of segmented electrodes near the exit of an aton-type Hall thruster on ion focusing acceleration

Energy Technology Data Exchange (ETDEWEB)

Yu, D.R.; Qing, S.W.; Liu, H.; Li, H. [Lab. of Plasma Propulsion, Harbin Institute of Technology (China)

2011-12-15

The effect of floating conductive electrodes near the channel exit of an Aton-type Hall thruster on ion focusing acceleration is studied by simulating the two-dimensional plasma flow with a fully kinetic Particle-in-Cell method for the gas flow rate j{sub a} ranged in 1{proportional_to}3 mg/s. Numerical results show that low-emissive electrodes can reduce plume divergence if the electrode length is less than 2 mm due to the low secondary electron emissive characteristic, but widen plume in all the gas flow rate range if the electrode length is greater than 2mm since the conductive property of segmented electrodes trends to make equipotential lines convex toward channel exit and is even parallel to the wall surface in the near-wall region. Further investigation predicts that the combination of high emissive dielectric wall and segmented low-emissive dielectric wall is a promising way to reduce plume divergence (copyright 2011 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim) (orig.)

Performance of nitrogen-doped graphene aerogel particle electrodes for electro-catalytic oxidation of simulated Bisphenol A wastewaters

Energy Technology Data Exchange (ETDEWEB)

Chen, Zhuang [Environmental Research Academy, North China Electric Power University, Beijing 102206 (China); Labortory of Environmental Remediation and Functional Material, Suzhou Research Academy of North China Electric Power University, Suzhou, Jiangsu, 215026 (China); Zhang, Yimei, E-mail: yimei.zhang1@gmail.com [Environmental Research Academy, North China Electric Power University, Beijing 102206 (China); Labortory of Environmental Remediation and Functional Material, Suzhou Research Academy of North China Electric Power University, Suzhou, Jiangsu, 215026 (China); Zhou, Lincheng; Zhu, Hao; Wan, Fei; Wang, Yue [Labortory of Environmental Remediation and Functional Material, Suzhou Research Academy of North China Electric Power University, Suzhou, Jiangsu, 215026 (China); Zhang, Dandan [Environmental Research Academy, North China Electric Power University, Beijing 102206 (China)

2017-06-15

Highlights: • The new three-dimensional electrode system with the nitrogen-doped graphene aerogels particle electrodes is developed. • The optimal technique parameters were explored. • Simulated BPA waste waters are effectively degraded. • Degradation pathway and intermediates are proposed. • Three-dimensional electrode system shows good and stable removal performance. - Abstract: The treatment of effluent containing Bisphenol A (BPA) was investigated experimentally using nitrogen-doped graphene aerogel (NGAs) as particle electrodes in a three-dimensional electrode reactor for the electrochemical treatment was studied. The effects of the cell voltage, pH, the ratio of NGAs mass to solution volume and repeated times on the removal efficiency were investigated. Compared with commercial carbon particle electrodes, the NGAs exhibited stronger activity to remove BPA simulated wastewater. For 15 mg L{sup −1} of BPA solution, the degradation rate of BPA exceeded 90% after treatment for only 30 min under the optimum conditions. The COD{sub Cr} removal rate of BPA was 85%. Moreover, in the process of reused 50 times, the degradation rate of BPA can be kept in more than 85%. The COD{sub Cr} removal rate was stable at about 73%. The intermediate products of electrochemical degradation of BPA were identified by liquid chromatography-mass spectrometry liquid chromatography (LC–MS), and a probable BPA degradation pathway was proposed. It was considered that ·OH radicals by water electrolysis could constantly attack the aromatic ring to form various intermediates such as hydroxylated-BPA, isopropylphenol, hydroquinone, phenol and butantetraol, maleic acid, oxalic acid. These compounds were eventually mineralized by electrolysis into CO{sub 2} and H{sub 2}O.

Deformation of Two-Dimensional Nonuniform-Membrane Red Blood Cells Simulated by a Lattice Boltzmann Model

International Nuclear Information System (INIS)

Hua-Bing, Li; Li, Jin; Bing, Qiu

2008-01-01

To study two-dimensional red blood cells deforming in a shear Bow with the membrane nonuniform on the rigidity and mass, the membrane is discretized into equilength segments. The fluid inside and outside the red blood cell is simulated by the D2Q9 lattice Boltzmann model and the hydrodynamic forces exerted on the membrane from the inner and outer of the red blood cell are calculated by a stress-integration method. Through the global deviation from the curvature of uniform-membrane, we find that when the membrane is nonuniform on the rigidity, the deviation first decreases with the time increases and implies that the terminal profile of the red blood cell is static. To a red blood cell with the mass nonuniform on the membrane, the deviation becomes more large, and the mass distribution affects the profile of the two sides of the flattened red blood cell in a shear flow. (fundamental areas of phenomenology(including applications))

Two-dimensional Simulations of Correlation Reflectometry in Fusion Plasmas

International Nuclear Information System (INIS)

Valeo, E.J.; Kramer, G.J.; Nazikian, R.

2001-01-01

A two-dimensional wave propagation code, developed specifically to simulate correlation reflectometry in large-scale fusion plasmas is described. The code makes use of separate computational methods in the vacuum, underdense and reflection regions of the plasma in order to obtain the high computational efficiency necessary for correlation analysis. Simulations of Tokamak Fusion Test Reactor (TFTR) plasma with internal transport barriers are presented and compared with one-dimensional full-wave simulations. It is shown that the two-dimensional simulations are remarkably similar to the results of the one-dimensional full-wave analysis for a wide range of turbulent correlation lengths. Implications for the interpretation of correlation reflectometer measurements in fusion plasma are discussed

Simulation of radiation effects on three-dimensional computer optical memories

International Nuclear Information System (INIS)

Moscovitch, M.; Emfietzoglou, D.

1997-01-01

A model was developed to simulate the effects of heavy charged-particle (HCP) radiation on the information stored in three-dimensional computer optical memories. The model is based on (i) the HCP track radial dose distribution, (ii) the spatial and temporal distribution of temperature in the track, (iii) the matrix-specific radiation-induced changes that will affect the response, and (iv) the kinetics of transition of photochromic molecules from the colored to the colorless isomeric form (bit flip). It is shown that information stored in a volume of several nanometers radius around the particle close-quote s track axis may be lost. The magnitude of the effect is dependent on the particle close-quote s track structure. copyright 1997 American Institute of Physics

Artificial cooling due to quiet injection in bounded plasma particle simulations

International Nuclear Information System (INIS)

Lawson, W.S.

1988-01-01

An explanation is proposed for an artificial cooling effect seen in electrostatic particle-in-cell plasma simulations. The effect hinges on heat transport from the trapped electrons to fluctuations of the electric field, which are kept at a sub-thermal level through the continuous ''quiet'' injection of passing electrons. Further simulations are done which test and support the explanation. copyright 1988 Academic Press, Inc

Simulation of neutron transport process, photons and charged particles within the Monte Carlo method

International Nuclear Information System (INIS)

Androsenko, A.A.; Androsenko, P.A.; Artamonov, S.N.; Bolonkina, G.V.; Lomtev, V.L.; Pupko, S.V.

1991-01-01

Description is given to the program system BRAND designed for the accurate solution of non-stationary transport equation of neutrons, photons and charged particles in the conditions of real three-dimensional geometry. An extensive set of local and non-local estimates provides an opportunity of calculating a great set of linear functionals normally being of interest in the calculation of reactors, radiation protection and experiment simulation. The process of particle interaction with substance is simulated on the basis of individual non-group data on each isotope of the composition. 24 refs

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.

Large eddy simulation modeling of particle-laden flows in complex terrain

Science.gov (United States)

Salesky, S.; Giometto, M. G.; Chamecki, M.; Lehning, M.; Parlange, M. B.

2017-12-01

The transport, deposition, and erosion of heavy particles over complex terrain in the atmospheric boundary layer is an important process for hydrology, air quality forecasting, biology, and geomorphology. However, in situ observations can be challenging in complex terrain due to spatial heterogeneity. Furthermore, there is a need to develop numerical tools that can accurately represent the physics of these multiphase flows over complex surfaces. We present a new numerical approach to accurately model the transport and deposition of heavy particles in complex terrain using large eddy simulation (LES). Particle transport is represented through solution of the advection-diffusion equation including terms that represent gravitational settling and inertia. The particle conservation equation is discretized in a cut-cell finite volume framework in order to accurately enforce mass conservation. Simulation results will be validated with experimental data, and numerical considerations required to enforce boundary conditions at the surface will be discussed. Applications will be presented in the context of snow deposition and transport, as well as urban dispersion.

Numerical simulation of two-dimensional late-stage coarsening for nucleation and growth

International Nuclear Information System (INIS)

Akaiwa, N.; Meiron, D.I.

1995-01-01

Numerical simulations of two-dimensional late-stage coarsening for nucleation and growth or Ostwald ripening are performed at area fractions 0.05 to 0.4 using the monopole and dipole approximations of a boundary integral formulation for the steady state diffusion equation. The simulations are performed using two different initial spatial distributions. One is a random spatial distribution, and the other is a random spatial distribution with depletion zones around the particles. We characterize the spatial correlations of particles by the radial distribution function, the pair correlation functions, and the structure function. Although the initial spatial correlations are different, we find time-independent scaled correlation functions in the late stage of coarsening. An important feature of the late-stage spatial correlations is that depletion zones exist around particles. A log-log plot of the structure function shows that the slope at small wave numbers is close to 4 and is -3 at very large wave numbers for all area fractions. At large wave numbers we observe oscillations in the structure function. We also confirm the cubic growth law of the average particle radius. The rate constant of the cubic growth law and the particle size distribution functions are also determined. We find qualitatively good agreement between experiments and the present simulations. In addition, the present results agree well with simulation results using the Cahn-Hilliard equation

Optical separation and controllable delivery of cells from particle and cell mixture

Directory of Open Access Journals (Sweden)

Li Yuchao

2015-11-01

Full Text Available Cell separation and delivery have recently gained significant attention in biological and biochemical studies. In thiswork, an optical method for separation and controllable delivery of cells by using an abruptly tapered fiber probe is reported. By launching a laser beam at the wavelength of 980 nm into the fiber, a mixture of cells with sizes of ~5 and ~3 μm and poly(methyl methacrylate particles with size of 5 μm are separated into three chains along the direction of propagation of light. The cell and particle chains are delivered in three dimensions over 600 μm distance. Experimental results are interpreted by numerical simulations. Optical forces and forward migration velocities of different particles and cells are calculated and discussed.

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

Collective ratchet effects and reversals for active matter particles on quasi-one-dimensional asymmetric substrates.

Science.gov (United States)

McDermott, Danielle; Olson Reichhardt, Cynthia J; Reichhardt, Charles

2016-10-19

Using computer simulations, we study a two-dimensional system of sterically interacting self-mobile run-and-tumble disk-shaped particles with an underlying periodic quasi-one-dimensional asymmetric substrate, and show that a rich variety of collective active ratchet behaviors arise as a function of particle density, activity, substrate period, and the maximum force exerted by the substrate. The net dc drift, or ratchet transport flux, is nonmonotonic since it increases with increased activity but is diminished by the onset of self-clustering of the active particles. Increasing the particle density decreases the ratchet transport flux for shallow substrates but increases the ratchet transport flux for deep substrates due to collective hopping events. At the highest particle densities, the ratchet motion is destroyed by a self-jamming effect. We show that it is possible to realize reversals of the direction of the net dc drift in the deep substrate limit when multiple rows of active particles can be confined in each substrate minimum, permitting emergent particle-like excitations to appear that experience an inverted effective substrate potential. We map out a phase diagram of the forward and reverse ratchet effects as a function of the particle density, activity, and substrate properties.

H5Part A Portable High Performance Parallel Data Interface for Particle Simulations

CERN Document Server

Adelmann, Andreas; Shalf, John M; Siegerist, Cristina

2005-01-01

Largest parallel particle simulations, in six dimensional phase space generate wast amont of data. It is also desirable to share data and data analysis tools such as ParViT (Particle Visualization Toolkit) among other groups who are working on particle-based accelerator simulations. We define a very simple file schema built on top of HDF5 (Hierarchical Data Format version 5) as well as an API that simplifies the reading/writing of the data to the HDF5 file format. HDF5 offers a self-describing machine-independent binary file format that supports scalable parallel I/O performance for MPI codes on a variety of supercomputing systems and works equally well on laptop computers. The API is available for C, C++, and Fortran codes. The file format will enable disparate research groups with very different simulation implementations to share data transparently and share data analysis tools. For instance, the common file format will enable groups that depend on completely different simulation implementations to share c...

Self-Assembly of DNA-Coated Particles: Experiment, Simulation and Theory

Science.gov (United States)

Song, Minseok

The bottom-up assembly of material architectures with tunable complexity, function, composition, and structure is a long sought goal in rational materials design. One promising approach aims to harnesses the programmability and specificity of DNA hybridization in order to direct the assembly of oligonucleotide-functionalized nano- and micro-particles by tailoring, in part, interparticle interactions. DNA-programmable assembly into three-dimensionally ordered structures has attracted extensive research interest owing to emergent applications in photonics, plasmonics and catalysis and potentially many other areas. Progress on the rational design of DNA-mediated interactions to create useful two-dimensional structures (e.g., structured films), on the other hand, has been rather slow. In this thesis, we establish strategies to engineer a diversity of 2D crystalline arrangements by designing and exploiting DNA-programmable interparticle interactions. We employ a combination of simulation, theory and experiments to predict and confirm accessibility of 2D structural diversity in an effort to establish a rational approach to 2D DNA-mediated particle assembly. We start with the experimental realization of 2D DNA-mediated assembly by decorating micron-sized silica particles with covalently attached single-stranded DNA through a two-step reaction. Subsequently, we elucidate sensitivity and ultimate controllability of DNA-mediated assembly---specifically the melting transition from dispersed singlet particles to aggregated or assembled structures---through control of the concentration of commonly employed nonionic surfactants. We relate the observed tunability to an apparent coupling with the critical micelle temperature in these systems. Also, both square and hexagonal 2D ordered particle arrangements are shown to evolve from disordered aggregates under appropriate annealing conditions defined based upon pre-established melting profiles. Subsequently, the controlled mixing of

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

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.

Particle-in-cell/accelerator code for space-charge dominated beam simulation

Energy Technology Data Exchange (ETDEWEB)

2012-05-08

Warp is a multidimensional discrete-particle beam simulation program designed to be applicable where the beam space-charge is non-negligible or dominant. It is being developed in a collaboration among LLNL, LBNL and the University of Maryland. It was originally designed and optimized for heave ion fusion accelerator physics studies, but has received use in a broader range of applications, including for example laser wakefield accelerators, e-cloud studies in high enery accelerators, particle traps and other areas. At present it incorporates 3-D, axisymmetric (r,z) planar (x-z) and transverse slice (x,y) descriptions, with both electrostatic and electro-magnetic fields, and a beam envelope model. The code is guilt atop the Python interpreter language.

Solar energetic particles: observational studies and magnetohydrodynamic simulation

International Nuclear Information System (INIS)

Masson, S.

2010-10-01

Solar activity manifests itself through highly dynamical events, such as flares and coronal mass ejections, which result in energy release by magnetic reconnection. This thesis focuses on two manifestations of this energy release: solar energetic particles and dynamics of magnetic reconnection. The first part of my work consists in the detailed temporal analysis of several electromagnetic signatures, produced by energetic particles in the solar atmosphere, with respect to the energetic particle flux at Earth. Using multi-instrument observations, I highlighted that particles can be accelerated by the flare to relativistic energies during a specific episode of acceleration in the impulsive phase. This showed that particles traveled a longer path length than the theoretical length generally assumed. Using in-situ measurements of magnetic field and plasma, I identified the interplanetary magnetic field for 10 particle events, and performing a velocity dispersion analysis I obtained the interplanetary length traveled by particles. I showed that the magnetic structure of the interplanetary medium play a crucial role in the association of the particle flux at Earth and the acceleration signatures of particles at the Sun. The second part of my work focuses on the dynamics of magnetic reconnection. Observationally, the best evidence for magnetic reconnection is the appearance of brightnesses at the solar surface. Performing the first data-driven 3 dimensional magneto-hydrodynamic (MHD) simulation of an observed event, I discovered that the evolution of brightnesses can be explained by the succession of two different reconnection regimes, induced by a new topological association where null-point separatrix lines are embedded in quasi-separatrix layers. This new topological association induces a change of field line connectivity, but also a continuous reconnection process, leading to an apparent slipping motion of reconnected field lines. From a MHD simulation I showed that

A Three-Dimensional Numerical Study of Gas-Particle Flow and Chemical Reactions in Circulating Fluidised Bed Reactors

DEFF Research Database (Denmark)

Hansen, Kim Granly

Three-dimensional Computational Fluid Dynamics (CFD) simulations of Circulating Fluidized Beds (CFB's) have been performed. The computations are performed using a 3D multiphase computational fluid dynamics code with an Eulerian description of both gas and particle phases. The turbulent motion...... implemented in the CFD code FLOTRACS-MP-3D. The decomposition reaction is studied in a 3D representation of a 0.254 m i.d. riser, which has been studied experimentally by Ouyang et al. (1993). Comparison between measured and simulated time-averaged ozone concentration at different elevations in the riser...

Three-dimensional random resistor-network model for solid oxide fuel cell composite electrodes

International Nuclear Information System (INIS)

Abbaspour, Ali; Luo Jingli; Nandakumar, K.

2010-01-01

A three-dimensional reconstruction of solid oxide fuel cell (SOFC) composite electrodes was developed to evaluate the performance and further investigate the effect of microstructure on the performance of SOFC electrodes. Porosity of the electrode is controlled by adding pore former particles (spheres) to the electrode and ignoring them in analysis step. To enhance connectivity between particles and increase the length of triple-phase boundary (TPB), sintering process is mimicked by enlarging particles to certain degree after settling them inside the packing. Geometrical characteristics such as length of TBP and active contact area as well as porosity can easily be calculated using the current model. Electrochemical process is simulated using resistor-network model and complete Butler-Volmer equation is used to deal with charge transfer process on TBP. The model shows that TPBs are not uniformly distributed across the electrode and location of TPBs as well as amount of electrochemical reaction is not uniform. Effects of electrode thickness, particle size ratio, electron and ion conductor conductivities and rate of electrochemical reaction on overall electrochemical performance of electrode are investigated.

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)

Three-dimensional simulation of large-scale structure in the universe

Energy Technology Data Exchange (ETDEWEB)

Centrella, J.; Melott, A.L.

1983-09-15

High and low density cloud-in-cell models were used to simulate the nonlinear growth of adiabatic perturbations in collisionless matter to demonstrate the development of a cellular structure in the universe. Account was taken of a short wvelength cutoff in collisionless matter, with a focus on resolving filaments and low density pancakes. The calculations were performed with a Friedmann-Robertson-Walker model, and the gravitational potential of dark matter was obtained through solution of the Poisson equation. The simulation began with z between 100-1000, and initial particle velocities were set at zero. Spherically symmetric voids were observed to form, then colide and interact. Sufficient particles were employed to avoid depletion during nonlinear collapse. No galaxies formed during the epoch studied, which has implications for the significance of dark, baryonic matter in the present universe.

A fast sorting algorithm for a hypersonic rarefied flow particle simulation on the connection machine

Science.gov (United States)

Dagum, Leonardo

1989-01-01

The data parallel implementation of a particle simulation for hypersonic rarefied flow described by Dagum associates a single parallel data element with each particle in the simulation. The simulated space is divided into discrete regions called cells containing a variable and constantly changing number of particles. The implementation requires a global sort of the parallel data elements so as to arrange them in an order that allows immediate access to the information associated with cells in the simulation. Described here is a very fast algorithm for performing the necessary ranking of the parallel data elements. The performance of the new algorithm is compared with that of the microcoded instruction for ranking on the Connection Machine.

Application of Simulated Three Dimensional CT Image in Orthognathic Surgery

Energy Technology Data Exchange (ETDEWEB)

Kim, Hyun Don; Park, Chang Seo [Dept. of Dental Radiology, College of Dentistry, Yensei University, Seoul (Korea, Republic of); Yoo, Sun Kook; Lee, Kyoung Sang [Dept. of Medical Engineering, College of Medicine, Yensei University, Seoul (Korea, Republic of)

1998-08-15

In orthodontics and orthognathic surgery, cephalogram has been routine practice in diagnosis and treatment evaluation of craniofacial deformity. But its inherent distortion of actual length and angles during projecting three dimensional object to two dimensional plane might cause errors in quantitative analysis of shape and size. Therefore, it is desirable that three dimensional object is diagnosed and evaluated three dimensionally and three dimensional CT image is best for three dimensional analysis. Development of clinic necessitates evaluation of result of treatment and comparison before and after surgery. It is desirable that patient that was diagnosed and planned by three dimensional computed tomography before surgery is evaluated by three dimensional computed tomography after surgery, too. But Because there is no standardized normal values in three dimension now and three dimensional Computed Tomography needs expensive equipment and because of its expenses and amount of exposure to radiation, limitations still remain to be solved in its application to routine practice. If postoperative three dimensional image is constructed by pre and postoperative lateral and postero-anterior cephalograms and preoperative three dimensional computed tomogram, pre and postoperative image will be compared and evaluated three dimensionally without three dimensional computed tomography after surgery and that will contribute to standardize normal values in three dimension. This study introduced new method that computer-simulated three dimensional image was constructed by preoperative three dimensional computed tomogram and pre and postoperative lateral and postero-anterior cephalograms, and for validation of new method, in four cases of dry skull that position of mandible was displaced and four patients of orthognathic surgery, computer-simulated three dimensional image and actual postoperative three dimensional image were compared. The results were as follows. 1. In four cases of

Application of Simulated Three Dimensional CT Image in Orthognathic Surgery

International Nuclear Information System (INIS)

Kim, Hyun Don; Park, Chang Seo; Yoo, Sun Kook; Lee, Kyoung Sang

1998-01-01

In orthodontics and orthognathic surgery, cephalogram has been routine practice in diagnosis and treatment evaluation of craniofacial deformity. But its inherent distortion of actual length and angles during projecting three dimensional object to two dimensional plane might cause errors in quantitative analysis of shape and size. Therefore, it is desirable that three dimensional object is diagnosed and evaluated three dimensionally and three dimensional CT image is best for three dimensional analysis. Development of clinic necessitates evaluation of result of treatment and comparison before and after surgery. It is desirable that patient that was diagnosed and planned by three dimensional computed tomography before surgery is evaluated by three dimensional computed tomography after surgery, too. But Because there is no standardized normal values in three dimension now and three dimensional Computed Tomography needs expensive equipment and because of its expenses and amount of exposure to radiation, limitations still remain to be solved in its application to routine practice. If postoperative three dimensional image is constructed by pre and postoperative lateral and postero-anterior cephalograms and preoperative three dimensional computed tomogram, pre and postoperative image will be compared and evaluated three dimensionally without three dimensional computed tomography after surgery and that will contribute to standardize normal values in three dimension. This study introduced new method that computer-simulated three dimensional image was constructed by preoperative three dimensional computed tomogram and pre and postoperative lateral and postero-anterior cephalograms, and for validation of new method, in four cases of dry skull that position of mandible was displaced and four patients of orthognathic surgery, computer-simulated three dimensional image and actual postoperative three dimensional image were compared. The results were as follows. 1. In four cases of

Particle-in-Cell Simulation Study on the Floating Potential of Spacecraft in the Low Earth Orbit

International Nuclear Information System (INIS)

Tang Daotan; Yang Shengsheng; Zheng Kuohai; Qin Xiaogang; Li Detian; Liu Qing; Zhao Chengxuan; Du Shanshan

2015-01-01

In order to further understand the characteristics of the floating potential of low earth orbit spacecraft, the effects of the electron current collection area, background electron temperature, photocurrent emission, spacecraft wake, and the shape of spacecraft on spacecraft floating potential were studied here by particle-in-cell simulation in the low earth orbit. The simulation results show that the electron current collection area and background electron temperature impact on the floating potential by changing the electron current collection of spacecraft. By increasing the electron current collection area or background electron temperature, the spacecraft will float at a lower electric potential with respect to the surrounding plasma. However, the spacecraft wake affects the floating potential by increasing the ion current collected by spacecraft. The emission of the photocurrent from the spacecraft surface, which compensates for the electrons collected from background plasma, causes the floating potential to increase. The shape of the spacecraft is also an important factor influencing the floating potential. (paper)

Three-dimensional simulation of the electromagnetic ion/ion beam instability: cross field diffusion

Directory of Open Access Journals (Sweden)

H. Kucharek

2000-01-01

Full Text Available In a system with at least one ignorable spatial dimension charged particles moving in fluctuating fields are tied to the magnetic field lines. Thus, in one-and two-dimensional simulations cross-field diffusion is inhibited and important physics may be lost. We have investigated cross-field diffusion in self-consistent 3-D magnetic turbulence by fully 3-dimensional hybrid simulation (macro-particle ions, massless electron fluid. The turbulence is generated by the electromagnetic ion/ion beam instability. A cold, low density, ion beam with a high velocity stream relative to the background plasma excites the right-hand resonant instability. Such ion beams may be important in the region of the Earth's foreshock. The field turbulence scatters the beam ions parallel as well as perpendicular to the magnetic field. We have determined the parallel and perpendicular diffusion coefficient for the beam ions in the turbulent wave field. The result compares favourably well (within a factor 2 with hard-sphere scattering theory for the cross-field diffusion coefficient. The cross-field diffusion coefficient is larger than that obtained in a static field with a Kolmogorov type spectrum and similar total fluctuation power. This is attributed to the resonant behaviour of the particles in the fluctuating field.

Study of asymmetrical electric discharges using particle simulation

International Nuclear Information System (INIS)

Alves, M.V.

1990-11-01

Asymmetrical electric discharges are been widely used in the microelectronic industry. The asymmetry in the electrode areas determines the magnitude of the plasma-to-electrode voltage, V sub(a), at the powered electrode, which determines the ion bombarding energy, a critical plasma processing parameter. Two many-particle simulation codes, called P D C 1 and P D S 1, were developed. These codes are electrostatic, one-dimensional (radial) and model (a bounded plasma between two infinite cylinders or two concentric spheres that can be connected to a RLC external circuit. Both codes consider asymmetrical electrodes areas. In order to simulate electrical discharges, Monte-Carlo simulation of electron-neutral and ion-neutral collisions were included. These codes were used to study the relationship between the voltage area ratio across the sheaths, V sub(a) / V sub(b), and the electrode area ratio A sub(b) / A sub(a). Simulation results agree with experimental results and also with the analytical model that includes local ionization near the electrodes, observed to occur in almost all our simulations. (author)

Computer simulation of the martensite transformation in a model two-dimensional body

International Nuclear Information System (INIS)

Chen, S.; Khachaturyan, A.G.; Morris, J.W. Jr.

1979-05-01

An analytical model of a martensitic transformation in an idealized body is constructed and used to carry out a computer simulation of the transformation in a pseudo-two-dimensional crystal. The reaction is assumed to proceed through the sequential transformation of elementary volumes (elementary martensitic particles, EMP) via the Bain strain. The elastic interaction between these volumes is computed and the transformation path chosen so as to minimize the total free energy. The model transformation shows interesting qualitative correspondencies with the known features of martensitic transformations in typical solids

Computer simulation of the martensite transformation in a model two-dimensional body

International Nuclear Information System (INIS)

Chen, S.; Khachaturyan, A.G.; Morris, J.W. Jr.

1979-06-01

An analytical model of a martensitic transformation in an idealized body is constructed and used to carry out a computer simulation of the transformation in a pseudo-two-dimensional crystal. The reaction is assumed to proceed through the sequential transformation of elementary volumes (elementary martensitic particles, EMP) via the Bain strain. The elastic interaction between these volumes is computed and the transformation path chosen so as to minimize the total free energy. The model transformation shows interesting qualitative correspondencies with the known features of martensitic transformations in typical solids

Simulations of reactive transport and precipitation with smoothed particle hydrodynamics

Science.gov (United States)

Tartakovsky, Alexandre M.; Meakin, Paul; Scheibe, Timothy D.; Eichler West, Rogene M.

2007-03-01

A numerical model based on smoothed particle hydrodynamics (SPH) was developed for reactive transport and mineral precipitation in fractured and porous materials. Because of its Lagrangian particle nature, SPH has several advantages for modeling Navier-Stokes flow and reactive transport including: (1) in a Lagrangian framework there is no non-linear term in the momentum conservation equation, so that accurate solutions can be obtained for momentum dominated flows and; (2) complicated physical and chemical processes such as surface growth due to precipitation/dissolution and chemical reactions are easy to implement. In addition, SPH simulations explicitly conserve mass and linear momentum. The SPH solution of the diffusion equation with fixed and moving reactive solid-fluid boundaries was compared with analytical solutions, Lattice Boltzmann [Q. Kang, D. Zhang, P. Lichtner, I. Tsimpanogiannis, Lattice Boltzmann model for crystal growth from supersaturated solution, Geophysical Research Letters, 31 (2004) L21604] simulations and diffusion limited aggregation (DLA) [P. Meakin, Fractals, scaling and far from equilibrium. Cambridge University Press, Cambridge, UK, 1998] model simulations. To illustrate the capabilities of the model, coupled three-dimensional flow, reactive transport and precipitation in a fracture aperture with a complex geometry were simulated.

Spacer geometry and particle deposition in spiral wound membrane feed channels

KAUST Repository

Radu, A.I.

2014-11-01

Deposition of microspheres mimicking bacterial cells was studied experimentally and with a numerical model in feed spacer membrane channels, as used in spiral wound nanofiltration (NF) and reverse osmosis (RO) membrane systems. In-situ microscopic observations in membrane fouling simulators revealed formation of specific particle deposition patterns for different diamond and ladder feed spacer orientations. A three-dimensional numerical model combining fluid flow with a Lagrangian approach for particle trajectory calculations could describe very well the in-situ observations on particle deposition in flow cells. Feed spacer geometry, positioning and cross-flow velocity sensitively influenced the particle transport and deposition patterns. The deposition patterns were not influenced by permeate production. This combined experimental-modeling approach could be used for feed spacer geometry optimization studies for reduced (bio)fouling. © 2014 Elsevier Ltd.

Particle simulation models and their application to controlled fusion

International Nuclear Information System (INIS)

Okuda, H.

1977-01-01

Plasma simulation models using particles are described which have been developed for studying the microscopic behavior of a confined plasma in a magnetic field. The first model is developed to investigate the anomalous diffusion of particles and energy due to low-frequency electrostatic microinstabilities in cylindrical and toroidal systems. The model makes use of the combination of eigenfunction expansion in one direction and the multipole expansion on a two-dimensional spatial grid for solving the Maxwell equations and for pushing particles. The second model is developed to study the neutral beam injection heating of a tokamak plasma taking into account the spatial variation of plasma parameters and the finite ion beam banana orbit. The self-consistent electric and magnetic fields are totally ignored in this model and the Fokker-Planck collisions on the beam ions due to background ions and electrons are built in through the Monte Carlo method

Comparison of GPU-Based Numerous Particles Simulation and Experiment

International Nuclear Information System (INIS)

Park, Sang Wook; Jun, Chul Woong; Sohn, Jeong Hyun; Lee, Jae Wook

2014-01-01

The dynamic behavior of numerous grains interacting with each other can be easily observed. In this study, this dynamic behavior was analyzed based on the contact between numerous grains. The discrete element method was used for analyzing the dynamic behavior of each particle and the neighboring-cell algorithm was employed for detecting their contact. The Hertzian and tangential sliding friction contact models were used for calculating the contact force acting between the particles. A GPU-based parallel program was developed for conducting the computer simulation and calculating the numerous contacts. The dam break experiment was performed to verify the simulation results. The reliability of the program was verified by comparing the results of the simulation with those of the experiment

3D Lagrangian Model of Particle Saltation in an Open Channel Flow with Emphasis on Particle-Particle Collisions

Science.gov (United States)

Moreno, P. A.; Bombardelli, F. A.

2012-12-01

Particles laying motionless at the bed of rivers, lakes and estuaries can be put into motion when the shear stress exerted by the flow on the particles exceeds the critical shear stress. When these particles start their motion they can either remain suspended by long periods of time (suspended load) or move close to the bed (bed load). Particles are transported as bed load in three different modes: Sliding, rolling and saltation. Saltation is usually described as the bouncing motion of sediment particles in a layer a few particle diameters thick. The amount of particles and the bed-load mode in which they move depend on the particle size and density, and the flow intensity, usually quantified by the shear velocity. The bottom shear stress in natural streams will most likely be large enough to set saltation as the most important bed-load transport mechanism among all three modes. Thus, studying the saltation process is crucial for the overall understanding of bed-load transport. Particularly, numerical simulations of this process have been providing important insight regarding the relative importance of the physical mechanisms involved in it. Several processes occur when particles are saltating near the bed: i) Particles collide with the bed, ii) they "fly" between collisions with the bed, as a result of their interaction with the fluid flow, iii) and they collide among themselves. These processes can be simulated using a three-dimensional Eulerian-Lagrangian model. In order to mimic these processes we have experimented with an averaged turbulent flow field represented by the logarithmic law of the wall, and with a more involved approach in which a computed turbulent velocity field for a flat plate was used as a surrogate of the three-dimensional turbulent conditions present close to stream beds. Since flat-plate and open-channel boundary layers are essentially different, a dynamic similarity analysis was performed showing that the highly-resolved three-dimensional

FAST: a three-dimensional time-dependent FEL simulation code

International Nuclear Information System (INIS)

Saldin, E.L.; Schneidmiller, E.A.; Yurkov, M.V.

1999-01-01

In this report we briefly describe the three-dimensional, time-dependent FEL simulation code FAST. The equations of motion of the particles and Maxwell's equations are solved simultaneously taking into account the slippage effect. Radiation fields are calculated using an integral solution of Maxwell's equations. A special technique has been developed for fast calculations of the radiation field, drastically reducing the required CPU time. As a result, the developed code allows one to use a personal computer for time-dependent simulations. The code allows one to simulate the radiation from the electron bunch of any transverse and longitudinal bunch shape; to simulate simultaneously an external seed with superimposed noise in the electron beam; to take into account energy spread in the electron beam and the space charge fields; and to simulate a high-gain, high-efficiency FEL amplifier with a tapered undulator. It is important to note that there are no significant memory limitations in the developed code and an electron bunch of any length can be simulated

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)

Three-dimensional ordering of cold ion beams in a storage ring: A molecular-dynamics simulation study

Energy Technology Data Exchange (ETDEWEB)

Yuri, Yosuke, E-mail: yuri.yosuke@jaea.go.jp [Takasaki Advanced Radiation Research Institute, Japan Atomic Energy Agency, 1233 Watanuki-machi Takasaki, Gunma 370-1292 Japan (Japan)

2015-06-29

Three-dimensional (3D) ordering of a charged-particle beams circulating in a storage ring is systematically studied with a molecular-dynamics simulation code. An ion beam can exhibit a 3D ordered configuration at ultralow temperature as a result of powerful 3D laser cooling. Various unique characteristics of the ordered beams, different from those of crystalline beams, are revealed in detail, such as the single-particle motion in the transverse and longitudinal directions, and the dependence of the tune depression and the Coulomb coupling constant on the operating points.

A one-dimensional transport code for the simulation of D-T burning tokamak plasma

International Nuclear Information System (INIS)

Tone, Tatsuzo; Maki, Koichi; Kasai, Masao; Nishida, Hidetsugu

1980-11-01

A one-dimensional transport code for D-T burning tokamak plasma has been developed, which simulates the spatial behavior of fuel ions(D, T), alpha particles, impurities, temperatures of ions and electrons, plasma current, neutrals, heating of alpha and injected beam particles. The basic transport equations are represented by one generalized equation so that the improvement of models and the addition of new equations may be easily made. A model of burn control using a variable toroidal field ripple is employed. This report describes in detail the simulation model, numerical method and the usage of the code. Some typical examples to which the code has been applied are presented. (author)

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.)

Numerical Simulations of Particle Deposition in Metal Foam Heat Exchangers

Science.gov (United States)

Sauret, Emilie; Saha, Suvash C.; Gu, Yuantong

2013-01-01

Australia is a high-potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage. However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two-dimensional (2D) numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Simulating variable source problems via post processing of individual particle tallies

International Nuclear Information System (INIS)

Bleuel, D.L.; Donahue, R.J.; Ludewigt, B.A.; Vujic, J.

2000-01-01

Monte Carlo is an extremely powerful method of simulating complex, three dimensional environments without excessive problem simplification. However, it is often time consuming to simulate models in which the source can be highly varied. Similarly difficult are optimization studies involving sources in which many input parameters are variable, such as particle energy, angle, and spatial distribution. Such studies are often approached using brute force methods or intelligent guesswork. One field in which these problems are often encountered is accelerator-driven Boron Neutron Capture Therapy (BNCT) for the treatment of cancers. Solving the reverse problem of determining the best neutron source for optimal BNCT treatment can be accomplished by separating the time-consuming particle-tracking process of a full Monte Carlo simulation from the calculation of the source weighting factors which is typically performed at the beginning of a Monte Carlo simulation. By post-processing these weighting factors on a recorded file of individual particle tally information, the effect of changing source variables can be realized in a matter of seconds, instead of requiring hours or days for additional complete simulations. By intelligent source biasing, any number of different source distributions can be calculated quickly from a single Monte Carlo simulation. The source description can be treated as variable and the effect of changing multiple interdependent source variables on the problem's solution can be determined. Though the focus of this study is on BNCT applications, this procedure may be applicable to any problem that involves a variable source

Micromagnetic simulations of spinel ferrite particles

International Nuclear Information System (INIS)

Dantas, Christine C.; Gama, Adriana M.

2010-01-01

This paper presents the results of simulations of the magnetization field ac response (at 2-12 GHz) of various submicron ferrite particles (cylindrical dots). The ferrites in the present simulations have the spinel structure, expressed here by M 1 - n Zn n Fe 2 O 4 (where M stands for a divalent metal), and the parameters chosen were the following: (a) for n=0: M={Fe, Mn, Co, Ni, Mg, Cu }; (b) for n=0.1: M = {Fe, Mg} (mixed ferrites). These runs represent full 3D micromagnetic (one-particle) ferrite simulations. We find evidences of confined spin waves in all simulations, as well as a complex behavior nearby the main resonance peak in the case of the M = {Mg, Cu} ferrites. A comparison of the n=0 and n=0.1 cases for fixed M reveals a significant change in the spectra in M = Mg ferrites, but only a minor change in the M=Fe case. An additional larger scale simulation of a 3 by 3 particle array was performed using similar conditions of the Fe 3 O 4 (magnetite; n=0, M = Fe) one-particle simulation. We find that the main resonance peak of the Fe 3 O 4 one-particle simulation is disfigured in the corresponding 3 by 3 particle simulation, indicating the extent to which dipolar interactions are able to affect the main resonance peak in that magnetic compound.

Partially linearized algorithms in gyrokinetic particle simulation

Energy Technology Data Exchange (ETDEWEB)

Dimits, A.M.; Lee, W.W.

1990-10-01

In this paper, particle simulation algorithms with time-varying weights for the gyrokinetic Vlasov-Poisson system have been developed. The primary purpose is to use them for the removal of the selected nonlinearities in the simulation of gradient-driven microturbulence so that the relative importance of the various nonlinear effects can be assessed. It is hoped that the use of these procedures will result in a better understanding of the transport mechanisms and scaling in tokamaks. Another application of these algorithms is for the improvement of the numerical properties of the simulation plasma. For instance, implementations of such algorithms (1) enable us to suppress the intrinsic numerical noise in the simulation, and (2) also make it possible to regulate the weights of the fast-moving particles and, in turn, to eliminate the associated high frequency oscillations. Examples of their application to drift-type instabilities in slab geometry are given. We note that the work reported here represents the first successful use of the weighted algorithms in particle codes for the nonlinear simulation of plasmas.

Partially linearized algorithms in gyrokinetic particle simulation

International Nuclear Information System (INIS)

Dimits, A.M.; Lee, W.W.

1990-10-01

In this paper, particle simulation algorithms with time-varying weights for the gyrokinetic Vlasov-Poisson system have been developed. The primary purpose is to use them for the removal of the selected nonlinearities in the simulation of gradient-driven microturbulence so that the relative importance of the various nonlinear effects can be assessed. It is hoped that the use of these procedures will result in a better understanding of the transport mechanisms and scaling in tokamaks. Another application of these algorithms is for the improvement of the numerical properties of the simulation plasma. For instance, implementations of such algorithms (1) enable us to suppress the intrinsic numerical noise in the simulation, and (2) also make it possible to regulate the weights of the fast-moving particles and, in turn, to eliminate the associated high frequency oscillations. Examples of their application to drift-type instabilities in slab geometry are given. We note that the work reported here represents the first successful use of the weighted algorithms in particle codes for the nonlinear simulation of plasmas

Three-dimensional simulation of viscous-flow agglomerate sintering.

Science.gov (United States)

Kirchhof, M J; Schmid, H -J; Peukert, W

2009-08-01

The viscous-flow sintering of different agglomerate particle morphologies is studied by three-dimensional computer simulations based on the concept of fractional volume of fluid. For a fundamental understanding of particle sintering characteristics, the neck growth kinetics in agglomerate chains and in doublets consisting of differently sized primary particles is investigated. Results show that different sintering contacts in agglomerates even during the first stages are not completely independent from each other, even though differences are small. The neck growth kinetics of differently sized primary particles is determined by the smaller one up to a size difference by a factor of approximately 2, whereas for larger size differences, the kinetics becomes faster. In particular, the agglomerate sintering kinetics is investigated for particle chains of different lengths and for different particle morphologies each having ten primary particles and nine initial sintering contacts. For agglomerate chains, the kinetics approximately can be normalized by using the radius of the fully coalesced sphere. In general, different agglomerate morphologies show equal kinetics during the first sintering stages, whereas during advanced stages, compact morphologies show significantly faster sintering progress than more open morphologies. Hence, the overall kinetics cannot be described by simply using constant morphology correction factors such as fractal dimension or mean coordination number which are used in common sintering models. However, for the first stages of viscous-flow agglomerate sintering, which are the most important for many particle processes, a sintering equation is presented. Although we use agglomerates consisting of spherical primary particles, our methodology can be applied to other aggregate geometries as well.

Simulation of plasma double-layer structures

International Nuclear Information System (INIS)

Borovsky, J.E.; Joyce, G.

1982-01-01

Electrostatic plasma double layers are numerically simulated by means of a magnetized 2 1/2-dimensional particle-in-cell method. The investigation of planar double layers indicates that these one-dimensional potential structures are susceptible to periodic disruption by instabilities in the low-potential plasmas. Only a slight increase in the double-layer thickness with an increase in its obliqueness to the magnetic field is observed. Weak magnetization results in the double-layer electric-field alignment of accelerated particles and strong magnetization results in their magnetic-field alignment. The numerial simulations of spatially periodic two-dimensional double layers also exhibit cyclical instability. A morphological invariance in two-dimensional double layers with respect to the degree of magnetization implies that the potential structures scale with Debye lengths rather than with gyroradii. Electron-beam excited electrostatic electron-cyclotron waves and (ion-beam driven) solitary waves are present in the plasmas adjacent to the double layers

A collision model in plasma particle simulations

International Nuclear Information System (INIS)

Ma Yanyun; Chang Wenwei; Yin Yan; Yue Zongwu; Cao Lihua; Liu Daqing

2000-01-01

In order to offset the collisional effects reduced by using finite-size particles, β particle clouds are used in particle simulation codes (β is the ratio of charge or mass of modeling particles to real ones). The method of impulse approximation (strait line orbit approximation) is used to analyze the scattering cross section of β particle clouds plasmas. The authors can obtain the relation of the value of a and β and scattering cross section (a is the radius of β particle cloud). By using this relation the authors can determine the value of a and β so that the collisional effects of the modeling system is correspondent with the real one. The authors can also adjust the values of a and β so that the authors can enhance or reduce the collisional effects fictitiously. The results of simulation are in good agreement with the theoretical ones

Simulation on three dimensional bubble formation using MARS

International Nuclear Information System (INIS)

Kunugi, Tomoaki

1997-01-01

This paper describes a numerical simulation on three-dimensional bubble formation by means of the MARS (Multi-interfaces Advection and Reconstruction Solver) developed by the author. The comparison between two-dimensional and three-dimensional simulation on an agglomeration of two bubbles is discussed. Moreover, some simulation results regarding a phase change phenomena such as a boiling and condensation in a two dimensional enclosure with heated and cooled walls are presented. (author)

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.

Numerical simulation of particle settling and cohesion in liquid

Energy Technology Data Exchange (ETDEWEB)

Johno, Y; Nakashima, K; Shigematsu, T; Ono, B [SASEBO National College of Technology, 1-1 Okishin, Sasebo, Nagasaki, 857-1193 (Japan); Satomi, M, E-mail: yjohno@post.cc.sasebo.ac.j [Sony Semiconductor Kyushu Corporation, Kikuchigun, Kumamoto (Japan)

2009-02-01

In this study, the motions of particles and particle clusters in liquid were numerically simulated. The particles of two sizes (Dp=40mum and 20mum) settle while repeating cohesion and dispersion, and finally the sediment of particles are formed at the bottom of a hexahedron container which is filled up with pure water. The flow field was solved with the Navier-Stokes equations and the particle motions were solved with the Lagrangian-type motion equations, where the interaction between fluid and particles due to drag forces were taken into account. The collision among particles was calculated using Distinct Element Method (DEM), and the effects of cohesive forces by van der Waals force acting on particle contact points were taken into account. Numerical simulations were performed under conditions in still flow and in shear flow. It was found that the simulation results enable us to know the state of the particle settling and the particle condensation.

Nonlinear Delta-f Particle Simulations of Collective Effects in High-Intensity Bunched Beams

CERN Document Server

Qin, Hong; Hudson, Stuart R; Startsev, Edward

2005-01-01

The collective effects in high-intensity 3D bunched beams are described self-consistently by the nonlinear Vlasov-Maxwell equations.* The nonlinear delta-f method,** a particle simulation method for solving the nonlinear Vlasov-Maxwell equations, is being used to study the collective effects in high-intensity 3D bunched beams. The delta-f method, as a nonlinear perturbative scheme, splits the distribution function into equilibrium and perturbed parts. The perturbed distribution function is represented as a weighted summation over discrete particles, where the particle orbits are advanced by equations of motion in the focusing field and self-consistent fields, and the particle weights are advanced by the coupling between the perturbed fields and the zero-order distribution function. The nonlinear delta-f method exhibits minimal noise and accuracy problems in comparison with standard particle-in-cell simulations. A self-consistent 3D kinetic equilibrium is first established for high intensity bunched beams. The...

Particle-in-cell simulations of asymmetric guide-field reconnection: quadrupolar structure of Hall magnetic field

Science.gov (United States)

Schmitz, R. G.; Alves, M. V.; Barbosa, M. V. G.

2017-12-01

One of the most important processes that occurs in Earth's magnetosphere is known as magnetic reconnection (MR). This process can be symmetric or asymmetric, depending basically on the plasma density and magnetic field in both sides of the current sheet. A good example of symmetric reconnection in terrestrial magnetosphere occurs in the magnetotail, where these quantities are similar on the north and south lobes. In the dayside magnetopause MR is asymmetric, since the plasma regimes and magnetic fields of magnetosheath and magnetosphere are quite different. Symmetric reconnection has some unique signatures. For example, the formation of a quadrupolar structure of Hall magnetic field and a bipolar Hall electric field that points to the center of the current sheet. The different particle motions in the presence of asymmetries change these signatures, causing the quadrupolar pattern to be distorted and forming a bipolar structure. Also, the bipolar Hall electric field is modified and gives rise to a single peak pointing toward the magnetosheat, considering an example of magnetopause reconnection. The presence of a guide-field can also distort the quadrupolar pattern, by giving a shear angle across the current sheet and altering the symmetric patterns, according to previous simulations and observations. Recently, a quadrupolar structure was observed in an asymmetric guide-field MR event using MMS (Magnetospheric Multiscale) mission data [Peng et al., JGR, 2017]. This event shows clearly that the density asymmetry and the guide-field were not sufficient to form signatures of asymmetric reconnection. Using the particle-in-cell code iPIC3D [Markidis et al, Mathematics and Computers in Simulation, 2010] with the MMS data from this event used to define input parameters, we found a quadrupolar structure of Hall magnetic field and a bipolar pattern of Hall electric field in ion scales, showing that our results are in an excellent agreement with the MMS observations. To our

Response to a small external force and fluctuations of a passive particle in a one-dimensional diffusive environment

Science.gov (United States)

Huveneers, François

2018-04-01

We investigate the long-time behavior of a passive particle evolving in a one-dimensional diffusive random environment, with diffusion constant D . We consider two cases: (a) The particle is pulled forward by a small external constant force and (b) there is no systematic bias. Theoretical arguments and numerical simulations provide evidence that the particle is eventually trapped by the environment. This is diagnosed in two ways: The asymptotic speed of the particle scales quadratically with the external force as it goes to zero, and the fluctuations scale diffusively in the unbiased environment, up to possible logarithmic corrections in both cases. Moreover, in the large D limit (homogenized regime), we find an important transient region giving rise to other, finite-size scalings, and we describe the crossover to the true asymptotic behavior.

IMPETUS - Interactive MultiPhysics Environment for Unified Simulations.

Science.gov (United States)

Ha, Vi Q; Lykotrafitis, George

2016-12-08

We introduce IMPETUS - Interactive MultiPhysics Environment for Unified Simulations, an object oriented, easy-to-use, high performance, C++ program for three-dimensional simulations of complex physical systems that can benefit a large variety of research areas, especially in cell mechanics. The program implements cross-communication between locally interacting particles and continuum models residing in the same physical space while a network facilitates long-range particle interactions. Message Passing Interface is used for inter-processor communication for all simulations. Copyright © 2016 Elsevier Ltd. All rights reserved.

High-Fidelity RF Gun Simulations with the Parallel 3D Finite Element Particle-In-Cell Code Pic3P

Energy Technology Data Exchange (ETDEWEB)

Candel, A; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Schussman, G.; Ko, K.; /SLAC

2009-06-19

SLAC's Advanced Computations Department (ACD) has developed the first parallel Finite Element 3D Particle-In-Cell (PIC) code, Pic3P, for simulations of RF guns and other space-charge dominated beam-cavity interactions. Pic3P solves the complete set of Maxwell-Lorentz equations and thus includes space charge, retardation and wakefield effects from first principles. Pic3P uses higher-order Finite Elementmethods on unstructured conformal meshes. A novel scheme for causal adaptive refinement and dynamic load balancing enable unprecedented simulation accuracy, aiding the design and operation of the next generation of accelerator facilities. Application to the Linac Coherent Light Source (LCLS) RF gun is presented.

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.

Helical structures in vertically aligned dust particle chains in a complex plasma

Science.gov (United States)

Hyde, Truell W.; Kong, Jie; Matthews, Lorin S.

2013-05-01

Self-assembly of structures from vertically aligned, charged dust particle bundles within a glass box placed on the lower, powered electrode of a Gaseous Electronics Conference rf reference cell were produced and examined experimentally. Self-organized formation of one-dimensional vertical chains, two-dimensional zigzag structures, and three-dimensional helical structures of triangular, quadrangular, pentagonal, hexagonal, and heptagonal symmetries are shown to occur. System evolution is shown to progress from a one-dimensional chain structure, through a zigzag transition to a two-dimensional, spindlelike structure, and then to various three-dimensional, helical structures exhibiting multiple symmetries. Stable configurations are found to be dependent upon the system confinement, γ2=ω0h/ω0v2 (where ω0h,v are the horizontal and vertical dust resonance frequencies), the total number of particles within a bundle, and the rf power. For clusters having fixed numbers of particles, the rf power at which structural phase transitions occur is repeatable and exhibits no observable hysteresis. The critical conditions for these structural phase transitions as well as the basic symmetry exhibited by the one-, two-, and three-dimensional structures that subsequently develop are in good agreement with the theoretically predicted configurations of minimum energy determined employing molecular dynamics simulations for charged dust particles confined in a prolate, spheroidal potential as presented theoretically by Kamimura and Ishihara [Kamimura and Ishihara, Phys. Rev. EPLEEE81063-651X10.1103/PhysRevE.85.016406 85, 016406 (2012)].

Particle-based simulation of charge transport in discrete-charge nano-scale systems: the electrostatic problem.

Science.gov (United States)

Berti, Claudio; Gillespie, Dirk; Eisenberg, Robert S; Fiegna, Claudio

2012-02-16

The fast and accurate computation of the electric forces that drive the motion of charged particles at the nanometer scale represents a computational challenge. For this kind of system, where the discrete nature of the charges cannot be neglected, boundary element methods (BEM) represent a better approach than finite differences/finite elements methods. In this article, we compare two different BEM approaches to a canonical electrostatic problem in a three-dimensional space with inhomogeneous dielectrics, emphasizing their suitability for particle-based simulations: the iterative method proposed by Hoyles et al. and the Induced Charge Computation introduced by Boda et al.

Monte Carlo Simulation for Particle Detectors

CERN Document Server

Pia, Maria Grazia

2012-01-01

Monte Carlo simulation is an essential component of experimental particle physics in all the phases of its life-cycle: the investigation of the physics reach of detector concepts, the design of facilities and detectors, the development and optimization of data reconstruction software, the data analysis for the production of physics results. This note briefly outlines some research topics related to Monte Carlo simulation, that are relevant to future experimental perspectives in particle physics. The focus is on physics aspects: conceptual progress beyond current particle transport schemes, the incorporation of materials science knowledge relevant to novel detection technologies, functionality to model radiation damage, the capability for multi-scale simulation, quantitative validation and uncertainty quantification to determine the predictive power of simulation. The R&D on simulation for future detectors would profit from cooperation within various components of the particle physics community, and synerg...

Scaling Properties of Particle Density Fields Formed in Simulated Turbulent Flows

Science.gov (United States)

Hogan, Robert C.; Cuzzi, Jeffrey N.; Dobrovolskis, Anthony R.; DeVincenzi, Donald (Technical Monitor)

1998-01-01

Direct numerical simulations (DNS) of particle concentrations in fully developed 3D turbulence were carried out in order to study the nonuniform structure of the particle density field. Three steady-state turbulent fluid fields with Taylor microscale Reynolds numbers (Re(sub lambda)) of 40, 80 and 140 were generated by solving the Navier-Stokes equations with pseudospectral methods. Large scale forcing was used to drive the turbulence and maintain temporal stationarity. The response of the particles to the fluid was parameterized by the particle Stokes number St, defined as the ratio of the particle's stopping time to the mean period of eddies on the Kolmogorov scale (eta). In this paper, we consider only passive particles optimally coupled to these eddies (St approx. = 1) because of their tendency to concentrate more than particles with lesser or greater St values. The trajectories of up to 70 million particles were tracked in the equilibrated turbulent flows until the particle concentration field reached a statistically stationary state. The nonuniform structure of the concentration fields was characterized by the multifractal singularity spectrum, f(alpha), derived from measures obtained after binning particles into cells ranging from 2(eta) to 15(eta) in size. We observed strong systematic variations of f(alpha) across this scale range in all three simulations and conclude that the particle concentration field is not statistically self similar across the scale range explored. However, spectra obtained at the 2(eta), 4(eta), and 8(eta) scales of each flow case were found to be qualitatively similar. This result suggests that the local structure of the particle concentration field may be flow-Independent. The singularity spectra found for 2n-sized cells were used to predict concentration distributions in good agreement with those obtained directly from the particle data. This Singularity spectrum has a shape similar to the analogous spectrum derived for the

Particle transport simulation for spaceborne, NaI gamma-ray spectrometers

International Nuclear Information System (INIS)

Dyer, C.S.; Truscott, P.R.; Sims, A.J.; Comber, C.; Hammond, N.D.A.

1988-11-01

Radioactivity induced in detectors by protons and secondary neutrons limits the sensitivity of spaceborne gamma-ray spectrometers. Three dimensional Monte Carlo transport codes have been employed to simulate particle transport of cosmic rays and inner-belt protons in various representations of the Gamma Ray Observatory Spacecraft and the Oriented Scintillation Spectrometer Experiment. Results are used to accurately quantify the contributions to the radioactive background, assess shielding options and examine the effect of detector and space-craft orientation in anisotropic trapped proton fluxes. (author)

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

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

Particle-in-cell simulations of multi-MeV pulsed X-ray induced air plasmas at low pressures

International Nuclear Information System (INIS)

Ribière, M.; D'Almeida, T.; Gaufridy de Dortan, F. de; Maulois, M.; Delbos, C.; Garrigues, A.; Cessenat, O.; Azaïs, B.

2016-01-01

A full kinetic modelling of the charge particles dynamics generated upon the irradiation of an air-filled cavity by a multi-MeV pulsed x-ray is performed. From the calculated radiative source generated by the ASTERIX generator, we calculated the electromagnetic fields generated by x-ray induced air plasmas in a metallic cavity at different pressures. Simulations are carried out based on a Particle-In-Cell interpolation method which uses 3D Maxwell-Vlasov calculations of the constitutive charged species densities of air plasmas at different pressures at equilibrium. The resulting electromagnetic fields within the cavity are calculated for different electron densities up to 4 × 10"1"0" cm"−"3. For each air pressure, we show electronic plasma waves formation followed by Landau damping. As electron density increases, the calculations exhibit space-charged neutralization and return current formation.

Numerical Simulation of Particle Flow Motion in a Two-Dimensional Modular Pebble-Bed Reactor with Discrete Element Method

Directory of Open Access Journals (Sweden)

Guodong Liu

2013-01-01

Full Text Available Modular pebble-bed nuclear reactor (MPBNR technology is promising due to its attractive features such as high fuel performance and inherent safety. Particle motion of fuel and graphite pebbles is highly associated with the performance of pebbled-bed modular nuclear reactor. To understand the mechanism of pebble’s motion in the reactor, we numerically studied the influence of number ratio of fuel and graphite pebbles, funnel angle of the reactor, height of guide ring on the distribution of pebble position, and velocity by means of discrete element method (DEM in a two-dimensional MPBNR. Velocity distributions at different areas of the reactor as well as mixing characteristics of fuel and graphite pebbles were investigated. Both fuel and graphite pebbles moved downward, and a uniform motion was formed in the column zone, while pebbles motion in the cone zone was accelerated due to the decrease of the cross sectional flow area. The number ratio of fuel and graphite pebbles and the height of guide ring had a minor influence on the velocity distribution of pebbles, while the variation of funnel angle had an obvious impact on the velocity distribution. Simulated results agreed well with the work in the literature.

Radio Evolution of Supernova Remnants Including Nonlinear Particle Acceleration: Insights from Hydrodynamic Simulations

Science.gov (United States)

Pavlović, Marko Z.; Urošević, Dejan; Arbutina, Bojan; Orlando, Salvatore; Maxted, Nigel; Filipović, Miroslav D.

2018-01-01

We present a model for the radio evolution of supernova remnants (SNRs) obtained by using three-dimensional hydrodynamic simulations coupled with nonlinear kinetic theory of cosmic-ray (CR) acceleration in SNRs. We model the radio evolution of SNRs on a global level by performing simulations for a wide range of the relevant physical parameters, such as the ambient density, supernova (SN) explosion energy, acceleration efficiency, and magnetic field amplification (MFA) efficiency. We attribute the observed spread of radio surface brightnesses for corresponding SNR diameters to the spread of these parameters. In addition to our simulations of Type Ia SNRs, we also considered SNR radio evolution in denser, nonuniform circumstellar environments modified by the progenitor star wind. These simulations start with the mass of the ejecta substantially higher than in the case of a Type Ia SN and presumably lower shock speed. The magnetic field is understandably seen as very important for the radio evolution of SNRs. In terms of MFA, we include both resonant and nonresonant modes in our large-scale simulations by implementing models obtained from first-principles, particle-in-cell simulations and nonlinear magnetohydrodynamical simulations. We test the quality and reliability of our models on a sample consisting of Galactic and extragalactic SNRs. Our simulations give Σ ‑ D slopes between ‑4 and ‑6 for the full Sedov regime. Recent empirical slopes obtained for the Galactic samples are around ‑5, while those for the extragalactic samples are around ‑4.

Three-dimensional visualization of forming Hepatitis C virus-like particles by electron-tomography

Energy Technology Data Exchange (ETDEWEB)

Badia-Martinez, Daniel; Peralta, Bibiana [Structural Biology Unit, CIC bioGUNE, CIBERehd, 48160 Derio (Spain); Andres, German; Guerra, Milagros [Electron Microscopy Unit, Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Campus Cantoblanco, 28049 Madrid (Spain); Gil-Carton, David [Structural Biology Unit, CIC bioGUNE, CIBERehd, 48160 Derio (Spain); Abrescia, Nicola G.A., E-mail: nabrescia@cicbiogune.es [Structural Biology Unit, CIC bioGUNE, CIBERehd, 48160 Derio (Spain); IKERBASQUE, Basque Foundation for Science, 48011 Bilbao (Spain)

2012-09-01

Hepatitis C virus infects almost 170 million people per year but its assembly pathway, architecture and the structures of its envelope proteins are poorly understood. Using electron tomography of plastic-embedded sections of insect cells, we have visualized the morphogenesis of recombinant Hepatitis C virus-like particles. Our data provide a three-dimensional sketch of viral assembly at the endoplasmic reticulum showing different budding stages and contiguity of buds. This latter phenomenon could play an important role during the assembly of wt-HCV and explain the size-heterogeneity of its particles.

Three-dimensional visualization of forming Hepatitis C virus-like particles by electron-tomography

International Nuclear Information System (INIS)

Badia-Martinez, Daniel; Peralta, Bibiana; Andrés, German; Guerra, Milagros; Gil-Carton, David; Abrescia, Nicola G.A.

2012-01-01

Hepatitis C virus infects almost 170 million people per year but its assembly pathway, architecture and the structures of its envelope proteins are poorly understood. Using electron tomography of plastic-embedded sections of insect cells, we have visualized the morphogenesis of recombinant Hepatitis C virus-like particles. Our data provide a three-dimensional sketch of viral assembly at the endoplasmic reticulum showing different budding stages and contiguity of buds. This latter phenomenon could play an important role during the assembly of wt-HCV and explain the size-heterogeneity of its particles.

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.

Data parallel sorting for particle simulation

Science.gov (United States)

Dagum, Leonardo

1992-01-01

Sorting on a parallel architecture is a communications intensive event which can incur a high penalty in applications where it is required. In the case of particle simulation, only integer sorting is necessary, and sequential implementations easily attain the minimum performance bound of O (N) for N particles. Parallel implementations, however, have to cope with the parallel sorting problem which, in addition to incurring a heavy communications cost, can make the minimun performance bound difficult to attain. This paper demonstrates how the sorting problem in a particle simulation can be reduced to a merging problem, and describes an efficient data parallel algorithm to solve this merging problem in a particle simulation. The new algorithm is shown to be optimal under conditions usual for particle simulation, and its fieldwise implementation on the Connection Machine is analyzed in detail. The new algorithm is about four times faster than a fieldwise implementation of radix sort on the Connection Machine.

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.

Visual interrogation of gyrokinetic particle simulations

International Nuclear Information System (INIS)

Jones, Chad; Ma, K-L; Sanderson, Allen; Myers, Lee Roy Jr

2007-01-01

Gyrokinetic particle simulations are critical to the study of anomalous energy transport associated with plasma microturbulence in magnetic confinement fusion experiments. The simulations are conducted on massively parallel computers and produce large quantities of particles, variables, and time steps, thus presenting a formidable challenge to data analysis tasks. We present two new visualization techniques for scientists to improve their understanding of the time-varying, multivariate particle data. One technique allows scientists to examine correlations in multivariate particle data with tightly coupled views of the data in both physical space and variable space, and to visually identify and track features of interest. The second technique, built into SCIRun, allows scientists to perform range-based queries over a series of time slices and visualize the resulting particles using glyphs. The ability to navigate the multiple dimensions of the particle data, as well as query individual or a collection of particles, enables scientists to not only validate their simulations but also discover new phenomena in their data

Electron behavior in ion beam neutralization in electric propulsion: full particle-in-cell simulation

International Nuclear Information System (INIS)

Usui, Hideyuki; Hashimoto, Akihiko; Miyake, Yohei

2013-01-01

By performing full Particle-In-Cell simulations, we examined the transient response of electrons released for the charge neutralization of a local ion beam emitted from an ion engine which is one of the electric propulsion systems. In the vicinity of the engine, the mixing process of electrons in the ion beam region is not so obvious because of large difference of dynamics between electrons and ions. A heavy ion beam emitted from a spacecraft propagates away from the engine and forms a positive potential region with respect to the background. Meanwhile electrons emitted for a neutralizer located near the ion engine are electrically attracted or accelerated to the core of the ion beam. Some electrons with the energy lower than the ion beam potential are trapped in the beam region and move along with the ion beam propagation with a multi-streaming structure in the beam potential region. Since the locations of the neutralizer and the ion beam exit are different, the above-mentioned bouncing motion of electrons is also observed in the direction of the beam diameter