On May 23, 2012, at 6:21 AM, <[email protected]>
<[email protected]> wrote:
> In fact, I just didn't setup F to work witth the spatial dimension, i.e. F
> is for me a chemistry solver and does not know about space.
OK. I was confused when you said that y was an array of CellVariable. Are there
different values of y for each point in space or just for each PDE or are there
just a set of values of y that are used for all PDEs for all points in space?
> \phi and y are like components and species (N\phi is the number of components
> and Ny is the number of species, and we have in general N\phi > 1 and always
> Ny > N\phi because of chemical reactions)
Oh. So is this a fair representation of what you want for N\phi = 2 and Ny = 3?
\frac{\partial{\phi_1}}{\partial t} = \nabla ( D(y_1) \nabla y_1) + \nabla (
D(y_2) \nabla y_2) + \nabla ( D(y_3) \nabla y_3)
\frac{\partial{\phi_2}}{\partial t} = \nabla ( D(y_1) \nabla y_1) + \nabla (
D(y_2) \nabla y_2) + \nabla ( D(y_3) \nabla y_3)
That can't be enough because then \phi_1 and \phi_2 would be identical.
Please give the simplest set of toy equations for your PDEs and for F that
represent what you want.
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
