1 / 2

# Numerical modeling of rock deformation: 12 FEM 2D diffusion application - PowerPoint PPT Presentation

Numerical modeling of rock deformation: 12 FEM 2D diffusion application. Stefan Schmalholz [email protected] NO E61 AS2009, Thursday 10-12, NO D 11. Contact metamorphic aureole.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

## PowerPoint Slideshow about ' Numerical modeling of rock deformation: 12 FEM 2D diffusion application' - junior

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

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

Stefan Schmalholz

NO E61

AS2009, Thursday 10-12, NO D 11

Numerical modeling of rock deformation: FEM 2D Elasticity. Stefan Schmalholz, ETH Zurich

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