MathGroup Archive 2006

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

Search the Archive

Solving a PDE with fixed endpoints

  • To: mathgroup at
  • Subject: [mg70756] Solving a PDE with fixed endpoints
  • From: "Charlie Brummitt" <cbrummitt at>
  • Date: Thu, 26 Oct 2006 02:38:58 -0400 (EDT)


I am trying to solve PDE's with fixed endpoints. That is, my function
is u[x,t], and I fix u[0,t] = A and u[L,t]=B, where A and B are

However, when I solve these PDE's out to, say, time t = 8,000, I find
that the endpoints don't always stay fixed. Instead, they sometimes
"drift" away from A and B. An example is shown in the following JPG

3D plot of times 7000 through 8000 and plot of u at the final time.
Note that the endpoints should be fixed at 0.1 on both ends, but the
left end (u(0,t)) drifts downward to -120.

Here is the initial condition, showing that the endpoints should be
fixed at 0.1 on both ends:

Should I be circumspect of the solution that NDSolve gives me if the
endpoint condition isn't satisfied correctly for all time? Is there
any way I could fix this?

I think this phenomenon of the end points drifting away from A, B
occurs whenever there is a constant term in the equation. That is,
du/dt = const. + (rest of equation). Thus, this constant term
increases the value of u at every spatial point, which makes it hard
(impossible?) to satisfy the endpoint conditions.

Thanks for any ideas/suggestions you may have,

  • Prev by Date: DownValues and Cases
  • Next by Date: Re: Defining two functions at once?
  • Previous by thread: Re: DownValues and Cases
  • Next by thread: MathModelica System Designer Professional now available