PLUTO
4.4-patch2
|
Compute magnetic field from vector potential. More...
#include "pluto.h"
Functions | |
void | VectorPotentialDiff (double *b, Data *d, int i, int j, int k, Grid *grid) |
The function VectorPotentialDiff() computes either staggered or cell-center magnetic field components by suitable central differencing of the vector potential, .
For staggered MHD the convention is the following:
Ax1[i,j,k] –> Ax1(0++)
Ax2[i,j,k] –> Ax2(+0+)
Ax3[i,j,k] –> Ax3(++0)
For cell-centered MHD vector potential is compute at the cell-center. In the case of staggered MHD, the position of A is edge-centered and it is shown below:
______________________ / /| / / | / z face / | / Ax | / / | / / Az ----------Ay----------- | | | | | | y | | | face | | | / Az x face Az / | | Ax | | / | | / | |/ ----------Ay-----------
change D. Mukherjee, A. Mignone (dipan) jan. mukhe rjee @unit o.it
Assign face- or cell-centered magnetic field by differentiating the vector potential.