{'fraction_non_alphanumeric': 0.037133097165991905, 'fraction_numerical': 0.021919281376518218, 'mean_word_length': 4.388955173001534, 'pattern_counts': {'":': 0, '<': 12, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
How flux feeding causes eruptions of solar magnetic flux ropes with the hyperbolic flux tube configuration?

Coronal magnetic flux ropes are generally considered to be the core structure of large-scale solar eruptions. Recent observations found that solar eruptions could be initiated by a sequence of"flux feeding,"during which chromospheric fibrils rise upward from below, and merge with a pre-existing prominence. Further theoretical study has confirmed that the flux feeding mechanism is efficient in causing the eruption of flux ropes that are wrapped by bald patch separatrix surfaces. But it is unclear how flux feeding influences coronal flux ropes that are wrapped by hyperbolic flux tubes (HFT), and whether it is able to cause the flux-rope eruption. In this paper, we use a 2.5-dimensional magnetohydrodynamic model to simulate the flux feeding processes in HFT configurations. It is found that flux feeding injects axial magnetic flux into the flux rope, whereas the poloidal flux of the rope is reduced after flux feeding. Flux feeding is able to cause the flux rope to erupt, provided that the injected axial flux is large enough so that the critical axial flux of the rope is reached. Otherwise, the flux rope system evolves to a stable equilibrium state after flux feeding, which might be even farther away from the onset of the eruption, indicating that flux feeding could stabilize the rope system with the HFT configuration in this circumstance.


Introduction
Large-scale solar eruptions include prominence/filament eruptions, flares, and coronal mass ejections (CMEs) (Benz 2008;Chen 2011;Parenti 2014;Liu 2020). They are capable of inflicting huge impacts on the solar-terrestrial system (Švestka 2001;Cheng et al. 2014;Shen et al. 2014;Lugaz et al. 2017;Gopalswamy et al. 2018). It is widely accepted that different kinds of large-scale solar eruptions are close related to each other: they are essentially different manifestations of the same eruptive process of a coronal magnetic flux rope system (Zhang et al. 2001;Vršnak et al. 2005;van Driel-Gesztelyi & Green 2015;Jiang et al. 2018;Liu et al. 2018;Yan et al. 2020). Therefore, it is of great significance to investigate how the eruption of coronal magnetic flux ropes is initiated. According to the magnetic topology, coronal flux ropes are classified into two types of configurations: if the flux rope sticks to the photosphere, with a bald patch separatrix surface (BPSS, Titov et al. 1993;Titov & Démoulin 1999;Gibson & Fan 2006) wrapping the flux rope, this is usually called the BPS configuration (Filippov 2013); for the flux rope system in which the rope is suspended in the corona and wrapped around by a hyperbolic flux tube (HFT), it is called the HFT configuration Aulanier et al. 2005;Chintzoglou et al. 2017).
Many theoretical analyses have been carried out to investigate the eruptive mechanism of coronal magnetic flux ropes. Apart from the well-known magnetohydrodynamic (MHD) instabilities (Romano et al. 2003;Török & Kliem 2003;Kliem & Török 2006;Fan & Gibson 2007;Aulanier et al. 2010;Guo et al. 2010;Savcheva et al. 2012), it was also suggested by many previous studies that catastrophes could be responsible for solar eruptions (e.g., Forbes & Isenberg 1991;Isenberg et al. 1993;Lin et al. 2001;Chen et al. 2007;Longcope & Forbes 2014;Kliem et al. 2014). For example, if either the axial flux or the poloidal flux of a flux rope exceeds the corresponding critical value, a catastrophic loss of equilibrium will occur, resulting in the eruption of the rope (Li & Hu 2002;Bobra et al. 2008;Su et al. 2011;Zhang et al. 2016;Zhuang et al. 2018). Alternatively, magnetic reconnection may play a dominant role in flux rope eruptions, such as break-out model (Antiochos et al. 1999;Sterling & Moore 2004), tether-cutting model (Moore et al. 2001;Inoue et al. 2015), and flux emergence model (Chen & Shibata 2000;Archontis & Hood 2008). Recently, Zhang et al. (2014) observed that chromospheric fibrils rise from below a quiescent prominence and merge with this prominence. This phenomenon is called "flux feeding" process. As observed by Zhang et al. (2014), the flux feeding processes occur 3 times, followed by the A&A proofs: manuscript no. ams eruption of the quiescent prominence. This implies that coronal eruptions could be initiated by flux feeding.
Since there lacks accurate measurement of the local conditions in the corona, the physical essence of flux feeding mechanism is unclear based on observational results alone. To solve this, Zhang et al. (2020) (hereafter Paper I) carried out numerical simulations to investigate the scenario of the flux feeding mechanism in BPS configurations. In their simulations, the rising chromospheric fibril is represented by a small flux rope emerging from below a pre-existing coronal magnetic flux rope with a BPS configuration. It was found that flux feeding processes only inject axial magnetic flux into the pre-existing flux rope, whose poloidal flux, however, remains unchanged after flux feeding. If the amount of the injected axial flux is large enough so that the total axial flux of the flux rope exceeds its critical value, the eruption of the rope is initiated by the upward catastrophe associated with the critical axial flux; otherwise, the flux rope remains in the stable BPS configuration after flux feeding. Different from the BPS configurations, HFT configurations are usually considered as the pre-eruptive states of coronal eruptions Aulanier et al. 2010;Filippov 2013), indicating that HFT configurations could be metastable or even unstable. Here come the questions: how does flux feeding influence the flux rope systems with an HFT configuration? Is flux feeding still able to cause eruptions? If flux feeding mechanism is still efficient, what are the similarities and differences between the BPS and the HFT cases? To resolve these questions, theoretical analyses about flux feeding in HFT configurations are needed, which will further shed light on the physical nature of the flux feeding mechanism.
In this paper, we carry out 2.5-dimensional numerical simulations to investigate the influence of flux feeding on flux rope system with an HFT configuration, especially how flux feeding causes the rope's eruption. The rest of the paper is arranged as follows: the simulating procedures are introduced in Sect. 2; the simulation results of a typical flux feeding process and the flux rope eruption caused by this process are presented in Sect. 3; the influence of flux feeding on the rope system and the physical scenario of flux feeding mechanism in HFT configurations are analyzed in Sect. 4. Finally, discussion and conclusion are given in Sect. 5.

Simulating procedures
For the 2.5-dimensional cases in our simulations, all the magnitudes satisfy ∂/∂z = 0, thus the magnetic field can be denoted as Here B z is the component of the magnetic field in z−dirention; ψ is the magnetic flux function, and the isolines of ψ correspond to the magnetic field lines projected in x − y plane. In this paper, the multi-step implicit scheme (Hu 1989) is used to simulate the evolution of coronal flux rope systems. The initial state in our simulation is obtained by numerical procedures: first we insert a flux rope into a potential background field from the lower base, and then by adjusting the properties of this rope, a flux rope system with an HFT configuration is eventually obtained. The basic equations and the detailed simulating procedures are introduced in Appendix A. The magnetic configuration of the initial state is shown in Fig. 1(a): a flux rope is suspended above a bundle of closed arcades, resulting in a X-type magnetic structure below the rope. The green curves in Fig. 1(a) mark the boundary of the flux rope and that of the arcades below the rope. The magnetic flux of these arcades per unit length along the z−direction is Φ a = 0.069 × 10 10 Mx cm −1 . The background field is a partially open bipolar field, with a negative and a positive surface magnetic charge located at the lower photosphere within −25 Mm < x < −15 Mm and 15 Mm < x < 25 Mm, respectively. B z is zero in the background field. Similar as that in Paper I, the emerging fibril in the scenario of flux feeding is represented by a small flux rope, and the pre-existing large flux rope is called "major rope" for simplicity in the rest of this paper. The emergence of the small rope is achieved by similar simulating procedures as those in Paper I: the small rope emerges from the central region right below the major rope at a constant speed; the emergence begins at t = 0 and ends at t = τ E = 60τ A , where τ A = 17.4 s is the characteristic Alfvén transit time. Thus the emerged part of the small rope at the base is located within − , and a is the radius of the small rope. Based on this, by adjusting ψ, B z , the velocities v x,y,z , the temperature T , and the density ρ at the base (y = 0, −x E x x E ), the emergence of the small flux rope is achieved by: Here v E = 2a/τ E is the emerging velocity; ψ i is the magnetic flux function of the initial state, and it is determined by numerical means (please see Eq. A.10 in Appendix A). It is noteworthy that ψ at the lower base is fixed at ψ i except during the emergence of the small rope, so that the normal component of the photospheric magnetic field is unchanged, implying that the lower base corresponds to the photosphere (Lin et al. 2003). The parameter C E determines the magnetic field strength of the small rope; its dimensionless values given in this paper are in the unit of 0.373 × 10 1 0 Mx cm −1 . In both the small and the major ropes, B z is positive and the component of the magnetic field in x-y plane is counterclockwise. Anomalous resistivity is used in our simulations, so as to restrict magnetic reconnection within the regions of current sheets: Here η m = 9.95 × 10 −2 , L 0 = 10 7 m, and v 0 = 128.57 km s −1 ; µ 0 is the vacuum magnetic permeability; the critical current density is j c = 4.45 × 10 −4 A m −2 .

Flux feeding process
With the procedures introduced above, we are able to simulate flux feeding processes in HFT configurations. The simulation results of a typical flux feeding process with C E = 2.23 are illustrated in Fig. 1(b)-1(i); different colors are used to depict the magnetic fields lines corresponding to different values of the flux function ψ, so as to clearly demonstrate the interaction between the major rope and emerging small rope in detail.
At the beginning of the flux feeding process, the small rope starts to emerge from the lower boundary, as marked by the brown curve in Fig. 1(b). The arcades below the major rope in the initial state are then pushed upward by the emerging small rope, so that the top boundary of these arcades reaches the lower boundary of the major rope (as shown by the green curves in Fig. 1(c)), and then these arcades reconnect with the magnetic fields of the major rope ( Fig. 1(d)). A comparison between Fig. 1(d) and Fig. 1(b) shows that the outer section of the major rope is peeled off by the reconnection. After all the arcades have reconnected with the major rope, the small rope itself also reaches the major rope. During the first half (0 ∼ τ E /2 = 30τ A ) of the flux feeding process, what emerges from the lower base is the top half of the small rope, in which the direction of B x is opposite to that in the bottom half of the major rope. Thus a horizontal current sheet forms at the interface between the two ropes. This current sheet can be clearly recognized in Fig. 2(a), in which the blue color depicts the distribution of the current density. As a result of the reconnection within this current sheet, the two ropes gradually merge together, as shown by, for example, the magnetic field lines plotted by the red and the pink colors in Fig. 1(e)-1(f). The magnetic cancelation occurring in this current sheet should decrease the local magnetic pressure below the major rope, so that the major rope gradually descends during the early period of the flux feeding process (see Fig. 1(f) and 1(g), also see the height profile of the major rope axis plotted in Fig. 2(i)). During the second half (30τ A ∼ 60τ A ) of the flux feeding process, however, it is the bottom half of the small rope that emerges from the lower base; the direction of the magnetic field lines in the small rope is now the same as that in the bottom half of the major rope; there should be no reconnection any more. Therefore, after t = 30τ A , magnetic flux is injected from the lower base, so that the local magnetic pressure accumulates below the major rope, which causes the major rope to gradually rise, as shown by Fig. 1(g)-1(h).
Eventually, the flux feeding process ends at t = 60τ A . The configuration of the resultant flux rope after flux feeding is illustrated in Fig. 1(i): the rope sticks to the photosphere, without any arcade below the rope.

Eruption caused by flux feeding
Further evolution of the resultant major rope after flux feeding indicates that the major rope is caused to erupt by this flux feeding process with C E = 2.23. As demonstrated by Fig. 2, the major rope keeps rising after flux feeding. Due to the line-tying effect from the photosphere, the rope does not separate from the lower base immediately, but its lower boundary sticks to the photosphere for a certain period ( Fig. 2(b)). As the rope rises, the lower section of the flux rope and the ambient background field are stretched. Eventually, as shown in Fig. 2(c)-2(d), the rope is fully detached from the lower base, and magnetic reconnection occurs within the vertical current sheet formed below the rope. The reconnection results in closed arcades below the rope, which can be clearly recongized in the lower section of Fig. 2(d). Moreover, as shown in Fig. 2(e), a plasmoid appears within the current sheet below the rope, and then rises and merges with the major rope ( Fig. 2(f)), which is similar as the observations and simulations in Gou et al. (2019). This process is followed by another similar one, as shown in Fig. 2(g)-2(h). These plasmoids should result from the tearing mode instability within the current sheet, which have been investigated in detail in many previous studies (e.g., Bárta et al. 2008;Shen et al. 2011;Ni et al. 2012;Huang et al. 2017).
The height of the major rope axis versus time is plotted in Fig. 2(i). As explicated in Sect. 3.1, due to the magnetic cancellation during the early period of the flux feeding process, the major flux rope gradually descends, which is accompanied by the contraction of the background arcades above the major rope (see Fig. 1(f) and 1(g)). This kind of downward movement before the onset of the eruption is often considered as a signature of "coronal implosion" (Hudson 2000), which has been reported in many previous observational studies (e.g., Ji et al. 2004;Veronig et al. A&A proofs: manuscript no. ams 2006; Joshi et al. 2009;Liu et al. 2009;Li et al. 2017). For example, Liu et al. (2009) observed the large-scale contraction of the overlying coronal loops before the onset of the flare, which was attributed to the prolonged preheating phase dominated by coronal thermal emissions. Our simulation results suggest another scenario of coronal implosions: internal magnetic cancellation during the early period of the flux feeding process in HFT configurations decreases the local magnetic pressure below the major flux rope, resulting in the downward movement of the rope before the onset of its eruption. For comparison, the height profile of the flux rope eruption caused by flux feeding in BPS configurations is monotonic, without any downward movement before the eruption (see Paper I).
It has been suggested in many previous studies that solar eruptions not only could be initiated by the gradual variations of photospheric magnetic fields, but also might in turn cause rapid magnetic field changes in the photosphere (Wang et al. 1994;Liu et al. 2005;Wang & Liu 2015;Toriumi & Wang 2019). The issue about this kind of photospheric response to the coronal eruptions is graphically called "tail wags the dog" problem (Aulanier 2016). This phenomenon has been observed in many previous studies (e.g., Wang et al. 2012;Liu et al. 2012;Petrie 2012;Sun et al. 2017;Castellanos Durán et al. 2018). For example, by measuring the various magnetic parameters of a compact region at the central polarity inversion line (PIL) of an active region, Sun et al. (2017) found that the horizontal component of the photospheric magnetic field in this compact region obviously increase after the flare. In our simulation results, Fig. 3 shows the photospheric distribution of the horizontal magnetic field, B x , at the central region right below the major rope (we note that B z is zero in the background field). In the initial state, as shown in Fig. 3(e), B x at the photosphere is negative, consistent with the arcades below the rope in the initial state. During the eruptive process of the resultant major rope after flux feeding, closed arcades forms right below the rising rope, which also results in the negative photospheric B x , as shown in Fig. 3(f). Moreover, as the rope rises, more and more field lines are closed by reconnection and pile up above the PIL. As a result of the reconnection-driven contraction of the flare loops, the horizontal photospheric magnetic field B x increases accordingly (see Fig. 3(f)-3(h)), which is consistent with the simulation results in Barczynski et al. (2019).
Comparing the photospheric B x at t = 200τ A with that in the initial state, the average strength of the photospheric B x after the eruption increases from 6.27 G to 8.50 G. Moreover, as shown in Fig. 3, the magnetic field strength is of the order of 10 G, indicating that coronal eruptions originating in quiescent regions could also result in the phenomenon of "tail wags the dog", that is coronal eruptions originating in quiescent regions could also cause variations of the photospheric magnetic fields.

Initiation of Eruptions
In order to investigate how the eruption is initiated, we first investigate the influence of flux feeding on the major flux rope. The magnetic properties of a flux rope could be characterized by its axial magnetic flux, Φ z , and its poloidal magnetic flux per unit length along the z−direction, Φ p . For the initial major rope shown in Fig. 1(a), the axial flux Φ z0 = 4.366 × 10 19 Mx, and the poloidal flux Φ p0 = 1.186 × 10 10 Mx cm −1 . Here we select the flux rope at t = 80τ A (Fig. 1(i)) as the resultant rope after flux feeding, whose axial flux Φ z = 7.115 × 10 19 Mx, and poloidal flux Φ p = 1.116 × 10 10 Mx cm −1 . The increased axial flux is injected by the flux feeding process, which is similar as the BPS cases investigated in Paper I. Different from the result in BPS cases, however, the poloidal flux is reduced by 0.070 × 10 10 Mx cm −1 after flux feeding.
Since the distribution of the coronal and the photospheric magnetic fields play a dominant role in triggering coronal flux rope eruptions (Yeates 2014;Yang et al. 2018;Thalmann et al. 2019;Xing et al. 2020), the influence of flux feeding processes on the major rope system should be sensitive to the scale of the magnetic field strength in the small rope. The magnetic parameters of the resultant flux ropes after the flux feeding processes with different C E are tabulated in Table 1. Larger C E implies stronger magnetic field strength in the small emerging rope, so that more axial flux is injected into the major rope. The deduced poloidal fluxes in different cases, however, are very close to each other: ∆Φ p = 0.069 ∼ 0.070 × 10 10 Mx cm −1 , which is almost the same as the magnetic flux of the arcades below the rope in the initial state, Φ a . This indicates that the reduced poloidal flux should be caused by the interaction between the major rope and the arcades below the rope with the HFT configuration, which  has been introduced in Sect. 3.1: as pushed by the rising small flux rope, these arcades reconnect with the major rope, so that the outer section of the major rope is peeled off. In the cases with differrent C E , the further evolution of the resultant rope also varies, as tabulated in the last column in Table 1. The eruptive and the non-eruptive cases are well separated: the resultant rope erupts if C E is no smaller than 2.23, and the axial fluxes of the resultant ropes in all the eruptive cases are larger than those in the non-eruptive cases. This indicates that there should exist a critical axial magnetic flux of the resultant rope, and the onset of its eruption should be initiated by the upward catastrophe associated with this critical axial flux (e.g., Su et al. 2011;Zhang et al. 2016). The value of the critical axial flux is of the order of 7.1 × 10 19 Mx. The simulation results of the critical non-eruptive case with C E = 2.22 are illustrated in Fig. 4: the resultant rope does not keep rising after flux feeding ( Fig. 4(b)-4(c)), but falls back to the photosphere (Fig. 4(d)-4(e)), resulting in an equilibrium state in which the rope sticks to the photosphere (Fig. 4(f)). Therefore, if the injected axial flux is not large enough so that the critical axial flux of the major rope is not reached, not only is the flux feeding process unable to cause the major rope to erupt, but also the major rope system is transformed from an HFT configuration to a stable BPS configuration. This result might also suggest a possible scenario for confined flares in HFT configurations.
As demonstrated in Fig. 2, obvious magnetic reconnection occurs in the vertical current sheet below the rising major flux rope. Here we simulate a specific case to distinguish between the role that reconnection plays and that upward catastrophe plays in initiating the flux rope eruption: C E is also 2.23, the same as the eruptive case shown in Fig. 2, but the magnetic reconnection in the current sheet below the rope is prohibited. To achieve this, similar simulating procedures as, for example, Zhang et al. (2016), are used here: after t = 80τ A , first, the resistivity is adjusted to zero, and second, reassign the flux function ψ along the entire vertical current sheet (if it exists) below the major rope with the initial value ψ c = ψ i (x = 0, y = 0) at each time step, so that ψ along the current sheet is invariant. With these procedures, both numerical and physical magnetic reconnections below the major rope are prohibited. The simulation results for A&A proofs: manuscript no. ams  Simulation results of the critical eruptive case with C E = 2.23, in which the magnetic reconnection below the major rope is prohibited. The meanings of the symbols are the same as those in Fig. 2. The red dotted line in panel(g) is the corresponding eruptive profile with reconnection, which is the same as that in Fig. 2 this specific case are illustrated in Fig. 5. After flux feeding, the upward catastrophe occurs and the resultant rope rises for a while (Fig. 5(a)-5(c)), but then stops rising, and eventually, the rope does not erupt but levitates in the corona (Fig. 5(d)-5(f)). Without magnetic reconnection, the closed background arcades above the major rope are always rooted at the photosphere. These closed arcades constrain the major rope, so that prevent the rope from further rising. For comparison, in the case with reconnection demonstrated in Fig. 2, background arcades above the rope reconnect at the vertical current sheet below the rope, so that the "tether" constraining the rising major rope is removed. As suggested by Green et al. (2018), different kinds of eruptive mechanisms can be classified into "trigger" and "driver" according to their predominant role in the eruption process. In our simulation results, the flux feeding process causes the upward catastrophe, which triggers the flux rope eruption. Magnetic reconnection as well as the further evolution of the upward catastrophe act as the drivers of the flux rope eruption.

Discussion and conclusion
In this paper, we investigate the influence of flux feeding on coronal magnetic flux ropes with an HFT configuration, especially how flux feeding causes the flux rope to erupt. During the flux feeding process, the small emerging flux rope as well as the arcades below the major rope with the HFT configuration reconnect with the pre-existing major flux rope, so that axial magnetic flux is injected into the major rope, whereas the poloidal magnetic flux of the major rope is reduced. Flux feeding is able to cause the major rope to erupt, provided that the amount of the axial flux inject by flux feeding is large enough so that the total axial flux of the major rope reaches its critical value, which is similar as that in the BPS cases investigated in Paper I. During the early period of the flux feeding process, the major rope gradually descends due to the magnetic cancellation, demonstrating a downward movement of the flux rope right before its eruption, which suggests an alternative theoretical scenario for coronal implosions. Moreover, our simulation reproduce the photospheric magnetic response to the eruption of the coronal flux rope: the horizontal component of the photospheric magnetic field increases around the PIL after the eruption, as is often observed in the wake of major solar eruptions.
Although the onset of the eruption in both the BPS and the HFT cases is dominated by the upward catastrophe associated with the critical axial flux, the influence of flux feeding on the flux ropes with an HFT configuration is different from those with a BPS configuration. With BPS, the major flux rope sticks to the photosphere, so that the small rope directly interacts with the major rope from the very beginning of the flux feeding process, and the major rope only reconnects with the small emerging rope during the flux feeding process; the non-eruptive resultant rope after flux feeding is still with a BPS configuration. For HFT cases, there exisits additional interaction between the major rope and the arcades below it in the initial state. Consequently, not only the poloidal flux of the major rope is reduced due to the reconnection with these arcades, but the HFT configuration collapses after flux feeding. The flux feeding in HFT configurations results in a dramatic change of topology, with the flux rope either erupting or falling back to a BPS configuration. Moreover, previous theoretical studies found that a flux rope could erupt if either its axial or poloidal magnetic flux exceeds the corresponding critical value. For BPS cases, each episode of flux feeding injects axial flux into the flux rope while maintaining its poloidal flux, therefore pushing the rope system one step closer toward the eruption, no matter how small the amount of the injected axial flux is. For HFT cases, however, if the injected axial flux is too little, the reduction in poloidal flux could push the rope system even farther away from the onset of the eruption, and, what's more, the rope system may also evolve from a metastable HFT configuration to a stable BPS configuration after flux feeding. Therefore, the flux feeding process in HFT configurations is not always in favor of the eruption of the major rope; it could even stabilize the rope system.
As introduced in Sect. 4, flux feeding cases with different C E are simulated in this paper. The emerging velocity of the small rope, however, is unchanged in different cases. Since the emerging velocity might also have some influence on the evolution of the resultant rope, we will try to simulate flux feeding cases with various emerging velocities in our future work, so as to investigate the role that the kinetic properties of flux feeding processes plays in initiating coronal flux rope eruptions.

Appendix A: Basic equations and initial preparations
Combining the 2.5-Dimensional form of the magnetic field (Eq. 1), we could write the MHD equations in the non-dimensional form as: Here ρ and T denote the density and the temperature; v x , v y , v z represent the x, y, and z−component of the velocity, respectively; γ = 5/3 is the polytropic index in our simulation; g is the normalized gravity; η is the resistivity. Here β 0 = 2µ 0 ρ 0 RT 0 L 2 0 /ψ 2 0 = 0.1 is the characteristic ratio of the gas pressure to the magnetic pressure, where ρ 0 = 3.34 × 10 −13 kg m −3 , T 0 = 10 6 K, L 0 = 10 7 m, and ψ 0 = 3.73 × 10 3 Wb m −1 are the characteristic values of density, temperature, length and magnetic flux function, respectively, which are also the calculating units in the simulation. For the other quantities, the corresponding characteristic values are v 0 = 128.57 km s −1 , t 0 = 77.8 s, B 0 = 3.37 × 10 −4 T, g 0 = 1.65 × 10 3 m s −2 . The numerical domain in our simulation is 0 < x < 200 Mm, 0 < y < 300 Mm, which is discretized into 400×600 uniform meshes. At the left side of the domain (x = 0), symmetric boundary condition is used. The radiation and the heat conduction in the energy equation are neglected.
In order to simulate flux feeding processes in HFT cases, we must first construct a typical coronal flux rope system with an HFT configuration. Here we select a partially open bipolar field as the background field, in which a negative and a positive surface magnetic charge located at the photosphere within −b < x < −a and a < x < b, respectively. The background field could then be obtained by the complex variable method (e.g., Hu et al. 1995;Zhang et al. 2017). The background magnetic field can be cast in the complex variable form where ω = x + iy, and There is a neutral current sheet located at (x = 0, y ≥ y N ). The magnetic flux function could then be calculated by: ψ(x, y) = Im f (ω)dω , (A.9) and the flux function at the lower base is where ψ c = πψ 0 ; the flux function at the neutral point y = y N can be calculated as: Letting B z equals 0 in the background field, the configuration of the background field could then be obtained, as shown in Fig. A.1(a). The initial corona is isothermal and static: T c ≡ T (0, x, y) = 1 × 10 6 K, ρ c ≡ ρ(0, x, y) = ρ 0 e −gy . (A.12) Except during the emergence of the small rope, the lower boundary is fixed: the flux function ψ is fixed at ψ i given by Eq. A.10; B z and the velocity is always 0; the density and the temperature are fixed at their initial values. Increment equivalent extrapolation is used at the right and top boundaries. With the initial condition and the background configuration, equations (A.1) to (A.5) are solved by the multi-step implicit scheme. First, we let a flux rope emerge from the lower base of the background field, as shown in Fig. A.1(b). By multiplying B z within the rope by a factor larger than 1, we increase the axial magnetic flux of the flux rope. The rope then rises, and magnetic reconnection occurs within the veritical current sheet formed below the rope, resulting in closed arcades below the rope. After that, adjust the axial flux of the rope again, so that the rope stops rising. Eventually, let the rope system relax to a equilibrium state with an HFT configuration, as shown in Fig. A.1(c). This is the initial state of our simulation. It is noteworthy that the radius of the flux rope in our simulation is finite, so that the thin-rope approximation is not satisfied. Under this circumstance, the initial state could only be obtained by numerical procedures.   (b) is an interim state to construct the initial state. Panel (c) is the obtained flux rope system with an HFT configuration, which is the initial state of our simulation.