Drexel dragonThe Math ForumDonate to the Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.math.mathematica

Topic: bad PDE solution...
Replies: 1   Last Post: Jan 14, 2013 11:27 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Fred Bartoli

Posts: 6
Registered: 7/1/09
bad PDE solution...
Posted: Jan 14, 2013 12:01 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


As a first simple test case of a more complicated pb I want to solve the
plane diffusion equation.
It seems NDSolve gives a bad answer to the simple step case but I can't
find what's happening...

(* Set up the equation: *)

eq1=D[h[t,x],x,x]- k1 D[h[t,x],t]==0

const=k1-> 10000

(*"solve" it...*)

sol = NDSolve[
{D[h[t, x], {x, 2}] - k1 D[h[t, x], t] == 0, h[0, x] == 0,
h[t, 1] == 0, h[t, 0] == 1 - Exp[-10^6*t]} /. const,
h, {t, 0, 10^-3}, {x, 0, 0.01}][[1]]

right one"]

(* Test the Erfc[x/(2Sqrt[t/k1])] solution *)
Release[Hold[D[h[t,x],x,x]- k1*D[h[t,x],t]]/.

I hope I'm doing something wrong, but I can find what...


Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum 1994-2015. All Rights Reserved.