Hello,
I would like to implement "partial derivative of u with respect to x +
partial derivative of v with respect to y = a constant" at boundaries.
Is it possible in Fipy? Thanks in advance.
or maybe something similar as
mesh = Grid2D(dx=dx,dy=dy,nx=nx,ny=ny)
u = CellVariable(mesh=mesh, hasOld=1, rank=0)
v = CellVariable(mesh=mesh, hasOld=1, rank=0)
(u.faceGrad[0] + v.faceGrad[1]).constrain(1., where=mesh.exteriorFaces)
Regards
Ronghai
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
