In the project kick-off meeting the proposed research plan was reviewed and we agreed that the work package WP3 (theoretical analysis) will have a wider scientific impact if the theoretical developments were done for matrices of arbitrary dimensions. The tasks in the work packages WP3 and WP4 were modified to accommodate this generality.
I trained in the theory of matrix functions, backward-error analysis, Padé approximations and in the development of numerically stable algorithms to compute matrix functions. A career development plan was put together to improve my career prospects in academia. Teaching experience, academic service and leadership are essential criteria in any academic job application. As my CV was weak in these areas, I did the following professional development activities.
1. Trained in university research and teaching (New Academic Programme) and public engagement (Royal Society).
2. Taught one-third of a Level 4 course at the School of Mathematics.
3. Assumed an administrative role (member of School's Athena SWAN committee).
4. Assumed an organisation role (Numerical Analysis seminars).
The first two work packages of the project are related to project management (WP1) and quality control (WP2). So there are no scientific results to be reported here. The next two work packages related to theoretical analysis (WP3) and software implementation (WP4) constitute the original research work of the project.
In the work package WP3, the properties of the wave-kernel functions and their inverses were established, the theory for the backward-error analysis was developed and the error estimates of approximations to compute the wave-kernel matrix functions were established. Based on these results, in work package WP4, a cost-optimal algorithm was developed. Additionally, MATLAB functions to compute and test the wave-kernel matrix functions were developed and this software is freely available in the GitHub repository:
https://github.com/nadukandi/wkm(se abrirá en una nueva ventana)A manuscript containing the theoretical results and the algorithms for computing the wave-kernel matrix functions was written and submitted to the SIAM Journal on Scientific Computing.
EPrint:
http://eprints.maths.manchester.ac.uk/2651/(se abrirá en una nueva ventana)Additionally, an algorithm was developed for the piecewise integration of particle trajectories on triangular meshes distributed over several computer nodes. High-performance software implementation of this algorithm has been done in C++ and using MPI/OpenMP directives. The software is hosted in the Bitbucket repository:
https://bitbucket.org/nadukandi/slpfem(se abrirá en una nueva ventana)The work package WP5 (simulation campaign) has been partially accomplished. From the perspective of the FastFlowSim project, a flow solver has two components. The first is the advective transport of flow data, e.g. position, velocity and passive intensive scalars like density of the flow. It is to this component of the flow solver that the theoretical, algorithmic and software developments of the FastFlowSim project contribute. The second component corresponds to the dispersion of the aforesaid flow data due to viscous processes and imposition of the incompressibility constraint. The strategy adopted here is to delegate the tasks associated to this component to open source numerical computation software, e.g. MUMPS (linear solvers) and BDDCML (domain decomposition preconditioners). The compilation of these open source software has been done and the interface to integrate the software developed in the FastFlowSim project has been prepared. The combined testing of the integrated software and its performance analysis on flow problems is pending. This delay is a result of the time allocated to the successful generalization of the theoretical developments in WP3 and in the training required to write a manuscript with appropriate mathematical rigour.
The last work package is on dissemination and communication (WP6). The theoretical results were presented at the SIAM UKIE annual meetings and the 27th Biennial Numerical Analysis Conference. I attended a workshop on “Computational Aspects of Time Dependent Electromagnetic Wave Problems in Complex Materials”, Brown University (USA). I presented my research findings and networked with world leading scientists working on wave-propagation problems, e.g. Prof. Thomas Hagstrom, Dr. Daniel Appelo and Dr. Vladimir Druskin. For the same purpose, I also visited Prof. Gilbert Strang (Mathematics, MIT) and Prof. Ram-Mohan Ramdas (Physics, WPI).
I successfully led research stations and Q&A sessions at two science festivals (Science Uncovered 2017 and Science Spectacular 2017) held in Manchester. Additionally, I presented my work at the STEM for BRITAIN 2018 poster competition held at the House of Commons, UK Parliament.
