Re: Integrating field

2012-11-30 Thread A. A. Eftekhari
Dear all, I have a question not too unrelated to this one. Let's say I have a 2D numpy array and I want to assign its value to a cell variable defined over a 2D mesh of the same dimensions. Is there any simple way to do that. n = 5 mesh = Grid2D(nx = n, ny = n) k = rand(n,n) phi = CellVariable(me

Re: Flux=0 constrains on both sizes cause problems with simple 1D diffusion problem

2012-11-08 Thread A. A. Eftekhari
Hi Daniel, In your script, you are solving an steady state problem. Add a transient term, and preferably sweep to get the solution. Perhaps sth like this: equation = (TransientTerm() == diffusion_term + source_term) timestep = steps = 10 for step in range(steps) res = 1.e6 phi.u

Re: Buckley Leverett equation

2012-09-14 Thread A. A. Eftekhari
ition. It is sort of crazy, but an upwind convection term > with a zero velocity will not interact with the boundary (I think, I > didn't think it through carefully). It sort of isn't a boundary value > problem at that point. It's confusing. > > Also I used $\fr

Buckley Leverett equation

2012-09-13 Thread A. A. Eftekhari
Hello all, I'm using fipy to solve flow in porous media problems, and I must say with the new version (3.0) there is no problem in solving single phase multi-component convection diffusion equation coupled with Darcy's law. My difficulty is with solving two-phase flow and specially Buckley-Leveret