Posted by ADSW1243 on
URL: http://forum.featool.com/Vector-Poisson-Equation-for-a-3D-Current-Source-tp338p361.html

The (hopefully I know I must be getting annoying!) final question I have is with regards to the coefficients J1 J2 J3.

From trying to do FEM in other methods I have a function (attached with its dependent) which, for a 3D coil with rectangular cross section, can map the current density to nodes (well kind of, honestly the function isn't great because it only works if nodes are symmetrically distributed through the coil and the surface nodes at the corners should definitely be assigned differently. The aim is sqrt(J1^2 + J2^2 + J3^2) = Constant everywhere in the coil, with rotational flow in the XY plane like the initial figure on this post in case you need to test this to answer my question).

Currently it is of the form: [Jx,Jy,Jz] = J_coeff(x,y,z,N_wire,I_mag,mu_o).

Given this setup I have two simple questions:

1. Do I need to split this into 3 functions which only output one of J1, J2, J3?
2. In my geometry for the initial problem, I have two sub-domains: an internal coil of the type discussed above, surrounded by a 'universe' box which forms a typical black-box setup. Currently if I understand the GUI correctly, I set Jx as a constant which seems to apply to both sub-domains. The expectation is that this will only apply to the internal sub-domain, but the values in the outer-sub-domain call the internal as a 3D source. How do I approach ensuring this is achieved? I thought perhaps a different equation is needed for each sub-domain, but I am unsure.




J_coeff.mcoilcorner.m