• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.292-300
• /
• 2024

Kinetic simulations have been performed on an Inertial Electrostatic Confinement Fusion (IECF) device. These simulations were performed using the particle-in-cell (PIC) method to analyze the behavior of ions in an IEC device and the effects of some parameters on the Confinement Time (CT). CT is an essential factor that significantly contributes to the IEC's performance as a nuclear fusion device. Using the PIC method, the geometry of a two-grided device with variable grid radius, the number of cathode grid rings, variable pressure and different dielectric thickness for the feed stalk was simulated. In this research, with the development of previous works, the interaction of particles was simulated and compared with previous results. The simulation results are in good agreement with the previous results. In these simulations, it was found that with the increase of the dielectric thickness of the feed stalk, the electric field was weakened and as a result, the confinement time was reduced. On the other hand, with the increase of the cathode radius, the confinement time increased. Using the results, an IEC device can be designed with higher efficiency and more optimal CT for ions.

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

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.6
• /
• pp.1034-1043
• /
• 2015

This is the continuation of a numerical study on vortex shedding from a blunt trailing-edge of a hydrofoil. In our previous work (Lee et al., 2015), numerical schemes for efficient computations were successfully implemented; i.e. multiple domains, the approximation of domain boundary conditions using cubic spline functions, and particle-based domain decomposition for better load balancing. In this study, numerical results through a hybrid particle-mesh method which adopts the Vortex-In-Cell (VIC) method and the Brinkman penalization model are further rigorously validated through comparison to experimental data at the Reynolds number of $2{\times}10^6$. The effects of changes in numerical parameters are also explored herein. We find that the present numerical method enables us to reasonably simulate vortex shedding phenomenon, as well as turbulent wakes of a hydrofoil.

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 2006.09b
• /
• pp.916-917
• /
• 2006

Solid Oxide Fuel Cell technology uses powder processes to produce electrodes with residual porosity by partially sintering a mixture of electronically and ionically conducting particles. We model porous fuel cell electrodes with 3D packings of monosized spherical particles. These packings are created by numerical sintering. Each particle-particle contact is characteristic for an ionic, electronic or electrochemical resistance. The numerical packing is then discretized into a resistor network which is solved by using Kirchhoff's current law to evaluate the electrode's electrochemical performance. We investigate in particular percolation effects in functionally graded electrodes as compared to other types of electrodes.

• Journal of Astronomy and Space Sciences
• /
• v.37 no.2
• /
• pp.85-94
• /
• 2020

Equatorial noise, also known magnetosonic waves (MSWs), are one of the frequently observed plasma waves in Earth's inner magnetosphere. Observations have shown that wave amplitudes maximize at the magnetic equator with a narrow extent in their latitudinal distribution. It has been understood that waves are generated from an equatorial source region and confined within a few degrees magnetic latitude. The present study investigates whether the MSW instability and saturation amplitudes maximize at the equator, given an energetic proton ring-like distribution derived from an observed wave event, and using linear instability analysis and particle-in-cell simulations with the plasma conditions at different latitudes along the dipole magnetic field line. The results show that waves initially grow fastest (i.e., with the largest growth rate) at high latitude (20°-25°), but consistent with observations, their saturation amplitudes maximize within ±10° latitude. On the other hand, the slope of the saturation amplitudes versus latitude revealed in the present study is not as steep as what the previous statistical observation results suggest. This may be indicative of some other factors not considered in the present analyses at play, such as background magnetic field and plasma inhomogeneities and the propagation effect.

• Journal of Astronomy and Space Sciences
• /
• v.37 no.1
• /
• pp.35-42
• /
• 2020

Particle-in-cell simulations were performed to understand the interaction of the solar wind with localized magnetic fields on the sunlit surface of the Moon. The results indicated a mini-magnetosphere was formed which had a thin magnetopause with the thickness of the electron skin depth. It was also found that the solar wind penetrated into the cavity of the magnetosphere intermittently rather than in a steady manner. The solar wind that moved around the magnetosphere was observed to hit the surface of the Moon, implying that it may be the cause of the lunar swirl formation on the surface.

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

• 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 computational fluids engineering
• /
• v.19 no.2
• /
• pp.72-78
• /
• 2014

A hybrid particle-mesh method as the combination between the Vortex-In-Cell (VIC) method and penalization method has been achieved in recent years. The VIC method, which is based on the vorticity-velocity formulation, offers particle-mesh algorithms to numerically simulate flows past a solid body. The penalization method is used to enforce boundary conditions at a body surface with a decoupling between body boundaries and computational grids. The main advantage of the hybrid particle-mesh method is an efficient implementation for solid boundaries of arbitrary complexity on Cartesian grids. However, a numerical simulation of flows in large domains is still not too easy. In this study, a multi-domain approach is thus proposed to further reduce computation cost and easily implement it. We validate the implementation by numerical simulations of an incompressible viscous flow around an impulsively started circular cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.5
• /
• pp.741-749
• /
• 2001

An inhomogeneous hyperbolic electric field is established among two-dimensional quadrupole electrodes to which an ac voltage is applied. Conditions under which charged particles are focused into a narrow axis region of the plug laminar flow are discussed. The aerodynamic forces influence the behavior of the charged particles in the quadrupole electric field. We derived the dimensionless equations of motion of a charged particle in the alternating quadrupole electric field, and discussed particle trajectories and focusing performance in terms of two dimensionless parameters, which are functions of particle size, operating pressure, and the amplitude and frequency of applied AC voltage, with the results of numerical simulations and experiments.