Science & Research
Research Highlights > Theory & Computational Astrophysics

General Info

Research

People/Visitors

Activities

Public Outreach

Links

Home
Job Openings
Internal Document
Web Mail

 

 

 

 

Magnetohydrodynamics by Dr. Hsien Shang

Zeus 2D Simulations of Wide-angle winds

We (Hsien Shang together with research assistant Mei-Yin Chou and Anthony Allen) are carrying out studies of MHD flows by Zeus2D. Simulations of wide-angle winds are done on their asymptotic regime and for long-term evolution to achieve time-steady state, and to further investigate the implications on the formation of outflows by interacting these wide-angle winds with ambient media. They adopted analytic construction for the MHD wind where the density, poloidal and toroidal velocity, and toroidal magnetic fields follow a simple power-law relationship from Li & Shu (1996). Due to the sharp drop-off in density and magnetic field strengths characterized by such winds, and large spatial coverage in the computational domains, establishing robust numerical results has been a challenging task. The numerical boundary and reflection effects arise very easily. To establish robust runs and performances based on local resources, they first established benchmark tests on various platforms of locally accessible ASIAA machines and with the use of a variety of optimization flags, the benchmark is summarized in Allen, Chou, and Shang (2002). Steady-state runs have also been established for super-Alfv\'enic winds running into vacuum from a narrow launch sphere of hidden zones in the hope to mimic the x-winds.

On the effects of wide-angle winds interacting with ambient environment, in progress are the simulations of winds running into surrounding medium of constant density, and later on more realistic configurations from Li and Shu (1997) toroid model sequences may be implemented. In the figures are simulations of winds of terminal speed of 300 km/s (~ 3 Va) and ambient densities from 10-16 and 10-19g/cm3. These results are plotted on four-panels to illustrate selected time frames of 100, 400, 700, and 1000 years of evolution. The appearance of the resultant density structures and shapes of shock fronts suggests strong connections to the lobe shapes of molecular outflows. Further work will continue to investigate such connections.

The ambient density is 2x10-16 gm/cm3 (top) and 2x10-19 gm/cm3 respectively. The running time from left to right, top to bottom of each picture are 100, 400, 700 and 1000 years accordingly. We divide the r-axis zone into 200 grids and z-axis zone into 1200 grids uniformly.

Topics of Magnetohydrodynamics:

Simulation of Magnetic Collapse
Zeus 2D Simulations of Wide-angle winds

<< Back to Theory Index
 
HOME INDEX MAIL CHINESE