Jan Verwer - PowerPoint PPT Presentation

Jan verwer
1 / 23

  • Uploaded on
  • Presentation posted in: General

Convergence and Component Splitting for the Crank-Nicolson Leap-Frog Scheme. Jan Verwer. Hairer-60 Conference, Geneva, June 2009. TexPoint fonts used in EMF. Read the TexPoint manual before you delete this box.: A A A A A. Crank-Nicolson Leap-Frog (CNLF). non-stiff. stiff.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

Download Presentation

Jan Verwer

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

Jan verwer

Convergence and Component Splitting

for the Crank-Nicolson Leap-Frog Scheme

Jan Verwer

Hairer-60 Conference, Geneva, June 2009

TexPoint fonts used in EMF.

Read the TexPoint manual before you delete this box.: AAAAA

Jan verwer

Crank-Nicolson Leap-Frog (CNLF)





Usually IMEX-Euler for

Jan verwer

Contents of this talk

CNLF is a two-step IMEX scheme. Used for PDEs in

CFD (method of lines). Non-stiff term then represents

convection and the stiff term diffusion + reactions.

This talk is about an alternative use of CNLF:

  • A splitting (convergence) condition justifying a

  • wider class of splittings than normally seen in CFD

  • As an example, component splitting for 1st order

  • Maxwell-type wave equations

  • Two numerical illustrations of component splitting

Jan verwer

Consistency of CNLF

We always think of semi-discrete systems

but suppress for convenience the spatial mesh size

Further, order terms like

are always supposed to be derived and valid for

Jan verwer

Consistency of CNLF

Just for convenience we neglect spatial errors.


Then the local truncation of CNLF satisfies



In CFD applications this splitting (convergence) condition is mostly satisfied!

Jan verwer

Consistency of CNLF

For the IMEX-Euler scheme

the splitting (convergence) condition features in the

same way. That is, if


uniformly in the spatial mesh size

Jan verwer

Convergence of CNLF

Hence, if

and assuming stability, CNLF with Euler start will converge with order two uniformly in the spatial mesh width!

Q: is this common splitting (convergence) condition also necessary for 2nd – order convergence?

Jan verwer

Numerical counter example

Semi-discrete 1st-order wave equation, with a splitting such that

is violated (splitting details later).

We let

1st order

(i) The common splitting condition is not necessary for 2nd order CNLF

convergence. What is the right condition?

(ii) But why only 1st order when

IMEX-Euler is used to start up?

2nd order

-o- : Exact (or CN) start

-*- : IMEX-Euler start

Jan verwer

A new splitting (convergence) condition

First the linear case:


  • Thm. Assume stability and condition (n). Then, uniformly in h,

  • IMEX-Euler is 1st-order convergent

  • CNLF with IMEX-Euler startis 1st-order convergent

  • CNLF with “exact start” is 2nd-order convergent

Proofs rest on local error cancellation of terms that cause order

reduction if is violated. The cancellation fails

at the first CNLF step when IMEX-Euler is used to compute .

Jan verwer

A new splitting (convergence) condition

The non-linear case:

The new condition reads

Jan verwer

Component splitting

Discussed for linear, semi-discrete 1st order wave equations



with S a diagonal matrix satisfying thegeneral ansatz

Jan verwer

The splitting condition

- The common splitting condition requires

- However



- The new splitting condition

is to be interpreted as a discrete spatial integration

which “removes” the factor

Jan verwer


- Stability analysis of IMEX methods normally requires

commuting operators. However,

which isnot true!

- All we can say is that

which regarding stability is necessary for

the LF part and sufficient for the CN part in CNLF

- Experience: runs are stable for the maximal

stable step size for the LF part

Jan verwer

Numerical illustration I

The component splitting

matrix S is chosen in the form

Jan verwer

Illustration I (piecewise uniform grid)

Splitting matrix S such that LF is applied

at the coarse grid and CN at the fine grid.

Factor 10 between coarse & fine grid!

Jan verwer

Illustration I (the splitting conditions)

Plots for time t = 0


Jan verwer

Illustration I (global errors)

Global errors at t = 0.25

Maximal step sizeτ = h with

hthe coarse grid size

1st order

--- : 2nd - order

-o- : CNLF with CN start

-*- : CNLF with IMEX-Euler start

-+- : CN

CNLF with CN start

gives 2nd order

The IMEX-Euler start

causes order reduction !!!


Jan verwer

Illustration I (uniform grid, random S)

Uniform grid and S randomly chosen as

Global errors at t = 0.25

Step sizeτ = h

--- : 2nd order

-o- : CNLF with CN start

-*- : CNLF with IMEX-Euler start

-+- : CN

Results are in line with

those on the non-uniform grid

Jan verwer

Illustration II

2D Maxwell

type problem

on unit square

U(x,y,t = 1)

U(x,y,t = 0)

Jan verwer

Illustration II

Strongly peaked 0.95 < d(x,y) ≤ 100.Through component splitting, we use CN near the peak (d ≥ 1) and LF else-where, to avoid the step size limitation for LF near the peak

A uniform staggered grid and 2nd order differencing with grid size h requires for LF

The following results at t = 1 are obtained with CNLF for

using only a very small amount of implicitly treated points

Jan verwer

Illustration II

CNLF is as accurate as CN

Jan verwer

Illustration II

nnz: number of nonzeros in linear system matrix (sparsity indicator)

Jan verwer


-- Component splitting tests confirm

the new CNLF convergence condition

-- Component splitting can be set up in the

same way for 3D Maxwell

-- But, how practical this is for real

applications, I don’t know yet

  • Login