[SciPy-User] SciPy and Diffusion Equation

[Jonathan Guyer]

FiPy can indeed handle problems like this in 1D, 2D & 3D and with non-linear coefficients and arbitrary combinations of terms and arbitrary combinations of equations on multiple solution variables.

See http://www.ctcms.nist.gov/fipy/examples/README.html for examples.

FiPy is not "under the scipy umbrella" but happily sits alongside it.

- Jon (a FiPy developer)

On Jan 16, 2011, at 2:35 PM, Johann Rohwer wrote:

> FiPy (http://www.ctcms.nist.gov/fipy/) is a general PDE solver written in 
> Python, based on the finite volume approach. It should be able to handle 
> diffusion problems and a number of examples are available from the website.
> 
> Johann
> 
> On 16/01/2011 21:29, Lorenzo Isella wrote:
>> Dear All,
>> I hope that this is not too off-topic.
>> I wonder if there is any tool under the scipy umbrella to deal with
>> diffusion problems. I have tackled these problems with proprietary
>> software in the past (e.g. Comsol Multiphysics), but I wonder if there
>> are other (simple) pythonic options.
>> To fix the ideas, consider c(x,t) as a time-dependent concentration
>> profile and the equation
>> 
>> \partial_t c= D\partial_x^2 c-\alpha c +beta\delta(x_0)
>> 
>> i.e. a standard diffusion equation with a source term beta\delta(x_0)
>> (what I mean is that the initial density distribution is a delta
>> centered about x_0, one can probably remove that term from the equation
>> while giving the initial density distribution at t=0).
>> However, it all boils down to a diffusion equation in an interval [0,L]
>> with an initial density profile c(x,t_0). In the case of the equation
>> above, I have an analytical solution, but in general I may want to add
>> sources, non-linear terms etc...so it should be considered as a starting
>> point.
>> Any suggestion is appreciated.
>> Cheers
>> 
>> Lorenzo
>> _______________________________________________
>> SciPy-User mailing list
>> SciPy-User at scipy.org
>> http://mail.scipy.org/mailman/listinfo/scipy-user
> 
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user