[1] In the framework of integrated numerical space weather prediction, we have developed a 3-D MHD simulation model of the solar surface-solar wind system. We report the construction method of the model and its first results. By implementing a grid system with angularly unstructured and increasing radial spacing, we realized a spherical grid that has no pole singularity and realized a fine grid size around the inner boundary and a wide-range grid up to a size of 1 AU simultaneously. The magnetic field at the inner boundary is specified by the observational data. In order to obtain the supersonic solar wind speed, parameterized source functions are introduced into the momentum and energy equations. These source functions decay exponentially in altitude as widely used in previous studies.The absolute values of the source functions are controlled so as to reflect the topology of the coronal magnetic field. They are increased inside the magnetic flux tube with subradial expansion and reduced inside the magnetic flux tube with overradial expansion. This adjustment aims to reproduce the variation of the solar wind speed according to the coronal magnetic structure. The simulation simultaneously reproduces the plasma-exit structure, the high-and low-temperature regions, the open and closed magnetic field regions in the corona, the fast and slow solar wind, and the sector structure in interplanetary space. It is confirmed from the comparison with observations that the MHD model successfully reproduces many features of both the fine solar coronal structure and the global solar wind structure.

The heliospheric structure ranging from the solar surface to the earth's orbit is self-consistently reproduced from a time-stationary three-dimensional (3D) magnetohydrodynamic (MHD) simulation. The simulation model incorporates gravity, Coriolis, and centrifugal forces into the momentum equation, and coronal heating and fieldaligned thermal conduction into the energy equation. The heating term in the present model has its peak at 2.8 solar radius (R s ) and exponentially falls to zero at greater distance from the solar surface. The absolute value of heating depends on the topology of the solar magnetic field so as to be in inverse proportion with the magnetic expansion factor. The results of the simulation simultaneously reproduce the plasma-exit structure on the solar surface, the high-temperature region in the corona, the open-and closed-magnetic-field structures in the corona, the fast and slow streams of the solar wind, and the sector structure in the heliosphere.

No abstract