Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

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

Search the Archive

Re: Boundary condition of diffusion equation in a sphere

  • To: mathgroup at smc.vnet.net
  • Subject: [mg126931] Re: Boundary condition of diffusion equation in a sphere
  • From: Alexei Boulbitch <Alexei.Boulbitch at iee.lu>
  • Date: Mon, 18 Jun 2012 05:45:29 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com

Hi there!

I am solving a diffusion equation in sphere, can anybody tell me how to write this boundary condition in mathematica

c[x,0]=1, 0<x<1

I tried c[x,t]/.0<x<1 and c[x/;0<x<1,t] but they don't work out.

Thank you so much if you could shed some light on it

You may do it this way:

ss = NDSolve[{x*D[u[t, x], t] == x*D[u[t, x], x, x] + 2 D[u[t, x], x],
     u[0, x] == 1, u[t, 0.1] == Cos[t], u[t, 1] == 1},
   u, {t, 0, 2}, {x, 0.1, 1}, MaxStepSize -> 0.001][[1]]


Check this after:

Plot3D[Evaluate[u[t, x] /. ss], {t, 0, 10}, {x, 0.1, 5},
 PlotRange -> All,
 AxesLabel -> Map[Style[#, 16, Italic, Red] &, {x, t, u}]]


You need more boundary conditions here than you have specified above. Namely, at u[t,0] I fixed it for Cos[t] and at the outer boundary of the sphere x=1 I fixed u=1. The consistency requirements for the boundary conditions are here too rigid due to the specific initial condition that you have. I started from x=0.1 instead of x=0 in order to avoid the divergence. It might be chosen closer to zero though. You may try this for the beginning.

Have fun, Alexei



Alexei BOULBITCH, Dr., habil.
IEE S.A.
ZAE Weiergewan,
11, rue Edmond Reuter,
L-5326 Contern, LUXEMBOURG

Office phone :  +352-2454-2566
Office fax:       +352-2454-3566
mobile phone:  +49 151 52 40 66 44

e-mail: alexei.boulbitch at iee.lu






  • Prev by Date: Re: Add scale bar to Image
  • Next by Date: Re: Structure of "identical" data not equal in size
  • Previous by thread: Boundary condition of diffusion equation in a sphere
  • Next by thread: modulo solving lacking domain?