On Wed, 13 Feb 2013 22:28:51 -0800 (PST), Sandeep Kumar <email@example.com> wrote:
> I am studying Cahn-Hilliard equation which is a fourth order > parabolic type equation and wanted to implement this equation > using numerical methods. I am trying to solve this with > Crank-Nikolson scheme, but am not able to formulate the > problem. The equation contains a fourth order derivative( > square of laplacian) and also a laplacian acting on a cubing > term. As per the scheme I have taken the mean of Ui's at n and > n+1 time scale. After the action of laplacian on the cubic > term, I get two terms. Now, I am not able to separate the n and > n+1 term to get the expression into a matrix form. With my > knowledge about CN scheme, we are supposed to separate the n > and n+1 terms. Or, is there another way of solving the semi > linear PDEs? > > The equation is: del u/ del t = - epsilonsquare* > laplacian^2*u+ laplacian(u^3)-lapacian(u)
This (AFAIK) is an initial-boundary value problem. Initial with respect to time and boundary with respect to space. CN can be used for IVP. For BVP you have to use other method, like FDM, FEM or other.