Just in case the figure is not showing : Le mardi 3 novembre 2020 à 18:49:19 UTC+1, Henri Colaux a écrit :
> > Dear devs, > > I am writing a code that simulates the diffusion of magnetisation from > proton to proton inside a solid. This process behaves like a standard > diffusion program (but with an equilibrium term) so your code have been > very helpful and I got some great results so far, but I stambled into some > problem whose solution is above my knowledge. > > I am trying to simulate a « drill-core » sort of structure with two > different components - denoted « domains », one being a small slice inside > a much more abundant one, as pictured in the following schematics : > > > I figured that this system is equivalent to producing a 2D array using > Grid2D as picture of the right of this figure, and ensuring that there is > no magnetisation transfer (i.e. no flux) taking place at the border of the > system. To account for the rotation centre and the cylindrical symmetry, I > have taken the effect of the radius directly in the differential > equation.Yet, for some reasons, it seems like the boundery conditions is > not set correctly in my code, which gave unwanted behaviours. > > I have included a minimum working example of the code I am trying to > simulate. First and formost, could you let me know if I do everything right > for what I am trying to simulate ? > > Cheers ! > > Henri Colaux > -- To unsubscribe from this group, send email to [email protected] View this message at https://list.nist.gov/fipy --- To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].
