<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=us-ascii">
<META content="MSHTML 6.00.6000.17080" name=GENERATOR><LINK 
href="BLOCKQUOTE{margin-Top: 0px; margin-Bottom: 0px; margin-Left: 2em}" 
rel=stylesheet></HEAD>
<BODY style="FONT-SIZE: 10pt; MARGIN: 10px; FONT-FAMILY: verdana">
<DIV><FONT face=Verdana size=2>Dear all,</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV>I want to numerically solve a strongly nonlinear fourth-order equation 
(2-D), which is used to describe the dynamics of a liquid film.</DIV>
<DIV>please find the form of the&nbsp;equation below or in the attachment, The 
thickness of the&nbsp;film,&nbsp; H(X,Y,T) is the&nbsp;function to be solved, C, 
G, A0 are constant parameters.</DIV>
<DIV>&nbsp;</DIV>
<DIV><IMG src="cid:__0@Foxmail.net"></DIV>
<DIV>&nbsp;</DIV>
<DIV>I wrote a PETSc programs for this problem, using CN method in time. 
However, it does not work well. In some cases, there are some negativevalues in 
the solution (Obviously, the thickness of the liquid film could&nbsp;not be 
negative.) While in some other cases,&nbsp;the solution remains&nbsp;as the 
initial condition, looking like the governing equation doesn't work at 
all.&nbsp;&nbsp;</DIV>
<DIV>&nbsp;</DIV>
<DIV>I think it is because that the solution of this equation may form infinite 
fradient when the equation is developed.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Some literatures reported that the similar problem was solved successfully 
using the ADI (Alternating Direction Implicit) method. Unfortunately, I found 
that the current PETSc doesn't support ADI.</DIV>
<DIV>&nbsp;</DIV>
<DIV>I also read a paper in which the similar equation was solved by using the 
CN scheme in time. Besides, the author employed the modified second-order upwind 
difference method to handle the nonlinear terms due to the inability of the 
centered differences in space.</DIV>
<DIV>&nbsp;</DIV>
<DIV>I prefer to use PETSc because&nbsp;this powerful toolkit enables easy 
parallel computation of PDE since I am not familar with MPI. Could anyone please 
give me some suggestions on how I can solve the above equation within the 
framework of PETSc?</DIV>
<DIV>&nbsp;</DIV>
<DIV>Thanks very much in advance!</DIV>
<DIV>&nbsp;</DIV>
<DIV>Feng-Chao Wang</DIV>
<DIV><FONT face=Verdana size=2></FONT>&nbsp;</DIV>
<DIV align=left><FONT face=Verdana color=#c0c0c0 size=2>2011-12-28 
</FONT></DIV><FONT face=Verdana size=2>
<HR style="WIDTH: 122px; HEIGHT: 2px" align=left SIZE=2>

<DIV><FONT face=Verdana color=#c0c0c0 
size=2></FONT>&nbsp;</DIV></FONT></BODY></HTML>