• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.03a
• /
• pp.852-855
• /
• 2010

Numerical simulations for large scale triaxial tests with large diameter rockfill materials are conducted using distinct element method. For generation of compacted assembly with specific grain size distribution and initial material porosity, the clump logic method and expansion of generated particles are adapted. With micro parameters which are chosen by calibration process, discrete particle modelling of triaxial test in case of other confining stress and cyclic loading condition were conducted. Also numerical simulations of fluid injection into particulate materials were conducted to observe cavity initiation and propagation using distinct element method. The fluid scheme solves the continuity and Navior-Stokes equations numerically, then derives pressure and velocity vectors for fixed grid by considering the existence of particles within the fluid cell.

• Geomechanics and Engineering
• /
• v.9 no.6
• /
• pp.815-827
• /
• 2015

The mechanical behavior of soil and soil-rock mixture is investigated via the discrete element method. A non-overlapping combination method of spheres is used to model convex polyhedron rock blocks of soil-rock mixture in the DEM simulations. The meso-mechanical parameters of soil and soil-rock interface in DEM simulations are obtained from the in-situ tests. Based on the Voronoi cell, a method representing volumtric strain of the sample at the particle scale is proposed. The numerical results indicate that the particle rotation, occlusion, dilatation and self-organizing force chains are a remarkable phenomena of the localization band for the soil and soil-rock mixture samples. The localization band in a soil-rock mixture is wider than that in the soil sample. The current research shows that the 3D discrete element method can effectively simulate the mechanical behavior of soil and soil-rock mixture.

• Korea-Australia Rheology Journal
• /
• v.18 no.4
• /
• pp.171-181
• /
• 2006

We present a short review for authors' previous work on direct numerical simulations for inertialess hard particle suspensions formulated either with a Newtonian fluid or with viscoelastic polymeric fluids to understand the microstructural evolution and the bulk material behavior. We employ two well-defined bi-periodic domain concepts such that a single cell problem with a small number of particles may represent a large number of repeated structures: one is the sliding bi-periodic frame for simple shear flow and the other is the extensional bi-periodic frame for planar elongational flow. For implicit treatment of hydrodynamic interaction between particle and fluid, we use the finite-element/fictitious-domain method similar to the distributed Lagrangian multiplier (DLM) method together with the rigid ring description. The bi-periodic boundary conditions can be effectively incorportated as constraint equations and implemented by Lagrangian multipliers. The bulk stress can be evaluated by simple boundary integrals of stresslets on the particle boundary in such formulations. Some 2-D example results are presented to show effects of the solid fraction and the particle configuration on the shear and elongational viscosity along with the micro-structural evolution for both particles and fluid. Effects of the fluid elasticity has been also presented.

• Journal of the Optical Society of Korea
• /
• v.13 no.1
• /
• pp.8-14
• /
• 2009

In a laser-plasma wakefield accelerator, the ponderomotive force of an ultrashort high intensity laser pulse excites a longitudinal wave or plasma bubble in a way similar to the excitation of a wake wave behind a boat as it propagates on the water surface. Electric fields inside the plasma bubble can be several orders of magnitude higher than those available in conventional RF-based particle accelerator facilities which are limited by material breakdown. Therefore, if an electron bunch is properly phase-locked with the bubble's acceleration field, it can gain relativistic energies within an extremely short distance. Here, in the bubble regime we show the generation of stable and reproducible sub GeV, and GeV-class electron beams. Supported by three-dimensional particle-in-cell simulations, our experimental results show the highest acceleration gradients produced so far. Simulations suggested that the plasma bubble elongation should be minimized in order to achieve higher electron beam energies.

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

For two-dimensional grid electrodes immersed in plasmas, sheath expansion due to negative high-voltage pulse applied to the electrode generates high-energy pseudowave. The high-energy pseudowave can be used as ion beam for ion implantation. To estimate ion dose due to high-energy pseudowave, investigation on sheath expansion of grid electroes is necessary. To investigate sheath expansion, an analytic model was developed by Vlasov equation and applying the 1-D sheath expansion model to 2-D. Because of lack of generalized 2-D Child-Langmuir current, model cannot give solvable equation. Instead, for a given grid electrode geometry, the model found the relations between ion distribution functions, Child-Langmuir currents, and sheath expansions. With these relations and particle-in-cell (PIC) simulations, for given grid electrode geometry, computation time was greatly reduced for various conditions such as electrode voltages, plasma densities, and ion species. The model was examined by PIC simulations and experiments, and they well agreed.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.49-57
• /
• 2014

In the present study, 2-D gas-kinetic flow solver on unstructured meshes was developed for flows from continuum to transitional regimes. The gas-kinetic BGK scheme is based on numerical solutions of the BGK simplification of the Boltzmann transport equation. In the initial reconstruction, the unstructured version of the linear interpolation is applied to compute left and right states along a cell interface. In the gas evolution step, the numerical fluxes are computed from the evaluation of the time-dependent gas distribution function around a cell interface. Two-dimensional compressible flow calculations were performed to verify the accuracy and robustness of the current gas-kinetic approach. Gas-kinetic BGK scheme was successfully applied to two-dimensional steady and unsteady flow simulations with strong contact discontinuities. Exemplary hypersonic viscous simulations have been conducted to analyze the performances of the gas-kinetic scheme. The computed results show fair agreement with other standard particle-based approaches for both continuum part and transitional part.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.49.1-49.1
• /
• 2019

Giant radio relics in the outskirts of galaxy clusters have been observed and they are interpreted as synchrotron emission from relativistic electrons accelerated via diffusive shock acceleration (DSA) in weak shocks of Ms < 3.0. In the DSA theory, the particle momentum should be greater than a few times the momentum of thermal protons to cross the shock transition and participate in the Fermi acceleration process. In the equilibrium, the momentum of thermal electrons is much smaller than the momentum of thermal protons, so electrons need to be pre-accelerated before they can go through DSA. To investigate such electron injection process, we study the electron pre-acceleration in weak quasi-perpendicular shocks (Ms = 2.0 - 3.0) in an ICM plasma (kT = 8.6 keV, beta = 100) through 2D particle-in-cell simulations. It is known that in quasi-perpendicular shocks, a substantial fraction of electrons could be reflected upstream, gain energy via shock drift acceleration (SDA), and generate oblique waves via the electron firehose instability (EFI), leading the energization of electrons through wave-particle interactions. We find that such kinetic processes are effective only in supercritical shocks above a critical Mach number, $Ms{\ast}{\sim}2.3$. In addition, even in shocks with Ms > 2.3, energized electrons may not reach high energies to be injected to DSA, because the oblique EFI alone fails to generate long-wavelength waves. Our results should have implications for the origin and nature of radio relics.

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

• The Bulletin of The Korean Astronomical Society
• /
• v.43 no.2
• /
• pp.42.1-42.1
• /
• 2018

Collisionless shocks with low sonic Mach numbers, M<4, are expected to accelerate cosmic ray (CR) protons via diffusive shock acceleration (DSA) in the intracluster medium (ICM). However, observational evidence for CR protons in the ICM has yet to be established. Performing particle-in-cell simulations, we study the injection of protons into DSA and the early development of a nonthermal particle population in weak shocks in high ${\beta}$ plasmas. Reflection of incident protons, self-excitation of plasma waves via CR-driven instabilities, and multiple cycles of shock drift acceleration are essential to the early acceleration of CR protons in supercritical quasi-parallel shocks. We find that only in ICM shocks with $M{\geq}2.3$, a sufficient fraction of incoming protons are reflected by the overshoot in the shock electric potential and magnetic mirror at locally perpendicular magnetic fields, leading to efficient excitation of magnetic waves via CR streaming instabilities and the injection into the DSA process. Since a significant fraction of ICM shocks have M < 2.3 CR proton acceleration in the ICM might be less efficient than previously expected. This may explain why the diffuse gamma-ray emission from galaxy clusters due to proton-proton collisions has not been detected so far.

• Journal of the Optical Society of Korea
• /
• v.13 no.1
• /
• pp.42-47
• /
• 2009

The propagation of an intense laser pulse through an upward density-gradient plasma in a self-modulated laser wakefield acceleration (SM-LWFA) is investigated by using particle-in-cell (PIC) simulations. In the fully relativistic and kinetic PIC simulations, the relativistic and kinetic effects including Landau damping enhance the electron dephasing. This electron dephasing is the most important factor for limiting the energy of accelerated electrons. However, the electron dephasing, which is enhanced by relativistic and kinetic effects in the homogeneous plasma, can be forestalled through the detuning process arising from the longitudinal density gradient. Simulation results show that the detuning process can effectively maintain the coherence of the laser wake wave in the spatiotemporal wakefield pattern, hence considerable energy enhancement is achievable. The spatiotemporal profiles are analyzed for the detailed study on the relativistic and kinetic effects. In this paper, the optimum slope of the density gradient for increasing electron energy is presented for various laser intensities.