Coagulation of dust aggregates plays an important role in the formation of planets and is of key importance to the evolution of protoplanetary disks (PPDs). Characteristics of dust, such as the diversity of particle size, porosity, charge, and the manner in which dust couples to turbulent gas, affect the collision outcome and the rate of dust growth. Here we present a numerical model of the evolution of the dust population within a PPD which incorporates all of these effects. The probability that any two particles collide depends on the particle charge, cross-sectional area and their relative velocity. The actual collision outcome is determined by a detailed collision model which takes into account the aggregate morphology, trajectory, orientation, and electrostatic forces acting between charged grains. The data obtained in this research reveal the characteristics of dust populations in different environments at the end of the hit-and-stick growth, which establishes the foundation for the onset of the next growth stage where bouncing, mass transfer and fragmentation become important. For a given level of turbulence, neutral and weakly charged particles collide more frequently and grow faster than highly charged particles. However, highly charged particles grow to a larger size before reaching the bouncing barrier, and exhibit a "Runaway" growth, in which a few large particles grow quickly by accreting smaller particles while the rest of the population grows very slowly. In general, highly charged aggregates have a more compact structure and are comprised of larger monomers than neutral/weakly charged aggregates. The differences in the particle structure/composition not only affect the threshold velocities for bouncing and fragmentation, but also change the scattering and absorption opacity of dust, influencing the appearance of PPDs.

In order to characterize the early growth of fine-grained dust rims (FGRs) that commonly surround chondrules, we simulate the growth of FGRs through direct accretion of monomers of various sizes onto the chondrule surfaces. Dust becomes charged to varying degrees in the radiative plasma environment of the solar nebula (SN), and the resulting electrostatic force alters the trajectories of colliding dust grains, influencing the structure of the dust rim as well as the time scale of rim formation. We compare the growth of FGRs in protoplanetary disks (PPD) with different turbulence strengths and plasma conditions to previous models which assumed neutral dust grains. We use a combination of a Monte Carlo method and an N-body code to simulate the collision of dust monomers with chondrules: a Monte Carlo algorithm is used to randomly select dust particles that will collide with the chondrule as well as determine the elapsed time interval between collisions; at close approach, the detailed collision process is modeled using an N-body algorithm, Aggregate Builder (AB), to determine the collision outcome, as well as any restructuring of the chondrule rim. The collisions are driven by Brownian motion and coupling to turbulent gas motion in the protoplanetary disk. The charge distribution of the dust rim is modeled, used to calculate the trajectories of dust grains, and then analyze the resulting morphology of the dust rim. In a weakly turbulent region, the decreased relative velocity between charged particles causes small grains to be repelled from the chondrule, causing dust rims to grow more slowly and be composed of larger monomers, which results in a more porous structure. In a highly turbulent region, the presence of charge mainly affects the porosity of the rim by causing dust particles to deviate from the extremities of the rim and reducing the amount of restructuring.

In order to characterize the early growth of fine-grained dust rims (FGRs) that commonly surround chondrules, we perform numerical simulations of dust accretion onto chondrule surfaces. We employ a Monte Carlo algorithm to simulate the collision of dust monomers having radii between 0.5 and 10 $\mu$m with chondrules whose radii are between 500 and 1000 $\mu$m, in 100-$\mu$m increments. The collisions are driven by Brownian motion and solar nebula turbulence. After each collision, the colliding particles either stick at the point of contact, roll or bounce. We limit accretion of dust monomers (and in some cases, dust aggregates) to a small patch of the chondrule surface, for computational expediency. We model the morphology of the dust rim and the trajectory of the dust particle, which are not considered in most of the previous works. Radial profiles of FGR porosity show that rims formed in weak turbulence are more porous (with a porosity of 60-74\%) than rims formed in stronger turbulence (with a porosity of 52-60\%). The lower end of each range corresponds to large chondrules and the upper end to small chondrules, meaning that the chondrule size also has an impact on FGR porosity. The thickness of FGRs depends linearly on chondrule radius, and the slope of this linear dependency increases with time, and decreases with the turbulence strength. The porosity of FGRs formed by dust aggregates is $\sim 20\%$ on average greater than that of FGRs formed by single monomers. In general, the relatively high porosities that we obtain are consistent with those calculated by previous authors from numerical simulations, as well as with initial FGR porosities inferred from laboratory measurements of rimmed chondrule samples and rimmed chondrule analogs.

Though coronal mass ejections (CMEs) are magnetized fully-ionized gases, a recent observational study of a CME collision event in 2008 November has suggested that their behavior in the heliosphere is like elastic balls, and their collision is probably super-elastic \citep{Shen_etal_2012}. If this is true, this finding has an obvious impact on the space weather forecasting because the direction and veliocity of CMEs may change. To verify it, we numerically study the event through three-dimensional MHD simulations. The nature of CMEs' collision is examined by comparing two cases. In one case the two CMEs collide as observed, but in the other, they do not. Results show that the collision leads to extra kinetic energy gain by 3%--4% of the initial kinetic energy of the two CMEs. It firmly proves that the collision of CMEs could be super-elastic.

Reliable measurements of the solar magnetic field are still restricted to the photosphere, and our present knowledge of the three-dimensional coronal magnetic field is largely based on extrapolation from photospheric magnetogram using physical models, e.g., the nonlinear force-free field (NLFFF) model as usually adopted. Most of the currently available NLFFF codes have been developed with computational volume like Cartesian box or spherical wedge while a global full-sphere extrapolation is still under developing. A high-performance global extrapolation code is in particular urgently needed considering that Solar Dynamics Observatory (SDO) can provide full-disk magnetogram with resolution up to $4096\times 4096$. In this work, we present a new parallelized code for global NLFFF extrapolation with the photosphere magnetogram as input. The method is based on magnetohydrodynamics relaxation approach, the CESE-MHD numerical scheme and a Yin-Yang spherical grid that is used to overcome the polar problems of the standard spherical grid. The code is validated by two full-sphere force-free solutions from Low & Lou's semi-analytic force-free field model. The code shows high accuracy and fast convergence, and can be ready for future practical application if combined with an adaptive mesh refinement technique.