Numerical simulation of cooling coronal flare plasma

We have simulated the cooling of coronal flare plasma (T/sub e/>107 K) using a numerical model of a vertical magnetic flux tube containing an idealized flare chromosphere, transition region, and corona. The model solves the set of one-dimensional, two-fluid hydrodynamic equations. The cooling of the flux tube is calculated for a specific case beginning with an initial atmosphere in hydrostatic equilibrium and a maximum temperature of about 18 x 106 K. The behavior of temperature, density, and velocity is calculated as a function of height as the system cools. Early in the cooling energy is transported by conduction into the transition region and chromosphere where it is radiated away. Later, the transition region-corona interface moves upward into the tube at velocities of about 20 km s-1, while the chromosphere cools and the coronal component cools by both conduction and radiation. Coronal downflow velocities of about 60 km s-1 are evident during this phase. The expected spectral line emission from the system in X-ray lines of Fe XXV, Fe XXIV, Fe XXII, O VIII, and O VII is also calculated and compared to recent observational results. Some observational results can be explained as a consequence of simple cooling of flare flux tubes. The expected spectral line emission from certain transition region lines is also briefly considered. The dependence of our results on flare size is discussed, and our results are compared with similar previous work