The objective of this project is to incorporate magnetic effects into the existing wave propagation code Hawen, which currently models wave oscillations in a variety of media. This extension will require implementing the linear induction equation and coupling it with the linear equation of motion already available in the code.
The project will begin with a familiarization phase, during which the successful candidate will gain an understanding of the team's computational framework and numerical methods. They will then develop and implement the appropriate numerical strategy, validate the implementation through representative benchmark problems, and document the new developments. In the final stage of the project, the code will be used to investigate the influence of magnetic fields on solar oscillations by comparing numerical predictions with observational data. The provisional schedule is as follows:
1. Familiarization with numerical framework in Hawen (structure of the benchmarks, code parallelism, I/O), the discretization method used in the code, and the existing solar oscillation equations currently available.
2. Implement an initial strategy which consists in a zero-th order term to add in the system of equations. Indeed, when ignoring the perturbation in magnetic field, the contribution of the Lorentz force in the equation of motion can be included as a zero-th order term, which avoid an increase in the number of unknowns. This step will help the successful candidate to familiarize with the code structure.
3. Implement the full coupling between the induction equation and the oscillation equation, and compare with the strategy developed in step 2.
4. Develop automated Python routines to post-process Hawen simulations into helioseismic observables (such as power spectrum, cross-covariance), enabling a quantitative investigation of the impact of considering the magnetic field. This step also includes a familiarization with observables and data processing performed in existing Python routines to process Dopplergrams to helioseismic observable.
