Research Outputs

Now showing 1 - 2 of 2
  • Publication
    Finite element approximation of the spectrum of the curl operator in a multiply connected domain
    (Springer, 2019)
    Alonso Rodríguez, Ana María
    ;
    ;
    Rodríguez, R.
    ;
    Valli, A.
    ;
    Venegas, P.
    In this paper we are concerned with two topics: the formulation and analysis of the eigenvalue problem for the curlcurl operator in a multiply connected domain and its numerical approximation by means of finite elements. We prove that the curlcurl operator is self-adjoint on suitable Hilbert spaces, all of them being contained in the space for which curlvv⋅nn=0curl⁡vv⋅nn=0 on the boundary. Additional constraints must be imposed when the physical domain is not topologically trivial: we show that a viable choice is the vanishing of the line integrals of vvvv on suitable homological cycles lying on the boundary. A saddle-point variational formulation is devised and analyzed, and a finite element numerical scheme is proposed. It is proved that eigenvalues and eigenfunctions are efficiently approximated and some numerical results are presented in order to assess the performance of the method.
  • Publication
    Assessment of two approximation methods for the inverse problem of electroencephalography
    (International journal of numerical analysis and modeling, 2016) ;
    Alonso-Rodríguez, A.
    ;
    Rodríguez, R.
    ;
    Valli, A.
    The goal of this paper is to compare two computational models for the inverse problem of electroencephalography: the localization of brain activity from measurements of the electric potential on the surface of the head. The source current is modeled as a dipole whose localization and polarization has to be determined. Two methods are considered for solving the corresponding forward problems: the so called subtraction approach and direct approach. The former is based on subtracting a fundamental solution, which has the same singular character of the actual solution, and solving computationally the resulting non-singular problem. Instead, the latter consists in solving directly the problem with singular data by means of an adaptive process based on an aposteriori error estimator, which allows creating meshes appropriately refined around the singularity. A set of experimental tests for both, the forward and the inverse problems, are reported. The main conclusion of these tests is that the direct approach combined with adaptivity is preferable when the localization of the dipole is close to an interface between brain tissues with different conductivities.