Transcript Slide 1

Numerical modeling of rock
deformation:
12 FEM 2D
diffusion application
Stefan Schmalholz
[email protected]
NO E61
AS2009, Thursday 10-12, NO D 11
Numerical modeling of rock deformation: FEM 2D Elasticity. Stefan Schmalholz, ETH Zurich
Contact metamorphic aureole
Use the Matlab file “ML_FEM_2D_Diff_numint” from the course page which includes a
finite element code that solves the 2D transient diffusion equation.
Set up a model as displayed below. Set the initial temperatures as indicated. Use a
constant value of 1 for kappa. Heat production is zero.
An intrusion (blue) with temperature 1 is emplaced in rocks with temperature 0. A mineral
reaction takes place in the rocks if the temperature is larger than 0.3. Determine the
area in the rocks in which the reaction takes place. Do not set any boundary
conditions and explain what this physically means for the investigated problem.
Determine the rock area where reactions take place for the case when the temperature at
the top boundary is set to 0. Compare and discuss the results.
Before you do the simulations guess the result. For what boundary conditions you expect
a larger contact aureole?
50
60
30
30
Tini=0
Tini=1
0
0
50
100
Numerical modeling of rock deformation: FEM 2D Elasticity. Stefan Schmalholz, ETH Zurich