The thermal stability of coronal loops - Numerical simulations

The radiative stability of thermally isolated coronal loops with free-flow boundary conditions by nonlinear numerical simulation are studied. A chromosphere-to-corona loop equilibrium (including the option of a deep chromosphere) is first established by following the nonlinear evolution from an initial isothermal state with rigid boundaries. The end conditions are then changed, to allow free flow and to fix the temperature, and investigate the response to nonisobaric perturbations. Within a family of loops of the same pressure, long hot loops to be stable and short cool loops to be unstable to the thermal chromospheric-expansion mode are found. The stable cases remain so, even when long chromospheric ends and/or gravity are added. In those cases which are unstable, the subsequent nonlinear evolution which exhibits swelling of the chromosphere until the entire loop becomes cool and dense are followed.