New advances in creating and testing initial configurations of thermomechanical simulations with Tobias Baumann & Arthur Bauville Updated repository (bitbucket.org/geomio/geomio) & new wiki pages (https://bitbucket.org/geomio/geomio/wiki/Home) Motivation How to create a 3D model geometry for thermomechanical simulations? It should incorporate geophysical constraints (grids, points, ...) geological interpretations (cross-sections, maps,...) Motivation Commercial software exists, but ... expensive steep learning curve lack of flexibility (i.e. API to own thermomechanical code) geomIO workflow geomIO workflow geomIO workflow Usually, we start with collecting data We merge all data that we find in maps, and create orthogonal cross-sections The generic mapping toolbox is a great tool for doing that geomIO workflow Interpreting the data We use the vectorgraphics editor (Inkscape) to interpret the data on orthogonal cross-sections 3D example Example Gorleben salt structure Structures were created with data from maps-servers (NIBIS and geotis) that allow you to visualize 2D cross-sections Creating a 3D structure from 2D cross-sections Adding geometry constraints Add a hidden layer and name it as #_<layername> Drag and drop image on the Inkscape canvas Creating a 3D structure from 2D cross-sections Defining the scale of the model Add a Reference layer with a Bezier line object Associate the object with an attribute CoordRef and coordinates as values 0,-3500,2000,-2000 Creating a 3D structure from 2D cross-sections Start interpreting The layername specifies the orientation (EW,NS,HZ) & coordinate normal to the drawing layer Creating a 3D structure from 2D cross-sections Adding a second layer The layername is adapted such that it defines the offset of the layers Creating a 3D structure from 2D cross-sections Adding a second layer ... we duplicate <kbd>ctrl+d</kbd> and modify the object in the new layer. 3D example: rabbit-and-cat tutorial (/tutorials/ Tuto3D_Basic) 3D Rabbit-and-cat tutorial We provide several tutorials in /tutorials. This one illustrates well how geomIO works. 3D example: rabbit-and-cat tutorial (/tutorials/ Tuto3D_Basic) geomIO user code clc, clear, close all %opt = geomIO_Options(); opt.DrawCoordRes = 10; opt.inputFileName = ['Input/RabbitAndCat.HZ.svg']; opt.writeParaview = true; opt.shiftPVobj = [0 0 1]; opt.readLaMEM = true; opt.NewLaMEM = false; opt.LaMEMinputFileName = './Input/RabbitAndCat.dat'; opt.writePolygons = true; [~, Volumes] = run_geomIO(opt); plotVolumes(Volumes,{'CatAndRabbit'}) Initial temperature configuration (tutorials/Example3D_Geodynamics) Basic temperature configuration Simplified plate model Input: SVG file with horizontal (HZ) contours of oceanic and continental lithosphere geomIO associates 1D temperature profiles to each geological unit Initial temperature configuration (tutorials/Example3D_Slab_temperature) Advanced temperature configuration Japanese subduction zones Input: Slab surface geometries (slab contours as HZ.SVG file) + age grid (e.g., MοΏ½ller et al.) Half-space cooling model (depth dependent) Approx. Voronoi diagram for efficient distance-to-surface calculations Gravity modelling with geomIO (/tutorials/Example2D_Gravity/) 2D Gravity model of salt structures Method: Talwani et al. (1959) Works with any 2D geomIO object Gravity modelling with geomIO (/tutorials/Benchmark3D_gravity/) 3D Gravity modelling Method: Talwani & Ewing (1960) works with any 3D geomIO object Summary Summary geomIO repository now contains a detailed wiki & tutorials paper in press (G3 technical report) Initial temperature conditions 2D & 3D Gravity modelling to test interpretations