We develop a new numerical scheme for ideal magnetohydrodynamic (MHD) simulations, which is robust against one- and multi-dimensional shocks, and is accurate for low Mach number flows and discontinuities. The scheme belongs to a family of the advection upstream splitting method employed in computational aerodynamics, and it splits the inviscid flux in MHD equations into advection, pressure, and magnetic tension parts, and then individually evaluates mass, pressure, and magnetic tension fluxes at the interface of a computational cell. The mass flux is designed to avoid numerical shock instability in multidimension, while preserving contact discontinuity. The pressure flux possesses a proper scaling for low Mach number flows, allowing reliable simulations of nearly incompressible flows. The magnetic tension flux is built to be consistent with the HLLD approximate Riemann solver to preserve rotational discontinuity. We demonstrate various benchmark tests to verify the novel performance of the scheme. Our results indicate that the scheme must be a promising tool to tackle astrophysical systems that include both low and high Mach number flows, as well as magnetic field inhomogeneities.

We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the HLLD approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibited a convergent picture of the early nonlinear growth of the MRI in a global model, and indicated that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.

We design a conservative finite difference scheme for ideal magnetohydrodynamic simulations that attains high-order accuracy, shock-capturing, and divergence-free condition of the magnetic field. The scheme interpolates pointwise physical variables from computational nodes to midpoints through a high-order nonlinear weighted average. The numerical flux is evaluated at the midpoint by a multi-state approximate Riemann solver for correct upwinding, and its spatial derivative is approximated by a high-order linear central difference to update the variables with designed order of accuracy and conservation. The magnetic and electric fields are defined at staggered grid points employed in the Constrained Transport (CT) method by Evans & Hawley (1988). We propose a new CT variant, in which the staggered electric field is evaluated so as to be consistent with the base one-dimensional Riemann solver and the staggered magnetic field is updated to be divergence-free as designed high-order finite difference representation. We demonstrate various benchmark tests to measure the performance of the present scheme. We discuss the effect of the choice of interpolation methods, Riemann solvers, and the treatment for the divergence-free condition on the quality of numerical solutions in detail.

We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD Runge-Kutta integration. We apply our spatially fifth- and seventh-order schemes to a suite of simulations of collisionless self-gravitating systems and electrostatic plasma simulations, including linear and nonlinear Landau damping in one dimension and Vlasov--Poisson simulations in a six-dimensional phase space. The high-order schemes achieve a significantly improved accuracy in comparison with the third-order positive-flux-conserved scheme adopted in our previous study. With the semi-Lagrangian time integration, the computational cost of our high-order schemes does not significantly increase, but remains roughly the same as that of the third-order scheme. Vlasov--Poisson simulations on $128^3 \times 128^3$ mesh grids have been successfully performed on a massively parallel computer.

Nonlinear evolution of magnetic reconnection is investigated by means of magnetohydrodynamic simulations including uniform resistivity, uniform viscosity, and anisotropic thermal conduction. When viscosity exceeds resistivity (the magnetic Prandtl number Prm > 1), the viscous dissipation dominates outflow dynamics and leads to the decrease in the plasma density inside a current sheet. The low-density current sheet supports the excitation of the vortex. The thickness of the vortex is broader than that of the current for Prm > 1. The broader vortex flow more efficiently carries the upstream magnetic flux toward the reconnection region, and consequently boosts the reconnection. The reconnection rate increases with viscosity provided that thermal conduction is fast enough to take away the thermal energy increased by the viscous dissipation (the fluid Prandtl number Pr < 1). The result suggests the need to control the Prandtl numbers for the reconnection against the conventional resistive model.

A large set of numerical simulations of magnetohydrodynamic (MHD) turbulence induced by the magnetorotational instability (MRI) is presented. Revisiting the previous survey conducted by Sano et al. (2004), we investigate the gas pressure dependence of the saturation level. In ideal MHD simulations, the gas pressure dependence is found to be very sensitive to the choice of a numerical scheme. This is because the numerical magnetic Prandtl number varies according to the scheme as well as the pressure, which considerably affects the results. The saturation level is more sensitive to the numerical magnetic Prandtl number than the pressure. In MHD simulations with explicit viscosity and resistivity, the saturation level increases with the physical magnetic Prandtl number, and it is almost independent of the gas pressure when the magnetic Prandtl number is constant. This is indicative of the incompressible turbulence saturated by the secondary tearing instability.

We present an extension of the multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

We present a new numerical scheme for solving the advection equation and its application to the Vlasov simulation. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions compared to other existing schemes with the same memory usage. The two-dimensional scheme can solve the solid body rotation problem of a gaussian profile with little numerical diffusion. This is a very important property for Vlasov simulations of magnetized plasma. The application of the scheme to the electromagnetic Vlasov simulation of collisionless shock waves is presented as a benchmark test.

We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, for better conservation of the information entropy. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions within reasonable usage of computational resources compared to other existing schemes. The two-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile for more than hundred rotation periods with little numerical diffusion. This is crucially important for Vlasov simulations of magnetized plasmas. Applications of the one- and two-dimensional schemes to electrostatic and electromagnetic Vlasov simulations are presented with some benchmark tests.

Using a model of particle acceleration and transport in solar flares, we investigate the height distribution of coronal electrons by focusing on the energy-dependent pitch-angle scattering. When pitch-angle scattering is not included, the peak heights of loop-top electrons are constant, regardless of their energy, owing to the continuous acceleration and compression of the electrons via shrinkage of magnetic loops. On the other hand, under pitch-angle scattering, the electron heights are energy dependent; intermediate energy electrons are at a higher altitude, whereas lower and higher energy electrons are at lower altitudes. This implies that the intermediate energy electrons are inhibited to follow the shrinking field lines to lower altitudes because pitch-angle scattering causes efficient precipitation of these electrons into the footpoint and their subsequent loss from the loop. This result is qualitatively consistent with the position of the above-the-loop-top hard X-ray (HXR) source that is located above coronal HXR loops emitted by lower energy electrons and microwaves emitted by higher energy electrons. Quantitative agreement with observations might be achieved by considering primary acceleration before the onset of loop shrinkage and additional pitch-angle scattering via wave-particle interactions.

Based on the drift-kinetic theory, we develop a model for particle acceleration and transport in solar flares. The model describes the evolution of the particle distribution function by means of a numerical simulation of the drift-kinetic Vlasov equation, which allows us to directly compare simulation results with observations within an actual parameter range of the solar corona. Using this model, we investigate the time evolution of the electron distribution in a flaring region. The simulation identifies two dominant mechanisms of electron acceleration. One is the betatron acceleration at the top of closed loops, which enhances the electron velocity perpendicular to the magnetic field line. The other is the inertia drift acceleration in open magnetic field lines, which produces antisunward electrons. The resulting velocity space distribution significantly deviates from an isotropic distribution. The former acceleration can be a generation mechanism of electrons that radiate loop-top nonthermal emissions, and the latter be of escaping electrons from the Sun that should be observed by in-situ measurements in interplanetary space and resulting radio bursts through plasma instabilities.

We present a multi-wavelength observation of a solar flare occurring on 2006 December 13 with Hinode, RHESSI, and the Nobeyama Radio Observatory, to study the electron acceleration site and mechanism. The Solar Optical Telescope (SOT) on board Hinode observed elongated flare ribbons, and RHESSI observed double-footpoint hard X-ray (HXR) sources appearing in part of the ribbons. A photospheric vector magnetogram obtained from SOT reveals that the HXR sources are located at the region where horizontal magnetic fields change the direction. The region is interpreted as the footpoint of magnetic separatrix. Microwave images taken with the Nobeyama Radioheliograph show a loop structure connecting the HXR sources. The brighter parts of the microwave intensity are located between the top and footpoints of the loop. We consider these observations as an evidence of the electron acceleration near the magnetic separatrix and injection parallel to the field line.

We analytically and numerically study the motion of electrons along a magnetic loop, to compare with the observation of the propagating feature of the non-thermal microwave source in the 1999 August 28 solar flare reported by Yokoyama et al. (2002). We model the electron motion with the Fokker-Planck equation and calculate the spatial distribution of the gyrosynchrotron radiation. We find that the microwave propagating feature does not correspond to the motion of electrons with a specific initial pitch angle. This apparent propagating feature is a consequence of the motion of an ensemble of electrons with different initial pitch angles, which have different time and position to produce strong radiation in the loop. We conclude that the non-thermal electrons in the 1999 August 28 flare were isotropically accelerated and then are injected into the loop.

The Solar Optical Telescope (SOT) on board Hinode satellite observed an X3.4 class flare on 2006 December 13. Typical two-ribbon structure was observed, not only in the chromospheric CaII H line but also in G-band and FeI 6302A line. The high-resolution, seeing-free images achieved by SOT revealed, for the first time, the sub-arcsec fine structures of the "white light" flare. The G-band flare ribbons on sunspot umbrae showed a sharp leading edge followed by a diffuse inside, as well as previously known core-halo structure. The underlying structures such as umbral dots, penumbral filaments and granules were visible in the flare ribbons. Assuming that the sharp leading edge was directly heated by particle beam and the diffuse parts were heated by radiative back-warming, we estimate the depth of the diffuse flare emission using the intensity profile of the flare ribbon. We found that the depth of the diffuse emission is about 100 km or less from the height of the source of radiative back-warming. The flare ribbons were also visible in the Stokes-V images of FeI 6302A, as a transient polarity reversal. This is probably related to "magnetic transient" reported in the literature. The intensity increase in Stokes-I images indicates that the FeI 6302A line was significantly deformed by the flare, which may cause such a magnetic transient.

We study the non-thermal emissions in a solar flare occurring on 2003 May 29 by using RHESSI hard X-ray (HXR) and Nobeyama microwave observations. This flare shows several typical behaviors of the HXR and microwave emissions: time delay of microwave peaks relative to HXR peaks, loop-top microwave and footpoint HXR sources, and a harder electron energy distribution inferred from the microwave spectrum than from the HXR spectrum. In addition, we found that the time profile of the spectral index of the higher-energy ($\gsim 100$ keV) HXRs is similar to that of the microwaves, and is delayed from that of the lower-energy ($\lsim 100$ keV) HXRs. We interpret these observations in terms of an electron transport model called {\TPP}. We numerically solved the spatially-homogeneous {\FP} equation to determine electron evolution in energy and pitch-angle space. By comparing the behaviors of the HXR and microwave emissions predicted by the model with the observations, we discuss the pitch-angle distribution of the electrons injected into the flare site. We found that the observed spectral variations can qualitatively be explained if the injected electrons have a pitch-angle distribution concentrated perpendicular to the magnetic field lines rather than isotropic distribution.