Re: Octave PDE initial value problem solvers

[Geordie McBain]

There are some good examples of solving nonlinear two-point boundary
value problems using orthogonal collocation in Octave in the book

        J. B. Rawlings & J. G. Ekerdt 2002 Chemical Reactor Analysis 
        and Design Fundamentals, Nob Hill.  ISBN 0-615-11884-4

I think you might even be able to find the Octave code on the Web
somewhere.

I don't recall whether they also treated 1-D parabolic initial--boundary
value problems, but it's easy enough to reduce those to a sequence of
the former if not, using finite differences in time.

Geordie McBain
www.aeromech.usyd.edu.au/~mcbain


On Fri, 2005-10-28 at 12:31 +0800, address@hidden wrote:
> Hello,
> 
> I need to put together a model that describes how moisture and salts
> move through a layer of unsaturated mud as it dries (by evaporation off
> the top, and by drainage out the bottom).  I'll start by trying to get a
> moisture transport model working, then will worry about coupling the
> salt equation.  At some point I'll also need to concern myself with
> volume changes as the mud dries out.  I'm going to start by trying to
> numerically solve the "Richards Equation", which is a nonlinear
> parabolic equation.  The usual way to solve this seems to be to use a
> Galerkin-type finite element scheme - something I need to read up on.
> 
> Anyway, I'm wondering if anybody knows of publicly available Octave code
> that might be useful to me?  I'll only be working on a 1D model.
> 
> As somebody who used to use Matlab (though in signal processing and
> control, not for solving equations) at uni, I've decided that using
> Octave should get me started faster then trying to use C++ or Delphi...
> 
> Thanks,
> Neil
> 
> 
> 
> -------------------------------------------------------------
> Octave is freely available under the terms of the GNU GPL.
> 
> Octave's home on the web:  http://www.octave.org
> How to fund new projects:  http://www.octave.org/funding.html
> Subscription information:  http://www.octave.org/archive.html
> -------------------------------------------------------------
> 



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------