Function to handle vector potential.
Update the cell-centered vector potential by integrating
Depending on the time-stepping algorithm, we advance the solution in time using the following increments
The cell-centered electric field E is computed by averaging values available at faces or edges, depending on the MHD formulation:
- staggered_mhd: we average the electric field values available at cell edges and previously computed . In this case the function does the job only if called from CT;
- Cell-cent. mhd: we average the face values using upwind fluxes. In this case the function does the job only if called from time stepping algorithms. The flux components of induction equation contain the emf:
,
,
,
Note that the vector potential is NOT used in the rest of the code and serves solely as a diagnostic tool.
- Parameters
-
[in,out] | d | pointer to PLUTO Data structure |
[in] | vp | pointer to an EMF structure, used only with CT |
[in] | state | pointer to State_1D structure containing the electric field components (cell-centered MHD) |
[in] | grid | pointer to an array of Grid structures |
- Attention
- In cylindrical coordinates: PLUTO treats cylindrical (but not polar) coordinates as Cartesian so that
. Nevertheless, the correct right-handed cylindrical system is defined by
. The fact that the
coordinate is missing causes an ambiguity when computing the third component of the EMF in 2.5 D:
which is (wrongly) called
but is, indeed,
. This is not a problem in the update of the magnetic field, since the contribution coming from
is also reversed:
as it should be. However, it is a problem in the definition of the vector potential which is consistently computed using its physical definition. Therefore, to correctly update the phi component of the vector potential we change sign to
(
).
- Authors
- A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)
P. Tzeferacos (petro.nosp@m.s.tz.nosp@m.efera.nosp@m.cos@.nosp@m.ph.un.nosp@m.ito..nosp@m.it)
- Date
- Sep 24, 2012
Definition in file vec_pot_update.c.