Truncation error in SPH
Truncation error in SPH
Truncation error in SPH
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>