3
$\begingroup$

Suppose I had the following problem:

$U_{tt}=U_{xx}+U_{yy}$ in $\Omega=[0,1]\times[0,1]$
$U(x,y,0)=f(x,y)$
$U_{t}(x,y,0)=g(x,y)$
$U=0$ on $\partial \Omega$

I know that there is an explicit finite difference scheme to solve this problem of the form:

$\frac{U(i,j,k+1) - 2U(i,j,k)+U(i,j,k-1)}{\Delta t^2} = \frac{U(i+1,j,k) - 2U(i,j,k)+U(i-1,j,k)}{\Delta x^2} + \frac{U(i,j+1,k) - 2U(i,j,k)+U(i,j-1,)}{\Delta y^2}$

by using a centered finite difference in time. My mind is thinking that it's possible to simply change the k's to k+1's to obtain an implicit scheme, but I haven't worked out the ensuing truncation error or stability analysis. My guess is that it's likely to be wrong... How do I derive an implicit scheme for this PDE?

Improve this question
$\endgroup$
1
  • 2
    $\begingroup$ Why do you want to use an implicit scheme for this problem? Please remember what the FAQ says: You should only ask practical, answerable questions based on actual problems that you face. $\endgroup$
    – David Ketcheson
    Jan 16, 2012 at 20:06

1 Answer 1

Reset to default
6
$\begingroup$

There is not much point using an implicit method for pure wave propagation because you have to resolve phase to have an accurate method. If you have a hyperbolic system in which some waves are very stiff (not interesting except for their influence on evolution of a slow manifold), you might want an implicit method. It is fairly problem-dependent whether you want an A- or L-stable method, a strong stability preserving method, or a symplectic method. Unfortunately, the order conditions are mutually exclusive such that it is not possible to have all attractive properties in one integration scheme. You can also consider IMEX methods.

Improve this answer
$\endgroup$
4
  • $\begingroup$ if you want to increase your timestep while preserving or decreasing your spatial resolution, wouldn't that also motivate the need to use an implicit method? $\endgroup$
    – Paul
    Jan 16, 2012 at 15:30
  • 2
    $\begingroup$ The point is that no general Eulerian method can resolve phase without resolving CFL ($c \Delta t \sim \Delta x$ where $c$ is the wave speed). If you want to exceed CFL, you have to give up accuracy for some components. If all your waves travel the same speed, you can't gain much. If you have stiff waves or dispersive waves, an implicit method can be an advantage. So the answer is problem-dependent, but your model problem does not offer much potential for implicit methods. $\endgroup$
    – Jed Brown
    Jan 16, 2012 at 15:49
  • 1
    $\begingroup$ Note that for linear standing wave problems (with sources that have few frequency components), you do not need a time discretization. Instead, you can solve the time-periodic problem directly, although implicit solvers become challenging at high frequency, see scicomp.stackexchange.com/questions/40/… and scicomp.stackexchange.com/questions/619/…. $\endgroup$
    – Jed Brown
    Jan 16, 2012 at 15:53
  • $\begingroup$ can you please look at my question (a bit related) with no boundary condition? scicomp.stackexchange.com/questions/41366/… $\endgroup$
    – NotaChoice
    May 19 at 16:47

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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