I'm working on a finite element simulation of the magnetostatic magnetic vector potential in 3D, with the following geometry:.
The internal structure is a hollow cored-rectangular current loop. The outer box is the spacial domain, set to be approximately 10x the characteristic length scale of the problem. As a result I expect magnetic flux to be wholly contained within this box. My equation to be solved is the vector Poisson equation in 3D: $$\nabla^2\vec A=-\mu\vec J$$
My question then is with regards to the boundary conditions upon the faces. For faces 11-16 (the outer box) I choose Dirichlet zero conditions which correspond to the flux-enclosed approximation mentioned.
Faces 1-10, i.e the faces of my coil I'm unsure of. Assuming A is non-divergent, I expect the magnetic vector potential to be continuous across boundary but I also expect the normal derivative of A to be discontinuous by the surface current density K (D.J. Griffiths (1999),"Introduction to Electrodynamics", Upper Saddle River, NJ, Prentice Hall). This corresponds to a tangential discontinuity in the magnetic flux due to some K. However I have noticed as I work through the literature that many papers seem to ignore this component and simply require A to be continuous e.g:
N. Demerdash, T. Nehl and F. Fouad, "Finite element formulation and analysis of three dimensional magnetic field problems," in IEEE Transactions on Magnetics, vol. 16, no. 5, pp. 1092-1094, September 1980. doi: 10.1109/TMAG.1980.1060817
N. A. Demerdash, F. A. Fouad, T. W. Nehl and O. A. Mohammed, "Three Dimensional Finite Element Vector Potential Formulation of Magnetic Fields in Electrical Apparatus," in IEEE Transactions on Power Apparatus and Systems, vol. PAS-100, no. 8, pp. 4104-4111, Aug. 1981. doi: 10.1109/TPAS.1981.317005
The part I am unsure of is how to find the surface current density (assuming this is important). I know for my rectangular-cross sectioned coil with winding height z, inner length L_i, outer length L_o with N turns carrying a current I, the average current density through a winding section is: $ j_{0} = \frac{IN}{z(L_{o} - L_{i})/2} $.
I also know the surface current density is also likely to vary with my faces in accordance with the geometries.