The objective of this paper is to explore the application of a six-component overset grid to solar wind simulation with a three-dimensional (3D) Solar-InterPlanetary Conservation Element/Solution Element MHD model. The essential focus of our numerical model is devoted to dealing with: (1) the singularity and mesh convergence near the poles via the use of the six-component grid system, (2) the ∇ · B constraint error via an easy-to-use cleaning procedure by a fast multigrid Poisson solver, (3) the Courant-Friedrichs-Levy number disparity via the Courant-number insensitive method, (4) the time integration by multiple time stepping, and (5) the time-dependent boundary condition at the subsonic region by limiting the mass flux escaping through the solar surface. In order to produce fast and slow plasma streams of the solar wind, we include the volumetric heating source terms and momentum addition by involving the topological effect of the magnetic field expansion factor f S and the minimum angular distance θ b (at the photosphere) between an open field foot point and its nearest coronal hole boundary. These considerations can help us easily code the existing program, conveniently carry out the parallel implementation, efficiently shorten the computation time, greatly enhance the accuracy of the numerical solution, and reasonably produce the structured solar wind. The numerical study for the 3D steady-state background solar wind during Carrington rotation 1911 from the Sun to Earth is chosen to show the above-mentioned merits. Our numerical results have demonstrated overall good agreements in the solar corona with the Large Angle and Spectrometric Coronagraph on board the Solar and Heliospheric Observatory satellite and at 1 AU with WIND observations.

In this paper, the spacetime conservation element and solution element (CESE ) method is applied to threedimensional magnetohydrodynamics ( MHD) equations in Cartesian coordinates for solar wind plasma, with the purpose of modeling the steady state solar atmospheric study. To illustrate this newly developed scheme we have studied two examples: (1) two-dimensional coronal dynamical structure with multipole magnetic fields and (2) threedimensional coronal dynamical structure, using measured solar surface magnetic fields and the empirical values of the plasma properties on the solar surface as the initial conditions for the set of MHD equations and then the relaxation method to achieve a quasiYsteady state. From these examples we have shown that the newly developed modified spacetime CESE scheme possesses the ability to model the Sun-Earth environment and other astrophysical flows.