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.

5.3. CONTRACTIVITY OF THE NONLINEAR CONTINUUM PROBLEM 125<br />

5.3 Contractivity of the Nonlinear Continuum Problem<br />

5.3.1 Non-contractivity of the fixed-strain split<br />

The coupled problem is contractive relative to the norms · N <strong>and</strong> · T , as shown in<br />

Chapter 3. We investigate whether the fixed-strain split is contractive, or not. We consider<br />

two arbitrary initial conditions as we did in Chapter 3. The fixed-strain split can be written<br />

as<br />

⎡<br />

⎤<br />

⎡<br />

⎤<br />

⎣ dun<br />

dpn ⎦ Apsn −→ ⎣ du∗<br />

dpn+1 ⎦ Ausn −→<br />

⎡<br />

⎣ dun+1<br />

⎤<br />

dpn+1 ⎦, where<br />

⎧<br />

A<br />

⎪⎨<br />

p sn : dm ˙ + Div dv = 0, δdεv ˙ = 0,<br />

Au sn : Div dσ = 0, dp = 0<br />

⎪⎩<br />

⇒ Div dσ ′ = 0.<br />

(5.27)<br />

Equation 5.27 has homogeneous boundary conditions with no source terms. Since the<br />

pressure in the mechanical problem is prescribed, the pressure is not affected by pertur-<br />

bations of the initial condition, thus dp = 0. When solving the flow problem, A p sn, we<br />

obtain<br />

d dχ 2<br />

N<br />

dt<br />

= ∂dχ2 N : dεe<br />

˙ +<br />

∂dεe<br />

∂dχ2 N · ˙<br />

∂dχ<br />

dξ +<br />

∂dξ<br />

2<br />

N dme<br />

˙<br />

∂dme<br />

<br />

= dσ<br />

Ω<br />

′ <br />

dme<br />

: dεe<br />

˙ − dκ · dξ ˙ − M − bdεe,v bdεe,v ˙<br />

ρf,0<br />

+ M<br />

<br />

dme<br />

− bεe,v dme<br />

˙ dΩ<br />

ρf,0 ρf,0<br />

<br />

<br />

= dσ : dε ˙<br />

dp<br />

<br />

+ dm ˙ dΩ − dσ ′ <br />

: dεp<br />

˙ + dκ · dξ ˙ dΩ<br />

=<br />

Ω<br />

ρf,0<br />

Ω <br />

D<br />

<br />

d p≥0 <br />

<br />

dσ : dε ˙<br />

−1<br />

− dv · k µdv<br />

Ω<br />

<br />

dm ˙<br />

−1<br />

= −dv · k µdv from Equation 5.27<br />

dΩ − D d p ≤ 0, (5.28)<br />

where v ∈ [H (div,Ω)] ndim. Note that we consider maximum plastic dissipation, D d p ≥ 0,<br />

when solving the flow problem.

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

Saved successfully!

Ooh no, something went wrong!