• Journal of The Korean Astronomical Society
• /
• v.54 no.3
• /
• pp.103-112
• /
• 2021

The intracluster medium (ICM) is expected to experience on average about three passages of weak shocks with low sonic Mach numbers, M ≲ 3, during the formation of galaxy clusters. Both protons and electrons could be accelerated to become high energy cosmic rays (CRs) at such ICM shocks via diffusive shock acceleration (DSA). We examine the effects of DSA by multiple shocks on the spectrum of accelerated CRs by including in situ injection/acceleration at each shock, followed by repeated re-acceleration at successive shocks in the test-particle regime. For simplicity, the accelerated particles are assumed to undergo adiabatic decompression without energy loss and escape from the system, before they encounter subsequent shocks. We show that in general the CR spectrum is flattened by multiple shock passages, compared to a single episode of DSA, and that the acceleration efficiency increases with successive shock passages. However, the decompression due to the expansion of shocks into the cluster outskirts may reduce the amplification and flattening of the CR spectrum by multiple shock passages. The final CR spectrum behind the last shock is determined by the accumulated effects of repeated re-acceleration by all previous shocks, but it is relatively insensitive to the ordering of the shock Mach numbers. Thus multiple passages of shocks may cause the slope of the CR spectrum to deviate from the canonical DSA power-law slope of the current shock.

• Journal of the Korea Society for Simulation
• /
• v.31 no.4
• /
• pp.1-10
• /
• 2022

This paper describes the development of a modular simulation environment to analyze the performance of guided weapons. For the reusability of the simulation environment, components of the guided weapon simulation are modularized based on MATLAB Simulink. And the module management system is implemented based on Excel and MATLAB GUI for simulation interface and module management. In addition, a method of integration for modules with different interface and a method to set up the guidance and control phase for designing a guided flight scenario are suggested. Finally, to verify the performance of the implemented simulation environment, it was compared with the existing simulation results.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.2
• /
• pp.51.1-51.1
• /
• 2021

Merger shocks with M_s < ~ 3 - 4 have been detected in galaxy clusters through radio observations of synchrotron radiations emitted from cosmic-ray (CR) electrons. The CR electrons are believed to be produced by the so-called diffusive shock acceleration (DSA) at the merger shocks. To describe the acceleration of electrons, the injection into DSA has to be understood. Recent studies have showed that electrons could be energized through stochastic shock drift acceleration (SSDA), a mechanism mediated by multi-scale plasma waves at shock transition zone. However, such preacceleration process seems to be effective only at the supercritical shocks with M_s > ~ 2.3, implying that further studies should be done to explain radio relics with weaker shocks. In this talk, we present the results obtained by fully kinetic 2D particle-in-cell (PIC) simulations, which include pre-existing suprathermal electrons possibly ejected from active galactic nuclei (AGNs) or produced by previous episodes of turbulence/shocks. The simulations indicate that the pre-existing electrons enhance the upstream plasma waves in shocks with M_s < ~ 2.3. However, the wavelength of such waves is not long enough to scatter off suprathermal electrons and energize them to the injection momentum for DSA. Hence, we conclude that preexciting suprathermal electrons alone would not solve the problem of electron acceleration at radio relic shocks.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.2
• /
• pp.55.2-55.2
• /
• 2019

The firehose instability is driven by a pressure anisotropy in a magnetized plasma when the temperature along the magnetic field is higher than the perpendicular temperature. Such condition occurs commonly in astrophysical and space environments, for instance, when there are beams aligned with the background magnetic field. Recently, it was argued that, in weak quasi-perpendicular shocks in the high-β intracluster medium (ICM), shock-reflected electrons propagating upstream cause the temperature anisotropy. This electron temperature anisotropy can trigger the electron firehose instability (EFI), which excites oblique waves in the shock foot. Scattering of electrons by these waves enables multiple cycles of shock drift acceleration (SDA) in the preshock region, leading to the electron injection to diffusive shock acceleration (DSA). In the study, the kinetic properties of the EFI are examined by the linear stability analysis based on the kinetic Vlasov-Maxwell theory and then further investigated by 2D Particle-in-Cell (PIC) simulations, especially focusing on those in high-β (β~100) plasmas. We then discuss the basic properties of the firehose instability, and the implication of our work on electron acceleration in ICM shock.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.2
• /
• pp.47.2-47.2
• /
• 2019

Observations indicate that turbulence in molecular clouds of the interstellar medium (ISM) is highly supersonic (M >> 1) and strongly magnetized (β ≈ 0.1), while in the intracluster medium (ICM) it is subsonic (M <~1) and weakly magnetized (β ≈ 100). Here, M is the turbulent Mach number and β is the ratio of the gas to magnetic pressures. Although magnetohydrodynamic (MHD) turbulence in such environments has been previously studied through numerical simulations, some of its properties as well as its consequences are not yet fully described. In this talk, we report a study of MHD turbulence in molecular clouds and the ICM using a newly developed code based the high-order accurate, WENO (Weighted Essentially Non-Oscillatory) scheme. The simulation results using the WENO code are generally in agreement with those presented in the previous studies with, for instance, a TVD code (Porter et al. 2015 &, Park & Ryu 2019), but reveal more detailed structures on small scales. We here present and compare the properties of simulated turbulences with WENO and TVD codes, such as the spatial distribution of density, the density probability distribution functions, and the power spectra of kinetic and magnetic energies. We also describe the populations of MHD shocks and the energy dissipation at the shocks. Finally, we discuss the implications of this study on star formation processes in the ISM and shock dissipation in the ICM.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.2
• /
• pp.54.3-54.3
• /
• 2019

During the formation of large-scale structures in the universe, shocks with the sonic Mach number Ms <~ 5 are naturally induced by supersonic flow motions of baryonic matter in the intracluster medium (ICM). Cosmic rays (CRs) are expected to be accelerated via diffusive shock acceleration (DSA) at these ICM shocks, although the existence of CR protons in the ICM remains to be confirmed through gamma-ray observations. Based on the results obtained from kinetic plasma simulations, we build an analytic DSA model for weak, quasi-parallel shocks in the test-particle regime. With our DSA model, the CR acceleration efficiency ranges ~ 0.001 - 0.02 in supercritical quasi-parallel shocks with sonic Mach number Ms ~ 2.25 - 5, and the acceleration would be negligible in subcritical shocks wth Ms <~ 2.25. Adopting our DSA model, we estimate gamma-ray and neutrino emissions from clusters of galaxies by performing cosmological hydrodynamic simulations. The estimated gamma-ray flux is below the Fermi-LAT upper limit. In addition, the possible neutrino emission due to the decay of charged pions in galaxy clusters would be about <~ 1% of the atmospheric neutrino intensity in the energy range of <~ 100 GeV. In this talk, we will discuss the implication of our results.

• The Bulletin of The Korean Astronomical Society
• /
• v.45 no.1
• /
• pp.37.1-37.1
• /
• 2020

Cosmic-ray protons (CRp) are efficiently produced at starburst galaxies (SBGs), where the star formation rate (SFR) rate is high. In this talk, we present estimates of gamma-ray and neutrino emissions from nearby SBGs, M82, NGC253, and Arp220. Inside the starburst nucleus (SBN), CRp are accelerated at supernova remnant (SNR) shocks as well as at stellar wind (SW) termination shocks, and their transport is governed by the advection due to starburst-driven wind and diffusion mediated by turbulence. We here model the momentum distributions of SNR and SW-produced CRp with single or a double power-law forms. We also employ two different diffusion models, where CRp are resonantly scattered off large-scale turbulence in SBN or self-excited waves driven by CR streaming instability. We then calculate gamma-ray/neutrino fluxes. The observed gamma-ray fluxes by Fermi-LAT, Veritas, and H.E.S.S are well reproduced with double power-law distribution for SNR-produced CRp and the CRp diffusion by self-excited turbulence. The estimated neutrino fluxes are <~10-3 of the atmospheric neutrino flux in the energy range of Eneutrino <~100 GeV and <~10-1 of the IceCube point source sensitivity in the energy range of Eneutrino >~60 TeV.

• The Bulletin of The Korean Astronomical Society
• /
• v.45 no.1
• /
• pp.52.2-52.2
• /
• 2020

At quasi-perpendicular shocks in the high-�� (��=Pgas/Pmag~100) intracluster medium (ICM), various microinstabilities occur by the temperature anisotropies and/or drift motions of plasma. In the downstream, the Alfvén ion cyclotron instability (AIC) due to the ion temperature anisotropy (Ti⊥>Ti║) is triggered by shock-reflected ions, the whistler instability (WI) is driven by the electron temperature anisotropy (Te⊥>Te║) as a consequence of the shock compression of magnetic fields, and the mirror instability is generated due to the ion and/or electron temperature anisotropy. At the shock foot, the modified two stream instability (MTSI) is possibly excited by the cross-field drift between ions and electrons. In the upstream, electron firehose instability (EFI) is driven by the electron temperature anisotropy or the relative drift between incoming and reflected electrons. These microinstabilities play important roles in the particle acceleration in ICM shocks, so understanding of the microinstabilities and the resultant plasma waves is essential. In this study, based on a linear stability analysis, the basic properties of the microinstabilities in ICM shocks and the ion/electron scale fluctuations are described. We then discuss the implication of our work on the electron pre-acceleration in ICM shocks.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.1
• /
• pp.36.2-36.2
• /
• 2021

In the outskirts of galaxy clusters, weak shocks with M_s < ~3 appear as radio relics where the synchrotron radiation is emitted from cosmic-ray (CR) electrons. To understand the production of CR electrons through the so-called diffusive shock acceleration (DSA), the electron injection into the DSA process at shocks in the hot intracluster medium (ICM) has to be described. However, the injection remains as an unsolved, outstanding problem. To explore this problem, 2D Particle-in-Cell (PIC) simulations were performed. In this talk, we present the electron preacceleration mechanism mediated by multi-scale plasma waves in the shock transition zone. In particular, we find that the electron preacceleration is effective only in the supercritical shocks, which have the sonic Mach number M_s > M_crit ≈ 2.3 in the high-beta (β~100) plasma of the ICM, because the Alfven ion cyclotron instability operates and hence multi-scale plasma waves are induced only in such supercritical shocks. Our findings will help to understand the nature of radio relics in galaxy clusters.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.25 no.4
• /
• pp.299-305
• /
• 2015

A statistical annoyance model to railway noise at platform was proposed by jury evaluation test performed in hearing laboratory. ITX-Saemaeul and Mugunghwa were chosen as the noise sources of the test, and announcement sound was included to simulate real situation. Logistic regression analysis produced %HA_LAB curve. Hosmer-Lemeshow test and classification accuracy test were used to verify the model's statistical significance. It was shown that the model which was generated from relatively small number of samples is statistically significant.