2021MNRAS.505.5238J

Thermal instability plays a crucial role in the formation of multiphase structures and their dynamics in the interstellar medium, and is a leading theory for cold cloud creation in various astrophysical environments. In this paper, we use 2D simulations to investigate thermal instability under the influence of various initial conditions and physical processes. We experiment with Gaussian random field (GRF) density perturbations of different initial power spectra. We also enrol thermal conduction and physical viscosity in isotropic hydrodynamic and anisotropic magnetohydrodynamic (MHD) simulations. We find that the initial GRF spectral index α has a dramatic impact on the pure hydrodynamic development of thermal instability, influencing the size, number, and motions of clouds. Cloud fragmentation happens due to two mechanisms: tearing and contraction rebound. In the runs with isotropic conduction and viscosity, the structures and dynamics of the clouds are dominated by evaporation and condensation flows in the non-linear regime, and the flow speed is regulated by viscosity. Cloud disruptions happen as a result of the Darrieus-Landau instability. However, at very late times, all individual clouds merge into one cold structure in all hydrodynamic runs. In the MHD case, the cloud structure is determined by both the initial perturbations and the initial magnetic field strength. In high-β runs, anisotropic conduction causes dense filaments to align with the local magnetic fields and the field direction can become reoriented. Strong magnetic fields suppress cross-field contraction and cold filaments can form along or perpendicular to the initial fields.