• No results found

The effect of elastic compressibility in geodynamic

N/A
N/A
Protected

Academic year: 2022

Share "The effect of elastic compressibility in geodynamic "

Copied!
1
0
0

Bezig met laden.... (Bekijk nu de volledige tekst)

Hele tekst

(1)

The effect of elastic compressibility in geodynamic

models using a benchmark of a sinking cylinder in a bounded medium

Walter Austmann

1

, Rob Govers

1

1

Department of Earth Sciences, Utrecht University, Utrecht, The Netherlands

E-mail: w.f.austmann@uu.nl

Introduction

Although it is commonly accepted that the Earth is compressible, this effect is neglect- ed in most numerical models. For purely viscous rheologies, this is justified, since the bulk viscosity of the Earth is negligi- ble. However, for visco-elastic rheologies, the elastic deformation allows volumetric changes. To study the effect of elastic com- pressibility, a well-defined benchmark in fluid mechanics is used: a sinking cylinder in a homogenous fluid. An analytical solu- tion exists in the case of an incompressible purely viscous fluid and is compared to our visco-elastic solution.

a

2R

g H

x y

Figure 1: Physical sketch of the model. In the reference model H=25cm; R=0.5cm;

a=0.25cm. After:DeFranco (2008)

Model Setup

A rigid cylinder is placed midway in a ho- mogeneous visco-elastic fluid. The fluid is bounded by impermeable walls. The tank is 25cm high and 1cm wide. The cylinder is placed at a distance of 1,5cm from the top of the domain. A sketch is given in Fig. 1 and the model parameters in Table 1.

This setup is solved using the lagrangian finite element package G-TECTON con- sisting of both triangular and quadrilateral elements. A remeshing is used to avoid nu- merical errors due to high deformations.

Benchmark description

In the analytical benchmark, the steady

state Stokes equations are solved for a rigid cylinder moving with a constant velocity

through a homogeneous fluid between two parallel impermeable walls. The fluid ex- tends infinitely in the direction of motion.

The resulting flow field (Fig. 2) is only non- uniform in a small confined region around the cylinder (~2 radius width).

From the flow field, the dimensionless drag- force is calculated. In our models, it is as- sumed that when the cylinder is in steady state, the dragforce balances the buoyancy force. This implies that the fluid is not self- gravitating.

Figure 2:

Streamfunction of the bench-

mark around a fixed cylinder

(Liu et. al, 1998)

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5

Velocity (cm s−1 )

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0

time (s)

Steady state solution: Champmartin (2007)

ν=0.49 ν=0.45 ν=0.40 ν=0.35 ν=0.30 ν=0.25 ν=0.20 ν=0.15 ν=0.10

0

5

10

15

20

25

y (cm)

0 1

x (cm)

0 1

x (cm)

0

5

10

15

20

25

y (cm)

0 1

x (cm)

0 1

x (cm)

0

5

10

15

20

25

y (cm)

0 1

x (cm)

0 1

x (cm)

−0.70

−0.56

−0.42

−0.28

−0.14 0.00 0.14 0.28 0.42 0.56 0.70

v x (cm s−1 )

0

5

10

15

20

25

y (cm)

0 1

x (cm)

0 1

x (cm)

0 1 0 1

0 1 0 1

0 1

0 1 −1.0

−0.8

−0.6

−0.4

−0.2 0.0 0.2 0.4 0.6 0.8 1.0

v y (cm s−1 )

0 1 0 1

Figure 3: Velocity field after 5 sec- onds. The cylinder has a velocity

close to steady state. Horizontal velocity on the left and vertical velocity on the right. Different

scale bars are used for the two components.

Figure 6: Velocity profiles for different Poisson ratios. A Poisson ratio of 0.5 represents an incompressible flow. Each model run reaches an almost constant velocity at the end of the run. The highest Poisson ratios are in good agreement with the analytical solution. Due to elastic effects, the initial velocity is large and it monotonically decreasing with time.

For Earth-like poisson ratios (i.e. n=0.3) the difference with the incompressible (analyti- cal) solution is 40%.

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

y (cm)

0.00 0.25 0.50 0.75 1.00

Velocity (cm s−1)

t=0.5 s t=1.0 s t=1.5 s t=2.0 s t=2.5 s t=3.0 s t=3.5 s t=4.0 s t=4.5 s t=5.0 s 0

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

y (cm)

0.00 0.25 0.50 0.75 1.00

Velocity (cm s−1)

Figure 4: Time evolution of the vertical velocity at x=0.5cm. For near steady state velocities, the vertical motion of the fluid approaches a linear profile. The gradient is dependent on the distance between the cyl- inder and the domain boundaries. Thus, as long as the distance travelled by the cylin- der is much smaller than the distance be- tween the cylinder and the boundary, the cylinder approaches a steady state velocity.

Analysis

In the incompressible solution, the downward mo- tion of the cylinder is balanced by the upward flow

around the cylinder. In the compressible solu- tion, a part of the motion of the cylinder is ac-

comodated by the compressibility of the fluid below. Thus the upward motion is smaller and

the dragforce is decreased. Compressibility in- creases therefore the sinking velocity of the

cylinder.

Discussion

The motion of tectonic plates is the effect of driving and resisting forces acting on the plate. The driving forces (e.g. slab pull, ridge push and continental lith- ospheric body forces) are relatively well constrained.

The resisting forces are much harder to quantify and are only implicit in the numerical models. By using an incompressible viscous rheology instead of a (more

realistic) compressible visco-elastic rheology, these resistive forces are overestimated. This causes an un- derestimation of the deformation rates (e.g. predicted plate velocities, sinking velocities of (detached) slabs or topographic uplift). In simple models this can be as large as 40% for Earth-like Poisson ratios.

Conclusions

• Compressible visco-elastic flow increases the steady state velocity by 40% in the case of a sinking cylin- der for Earth-like Poisson ratios.

• The motion of the fluid is no longer confined to a small region around the cylinder

• The boundaries of fluid have large effect on the mo- tion: the closer a boundary, the larger the resistance on the motion

• Steady state results are achieved when the cylinder is initially (1) close to the top domain to damp out the visco-elastic overshoot and (2) far from the bottom so that the cylinder is now slowed down after obtaining the steady state velocity.

Results

Material Cylinder Fluid

Viscosity (Pa s) 1.64*10

20

1.64 Young Modulus (Pa) 3.5*10

13

350

Poisson ratio 0.4 0.4

Density (kg m

-3

) 7161 0

Table 1: model parameter of the reference model

−0.10

−0.05 0.00

ε yy (s−1 )

0.2 0.4

Poisson ratio

−0.10

−0.05 0.00

ε yy (s−1 )

0.2 0.4

Poisson ratio

.

0.0 0.1 0.2 0.3

∆ V (cm s−1 )

0.0 0.3 0.6 0.9

Aspect ratio 0.0

0.1 0.2 0.3

∆ V (cm s−1 )

0.0 0.3 0.6 0.9

Aspect ratio 0.00

0.25 0.50 0.75 1.00

∆ V (cm s−1 )

0.0 0.2 0.4 0.6

Poisson ratio 0.00

0.25 0.50 0.75 1.00

∆ V (cm s−1 )

0.0 0.2 0.4 0.6

Poisson ratio

Figure 5: For an increasing Poisson ratio (1) the flow becomes more confined (left: the av- erage vertical strain rate in the region below the cylinder approaches zero) and (2) the ver- tical velocity of the cylinder approaches the analytical solution (n=0.5) (middle)

Vertical velocity differences decrease for increasing a/R

Referenties

GERELATEERDE DOCUMENTEN

Kees Lokhorst (l) en Ferry Leenstra helpen ondernemers hun ideeën voor een vleeskuiken- houderij met toekomst, uit te werken.

Op basis van de gemaakte keuzes in stap 1 en stap 2 kan het nieuwe niveau voor de TW voor alle steekjaren na revisie worden berekend.. Berekenen van het nieuwe niveau

On the other hand, to present preliminary work on a classification problem of fire identification in forests and a counting problem of estimating number of houses in urban areas, based

In het bij- zonder gaat mijn dank uit naar André Jansen en Henk Mulder voor de assistentie tijdens de..

Lorsqu'il fut procédé, en vue du placement d'une chaufferie, au ereasement de tranchées dans la partie occidentale de la collégiale Saint-Feuillen à Fosse, et

• The final author version and the galley proof are versions of the publication after peer review.. • The final published version features the final layout of the paper including

plantarum 423-neg pGKV-plaAPrV2mChB is cultured with a manganese concentration between 295.00 and 12.76 µM, transcription appears to reach a relatively higher

Hoofdstuk 4 Beschrijving van de sporen, structuren en vondsten In totaal werden er 62 individuele grondsporen geregistreerd, waarvan 18 exemplaren van natuurlijke aard