When central finite differencing gives complex values for a real solution!

References

• Sachenko , P , Schneibel , JH and Zhang , W . 2002 . Kinetics of thermal grain boundary grooving for changing dihedral angles . Phil. Mag. , A82 : 815 – 829 .
   Google Scholar
• Zhang , Z , Sigle , W , Phillipp , F and Rühle , M . 2003 . Direct atom-resolved imaging of oxides and their grain boundaries . Science , 302 ( 5646 ) : 846 – 849 .
   PubMed Web of Science ®Google Scholar
• Mullins , WW . 1957 . Theory of thermal grooving . J. Appl. Phys. , 28 : 333 – 339 .
   Web of Science ®Google Scholar
• Broadbridge , P . 1989 . Exact solvability of the Mullins nonlinear diffusion model of groove development . J. Math. Phys. , 30 : 1648 – 51 .
   Web of Science ®Google Scholar
• Broadbridge , P and Tritscher , P . 1994 . An integrable fourth-order nonlinear evolution equation applied to thermal grooving of metal surfaces . IMA J. Appl. Math. , 53 : 249 – 265 .
   Web of Science ®Google Scholar
• Tritscher , P and Broadbridge , P . 1995 . Grain boundary grooving by surface diffusion: An analytic nonlinear model for a symmetric groove . Proc. Roy. Soc Lond. A , 450 : 569 – 587 .
   Web of Science ®Google Scholar
• Fornberg , B . 2006 . A pseudospectral fictitious point method for high order initial-boundary value problems . SIAM J. Sci. Comp. , 28 : 1716 – 1729 .
   Web of Science ®Google Scholar
• Fujita , H . 1954 . The exact pattern of a concentration-dependent diffusion in a semi-infinite medium. 3 . Textile Res. J. , 24 : 234 – 240 .
   Google Scholar
• Mathworks Inc., MATLAB and SIMULINK Documentation, 2010; available at http://www.mathworks.com/products/
   Google Scholar
• Ha , Y , Kim , Y-J and Myers , TG . 2008 . On the numerical solution of a driven thin film equation . J. Comput. Phys. , 227 ( 15 ) : 7246 – 7263 .
   Web of Science ®Google Scholar
• Broadbridge , P and Goard , JM . 2010 . Grain boundary evolution with time dependent material properties . J. Eng. Math. , 66 parts 1–3 : 87 – 102 .
   Web of Science ®Google Scholar
• Broadbridge , P , Triadis , D and Hill , JM . 2009 . Infiltration from supply at constant water content: an integrable model . J. Eng. Math. , 64 part 2 : 193 – 206 .
   Web of Science ®Google Scholar
• Triadis , D and Broadbridge , P . 2010 . Analytical model of infiltration under constant-concentration boundary conditions . Water Resour. Res. , 46 part 3 : W03526
   Web of Science ®Google Scholar
• White , I and Broadbridge , P . 1988 . Constant Rate Rainfall Infiltration: A Versatile Non-linear Model. 2. Applications of Solutions . Water Resour. Res. , 24 : 155 – 162 .
   Web of Science ®Google Scholar
• Miller Jr , W . 1977 . Symmetry and Separation of Variables , Reading, MA : Addison-Wesley .
   Google Scholar
• Maplesoft Inc, Maplesoft Documentation Center. Available at: http://www.maplesoft.com/
   Google Scholar
• Childs , EC . 1969 . An Introduction to the Physical Basis of Soil Water Phenomena , New York : John Wiley & Sons .
   Google Scholar
• Ilic , S , Hackl , K and Gilbert , R . 2009 . Application of the multiscale FEM to the modeling of cancellous bone . Biomech. Model. Mechanobiol. , 9 part 1 : 87 – 102 .
   Web of Science ®Google Scholar
• Bohner , M and Peterson , A . 2001 . Dynamic Equations on Time Scales: An Introduction with Applications , Boston : Birkhuser .
   Google Scholar
• Hydon , PE and Mansfield , EL . 2004 . A variational complex for difference equations . Found. Comput. Math. , 4 part 2 : 187 – 217 .
   Web of Science ®Google Scholar
• Becker , J , Grün , G , Lenz , M and Rumpf , M . 2002 . Numerical methods for fourth order nonlinear degenerate diffusion problems . Appl. Math. , 47 part 6 : 517 – 543 .
   Google Scholar
• Wolfram Inc., Wolfram Mathematica 7 Documentation, 2009; available at http://reference.wolfram.com/mathematica/guide/
   Google Scholar