On Mon, Nov 18, 2013 at 4:40 PM, olivier atteia <olivier.att...@ipb.fr>wrote:
> Dear colleague, > thanks a lot for creating FiPy > I begun to use it and it has been very helpfull for me for the transport > questions (i will surely cite fipy in my next papaer!) > > however, i am not a specialist of numerical solvers and i might want to go > too far, but i want to solve the richards equations : unsaturated water > flow, like > > Hi Olivier, I would not use the "var" argument when instantiating the terms since you are only really solving for one variable, "S", in this example. The next step is to reformulate the diffusion term so FiPy understands it. The form of the equation should be \frac{\partial S}{\partial t} = \nabla \cdot \left( K \frac{\partial H}{\partial S} \nabla S \right) So, basically, you need to find find the derivative of H with respect to S and then the diffusion coefficient is just "K * Hprime". Hope that helps. -- Daniel Wheeler
<<image/png>>
_______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]