PLUTO Test Problems  4.4-patch2
init.c File Reference

Resistive RMHD rotor test problem in 2 dimensions. More...

Detailed Description

The 2D rotor problem [dZBL03][Mig12] consists of a rapidly spinning disk embedded in a uniform background medium treaded by a costant magnetic field. The initial conditin reads as

\[ \left(\rho, p, v_x, v_y, B_x\right) = \left\{\begin{array}{ll} \left(10, 1,-\omega y,\,\omega x,\, 1\right) & \quad\mathrm{for}\quad r < 0.1 \\ \noalign{\medskip} \left(1,\, 1,\, 0,\, 0,\, 1\right) & \quad\mathrm{otherwise} \end{array}\right. \]

where $\omega = 8.5$ is the angular frequency of rotation. The electric field is initialized to its ideal value while remaining variables are set to zero. The computational domain is the unit square and outflow boundary conditions are imposed everywhere.

A list of tested configurations is given in the following table:

Conf.DIM divB eta
#01 2 (np=1) GLM 1.e-6
#02 2 (np=8) GLM 1.e-6
#03 2 (np=8) GLM 1.e-1
#04 2 (np=1) CT/ARITHMETIC 1.e-6
#05 2 (np=8) CT/ARITHMETIC 1.e-6
#06 3 (np=8) CT/MAXWELL 1.e-6
#07 3 (np=8) CT/MAXWELL 1.e-6

The final solution is computed at t = 0.3. Conf. #06 and #07 are 2D solutions computed in the XY plane and YZ plane. Used to check invariance under rotations.

Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Aug 31, 2019

Reference: