“…It also has the advantage to investigate the transient fast flow with designated injection interval by changing the boundary conditions, which remains a challenge for standalone MHD codes. Similar to our previous work (e.g., Wang et al., 2021; Yang et al., 2011), we set up the initial and boundary conditions using the empirical magnetic field (Tsyganenko, 1989) and plasma (Tsyganenko & Mukai, 2003) models driven by typical solar wind and geomagnetic conditions, that is, IMF B Z = –3.0 nT, solar wind number density N = 5.0 cm −3 , velocity V = 400 km/s, Kp = 1.0, cross‐polar‐cap potential drops = 50 kV. As most bursty bulk flows (BBFs) are found between the mid‐tail region and the near‐Earth region (−30 < X < −10 Re) and with the velocity greater than 300 km/s (e.g., Angelopoulos et al., 1994; Ohtani et al., 2004), we inject the bubble from the tail at −30 Re following 30 min of adiabatic convection.…”