{"id":93219,"date":"2026-10-15T21:59:25","date_gmt":"2026-10-15T19:59:25","guid":{"rendered":"https:\/\/fusion.bsc.es\/?p=93219"},"modified":"2026-10-15T21:59:25","modified_gmt":"2026-10-15T19:59:25","slug":"my-masters-thesis-at-bsc-simulating-electromagnetic-wave-propagation-for-fusion","status":"publish","type":"post","link":"https:\/\/fusion.bsc.es\/index.php\/2026\/10\/15\/my-masters-thesis-at-bsc-simulating-electromagnetic-wave-propagation-for-fusion\/","title":{"rendered":"My Master\u2019s Thesis at BSC: Simulating Electromagnetic Wave Propagation for Fusion"},"content":{"rendered":"<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-93229\" src=\"https:\/\/fusion.bsc.es\/wp-content\/uploads\/1774980330751-e1789732567885-400x226.jpg\" alt=\"\" width=\"750\" height=\"423\" srcset=\"https:\/\/fusion.bsc.es\/wp-content\/uploads\/1774980330751-e1789732567885-400x226.jpg 400w, https:\/\/fusion.bsc.es\/wp-content\/uploads\/1774980330751-e1789732567885.jpg 800w\" sizes=\"auto, (max-width: 750px) 100vw, 750px\" \/><\/p>\n<p style=\"text-align: justify;\">Hello everyone! My name is <strong>Adam Canarutto<\/strong> and over the past six months, <strong>I had the opportunity to carry out my Master&#8217;s thesis within the Fusion Group at the <a href=\"http:\/\/bsc.es\">Barcelona Supercomputing Center (BSC)<\/a><\/strong>. My work focused on a problem at the intersection of plasma physics, numerical modelling and high-performance computing: the simulation of electromagnetic wave propagation for Ion Cyclotron Resonance Heating (ICRH).<\/p>\n<p style=\"text-align: justify;\">Throughout my project, I had the privilege of being guided by <strong>Alejandro Soba<\/strong>, <strong>Ezequiel Goldberg<\/strong>, and <strong>Hern\u00e1n Domingo<\/strong>, and I am incredibly grateful for their support. From the very beginning, they were always there to guide me, answer my questions, and help me overcome any difficulties along the way, while creating a friendly and enjoyable working environment.<\/p>\n<p><!--more--><\/p>\n<p class=\"isSelectedEnd\" style=\"text-align: justify;\">Understanding wave propagation in a magnetised plasma is a challenging numerical problem. The plasma response is anisotropic and depends on parameters such as the magnetic field, density, particle species and wave frequency. At the same time, realistic fusion simulations require large computational resources and therefore efficient parallel numerical methods.<\/p>\n<p class=\"isSelectedEnd\" style=\"text-align: justify;\">My work was carried out using <a href=\"https:\/\/www.bsc.es\/research-development\/research-areas\/engineering-simulations\/alya-high-performance-computational\"><strong>Alya<\/strong><\/a>, BSC&#8217;s multiphysics simulation code designed for large-scale parallel computing, and specifically its <strong>EMWAVE<\/strong> module for electromagnetic-wave simulations based on the finite element method. The plasma description used in the project relies on the cold-plasma approximation, which provides a useful foundation for studying wave propagation before moving towards more complete plasma models.<\/p>\n<p class=\"isSelectedEnd\" style=\"text-align: justify;\">A significant part of the thesis involved extending and testing EMWAVE&#8217;s numerical capabilities. I implemented and verified <strong>quadratic triangular and quadrilateral finite elements<\/strong>, introduced <strong>electrical conductivity<\/strong> into the electromagnetic formulation, and extended the available outputs to include the magnetic flux density <strong>B<\/strong>, magnetic field <strong>H<\/strong>, and current density <strong>J<\/strong>, in addition to the electric field. This required working directly with Alya&#8217;s Fortran code, understanding its finite-element formulation, debugging the implementation and designing numerical tests to verify the new features.<\/p>\n<p class=\"isSelectedEnd\" style=\"text-align: justify;\"><strong>Verification was a central part of the project.<\/strong> Numerical results were first tested using controlled electromagnetic problems and then compared with results from <strong>ERMES<\/strong>, an independent electromagnetic solver. The agreement between the two codes for the electromagnetic fields and derived quantities provided confidence in the implementation and highlighted the importance of independent verification when developing scientific software.<\/p>\n<p class=\"isSelectedEnd\" style=\"text-align: justify;\">I also moved from local tests to large-scale simulations on <a href=\"https:\/\/www.bsc.es\/marenostrum\/marenostrum\"><strong>MareNostrum 5<\/strong><\/a>, BSC&#8217;s supercomputer. Parallel scalability studies allowed me to investigate how the solver behaves as the computational workload is distributed across increasing numbers of processors. This provided practical experience with parallel efficiency, workload distribution, memory requirements and the challenges of running large HPC simulations. It also demonstrated that, for realistic fusion applications, reproducing the correct physics is only part of the problem: numerical methods must also exploit modern supercomputing architectures efficiently.<\/p>\n<figure id=\"attachment_93224\" aria-describedby=\"caption-attachment-93224\" style=\"width: 740px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-93224\" src=\"https:\/\/fusion.bsc.es\/wp-content\/uploads\/scaling_test_tol1e6-1024x634.png\" alt=\"\" width=\"750\" height=\"465\" srcset=\"https:\/\/fusion.bsc.es\/wp-content\/uploads\/scaling_test_tol1e6-1024x634.png 1024w, https:\/\/fusion.bsc.es\/wp-content\/uploads\/scaling_test_tol1e6-400x248.png 400w, https:\/\/fusion.bsc.es\/wp-content\/uploads\/scaling_test_tol1e6-1536x952.png 1536w, https:\/\/fusion.bsc.es\/wp-content\/uploads\/scaling_test_tol1e6.png 1590w\" sizes=\"auto, (max-width: 750px) 100vw, 750px\" \/><figcaption id=\"caption-attachment-93224\" class=\"wp-caption-text\">Scalability analysis performed on Marenostrum 5. The graphs shows the runtime, speedup and parallel efficiency of the execution. Nodes are just a division of MPI workers, with each node having a maximum of 112 processes. These results clearly indicate room for improvement in the EMWAVE module.<\/figcaption><\/figure>\n<p class=\"isSelectedEnd\" style=\"text-align: justify;\">In the final stage of the thesis, I explored coupling Alya with <a href=\"https:\/\/snl-dakota.github.io\/docs\/6.19.0\/users\/usingdakota\/introduction\/aboutdakota.html\"><strong>DAKOTA<\/strong><\/a>, a toolkit for parameter studies, optimisation and uncertainty quantification. I used this coupling to automate frequency-sensitivity studies and examine how the electromagnetic response changes with excitation frequency. These studies also showed the importance of ensuring sufficient numerical resolution before drawing physical conclusions from parameter scans. A more detailed explanations of the results can be found in my <a href=\"https:\/\/fusion.bsc.es\/wp-content\/uploads\/High_Performance_Finite_Element_Simulation_of_Electromagnetic_Waves_for_Ion_Cyclotron_Resonance_Heating_compressed.pdf\">Master&#8217;s Thesis<\/a>.<\/p>\n<p>The thesis concluded with its presentation in <strong>September 2026<\/strong>, but the work forms part of a longer-term development of EMWAVE. Future directions include more realistic geometries and, eventually, moving beyond the cold-plasma approximation towards <strong>hot-plasma models<\/strong> capable of describing the wave-particle interactions responsible for plasma heating.<\/p>\n<p class=\"isSelectedEnd\" style=\"text-align: justify;\">Beyond the technical developments, the internship gave me the opportunity to experience scientific research within a large computational organisation. Coming from a background in plasma and fusion physics, I had to connect physical concepts with finite-element methods, scientific software development and HPC. Working on a large research code also taught me the importance of careful testing, integration and reproducibility in a parallel computing environment.<\/p>\n<p style=\"text-align: justify;\">For me, the project brought together three areas I would like to continue exploring: <strong>fusion plasma physics, numerical modelling and high-performance scientific computing<\/strong>. It was also a memorable opportunity to spend six months in Barcelona, surrounded by an international research community and supported by colleagues who made the experience valuable both scientifically and personally. I am very grateful to my supervisors and to the BSC Fusion Group for the opportunity and for everything I learned during these six months. Finally, I would also like to acknowledge the financial support provided by FuseNet and the Erasmus+ programme, which contributed to making this research stay at BSC possible.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Hello everyone! My name is Adam Canarutto and over the past six months, I had the opportunity to carry out my Master&#8217;s thesis within the Fusion Group at the Barcelona Supercomputing Center (BSC). My work focused on a problem at the intersection of plasma physics, numerical modelling and high-performance computing: the simulation of electromagnetic wave &#8230; <a title=\"My Master\u2019s Thesis at BSC: Simulating Electromagnetic Wave Propagation for Fusion\" class=\"read-more\" href=\"https:\/\/fusion.bsc.es\/index.php\/2026\/10\/15\/my-masters-thesis-at-bsc-simulating-electromagnetic-wave-propagation-for-fusion\/\" aria-label=\"Read more about My Master\u2019s Thesis at BSC: Simulating Electromagnetic Wave Propagation for Fusion\">Read more<\/a><\/p>\n","protected":false},"author":47,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":"","_members_access_role":[],"_members_access_error":""},"categories":[9,1],"tags":[],"class_list":["post-93219","post","type-post","status-publish","format-standard","hentry","category-group","category-news"],"_links":{"self":[{"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/posts\/93219","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/users\/47"}],"replies":[{"embeddable":true,"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/comments?post=93219"}],"version-history":[{"count":14,"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/posts\/93219\/revisions"}],"predecessor-version":[{"id":93274,"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/posts\/93219\/revisions\/93274"}],"wp:attachment":[{"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/media?parent=93219"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/categories?post=93219"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/fusion.bsc.es\/index.php\/wp-json\/wp\/v2\/tags?post=93219"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}