• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.568-574
• /
• 2020

The feasibility of the particle-in-cell (PIC) method is assessed to simulate the non-collective phenomena like non-collective Thomson scattering (TS). The non-collective TS in the laser-plasma interaction, which is related to the single-particle behavior, is simulated through a 2D relativistic PIC code (XOOPIC). For this simulation, a non-collective TS is emitted from a 50-50 DT plasma with electron density and temperature of n_e = 3.00 × 10¹³ cm^-3 and T_e = 1000 eV, typical for the edge plasma at ITER measured by ETS system, respectively. The wavelength, intensity, and FWHM of the laser applied in the ETS system are λ_i,0 = 1.064 × 10^-4 cm, I_i = 2.24 × 10¹⁷ erg=s·㎠, and 12.00 ns, respectively. The electron density and temperature predicted by the PIC simulation, obtained from the TS scattered wave, are n_e,TS = 2.91 × 10¹³ cm^-3 and T_e,TS = 1089 eV, respectively, which are in accordance with the input values of the simulated plasma. The obtained results indicate that the ambiguities rising due to the contradiction between the PIC statistical collective mechanism caused by the super-particle concept and the non-collective nature of TS are resolved. The ability and validity to use PIC method to study the non-collective regimes are verified.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.242.1-242.1
• /
• 2014

The electron cyclotron resonance plasma source with a belt-type magnet assembly (BMA) is designed for effective plasma confinements. For characterizing the plasma source, the plasma parameters are measured by Langmuir probe. However, the plasma parameters and the motion of charged particles near the ECR zone are not easy to diagnostics, because of the high plasma density and temperature. Thus, as an alternative method, the electromagnetic simulation of the plasma source has been performed by using three-dimensional particle-in-cell and Monte Carlo collisional (PIC-MCC) simulation codes. For considering the limitation of simulation resources and time, the periodic boundary condition is applied and the coulomb collision is neglected. In this paper, we present the results of 3D PIC simulations of ECR plasmas with BMA and we compare them with the experimental results.

• Journal of Korea Society of Waste Management
• /
• v.34 no.7
• /
• pp.685-696
• /
• 2017

This study focuses on computational particle fluid dynamics (CPFD) modeling for the fast pyrolysis of biomass in a conical spouted bed reactor. The CPFD simulation was conducted to understand the hydrodynamics, heat transfer, and biomass fast pyrolysis reaction of the conical spouted bed reactor and the multiphase-particle in cell (MP-PIC) model was used to investigate the fast pyrolysis of biomass in a conical spouted bed reactor. A two-stage semi-global kinetics model was applied to model the fast pyrolysis reaction of biomass and the commercial code (Barracuda) was used in simulations. The temperature of solid particles in a conical spouted bed reactor showed a uniform temperature distribution along the reactor height. The yield of fast pyrolysis products from the simulation was compared with the experimental data; the yield of fast pyrolysis products was 74.1wt.% tar, 17.4wt.% gas, and 8.5wt.% char. The comparison of experimental measurements and model predictions shows the model's accuracy. The CPFD simulation results had great potential to aid the future design and optimization of the fast pyrolysis process for biomass.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.557-560
• /
• 2008

A two-dimensional particle-in-cell(PIC) simulation with a Monte-Carlo Collision(MCC) has been developed to investigate the discharge characteristics of the acceleration channel of a HET. The dynamics of electrons and ions are treated with PIC method at the time scale of electrons in order to investigate the particle transport. The densities of charged particles are coupled with Poisson's equation. Xenon neutrals are injected from the anode and experience elastic, excitation, and ionization collisions with electrons, and are scattered by ions. These collisions are simulated by using an MCC model. The effects of control parameters such as magnetic field profile, electron current density, and the applied voltage have been investigated. The secondary electron emission on the dielectric surface is also considered.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.11
• /
• pp.932-938
• /
• 2017

Particle-in-cell method which blends Eulerian grids and Lagrangian particle is utilized to solve simplified hall-effect thruster. Since this study individually tracks not only neutrons and ions but also electrons, message passing interface(mpi) scheme is adopted for parallel computer cluster. Helical movement of an electron cloud in constant magnetic field is validated comparing with an exact solution. A plasma in radial magnetic field and axial electric field in a reaction cylinder is established. Electrons do double helix movement and are well anchored in a cylinder. Ionization of neutrons by impact with high-speed electrons generates ion particles. They are accelerated by axial electric field, which forms a plume of a plasma-effect thruster.

• The Bulletin of The Korean Astronomical Society
• /
• v.37 no.1
• /
• pp.67.1-67.1
• /
• 2012

Shocks are evolved when the relativistic jets in active galactic nuclei (AGNs), black hole binaries, supernova remnants (SNR) and gamma-ray bursts (GRBs) interact with the surrounding medium. The high energy particles are believed to be accelerated by the diffusive shock acceleration and the strong magnetic field is generated by Weibel instability in the shock. When ultrarelativistic electrons with strong magnetic field cool by the synchrotron emission, the radiation is observed in gamma-ray burst and the near-equipartitioned magnetic field in the external shock delays the afterglow emission. In this paper, we performed the 3D particle-in-cell (PIC) simulations to understand the characteristics of these relativistic shock and particle acceleration. Forward and reverse shocks are shaped while the unmagnetized injecting jet interacts with the unmagnetized ambient medium. Both upstream and downstream become thermalized and the particle accelerations are shown in each transition region of the shock structures.

• Journal of Astronomy and Space Sciences
• /
• v.32 no.1
• /
• pp.45-50
• /
• 2015

We have developed a 2.5-dimensional electromagnetic particle simulation code using the particle-in-cell (PIC) method to investigate electromagnetic phenomena that occur in space plasmas. Our code is based on the leap-frog method and the centered difference method for integration and differentiation of the governing equations. We adopted the relativistic Buneman-Boris method to solve the Lorentz force equation and the Esirkepov method to calculate the current density while maintaining charge conservation. Using the developed code, we performed test simulations for electron two-stream instability and electron temperature anisotropy induced instability with the same initial parameters as used in previously reported studies. The test simulation results are almost identical with those of the previous papers.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2012.02a
• /
• pp.528-528
• /
• 2012

Recently, atmospheric pressure plasmas attract lots of interests for the useful applications such as surface modification and bio-medical treatment. In this study, a particle-in-cell Monte Carlo collision (PIC-MCC) simulation was adopted to investigate the discharge characteristics of a planar micro dielectric barrier discharge (DBD) with a driving frequency from 13.56 MHz to 162.72 MHz and with a gap distance of 80 micrometers. The variation of frequency, in the change in the electron energy probability function (EEPF). Through the relation between the ion trajectories and the frequency, results in the change of EEPFs is achievable with the turning point of frequency mode. Therefore, it is possible to categorize the efficient operation range of DBDs for its applications by controlling the interactions between plasmas and neutral gas for the generation of preferable radicals.

• The Bulletin of The Korean Astronomical Society
• /
• v.36 no.2
• /
• pp.109.1-109.1
• /
• 2011

Shock waves are ubiquitous in astrophysical environments. In particular, shocks formed by merger of subclumps, infall of matter and internal flow motion in intracluster media (ICMs) and cluster outskirts are relatively weak with Mach number M ${\lesssim}$ a few. At such weak shocks, it has been believed that the diffusive shock acceleration (DSA) of cosmic rays is rather inefficient. Yet, the presence of nonthermal phenomena, such as radio halos and relics, suggests that contrary to the expectation, DSA as well as magnetic field amplification should operate at weak shocks in cluster environments. We recently initiated a study of weak, collisionless, astrophysical shocks using a PIC(Particle-in-Cell) code. The PIC code describes the motion of electron and ion particles under the electromagnetic field which is represented in grid zones. Here, we present a preliminary work of one-dimensional simulations. We show how shocks are set up as the turbulent electromagnetic field is developed in the shock transition layer, and discuss the implication on DSA and magnetic field amplification.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.61 no.12
• /
• pp.1880-1885
• /
• 2012

Nonlinear dynamical behaviors in thermionic low pressure discharge are investigated using a particle-in-cell(PIC) simulation. An electrostatic PIC code is developed to model the plasma discharge system including the kinetic effects. The elastic collision, excitation collision, ionization collision, and electron-ion recombination collision are considered in this code. The generated electrons and ions are traced to analyze physical characteristics of the plasma. The simulation results show that the nonlinear oscillation structures are observed for cold plasma in the system and the similar structures are observed for warm plasma with a shift in values of the bifurcation parameter. The detailed oscillation process can be subdivided into three distinct mode; anode-glow, temperature-limited, and double-layer modes.