06.08.2013 Views

Sequential Methods for Coupled Geomechanics and Multiphase Flow

Sequential Methods for Coupled Geomechanics and Multiphase Flow

Sequential Methods for Coupled Geomechanics and Multiphase Flow

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

4.1. ERROR ESTIMATION 83<br />

displacement between the fully coupled method <strong>and</strong> drained split are written as<br />

⎡<br />

⎣ efs,u<br />

where efs =<br />

efs,p<br />

⎤<br />

⎦<br />

n+1,k+1<br />

⎡ ⎤ ⎡ ⎤⎡<br />

K 0<br />

= ⎣ ⎦<br />

L F<br />

⎜ 0 Lt<br />

⎝⎣<br />

⎦⎣<br />

0 0<br />

efs,u<br />

⎤<br />

⎦<br />

efs,p<br />

⎡<br />

K<br />

= ⎣<br />

−1 −F<br />

0<br />

−1LK−1 F−1 ⎤⎛⎡<br />

⎤⎡<br />

⎦⎜<br />

0 Lt<br />

⎝⎣<br />

⎦<br />

0 0<br />

⎡ ⎤⎡<br />

0 0<br />

+ ⎣ ⎦⎣<br />

L Q<br />

efs,u<br />

⎤n⎞<br />

⎦ ⎠<br />

efs,p<br />

⎡<br />

⎤⎡<br />

⎤<br />

=<br />

−1 ⎛<br />

⎣ 0 K−1Lt 0 −F−1LK−1Lt ⎡<br />

⎦<br />

<br />

D<br />

⎣ efs,u<br />

efs,p<br />

e n+1,niter<br />

fs<br />

⎤<br />

⎦, <strong>and</strong> e n fs<br />

⎣ efs,u<br />

efs,p<br />

⎦<br />

n+1,k<br />

n+1,k<br />

⎣ efs,u<br />

⎡ ⎤ ⎡<br />

0 0<br />

+ ⎣ ⎦<br />

L Q<br />

efs,p<br />

⎡<br />

+ ⎣<br />

= en,niter<br />

fs . Then e n+1,niter<br />

fs<br />

= D niter<br />

niter <br />

n+1,0<br />

e +<br />

l=1<br />

⎤<br />

⎦<br />

D l−1 He n,niter<br />

fs<br />

= D niter n+1<br />

(xf − xnf ) + (Dniter niter<br />

+<br />

where en+1,0 is written as (x n+1<br />

f −xn f )+(xn f −xn,niter s<br />

in terms of the initial error e 0,0<br />

fs , Dniter , <strong>and</strong> S as<br />

e n+1,niter<br />

fs<br />

=<br />

n<br />

l=0<br />

D niter S l (x n+1−l<br />

f<br />

n+1,k<br />

0 0<br />

F −1 L F −1 Q<br />

⎣ efs,u<br />

efs,p<br />

⎤⎡<br />

⎦<br />

<br />

H<br />

can be written as<br />

l=1<br />

<br />

S<br />

⎤<br />

⎦<br />

⎣ efs,u<br />

efs,p<br />

<br />

D l−1 H) e n,niter<br />

fs ,<br />

n ⎞<br />

⎟<br />

⎠<br />

⎤<br />

⎦<br />

n<br />

,<br />

(4.4)<br />

(4.5)<br />

). By recursion, e n+1,niter is expressed<br />

− x n−l<br />

f )+S n+1 e 0,0<br />

fs<br />

. (4.6)

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

Saved successfully!

Ooh no, something went wrong!