Numerical simulations of passive scalar advection in MHD turbulence using OpenFOAM

One of the biggest challenges has ever faced in human history is to be self-sufficient in terms of energy. In the context of the global warming, it is important to find sustainable solutions for our growing need on energy without putting at risk the lives of our generation and the ones to come. Nuclear fusion is one of the possible solutions of a potential clean and unlimited energy in the near future. In order to achieve the nuclear fusion on the Earth, the plasma has to be controlled in a magnetically confined region. One of the simple but useful model to understand the plasma dynamics and its behavior is magnetohydrodynamics (MHD). Computational fluids dynamic (CFD) is used to solve the complex MHD set of equations which governs the plasma. This requires high computational resources, which can be time and economically expensive. In order to tackle the problem, free-licensed software, OpenFOAM is applied. OpenFOAM has a module which solves MHD, mhdFOAM. The aim of the project is to implement the passive scalar model which is a simple description for temperature advection in the present module of OpenFOAM. In the present work, the implementation of the passive scalar model has been achieved. The simulations of the passive scalar advection in the MHD have been performed using Marconi machine which is one of powerful European supercomputers. The simulation results show a correlation between the passive scalar and the velocity field which is reasonable to consider the implementation is correctly achieved