Re: Problem with string script input Gmsh

2016-07-27 Thread Michael Ghijs
UPDATE: both files (.gmsh and .msh) are present, Popen() could not find them though. At line 226-229 in gmshMesh.py: while True: print ["gmsh", geoFile] + gmshFlags + ["-o", mshFile] p = Popen(["gmsh", geoFile] + gmshFlags + ["-o",

Re: Problem with string script input Gmsh

2016-07-27 Thread Guyer, Jonathan E. Dr. (Fed)
I think it's possible that Popen can't find gmsh, rather than not being able to find your .geo or .msh file. I'm not in a position to diagnose right now, but hopefully somebody familiar with Windows can help you diagnose. My approach would be to write some simple scripts, separate from FiPy, to

Re: Complex PDE Coupling; Extending to a Further PDE

2016-07-27 Thread Guyer, Jonathan E. Dr. (Fed)
Q1: You should be able to solve some equations coupled and others uncoupled. You'll need to sweep them all together. Better would be to figure out how to couple them all. I don't recommend doing them uncoupled, I just couldn't quickly get it to work. Q2: In general, the coupled forms should

Re: How to impose faceGrad.constrain with a numpy array ?

2016-07-27 Thread Guyer, Jonathan E. Dr. (Fed)
Try _Az.faceGrad.constrain( _By.arithmeticFaceValue * mesh._orientedFaceNormals, where = mesh.facesLeft) _Az.faceGrad.constrain(-_By.arithmeticFaceValue * mesh._orientedFaceNormals, where = mesh.facesRight) _Az.faceGrad.constrain(-_Bx.arithmeticFaceValue * mesh._orientedFaceNormals,

away for awhile

2016-07-27 Thread Guyer, Jonathan E. Dr. (Fed)
I'm going to be on leave without any email until the middle of August. Daniel should be back in the office next week. Replies from us may be slow for awhile. ___ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY:

Problem with string script input Gmsh

2016-07-27 Thread Michael Ghijs
Hi everyone, I am experiencing /troubles /with the *fipy.Gmsh2D function with* *string*-wise script *input* in the /circle diffusion example/. Please find the script in attachment (and note that I applied a little trick to circumvent the

Re: Diffusion-Convection-Reactive Chemisty

2016-07-27 Thread Guyer, Jonathan E. Dr. (Fed)
I recommend you work through and check your units in all equations. I don't have any intuition for what's a big value and what's small when we're talking kmol. There aren't any units at all specified for D, k, etc. Is your solution converged when it finishes `while res > 1.e-2`? How did you