• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.275-279
• /
• 2001

The study of incompressible magnetohydrodynamic (MHD) turbulence gives useful insights on many astrophysical problems. We describe a pseudo-spectral MHD code suitable for the study of incompressible turbulence. We review our recent' works on direct three-dimensional numerical simulations for MHD turbulence in a periodic box. In those works, we use a pseudo-spectral code to solve the incompressible MHD equations. We first discuss the structure and properties of turbulence as functions of scale. The results are consistent with the scaling law recently proposed by Goldreich & Sridhar. The scaling law is based on the concept of scale-dependent isotropy: smaller eddies are more elongated than larger ones along magnetic field lines. This scaling law substantially changes our views on MHD turbulence. For example, as noted by Lazarian & Vishniac, the scaling law can provide a fast reconnection rate. We further discuss how the study of incompressible MHD turbulence can help us to understand physical processes in interstellar medium (ISM) by considering imbalanced cascade and viscous damped turbulence.

• Journal of The Korean Astronomical Society
• /
• v.40 no.4
• /
• pp.165-169
• /
• 2007

The numerical simulations associated with the interaction of High Velocity Clouds (HVC) with the Magnetized Galactic Interstellar Medium (ISM) are a powerful tool to describe the evolution of the interaction of these objects in our Galaxy. In this work we present a new project referred to as Theoretical Virtual i Observatories. It is oriented toward to perform numerical simulations in real time through a Web page. This is a powerful astrophysical computational tool that consists of an intuitive graphical user interface (GUI) and a database produced by numerical calculations. In this Website the user can make use of the existing numerical simulations from the database or run a new simulation introducing initial conditions such as temperatures, densities, velocities, and magnetic field intensities for both the ISM and HVC. The prototype is programmed using Linux, Apache, MySQL, and PHP (LAMP), based on the open source philosophy. All simulations were performed with the MHD code ZEUS-3D, which solves the ideal MHD equations by finite differences on a fixed Eulerian mesh. Finally, we present typical results that can be obtained with this tool.

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

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.341-343
• /
• 2001

Here were continue the MHD study started by Santillan et al (1999) for the interaction of high-velocity clouds (HVCs) with the magnetized thick gaseous disk of our Galaxy. We use the MHD code ZEUS-3D and perform 3D-numerical simulations of this interaction, and study the formation of head-tail structures in HVCs. Our results show that clouds located above 2 kpc from mindplane present velocity and column density gradients with a cometary structure that is similar to those observed in 21 cm emission

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

High-velocity clouds are flows of neutral hydrogen, located at high galactic latitudes, with large velocities ($[VLSR]{\ge} 100 km/s$) that do not match a simple model of circular rotation for our Galaxy. Numerical simulations have been performed for the last 20 years to study the details of their evolution, and their possible interaction with the Galactic disk. Here we present a brief review of the models that have been already published, and describe newly performed three-dimensional magnetohydrodynamic simulations.

• 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.43 no.1
• /
• pp.62.1-62.1
• /
• 2018

Galaxies in the cluster environment interact with the intracluster medium (ICM), losing the interstellar medium (ISM) and alternating their evolution. Observational evidences of the extraplanar ISM stripped by the ICM's ram pressure are prevalent in HI imaging studies of cluster galaxies. However, current theoretical understanding of the ram pressure stripping (or ICM-ISM interaction in general) is still limited mainly due to the lack of numerical resolution at ISM scales in large-scale simulations. Especially, self-consistent modeling of the turbulent, multiphase ISM is critical to understand star formation in galaxies interacting with the ICM. To achieve this goal, we utilize the TIGRESS simulation suite, simulating a local patch of galactic disks with high resolution to resolve key physical processes in the ISM, including cooling/heating, self-gravity, MHD, star formation, and supernova feedback. We then expose the ISM disk to ICM flows and investigate the evolution of star formation rate and the properties of the ISM. By exploring ICM parameter space, we discuss an implication of the simple ram pressure stripping condition (so called the Gunn-Gott condition) to the realistic ISM.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.181-190
• /
• 2001

The advent of robust, reliable and accurate higher order Godunov schemes for many of the systems of equations of interest in computational astrophysics has made it important to understand how to solve them in multi-scale fashion. This is so because the physics associated with astrophysical phenomena evolves in multi-scale fashion and we wish to arrive at a multi-scale simulational capability to represent the physics. Because astrophysical systems have magnetic fields, multi-scale magnetohydrodynamics (MHD) is of especial interest. In this paper we first discuss general issues in adaptive mesh refinement (AMR), We then focus on the important issues in carrying out divergence-free AMR-MHD and catalogue the progress we have made in that area. We show that AMR methods lend themselves to easy parallelization. We then discuss applications of the RIEMANN framework for AMR-MHD to problems in computational astophysics.

• The Bulletin of The Korean Astronomical Society
• /
• v.38 no.1
• /
• pp.47.1-47.1
• /
• 2013

Supersonic turbulence is believed to decay rapidly within a flow crossing time irrespective of the degree of magnetization. However, this consensus of decaying magnetohydrodynamic (MHD) turbulence relies on local isothermal simulations, which are unable to investigate the role of global magnetic fields and structures. Utilizing three-dimensional MHD simulations including interstellar cooling and heating, we investigate decaying MHD turbulence within cold neutral medium sheets embedded in warm neutral medium. Early evolution is consistent with previous studies characterized rapid decay of turbulence with the decaying time shorter than a flow crossing time and power-law temporal decay of turbulent kinetic energy with slope of -1. If initial magnetic fields are strong and perpendicular to the sheet, however long term evolutions of kinetic energy shows that a significant amount of turbulent energy still remains even after ten flow crossing times, and decaying rate is reduced as field strengths increase. We analyse power spectra of remaining turbulence to show that incompressible, in-plane motions dominate.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.333-335
• /
• 2001

We study the properties of supernova (SN) driven interstellar turbulence with a numerical magnetohydrodynamic (MHD) model. Calculations were done using the RIEMANN framework for MHD, which is highly suited for astrophysical flows because it tracks shocks using a Riemann solver and ensures pressure positivity and a divergence-free magnetic field. We start our simulations with a uniform density threaded by a uniform magnetic field. A simplified radiative cooling curve and a constant heating rate are also included. In this radiatively-cooling magnetized medium, we explode SNe one at a time at randomly chosen positions with SN explosion rates equal to and 12 times higher than the Galactic value. The evolution of the system is basically determined by the input energy of SN explosions and the output energy of radiative cooling. We follow the simulations to the point where the total energy of the system, as well as thermal, kinetic, and magnetic energy individually, has reached a quasi-stationary value. From the numerical experiments, we find that: i) both thermal and dynamical processes are important in determining the phases of the interstellar medium, and ii) the power index n of the $B-p^n$ relation is consistent with observed values.