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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

5.5. ERROR ESTIMATION 133<br />

Using the fully coupled method, the problem is written as<br />

⎡<br />

F L<br />

⎣<br />

−Lt ⎤⎡<br />

⎦⎣<br />

K<br />

<br />

A<br />

p<br />

⎤n+1<br />

⎡ ⎤⎡<br />

⎦<br />

u<br />

L Q<br />

− ⎣ ⎦⎣<br />

0 0<br />

<br />

B<br />

p<br />

⎤<br />

⎦<br />

u<br />

⎡<br />

x = ⎣ p<br />

⎤<br />

⎦, F = Q + ∆tT,<br />

u<br />

<strong>Sequential</strong> methods decompose the matrix A into<br />

⎡<br />

F + R 0<br />

⎣<br />

−Lt ⎤⎡<br />

⎦⎣<br />

K<br />

p<br />

⎤<br />

⎦<br />

u<br />

⎡ ⎤⎡<br />

L Q<br />

− ⎣ ⎦⎣<br />

0 0<br />

p<br />

⎤n<br />

⎡<br />

⎦ =<br />

u<br />

n+1,k+1<br />

⎣ fp<br />

fu<br />

n<br />

=<br />

⎡<br />

⎣ fp<br />

fu<br />

⎤<br />

⎦<br />

<br />

f<br />

⎡ ⎤⎡<br />

R −L<br />

− ⎣ ⎦⎣<br />

0 0<br />

p<br />

⎤<br />

⎦<br />

u<br />

⎤<br />

⎦<br />

n+1<br />

n+1<br />

n+1,k<br />

, (5.56)<br />

, (5.57)<br />

where R is the relaxation matrix in a sequential method. R depends on the specific sequen-<br />

tial method used. Subtracting Equation 5.57 from Equation 5.56, the errors of pressure <strong>and</strong><br />

displacement are written as<br />

⎡<br />

⎣ efs,p<br />

efs,u<br />

⎤<br />

⎦<br />

n+1,k+1<br />

=<br />

=<br />

⎡<br />

F + R 0<br />

⎣<br />

−Lt ⎤<br />

⎦<br />

K<br />

⎡ ⎤ ⎡<br />

L Q<br />

+ ⎣ ⎦<br />

0 0<br />

⎡ ⎤⎡<br />

⎜ R −L<br />

⎝⎣<br />

⎦<br />

0 0<br />

⎤n⎞<br />

−1 ⎛<br />

⎣ efs,p<br />

efs,u<br />

⎦<br />

⎠<br />

⎛<br />

⎣ efs,p<br />

efs,u<br />

⎤<br />

⎦<br />

n+1,k<br />

⎡<br />

(F + R)<br />

⎣<br />

−1<br />

K<br />

0<br />

−1Lt (F + R) −1 K−1 <br />

⎤⎜⎡<br />

⎤⎡<br />

⎜<br />

⎦⎜<br />

R −L<br />

⎜⎣<br />

⎦⎣<br />

⎜<br />

⎝ 0 0<br />

<br />

M<br />

N<br />

efs,p<br />

⎤n+1,k<br />

⎦<br />

efs,u<br />

⎡ ⎤ ⎡<br />

L Q<br />

+ ⎣ ⎦ ⎣<br />

0 0<br />

efs,p<br />

⎤n⎞<br />

⎦ ⎠ . (5.58)<br />

efs,u

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

Saved successfully!

Ooh no, something went wrong!