I would consult the FiPy online docs, but with the shutdown I'm having a hard
time.
I have simple analytic forms for a 2D vector flow field that I want to use in
terms such as
ExplicitUpwindConvection etc.
What is the best way to do this? I'm thinking something like this:
def V(x,y):
"" example dummy function """
return [x,y] # is this the right return type? I doubt it.
mesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny)
x = mesh.cellCenters[0]
y = mesh.cellCenters[1]
vel = CellVariable(name = "velocity field",
mesh = mesh,
rank = 1)
vfield = V(x,y)
vel.setValue(vfield)
... ExplicitUpwindConvectionTerm(coeff=vel) ...
# What, exactly, does coeff expect as legal input?
Thanks
Richard
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
