MathGroup Archive 2011

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Numerical Solver for 2D/3D (Reaction-)Diffusion equation

  • To: mathgroup at smc.vnet.net
  • Subject: [mg115346] Re: Numerical Solver for 2D/3D (Reaction-)Diffusion equation
  • From: Oliver Ruebenkoenig <ruebenko at wolfram.com>
  • Date: Sat, 8 Jan 2011 03:39:20 -0500 (EST)

On Tue, 4 Jan 2011, Alois wrote:

> I'm trying to solve numerically a 2d/3d reaction-diffusion equation
> using NDSolve. The system can be solved for one spatial dimension (+
> time dimension), but it does not converge for two spatial dimensions
> (x,y,t). Actually, I'm trying to use polar coordinates (r, phi, t),
> because the problem has a point source.
>
> I'm using NDsolve with the following method:
> Method -> {"MethodOfLines", "SpatialDiscretization" ->
> {"TensorProductGrid", "Coordinates" -> myGrid}
>
> which fails. Can you recommend a better method to solve the problem in
> 2D (and eventually in 3D) ?
>
> Kind regards,
>   Alois
>
>
>
>

Alois,

things like


eqn = \[Rho]*cp*D[u[t, x, y], t] - \[Lambda]*
     D[u[t, x, y], x, x] + -\[Lambda]*D[u[t, x, y], y, y] ==
   f + r*u[t, x, y]*If[x == p1X, y == p1Y, 1, 0]

NDSolve[Join[{eqn},
    bc, ic], u, {t, 0, tEnd}, {x, 0, xEnd}, {y, 0, yEnd},
   Method -> {"MethodOfLines",
     "SpatialDiscretization" -> {"TensorProductGrid",
        MaxPoints -> {35, 35}, MinPoints -> {35, 35}
       }}
   ];

work in principal. Now, the exact details depend on the equation 
parameters. Note that the Max and Min points need to be chosen 
such that the tensor grid generates a point {p1X,p1Y}.

Hth,
Oliver


  • Prev by Date: Mathematica will be a new feature of HTML6
  • Next by Date: Re: Reduce in Mathematica 5 vs Mathematica 8 (2nd problem)
  • Previous by thread: Re: Numerical Solver for 2D/3D (Reaction-)Diffusion equation
  • Next by thread: Re: Best Way to Copy from Mathematica to Powerpoint?