On 07/22/2014 08:08 PM, Joseph Batley wrote:
Ah, thanks you Robert, It appears I did not quite understand the conversion from strong to weak. I have made the changes and attached the problem file.
I have found a thesis "THEORY BASED DESIGN AND OPTIMIZATION OF MATERIALS FOR SPINTRONICS APPLICATIONS", where the same equations are solved by FEM, and it seems to me that what is in the attached file (*) is ok. Have a look at it as well, please.
(*) updates:
- fixed dw_surface_flux material parameter shape (needs a tensor - compile sfepy with -DDEBUG_FMF in site_cfg.py to catch that)
- fixed the Jc boundary condition (not sure about sign...)
- fixed/simplified interface_FM, interface_NM region definitions - now their faces have opposite normals (is that ok?)
- added more things to output for debugging
It looks promising, the post process shows \mu_{c} = \mu_{\uparrow} + \mu_{\downarrow} looks correct so the normal charge current side works, the \mu_{s} = \mu_{\uparrow} - \mu_{\downarrow} which is the spin accumulation however should simply show a peak around the interface which decays with different lengths in each domain.
After the above changes I still do not see a peak. Instead, mu_up and mu_down are the same, so \mu_{s} is zero.
You have more debugging ahead :)
r. PS: you may need to tweak the mesh filename. PPS: one.val material is not needed, you can use simply dw_laplace.i.Omega(d, mu_up)...
And yes you're correct, the equations oars the same in both domains, it is only the material parameters that change.