codes

If you would like to use these codes, please contact me via anna.guseva at obspm.fr

Magnetohydrodynamic Taylor-Couette flow code

MPI-parallelized code solving coupled Navier-Stokes and induction equations for a flow between two rotating cylinders, developed in collaboration with A.P. Willis (U. Sheffield, UK). The code has spectral discretization in the axial and azimuthal directions, and the radial coordinate is discretized using finite differences. The nonlinear terms are evaluated in the physical space and are de-aliased using the 3/2 rule; more details of the numerical method can be found in (Guseva et al., 2015). The code showed excellent scalability in high-performance computing centers (Jülich Supercomputing Centre), and allowed to model low-resistivity fluids, unachievable numerically before; see (Guseva et al., 2016) (Gellert et al., 2016) (Guseva et al., 2017) (Guseva et al., 2017) (Mamatsashvili et al., 2018) (Guseva & Tobias, 2023) for works using this code.

Channel flow code with control

MPI-parallelized high-performance code for channel flow in Fortran 90, building upon an older formulation by Flores and Jimenez (2006), with the opposition control strategy. The code solves coupled equations for vorticity and wall-normal velocity using pseudo-spectral direct numerical simulations, with Fourier-Compact-Finite-Differences spatial discretization, Runge-Kutta timestepping scheme and constant mass flux. The opposition control strategy is implemented as a proportionality between the boundary conditions on velocity at the wall and velocity at a given location above it. See (Guseva & Jiménez, 2022) for more details. The code is available here.

Linear stability of turbulent channel flow

A set of subroutines to determine linear stability of a turbulent mean velocity profile in a channel flow, with the opposition control strategy. The code solves an eigenvalue problem for coupled equations for vorticity and wall-normal velocity using pseudo-spectral Fourier-Chebyshev spatial discretization. The opposition control strategy is implemented as a proportionality between the boundary conditions on velocity at the wall and above it. See (Guseva & Jiménez, 2022) and the references therein.