Numerical Simulation of Longitudinal Oscillation of Filament Based on Smoothed Particle Hydrodynamics (SPH) Method
Yu Xiang Liu,Hong Fu Qiang,Xue Ren Wang,Du Dou Wang,Yan Chao Wang,Lin Tao Zhang
DOI: https://doi.org/10.1007/s11207-022-01999-2
IF: 2.961
2022-06-01
Solar Physics
Abstract:The oscillations of solar coronal filaments have received significant interest, and the restoring force and damping mechanism are the focus of discussion. We study the basic physical principles of the multi-factor, coupled, longitudinal oscillations of a filament, including its dominant restoring force and the damping mechanism. Using the Smoothed Particle Hydrodynamics (SPH) method we have performed for the first time a one-dimensional (1D) hydrodynamic numerical simulation of the longitudinal oscillation of the prominence. We modeled the process of mass uplift, that is the chromosphere mass enters the corona directly with the uplifted magnetic-flux tube. We studied multiple parameters, including the mass density, the initial perturbation velocity, and the geometry of the magnetic depression. We found that a sudden velocity disturbance can result in the oscillation of the prominence mass along the magnetic dip. Our numerical simulation analysis and theoretical study show that the oscillation period can be expressed as (where is the dip depth of the magnetic tube, is half the length of the prominence magnetic dip, and is the initial velocity disturbance), which is slightly different from the previous linear pendulum theory. The particle-collision viscosity also affects the longitudinal-oscillation period. In a non-adiabatic process, the oscillation decays with time. Radiant cooling is the primary factor that causes damping. Plasma particles' collisional viscosity is also a cause of damping, which primarily occurs once the prominence moves to the farthest point from the equilibrium state. The proportional law of the attenuation time scale is deduced, that is (where is the prominence length, and is the prominence density). Compared to previous simulations, plays a more important role. The greater the amplitude, the greater the oscillation damping. The advantage of the SPH method used in this article is that it can avoid the accuracy loss caused by mesh distortion in extremely large deformations, and it can avoid the interface problem of the Euler mesh in the Euler method. Therefore, it is especially suitable for solving dynamic large deformation problems such as prominence oscillation and explosion.
astronomy & astrophysics