AGB Stellar Winds

I was testing the program on bluehive2.

This is the wind vprofile when it reaches steady state. However, it is very weird that when I use color plot, it shows that the lowest velocity is 0.005691 which is exactly the boundary condition I specified. Besides that, a considerable region is below unity, but in the lineout plot, the lowest velocity is above unity.

The boundary condition would agree with the color plot, but why there is a difference?

The asymptotic behavior of the wind is very close to Parker's solution (

comparing the final velocity, the relative diff is only 5%

), the global quantitative difference has not been determined yet since there maybe a mistake in the lineout plot.

They are not very symmetric…

Attachments (2)

Download all attachments as: .zip

Comments

No comments.