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

MathGroup Archive 2009

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

Search the Archive

Re: question about NDSolve for Diffusion equation

  • To: mathgroup at smc.vnet.net
  • Subject: [mg102840] Re: [mg102823] question about NDSolve for Diffusion equation
  • From: "Ingolf Dahl" <ingolf.dahl at telia.com>
  • Date: Sat, 29 Aug 2009 06:31:07 -0400 (EDT)
  • References: <200908280446.AAA28981@smc.vnet.net>

First things first:
Is this really the diffusion equation? With the second derivative in t it
looks more like the wave equation.
But then, look in the MathGroup archive, less than one year back, for the
threads "Diffusion Model using NDSolve - Advice needed", "PDE heat equation
(inconsisten problem)" and "Partial differential equation with evolving
boundary conditions". I think you have stumbled on the same problems as
discussed there, not a mathematical necessity but a bug-like feature of the
solution method used by Mathematica. I hope you will find something
interesting or useful for your problem there. I think your letter is a FAQ,
and I think you should avoid the sharp change of C[0,t] near t=0. You might
try with a ramp in a short time interval. There are also other methods
available. 


Best regards

Ingolf Dahl


> -----Original Message-----
> From: Nam Viet Nguyen [mailto:teivn109 at yahoo.com] 
> Sent: den 28 augusti 2009 06:46
> To: mathgroup at smc.vnet.net
> Subject: [mg102823] question about NDSolve for Diffusion equation
> 
> Hello,
> I am trying to use the NDSolve of Mathematica to solve the 
> diffusion equation which has diffusivity dependence to 
> concentration. My work is dealing with phosphorous diffusion 
> into silicon.
>  
> The equation is
> d''C/dt^2 = d( C^2 *dC/dx)/dx
> where C= C(x,t).
> In my problem, the initial conditions should be : C(x,0)=0, 
> C(0,0) = 0, C(0,t) = Cs if t > 0. 
> If I used the NDSolve, there is a error message, the Initial 
> conditions are not consistent.
> I dont known how to go further as I have checked again and 
> again the initial conditions are fine. 
> Your helps anf advice are gratefully acknowledged,
> 
> Nam Nguyen.
> Stuttgart, Germany
> 
> 



  • Prev by Date: Problems encountered with Mathematica
  • Next by Date: problem with integrating an interpolated list
  • Previous by thread: question about NDSolve for Diffusion equation
  • Next by thread: Viewing .mx files (Question)