Lagrangian MHD in 2D and 3D

Summary form only given, as follows. The cell-centered diffusion differencing scheme presented by Morel et al. (1992) has been applied to magnetic diffusion associated with Lagrangian hydrodynamic codes. Thus, the method applies to non-orthogonal meshes. Although the present application involves structured meshes, the method applies equally well to unstructured meshes. Morel´s example of application is to 2D diffusion using Ficke´s law. Thus, a volume integral approach is applied to the divergence operator. In 2D magnetic diffusion symmetry allows the use of an area integral approach involving the field components normal to the area, e.g. A-theta and B-theta. Instead of a divergence of a term proportional to the field gradient a curl of a term proportional to the curl of the field is used. An essential fact that allows this procedure is that the variable theta is ignorable. A benefit of this approach is that the solenoidal property of the magnetic field is automatic. In the case of 3D it is necessary to return to the volumetric integral approach and to use rectangular components of the vector potential. Successful benchmarks have been run in comparison with the 1D code RAVEN. A typical example is that of a metal cylinder being compressed by a magnetic field applied at the outer boundary. So far, the 3D diffusion model has been tested in the orthogonal case and found to preserve the linear, homogeneous solution.