No Arabic abstract
The Coronal Multichannel Polarimeter (CoMP) routinely performs coronal polarimetric measurements using the Fe XIII 10747 $AA$ and 10798 $AA$ lines, which are sensitive to the coronal magnetic field. However, inverting such polarimetric measurements into magnetic field data is a difficult task because the corona is optically thin at these wavelengths and the observed signal is therefore the integrated emission of all the plasma along the line of sight. To overcome this difficulty, we take on a new approach that combines a parameterized 3D magnetic field model with forward modeling of the polarization signal. For that purpose, we develop a new, fast and efficient, optimization method for model-data fitting: the Radial-basis-functions Optimization Approximation Method (ROAM). Model-data fitting is achieved by optimizing a user-specified log-likelihood function that quantifies the differences between the observed polarization signal and its synthetic/predicted analogue. Speed and efficiency are obtained by combining sparse evaluation of the magnetic model with radial-basis-function (RBF) decomposition of the log-likelihood function. The RBF decomposition provides an analytical expression for the log-likelihood function that is used to inexpensively estimate the set of parameter values optimizing it. We test and validate ROAM on a synthetic test bed of a coronal magnetic flux rope and show that it performs well with a significantly sparse sample of the parameter space. We conclude that our optimization method is well-suited for fast and efficient model-data fitting and can be exploited for converting coronal polarimetric measurements, such as the ones provided by CoMP, into coronal magnetic field data.
We introduce and investigate matrix approximation by decomposition into a sum of radial basis function (RBF) components. An RBF component is a generalization of the outer product between a pair of vectors, where an RBF function replaces the scalar multiplication between individual vector elements. Even though the RBF functions are positive definite, the summation across components is not restricted to convex combinations and allows us to compute the decomposition for any real matrix that is not necessarily symmetric or positive definite. We formulate the problem of seeking such a decomposition as an optimization problem with a nonlinear and non-convex loss function. Several mode
We present an investigation of a coronal cavity observed above the western limb in the coronal red line Fe X 6374 {AA} using a telescope of Peking University and in the green line Fe XIV 5303 {AA} using a telescope of Yunnan Observatories, Chinese Academy of Sciences during the total solar eclipse on 2017 August 21. A series of magnetic field models are constructed based on the magnetograms taken by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO) one week before the eclipse. The model field lines are then compared with coronal structures seen in images taken by the Atmospheric Imaging Assembly on board SDO and in our coronal red line images. The best-fit model consists of a flux rope with a twist angle of 3.1$pi$, which is consistent with the most probable value of the total twist angle of interplanetary flux ropes observed at 1 AU. Linear polarization of the Fe XIII 10747 {AA} line calculated from this model shows a lagomorphic signature that is also observed by the Coronal Multichannel Polarimeter of the High Altitude Observatory. We also find a ring-shaped structure in the line-of-sight velocity of Fe XIII 10747 {AA}, which implies hot plasma flows along a helical magnetic field structure, in the cavity. These results suggest that the magnetic structure of the cavity is a highly twisted flux rope, which may erupt eventually. The temperature structure of the cavity has also been investigated using the intensity ratio of Fe XIII 10747 {AA} and Fe X 6374 {AA}.
Coronal waves are large-scale disturbances often driven by coronal mass ejections (CMEs). We investigate a spectacular wave event on 7 March 2012, which is associated with an X5.4 flare (SOL2012-03-07). By using a running center-median (RCM) filtering method for the detection of temporal variations in extreme ultraviolet (EUV) images, we enhance the EUV disturbance observed by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) and the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) onboard the PRoject for Onboard Autonomy 2 (PROBA2). In coronagraph images, a halo front is observed to be the upper counterpart of the EUV disturbance. Based on the EUV and coronagraph images observed from three different perspectives, we have made three-dimensional (3D) reconstructions of the wave surfaces using a new mask-fitting method. The reconstructions are compared with those obtained from forward-fitting methods. We show that the mask fitting method can reflect the inhomogeneous coronal medium by capturing the concave shape of the shock wave front. Subsequently, we trace the developing concave structure and derive the deprojected wave kinematics. The speed of the 3D-wave nose increases from a low value below a few hundred $mathrm{km,s^{-1}}$ to a maximum value of about 3800 $mathrm{km,s^{-1}}$, and then slowly decreases afterwards. The concave structure starts to decelerate earlier and has significantly lower speeds than those of the wave nose. We also find that the 3D-wave in the extended corona has a much higher speed than the speed of EUV disturbances across the solar disk.
We propose and test the first Reduced Radial Basis Function Method (R$^2$BFM) for solving parametric partial differential equations on irregular domains. The two major ingredients are a stable Radial Basis Function (RBF) solver that has an optimized set of centers chosen through a reduced-basis-type greedy algorithm, and a collocation-based model reduction approach that systematically generates a reduced-order approximation whose dimension is orders of magnitude smaller than the total number of RBF centers. The resulting algorithm is efficient and accurate as demonstrated through two- and three-dimensional test problems.
Recently, collocation based radial basis function (RBF) partition of unity methods (PUM) for solving partial differential equations have been formulated and investigated numerically and theoretically. When combined with stable evaluation methods such as the RBF-QR method, high order convergence rates can be achieved and sustained under refinement. However, some numerical issues remain. The method is sensitive to the node layout, and condition numbers increase with the refinement level. Here, we propose a modified formulation based on least squares approximation. We show that the sensitivity to node layout is removed and that conditioning can be controlled through oversampling. We derive theoretical error estimates both for the collocation and least squares RBF-PUM. Numerical experiments are performed for the Poisson equation in two and three space dimensions for regular and irregular geometries. The convergence experiments confirm the theoretical estimates, and the least squares formulation is shown to be 5-10 times faster than the collocation formulation for the same accuracy.