• Earthquakes and Structures
• /
• v.3 no.3_4
• /
• pp.341-363
• /
• 2012

The spectral representation method is a quick and versatile tool for the generation of spatially variable, response-spectrum-compatible simulations to be used in the nonlinear seismic response evaluation of extended structures, such as bridges. However, just as recorded data, these simulated accelerations require processing, but, unlike recorded data, the reasons for their processing are purely numerical. Hence, the criteria for the processing of acceleration simulations need to be tied to the effect of processing on the structural response. This paper presents a framework for processing acceleration simulations that is based on seismological approaches for processing recorded data, but establishes the corner frequency of the high-pass filter by minimizing the effect of processing on the response of the structural system, for the response evaluation of which the ground motions were generated. The proposed two-step criterion selects the filter corner frequency by considering both the dynamic and the pseudo-static response of the systems. First, it ensures that the linear/nonlinear dynamic structural response induced by the processed simulations captures the characteristics of the system's dynamic response caused by the unprocessed simulations, the frequency content of which is fully compatible with the target response spectrum. Second, it examines the adequacy of the selected estimate for the filter corner frequency by evaluating the pseudo-static response of the system subjected to spatially variable excitations. It is noted that the first step of this two-fold criterion suffices for the establishment of the corner frequency for the processing of acceleration time series generated at a single ground-surface location to be used in the seismic response evaluation of, e.g. a building structure. Furthermore, the concept also applies for the processing of acceleration time series generated by means of any approach that does not provide physical considerations for the selection of the corner frequency of the high-pass filter.

• Journal of The Korean Astronomical Society
• /
• v.36 no.1
• /
• pp.1-12
• /
• 2003

Nonthermal particles can be produced due to incomplete thermalization at collisionless shocks and further accelerated to very high energies via diffusive shock acceleration. In a previous study we explored the cosmic ray (CR) acceleration at cosmic shocks through numerical simulations of CR modified, quasi-parallel shocks in 1D plane-parallel geometry with the physical parameters relevant for the shocks emerging in the large scale structure formation of the universe (Kang & Jones 2002). Specifically we considered pancake shocks driven by accretion flows with $U_o = 1500 km\;s^{-l}$ and the preshock gas temperature of $T_o = 10^4 - 10^8K$. In order to consider the CR acceleration at shocks with a broader range of physical properties, in this contribution we present additional simulations with accretion flows with $U_o = 75 - 1500 km\;s^{-l}$ and $T_o = 10^4K$. We also compare the new simulation results with those reported in the previous study. For a given Mach number, shocks with higher speeds accelerate CRs faster with a greater number of particles, since the acceleration time scale is $t_{acc}\;{\propto}\;U_o^{-2}$. However, two shocks with a same Mach number but with different shock speeds evolve qualitatively similarly when the results are presented in terms of diffusion length and time scales. Therefore, the time asymptotic value for the fraction of shock kinetic energy transferred to CRs is mainly controlled by shock Mach number rather than shock speed. Although the CR acceleration efficiency depends weakly on a well-constrained injection parameter, $\epsilon$, and on shock speed for low shock Mach numbers, the dependence disappears for high shock Mach numbers. We present the 'CR energy ratio', ${\phi}(M_s)$, for a wide range of shock parameters and for $\epsilon$ = 0.2 - 0.3 at terminal time of our simulations. We suggest that these values can be considered as time-asymptotic values for the CR acceleration efficiency, since the time-dependent evolution of CR modified shocks has become approximately self-similar before the terminal time.

• Publications of The Korean Astronomical Society
• /
• v.30 no.2
• /
• pp.545-548
• /
• 2015

Most high energy cosmic rays (CRs) are thought to be produced by diffusive shock acceleration (DSA) in supernova remnants (SNRs) within the Galaxy. Plasma and MHD simulations have shown that the self-excitation of MHD waves and amplification of magnetic fields via plasma instabilities are an integral part of DSA for strong collisionless shocks. In this study we explore how plasma processes such as plasma instabilities and wave-particle interactions can affect the energy spectra of CR protons and electrons, using time-dependent DSA simulations of SNR shocks. We demonstrate that the time-dependent evolution of the shock dynamics, the self-amplified magnetic fields and $Alfv{\acute{e}nic$ drift govern the highest energy end of the CR energy spectra. As a result, the spectral cutoffs in nonthermal X-ray and ${\gamma}$-ray radiation spectra are regulated by the evolution of the highest energy particles, which are injected at the early phase of SNRs. We also find that the maximum energy of CR protons can be boosted significantly only if the scale height of the magnetic field precursor is long enough to contain the diffusion lengths of the particles of interests. Thus, detailed understandings of nonlinear wave-particle interactions and time-dependent DSA simulations are crucial for understanding the nonthermal radiation from CR acceleration sources.

• Journal of The Korean Astronomical Society
• /
• v.37 no.5
• /
• pp.405-412
• /
• 2004

Cosmological shocks form as an inevitable consequence of gravitational collapse during the large scale structure formation and cosmic-rays (CRs) are known to be accelerated at collisionless shocks via diffusive shock acceleration (DSA). We have calculated the evolution of CR modified shocks for a wide range of shock Mach numbers and shock speeds through numerical simulations of DSA in 1D quasi-parallel plane shocks. The simulations include thermal leakage injection of seed CRs, as well as pre-existing, upstream CR populations. Bohm-like diffusion is assumed. We show that CR modified shocks evolve to time-asymptotic states by the time injected particles are accelerated to moderately relativistic energies (p/mc $\ge$ 1), and that two shocks with the same Mach number, but with different shock speeds, evolve qualitatively similarly when the results are presented in terms of a characteristic diffusion length and diffusion time. We find that $10^{-4} - 10^{-3}$ of the particles passed through the shock are accelerated to form the CR population, and the injection rate is higher for shocks with higher Mach number. The CR acceleration efficiency increases with shock Mach number, but it asymptotes to ${\~}50\%$ in high Mach number shocks, regardless of the injection rate and upstream CR pressure. On the other hand, in moderate strength shocks ($M_s {\le} 5$), the pre-existing CRs increase the overall CR energy. We conclude that the CR acceleration at cosmological shocks is efficient enough to lead to significant nonlinear modifications to the shock structures.

• Journal of The Korean Astronomical Society
• /
• v.36 no.3
• /
• pp.111-121
• /
• 2003

In order to explore the cosmic ray acceleration at the cosmological shocks, we have performed numerical simulations of one-dimensional, plane-parallel, cosmic ray (CR) modified shocks with the newly developed CRASH (Cosmic Ray Amr SHock) numerical code. Based on the hypothesis that strong Alfven waves are self-generated by streaming CRs, the Bohm diffusion model for CRs is adopted. The code includes a plasma-physics-based 'injection' model that transfers a small proportion of the thermal proton flux through the shock into low energy CRs for acceleration there. We found that, for strong accretion shocks with Mach numbers greater than 10, CRs can absorb most of shock kinetic energy and the accretion shock speed is reduced up to $20\%$, compared to pure gas dynamic shocks. Although the amount of kinetic energy passed through accretion shocks is small, since they propagate into the low density intergalactic medium, they might possibly provide acceleration sites for ultra-high energy cosmic rays of $E\ll10^{18}eV$. For internal/merger shocks with Mach numbers less than 3, however, the energy transfer to CRs is only about $10-20\%$ and so nonlinear feedback due to the CR pressure is insignificant. Considering that intracluster medium (ICM) can be shocked repeatedly, however, the CRs generated by these weak shocks could be sufficient to explain the observed non-thermal signatures from clusters of galaxies.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.4
• /
• pp.719-728
• /
• 2017

This paper proposes a synthesis of new guidance law to generate an evasive maneuver against enemy's missile interception while considering its impact angle, acceleration, and field-of-view constraints. The first component of the synthesis is a new function of repulsive Artificial Potential Field to generate the evasive maneuver as a real-time dynamic obstacle avoidance. The terminal impact angle and terminal acceleration constraints compliance are based on Time-to-Go Polynomial Guidance as the second component. The last component is the Logarithmic Barrier Function to satisfy the field-of-view limitation constraint by compensating the excessive total acceleration command. These three components are synthesized into a new guidance law, which involves three design parameter gains. Parameter study and numerical simulations are delivered to demonstrate the performance of the proposed repulsive function and guidance law. Finally, the guidance law simulations effectively achieve the zero terminal miss distance, while satisfying an evasive maneuver against intercept missile, considering impact angle, acceleration, and field-of-view limitation constraints simultaneously.

• Journal of The Korean Astronomical Society
• /
• v.37 no.4
• /
• pp.249-255
• /
• 2004

Here we present a linear stability analysis and an MHD 2D model for the Parker-Jeans instability in the Galactic gaseous disk. The magnetic field is assumed parallel to a Galactic spiral arm, and the gaseous disk is modelled as a multi-component, magnetized, and isothermal gas layer. The model employs the observed vertical stratifications for the gas density and the gravitational acceleration in the Solar neighborhood, and the self-gravity of the gas is also included. By solving Poisson's equation for the gas density stratification, we determine the vertical acceleration due to self-gravity as a function of z. Subtracting it from the observed gravitational acceleration, we separate the total acceleration into self and external gravities. The linear stability analysis provides the corresponding dispersion relations. The time and length scales of the fastest growing mode of the Parker-Jeans instability are about 40 Myr and 3.3 kpc, respectively. In order to confirm the linear stability analysis, we have performed two-dimensional MHD simulations. These show that the Parker-Jeans instability under the self and external gravities evolves into a quasi-equilibrium state, creating condensations on the northern and southern sides of the plane, in an alternate manner.

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.5
• /
• pp.1243-1250
• /
• 1995

For the numerical solution of frictional dynamic contact problems, correct contact points and displacements are determined by iteratively reducing the displacement error vector monotonically toward zero And spurious oscillations are prevented from the solution by enforcing the velocity and acceleration compatibilities of the contact points with the corresponding error vectors. Numerical simulations are conducted to demonstrate the accuracy of the solution and the necessity of the velocity and acceleration compatibilities on the contact surface.