Hello, I'm working in my PhD thesis and at the moment I have to solve the Richards' Equation in 1D.
C(h)*Dh/Dt = D/Dx[K(h)*Dh/Dx] + DK/Dx
with some boundary and initial condition.
I usually work in matlab and I need some point to start. If some of you know a matlab code that solve this kind of problem please tell me something.
Thanks a lot, and greetings.
Hector
Re: Richard's equation
A google search gives me hundreds of possibilities. If searched for "Richards equation matlab " and got 55000 matching pages. If you want to use finite differences then you can look at http://www.usyd.edu.au/su/agric/acpa/software/matlab/matlab.htm. There are several freely available papers that explain the issues involved if you use finite elements . I think it's a good idea to just derive a finite element formulation of the equation on your own and then code it up. That way you will understand why poeple use mixed methods for the problem and how they are implemented.
-- Biswajit