Hi everyone,
I am working with RDE system. I tried to search on fipy documentation and
examples about how to write source term in vector form. For example, the code
bellow is to solve 2 RDEs on case 1D.
from fipy import *m = Grid1D(nx=100, Lx=1.)v = CellVariable(mesh=m,
hasOld=True, value=[[0.5], [0.5]], elementshape=(2,))
v.constrain([[0], [1]], where=m.facesLeft)v.constrain([[1], [0]],
where=m.facesRight)M = [[1,0],[0,1]]N = [[1,0],[0,1]]Q = (-v[0] * v[1],
v[0]**2) eqn = TransientTerm(N) == DiffusionTerm([M]) + (Q)vi =
Viewer((v[0],v[1]))
for t in range(20): v.updateOld() eqn.solve(var=v, dt=1.e-3) if
__name__ == '__main__': vi.plot()if __name__ == '__main__':
raw_input("Press <return> to proceed...")
Can anybody tell me how to write the source term 'Q'? The point why I don't
define each variable separately is because I want to work with big RDE system.
Thanks,DhisaVanderbilt
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
