On Jun 6, 2007, at 6:34 AM, Ralf Banisch wrote:
can I use fipy to solve the following set of PDE's:
da/dt = s(a^2/b + b_a) - a*r_a + D_a*d^2a/dx^2
db/dt = s*a^2 - b*r_b + D_b*d^2b/dx^2
where a and b are the unknown variables and s, r_a, r_b, D_a and
D_b are given constants. a and b are one-dimensional functions a = a
(x,t) in this example. If possible, I would like to use fipy for
solving the set above in the case of two-dimensional functions a = a
(x,y,t) on a two-dimensional grid as well.
FiPy is well suited to solving this sort of problem. Our user manual
contains specific guidance on how to represent and solve these sorts
of equations.
In particular, please see FAQ 6.1 and Example 7.1.
If you have not already downloaded it, you can find our manual at
http://www.ctcms.nist.gov/fipy/download/fipy.pdf
--
Jonathan E. Guyer, PhD
Metallurgy Division
National Institute of Standards and Technology
<http://www.metallurgy.nist.gov/>
