16.08.2013 Views

Truncation error in SPH

Truncation error in SPH

Truncation error in SPH

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Non-uniform spac<strong>in</strong>g<br />

∑<br />

∂A<br />

− A W<br />

b bΔ xb<br />

=<br />

∂<br />

′<br />

b x = x<br />

2<br />

( ˆ )<br />

h<br />

3<br />

− A′ 1<br />

ˆ<br />

a Wds − + A′′′ a s W ′ ds +<br />

6<br />

3 4<br />

1 ⎡ ⎛⎛ Δx ⎞ ⎞ Aa<br />

⎛ 2 1 ⎞ ⎛⎛ Δx<br />

⎞ ⎞⎤<br />

− ⎢Aaδ O ⎜ + δ + O ⎥<br />

h ⎜⎜ h<br />

⎟ ⎟<br />

⎜ ⎟<br />

2<br />

⎜<br />

12<br />

⎟ ⎜⎜ h<br />

⎟<br />

⎢ ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ ⎟<br />

⎣ ⎝ ⎠ ⎝ ⎠⎥⎦<br />

3<br />

⎡ ⎛⎛ Δx<br />

⎞ ⎞⎤<br />

A′ aδO<br />

⎜<br />

⎜⎜ ⎟<br />

h<br />

⎟ ⎟<br />

4<br />

⎡ A′′ a ⎛ Δx ⎞ A′′<br />

⎛<br />

a ⎛ Δx<br />

⎞ ⎞⎤<br />

First <strong>SPH</strong>ERIC Workshop, Rome, 2006<br />

x<br />

− ⎢ ⎥<br />

⎢⎣ ⎝⎝ ⎠ ⎠⎥⎦<br />

a<br />

∫ ∫ …<br />

− h ⎢ δO<br />

O<br />

2<br />

⎜ +<br />

h<br />

⎟ ⎜ ⎟⎥<br />

2 ⎜⎜ h<br />

⎟ −…<br />

⎢ ⎝ ⎠ ⎝ ⎠ ⎟<br />

⎣ ⎝ ⎠⎦⎥<br />

smooth<strong>in</strong>g <strong>error</strong><br />

discretisation<br />

<strong>error</strong>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!