Hello, Daniel, On Mon, Feb 9, 2009 at 6:29 PM, Daniel Wheeler <[email protected]> wrote: > > If I replace the var**3 with h**3 the system seems stable, but uninteresting. > > ##eq = TransientTerm() == -DiffusionTerm((var**3, 1)) + > DiffusionTerm(-exp(-var) * var**3 + 3 * b / var) > eq = TransientTerm() == -DiffusionTerm((h**3, 1)) + > DiffusionTerm(-exp(-var) * h**3 + 3 * b * h**3 / var**4) > > In the email you wrote above you have: > > \partial_t \phi=-\nabla\{\phi^3\nabla[\triangle \phi > -\partial_\phi f(\phi)]\} > > where f(\phi) = \frac{b}{2\phi^2}-e^{-\phi} > > Why did you replace the phi^3 with a constant?
Yes, the above equations are correct. I just changed the name of "\phi" variable with "h" in the code, I dont't hink that I replaced phi^3 with h. In the article: http://prola.aps.org/abstract/PRL/v87/i1/e016104 eq. 3 is exact the above equation, and fig. 5 represents 1D simulations of that equation. What I am trying to do is to simulate 2D case. Thank you, Ionut
