On Mon, Feb 9, 2009 at 9:13 PM, Daniel Wheeler
<[email protected]> wrote:
>
> In the code you sent h**3 is being used as one of the fourth order
> diffusion coefficients. In the code the value of h is fixed to be 5.
> When you wrote the equations with tex, in one of the earlier emails,
> the diffusion coefficient is phi**3, which is the solution variable.
> Have I misunderstood something?
In the previous post I gave also the link to an article which should clarify
everything (equations, 1d case, ...). Maybe you can find some time to
take a look.
> The tex equation that you wrote is:
>
> \partial_t \phi=-\nabla\{\phi^3\nabla[\triangle \phi
> -\partial_\phi f(\phi)]\}
>
> The coded equation that you wrote is:
>
> eq = TransientTerm() == -DiffusionTerm((h**3, 1)) +
> (fDerivative.getFaceGrad() * h**3).getDivergence()
>
> and h is set to 5.0. These are different, hence my confusion.
ok, I am sorry for my mistake. Anyway, I still trying to get the code running.
Thanks,
Ionut
