Unsteady, thermally conductive coronal flow

A numerical algorithm using implicit time-differencing is applied to the solar wind equations allowing, for the first time, solutions including thermal conduction to be found for time-dependent flow traversing the subsonic to supersonic velocity transition region. Sample solutions are shown that demonstate the distinctive differences introduced by including thermal conduction, in comparison to the commonly available solutions assuming polytropic flow. Also, it is found that steady solutions are produced at least as quickly using a time-dependent relaxation to the steady state as when solving the steady-state equations.