Re: partial differential equation of diffusion
- To: mathgroup at smc.vnet.net
- Subject: [mg88377] Re: partial differential equation of diffusion
- From: "Ali YILDIRIM" <ayildirim10 at hotmail.com>
- Date: Sat, 3 May 2008 06:17:34 -0400 (EDT)
- References: <BLU129-DS716E7490C4367AC35F1EEA8DB0@phx.gbl> <22d35c5a0805012309s10c21a53me6a645fcb0be5a64@mail.gmail.com>
Dear Jean Marc Gulliet, >>When I put my equation as D[M[r, t], t] == a*(D[M[r, t], r, r] + >>(2/r)*D[M[r, t], r]) - k*M[r, t] in to mathematica >and initial and boundary condition as M[r, 0] == M0, D[M[0, t], t] == 0, >M[R, t] == Ms >It can not solve. The mathematica notebook is given in attachment. >Best regards > > dM/d t =a[d^2 M//dr^2 +(2/r) (dM/dr)]-kM >> > The differentials are partial differentials. >> > a and k are constants. >>> >> > Do you mean that you do not know how to multiply expressions with > Mathematica? If this is the case you should familiarize yourself with > the basis of Mathematica syntax, this will really help you and be a > time saver in the future. > > D[M[r, t], t] == > a*(D[M[r, t], r, r] + (2/r)*D[M[r, t], r]) - k*M[r, t] > >> >Best Regards.. >> > >> > >> Ã?Ä?r.Gör. Ali YILDIRIM >> Gaziantep Ã?niversitesi >> Nizip Meslek Yüksekokulu >> Gıda Teknolojisi Programı >> 27700-Nizip-Gaziantep/Türkiye >> ayildirim at gantep.edu.tr >> > -----------!!!----!!!! >> > Not sure for what this D stands here. I discarded it in the example >> > below. Also you have a spurious backslash. So you could write your >> > equation as >> > >> > D[M[r, t], t] == D[M[r, t], r, r] + 2/r D[M[r, t], r] - k M[r, t] >> > >> >> Boundary conditions are: >> >> >> >> 1-) at t=0, M= Mo >> > >> > M[r, 0] == M0 >> > >> >> 2-) r = 0, dM/dr=0 >> > >> > D[M[0, t], t] == 0 >> > >> >> 3-) at r = R, M = Ms >> > >> > M[R, t] == Ms >> > >> > Then you can use DSolve or NDSolve as in >> > >> > DSolve[{D[M[r, t], t] == >> > D[M[r, t], r, r] + 2/r D[M[r, t], r] - k M[r, t], M[r, 0] == M0, >> > D[M[0, t], t] == 0, M[R, t] == Ms}, M, {r, t}] >> > >> >> >> > > > Regards, > -- Jean-Marc >