1 / 2

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

Numerical modeling of rock deformation: 12 FEM 2D diffusion application. Stefan Schmalholz schmalholz@erdw.ethz.ch NO E61 AS2009, Thursday 10-12, NO D 11. Contact metamorphic aureole.

junior
Download Presentation

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

An Image/Link below is provided (as is) to download presentation 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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Numerical modeling of rock deformation:12 FEM 2D diffusion application Stefan Schmalholz schmalholz@erdw.ethz.ch NO E61 AS2009, Thursday 10-12, NO D 11 Numerical modeling of rock deformation: FEM 2D Elasticity. Stefan Schmalholz, ETH Zurich

  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 Tini=0 Tini=1 0 0 50 100 Numerical modeling of rock deformation: FEM 2D Elasticity. Stefan Schmalholz, ETH Zurich

More Related