MathGroup Archive 2006

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

Search the Archive

Re: Reaction-diffusion PDEs

  • To: mathgroup at smc.vnet.net
  • Subject: [mg65676] Re: [mg65660] Reaction-diffusion PDEs
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Thu, 13 Apr 2006 02:19:52 -0400 (EDT)
  • Reply-to: hanlonr at cox.net
  • Sender: owner-wri-mathgroup at wolfram.com

soln=NDSolve[{
          D[u[x,t],t]==0.6u[x,t]-v[x,t]-
              u[x,t]^3+0.0002D[u[x,t],{x,2}],
          D[v[x,t],t]==1.5u[x,t]-2v[x,t]+
              0.01D[v[x,t],{x,2}],
          u[x,0]==Sin[2Pi*x],
          v[x,0]==Sin[4Pi*x],
          u[0,t]==0,v[0,t]==0,
          u[1,t]==0,v[1,t]==0},
        {u[x,t],v[x,t]},
        {x,0,1},{t,0,30},MaxStepSize->0.01][[1]];

Plot3D[Evaluate[u[x,t]/.soln],{x,0,1},{t,0,30}];

Plot3D[Evaluate[v[x,t]/.soln],{x,0,1},{t,0,30}];


Bob Hanlon

> 
> From: Skirmantas <skirmantas.janusonis at yale.edu>
To: mathgroup at smc.vnet.net
> Subject: [mg65676] [mg65660] Reaction-diffusion PDEs
> 
> Hello,
> 
> I'd like to numerically solve the following reaction-diffusion system:
> 
> {D[u[x,t],t]==0.6u[x,t]-v[x,t]-u[x,t]^3+0.0002D[u[x,t],x,x], D[v[x,t],t]
==1.5u[x,t]-2v[x,t]+0.01D[v[x,t],x,x], u[x,0]==Sin[2*Pi*x], v[x,0]==Sin
[4*Pi*x]}, where 0<=x<=1, 0<=t<=30, assuming zero-flux boundary 
conditions (at x=0 and x=1). I can set up my own discretization and iteration 
with the For operator, but I suspect NDSolve may be able to to do it 
automatically. I'd appreciate any advice.
> 
> Skirmantas
> 
> 


  • Prev by Date: Re: Windowed standard deviation
  • Next by Date: Re: FindFit / NonlinearFit Problems
  • Previous by thread: Reaction-diffusion PDEs
  • Next by thread: A conditional random number generation problem (please help me!)