Clarify the role of visco-resistive heating in saturated SMAP dynamics

Clarify whether visco-resistive heating in the energy-conserving full-MHD model affects the numerical stability and quasi-linearly saturated dynamics of slow magnetoacoustic pulsations.

Background

The model converts dissipated wave energy into thermal pressure through visco-resistive heating. Removing these heating terms leaves the first few unsaturated SMAP pulses largely unchanged, but the modified simulation subsequently develops a rapidly diverging fluctuation energy and terminates abnormally.

Because disabling heating also violates energy conservation, the authors cannot determine whether the late-time behavior reflects a physical role of heating, a numerical instability, or both. They explicitly identify the effect on numerical stability and saturated dynamics as unresolved.