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 » sci.math.* » sci.math.num-analysis.independent

Topic: Crank Nikolson scheme for semi linear parabolic type equation
Replies: 1   Last Post: Feb 14, 2013 6:51 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View  
AMX

Posts: 35
Registered: 8/22/09
Re: Crank Nikolson scheme for semi linear parabolic type equation
Posted: Feb 14, 2013 6:51 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

On Wed, 13 Feb 2013 22:28:51 -0800 (PST),
Sandeep Kumar <searchsandy1712@gmail.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)


> boundary
> conditions Ux(0,t)=Ux(l,t)=Uxxx(0,t)=Uxxx(l,t)=0;
>


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.

For 1D (space) your problem can be written:

?u(x,t)/?t = - ?^2 ?^2u(x,t)/?x^2 + ?(u^3(x,t))/?x - ?u(x,t)/?x

After spatial discretisation RHS will be some expression of
discrete u(x_i,t).

Beside boundary conditions you have to specify initial condition
u(x,0) = f(x)

AMX

--
adres w rot13
Nyrxfnaqre Znghfmnx r-znk@b2.cy



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

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.