Numerical modeling of rock deformation: 12 FEM 2D diffusion application

1
Numerical modeling of rock deformation12 FEM
2D diffusion application

• Stefan Schmalholz
• schmalholz_at_erdw.ethz.ch
• NO E61
• AS2009, Thursday 10-12, NO D 11

2
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
Tini1
Tini0
0
0
50
100