5
$\begingroup$

Is the Crank-Nicolson method appropriate for solving a system of nonlinear parabolic PDEs like $\partial u/\partial t - a\Delta u + u^4 = 0$ ?

I tried to apply this method for solving such system but the solution was oscillating (maybe because of a small value of the coefficient of the time derivative) and the implicit Euler method calculates a correct solution.

Improve this question
$\endgroup$

2 Answers 2

Reset to default
10
$\begingroup$

This phenomenon is often called "ringing" and plagues methods that are not $L$-stable. This can be seen in this motivating example from Hairer & Wanner (1999) "Stiff differential equations solved by Radau methods". Consider the equation

$$ \dot y = -50 (y - \cos t) $$

and apply explicit Euler with time step near the stability limit, implicit midpoint (or equivalently for this problem, trapezoid rule, aka. Crank-Nicolson), and implicit Euler. The result, shown below, illustrates that the $A$-stable (but not $L$-stable) implicit midpoint method produces a poor-quality solution.

Figure 1 of Hairer and Wanner (1999)

To avoid this problem when solving stiff systems, you should use an $L$-stable method such as BDF-2, a suitable DIRK, or a Radau method. See Hairer and Wanner's second volume for extensive discussion of this topic.

Improve this answer
$\endgroup$
2
  • $\begingroup$ Do you mean methods for solving ODEs? How these methods can be applied for solving PDEs? Is it required to use the method of lines? $\endgroup$
    – jokersobak
    Commented Nov 3, 2014 at 8:38
  • $\begingroup$ You can discretize in space or time first, or both simultaneously (as in space-time). I think it's typically more useful in terms of software interfaces to use method of lines (discretize in space first), though space-time discretizations have their place when you want local adaptivity simultaneously in space and time. $\endgroup$
    – Jed Brown
    Commented Nov 4, 2014 at 0:35
1
$\begingroup$

Theoretically, if the implicit Euler method works for this equation, Crank--Nicolson scheme should also work. let $\tau$ be the step in time and if we only consider the temporal discretization, the linearized Crank--Nicolson scheme is given by $$ \frac{u^n-u^{n-1}}{\tau} - \frac{a}{2} ( \Delta u^n+ \Delta u^{n-1} ) = -(\frac{3 u^{n-1}-u^{n-2}}{2})^4 .$$

The above scheme is okay under the condition that the regularity of the equation is good. Moreover, spatial approximation does not affect the stability of the scheme and you can use both finite difference and finite element methods to solve it.

Improve this answer
$\endgroup$
5
  • 1
    $\begingroup$ Putting aside the lack of $L$-stability, extrapolation of the reaction term as you propose often results in instability. Even for split implicit methods, convergence is often lost for moderately large time steps, see, e.g., Ropp, Shadid, and Ober (2004). $\endgroup$
    – Jed Brown
    Commented Oct 22, 2014 at 16:50
  • $\begingroup$ @mathsgao Why did you write the term $-\left(\frac{3u^{n-1}-u^{n-2}}{2}\right)^4$ ? What does it mean? $\endgroup$
    – jokersobak
    Commented Oct 23, 2014 at 1:24
  • $\begingroup$ @mathsgao Why it is not the following scheme? $\frac{u^n - u^{n-1}}{\tau} - \frac{a}{2}(\Delta u^n + \Delta u^{n-1}) + \frac{1}{2}((u^n)^4 + (u^{n-1})^4)$ $\endgroup$
    – jokersobak
    Commented Oct 23, 2014 at 4:32
  • $\begingroup$ Sure, you can use it. But it is nonlinear, i.e., $F(u^n)=0$ has to be solved at each time step. $\endgroup$
    – mathsgao
    Commented Oct 23, 2014 at 15:59
  • $\begingroup$ If the regularity of the solution is very good, C-N scheme with linearization should be stable. However, if the solution of the PDE blows up, C-N scheme maybe not good. @jokersobak, would you like to test these schemes with an artificial example which admits smooth solution? $\endgroup$
    – mathsgao
    Commented Oct 23, 2014 at 16:07

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.