On the Modeling of Kink Oscillations in Fine-Structured Coronal Loops with Field-aligned Nonlinear Longitudinal Disturbances

This record contains the ideal-MHD simulation data, input decks, and analysis scripts underlying the study of standing kink oscillations in fine-structured, cool coronal loop strands subject to field-aligned nonlinear longitudinal disturbances. Two-dimensional ideal-MHD simulations were performed with the ATHENA code (Gardiner & Stone 2005) in the x–z plane, using a linearized Riemann solver and constrained transport, on a 300 × 400 uniform grid covering (0, 0.91L) × (−0.5L, 0.5L) with L = 100 Mm and open boundaries. An isobaric strand of width a = 0.4 Mm and density contrast d = 5 or d = 20 is embedded in a loop whose footpoints are anchored in a dense chromospheric-like layer (dp = 10³, sp = 2 Mm) at x = 0.10L and x = 0.81L. Kink oscillations are excited impulsively by a Gaussian Vz pulse of amplitude 0.15 vA0 and width 10 Mm applied at (x0, z0) = (0.455L, −0.1L). Field-aligned flows of 92 km/s (matching Hinode/SOT observations) are imposed as an initial condition, not as a driven steady state, in three geometries — bounded (strand interior only), unbounded (interior plus ambient corona), and external (ambient corona only) — in both uniform and longitudinally varying (Gaussian) forms, together with runs sampling the internal/external amplitude ratio Vout/Vin. The archive provides: ATHENA input files for every run; simulation output and reduced time series of mass density, longitudinal velocity Vx and transverse velocity Vz; derived quantities from the Gaussian strand fits (centre position, half-width) and the damped-sine fits (amplitude, period P, damping time τ, quality factor τ/P); and the Python and IDL scripts used for the fitting and for generating Figures 1–12 of the associated paper.