Report

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