One other thing : I tried the « constrain » method with either :
Pol.faceGrad.constrain(((0,),(0,)), where=gr.exteriorFaces) or Pol.faceGrad.constrain(((0,),(0,)), where=gr.facesTop) Pol.faceGrad.constrain(((0,),(0,)), where=gr.facesBottom) Pol.faceGrad.constrain(((0,),(0,)), where=gr.facesLeft) Pol.faceGrad.constrain(((0,),(0,)), where=gr.facesRight) but it does not seem to change anything. Could anyone tell me what I am doing wrong ? Thank you ! H. Le mardi 3 novembre 2020 à 18:50:29 UTC+1, Henri Colaux a écrit : > 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].
